1
|
Liu H, Song D, Zhang Y, Yang S, Luo R, Chen HF. Extensive tests and evaluation of the CHARMM36IDPSFF force field for intrinsically disordered proteins and folded proteins. Phys Chem Chem Phys 2019; 21:21918-21931. [PMID: 31552948 DOI: 10.1039/c9cp03434j] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Intrinsically disordered proteins (IDPs) have received increasing attention in recent studies due to their structural heterogeneity and critical biological functions. To fully understand the structural properties and determine accurate ensembles of IDPs, molecular dynamics (MD) simulation was widely used to sample diverse conformations and reveal the structural dynamics. However, the classical state-of-the-art force fields perform well for folded proteins while being unsatisfactory for the simulations of disordered proteins reported in many previous studies. Thus, improved force fields were developed to precisely describe both folded proteins and disordered proteins. Preliminary tests show that our newly developed CHARMM36IDPSFF (C36IDPSFF) force field can well reproduce the experimental observables of several disordered proteins, but more tests on different types of proteins are needed to further evaluate the performance of C36IDPSFF. Here, we extensively simulate short peptides, disordered proteins, and fast-folding proteins as well as folded proteins, and compare the simulated results with the experimental observables. The simulation results show that C36IDPSFF could substantially reproduce the experimental observables for most of the tested proteins but some limitations are also found in the radius of gyration of large disordered proteins and the stability of fast-folding proteins. This force field will facilitate large scale studies of protein structural dynamics and functions using MD simulations.
Collapse
Affiliation(s)
- Hao Liu
- State Key Laboratory of Microbial Metabolism, Department of Bioinformatics and Biostatistics, SJTU-Yale Joint Center for Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai, 200240, China.
| | | | | | | | | | | |
Collapse
|
2
|
Chatfield DC, Augsten A, D'Cunha C, Wong SE. Methyl dynamics in crystalline amino acids: MD and NMR. J Comput Chem 2003; 24:1052-8. [PMID: 12759905 DOI: 10.1002/jcc.10263] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Correlation times for rotation of deuterated methyls in crystalline leucine, valine, and cyclo-L-alanyl-L-alanine are calculated with molecular dynamics and compared with NMR data. The simulations distinguish between methyls having different steric environments in the crystal, yielding correlation times differing by a factor of up to 30 for methyls within a given crystal. MD and NMR correlation times agree to within a factor of 2. However, averaging over nonequivalent methyls can yield correlation functions that, although actually multiexponential, are well fit by single exponentials. This may have significance for interpreting NMR data; previous NMR data did not distinguish between the methyls in these crystals. Adiabatic rotational barriers calculated with the X-ray structure differ from effective barriers during simulation by up to +/-1 kcal/mol; the difference indicates that dynamical effects have a significant role in determining rotational correlation times.
Collapse
Affiliation(s)
- David C Chatfield
- Chemistry Department, Florida International University, Miami, Florida 33199, USA.
| | | | | | | |
Collapse
|
3
|
Dölle A. Reorientational Dynamics of the Model Compound 1,2,3,4-Tetrahydro-5,6-dimethyl- 1,4-methanonaphthalene in Neat Liquid from Temperature-Dependent 13C Nuclear Magnetic Relaxation Data: Spectral Densities and Correlation Functions. J Phys Chem A 2002. [DOI: 10.1021/jp0208438] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andreas Dölle
- Institut für Physikalische Chemie, Rheinisch-Westfälische Technische Hochschule, 52056 Aachen, Germany
| |
Collapse
|
4
|
Pfeiffer S, Fushman D, Cowburn D. Simulated and NMR-derived backbone dynamics of a protein with significant flexibility: a comparison of spectral densities for the betaARK1 PH domain. J Am Chem Soc 2001; 123:3021-36. [PMID: 11457013 DOI: 10.1021/ja0031117] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A 7.6 ns molecular dynamics trajectory of the betaARK1 PH domain in explicit water with appropriate ions was calculated at 300 K. Spectral densities at omega = 0, omega(N), and 0.87omega(H) and the model-free parameters were evaluated from the experimental as well as the simulated data, taking the anisotropic overall motion of the protein into account. Experimental and simulated spectral densities are in reasonable general agreement for NH bond vectors, where the corresponding motions have converged within the simulation time. A sufficient sampling of the motions for NH bonds within flexible parts of the protein requires a longer simulation time. The simulated spectral densities J(0) and J(omega(N)) are, on average, 4.5% and 16% lower than the experimental data; the corresponding numbers for the core residues are about 6%; the high-frequency spectral densities J(0.87omega(H)) are lower by, on average, 16% (21% for the core). The simulated order parameters, S(2), are also lower, although the overall disagreement between the simulation and experiment is less pronounced: 1% for all residues and 6% for the core. The observed systematic decrease of simulated spectral density and the order parameters compared to the experimental data can be partially attributed to the ultrafast librational motion of the NH bonds with respect to their peptide plane, which was analyzed in detail. This systematic difference is most pronounced for J(0.87omega(H)), which appears to be most sensitive to the slow, subnanosecond time scale of internal motion, whereas J(0) and J(omega(N)) are dominated by the overall rotational tumbling of the protein. Similar discrepancies are observed between the experimentally measured (15)N relaxation parameters (R(1), R(2), NOE) and their values calculated from the simulated spectral densities. The analysis of spectral densities provides additional information regarding the comparison of the simulated and experimental data, not available from the model-free analysis.
Collapse
Affiliation(s)
- S Pfeiffer
- The Rockefeller University, 1230 York Avenue, New York, NY 10021-6399, USA
| | | | | |
Collapse
|
5
|
Wu G, MacKenzie R, Durda PJ, Tsang P. The binding of a glycoprotein 120 V3 loop peptide to HIV-1 neutralizing antibodies. Structural implications. J Biol Chem 2000; 275:36645-52. [PMID: 10967109 DOI: 10.1074/jbc.m005369200] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The structural and antigenic properties of a peptide ("CRK") derived from the V3 loop of HIV-1 gp120 protein were studied using NMR and SPR techniques. The sequence of CRK corresponds to the central portion of the V3 loop containing the highly conserved "GPGR" residue sequence. Although the biological significance of this conserved sequence is unknown, the adoption of conserved secondary structure (type II beta-turn) in this region has been proposed. The tendency of CRK (while free or conjugated to protein), to adopt such structure and the influence of such structure upon CRK antigenicity were investigated by NMR and SPR, respectively. Regardless of conjugation, CRK is conformationally averaged in solution but a weak tendency of the CRK "GPGR" residues to adopt a beta-turn conformation was observed after conjugation. The influence of GPGR structure upon CRK antigenicity was investigated by measuring the affinities of two cognate antibodies: "5023A" and "5025A," for CRK, protein-conjugated CRK and gp120 protein. Each antibody bound to all the antigens with nearly the same affinity. From these data, it appears that: (a) antibody binding most likely involves an induced fit of the peptide and (b) the gp120 V3 loop is probably conformationally heterogeneous. Since 5023A and 5025A are HIV-1 neutralizing antibodies, neutralization in these cases appears to be independent of adopted GPGR beta-turn structure.
Collapse
Affiliation(s)
- G Wu
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221-0172, USA
| | | | | | | |
Collapse
|
6
|
Chatfield DC, Wong SE. Methyl Motional Parameters in Crystalline l-Alanine: Molecular Dynamics Simulation and NMR. J Phys Chem B 2000. [DOI: 10.1021/jp0018089] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- David C. Chatfield
- Chemistry Department, Florida International University, Miami, Florida 33199
| | - Sergio E. Wong
- Chemistry Department, Florida International University, Miami, Florida 33199
| |
Collapse
|
7
|
Abstract
NMR relaxation experiments can provide information on overall and internal motions in proteins. This review consists of a concise report on the evolution of the theories for nuclear relaxation followed by an overview of mathematical models for internal motions in proteins. Next, the method of spectral density mapping with recent developments is reviewed. This is followed by a discussion of pulse sequences for relaxation experiments. Finally, we review recent studies correlating relaxation parameters, in particular spectral density functions, with structural features of proteins and with results of molecular dynamics simulations.
Collapse
Affiliation(s)
- K T Dayie
- Program in Higher Degrees in Biophysics, Harvard University, Boston, MA 02115, USA
| | | | | |
Collapse
|
8
|
Ishima R, Yamasaki K, Saito M, Nagayama K. Spectral densities of nitrogen nuclei in Escherichia coli ribonuclease HI obtained by 15N NMR relaxation and molecular dynamics. JOURNAL OF BIOMOLECULAR NMR 1995; 6:217-220. [PMID: 8589610 DOI: 10.1007/bf00211786] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Spectral densities of the 15N amide in Escherichia coli ribonuclease HI, obtained from NMR relaxation experiments, were compared with those calculated using a molecular dynamics (MD) simulation. All calculations and comparisons assumed that the auto-correlation function describing the internal motions of the molecule was independent of the auto-correlation function associated with overall rotational diffusion. Comparisons were limited to those residues for which the auto-correlation function of internal motions rapidly relaxed and reached a steady state within 205 ps. The results show the importance of frequency components as well as amplitudes of internal motions in order to obtain a meaningful comparison of MD simulations with NMR data.
Collapse
Affiliation(s)
- R Ishima
- Nagayama Protein Array Project, ERATO, JRDC, Tsukuba Research Consortium, Pilot Laboratory, Japan
| | | | | | | |
Collapse
|
9
|
Fadel AR, Jin DQ, Montelione GT, Levy RM. Crankshaft motions of the polypeptide backbone in molecular dynamics simulations of human type-alpha transforming growth factor. JOURNAL OF BIOMOLECULAR NMR 1995; 6:221-226. [PMID: 8589611 DOI: 10.1007/bf00211787] [Citation(s) in RCA: 54] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Order parameters for the backbone N-H and C alpha-H bond vectors have been calculated from a 150 ps molecular dynamics (MD) simulation of human type-alpha transforming growth factor in H2O solvent. Two kinds of 'crankshaft motions' of the polypeptide backbone are observed in this MD trajectory. The first involves small-amplitude rocking of the rigid peptide bond due to correlated changes in the backbone dihedral angles psi i-1 and phi i. These high-frequency 'librational crankshaft' motions are correlated with systematically smaller values of motional order parameters for backbone N-H bond vectors compared to C alpha-H bond vectors. In addition, infrequent 'crankshaft flips' of the peptide bond from one local minimum to another are observed for several amino acid residues. These MD simulations demonstrate that comparisons of N-H and C alpha-H order parameters provide a useful approach for identifying crankshaft librational motions in proteins.
Collapse
Affiliation(s)
- A R Fadel
- Department of Chemistry, Wright-Rieman Laboratories, Rutgers University, Piscataway, NJ, USA
| | | | | | | |
Collapse
|
10
|
Beger RD, Balasubramanian S, Bennett SE, Mosbaugh DW, Bolton PH. Tertiary structure of uracil-DNA glycosylase inhibitor protein. J Biol Chem 1995; 270:16840-7. [PMID: 7622499 DOI: 10.1074/jbc.270.28.16840] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
The Bacillus subtilis bacteriophage PBS2 uracil-DNA glycosylase inhibitor (Ugi) is an acidic protein of 84 amino acids that inactivates uracil-DNA glycosylase from diverse organisms. The secondary structure of Ugi consists of five anti-parallel beta-strands and two alpha-helices (Balasubramanian, S., Beger, R.D., Bennett, S.E., Mosbaugh, D.W., and Bolton, P.H. (1995) J. Biol. Chem. 270, 296-303). The tertiary structure of Ugi has been determined by solution state multidimensional nuclear magnetic resonance. The Ugi structure contains an area of highly negative electrostatic potential produced by the close proximity of a number of acidic residues. The unfavorable interactions between these acidic residues are apparently accommodated by the stability of the beta-strands. This negatively charged region is likely to play an important role in the binding of Ugi to uracil-DNA glycosylase.
Collapse
Affiliation(s)
- R D Beger
- Chemistry Department, Wesleyan University, Middletown, Connecticut 06459, USA
| | | | | | | | | |
Collapse
|
11
|
Yamasaki K, Saito M, Oobatake M, Kanaya S. Characterization of the internal motions of Escherichia coli ribonuclease HI by a combination of 15N-NMR relaxation analysis and molecular dynamics simulation: examination of dynamic models. Biochemistry 1995; 34:6587-601. [PMID: 7756290 DOI: 10.1021/bi00020a003] [Citation(s) in RCA: 42] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
The backbone dynamics of Escherichia coli ribonuclease HI (RNase HI) in the picosecond to nanosecond time scale were characterized by a combination of measurements of 15N-NMR relaxation (T1, T2, and NOE), analyzed by a model-free approach, and molecular dynamics (MD) simulation in water. The MD simulations in water were carried out with long-range Coulomb interactions to avoid the artificial fluctuation caused by the cutoff approximation. The model-free analysis of the 15N-NMR relaxation indicated that RNase HI has a rotational correlation time of 10.9 ns at 27 degrees C. The generalized order parameter (S2) for the internal motions varied from 0.15 to 1.0, with an average value of 0.85, which is much larger than that of the RNase H domain of HIV-1 reverse transcriptase (0.78). Large internal motions (small order parameters) were observed in the N-terminal region (Leu2-Lys3), the loop between beta-strands A and B (Cys13-Gly15), the turn between alpha-helix I and beta-strand D (Glu61, His62), the loop between beta-strand D and alpha-helix II (Asp70-Tyr71), the loop between alpha-helices III and IV (Ala93-Lys96), the loop between beta-strand E and alpha-helix V (Gly123-His127), and the C-terminal region (Gln152-Val155). The effective correlation time observed in these regions varied from 0.45 ns (Glu61, Lys96) to 2.2 ns (Leu14). The order parameters calculated from the MD agreed well with those from the NMR experiment, with a few exceptions. The distributions of most of the backbone N-H vectors obtained by MD are approximately consistent with the diffusion-in-a-cone model. These distributions, however, were elliptic, with a long axis perpendicular to the plane defined by the N-H and N-C alpha vectors. Distributions supporting the axial fluctuation model or the jump-between-two-cones model were also observed in the MD simulation.
Collapse
Affiliation(s)
- K Yamasaki
- Protein Engineering Research Institute, Osaka, Japan
| | | | | | | |
Collapse
|
12
|
Ishima R, Nagayama K. Protein backbone dynamics revealed by quasi spectral density function analysis of amide N-15 nuclei. Biochemistry 1995; 34:3162-71. [PMID: 7880811 DOI: 10.1021/bi00010a005] [Citation(s) in RCA: 118] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Spectral density functions J(0), J(omega N), and J(omega H + omega N) of individual amide N-15 nuclei in proteins were approximated by a quasi spectral density function (QSDF). Using this function, the backbone dynamics were analyzed for seven protein systems on which data have been published. We defined J(0; omega N) as the difference between the J(0) and the J(omega N) values, which describes motions slower than 50 (or 60) MHz, and J(omega N; omega H+N) as the difference between the J(omega N) and the J(omega H + omega N) values, which describes motions slower than 450 (or 540) MHz. The QSDF analysis can easily extract the J(0; omega N) of protein backbones, which have often some relation to biologically relevant reactions. Flexible N-terminal regions in eglin c and glucose permease IIA and a loop region in eglin c showed smaller values of both the J(0; omega N) and the J(omega N; omega H+N) as compared with the other regions, indicating increases in motions faster than nanosecond. The values of the J(0; omega N) for the backbone of the FK506 binding protein showed a large variation in the apoprotein but fell in a very narrow range after the binding of FK506. Characteristic increase or decrease in the values of J(0) and J(omega N) was observed in two or three residues located between secondary structures.
Collapse
Affiliation(s)
- R Ishima
- Nagayama Protein Array Project, ERATO, JRDC, Tsukuba Research Consortium, Japan
| | | |
Collapse
|