1
|
Robinson MK, Monroe JI, Shell MS. Are AMBER Force Fields and Implicit Solvation Models Additive? A Folding Study with a Balanced Peptide Test Set. J Chem Theory Comput 2016; 12:5631-5642. [DOI: 10.1021/acs.jctc.6b00788] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Melina K. Robinson
- Department
of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Jacob I. Monroe
- Department
of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - M. Scott Shell
- Department
of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, United States
| |
Collapse
|
2
|
Pepin R, Laszlo KJ, Marek A, Peng B, Bush MF, Lavanant H, Afonso C, Tureček F. Toward a Rational Design of Highly Folded Peptide Cation Conformations. 3D Gas-Phase Ion Structures and Ion Mobility Characterization. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2016; 27:1647-60. [PMID: 27400696 PMCID: PMC5031493 DOI: 10.1007/s13361-016-1437-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/03/2016] [Revised: 05/27/2016] [Accepted: 06/16/2016] [Indexed: 05/18/2023]
Abstract
Heptapeptide ions containing combinations of polar Lys, Arg, and Asp residues with non-polar Leu, Pro, Ala, and Gly residues were designed to study polar effects on gas-phase ion conformations. Doubly and triply charged ions were studied by ion mobility mass spectrometry and electron structure theory using correlated ab initio and density functional theory methods and found to exhibit tightly folded 3D structures in the gas phase. Manipulation of the basic residue positions in LKGPADR, LRGPADK, KLGPADR, and RLGPADK resulted in only minor changes in the ion collision cross sections in helium. Replacement of the Pro residue with Leu resulted in only marginally larger collision cross sections for the doubly and triply charged ions. Disruption of zwitterionic interactions in doubly charged ions was performed by converting the C-terminal and Asp carboxyl groups to methyl esters. This resulted in very minor changes in the collision cross sections of doubly charged ions and even slightly diminished collision cross sections in most triply charged ions. The experimental collision cross sections were related to those calculated for structures of lowest free energy ion conformers that were obtained by extensive search of the conformational space and fully optimized by density functional theory calculations. The predominant factors that affected ion structures and collision cross sections were due to attractive hydrogen bonding interactions and internal solvation of the charged groups that overcompensated their Coulomb repulsion. Structure features typically assigned to the Pro residue and zwitterionic COO-charged group interactions were only secondary in affecting the structures and collision cross sections of these gas-phase peptide ions. Graphical Abstract ᅟ.
Collapse
Affiliation(s)
- Robert Pepin
- Department of Chemistry, Bagley Hall, University of Washington, Seattle, WA, USA
| | - Kenneth J Laszlo
- Department of Chemistry, Bagley Hall, University of Washington, Seattle, WA, USA
| | - Aleš Marek
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic
| | - Bo Peng
- Department of Chemistry, Bagley Hall, University of Washington, Seattle, WA, USA
| | - Matthew F Bush
- Department of Chemistry, Bagley Hall, University of Washington, Seattle, WA, USA
| | - Helène Lavanant
- Laboratoire COBRA CNRS UMR 6014 & FR 3038, Université de Rouen, INSA de Rouen, Mont St Aignan Cedex, France
| | - Carlos Afonso
- Laboratoire COBRA CNRS UMR 6014 & FR 3038, Université de Rouen, INSA de Rouen, Mont St Aignan Cedex, France
| | - František Tureček
- Department of Chemistry, Bagley Hall, University of Washington, Seattle, WA, USA.
| |
Collapse
|
3
|
Duan LL, Gao Y, Ji CG, Mei Y, Zhang QG, Tang B, Zhang JZH. Energetics of protein backbone hydrogen bonds and their local electrostatic environment. Sci China Chem 2014. [DOI: 10.1007/s11426-014-5246-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
4
|
Duan LL, Gao Y, Mei Y, Zhang QG, Tang B, Zhang JZH. Folding of a helix is critically stabilized by polarization of backbone hydrogen bonds: study in explicit water. J Phys Chem B 2012; 116:3430-5. [PMID: 22369598 DOI: 10.1021/jp212516g] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Multiple single-trajectory molecular dynamics (MD) simulation at room temperature (300 K) in explicit water was carried out to study the folding dynamics of an α-helix (PDB 2I9M ) using a polarized charge scheme that includes electronic polarization of backbone hydrogen bonds. Starting from an extended conformation, the 17-residue peptide was successfully folded into the native structure (α-helix) between 80 and 130 ns with a root-mean-square deviation of ~1.0 Å. Analysis of the time-dependent trajectories revealed that helix formation of the peptide started at the terminals and progressed toward the center of the peptide. For comparison, MD trajectories generated under various versions of standard AMBER force fields failed to show any significant or stable helix formation in our simulation. Our result shows clear evidence that the electronic polarization of backbone hydrogen bonds energetically stabilizes the helix formation and is critical to the stable folding of the short helix structure.
Collapse
Affiliation(s)
- Li L Duan
- State Key Laboratory of Precision Spectroscopy, Department of Physics, and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | | | | | | | | | | |
Collapse
|
5
|
Gao Y, Lu X, Duan LL, Zhang JZH, Mei Y. Polarization of intraprotein hydrogen bond is critical to thermal stability of short helix. J Phys Chem B 2011; 116:549-54. [PMID: 22126129 DOI: 10.1021/jp208953x] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Simulation result for protein folding/unfolding is highly dependent on the accuracy of the force field employed. Even for the simplest structure of protein such as a short helix, simulations using the existing force fields often fail to produce the correct structural/thermodynamic properties of the protein. Recent research indicated that lack of polarization is at least partially responsible for the failure to successfully fold a short helix. In this work, we develop a simple formula-based atomic charge polarization model for intraprotein (backbone) hydrogen bonding based on the existing AMBER force field to study the thermal stability of a short helix (2I9M) by replica exchange molecular dynamics simulation. By comparison of the simulation results with those obtained by employing the standard AMBER03 force field, the formula-based atomic charge polarization model gave the helix melting curve in close agreement with the NMR experiment. However, in simulations using the standard AMBER force field, the helix was thermally unstable at the temperature of the NMR experiment, with a melting temperature almost below the freezing point. The difference in observed thermal stability from these two simulations is the effect of backbone intraprotein polarization, which was included in the formula-based atomic charge polarization model. The polarization of backbone hydrogen bonding thus plays a critical role in the thermal stability of helix or more general protein structures.
Collapse
Affiliation(s)
- Ya Gao
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China
| | | | | | | | | |
Collapse
|
6
|
Duan LL, Mei Y, Zhang D, Zhang QG, Zhang JZH. Folding of a helix at room temperature is critically aided by electrostatic polarization of intraprotein hydrogen bonds. J Am Chem Soc 2010; 132:11159-64. [PMID: 20698682 DOI: 10.1021/ja102735g] [Citation(s) in RCA: 81] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report direct folding of a 17-residue helix protein (pdb:2I9M) by standard molecular dynamics simulation (single trajectory) at room temperature with implicit solvent. Starting from a fully extended structure, 2I9M successfully folds into the native conformation within 16 ns using adaptive hydrogen bond-specific charges to take into account the electrostatic polarization effect. Cluster analysis shows that conformations in the native state cluster have the highest population (78.4%) among all sampled conformations. Folding snapshots and the secondary structure analysis demonstrate that the folding of 2I9M begins at terminals and progresses toward the center. A plot of the free energy landscape indicates that there is no significant free energy barrier during folding, which explains the observed fast folding speed. For comparison, exactly the same molecular dynamics simulation but carried out under existing AMBER charges failed to fold 2I9M into native-like structures. The current study demonstrates that electrostatic polarization of intraprotein hydrogen bonding, which stabilizes the helix, is critical to the successful folding of 2I9m.
Collapse
Affiliation(s)
- Li L Duan
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | | | | | | | | |
Collapse
|
7
|
Xiang Y, Duan L, Zhang JZH. Folding dynamics of a small protein at room temperature via simulated coherent two-dimensional infrared spectroscopy. Phys Chem Chem Phys 2010; 12:15681-8. [PMID: 20676442 DOI: 10.1039/c0cp00375a] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Understanding protein folding is of fundamental and practical importance in chemistry and biology. Despite the great success that has been made in tackling this problem, a detailed knowledge of how the elementary processes such as hydrogen-bond formation occur during protein folding has remained largely elusive. Using the combined power of molecular dynamics simulation with electrostatic polarization and coherent two-dimensional infrared spectroscopy, we are able to delineate the order of the hydrogen-bond formation event of a 17-residue peptide during its folding from an extended state to the native α-helix state. The folding is carried out by a single trajectory room-temperature molecular dynamics simulation that includes the polarization effect of hydrogen bonding, which is critical to the successful folding of the peptide. The onset and evolution of the isotope-labeled amide I vibration diagonal and cross peaks on the simulated 2DIR spectra allow us to build a structure-spectrum connection, and thus provide a microscopic picture of the helix folding process.
Collapse
Affiliation(s)
- Yun Xiang
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Sciences, East China Normal University, Shanghai 200062, China.
| | | | | |
Collapse
|