Free energy surfaces of β‐hairpin and α‐helical peptides generated by replica exchange molecular dynamics with the AGBNP implicit solvent model

Citations Over TimeTop 10% of 2004 papers

Abstract

We have studied the potential of mean force of two peptides, one known to adopt a beta-hairpin and the other an alpha-helical conformation in solution. These peptides are, respectively, residues 41-56 of the C-terminus (GEWTYDDATKTFTVTE) of the B1 domain of protein G and the 13 residue C-peptide (KETAAAKFERQHM) of ribonuclease A. Extensive canonical ensemble sampling has been performed using a parallel replica exchange method. The effective potential employed in this work consists of the OPLS all-atom force field (OPLS-AA) and an analytical generalized Born (AGB) implicit solvent model including a novel nonpolar solvation free energy estimator (NP). An additional dielectric screening parameter has been incorporated into the AGBNP model. In the case of the beta-hairpin, the nonpolar solvation free energy estimator provides the necessary effective interactions for the collapse of the hydrophobic core (W43, Y45, F52, and V54), which the more commonly used surface-area-dependent nonpolar model does not provide. For both the beta-hairpin and the alpha-helix, increased dielectric screening reduces the stability of incorrectly formed salt bridges, which tend to disrupt the formation of the hairpin and helix, respectively. The fraction of beta-hairpin and alpha-helix content we obtained using the AGBNP model agrees well with experimental results. The thermodynamic stability of the beta-hairpin from protein G and the alpha-helical C-peptide from ribonuclease A as modeled with the OPLS-AA/AGBNP effective potential reflects the balance between the nonpolar effective potential terms, which drive compaction, and the polar and hydrogen bonding terms, which promote secondary structure formation.

Related Papers

• → Implementation of the CHARMM Force Field in GROMACS: Analysis of Protein Stability Effects from Correction Maps, Virtual Interaction Sites, and Water Models(2010)1,250 cited
• → First-Shell Solvation of Ion Pairs: Correction of Systematic Errors in Implicit Solvent Models(2004)78 cited
• → Assessment of the Transferability of a Protein Force Field for the Simulation of Peptide-Surface Interactions(2010)64 cited
• → Consistent free energy landscapes and thermodynamic properties of small proteins based on a single all-atom force field employing an implicit solvation(2007)26 cited
• → Direct folding simulation of α‐helices and β‐hairpins based on a single all‐atom force field with an implicit solvation model(2006)35 cited