A test on peptide stability of AMBER force fields with implicit solvation

We used replica exchange molecular dynamics (REMD) simulations to evaluate four different AMBER force fields and three different implicit solvent models. Our aim was to determine if these physics-based models captured the correct secondary structures of two alpha-helical and two beta-peptides: the 14-mer EK helix of Baldwin and co-workers, the C-terminal helix of ribonuclease, the 16-mer C-terminal hairpin of protein G, and the trpzip2 miniprotein. The different models gave different results, but generally we found that AMBER ff96 plus the implicit solvent model of Onufriev, Bashford, and Case gave reasonable structures, and is fairly well-balanced between helix and sheet. We also observed differences in the strength of ion pairing in the solvent models, we but found that the native secondary structures were retained even when salt bridges were prevented in the conformational sampling. Overall, this work indicates that some of these all-atom physics-based force fields may be good starting points for protein folding and protein structure prediction.     

On the influence of charged side chains on the folding-unfolding equilibrium of beta-peptides: a molecular dynamics simulation study

The influence of charged side chains on the folding-unfolding equilibrium of beta-peptides was investigated by means of molecular dynamics simulations. Four different peptides containing only negatively charged side chains, positively charged side chains, both types of charged side chains (with the ability to form stabilizing salt bridges) or no charged side chains were studied under various conditions (different simulation temperatures, starting structures and solvent environment). The NMR solution structure in methanol of one of the peptides (A) has already been published; the synthesis and NMR analysis of another peptide (B) is described here. The other peptides (C and D) studied herein have hitherto not been synthesized. All four peptides A-D are expected to adopt a left-handed 3(14)-helix in solution as well as in the simulations. The resulting ensembles of structures were analyzed in terms of conformational space sampled by the peptides, folding behavior, structural properties such as hydrogen bonding, side chain-side chain and side chain-backbone interactions and in terms of the level of agreement with the NMR data available for two of the peptides. It was found that the presence of charged side chains significantly slows down the folding process in methanol solution due to the stabilization of intermediate conformers with side chain-backbone interactions. In water, where the solvent competes with the solute-solute polar interactions, the folding process to the 3(14)-helix is faster in the simulations.     

Beta‐hairpin restraint potentials for calculations of potentials of mean force as a function of beta‐hairpin tilt,rotation, and distance

We have developed a set of restraint potentials for β‐hairpin tilt relative to the membrane normal, β‐hairpin rotation around the β‐hairpin axis, and hairpin–hairpin distance. Such restraint potentials enable us to characterize the molecular basis of specific β‐hairpin tilt and rotation in membranes and hairpin–hairpin interactions at the atomic level by sampling their conformational space along these degrees of freedom, i.e., reaction coordinates, during molecular dynamics simulations. We illustrate the efficacy of the β‐hairpin restraint potentials by calculating the potentials of mean force (PMFs) as a function of tilt and rotation angles of protegrin‐1 (PG‐1), a β‐hairpin antimicrobial peptide, in an implicit membrane model. The peptide association in the membrane is also examined by calculating the PMFs as a function of distance between two PG‐1 peptides in various dimer interfaces. These novel restraint potentials are found to perform well in each of these cases and are expected to be a useful means to study the microscopic driving forces of insertion, tilting, and rotation of β‐hairpin peptides in membranes as well as their association in aqueous solvent or membrane environments particularly when combined with explicit solvent models. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Establishing effective simulation protocols for beta- and alpha/beta-peptides. II. Molecular mechanical (MM) model for a cyclic beta-residue

All-atom molecular mechanical (MM) force field parameters are developed for a cyclic beta-amino acid, amino-cyclo-pentane-carboxylic acid (ACPC), using a multi-objective evolutionary algorithm. The MM model is benchmarked using several short, ACPC-containing alpha/beta-peptides in water and methanol with SCC-DFTB (self consistent charge-density functional tight binding)/MM simulations as the reference. Satisfactory agreements are found between the MM and SCC-DFTB/MM results regarding the distribution of key dihedral angles for the tetra-alpha/beta-peptide in water. For the octa-alpha/beta-peptide in methanol, the MM and SCC-DFTB/MM simulations predict the 11- and 14/15-helical form as the more stable conformation, respectively; however, the two helical forms are very close in energy (2-4 kcal/mol) at both theoretical levels, which is also the conclusion from recent NMR experiments. As the first application, the MM model is applied to an alpha/beta-pentadeca-peptide in water with both explicit and implicit solvent models. The stability of the peptide is sensitive to the starting configuration in the explicit solvent simulations due to their limited length ( approximately 10-40 ns). Multiple ( approximately 20 x 20 ns) implicit solvent simulations consistently show that the 14/15-helix is the predominant conformation of this peptide, although substantially different conformations are also accessible. The calculated nuclear Overhauser effect (NOE) values averaged over different trajectories are consistent with experimental data, which emphasizes the importance of considering conformational heterogeneity in such comparisons for highly dynamical peptides.     

Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment

A pathogenetic feature of Alzhemier disease is the aggregation of monomeric beta-amyloid proteins (Abeta) to form oligomers. Usually these oligomers of long peptides aggregate on time scales of microseconds or longer, making computational studies using atomistic molecular dynamics models prohibitively expensive and making it essential to develop computational models that are cheaper and at the same time faithful to physical features of the process. We benchmark the ability of our implicit solvent model to describe equilibrium and dynamic properties of monomeric Abeta(10-35) using all-atom Langevin dynamics (LD) simulations, since Alphabeta(10-35) is the only fragment whose monomeric properties have been measured. The accuracy of the implicit solvent model is tested by comparing its predictions with experiment and with those from a new explicit water MD simulation, (performed using CHARMM and the TIP3P water model) which is approximately 200 times slower than the implicit water simulations. The dependence on force field is investigated by running multiple trajectories for Alphabeta(10-35) using the CHARMM, OPLS-aal, and GS-AMBER94 force fields, whereas the convergence to equilibrium is tested for each force field by beginning separate trajectories from the native NMR structure, a completely stretched structure, and from unfolded initial structures. The NMR order parameter, S2, is computed for each trajectory and is compared with experimental data to assess the best choice for treating aggregates of Alphabeta. The computed order parameters vary significantly with force field. Explicit and implicit solvent simulations using the CHARMM force fields display excellent agreement with each other and once again support the accuracy of the implicit solvent model. Alphabeta(10-35) exhibits great flexibility, consistent with experiment data for the monomer in solution, while maintaining a general strand-loop-strand motif with a solvent-exposed hydrophobic patch that is believed to be important for aggregation. Finally, equilibration of the peptide structure requires an implicit solvent LD simulation as long as 30 ns.     

Comparison of solvation-effect methods for the simulation of peptide interactions with a hydrophobic surface

In this study we investigated the interaction behavior between thirteen different small peptides and a hydrophobic surface using three progressively more complex methods of representing solvation effects: a united-atom implicit solvation method [CHARMM 19 force field (C19) with Analytical Continuum Electrostatics (ACE)], an all-atom implicit solvation method (C22 with GBMV), and an all-atom explicit solvation method (C22 with TIP3P). The adsorption behavior of each peptide was characterized by the calculation of the potential of mean force as a function of peptide-surface separation distance. The results from the C22/TIP3P model suggest that hydrophobic peptides exhibit relatively strong adsorption behavior, polar and positively-charged peptides exhibit negligible to relatively weak favorable interactions with the surface, and negatively-charged peptides strongly resist adsorption. Compared to the TIP3P model, the ACE and GBMV implicit solvent models predict much stronger attractions for the hydrophobic peptides as well as stronger repulsions for the negatively-charged peptides on the CH(3)-SAM surface. These comparisons provide a basis from which each of these implicit solvation methods may be reparameterized to provide closer agreement with explicitly represented solvation in simulations of peptide and protein adsorption to functionalized surfaces.     

Molecular dynamics study of peptides in implicit water: ab initio folding of beta-hairpin, beta-sheet, and beta beta alpha-motif

In this communication, we have demonstrated that molecular dynamics simulations using a GB implicit solvation model with the all-atom based force field (CHARMM19) can describe the spontaneous folding of small peptides in aqueous solution. The native structures of peptides with various structural motifs (beta-hairpin, beta-sheet, and betabetaalpha-moiety) were successfully predicted within reasonable time scales by MD simulations at moderately elevated temperatures. It is expected that the present simulations provide further insight into mechanism/pathways of the peptide folding.     

Hydrazino peptides as foldamers: an extension of the beta-peptide concept.

Replacing the C(beta) atoms in the beta-amino acid constituents of beta-peptides by nitrogen atoms leads to hydrazino peptides. A systematic conformation analysis of blocked hydrazino peptide oligomers of the general type I at the HF/6-31G, MP2/6-31G, and DFT/B3LYP/6-31G levels of ab initio MO theory and on the basis of molecular mechanics reveals a wide variety of secondary structures, as for instance various helices and sheet- and turnlike conformers. Some of them are closely related to secondary structure types found in beta-peptides; others represent novel types. Thus, a very stable, novel helix with 14-membered hydrogen-bonded pseudocycles, which occupies a conformation space different from that of helices with 14-membered rings found among the most stable conformers in beta-peptides, is indicated. The most important secondary structure elements are characterized by interactions between peptidic NH and CO groups. The additional hydrazino N(alpha)H group takes part in special structuring effects but is of lesser importance for secondary structure formation. The influence of environmental effects on the existence and stability of the various structure types is discussed. Due to the wide variety of structural possibilities, hydrazino peptides might be a useful tool for peptide and protein design.     

Implementation and application of helix-helix distance and crossing angle restraint potentials

Based on the definition of helix-helix distance and crossing angle introduced by Chothia et al. (J Mol Biol 1981, 145, 215), we have developed the restraint potentials by which the distance and crossing angle of two selected helices can be maintained around target values during molecular dynamics simulations. A series of assessments show that calculated restraint forces are numerically accurate. Since the restraint forces are only exerted on atoms which define the helical principal axes, each helix can rotate along its helical axis, depending on the helix-helix intermolecular interactions. Such a restraint potential enables us to characterize the helix-helix interactions at atomic details by sampling their conformational space around specific distance and crossing angle with (restraint) force-dependent fluctuations. Its efficacy is illustrated by calculating the potential of mean force as a function of helix-helix distance between two transmembrane helical peptides in an implicit membrane model.     

Coarse-grained protein molecular dynamics simulations

A limiting factor in biological science is the time-scale gap between experimental and computational trajectories. At this point, all-atom explicit solvent molecular dynamics (MD) are clearly too expensive to explore long-range protein motions and extract accurate thermodynamics of proteins in isolated or multimeric forms. To reach the appropriate time scale, we must then resort to coarse graining. Here we couple the coarse-grained OPEP model, which has already been used with activated methods, to MD simulations. Two test cases are studied: the stability of three proteins around their experimental structures and the aggregation mechanisms of the Alzheimer's Abeta16-22 peptides. We find that coarse-grained isolated proteins are stable at room temperature within 50 ns time scale. Based on two 220 ns trajectories starting from disordered chains, we find that four Abeta16-22 peptides can form a three-stranded beta sheet. We also demonstrate that the reptation move of one chain over the others, first observed using the activation-relaxation technique, is a kinetically important mechanism during aggregation. These results show that MD-OPEP is a particularly appropriate tool to study qualitatively the dynamics of long biological processes and the thermodynamics of molecular assemblies.     

Site-specific hydrogen-bonding interaction between N-acetylproline amide and protic solvent molecules: comparisons of IR and VCD measurements with MD simulations

The effects of solute-solvent interactions on solution structures of small peptides have been paid a great deal of attention. To study the effect of hydrogen-bonding interactions on peptide solution structures, we measured the amide I IR and VCD spectra of N-acetylproline amide (AP) in various protic solvents, i.e., D2O, MeOD, EtOD, and PrOD, and directly compared them with theoretically simulated ones. The numbers of protic solvent molecules hydrogen-bonded to the two peptide bonds in the AP were quantitatively determined by carrying out the molecular dynamics (MD) simulations and then compared with the spectral analyses of the experimentally measured amide I bands. The two peptides in the AP have different propensities of forming H-bonds with protic solvent molecules, and the H-bond population distribution is found to be strongly site-specific and solvent-dependent. However, it is found that adoption of the polyproline II (PII) conformation by AP in protic solvents does not strongly depend on the hydrogen bond network-forming ability of protic solvents nor on the solvent polarity. We present a brief discussion on the validity as well as limitation of the currently available force field parameters used for the present MD simulation study.     

Role of ions on structure and stability of a synthetic gramicidin ion channel in solution. A molecular dynamics study

We performed a molecular dynamics (MD) simulation to the investigate structure and stability of a synthetic gramicidin-like peptide in solution with and without ions. The starting structures of the MD simulations were taken from two recently solved NMR structures of this peptide in isotropic solution, which forms stable monomers or dimers in the presence or absence of ions, respectively. The monomeric structure is channel-like and is assumed to be stabilized by the presence of two Cs(+) ions bound in the channel, each one close to one channel entrance. In our MD simulations, we observed how the Cs(+) ions bind in the channel formed by the monomeric gramicidin-like peptide using implicit solvent and explicit ions with a concentration of 2 M. MD simulations were performed with and without explicit ions but with an implicit solvent model defined by the generalized Born approximation, which was used to mimic the dielectric properties of the solvent and to speed up the computations.     

Effects of surface interactions on peptide aggregate morphology

The formation of peptide aggregates mediated by an attractive surface is investigated using replica exchange molecular dynamics simulations with a coarse-grained peptide representation. In the absence of a surface, the peptides exhibit a range of aggregate morphologies, including amorphous aggregates, β-barrels and multi-layered fibrils, depending on the chiral stiffness of the chain (a measure of its β-sheet propensity). In contrast, aggregate morphology in the presence of an attractive surface depends more on surface attraction than on peptide chain stiffness, with the surface favoring fibrillar structures. Peptide-peptide interactions couple to peptide-surface interactions cooperatively to affect the assembly process both qualitatively (in terms of aggregate morphology) and quantitatively (in terms of transition temperature and transition sharpness). The frequency of ordered fibrillar aggregates, the surface binding transition temperature, and the sharpness of the binding transition all increase with both surface attraction and chain stiffness.     

Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field

The efficient and accurate characterization of solvent effects is a key element in the theoretical and computational study of biological problems. Implicit solvent models, particularly generalized Born (GB) continuum electrostatics, have emerged as an attractive tool to study the structure and dynamics of biomolecules in various environments. Despite recent advances in this methodology, there remain limitations in the parametrization of many of these models. In the present work, we demonstrate that it is possible to achieve a balanced implicit solvent force field by further optimizing the input atomic radii in combination with adjusting the protein backbone torsional energetics. This parameter optimization is guided by the potentials of mean force (PMFs) between amino acid polar groups, calculated from explicit solvent free energy simulations, and by conformational equilibria of short peptides, obtained from extensive folding and unfolding replica exchange molecular dynamics (REX-MD) simulations. Through the application of this protocol, the delicate balance between the competing solvation forces and intramolecular forces appears to be better captured, and correct conformational equilibria for a range of both helical and beta-hairpin peptides are obtained. The same optimized force field also successfully folds both beta-hairpin trpzip2 and mini-protein Trp-Cage, indicating that it is quite robust. Such a balanced, physics-based force field will be highly applicable to a range of biological problems including protein folding and protein structural dynamics.     

The link between sequence and conformation in protein structures appears to be stereochemically established

In search of the link between sequence and conformation in protein structures, we perform molecular dynamics analysis of the effect of stereochemical mutation in end-protected octa-alanine Ac-Ala8-NHMe from poly-L to an alternating-L,D structure. The mutation has a dramatic effect, transforming the peptide from a condition of extreme sensitivity to one of extreme insensitivity to solvent. Examining the molecular folds of poly-L and alternating-L,D structure in atomistic detail, we find them to differ in the relationship between peptide dipolar interactions at the local and nonlocal levels, either conflicting or harmonious depending upon the chain stereochemistry. The stereochemical transformation of interpeptide electrostatics from a condition of conflict to one of harmony explains the long-standing puzzle of why poly-L and alternating-L,D peptides strongly differ in properties such as "stiffness" and solvent sensitivity. Furthermore, it is possible that poly-L stereochemistry is also the fulcrum of protein sensitivity to the effects of amino acid side-chain structures via dielectric arbitrations in interpeptide electrostatics. Indeed the evidence is accumulating that the amino acid side chains differing in alpha-helix and beta-sheet propensities also differ in their desolvating effects in the adjacent and nearest-neighbor peptides and thus possibly in the solvent screening of peptide dipolar interactions.     

Conformational behavior of comb-like polyelectrolytes in selective solvent: computer simulation study

In this work, we present molecular dynamics simulations of comb-like polyelectrolytes in selective solvent. The studied polymers have a neutral backbone and polyelectrolyte side chains. The solvent is poor for the backbone and the theta solvent for the side chains. The polymers are modeled on a coarse-grained level with implicit solvent. The simulations show that the comb-like polyelectrolytes tend to form intramolecular self-organized structures of the pearl necklace type. This type of conformational behavior has been predicted by Borisov and Zhulina (Borisov, O. V.; Zhulina, E. B. Macromolecules 2005, 38, 2506) for neutral comb-like copolymers in selective solvent. The present study shows that comb-like polyelectrolytes in selective solvent exhibit the same type of behavior; however, it can be controlled by one additional parameter, the degree of dissociation of the grafts. The local conformational characteristics are studied using the ensemble-averaged bond angle cosines as functions of monomer position in the chain, which reveal structural details invisible by other means.     

Potential of mean force of water–proton bath and molecular dynamic simulation of proteins at constant pH

An advanced implicit solvent model of water–proton bath for protein simulations at constant pH is presented. The implicit water–proton bath model approximates the potential of mean force of a protein in water solvent in a presence of hydrogen ions. Accurate and fast computational implementation of the implicit water–proton bath model is developed using the continuum electrostatic Poisson equation model for calculation of ionization equilibrium and the corrected MSR6 generalized Born model for calculation of the electrostatic atom–atom interactions and forces. Molecular dynamics (MD) method for protein simulation in the potential of mean force of water–proton bath is developed and tested on three proteins. The model allows to run MD simulations of proteins at constant pH, to calculate pH‐dependent properties and free energies of protein conformations. The obtained results indicate that the developed implicit model of water–proton bath provides an efficient way to study thermodynamics of biomolecular systems as a function of pH, pH‐dependent ionization‐conformation coupling, and proton transfer events. © 2012 Wiley Periodicals, Inc.