Secondary Structure Propensities in Peptide Folding Simulations: A Systematic Comparison of Molecular Mechanics Interaction Schemes

We present a systematic study directed toward the secondary structure propensity and sampling behavior in peptide folding simulations with eight different molecular dynamics force-field variants in explicit solvent. We report on the combinational result of force field, water model, and electrostatic interaction schemes and compare to available experimental characterization of five studied model peptides in terms of reproduced structure and dynamics. The total simulation time exceeded 18 μs and included simulations that started from both folded and extended conformations. Despite remaining sampling issues, a number of distinct trends in the folding behavior of the peptides emerged. Pronounced differences in the propensity of finding prominent secondary structure motifs in the different applied force fields suggest that problems point in particular to the balance of the relative stabilities of helical and extended conformations.     

Why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution when comparing with the GROMOS 43A1 force field?

The optimal combination of force field and water model is an essential problem that is able to increase molecular dynamics simulation quality for different types of proteins and peptides. In this work, an attempt has been made to explore the problem by studying H1 peptide using four different models based on different force fields, water models and electrostatic schemes. The driving force for H1 peptide conformation transition and the reason why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution while the GROMOS 43A1 force field can do were investigated by temperature replica exchange molecular dynamics simulation (T-REMD). The simulation using the GROMOS 43A1 force field preferred to adopt a β-hairpin structure, which was in good agreement with the several other simulations and the experimental evidences. However, the simulation using the OPLS-AA force field has a significant difference from the simulations with the GROMOS 43A1 force field simulation. The results show that the driving force in H1 peptide conformation transition is solvent exposure of its hydrophobic residues. However, the subtle balances between residue-residue interactions and residue-solvent interaction are disrupted by using the OPLS-AA force field, which induced the reduction in the number of residue-residue contact. Similar solvent exposure of the hydrophobic residues is observed for all the conformations sampled using the OPLS-AA force field. For H1 peptide which exhibits large solvent exposure of the hydrophobic residues, the GROMOS 43A1 force field with the SPC water model can provide more accurate results.     

Toward a predictive understanding of slow methyl group dynamics in proteins

The development of the most recent generation of molecular mechanics force fields promises an increasingly predictive understanding of the protein dynamics-function relationship. Based on extensive validation against various types of experimental data, the AMBER force field ff99SB was benchmarked in recent years as a favorable force field for protein simulations. Recent improvements of the side chain and backbone potentials, made by different groups, led to the ff99SB-ILDN and ff99SBnmr1 force fields, respectively. The combination of these potentials into a unified force field, termed ff99SBnmr1-ILDN, was used in this study to perform a microsecond time scale molecular dynamics simulation of free ubiquitin in explicit solvent for validation against an extensive set of experimental NMR methyl group residual dipolar couplings. Our results show a high level of consistency between the experimental data and the values predicted from the molecular dynamics trajectory reflecting a systematically improved performance of ff99SBnmr1-ILDN over the original ff99SB force field. Moreover, the unconstrained ff99SBnmr1-ILDN MD ensemble achieves a similar level of agreement as the recently introduced EROS ensemble, which was constructed based on a large body of NMR data as constraints, including the methyl residual dipolar couplings. This suggests that ff99SBnmr1-ILDN provides a high-quality representation of the motions of methyl-bearing protein side chains, which are sensitive probes of protein-protein and protein-ligand interactions.     

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.     

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

Effects of force fields on the conformational and dynamic properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1‐40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER‐ff99sb‐ILDN, AMBER‐ff99sb*‐ILDN, AMBER‐ff99sb‐NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER‐ff99sb‐ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α‐helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER‐ff99sb‐NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER‐ff99sb‐NMR force field, the others tended to under estimate the expected amount of β‐sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER‐ff99sb‐NMR, reproduce a theoretically expected β‐sheet‐turn‐β‐sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C‐terminal hydrophobic cores from residues 17‐21 and 30‐36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.     

Improved side‐chain torsion potentials for the Amber ff99SB protein force field

Recent advances in hardware and software have enabled increasingly long molecular dynamics (MD) simulations of biomolecules, exposing certain limitations in the accuracy of the force fields used for such simulations and spurring efforts to refine these force fields. Recent modifications to the Amber and CHARMM protein force fields, for example, have improved the backbone torsion potentials, remedying deficiencies in earlier versions. Here, we further advance simulation accuracy by improving the amino acid side‐chain torsion potentials of the Amber ff99SB force field. First, we used simulations of model alpha‐helical systems to identify the four residue types whose rotamer distribution differed the most from expectations based on Protein Data Bank statistics. Second, we optimized the side‐chain torsion potentials of these residues to match new, high‐level quantum‐mechanical calculations. Finally, we used microsecond‐timescale MD simulations in explicit solvent to validate the resulting force field against a large set of experimental NMR measurements that directly probe side‐chain conformations. The new force field, which we have termed Amber ff99SB‐ILDN, exhibits considerably better agreement with the NMR data. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Are Current Molecular Dynamics Force Fields too Helical?

Accurate force fields are essential for the success of molecular dynamics simulations. In apparent contrast to the conformational preferences of most force fields, recent NMR experiments suggest that short polyalanine peptides in water populate the polyproline II structure almost exclusively. To investigate this apparent contradiction, with its ramifications for the assessment of molecular force fields and the structure of unfolded proteins, we performed extensive simulations of Ala₅ in water (∼5 μs total time), using twelve different force fields and three different peptide terminal groups. Using either empirical or density-functional-based Karplus relations for the J-couplings, we find that most current force fields do overpopulate the α-region, with quantitative results depending on the choice of Karplus relation and on the peptide termini. Even after reweighting to match experiment, we find that Ala₅ retains significant α- and β-populations. In fact, several force fields match the experimental data well before reweighting and have a significant helical population. We conclude that radical changes to the best current force fields are not necessary, based on the NMR data. Nevertheless, experiments on short peptides open the way toward the systematic improvement of current simulation models.     

Comparison of force fields for Alzheimer's A : A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three‐dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN‐NMR and CHARMM22*, in their ability to model Aβ₄₂, an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN‐NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.     

Free-energy landscape of the GB1 hairpin in all-atom explicit solvent simulations with different force fields: Similarities and differences

Although it is now possible to fold peptides and miniproteins in molecular dynamics simulations, it is well appreciated that force fields are not all transferable to different proteins. Here, we investigate the influence of the protein force field and the solvent model on the folding energy landscape of a prototypical two‐state folder, the GB1 hairpin. We use extensive replica‐exchange molecular dynamics simulations to characterize the free‐energy surface as a function of temperature. Most of these force fields appear similar at a global level, giving a fraction folded at 300 K between 0.2 and 0.8 in all cases, which is a difference in stability of 2.8 kT, and are generally consistent with experimental data at this temperature. The most significant differences appear in the unfolded state, where there are different residual secondary structures which are populated, and the overall dimensions of the unfolded states, which in most of the force fields are too collapsed relative to experimental Förster Resonance Energy Transfer (FRET) data. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Test of molecular dynamics force fields in gramicidin A

The force fields commonly used in molecular dynamics simulations of proteins are optimized under bulk conditions. Whether the same force fields can be used in simulations of membrane proteins is not well established, although they are increasingly being used for such purposes. Here we consider ion permeation in the gramicidin A channel as a test of the AMBER force field in a membrane environment. The potentials of mean force for potassium ions are calculated along the channel axis and compared with the one deduced from the experimental conductance data. The calculated result indicates a rather large central barrier similar to those obtained from other force fields, which are incompatible with the conductance data. We suggest that lack of polarizability is the most likely cause of this problem, and, therefore, urge development of polarizable force fields for simulations of membrane proteins.     

Structure and dynamics of two β-peptides in solution from molecular dynamics simulations validated against experiment

We have studied two different beta-peptides in methanol using explicit solvent molecular dynamics simulations and the GROMOS 53A6 force field: a heptapeptide (peptide 1) expected to form a left-handed 3(14)-helix, and a hexapeptide (peptide 2) expected to form a beta-hairpin in solution. Our analysis has focused on identifying and analyzing the stability of the dominant secondary structure conformations adopted by the peptides, as well as on comparing the experimental NOE distance upper bounds and 3J-coupling values with their counterparts calculated on the basis of the simulated ensembles. Moreover, we have critically compared the present results with the analogous results obtained with the GROMOS 45A3 (peptide 1) and 43A1 (peptide 2) force fields. We conclude that within the limits of conformational sampling employed here, the GROMOS 53A6 force field satisfactorily reproduces experimental findings regarding the behavior of short beta-peptides, with accuracy that is comparable to but not exceeding that of the previous versions of the force field. GCE legend Conformational clustering analysis of the simulated ensemble of a ss-hexapeptide with two different simulation setups (a and b). The central members of all of the clusters populating more than 5% of all of the structures are shown, together with the most dominant hydrogen bonds and the corresponding percentages of cluster members containing them.     

Force field influences in beta-hairpin folding simulations

All-atom force fields are now routinely used for more detailed understanding of protein folding mechanisms. However, it has been pointed out that use of all-atom force fields does not guarantee more accurate representations of proteins; in fact, sometimes it even leads to biased structural distributions. Indeed, several issues remain to be solved in force field developments, such as accurate treatment of implicit solvation for efficient conformational sampling and proper treatment of backbone interactions for secondary structure propensities. In this study, we first investigate the quality of several recently improved backbone interaction schemes in AMBER for folding simulations of a beta-hairpin peptide, and further study their influences on the peptide's folding mechanism. Due to the significant number of simulations needed for a thorough analysis of tested force fields, the implicit Poisson-Boltzmann solvent was used in all simulations. The chosen implicit solvent was found to be reasonable for studies of secondary structures based on a set of simulations of both alpha-helical and beta-hairpin peptides with the TIP3P explicit solvent as benchmark. Replica exchange molecular dynamics was also utilized for further efficient conformational sampling. Among the tested AMBER force fields, ff03 and a revised ff99 force field were found to produce structural and thermodynamic data in comparably good agreement with the experiment. However, detailed folding pathways, such as the order of backbone hydrogen bond zipping and the existence of intermediate states, are different between the two force fields, leading to force field-dependent folding mechanisms.     

Structural equilibrium of DNA represented with different force fields.

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.     

Thermal unfolding simulations of apo-calmodulin using leap-dynamics

The simulation method leap-dynamics (LD) has been applied to protein thermal unfolding simulations to investigate domain-specific unfolding behavior. Thermal unfolding simulations of the 148-residue protein apo-calmodulin with implicit solvent were performed at temperatures 290 K, 325 K, and 360 K and compared with the corresponding molecular dynamics trajectories in terms of a number of calculated conformational parameters. The main experimental results of unfolding are reproduced in showing the lower stability of the C-domain: at 290 K, both the N- and C-domains are essentially stable; at 325 K, the C-domain unfolds, whereas the N-domain remains folded; and at 360 K, both domains unfold extensively. This behavior could not be reproduced by molecular dynamics simulations alone under the same conditions. These results show an encouraging degree of convergence between experiment and LD simulation. The simulations are able to describe the overall plasticity of the apo-calmodulin structure and to reveal details such as reversible folding/unfolding events within single helices. The results show that by using the combined application of a fast and efficient sampling routine with a detailed molecular dynamics force field, unfolding simulations of proteins at atomic resolution are within the scope of current computational power.     

Interchain hydrophobic clustering promotes rigidity in HIV-1 protease flap dynamics: new insights from Molecular Dynamics

The dynamics of HIV-1 protease (HIV-pr), a drug target for HIV infection, has been studied extensively by both computational and experimental methods. The flap dynamics of HIV-pr is considered to be more important for better ligand binding and enzymatic actions. Moreover, it has been demonstrated that the drug-induced mutations can change the flap dynamics of HIV-pr affecting the binding affinity of the ligands. Therefore, detailed understanding of flap dynamics is essential for designing better inhibitors. Previous computational investigations observed significant variation in the flap opening in nanosecond time scale indicating that the dynamics is highly sensitive to the simulation protocols. To understand the sensitivity of the flap dynamics on the force field and simulation protocol, molecular dynamics simulations of HIV-pr have been performed with two different AMBER force fields, ff99 and ff02. Two different trajectories (20?ns each) were obtained using the ff99 and ff02 force field. The results showed polarizable force field (ff02) make the flap tighter than the nonpolarizable force field (ff99). Some polar interactions and hydrogen bonds involving flap residues were found to be stronger with ff02 force field. The formation of interchain hydrophobic cluster (between flap tip of one chain and active site wall of another chain) was found to be dominant in the semi-open structures obtained from the simulations irrespective of the force field. It is proposed that an inhibitor, which will promote this interchain hydrophobic clustering, may make the flaps more rigid, and presumably the effect of mutation would be small on ligand binding.     

Molecular dynamics simulations for the prediction of the dielectric spectra of alcohols,glycols and monoethanolamine

The response of molecular systems to electromagnetic radiation in the microwave region (0.3–300 GHz) has been principally studied experimentally, using broadband dielectric spectroscopy. However, relaxation times corresponding to reorganisation of molecular dipoles due to their interaction with electromagnetic radiation at microwave frequencies are within the scope of modern molecular simulations. In this work, fluctuations of the total dipole moment of a molecular system, obtained through molecular dynamics simulations, are used to determine the dielectric spectra of water, a series of alcohols and glycols, and monoethanolamine. Although the force fields employed in this study have principally been developed to describe thermodynamic properties, most them give fairly good predictions of this dynamical property for these systems. However, the inaccuracy of some models and the long simulation times required for the accurate estimation of the static dielectric constant can sometimes be problematic. We show that the use of the experimental value for the static dielectric constant in the calculations, instead of the one predicted by the different models, yields satisfactory results for the dielectric spectra, and hence the heat absorbed from microwaves, avoiding the need for extraordinarily long simulations or re-calibration of molecular models.     

Computational investigations of structural changes resulting from point mutations in a collagen-like peptide

The results of 0.5-1.0 ns molecular dynamics simulations of the collagen-like peptides [(POG)4(POA)(POG)4]3 and [(POG)9]3 (POG: proline-hydroxyproline-glycine) are presented. All simulations were performed using the AMBER-94 molecular mechanical force field with a shell of TIP3P waters surrounding the peptides. The initial geometries for the collagen-like peptides included an x-ray crystallographic structure, a computer-generated structure, a [(POG)9]3 structure modeled from the x-ray structure, and the x-ray structure with crystallographic waters replaced with a shell of modeled TIP3P waters. We examined the molecular dynamics peptide residue rms deviation fluctuations, dihedral angles, molecular and chain end-to-end distances, helical parameters, and peptide-peptide and peptide-solvent hydrogen-bonding patterns. Our molecular dynamics simulations of [(POG)4(POA)(POG)4]3 show average structures and internal coordinates similar to the x-ray crystallographic structure. Our results demonstrate that molecular dynamics can be used to reproduce the experimental structures of collagen-like peptides. We have demonstrated the feasibility of using the AMBER-94 molecular mechanical force field, which was parameterized to model nucleic acids and globular proteins, for fibril proteins. We provide a new interpretation of peptide-solvent hydrogen bonding and a peptide-peptide hydrogen bonding pattern not previously reported in x-ray studies. Last, we report on the differences; in particular with respect to main-chain dihedral angles and hydrogen bonding, between the native and mutant collagen-like peptides.