From c9a55080a79dfe508962ef2f65156119828f9559 Mon Sep 17 00:00:00 2001 From: David van der Spoel Date: Wed, 12 Feb 2020 14:53:33 +0100 Subject: [PATCH] Fixed bug in orires code that was checking too strict. Added one test to see that input is accepted. Fixes #3375 Change-Id: I0f64b386477b9a72f18c95fd977ac8cf5dbeeccd --- docs/release-notes/2020/2020.1.rst | 8 + src/gromacs/listed_forces/orires.cpp | 4 +- src/programs/mdrun/tests/CMakeLists.txt | 3 +- src/programs/mdrun/tests/orires.cpp | 93 ++ src/testutils/simulationdatabase/orires_1lvz.gro | 218 +++ src/testutils/simulationdatabase/orires_1lvz.ndx | 23 + src/testutils/simulationdatabase/orires_1lvz.top | 1820 ++++++++++++++++++++++ 7 files changed, 2166 insertions(+), 3 deletions(-) create mode 100644 src/programs/mdrun/tests/orires.cpp create mode 100644 src/testutils/simulationdatabase/orires_1lvz.gro create mode 100644 src/testutils/simulationdatabase/orires_1lvz.ndx create mode 100644 src/testutils/simulationdatabase/orires_1lvz.top diff --git a/docs/release-notes/2020/2020.1.rst b/docs/release-notes/2020/2020.1.rst index 0c7c3b7f36..7aaf378c74 100644 --- a/docs/release-notes/2020/2020.1.rst +++ b/docs/release-notes/2020/2020.1.rst @@ -72,6 +72,14 @@ the new modular simulator approach) could not be read. :issue:`3377` +Avoid overzealous program abort with orientation restraints +""""""""""""""""""""""""""""""""""""""""""""""""""""""""""" + +It could happen that mdrun would abort on checking orientation restraints in multiple +molecules even though no restraints where applied to them. + +:issue:`3375` + Fixes for ``gmx`` tools ^^^^^^^^^^^^^^^^^^^^^^^ diff --git a/src/gromacs/listed_forces/orires.cpp b/src/gromacs/listed_forces/orires.cpp index c7d326571d..2867c07a15 100644 --- a/src/gromacs/listed_forces/orires.cpp +++ b/src/gromacs/listed_forces/orires.cpp @@ -3,7 +3,7 @@ * * Copyright (c) 1991-2000, University of Groningen, The Netherlands. * Copyright (c) 2001-2004, The GROMACS development team. - * Copyright (c) 2013,2014,2015,2016,2017,2018,2019, by the GROMACS development team, led by + * Copyright (c) 2013,2014,2015,2016,2017,2018,2020, by the GROMACS development team, led by * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl, * and including many others, as listed in the AUTHORS file in the * top-level source directory and at http://www.gromacs.org. @@ -122,7 +122,7 @@ void init_orires(FILE* fplog, int nmol; while (const InteractionLists* il = gmx_mtop_ilistloop_next(iloop, &nmol)) { - if (nmol > 1) + if (nmol > 1 && (*il)[F_ORIRES].size() > 0) { gmx_fatal(FARGS, "Found %d copies of a molecule with orientation restrains while the current " diff --git a/src/programs/mdrun/tests/CMakeLists.txt b/src/programs/mdrun/tests/CMakeLists.txt index c6bda33065..f2e7e2d8bb 100644 --- a/src/programs/mdrun/tests/CMakeLists.txt +++ b/src/programs/mdrun/tests/CMakeLists.txt @@ -1,7 +1,7 @@ # # This file is part of the GROMACS molecular simulation package. # -# Copyright (c) 2013,2014,2015,2016,2017,2018,2019, by the GROMACS development team, led by +# Copyright (c) 2013,2014,2015,2016,2017,2018,2019,2020, by the GROMACS development team, led by # Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl, # and including many others, as listed in the AUTHORS file in the # top-level source directory and at http://www.gromacs.org. @@ -58,6 +58,7 @@ gmx_add_gtest_executable( initialconstraints.cpp interactiveMD.cpp outputfiles.cpp + orires.cpp pmetest.cpp simulator.cpp swapcoords.cpp diff --git a/src/programs/mdrun/tests/orires.cpp b/src/programs/mdrun/tests/orires.cpp new file mode 100644 index 0000000000..0705af96be --- /dev/null +++ b/src/programs/mdrun/tests/orires.cpp @@ -0,0 +1,93 @@ +/* + * This file is part of the GROMACS molecular simulation package. + * + * Copyright (c) 2020, by the GROMACS development team, led by + * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl, + * and including many others, as listed in the AUTHORS file in the + * top-level source directory and at http://www.gromacs.org. + * + * GROMACS is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public License + * as published by the Free Software Foundation; either version 2.1 + * of the License, or (at your option) any later version. + * + * GROMACS is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with GROMACS; if not, see + * http://www.gnu.org/licenses, or write to the Free Software Foundation, + * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. + * + * If you want to redistribute modifications to GROMACS, please + * consider that scientific software is very special. Version + * control is crucial - bugs must be traceable. We will be happy to + * consider code for inclusion in the official distribution, but + * derived work must not be called official GROMACS. Details are found + * in the README & COPYING files - if they are missing, get the + * official version at http://www.gromacs.org. + * + * To help us fund GROMACS development, we humbly ask that you cite + * the research papers on the package. Check out http://www.gromacs.org. + */ + +/*! \internal \file + * \brief + * Test for MD with orientation restraints + * + * \author David van der Spoel + * \ingroup module_mdrun_integration_tests + */ +#include "gmxpre.h" + +#include "moduletest.h" + +namespace gmx +{ +namespace test +{ + +class OriresTestFixture : public MdrunTestFixture +{ +protected: + OriresTestFixture(); + ~OriresTestFixture() override; +}; + +OriresTestFixture::OriresTestFixture() {} + +OriresTestFixture::~OriresTestFixture() {} + +//! Test fixture for mdrun with orires +typedef gmx::test::OriresTestFixture OriresTest; + +/* Check whether the orires function works. */ +TEST_F(OriresTest, OriresCanRun) +{ + runner_.useTopGroAndNdxFromDatabase("orires_1lvz"); + const std::string mdpContents = R"( + dt = 0.002 + nsteps = 10 + tcoupl = Berendsen + tc-grps = System + tau-t = 0.5 + ref-t = 300 + constraints = h-bonds + cutoff-scheme = Verlet + orire = Yes + orire-fitgrp = backbone + )"; + runner_.useStringAsMdpFile(mdpContents); + + EXPECT_EQ(0, runner_.callGrompp()); + + ::gmx::test::CommandLine oriresCaller; + + // Do an mdrun with ORIRES enabled + ASSERT_EQ(0, runner_.callMdrun(oriresCaller)); +} + +} // namespace test +} // namespace gmx diff --git a/src/testutils/simulationdatabase/orires_1lvz.gro b/src/testutils/simulationdatabase/orires_1lvz.gro new file mode 100644 index 0000000000..40ad092f36 --- /dev/null +++ b/src/testutils/simulationdatabase/orires_1lvz.gro @@ -0,0 +1,218 @@ +GUANINE NUCLEOTIDE-BINDING PROTEIN G(T), ALPHA-1 + 215 + 1ILE N 1 2.423 2.571 1.751 0.0936 -0.3219 -0.0355 + 1ILE H1 2 2.471 2.606 1.669 -2.9915 0.2125 -1.7758 + 1ILE H2 3 2.346 2.637 1.754 0.0150 -0.2794 -1.7958 + 1ILE H3 4 2.488 2.574 1.828 -1.8243 0.9739 1.6753 + 1ILE CA 5 2.390 2.431 1.722 -0.0040 -0.3102 0.0167 + 1ILE HA 6 2.323 2.434 1.636 -0.2465 -1.2965 0.1610 + 1ILE CB 7 2.514 2.348 1.691 0.1548 -0.0167 -0.1390 + 1ILE HB 8 2.560 2.401 1.607 3.5618 -0.2132 1.4412 + 1ILE CG2 9 2.614 2.337 1.805 0.0404 -1.1364 -0.1327 + 1ILE HG21 10 2.712 2.313 1.764 -1.2409 -2.3870 -2.6543 + 1ILE HG22 11 2.622 2.425 1.870 -1.9742 0.2679 -1.6710 + 1ILE HG23 12 2.585 2.259 1.876 1.0656 -1.6613 -0.2640 + 1ILE CG1 13 2.466 2.213 1.638 -0.1684 -0.4516 1.2294 + 1ILE HG11 14 2.393 2.229 1.559 2.4089 2.0711 -0.8483 + 1ILE HG12 15 2.429 2.151 1.720 -1.2108 -0.2315 0.9452 + 1ILE CD 16 2.588 2.138 1.585 -0.4056 -0.1240 0.1927 + 1ILE HD1 17 2.656 2.208 1.536 2.6377 -2.2144 1.1352 + 1ILE HD2 18 2.636 2.088 1.669 0.0941 -0.2702 -0.1801 + 1ILE HD3 19 2.542 2.076 1.508 -0.8460 1.5249 -0.9365 + 1ILE C 20 2.302 2.367 1.828 -0.1261 -0.0840 0.0509 + 1ILE O 21 2.217 2.288 1.788 -0.1254 -0.2860 0.4413 + 2ARG N 22 2.313 2.403 1.957 0.3813 -0.2249 0.0459 + 2ARG H 23 2.368 2.486 1.972 0.5336 0.3645 -3.1220 + 2ARG CA 24 2.237 2.338 2.061 0.1502 0.1391 0.1040 + 2ARG HA 25 2.259 2.231 2.065 -2.6348 -0.5011 1.6122 + 2ARG CB 26 2.286 2.415 2.183 -0.6640 0.4800 0.2203 + 2ARG HB1 27 2.248 2.516 2.168 -1.9093 -0.1917 -1.4780 + 2ARG HB2 28 2.395 2.415 2.185 -0.6732 0.8448 1.2487 + 2ARG CG 29 2.247 2.355 2.318 0.1107 -0.1795 0.1531 + 2ARG HG1 30 2.276 2.250 2.327 -0.1239 -0.0974 2.6200 + 2ARG HG2 31 2.139 2.355 2.332 -0.3872 -1.6470 -2.7549 + 2ARG CD 32 2.306 2.435 2.433 0.4397 -0.4298 0.1617 + 2ARG HD1 33 2.408 2.466 2.407 -0.0073 1.8954 0.9116 + 2ARG HD2 34 2.304 2.379 2.526 1.7767 -2.1717 -0.7954 + 2ARG NE 35 2.221 2.552 2.456 0.6417 -0.2564 0.0480 + 2ARG HE 36 2.149 2.567 2.387 0.8932 1.5191 0.1445 + 2ARG CZ 37 2.237 2.639 2.557 -0.2182 0.1165 -0.1247 + 2ARG NH1 38 2.325 2.632 2.658 -0.2428 -0.1092 -0.1168 + 2ARG HH11 39 2.403 2.568 2.659 0.0097 0.1379 2.7157 + 2ARG HH12 40 2.307 2.692 2.737 -0.0540 1.5964 -1.3081 + 2ARG NH2 41 2.150 2.740 2.566 0.7550 1.0205 -0.6821 + 2ARG HH21 42 2.092 2.756 2.485 -0.5254 -0.5336 -0.1216 + 2ARG HH22 43 2.168 2.809 2.638 3.1877 0.7170 -0.9288 + 2ARG C 44 2.085 2.344 2.050 0.1568 -0.7156 -0.5178 + 2ARG O 45 2.018 2.254 2.100 0.1123 -0.5463 -0.2711 + 3GLU N 46 2.030 2.447 1.986 -0.0863 -0.6508 -0.2018 + 3GLU H 47 2.107 2.506 1.957 2.5076 -2.7563 1.8233 + 3GLU CA 48 1.892 2.459 1.943 -0.0791 0.2737 0.0090 + 3GLU HA 49 1.837 2.448 2.037 -0.0446 -2.3690 -0.2034 + 3GLU CB 50 1.865 2.601 1.894 -0.2581 0.1853 -0.1468 + 3GLU HB1 51 1.917 2.622 1.801 -1.4334 0.8957 -0.6630 + 3GLU HB2 52 1.907 2.668 1.969 0.3663 1.1113 -1.2777 + 3GLU CG 53 1.716 2.632 1.880 -0.2972 0.0131 -0.1084 + 3GLU HG1 54 1.660 2.592 1.965 -0.5576 0.8787 0.1337 + 3GLU HG2 55 1.681 2.577 1.793 1.2727 -2.5557 0.7537 + 3GLU CD 56 1.682 2.778 1.861 0.2056 0.1281 -0.1580 + 3GLU OE1 57 1.769 2.867 1.852 0.0423 0.3023 -0.0305 + 3GLU OE2 58 1.562 2.800 1.831 0.1337 -0.3582 -0.2334 + 3GLU C 59 1.852 2.355 1.840 0.0033 0.3192 -0.0686 + 3GLU O 60 1.749 2.291 1.855 0.2870 -0.1890 -0.2784 + 4ASN N 61 1.934 2.319 1.741 0.0582 -0.4105 0.2369 + 4ASN H 62 2.011 2.383 1.724 -0.7026 0.2943 -0.6897 + 4ASN CA 63 1.924 2.204 1.654 0.6640 -0.3126 0.0317 + 4ASN HA 64 1.832 2.201 1.597 -0.0998 -0.7187 1.2592 + 4ASN CB 65 2.035 2.216 1.550 0.7497 -0.1563 0.1416 + 4ASN HB1 66 2.138 2.213 1.585 0.6873 -2.0392 0.2656 + 4ASN HB2 67 2.020 2.315 1.507 0.5996 -0.2743 -0.0752 + 4ASN CG 68 2.028 2.114 1.438 0.6663 0.3106 -0.2822 + 4ASN OD1 69 2.086 2.007 1.454 0.4121 0.2228 0.0416 + 4ASN ND2 70 1.940 2.128 1.338 0.0395 -0.1433 0.2007 + 4ASN HD21 71 1.873 2.204 1.335 0.1268 -0.0060 1.2169 + 4ASN HD22 72 1.942 2.064 1.260 0.8502 -0.3474 0.3858 + 4ASN C 73 1.927 2.076 1.737 0.0386 -0.2349 0.1798 + 4ASN O 74 1.840 1.991 1.717 -0.0458 -0.3335 0.9355 + 5LEU N 75 2.021 2.063 1.831 -0.0687 0.6090 0.4036 + 5LEU H 76 2.094 2.132 1.831 -0.4774 1.1622 -2.5375 + 5LEU CA 77 2.032 1.944 1.913 -0.0413 0.5320 0.2883 + 5LEU HA 78 2.027 1.857 1.848 -2.4141 1.5946 -1.0762 + 5LEU CB 79 2.166 1.955 1.987 0.1109 -0.0905 0.1059 + 5LEU HB1 80 2.170 1.872 2.058 2.6174 0.6698 0.9454 + 5LEU HB2 81 2.162 2.047 2.045 0.8089 0.5599 -0.8616 + 5LEU CG 82 2.288 1.952 1.895 0.0512 -0.2565 0.0305 + 5LEU HG 83 2.285 2.039 1.829 -0.4586 0.4334 0.9230 + 5LEU CD1 84 2.408 1.964 1.988 0.0055 0.1775 0.0372 + 5LEU HD11 85 2.407 1.877 2.053 0.0883 0.1450 -0.0054 + 5LEU HD12 86 2.501 1.962 1.932 -0.4579 -1.7524 -0.7446 + 5LEU HD13 87 2.402 2.060 2.041 1.9484 -0.0759 0.8026 + 5LEU CD2 88 2.309 1.825 1.814 0.4229 -0.0957 -0.1266 + 5LEU HD21 89 2.222 1.814 1.748 -1.5139 2.0474 1.8761 + 5LEU HD22 90 2.406 1.834 1.765 0.0073 -2.6964 -1.5970 + 5LEU HD23 91 2.313 1.735 1.875 0.0387 -1.7936 -2.4350 + 5LEU C 92 1.915 1.936 2.010 -0.2894 0.1139 -0.0464 + 5LEU O 93 1.876 1.826 2.049 -0.1455 -0.1880 -0.7352 + 6LYS N 94 1.842 2.045 2.034 -0.6901 0.0243 -0.8226 + 6LYS H 95 1.882 2.132 2.004 0.8199 -0.9562 -1.7789 + 6LYS CA 96 1.719 2.046 2.110 -0.2719 0.1430 -0.1383 + 6LYS HA 97 1.741 1.993 2.203 0.9730 -0.6169 -0.8467 + 6LYS CB 98 1.672 2.187 2.146 -0.2506 0.0636 0.2090 + 6LYS HB1 99 1.625 2.232 2.059 0.9820 2.6593 0.7788 + 6LYS HB2 100 1.755 2.252 2.173 0.1270 -0.8653 1.3155 + 6LYS CG 101 1.583 2.188 2.270 -0.3322 -0.0661 0.1524 + 6LYS HG1 102 1.651 2.161 2.351 2.0085 2.7494 -0.6671 + 6LYS HG2 103 1.497 2.122 2.264 0.2178 -0.8095 0.2702 + 6LYS CD 104 1.541 2.333 2.293 0.4968 0.1971 0.0589 + 6LYS HD1 105 1.522 2.376 2.195 1.7394 0.6881 0.0130 + 6LYS HD2 106 1.620 2.389 2.343 0.0016 0.7250 0.2662 + 6LYS CE 107 1.416 2.344 2.380 0.7798 0.1243 0.4757 + 6LYS HE1 108 1.340 2.270 2.355 -1.1187 1.9929 0.4302 + 6LYS HE2 109 1.363 2.437 2.361 -1.3895 -1.1540 -0.0331 + 6LYS NZ 110 1.441 2.318 2.523 0.6853 0.4103 0.5445 + 6LYS HZ1 111 1.361 2.279 2.570 -0.0921 1.4011 0.0796 + 6LYS HZ2 112 1.513 2.248 2.533 0.3965 0.2699 1.8626 + 6LYS HZ3 113 1.474 2.400 2.572 0.9530 0.7453 -0.1836 + 6LYS C 114 1.609 1.980 2.028 -0.0506 0.3156 -0.5774 + 6LYS O 115 1.522 1.915 2.087 -0.0003 0.4952 -0.3006 + 7ASP N 116 1.606 1.984 1.895 0.1681 -0.2273 -0.5995 + 7ASP H 117 1.664 2.055 1.852 1.5543 -1.5568 -1.0282 + 7ASP CA 118 1.509 1.911 1.816 -0.1619 0.2418 -0.6279 + 7ASP HA 119 1.413 1.911 1.868 -0.3314 1.7001 -0.8973 + 7ASP CB 120 1.487 1.969 1.677 0.4599 -0.5002 -1.0444 + 7ASP HB1 121 1.422 1.905 1.618 2.1452 -2.1706 -1.1680 + 7ASP HB2 122 1.582 1.983 1.626 1.0920 1.4900 0.5546 + 7ASP CG 123 1.437 2.112 1.663 0.1593 -0.4950 0.0059 + 7ASP OD1 124 1.373 2.166 1.756 -0.4294 0.3704 -0.8781 + 7ASP OD2 125 1.451 2.169 1.553 0.1888 0.3519 0.4321 + 7ASP C 126 1.545 1.765 1.792 0.3318 0.3579 -0.5897 + 7ASP O 127 1.461 1.686 1.750 -0.0629 0.1795 0.5235 + 8SER N 128 1.670 1.724 1.815 0.2385 0.1590 -0.4441 + 8SER H 129 1.747 1.787 1.832 1.7195 -1.6506 -0.0484 + 8SER CA 130 1.707 1.584 1.799 -0.1400 -0.0362 0.3462 + 8SER HA 131 1.666 1.564 1.700 2.3315 -0.8831 -0.5812 + 8SER CB 132 1.859 1.565 1.798 -0.0427 0.6957 0.2654 + 8SER HB1 133 1.905 1.636 1.730 -0.3410 -0.1665 -0.8636 + 8SER HB2 134 1.879 1.464 1.763 0.7777 0.1943 2.0685 + 8SER OG 135 1.909 1.562 1.930 0.2456 0.6083 0.1532 + 8SER HG 136 1.910 1.652 1.963 -0.0997 1.1774 -1.3071 + 8SER C 137 1.634 1.498 1.900 -0.3761 -0.0841 0.1342 + 8SER O 138 1.618 1.379 1.872 -0.1086 -0.0425 -0.1938 + 9GLY N 139 1.603 1.546 2.021 -0.8083 -0.0860 0.0251 + 9GLY H 140 1.649 1.632 2.046 1.0775 -1.0793 0.1291 + 9GLY CA 141 1.559 1.466 2.134 -0.2229 -0.2196 0.1597 + 9GLY HA1 142 1.472 1.410 2.100 0.9214 0.4124 -4.5132 + 9GLY HA2 143 1.531 1.537 2.212 1.0939 0.5737 -0.0511 + 9GLY C 144 1.665 1.373 2.190 0.1423 0.1523 0.0935 + 9GLY O 145 1.634 1.303 2.287 -0.4717 0.4476 0.1139 + 10LEU N 146 1.786 1.358 2.135 -0.0857 0.0788 -0.3970 + 10LEU H 147 1.810 1.425 2.063 0.9925 -0.8606 -0.9407 + 10LEU CA 148 1.897 1.277 2.180 -0.2014 0.1010 -0.0670 + 10LEU HA 149 1.856 1.189 2.229 -0.3069 -1.5412 -2.9031 + 10LEU CB 150 1.971 1.232 2.055 -0.2947 0.0573 -0.1062 + 10LEU HB1 151 2.052 1.171 2.097 -0.0778 0.7093 0.4427 + 10LEU HB2 152 2.010 1.320 2.004 -0.2353 0.7669 1.1344 + 10LEU CG 153 1.884 1.156 1.956 -0.1920 0.5272 -0.5642 + 10LEU HG 154 1.806 1.212 1.904 -0.3825 -0.6450 -1.5878 + 10LEU CD1 155 1.978 1.101 1.849 -0.0609 0.4181 -0.3931 + 10LEU HD11 156 1.924 1.056 1.766 0.0468 -1.4596 0.5037 + 10LEU HD12 157 2.043 1.175 1.802 1.6876 -0.3308 0.7317 + 10LEU HD13 158 2.042 1.023 1.891 1.0543 1.5004 -0.0242 + 10LEU CD2 159 1.808 1.046 2.029 0.7995 0.0120 -0.2899 + 10LEU HD21 160 1.860 0.977 2.096 -1.4402 -1.0796 0.4376 + 10LEU HD22 161 1.719 1.077 2.085 -1.6181 -2.6169 -2.3914 + 10LEU HD23 162 1.761 0.986 1.951 -0.2814 -0.8529 0.9735 + 10LEU C 163 1.986 1.354 2.278 -0.1641 -0.1064 0.0627 + 10LEU O 164 2.078 1.295 2.334 -0.1827 -0.0007 0.2057 + 11PHE N 165 1.966 1.484 2.299 -0.1651 -0.0363 -0.3616 + 11PHE H 166 1.881 1.529 2.267 -1.5975 -1.4086 1.3043 + 11PHE CA 167 2.063 1.574 2.359 -0.4167 -0.1086 0.1618 + 11PHE HA 168 2.146 1.515 2.399 -0.9779 -0.5706 0.6735 + 11PHE CB 169 2.128 1.654 2.247 -0.7180 -0.2172 -0.0927 + 11PHE HB1 170 2.145 1.758 2.274 -0.0255 -0.3036 -0.1819 + 11PHE HB2 171 2.053 1.660 2.168 0.1930 0.2958 -0.9555 + 11PHE CG 172 2.257 1.600 2.189 -0.3828 0.5005 -0.0355 + 11PHE CD1 173 2.257 1.503 2.087 0.6644 0.5474 -0.0918 + 11PHE HD1 174 2.159 1.478 2.048 -1.1761 4.0170 1.7102 + 11PHE CE1 175 2.377 1.453 2.036 0.3452 -0.9590 0.5803 + 11PHE HE1 176 2.376 1.385 1.951 0.5811 -0.1294 -0.1025 + 11PHE CZ 177 2.497 1.479 2.103 0.4084 0.3115 0.0073 + 11PHE HZ 178 2.589 1.435 2.067 -0.2225 -1.0535 -0.0223 + 11PHE CE2 179 2.495 1.585 2.194 0.0892 -0.4109 0.8548 + 11PHE HE2 180 2.591 1.614 2.233 -0.9734 0.2875 3.1320 + 11PHE CD2 181 2.378 1.647 2.240 -0.1009 0.1958 -0.4170 + 11PHE HD2 182 2.371 1.726 2.313 -0.8800 -0.3576 0.1153 + 11PHE C 183 2.012 1.664 2.470 0.1919 0.0750 0.2944 + 11PHE OC1 184 1.903 1.721 2.453 0.1210 -0.0696 0.2605 + 11PHE OC2 185 2.078 1.675 2.576 -0.3707 0.4761 0.6025 + 12SOL OW 186 3.873 0.044 3.044 -0.2983 0.3579 -0.0212 + 12SOL HW1 187 3.929 0.094 2.984 2.4308 -0.0923 1.9919 + 12SOL HW2 188 3.882 0.088 3.128 0.1698 -2.0064 1.2441 + 13SOL OW 189 3.300 3.670 3.302 -0.2897 -0.5158 -0.1517 + 13SOL HW1 190 3.229 3.716 3.346 -0.9753 -1.2956 -0.4340 + 13SOL HW2 191 3.334 3.609 3.368 0.5166 0.2449 0.1700 + 14SOL OW 192 3.835 3.670 1.359 0.5733 -0.4605 -0.0371 + 14SOL HW1 193 3.783 3.614 1.302 1.2593 -2.0638 0.8277 + 14SOL HW2 194 3.915 3.620 1.376 1.1344 0.7012 0.9362 + 15SOL OW 195 0.782 2.741 1.583 0.3241 0.3358 0.3704 + 15SOL HW1 196 0.759 2.784 1.666 -1.5933 -0.1989 0.1555 + 15SOL HW2 197 0.713 2.677 1.569 -0.0345 1.1794 -2.2267 + 16SOL OW 198 1.556 1.656 0.760 0.0766 0.1757 -0.0926 + 16SOL HW1 199 1.602 1.738 0.744 1.4009 -0.7885 -1.6586 + 16SOL HW2 200 1.626 1.592 0.772 -0.8985 -0.9048 -0.0060 + 17SOL OW 201 3.141 3.132 3.392 0.0944 0.2707 0.1807 + 17SOL HW1 202 3.169 3.112 3.481 1.0627 -0.8159 -0.3515 + 17SOL HW2 203 3.073 3.068 3.373 -0.1843 0.5675 0.1572 + 18SOL OW 204 0.105 1.434 0.498 -0.4822 -0.4474 -0.5695 + 18SOL HW1 205 0.086 1.354 0.546 0.5427 0.0586 0.7387 + 18SOL HW2 206 0.124 1.405 0.409 -0.4729 -1.6162 -0.2027 + 19SOL OW 207 2.107 2.873 3.366 -0.2955 0.3323 -0.4196 + 19SOL HW1 208 2.126 2.797 3.421 1.6975 0.3369 -1.0365 + 19SOL HW2 209 2.096 2.945 3.429 -0.8651 -0.1900 0.0982 + 20SOL OW 210 0.224 1.018 0.877 0.0313 0.0274 -0.5438 + 20SOL HW1 211 0.137 0.989 0.853 -0.8808 2.0702 0.0306 + 20SOL HW2 212 0.231 0.998 0.971 0.3370 -0.4605 -0.6666 + 21SOL OW 213 2.917 1.502 3.799 0.1262 0.5285 0.2875 + 21SOL HW1 214 2.887 1.585 3.762 -0.1020 0.2306 -0.2081 + 21SOL HW2 215 2.963 1.458 3.728 -1.6917 -0.7588 -0.1707 + 3.98279 3.98279 3.98279 diff --git a/src/testutils/simulationdatabase/orires_1lvz.ndx b/src/testutils/simulationdatabase/orires_1lvz.ndx new file mode 100644 index 0000000000..590c396d67 --- /dev/null +++ b/src/testutils/simulationdatabase/orires_1lvz.ndx @@ -0,0 +1,23 @@ +[ System ] + 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 + 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 + 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 + 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 + 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 + 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 + 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 + 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 + 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 + 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 + 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 + 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 + 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 + 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 + 211 212 213 214 215 +[ Backbone ] + 1 5 20 22 24 44 46 48 59 61 63 73 75 77 92 + 94 96 114 116 118 126 128 130 137 139 141 144 146 148 163 + 165 167 183 +[ SOL ] + 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 + 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 diff --git a/src/testutils/simulationdatabase/orires_1lvz.top b/src/testutils/simulationdatabase/orires_1lvz.top new file mode 100644 index 0000000000..0748a4399f --- /dev/null +++ b/src/testutils/simulationdatabase/orires_1lvz.top @@ -0,0 +1,1820 @@ +; +; File '1lvz.top' was generated +; By user: anna (502) +; On host: phy-freki.physics.uu.se +; At date: Mon Feb 3 17:42:15 2020 +; +; This is a standalone topology file +; +; Created by: +; :-) GROMACS - gmx pdb2gmx, 2019.4 (-: +; +; Executable: /usr/local/gromacs/bin/gmx +; Data prefix: /usr/local/gromacs +; Working dir: /Users/anna/Git/NMRreader +; Command line: +; gmx -quiet pdb2gmx -f 1lvz.pdb -ignh -ff amber99sb-ildn -water tip3p -p 1lvz.top -o 1lvz.gro +; Force field was read from the standard GROMACS share directory. +; + +; Include forcefield parameters +#include "amber99sb-ildn.ff/forcefield.itp" + +[ moleculetype ] +; Name nrexcl +Protein_chain_A 3 + +[ atoms ] +; nr type resnr residue atom cgnr charge mass typeB chargeB massB +; residue 1 ILE rtp NILE q +1.0 + 1 N3 1 ILE N 1 0.0311 14.01 + 2 H 1 ILE H1 2 0.2329 1.008 + 3 H 1 ILE H2 3 0.2329 1.008 + 4 H 1 ILE H3 4 0.2329 1.008 + 5 CT 1 ILE CA 5 0.0257 12.01 + 6 HP 1 ILE HA 6 0.1031 1.008 + 7 CT 1 ILE CB 7 0.1885 12.01 + 8 HC 1 ILE HB 8 0.0213 1.008 + 9 CT 1 ILE CG2 9 -0.372 12.01 + 10 HC 1 ILE HG21 10 0.0947 1.008 + 11 HC 1 ILE HG22 11 0.0947 1.008 + 12 HC 1 ILE HG23 12 0.0947 1.008 + 13 CT 1 ILE CG1 13 -0.0387 12.01 + 14 HC 1 ILE HG11 14 0.0201 1.008 + 15 HC 1 ILE HG12 15 0.0201 1.008 + 16 CT 1 ILE CD 16 -0.0908 12.01 + 17 HC 1 ILE HD1 17 0.0226 1.008 + 18 HC 1 ILE HD2 18 0.0226 1.008 + 19 HC 1 ILE HD3 19 0.0226 1.008 + 20 C 1 ILE C 20 0.6123 12.01 + 21 O 1 ILE O 21 -0.5713 16 ; qtot 1 +; residue 2 ARG rtp ARG q +1.0 + 22 N 2 ARG N 22 -0.3479 14.01 + 23 H 2 ARG H 23 0.2747 1.008 + 24 CT 2 ARG CA 24 -0.2637 12.01 + 25 H1 2 ARG HA 25 0.156 1.008 + 26 CT 2 ARG CB 26 -0.0007 12.01 + 27 HC 2 ARG HB1 27 0.0327 1.008 + 28 HC 2 ARG HB2 28 0.0327 1.008 + 29 CT 2 ARG CG 29 0.039 12.01 + 30 HC 2 ARG HG1 30 0.0285 1.008 + 31 HC 2 ARG HG2 31 0.0285 1.008 + 32 CT 2 ARG CD 32 0.0486 12.01 + 33 H1 2 ARG HD1 33 0.0687 1.008 + 34 H1 2 ARG HD2 34 0.0687 1.008 + 35 N2 2 ARG NE 35 -0.5295 14.01 + 36 H 2 ARG HE 36 0.3456 1.008 + 37 CA 2 ARG CZ 37 0.8076 12.01 + 38 N2 2 ARG NH1 38 -0.8627 14.01 + 39 H 2 ARG HH11 39 0.4478 1.008 + 40 H 2 ARG HH12 40 0.4478 1.008 + 41 N2 2 ARG NH2 41 -0.8627 14.01 + 42 H 2 ARG HH21 42 0.4478 1.008 + 43 H 2 ARG HH22 43 0.4478 1.008 + 44 C 2 ARG C 44 0.7341 12.01 + 45 O 2 ARG O 45 -0.5894 16 ; qtot 2 +; residue 3 GLU rtp GLU q -1.0 + 46 N 3 GLU N 46 -0.5163 14.01 + 47 H 3 GLU H 47 0.2936 1.008 + 48 CT 3 GLU CA 48 0.0397 12.01 + 49 H1 3 GLU HA 49 0.1105 1.008 + 50 CT 3 GLU CB 50 0.056 12.01 + 51 HC 3 GLU HB1 51 -0.0173 1.008 + 52 HC 3 GLU HB2 52 -0.0173 1.008 + 53 CT 3 GLU CG 53 0.0136 12.01 + 54 HC 3 GLU HG1 54 -0.0425 1.008 + 55 HC 3 GLU HG2 55 -0.0425 1.008 + 56 C 3 GLU CD 56 0.8054 12.01 + 57 O2 3 GLU OE1 57 -0.8188 16 + 58 O2 3 GLU OE2 58 -0.8188 16 + 59 C 3 GLU C 59 0.5366 12.01 + 60 O 3 GLU O 60 -0.5819 16 ; qtot 1 +; residue 4 ASN rtp ASN q 0.0 + 61 N 4 ASN N 61 -0.4157 14.01 + 62 H 4 ASN H 62 0.2719 1.008 + 63 CT 4 ASN CA 63 0.0143 12.01 + 64 H1 4 ASN HA 64 0.1048 1.008 + 65 CT 4 ASN CB 65 -0.2041 12.01 + 66 HC 4 ASN HB1 66 0.0797 1.008 + 67 HC 4 ASN HB2 67 0.0797 1.008 + 68 C 4 ASN CG 68 0.713 12.01 + 69 O 4 ASN OD1 69 -0.5931 16 + 70 N 4 ASN ND2 70 -0.9191 14.01 + 71 H 4 ASN HD21 71 0.4196 1.008 + 72 H 4 ASN HD22 72 0.4196 1.008 + 73 C 4 ASN C 73 0.5973 12.01 + 74 O 4 ASN O 74 -0.5679 16 ; qtot 1 +; residue 5 LEU rtp LEU q 0.0 + 75 N 5 LEU N 75 -0.4157 14.01 + 76 H 5 LEU H 76 0.2719 1.008 + 77 CT 5 LEU CA 77 -0.0518 12.01 + 78 H1 5 LEU HA 78 0.0922 1.008 + 79 CT 5 LEU CB 79 -0.1102 12.01 + 80 HC 5 LEU HB1 80 0.0457 1.008 + 81 HC 5 LEU HB2 81 0.0457 1.008 + 82 CT 5 LEU CG 82 0.3531 12.01 + 83 HC 5 LEU HG 83 -0.0361 1.008 + 84 CT 5 LEU CD1 84 -0.4121 12.01 + 85 HC 5 LEU HD11 85 0.1 1.008 + 86 HC 5 LEU HD12 86 0.1 1.008 + 87 HC 5 LEU HD13 87 0.1 1.008 + 88 CT 5 LEU CD2 88 -0.4121 12.01 + 89 HC 5 LEU HD21 89 0.1 1.008 + 90 HC 5 LEU HD22 90 0.1 1.008 + 91 HC 5 LEU HD23 91 0.1 1.008 + 92 C 5 LEU C 92 0.5973 12.01 + 93 O 5 LEU O 93 -0.5679 16 ; qtot 1 +; residue 6 LYS rtp LYS q +1.0 + 94 N 6 LYS N 94 -0.3479 14.01 + 95 H 6 LYS H 95 0.2747 1.008 + 96 CT 6 LYS CA 96 -0.24 12.01 + 97 H1 6 LYS HA 97 0.1426 1.008 + 98 CT 6 LYS CB 98 -0.0094 12.01 + 99 HC 6 LYS HB1 99 0.0362 1.008 + 100 HC 6 LYS HB2 100 0.0362 1.008 + 101 CT 6 LYS CG 101 0.0187 12.01 + 102 HC 6 LYS HG1 102 0.0103 1.008 + 103 HC 6 LYS HG2 103 0.0103 1.008 + 104 CT 6 LYS CD 104 -0.0479 12.01 + 105 HC 6 LYS HD1 105 0.0621 1.008 + 106 HC 6 LYS HD2 106 0.0621 1.008 + 107 CT 6 LYS CE 107 -0.0143 12.01 + 108 HP 6 LYS HE1 108 0.1135 1.008 + 109 HP 6 LYS HE2 109 0.1135 1.008 + 110 N3 6 LYS NZ 110 -0.3854 14.01 + 111 H 6 LYS HZ1 111 0.34 1.008 + 112 H 6 LYS HZ2 112 0.34 1.008 + 113 H 6 LYS HZ3 113 0.34 1.008 + 114 C 6 LYS C 114 0.7341 12.01 + 115 O 6 LYS O 115 -0.5894 16 ; qtot 2 +; residue 7 ASP rtp ASP q -1.0 + 116 N 7 ASP N 116 -0.5163 14.01 + 117 H 7 ASP H 117 0.2936 1.008 + 118 CT 7 ASP CA 118 0.0381 12.01 + 119 H1 7 ASP HA 119 0.088 1.008 + 120 CT 7 ASP CB 120 -0.0303 12.01 + 121 HC 7 ASP HB1 121 -0.0122 1.008 + 122 HC 7 ASP HB2 122 -0.0122 1.008 + 123 C 7 ASP CG 123 0.7994 12.01 + 124 O2 7 ASP OD1 124 -0.8014 16 + 125 O2 7 ASP OD2 125 -0.8014 16 + 126 C 7 ASP C 126 0.5366 12.01 + 127 O 7 ASP O 127 -0.5819 16 ; qtot 1 +; residue 8 SER rtp SER q 0.0 + 128 N 8 SER N 128 -0.4157 14.01 + 129 H 8 SER H 129 0.2719 1.008 + 130 CT 8 SER CA 130 -0.0249 12.01 + 131 H1 8 SER HA 131 0.0843 1.008 + 132 CT 8 SER CB 132 0.2117 12.01 + 133 H1 8 SER HB1 133 0.0352 1.008 + 134 H1 8 SER HB2 134 0.0352 1.008 + 135 OH 8 SER OG 135 -0.6546 16 + 136 HO 8 SER HG 136 0.4275 1.008 + 137 C 8 SER C 137 0.5973 12.01 + 138 O 8 SER O 138 -0.5679 16 ; qtot 1 +; residue 9 GLY rtp GLY q 0.0 + 139 N 9 GLY N 139 -0.4157 14.01 + 140 H 9 GLY H 140 0.2719 1.008 + 141 CT 9 GLY CA 141 -0.0252 12.01 + 142 H1 9 GLY HA1 142 0.0698 1.008 + 143 H1 9 GLY HA2 143 0.0698 1.008 + 144 C 9 GLY C 144 0.5973 12.01 + 145 O 9 GLY O 145 -0.5679 16 ; qtot 1 +; residue 10 LEU rtp LEU q 0.0 + 146 N 10 LEU N 146 -0.4157 14.01 + 147 H 10 LEU H 147 0.2719 1.008 + 148 CT 10 LEU CA 148 -0.0518 12.01 + 149 H1 10 LEU HA 149 0.0922 1.008 + 150 CT 10 LEU CB 150 -0.1102 12.01 + 151 HC 10 LEU HB1 151 0.0457 1.008 + 152 HC 10 LEU HB2 152 0.0457 1.008 + 153 CT 10 LEU CG 153 0.3531 12.01 + 154 HC 10 LEU HG 154 -0.0361 1.008 + 155 CT 10 LEU CD1 155 -0.4121 12.01 + 156 HC 10 LEU HD11 156 0.1 1.008 + 157 HC 10 LEU HD12 157 0.1 1.008 + 158 HC 10 LEU HD13 158 0.1 1.008 + 159 CT 10 LEU CD2 159 -0.4121 12.01 + 160 HC 10 LEU HD21 160 0.1 1.008 + 161 HC 10 LEU HD22 161 0.1 1.008 + 162 HC 10 LEU HD23 162 0.1 1.008 + 163 C 10 LEU C 163 0.5973 12.01 + 164 O 10 LEU O 164 -0.5679 16 ; qtot 1 +; residue 11 PHE rtp CPHE q -1.0 + 165 N 11 PHE N 165 -0.3821 14.01 + 166 H 11 PHE H 166 0.2681 1.008 + 167 CT 11 PHE CA 167 -0.1825 12.01 + 168 H1 11 PHE HA 168 0.1098 1.008 + 169 CT 11 PHE CB 169 -0.0959 12.01 + 170 HC 11 PHE HB1 170 0.0443 1.008 + 171 HC 11 PHE HB2 171 0.0443 1.008 + 172 CA 11 PHE CG 172 0.0552 12.01 + 173 CA 11 PHE CD1 173 -0.13 12.01 + 174 HA 11 PHE HD1 174 0.1408 1.008 + 175 CA 11 PHE CE1 175 -0.1847 12.01 + 176 HA 11 PHE HE1 176 0.1461 1.008 + 177 CA 11 PHE CZ 177 -0.0944 12.01 + 178 HA 11 PHE HZ 178 0.128 1.008 + 179 CA 11 PHE CE2 179 -0.1847 12.01 + 180 HA 11 PHE HE2 180 0.1461 1.008 + 181 CA 11 PHE CD2 181 -0.13 12.01 + 182 HA 11 PHE HD2 182 0.1408 1.008 + 183 C 11 PHE C 183 0.766 12.01 + 184 O2 11 PHE OC1 184 -0.8026 16 + 185 O2 11 PHE OC2 185 -0.8026 16 ; qtot 0 + +[ bonds ] +; ai aj funct c0 c1 c2 c3 + 1 2 1 + 1 3 1 + 1 4 1 + 1 5 1 + 5 6 1 + 5 7 1 + 5 20 1 + 7 8 1 + 7 9 1 + 7 13 1 + 9 10 1 + 9 11 1 + 9 12 1 + 13 14 1 + 13 15 1 + 13 16 1 + 16 17 1 + 16 18 1 + 16 19 1 + 20 21 1 + 20 22 1 + 22 23 1 + 22 24 1 + 24 25 1 + 24 26 1 + 24 44 1 + 26 27 1 + 26 28 1 + 26 29 1 + 29 30 1 + 29 31 1 + 29 32 1 + 32 33 1 + 32 34 1 + 32 35 1 + 35 36 1 + 35 37 1 + 37 38 1 + 37 41 1 + 38 39 1 + 38 40 1 + 41 42 1 + 41 43 1 + 44 45 1 + 44 46 1 + 46 47 1 + 46 48 1 + 48 49 1 + 48 50 1 + 48 59 1 + 50 51 1 + 50 52 1 + 50 53 1 + 53 54 1 + 53 55 1 + 53 56 1 + 56 57 1 + 56 58 1 + 59 60 1 + 59 61 1 + 61 62 1 + 61 63 1 + 63 64 1 + 63 65 1 + 63 73 1 + 65 66 1 + 65 67 1 + 65 68 1 + 68 69 1 + 68 70 1 + 70 71 1 + 70 72 1 + 73 74 1 + 73 75 1 + 75 76 1 + 75 77 1 + 77 78 1 + 77 79 1 + 77 92 1 + 79 80 1 + 79 81 1 + 79 82 1 + 82 83 1 + 82 84 1 + 82 88 1 + 84 85 1 + 84 86 1 + 84 87 1 + 88 89 1 + 88 90 1 + 88 91 1 + 92 93 1 + 92 94 1 + 94 95 1 + 94 96 1 + 96 97 1 + 96 98 1 + 96 114 1 + 98 99 1 + 98 100 1 + 98 101 1 + 101 102 1 + 101 103 1 + 101 104 1 + 104 105 1 + 104 106 1 + 104 107 1 + 107 108 1 + 107 109 1 + 107 110 1 + 110 111 1 + 110 112 1 + 110 113 1 + 114 115 1 + 114 116 1 + 116 117 1 + 116 118 1 + 118 119 1 + 118 120 1 + 118 126 1 + 120 121 1 + 120 122 1 + 120 123 1 + 123 124 1 + 123 125 1 + 126 127 1 + 126 128 1 + 128 129 1 + 128 130 1 + 130 131 1 + 130 132 1 + 130 137 1 + 132 133 1 + 132 134 1 + 132 135 1 + 135 136 1 + 137 138 1 + 137 139 1 + 139 140 1 + 139 141 1 + 141 142 1 + 141 143 1 + 141 144 1 + 144 145 1 + 144 146 1 + 146 147 1 + 146 148 1 + 148 149 1 + 148 150 1 + 148 163 1 + 150 151 1 + 150 152 1 + 150 153 1 + 153 154 1 + 153 155 1 + 153 159 1 + 155 156 1 + 155 157 1 + 155 158 1 + 159 160 1 + 159 161 1 + 159 162 1 + 163 164 1 + 163 165 1 + 165 166 1 + 165 167 1 + 167 168 1 + 167 169 1 + 167 183 1 + 169 170 1 + 169 171 1 + 169 172 1 + 172 173 1 + 172 181 1 + 173 174 1 + 173 175 1 + 175 176 1 + 175 177 1 + 177 178 1 + 177 179 1 + 179 180 1 + 179 181 1 + 181 182 1 + 183 184 1 + 183 185 1 + +[ pairs ] +; ai aj funct c0 c1 c2 c3 + 1 8 1 + 1 9 1 + 1 13 1 + 1 21 1 + 1 22 1 + 2 6 1 + 2 7 1 + 2 20 1 + 3 6 1 + 3 7 1 + 3 20 1 + 4 6 1 + 4 7 1 + 4 20 1 + 5 10 1 + 5 11 1 + 5 12 1 + 5 14 1 + 5 15 1 + 5 16 1 + 5 23 1 + 5 24 1 + 6 8 1 + 6 9 1 + 6 13 1 + 6 21 1 + 6 22 1 + 7 17 1 + 7 18 1 + 7 19 1 + 7 21 1 + 7 22 1 + 8 10 1 + 8 11 1 + 8 12 1 + 8 14 1 + 8 15 1 + 8 16 1 + 8 20 1 + 9 14 1 + 9 15 1 + 9 16 1 + 9 20 1 + 10 13 1 + 11 13 1 + 12 13 1 + 13 20 1 + 14 17 1 + 14 18 1 + 14 19 1 + 15 17 1 + 15 18 1 + 15 19 1 + 20 25 1 + 20 26 1 + 20 44 1 + 21 23 1 + 21 24 1 + 22 27 1 + 22 28 1 + 22 29 1 + 22 45 1 + 22 46 1 + 23 25 1 + 23 26 1 + 23 44 1 + 24 30 1 + 24 31 1 + 24 32 1 + 24 47 1 + 24 48 1 + 25 27 1 + 25 28 1 + 25 29 1 + 25 45 1 + 25 46 1 + 26 33 1 + 26 34 1 + 26 35 1 + 26 45 1 + 26 46 1 + 27 30 1 + 27 31 1 + 27 32 1 + 27 44 1 + 28 30 1 + 28 31 1 + 28 32 1 + 28 44 1 + 29 36 1 + 29 37 1 + 29 44 1 + 30 33 1 + 30 34 1 + 30 35 1 + 31 33 1 + 31 34 1 + 31 35 1 + 32 38 1 + 32 41 1 + 33 36 1 + 33 37 1 + 34 36 1 + 34 37 1 + 35 39 1 + 35 40 1 + 35 42 1 + 35 43 1 + 36 38 1 + 36 41 1 + 38 42 1 + 38 43 1 + 39 41 1 + 40 41 1 + 44 49 1 + 44 50 1 + 44 59 1 + 45 47 1 + 45 48 1 + 46 51 1 + 46 52 1 + 46 53 1 + 46 60 1 + 46 61 1 + 47 49 1 + 47 50 1 + 47 59 1 + 48 54 1 + 48 55 1 + 48 56 1 + 48 62 1 + 48 63 1 + 49 51 1 + 49 52 1 + 49 53 1 + 49 60 1 + 49 61 1 + 50 57 1 + 50 58 1 + 50 60 1 + 50 61 1 + 51 54 1 + 51 55 1 + 51 56 1 + 51 59 1 + 52 54 1 + 52 55 1 + 52 56 1 + 52 59 1 + 53 59 1 + 54 57 1 + 54 58 1 + 55 57 1 + 55 58 1 + 59 64 1 + 59 65 1 + 59 73 1 + 60 62 1 + 60 63 1 + 61 66 1 + 61 67 1 + 61 68 1 + 61 74 1 + 61 75 1 + 62 64 1 + 62 65 1 + 62 73 1 + 63 69 1 + 63 70 1 + 63 76 1 + 63 77 1 + 64 66 1 + 64 67 1 + 64 68 1 + 64 74 1 + 64 75 1 + 65 71 1 + 65 72 1 + 65 74 1 + 65 75 1 + 66 69 1 + 66 70 1 + 66 73 1 + 67 69 1 + 67 70 1 + 67 73 1 + 68 73 1 + 69 71 1 + 69 72 1 + 73 78 1 + 73 79 1 + 73 92 1 + 74 76 1 + 74 77 1 + 75 80 1 + 75 81 1 + 75 82 1 + 75 93 1 + 75 94 1 + 76 78 1 + 76 79 1 + 76 92 1 + 77 83 1 + 77 84 1 + 77 88 1 + 77 95 1 + 77 96 1 + 78 80 1 + 78 81 1 + 78 82 1 + 78 93 1 + 78 94 1 + 79 85 1 + 79 86 1 + 79 87 1 + 79 89 1 + 79 90 1 + 79 91 1 + 79 93 1 + 79 94 1 + 80 83 1 + 80 84 1 + 80 88 1 + 80 92 1 + 81 83 1 + 81 84 1 + 81 88 1 + 81 92 1 + 82 92 1 + 83 85 1 + 83 86 1 + 83 87 1 + 83 89 1 + 83 90 1 + 83 91 1 + 84 89 1 + 84 90 1 + 84 91 1 + 85 88 1 + 86 88 1 + 87 88 1 + 92 97 1 + 92 98 1 + 92 114 1 + 93 95 1 + 93 96 1 + 94 99 1 + 94 100 1 + 94 101 1 + 94 115 1 + 94 116 1 + 95 97 1 + 95 98 1 + 95 114 1 + 96 102 1 + 96 103 1 + 96 104 1 + 96 117 1 + 96 118 1 + 97 99 1 + 97 100 1 + 97 101 1 + 97 115 1 + 97 116 1 + 98 105 1 + 98 106 1 + 98 107 1 + 98 115 1 + 98 116 1 + 99 102 1 + 99 103 1 + 99 104 1 + 99 114 1 + 100 102 1 + 100 103 1 + 100 104 1 + 100 114 1 + 101 108 1 + 101 109 1 + 101 110 1 + 101 114 1 + 102 105 1 + 102 106 1 + 102 107 1 + 103 105 1 + 103 106 1 + 103 107 1 + 104 111 1 + 104 112 1 + 104 113 1 + 105 108 1 + 105 109 1 + 105 110 1 + 106 108 1 + 106 109 1 + 106 110 1 + 108 111 1 + 108 112 1 + 108 113 1 + 109 111 1 + 109 112 1 + 109 113 1 + 114 119 1 + 114 120 1 + 114 126 1 + 115 117 1 + 115 118 1 + 116 121 1 + 116 122 1 + 116 123 1 + 116 127 1 + 116 128 1 + 117 119 1 + 117 120 1 + 117 126 1 + 118 124 1 + 118 125 1 + 118 129 1 + 118 130 1 + 119 121 1 + 119 122 1 + 119 123 1 + 119 127 1 + 119 128 1 + 120 127 1 + 120 128 1 + 121 124 1 + 121 125 1 + 121 126 1 + 122 124 1 + 122 125 1 + 122 126 1 + 123 126 1 + 126 131 1 + 126 132 1 + 126 137 1 + 127 129 1 + 127 130 1 + 128 133 1 + 128 134 1 + 128 135 1 + 128 138 1 + 128 139 1 + 129 131 1 + 129 132 1 + 129 137 1 + 130 136 1 + 130 140 1 + 130 141 1 + 131 133 1 + 131 134 1 + 131 135 1 + 131 138 1 + 131 139 1 + 132 138 1 + 132 139 1 + 133 136 1 + 133 137 1 + 134 136 1 + 134 137 1 + 135 137 1 + 137 142 1 + 137 143 1 + 137 144 1 + 138 140 1 + 138 141 1 + 139 145 1 + 139 146 1 + 140 142 1 + 140 143 1 + 140 144 1 + 141 147 1 + 141 148 1 + 142 145 1 + 142 146 1 + 143 145 1 + 143 146 1 + 144 149 1 + 144 150 1 + 144 163 1 + 145 147 1 + 145 148 1 + 146 151 1 + 146 152 1 + 146 153 1 + 146 164 1 + 146 165 1 + 147 149 1 + 147 150 1 + 147 163 1 + 148 154 1 + 148 155 1 + 148 159 1 + 148 166 1 + 148 167 1 + 149 151 1 + 149 152 1 + 149 153 1 + 149 164 1 + 149 165 1 + 150 156 1 + 150 157 1 + 150 158 1 + 150 160 1 + 150 161 1 + 150 162 1 + 150 164 1 + 150 165 1 + 151 154 1 + 151 155 1 + 151 159 1 + 151 163 1 + 152 154 1 + 152 155 1 + 152 159 1 + 152 163 1 + 153 163 1 + 154 156 1 + 154 157 1 + 154 158 1 + 154 160 1 + 154 161 1 + 154 162 1 + 155 160 1 + 155 161 1 + 155 162 1 + 156 159 1 + 157 159 1 + 158 159 1 + 163 168 1 + 163 169 1 + 163 183 1 + 164 166 1 + 164 167 1 + 165 170 1 + 165 171 1 + 165 172 1 + 165 184 1 + 165 185 1 + 166 168 1 + 166 169 1 + 166 183 1 + 167 173 1 + 167 181 1 + 168 170 1 + 168 171 1 + 168 172 1 + 168 184 1 + 168 185 1 + 169 174 1 + 169 175 1 + 169 179 1 + 169 182 1 + 169 184 1 + 169 185 1 + 170 173 1 + 170 181 1 + 170 183 1 + 171 173 1 + 171 181 1 + 171 183 1 + 172 176 1 + 172 177 1 + 172 180 1 + 172 183 1 + 173 178 1 + 173 179 1 + 173 182 1 + 174 176 1 + 174 177 1 + 174 181 1 + 175 180 1 + 175 181 1 + 176 178 1 + 176 179 1 + 177 182 1 + 178 180 1 + 178 181 1 + 180 182 1 + +[ angles ] +; ai aj ak funct c0 c1 c2 c3 + 2 1 3 1 + 2 1 4 1 + 2 1 5 1 + 3 1 4 1 + 3 1 5 1 + 4 1 5 1 + 1 5 6 1 + 1 5 7 1 + 1 5 20 1 + 6 5 7 1 + 6 5 20 1 + 7 5 20 1 + 5 7 8 1 + 5 7 9 1 + 5 7 13 1 + 8 7 9 1 + 8 7 13 1 + 9 7 13 1 + 7 9 10 1 + 7 9 11 1 + 7 9 12 1 + 10 9 11 1 + 10 9 12 1 + 11 9 12 1 + 7 13 14 1 + 7 13 15 1 + 7 13 16 1 + 14 13 15 1 + 14 13 16 1 + 15 13 16 1 + 13 16 17 1 + 13 16 18 1 + 13 16 19 1 + 17 16 18 1 + 17 16 19 1 + 18 16 19 1 + 5 20 21 1 + 5 20 22 1 + 21 20 22 1 + 20 22 23 1 + 20 22 24 1 + 23 22 24 1 + 22 24 25 1 + 22 24 26 1 + 22 24 44 1 + 25 24 26 1 + 25 24 44 1 + 26 24 44 1 + 24 26 27 1 + 24 26 28 1 + 24 26 29 1 + 27 26 28 1 + 27 26 29 1 + 28 26 29 1 + 26 29 30 1 + 26 29 31 1 + 26 29 32 1 + 30 29 31 1 + 30 29 32 1 + 31 29 32 1 + 29 32 33 1 + 29 32 34 1 + 29 32 35 1 + 33 32 34 1 + 33 32 35 1 + 34 32 35 1 + 32 35 36 1 + 32 35 37 1 + 36 35 37 1 + 35 37 38 1 + 35 37 41 1 + 38 37 41 1 + 37 38 39 1 + 37 38 40 1 + 39 38 40 1 + 37 41 42 1 + 37 41 43 1 + 42 41 43 1 + 24 44 45 1 + 24 44 46 1 + 45 44 46 1 + 44 46 47 1 + 44 46 48 1 + 47 46 48 1 + 46 48 49 1 + 46 48 50 1 + 46 48 59 1 + 49 48 50 1 + 49 48 59 1 + 50 48 59 1 + 48 50 51 1 + 48 50 52 1 + 48 50 53 1 + 51 50 52 1 + 51 50 53 1 + 52 50 53 1 + 50 53 54 1 + 50 53 55 1 + 50 53 56 1 + 54 53 55 1 + 54 53 56 1 + 55 53 56 1 + 53 56 57 1 + 53 56 58 1 + 57 56 58 1 + 48 59 60 1 + 48 59 61 1 + 60 59 61 1 + 59 61 62 1 + 59 61 63 1 + 62 61 63 1 + 61 63 64 1 + 61 63 65 1 + 61 63 73 1 + 64 63 65 1 + 64 63 73 1 + 65 63 73 1 + 63 65 66 1 + 63 65 67 1 + 63 65 68 1 + 66 65 67 1 + 66 65 68 1 + 67 65 68 1 + 65 68 69 1 + 65 68 70 1 + 69 68 70 1 + 68 70 71 1 + 68 70 72 1 + 71 70 72 1 + 63 73 74 1 + 63 73 75 1 + 74 73 75 1 + 73 75 76 1 + 73 75 77 1 + 76 75 77 1 + 75 77 78 1 + 75 77 79 1 + 75 77 92 1 + 78 77 79 1 + 78 77 92 1 + 79 77 92 1 + 77 79 80 1 + 77 79 81 1 + 77 79 82 1 + 80 79 81 1 + 80 79 82 1 + 81 79 82 1 + 79 82 83 1 + 79 82 84 1 + 79 82 88 1 + 83 82 84 1 + 83 82 88 1 + 84 82 88 1 + 82 84 85 1 + 82 84 86 1 + 82 84 87 1 + 85 84 86 1 + 85 84 87 1 + 86 84 87 1 + 82 88 89 1 + 82 88 90 1 + 82 88 91 1 + 89 88 90 1 + 89 88 91 1 + 90 88 91 1 + 77 92 93 1 + 77 92 94 1 + 93 92 94 1 + 92 94 95 1 + 92 94 96 1 + 95 94 96 1 + 94 96 97 1 + 94 96 98 1 + 94 96 114 1 + 97 96 98 1 + 97 96 114 1 + 98 96 114 1 + 96 98 99 1 + 96 98 100 1 + 96 98 101 1 + 99 98 100 1 + 99 98 101 1 + 100 98 101 1 + 98 101 102 1 + 98 101 103 1 + 98 101 104 1 + 102 101 103 1 + 102 101 104 1 + 103 101 104 1 + 101 104 105 1 + 101 104 106 1 + 101 104 107 1 + 105 104 106 1 + 105 104 107 1 + 106 104 107 1 + 104 107 108 1 + 104 107 109 1 + 104 107 110 1 + 108 107 109 1 + 108 107 110 1 + 109 107 110 1 + 107 110 111 1 + 107 110 112 1 + 107 110 113 1 + 111 110 112 1 + 111 110 113 1 + 112 110 113 1 + 96 114 115 1 + 96 114 116 1 + 115 114 116 1 + 114 116 117 1 + 114 116 118 1 + 117 116 118 1 + 116 118 119 1 + 116 118 120 1 + 116 118 126 1 + 119 118 120 1 + 119 118 126 1 + 120 118 126 1 + 118 120 121 1 + 118 120 122 1 + 118 120 123 1 + 121 120 122 1 + 121 120 123 1 + 122 120 123 1 + 120 123 124 1 + 120 123 125 1 + 124 123 125 1 + 118 126 127 1 + 118 126 128 1 + 127 126 128 1 + 126 128 129 1 + 126 128 130 1 + 129 128 130 1 + 128 130 131 1 + 128 130 132 1 + 128 130 137 1 + 131 130 132 1 + 131 130 137 1 + 132 130 137 1 + 130 132 133 1 + 130 132 134 1 + 130 132 135 1 + 133 132 134 1 + 133 132 135 1 + 134 132 135 1 + 132 135 136 1 + 130 137 138 1 + 130 137 139 1 + 138 137 139 1 + 137 139 140 1 + 137 139 141 1 + 140 139 141 1 + 139 141 142 1 + 139 141 143 1 + 139 141 144 1 + 142 141 143 1 + 142 141 144 1 + 143 141 144 1 + 141 144 145 1 + 141 144 146 1 + 145 144 146 1 + 144 146 147 1 + 144 146 148 1 + 147 146 148 1 + 146 148 149 1 + 146 148 150 1 + 146 148 163 1 + 149 148 150 1 + 149 148 163 1 + 150 148 163 1 + 148 150 151 1 + 148 150 152 1 + 148 150 153 1 + 151 150 152 1 + 151 150 153 1 + 152 150 153 1 + 150 153 154 1 + 150 153 155 1 + 150 153 159 1 + 154 153 155 1 + 154 153 159 1 + 155 153 159 1 + 153 155 156 1 + 153 155 157 1 + 153 155 158 1 + 156 155 157 1 + 156 155 158 1 + 157 155 158 1 + 153 159 160 1 + 153 159 161 1 + 153 159 162 1 + 160 159 161 1 + 160 159 162 1 + 161 159 162 1 + 148 163 164 1 + 148 163 165 1 + 164 163 165 1 + 163 165 166 1 + 163 165 167 1 + 166 165 167 1 + 165 167 168 1 + 165 167 169 1 + 165 167 183 1 + 168 167 169 1 + 168 167 183 1 + 169 167 183 1 + 167 169 170 1 + 167 169 171 1 + 167 169 172 1 + 170 169 171 1 + 170 169 172 1 + 171 169 172 1 + 169 172 173 1 + 169 172 181 1 + 173 172 181 1 + 172 173 174 1 + 172 173 175 1 + 174 173 175 1 + 173 175 176 1 + 173 175 177 1 + 176 175 177 1 + 175 177 178 1 + 175 177 179 1 + 178 177 179 1 + 177 179 180 1 + 177 179 181 1 + 180 179 181 1 + 172 181 179 1 + 172 181 182 1 + 179 181 182 1 + 167 183 184 1 + 167 183 185 1 + 184 183 185 1 + +[ dihedrals ] +; ai aj ak al funct c0 c1 c2 c3 c4 c5 + 2 1 5 6 9 + 2 1 5 7 9 + 2 1 5 20 9 + 3 1 5 6 9 + 3 1 5 7 9 + 3 1 5 20 9 + 4 1 5 6 9 + 4 1 5 7 9 + 4 1 5 20 9 + 1 5 7 9 9 torsion_ILE_N_CA_CB_CG2_mult1 + 1 5 7 9 9 torsion_ILE_N_CA_CB_CG2_mult2 + 1 5 7 8 9 + 1 5 7 13 9 + 6 5 7 8 9 + 6 5 7 9 9 + 6 5 7 13 9 + 20 5 7 8 9 + 20 5 7 9 9 + 20 5 7 13 9 + 1 5 20 21 9 + 1 5 20 22 9 + 6 5 20 21 9 + 6 5 20 22 9 + 7 5 20 21 9 + 7 5 20 22 9 + 5 7 9 10 9 + 5 7 9 11 9 + 5 7 9 12 9 + 8 7 9 10 9 + 8 7 9 11 9 + 8 7 9 12 9 + 13 7 9 10 9 + 13 7 9 11 9 + 13 7 9 12 9 + 5 7 13 14 9 + 5 7 13 15 9 + 5 7 13 16 9 + 8 7 13 14 9 + 8 7 13 15 9 + 8 7 13 16 9 + 9 7 13 14 9 + 9 7 13 15 9 + 9 7 13 16 9 + 7 13 16 17 9 + 7 13 16 18 9 + 7 13 16 19 9 + 14 13 16 17 9 + 14 13 16 18 9 + 14 13 16 19 9 + 15 13 16 17 9 + 15 13 16 18 9 + 15 13 16 19 9 + 5 20 22 23 9 + 5 20 22 24 9 + 21 20 22 23 9 + 21 20 22 24 9 + 20 22 24 25 9 + 20 22 24 26 9 + 20 22 24 44 9 + 23 22 24 25 9 + 23 22 24 26 9 + 23 22 24 44 9 + 22 24 26 27 9 + 22 24 26 28 9 + 22 24 26 29 9 + 25 24 26 27 9 + 25 24 26 28 9 + 25 24 26 29 9 + 44 24 26 27 9 + 44 24 26 28 9 + 44 24 26 29 9 + 22 24 44 45 9 + 22 24 44 46 9 + 25 24 44 45 9 + 25 24 44 46 9 + 26 24 44 45 9 + 26 24 44 46 9 + 24 26 29 30 9 + 24 26 29 31 9 + 24 26 29 32 9 + 27 26 29 30 9 + 27 26 29 31 9 + 27 26 29 32 9 + 28 26 29 30 9 + 28 26 29 31 9 + 28 26 29 32 9 + 26 29 32 33 9 + 26 29 32 34 9 + 26 29 32 35 9 + 30 29 32 33 9 + 30 29 32 34 9 + 30 29 32 35 9 + 31 29 32 33 9 + 31 29 32 34 9 + 31 29 32 35 9 + 29 32 35 36 9 + 29 32 35 37 9 + 33 32 35 36 9 + 33 32 35 37 9 + 34 32 35 36 9 + 34 32 35 37 9 + 32 35 37 38 9 + 32 35 37 41 9 + 36 35 37 38 9 + 36 35 37 41 9 + 35 37 38 39 9 + 35 37 38 40 9 + 41 37 38 39 9 + 41 37 38 40 9 + 35 37 41 42 9 + 35 37 41 43 9 + 38 37 41 42 9 + 38 37 41 43 9 + 24 44 46 47 9 + 24 44 46 48 9 + 45 44 46 47 9 + 45 44 46 48 9 + 44 46 48 49 9 + 44 46 48 50 9 + 44 46 48 59 9 + 47 46 48 49 9 + 47 46 48 50 9 + 47 46 48 59 9 + 46 48 50 51 9 + 46 48 50 52 9 + 46 48 50 53 9 + 49 48 50 51 9 + 49 48 50 52 9 + 49 48 50 53 9 + 59 48 50 51 9 + 59 48 50 52 9 + 59 48 50 53 9 + 46 48 59 60 9 + 46 48 59 61 9 + 49 48 59 60 9 + 49 48 59 61 9 + 50 48 59 60 9 + 50 48 59 61 9 + 48 50 53 54 9 + 48 50 53 55 9 + 48 50 53 56 9 + 51 50 53 54 9 + 51 50 53 55 9 + 51 50 53 56 9 + 52 50 53 54 9 + 52 50 53 55 9 + 52 50 53 56 9 + 50 53 56 57 9 + 50 53 56 58 9 + 54 53 56 57 9 + 54 53 56 58 9 + 55 53 56 57 9 + 55 53 56 58 9 + 48 59 61 62 9 + 48 59 61 63 9 + 60 59 61 62 9 + 60 59 61 63 9 + 59 61 63 64 9 + 59 61 63 65 9 + 59 61 63 73 9 + 62 61 63 64 9 + 62 61 63 65 9 + 62 61 63 73 9 + 73 63 65 68 9 torsion_ASN_C_CA_CB_CG_mult1 + 73 63 65 68 9 torsion_ASN_C_CA_CB_CG_mult2 + 73 63 65 68 9 torsion_ASN_C_CA_CB_CG_mult3 + 73 63 65 68 9 torsion_ASN_C_CA_CB_CG_mult4 + 73 63 65 68 9 torsion_ASN_C_CA_CB_CG_mult5 + 73 63 65 68 9 torsion_ASN_C_CA_CB_CG_mult6 + 61 63 65 66 9 + 61 63 65 67 9 + 61 63 65 68 9 + 64 63 65 66 9 + 64 63 65 67 9 + 64 63 65 68 9 + 73 63 65 66 9 + 73 63 65 67 9 + 61 63 73 74 9 + 61 63 73 75 9 + 64 63 73 74 9 + 64 63 73 75 9 + 65 63 73 74 9 + 65 63 73 75 9 + 63 65 68 70 9 torsion_ASN_CA_CB_CG_ND2_mult1 + 63 65 68 70 9 torsion_ASN_CA_CB_CG_ND2_mult2 + 63 65 68 70 9 torsion_ASN_CA_CB_CG_ND2_mult3 + 63 65 68 70 9 torsion_ASN_CA_CB_CG_ND2_mult4 + 63 65 68 70 9 torsion_ASN_CA_CB_CG_ND2_mult5 + 63 65 68 70 9 torsion_ASN_CA_CB_CG_ND2_mult6 + 63 65 68 69 9 + 66 65 68 69 9 + 66 65 68 70 9 + 67 65 68 69 9 + 67 65 68 70 9 + 65 68 70 71 9 + 65 68 70 72 9 + 69 68 70 71 9 + 69 68 70 72 9 + 63 73 75 76 9 + 63 73 75 77 9 + 74 73 75 76 9 + 74 73 75 77 9 + 73 75 77 78 9 + 73 75 77 79 9 + 73 75 77 92 9 + 76 75 77 78 9 + 76 75 77 79 9 + 76 75 77 92 9 + 92 77 79 82 9 torsion_LEU_C_CA_CB_CG_mult1 + 92 77 79 82 9 torsion_LEU_C_CA_CB_CG_mult2 + 92 77 79 82 9 torsion_LEU_C_CA_CB_CG_mult3 + 75 77 79 80 9 + 75 77 79 81 9 + 75 77 79 82 9 + 78 77 79 80 9 + 78 77 79 81 9 + 78 77 79 82 9 + 92 77 79 80 9 + 92 77 79 81 9 + 75 77 92 93 9 + 75 77 92 94 9 + 78 77 92 93 9 + 78 77 92 94 9 + 79 77 92 93 9 + 79 77 92 94 9 + 77 79 82 83 9 + 77 79 82 84 9 + 77 79 82 88 9 + 80 79 82 83 9 + 80 79 82 84 9 + 80 79 82 88 9 + 81 79 82 83 9 + 81 79 82 84 9 + 81 79 82 88 9 + 79 82 84 85 9 + 79 82 84 86 9 + 79 82 84 87 9 + 83 82 84 85 9 + 83 82 84 86 9 + 83 82 84 87 9 + 88 82 84 85 9 + 88 82 84 86 9 + 88 82 84 87 9 + 79 82 88 89 9 + 79 82 88 90 9 + 79 82 88 91 9 + 83 82 88 89 9 + 83 82 88 90 9 + 83 82 88 91 9 + 84 82 88 89 9 + 84 82 88 90 9 + 84 82 88 91 9 + 77 92 94 95 9 + 77 92 94 96 9 + 93 92 94 95 9 + 93 92 94 96 9 + 92 94 96 97 9 + 92 94 96 98 9 + 92 94 96 114 9 + 95 94 96 97 9 + 95 94 96 98 9 + 95 94 96 114 9 + 94 96 98 99 9 + 94 96 98 100 9 + 94 96 98 101 9 + 97 96 98 99 9 + 97 96 98 100 9 + 97 96 98 101 9 + 114 96 98 99 9 + 114 96 98 100 9 + 114 96 98 101 9 + 94 96 114 115 9 + 94 96 114 116 9 + 97 96 114 115 9 + 97 96 114 116 9 + 98 96 114 115 9 + 98 96 114 116 9 + 96 98 101 102 9 + 96 98 101 103 9 + 96 98 101 104 9 + 99 98 101 102 9 + 99 98 101 103 9 + 99 98 101 104 9 + 100 98 101 102 9 + 100 98 101 103 9 + 100 98 101 104 9 + 98 101 104 105 9 + 98 101 104 106 9 + 98 101 104 107 9 + 102 101 104 105 9 + 102 101 104 106 9 + 102 101 104 107 9 + 103 101 104 105 9 + 103 101 104 106 9 + 103 101 104 107 9 + 101 104 107 108 9 + 101 104 107 109 9 + 101 104 107 110 9 + 105 104 107 108 9 + 105 104 107 109 9 + 105 104 107 110 9 + 106 104 107 108 9 + 106 104 107 109 9 + 106 104 107 110 9 + 104 107 110 111 9 + 104 107 110 112 9 + 104 107 110 113 9 + 108 107 110 111 9 + 108 107 110 112 9 + 108 107 110 113 9 + 109 107 110 111 9 + 109 107 110 112 9 + 109 107 110 113 9 + 96 114 116 117 9 + 96 114 116 118 9 + 115 114 116 117 9 + 115 114 116 118 9 + 114 116 118 119 9 + 114 116 118 120 9 + 114 116 118 126 9 + 117 116 118 119 9 + 117 116 118 120 9 + 117 116 118 126 9 + 116 118 120 123 9 torsion_ASP_N_CA_CB_CG_mult1 + 116 118 120 123 9 torsion_ASP_N_CA_CB_CG_mult2 + 116 118 120 123 9 torsion_ASP_N_CA_CB_CG_mult3 + 116 118 120 123 9 torsion_ASP_N_CA_CB_CG_mult4 + 116 118 120 123 9 torsion_ASP_N_CA_CB_CG_mult5 + 116 118 120 123 9 torsion_ASP_N_CA_CB_CG_mult6 + 116 118 120 121 9 + 116 118 120 122 9 + 119 118 120 121 9 + 119 118 120 122 9 + 119 118 120 123 9 + 126 118 120 121 9 + 126 118 120 122 9 + 126 118 120 123 9 + 116 118 126 127 9 + 116 118 126 128 9 + 119 118 126 127 9 + 119 118 126 128 9 + 120 118 126 127 9 + 120 118 126 128 9 + 118 120 123 124 9 torsion_ASP_CA_CB_CG_OD_mult1 + 118 120 123 124 9 torsion_ASP_CA_CB_CG_OD_mult2 + 118 120 123 124 9 torsion_ASP_CA_CB_CG_OD_mult3 + 118 120 123 124 9 torsion_ASP_CA_CB_CG_OD_mult4 + 118 120 123 124 9 torsion_ASP_CA_CB_CG_OD_mult5 + 118 120 123 124 9 torsion_ASP_CA_CB_CG_OD_mult6 + 118 120 123 125 9 torsion_ASP_CA_CB_CG_OD_mult1 + 118 120 123 125 9 torsion_ASP_CA_CB_CG_OD_mult2 + 118 120 123 125 9 torsion_ASP_CA_CB_CG_OD_mult3 + 118 120 123 125 9 torsion_ASP_CA_CB_CG_OD_mult4 + 118 120 123 125 9 torsion_ASP_CA_CB_CG_OD_mult5 + 118 120 123 125 9 torsion_ASP_CA_CB_CG_OD_mult6 + 121 120 123 124 9 + 121 120 123 125 9 + 122 120 123 124 9 + 122 120 123 125 9 + 118 126 128 129 9 + 118 126 128 130 9 + 127 126 128 129 9 + 127 126 128 130 9 + 126 128 130 131 9 + 126 128 130 132 9 + 126 128 130 137 9 + 129 128 130 131 9 + 129 128 130 132 9 + 129 128 130 137 9 + 128 130 132 133 9 + 128 130 132 134 9 + 128 130 132 135 9 + 131 130 132 133 9 + 131 130 132 134 9 + 131 130 132 135 9 + 137 130 132 133 9 + 137 130 132 134 9 + 137 130 132 135 9 + 128 130 137 138 9 + 128 130 137 139 9 + 131 130 137 138 9 + 131 130 137 139 9 + 132 130 137 138 9 + 132 130 137 139 9 + 130 132 135 136 9 + 133 132 135 136 9 + 134 132 135 136 9 + 130 137 139 140 9 + 130 137 139 141 9 + 138 137 139 140 9 + 138 137 139 141 9 + 137 139 141 142 9 + 137 139 141 143 9 + 137 139 141 144 9 + 140 139 141 142 9 + 140 139 141 143 9 + 140 139 141 144 9 + 139 141 144 145 9 + 139 141 144 146 9 + 142 141 144 145 9 + 142 141 144 146 9 + 143 141 144 145 9 + 143 141 144 146 9 + 141 144 146 147 9 + 141 144 146 148 9 + 145 144 146 147 9 + 145 144 146 148 9 + 144 146 148 149 9 + 144 146 148 150 9 + 144 146 148 163 9 + 147 146 148 149 9 + 147 146 148 150 9 + 147 146 148 163 9 + 163 148 150 153 9 torsion_LEU_C_CA_CB_CG_mult1 + 163 148 150 153 9 torsion_LEU_C_CA_CB_CG_mult2 + 163 148 150 153 9 torsion_LEU_C_CA_CB_CG_mult3 + 146 148 150 151 9 + 146 148 150 152 9 + 146 148 150 153 9 + 149 148 150 151 9 + 149 148 150 152 9 + 149 148 150 153 9 + 163 148 150 151 9 + 163 148 150 152 9 + 146 148 163 164 9 + 146 148 163 165 9 + 149 148 163 164 9 + 149 148 163 165 9 + 150 148 163 164 9 + 150 148 163 165 9 + 148 150 153 154 9 + 148 150 153 155 9 + 148 150 153 159 9 + 151 150 153 154 9 + 151 150 153 155 9 + 151 150 153 159 9 + 152 150 153 154 9 + 152 150 153 155 9 + 152 150 153 159 9 + 150 153 155 156 9 + 150 153 155 157 9 + 150 153 155 158 9 + 154 153 155 156 9 + 154 153 155 157 9 + 154 153 155 158 9 + 159 153 155 156 9 + 159 153 155 157 9 + 159 153 155 158 9 + 150 153 159 160 9 + 150 153 159 161 9 + 150 153 159 162 9 + 154 153 159 160 9 + 154 153 159 161 9 + 154 153 159 162 9 + 155 153 159 160 9 + 155 153 159 161 9 + 155 153 159 162 9 + 148 163 165 166 9 + 148 163 165 167 9 + 164 163 165 166 9 + 164 163 165 167 9 + 163 165 167 168 9 + 163 165 167 169 9 + 163 165 167 183 9 + 166 165 167 168 9 + 166 165 167 169 9 + 166 165 167 183 9 + 165 167 169 170 9 + 165 167 169 171 9 + 165 167 169 172 9 + 168 167 169 170 9 + 168 167 169 171 9 + 168 167 169 172 9 + 183 167 169 170 9 + 183 167 169 171 9 + 183 167 169 172 9 + 165 167 183 184 9 + 165 167 183 185 9 + 168 167 183 184 9 + 168 167 183 185 9 + 169 167 183 184 9 + 169 167 183 185 9 + 167 169 172 173 9 + 167 169 172 181 9 + 170 169 172 173 9 + 170 169 172 181 9 + 171 169 172 173 9 + 171 169 172 181 9 + 169 172 173 174 9 + 169 172 173 175 9 + 181 172 173 174 9 + 181 172 173 175 9 + 169 172 181 179 9 + 169 172 181 182 9 + 173 172 181 179 9 + 173 172 181 182 9 + 172 173 175 176 9 + 172 173 175 177 9 + 174 173 175 176 9 + 174 173 175 177 9 + 173 175 177 178 9 + 173 175 177 179 9 + 176 175 177 178 9 + 176 175 177 179 9 + 175 177 179 180 9 + 175 177 179 181 9 + 178 177 179 180 9 + 178 177 179 181 9 + 177 179 181 172 9 + 177 179 181 182 9 + 180 179 181 172 9 + 180 179 181 182 9 + +[ dihedrals ] +; ai aj ak al funct c0 c1 c2 c3 + 5 22 20 21 4 + 20 24 22 23 4 + 24 46 44 45 4 + 32 37 35 36 4 + 35 38 37 41 4 + 37 39 38 40 4 + 37 42 41 43 4 + 44 48 46 47 4 + 48 61 59 60 4 + 53 57 56 58 4 + 59 63 61 62 4 + 63 75 73 74 4 + 65 70 68 69 4 + 68 71 70 72 4 + 73 77 75 76 4 + 77 94 92 93 4 + 92 96 94 95 4 + 96 116 114 115 4 + 114 118 116 117 4 + 118 128 126 127 4 + 120 124 123 125 4 + 126 130 128 129 4 + 130 139 137 138 4 + 137 141 139 140 4 + 141 146 144 145 4 + 144 148 146 147 4 + 148 165 163 164 4 + 163 167 165 166 4 + 167 184 183 185 4 + 169 172 181 173 4 + 172 175 173 174 4 + 172 179 181 182 4 + 173 177 175 176 4 + 175 179 177 178 4 + 177 181 179 180 4 + +[ orientation_restraints ] +; ai aj type exp. label alpha const. obs. weight +; Hz nm^3 Hz Hz^-2 + + 46 47 1 1 1 3 6.083 1.1 1.0 + 61 62 1 1 2 3 6.083 -0.4 1.0 + 75 76 1 1 3 3 6.083 0.7 1.0 + 94 95 1 1 4 3 6.083 2.0 1.0 + 116 117 1 1 5 3 6.083 -0.1 1.0 + 128 129 1 1 6 3 6.083 -0.1 1.0 + 139 140 1 1 7 3 6.083 2.4 1.0 + 146 147 1 1 8 3 6.083 -0.3 1.0 + 165 166 1 1 9 3 6.083 -0.7 1.0 + +; Include water topology +#include "amber99sb-ildn.ff/tip3p.itp" + +#ifdef POSRES_WATER +; Position restraint for each water oxygen +[ position_restraints ] +; i funct fcx fcy fcz + 1 1 1000 1000 1000 +#endif + +; Include topology for ions +#include "amber99sb-ildn.ff/ions.itp" + +[ system ] +; Name +GUANINE NUCLEOTIDE-BINDING PROTEIN G(T), ALPHA-1; 3 SUBUNIT in water + +[ molecules ] +; Compound #mols +Protein_chain_A 1 +SOL 10 -- 2.11.4.GIT