Protein Simulations with an Optimized Water Model: Cooperative Helix Formation and Temperature-Induced Unfolded State Collapse

被引:211
作者
Best, Robert B. [1 ]
Mittal, Jeetain [2 ]
机构
[1] Univ Cambridge, Dept Chem, Cambridge CB2 1EW, England
[2] Lehigh Univ, Dept Chem Engn, Bethlehem, PA 18015 USA
关键词
ALANINE-BASED PEPTIDES; VARYING CHAIN LENGTHS; AMBER FORCE-FIELDS; COIL TRANSITION; EMPIRICAL PARAMETERS; FOLDING PROBLEM; PRESSURE DENATURATION; MOLECULAR SIMULATION; DIPOLAR COUPLINGS; DYNAMICS;
D O I
10.1021/jp108618d
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
A recognized shortcoming in current protein simulations is that most force fields are parametrized with relatively primitive three-site water models. Since the deficiencies of the common three-site water models in reproducing the phase diagram of water are well-known, an improved description of the solvent will be required, for example, to study proteins in molecular simulations at thermodynamic conditions other than standard temperature and pressure. Here, we combine a protein force field derived from Amber ff03 together with the highly optimized TIP4P/2005 water model, with a small backbone modification to match the population of helical states obtained with the new water model to experiment. Remarkably, we find that the resulting force field, Amber ff03w, produces a more cooperative helix coil transition, compared with the similarly "backbone-corrected" Amber ff03* model with TIP3P water, with calculated helix propagation parameters in good agreement with the experiment. The radius of gyration for nonhelical conformations is significantly larger for Amber ff03w than for Amber ff03* and shows a collapse with increasing temperature as found in single-molecule experiments on longer proteins. The origin of the collapse appears to be a more favorable enthalpic component of the peptide solvent interaction and is correlated with increasing turn formation, in accord with the experiment. In addition to this enhanced cooperativity, we verify that, with the new force field, replica exchange folding simulations of the GB1 hairpin and Trp cage result in folded structures, starting from completely unfolded initial conditions; simulations of folded proteins are also stable. These results together suggest that Amber ff03w (with TIP4P/2005) will be well suited for studying protein folding and properties of unfolded state and intrinsically disordered proteins over a wide range of thermodynamic conditions.
引用
收藏
页码:14916 / 14923
页数:8
相关论文
共 82 条
[1]   A general purpose model for the condensed phases of water: TIP4P/2005 [J].
Abascal, JLF ;
Vega, C .
JOURNAL OF CHEMICAL PHYSICS, 2005, 123 (23)
[2]   Properties of ices at 0 K: A test of water models [J].
Aragones, J. L. ;
Noya, E. G. ;
Abascal, J. L. F. ;
Vega, C. .
JOURNAL OF CHEMICAL PHYSICS, 2007, 127 (15)
[3]   Assessing the thermodynamic signatures of hydrophobic hydration for several common water models [J].
Ashbaugh, Henry S. ;
Collett, Nicholas J. ;
Hatch, Harold W. ;
Staton, Jennifer A. .
JOURNAL OF CHEMICAL PHYSICS, 2010, 132 (12)
[4]   Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities [J].
Avbelj, F ;
Luo, PZ ;
Baldwin, RL .
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, 2000, 97 (20) :10786-10791
[5]   Dipolar couplings in macromolecular structure determination [J].
Bax, A ;
Kontaxis, G ;
Tjandra, N .
NUCLEAR MAGNETIC RESONANCE OF BIOLOGICAL MACROMOLECULES, PT B, 2001, 339 :127-174
[6]  
Berendsen H., 1981, INTERMOL FORCES
[7]   GROMACS - A MESSAGE-PASSING PARALLEL MOLECULAR-DYNAMICS IMPLEMENTATION [J].
BERENDSEN, HJC ;
VANDERSPOEL, D ;
VANDRUNEN, R .
COMPUTER PHYSICS COMMUNICATIONS, 1995, 91 (1-3) :43-56
[8]   THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS [J].
BERENDSEN, HJC ;
GRIGERA, JR ;
STRAATSMA, TP .
JOURNAL OF PHYSICAL CHEMISTRY, 1987, 91 (24) :6269-6271
[9]   Are current molecular dynamics force fields too helical? [J].
Best, Robert B. ;
Buchete, Nicolae-Viorel ;
Hummer, Gerhard .
BIOPHYSICAL JOURNAL, 2008, 95 (01) :L7-L9
[10]   Balance between a and β Structures in Ab Initio Protein Folding [J].
Best, Robert B. ;
Mittal, Jeetain .
JOURNAL OF PHYSICAL CHEMISTRY B, 2010, 114 (26) :8790-8798