Experimental Verification of Force Fields for Molecular Dynamics Simulations Using Gly-Pro-Gly-Gly

被引:37
作者
Aliev, Abil E. [1 ]
Courtier-Murias, Denis [1 ]
机构
[1] UCL, Dept Chem, London WC1H 0AJ, England
基金
英国工程与自然科学研究理事会;
关键词
NMR COUPLING-CONSTANTS; CHEMICAL-SHIFTS; SUBSTITUENT ELECTRONEGATIVITIES; CONFORMATIONAL PREFERENCES; NUCLEIC-ACIDS; PROTEIN; WATER; BACKBONE; PROLINE; MECHANICS;
D O I
10.1021/jp101581h
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
Experimental NMR verification of MD simulations using 12 different force fields (AMBER, CHARMM, GROMOS, and OPLS-AA) and 5 different water models has been undertaken to identify reliable MD protocols for structure and dynamics elucidations of small open chain peptides containing Gly and Pro. A conformationally flexible tetrapeptide Gly-Pro-Gly-Gly was selected for NMR (3)J-coupling, chemical shift, and internuclear distance measurements, followed by their calculations using 2 us long MD simulations in water, In addition, Ramachandran population maps for Pro-2 and Gly-3 residues of GPGG obtained from MD simulations were used for detailed comparisons with similar maps from the protein data bank (PDB) for large number of Gly and Pro residues in proteins. The MD simulations revealed strong dependence of the populations and geometries of preferred backbone and side chain conformations, as well as the time scales of the peptide torsional transitions on the force field used. On the basis of the analysis of the measured and calculated data, AMBER99SB is identified as the most reliable force field for reproducing NMR measured parameters, which are dependent on the peptide backbone and the Pro side chain geometries and dynamics. Ramachandran maps showing the dependence of conformational populations as a function of backbone phi/psi angles for Pro-2 and Gly-3 residues of GPGG from MD simulations using AMBER99SB. AMBER03, and CHARMM were found to resemble similar maps for Gly and Pro residues from the PDB survey. Three force fields (AMBER99, AMBER99 phi, and AMBER94) showed the least satisfactory agreement with both the solution NMR and the PDB survey data. The poor performance of these force fields is attributed to their propensity to overstabilize helical peptide backbone conformations at the Pro-2 and Gly-3 residues. On the basis of the similarity of the MD and PDB Ramachandran plots, the following sequence of transitions is suggested for the Gly backbone conformation: alpha(L) reversible arrow beta(pR) reversible arrow beta(s) reversible arrow beta(p) reversible arrow alpha where backbone secondary structures alpha(L), and alpha are associated with helices and turns, beta(P) and beta(PR) correspond to the left- and right-handed polyproline II structures and beta(s) denotes the fully stretched backbone conformation. Compared to the force field dependence, less significant, but noteworthy, variations in the populations of the peptide backbone conformations were observed. For different solvent models considered, a correlation was noted between the number of torsional transitions in GPGG and the water self-diffusion coefficient on using TIP3P, TIP4P, and TIP5P models. In addition to MD results, we also report DFT derived Karplus relationships for Gly and Pro residues using B972 and B3LYP functionals.
引用
收藏
页码:12358 / 12375
页数:18
相关论文
共 102 条
[71]   A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force-field parameter sets 53A5 and 53A6 [J].
Oostenbrink, C ;
Villa, A ;
Mark, AE ;
Van Gunsteren, WF .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 2004, 25 (13) :1656-1676
[72]   A NEW ANALYSIS OF PROTON CHEMICAL-SHIFTS IN PROTEINS [J].
OSAPAY, K ;
CASE, DA .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1991, 113 (25) :9436-9444
[73]   A new set of molecular mechanics parameters for hydroxyproline and its use in molecular dynamics simulations of collagen-like peptides [J].
Park, S ;
Radmer, RJ ;
Klein, TE ;
Pande, VS .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 2005, 26 (15) :1612-1616
[74]   POLYMORPHIC TRANSITIONS IN SINGLE-CRYSTALS - A NEW MOLECULAR-DYNAMICS METHOD [J].
PARRINELLO, M ;
RAHMAN, A .
JOURNAL OF APPLIED PHYSICS, 1981, 52 (12) :7182-7190
[75]  
Ponder JW, 2003, ADV PROTEIN CHEM, V66, P27
[76]  
Press W. H., 1996, NUMERICAL RECIPES FO
[77]   Modern protein force fields behave comparably in molecular dynamics simulations [J].
Price, DJ ;
Brooks, CL .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 2002, 23 (11) :1045-1057
[78]   Force Field-Dependant Structural Divergence Revealed During Long Time Simulations of Calbindin d9k [J].
Project, Elad ;
Nachliel, Esther ;
Gutman, Menachem .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 2010, 31 (09) :1864-1872
[79]   STEREOCHEMICAL CRITERIA FOR POLYPEPTIDE AND PROTEIN CHAIN CONFORMATIONS .2. ALLOWED CONFORMATIONS FOR A PAIR OF PEPTIDE UNITS [J].
RAMAKRISHNAN, C ;
RAMACHANDRAN, GN .
BIOPHYSICAL JOURNAL, 1965, 5 (06) :909-+
[80]   Development and validation of reliable quantum mechanical approaches for the study of free radicals in solution [J].
Rega, N ;
Cossi, M ;
Barone, V .
JOURNAL OF CHEMICAL PHYSICS, 1996, 105 (24) :11060-11067