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 条
[71]   PARAMETERS OF HELIX-COIL TRANSITION THEORY FOR ALANINE-BASED PEPTIDES OF VARYING CHAIN LENGTHS IN WATER [J].
SCHOLTZ, JM ;
QIAN, H ;
YORK, EJ ;
STEWART, JM ;
BALDWIN, RL .
BIOPOLYMERS, 1991, 31 (13) :1463-1470
[72]   DISTRIBUTION OF HELICITY WITHIN THE MODEL PEPTIDE ACETYL(AAQAA)(3)AMIDE [J].
SHALONGO, W ;
DUGAD, L ;
STELLWAGEN, E .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1994, 116 (18) :8288-8293
[73]   A test on peptide stability of AMBER force fields with implicit solvation [J].
Shell, M. Scott ;
Ritterson, Ryan ;
Dill, Ken A. .
JOURNAL OF PHYSICAL CHEMISTRY B, 2008, 112 (22) :6878-6886
[74]   SPARTA plus : a modest improvement in empirical NMR chemical shift prediction by means of an artificial neural network [J].
Shen, Yang ;
Bax, Ad .
JOURNAL OF BIOMOLECULAR NMR, 2010, 48 (01) :13-22
[75]   Exploring the helix-coil transition via all-atom equilibrium ensemble simulations [J].
Sorin, EJ ;
Pande, VS .
BIOPHYSICAL JOURNAL, 2005, 88 (04) :2472-2493
[76]  
Tanford C., 1980, HYDROPHOBIC EFFECT F, V2
[77]  
Thompson P. A., 2000, BIOCHEMISTRY-US, V104, P9200
[78]   The helix-coil kinetics of a heteropeptide [J].
Thompson, PA ;
Muñoz, V ;
Jas, GS ;
Henry, ER ;
Eaton, WA ;
Hofrichter, J .
JOURNAL OF PHYSICAL CHEMISTRY B, 2000, 104 (02) :378-389
[79]   A systematic study of water models for molecular simulation: Derivation of water models optimized for use with a reaction field [J].
van der Spoel, D ;
van Maaren, PJ ;
Berendsen, HJC .
JOURNAL OF CHEMICAL PHYSICS, 1998, 108 (24) :10220-10230
[80]   Vapor-liquid equilibria from the triple point up to the critical point for the new generation of TIP4P-like models: TIP4P/Ew, TIP4P/2005, and TIP4P/ice [J].
Vega, C. ;
Abascal, J. L. F. ;
Nezbeda, I. .
JOURNAL OF CHEMICAL PHYSICS, 2006, 125 (03)