Molecular dynamics simulation of α-melanocyte stimulating hormone in a water-membrane model interface

The conformation of the tridecapeptide α-melanocyte stimulating hormone in the presence of a double water-membrane interface was studied by molecular dynamics simulation, using the computational package THOR. In this program the solvent is represented by a continuous medium with dielectric constant ɛ, and the interface between different media is simulated by a surface of discontinuity of the dielectric constant. The electrostatic image method was used to write down the terms, added to the force field, that describe the polarisation effects induced in the interface by the atomic charges. The program was further improved by the introduction of a second surface, parallel to the first one, to mimic the membrane. A conformational search using the software Prelude was employed to find an initial geometry for the peptide in water. The molecular dynamics simulation performed during 10 ns showed that the peptide structure is flexible in water, without stabilisation of any preferential conformation. In the presence of the model membrane, the peptide moved to the medium representing the interior of the membrane. Inside the low dielectric constant medium, the structure of the peptide showed a turn in the central sequence of amino acids and a packed conformation remained stabilised during more than 7.0 ns of simulation. Received: 27 November 1998 / Revised version: 11 March 1999 / Accepted: 8 April 1999     

Conformational study of a nine residue fragment of the antigenic loop of foot-and-mouth disease virus.

The nine-residue peptide Ac-TASARGDLA-NHMe was selected as model peptide in order to understand the conformational features of the antigenic loop of foot-and-mouth disease virus (FMDV). A throughout exploration of the conformational space has been carried out by means of molecular dynamics (MD) and energy minimization. The calculations have been carried out using the AMBER force field. Solvent effects have been included by an effective dielectric constant of epsilon = 4r. The lowest energy conformation presents a secondary structure constituted by an alpha-helix at the N-terminal end followed by two gamma-turns in the central region. The rest of the accessible minima found present also a high tendency to form gamma-turns. Finally, a 100 ps MD trajectory calculation at 298 K suggest a stability of the secondary structure elements of the lowest energy conformation.     

Simulation of molecular crowding effects on an Alzheimer's α-amyloid peptide

Fibril formation by the Alzheimer's β-amyloid (Aβ) peptide in brain tissue is integral to the Alzheimer's disease pathology. Understanding the conformational properties and the mechanisms triggering aggregation of the Aβ peptides, at an atomic level of detail, is of crucial importance for the design of effective therapeutic agents against this disease. In this work, the conformational transitions and dynamic properties of an amyloidogenic peptide fragment (Aβ_10-35) were studied by molecular dynamics simulations in systems modeling infinite dilution and the presence of macromolecular crowding agents (CA). The model system consists of the peptide described with an atomistic force field, the CA represented by inert, quasi-hard spheres and a continuum solvent model. This combined model allowed the simulations to be extended to 100 ns each. Simulations were carried out starting from a completely extended structure, a β-strand structure, and four nuclear magnetic resonance structures in dilute aqueous solution. For all structures, two additional simulations were performed that included the inert CA in the solution and occupied approx 30 and 40% of the volume, respectively. For two of the nuclear magnetic resonance structures, additional simulations were carried out with 35% volume fraction of CA to further examine the diffusive behavior of the peptide. The peptide adopted a collapsed coil conformation in all simulations. The results of the simulations in dilute solution showed reasonable qualitative agreement with experimental and other simulation results, whereas the presence of volume excluding agents resulted in some distinct changes in properties (e.g., an increase in the appearance of transient β-structure or decreases in diffusivity with increasing CA concentration). At the same time, internal motion such as order parameter or atomic root mean square fluctuations showed less systematic responses to volume exclusion.     

Accurate prediction of protein dihedral angles through conditional random field

Identifying local conformational changes induced by subtle differences on amino acid sequences is critical in exploring the functional variations of the proteins. In this study, we designed a computational scheme to predict the dihedral angle variations for different amino acid sequences by using conditional random field. This computational tool achieved an accuracy of 87% and 84% in 10-fold cross validation in a large data set for φ and Ψ, respectively. The prediction accuracies of φ and Ψ are positively correlated to each other for most of the 20 types of amino acids. Helical amino acids can achieve higher prediction accuracy in general, while amino acids in beet sheet have higher accuracy at specific angular regions. The prediction accuracy of φ is negatively correlated with amino acid flexibility represented by Vihinen Index. The prediction accuracy of φ can also be negatively correlated with angle distribution dispersion.     

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.     

Modeling of the structure in aqueous solution of the exopolysaccharide produced by Lactobacillus helveticus 766.

A method is described for constructing a conformational model in water of a heteropolysaccharide built up from repeating units, and is applied to the exopolysaccharide produced by Lactobacillus helveticus 766. The molecular modeling method is based on energy minima, obtained from molecular mechanics calculations of each of the constituting disaccharide fragments of the repeating unit in vacuo, as starting points. Subsequently, adaptive umbrella sampling of the potential of mean force is applied to extract rotamer populations of glycosidic dihedral angles of oligosaccharide fragments in solution. From these analyses, the most probable conformations are constructed for the hexasaccharide-repeating unit of the polysaccharide. After exploring the conformational space of each of the individual structures by molecular dynamics simulations, the different repeating unit conformations are used as building blocks for the generation of oligo- and polysaccharide models, by using a polysaccharide building program. The created models of the exopolysaccharide produced by L. helveticus 766 exhibit a flexible twisted secondary structure and tend to adopt a random coil conformation as tertiary structure.     

Modeling of loops in protein structures

Comparative protein structure prediction is limited mostly by the errors in alignment and loop modeling. We describe here a new automated modeling technique that significantly improves the accuracy of loop predictions in protein structures. The positions of all nonhydrogen atoms of the loop are optimized in a fixed environment with respect to a pseudo energy function. The energy is a sum of many spatial restraints that include the bond length, bond angle, and improper dihedral angle terms from the CHARMM-22 force field, statistical preferences for the main-chain and side-chain dihedral angles, and statistical preferences for nonbonded atomic contacts that depend on the two atom types, their distance through space, and separation in sequence. The energy function is optimized with the method of conjugate gradients combined with molecular dynamics and simulated annealing. Typically, the predicted loop conformation corresponds to the lowest energy conformation among 500 independent optimizations. Predictions were made for 40 loops of known structure at each length from 1 to 14 residues. The accuracy of loop predictions is evaluated as a function of thoroughness of conformational sampling, loop length, and structural properties of native loops. When accuracy is measured by local superposition of the model on the native loop, 100, 90, and 30% of 4-, 8-, and 12-residue loop predictions, respectively, had <2 A RMSD error for the mainchain N, C(alpha), C, and O atoms; the average accuracies were 0.59 +/- 0.05, 1.16 +/- 0.10, and 2.61 +/- 0.16 A, respectively. To simulate real comparative modeling problems, the method was also evaluated by predicting loops of known structure in only approximately correct environments with errors typical of comparative modeling without misalignment. When the RMSD distortion of the main-chain stem atoms is 2.5 A, the average loop prediction error increased by 180, 25, and 3% for 4-, 8-, and 12-residue loops, respectively. The accuracy of the lowest energy prediction for a given loop can be estimated from the structural variability among a number of low energy predictions. The relative value of the present method is gauged by (1) comparing it with one of the most successful previously described methods, and (2) describing its accuracy in recent blind predictions of protein structure. Finally, it is shown that the average accuracy of prediction is limited primarily by the accuracy of the energy function rather than by the extent of conformational sampling.     

Protein Structure Prediction with a Combined Solvation Free Energy-Molecular Mechanics Force Field

Abstract

Models of protein structure are frequently used to determine the physical characteristics of a protein when the crystal structure is not available. We developed a procedure to optimize such models, by use of a combined solvation free energy and molecular mechanics force field. Appropriately chosen atomic solvation parameters were defined using the criterion that the resulting protein model should deviate least from the crystal structure upon a forty picosecond molecular dynamics simulation carried out using the combined force field. Several tests were performed to refine the set of atomic solvation parameters which best complement the molecular mechanics forces. Four sets of parameters from the literature were tested and an empirically optimized set is proposed. The parameters are defined on a well characterized small molecule (alanyl dipeptide) and on the highly refined crystal structure of rat trypsin, and then tested on a second highly refined crystal structure of α-lytic protease. The new set of atomic solvation parameters provides a significant improvement over molecular mechanics alone in energy minimization of protein structures. This combined force field also has advantages over the use of explicit solvent as it is possible to take solvent effects into account during energetic conformational searching when modeling a homologous protein structure from a known crystal structure.     

Effects of pH and temperature on the structural and thermodynamic character of alpha-syn12 peptide in aqueous solution

The structural and thermodynamic characters of alpha-syn12 peptide in aqueous solution at different pH and temperatures have been investigated through temperature replica exchange molecular dynamics (T-REMD) simulations with GROMOS 43A1 force field. The two independent T-REMD simulations were completed at pH = 7.0 and 10.0, respectively. Each replica was run for 300 ns. The structural and thermodynamic characters of alpha syn12 peptide were studied based on the distributions of backbone dihedral angles, the free energy surface, and the stability of different type structure and the favorite conformations of the peptide. The results showed that the simulation at pH = 10.0 produced more sampling in alpha region than the simulation at pH = 7.0. The temperature changes from 283 K to 308 K result in negligible effects on the distributions of backbone dihedral angle. The beta hairpin conformation with Turn(9-6) and four hydrogen bonds (HB(4-11), HB(6-9), HB(9-6) and HB(11-4)) is the lowest free energy state in the simulation at pH = 7.0. However, for the simulation at pH = 10.0, the lowest free energy state corresponds to a structure with Turn(9-6) and two hydrogen bonds (HB(6-10) and HB(10-6)) induced by an overly strong residue-residue interaction effect between lysine residues. For the simulation at pH = 7.0, the free energy change of the alpha-syn12 peptide from the unfolded state to the beta hairpin state was in good agreement with the experiments and molecular dynamics simulation results for the other beta-peptides, the beta hairpin state of the alpha-syn12 peptide included the conformations that not only the Turn(9-6) is formed, but also the terminus are closed together in space. However, the subtle balances between lysine-lysine interactions and lysine-solvent interaction are disrupted in the simulation at pH = 10.0, which induced the assembly of lysine residues, the beta hairpin conformation is destabilized by the deprotonation of the Lys side chain. This study can help us to understand the conformation changes and the thermodynamic character of alpha;-syn12 peptide at atomic level induced by changing pH and temperature, which is propitious to reveal the nosogenesis of Parkinson disease. In our knowledge, this is the first report to study the influence of pH and temperature on isolated alpha-syn12 peptide in water by T-REMD.     

Interplay of charge distribution and conformation in peptides: comparison of theory and experiment

We assessed the correlation between charge distribution and conformation of flexible peptides by comparing the theoretically calculated potentiometric-titration curves of two model peptides, Ac-Lys5-NHMe (a model of poly-L-lysine) and Ac-Lys-Ala11-Lys-Gly2-Tyr-NH2 (P1) in water and methanol, with the experimental curves. The calculation procedure consisted of three steps: (i) global conformational search of the peptide under study using the electrostatically driven Monte Carlo (EDMC) method with the empirical conformational energy program for peptides (ECEPP)/3 force field plus the surface-hydration (SRFOPT) or the generalized Born surface area (GBSA) solvation model as well as a molecular dynamics method with the assisted model building and energy refinement (AMBER)99/GBSA force field; (ii) reevaluation of the energy in the pH range considered by using the modified Poisson-Boltzmann approach and taking into account all possible protonation microstates of each conformation, and (iii) calculation of the average degree of protonation of the peptide at a given pH value by Boltzmann averaging over conformations. For Ac-Lys5-NHMe, the computed titration curve agrees qualitatively with the experimental curve of poly-L-lysine in 95% methanol. The experimental titration curves of peptide P1 in water and methanol indicate a remarkable downshift of the first pK(a) value compared to the values for reference compounds (n-butylamine and phenol, respectively), suggesting the presence of a hydrogen bond between the tyrosine hydroxyl oxygen and the H(epsilon) proton of a protonated lysine side chain. The theoretical titration curves agree well with the experimental curves, if conformations with such hydrogen bonds constitute a significant part of the ensemble; otherwise, the theory predicts too small a downward pH shift.     

On the Bioactive Conformation of a Small Peptide and its Set of Thermodynamically Accessible Conformations

Abstract

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics framework using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.     

The atomic structure of crystalline porcine pancreatic elastase at 2.5 A resolution: comparisons with the structure of alpha-chymotrypsin

Three isomorphous heavy-atom derivatives have been used to calculate a 2.5 Å resolution electron density map of tosyl-elastase at pH 5.0, from which an accurate atomic model has been constructed. Atomic co-ordinates measured from this model have been refined using model building, real-space refinement and energy minimization programs. The three-dimensional conformation of the polypeptide chain is described in terms of conformational angles, hydrogen-bonding networks and the environment of different types of amino acid side-chain.Difference Fourier calculation of the high resolution structure of native elastase at pH 5.0 shows it to be virtually identical to that of the tosyl derivative, except near the tosyl group. The conformation of the catalytically important residues in native elastase is very similar to that of native α-chymotrypsin, except for the orientation of the active centre serine oxygen. The significance of important structural similarities and differences between these two enzymes is discussed.Elastase contains 25 internal water molecules which play an important role in stabilizing the active conformation of the enzyme. Many of these water molecules are in identical positions to those found in the interior of α-chymotrypsin     

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

Molecular dynamics method for proteins with ionization-conformation coupling and equilibrium titration

A new realization of the constant pH molecular dynamics simulation is proposed in mean force potential of proton reservoir in constant titration conditions. The MD-pH-ET method is used for a protein in the most probable ionization microstate, which is the state with minimal energy of the current conformation taking into account the correcting ionization potential of mean force considering an equilibrium ensemble of ionization states. The new MD-pH-ET method allows one to carry out an optimization of protein structure and the total free energy of a protein in the aqueous solution at constant pH and to calculate the pH-dependent properties. The MD-pH-ET method possesses the following unique features: (1) it uses the most precise and computationally effective realization of calculation of electrostatic energy of a protein in aqueous solution, the model of continuous dielectric media with Poisson equation, and the generalized Born method with “ideal” Born atomic radii; (2) it uses the same model of the potential energy surface in the ionization-conformation phase space both for calculating the potential energy of the protein and atomic forces and for determining the most probable ionization states; (3) it calculates the total free energy of the protein in the aqueous solution in a proton reservoir under the conditions of equilibrium titration. The workability of the new method MD-pH-ET is demonstrated for a molecule of the BPTI protein.     

Flap opening dynamics in HIV-1 protease explored with a coarse-grained model

We present a one-bead coarse-grained model that enables dynamical simulations of proteins on the time scale of tens of microseconds. The parameterization of the force field includes accurate conformational terms that allow for fast and reliable exploration of the configurational space. The model is applied to the dynamics of flap opening in HIV-1 protease. The experimental structure of the recently crystallized semi-open conformation of HIV-1 protease is well reproduced in the simulation, which supports the accuracy of our model. Thanks to very long simulations and extensive sampling of opening and closing events, we also investigate the thermodynamics and kinetics of the opening process. We have shown that the effect of the solvent slows down the dynamics to the experimentally observed time scales. The model is found to be reliable for application to substrate docking simulations, which are currently in progress.     

A small tripeptide AFA undergoes two state cooperative conformational transitions: implications for conformational biases in unfolded states

It is important to understand the conformational biases that are present in unfolded states to understand protein folding. In this context, it is surprising that even a short tripeptide like AFA samples folded/ordered conformation as demonstrated recently by NMR experiments of the peptide in aqueous solution at 280 K. In this paper, we present molecular dynamics simulation of the peptide in explicit water using OPLS-AA/L all-atom force field. The results are in overall agreement with NMR results and provide some further insights. The peptide samples turn and extended conformational forms corresponding to minima in free energy landscape. Frequent transitions between the minima are observed due to modest free energy barriers. The turn conformation seems to be stabilized by hydrophobic interactions and possibly by bridging water molecules between backbone donors and acceptors. Thus the peptide does not sample conformations randomly, but samples well defined conformations. The peptide served as a model for folding-unfolding equilibrium in the context of peptide folding. Further, implications for drug design are also discussed.     

Controlling the secondary-structure-forming tendencies of proteins by a backbone torsion-energy term

We examined a new backbone torsion-energy term proposed by us in the force field for protein systems. This torsion-energy term is represented by a double Fourier series in two variables, namely the backbone dihedral angles φ and ψ. It gives a natural representation of the torsion energy in the Ramachandran space in the sense that any two-dimensional energy surface periodic in both φ and ψ can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by modifying the torsion-energy surface. For instance, we can increase or decrease the α-helix-forming-tendencies by lowering or raising the torsion-energy surface in the α-helix region and likewise increase or decrease the β-sheet-forming tendencies by lowering or raising the surface in the β-sheet region in the Ramachandran space. We applied this torsion-energy modification method to six force fields, AMBER parm94, AMBER parm96, AMBER parm99, CHARMM27, OPLS-AA and OPLS-AA/L, and demonstrated that our modifications of the torsion-energy terms resulted in the expected changes of secondary-structure-forming tendencies by performing folding simulations of α-helical and β-hairpin peptides.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

Conformation of peptides in lipid membranes studied by x-ray grazing incidence scattering

Although the antimicrobial, fungal peptide alamethicin has been extensively studied, the conformation of the peptide and the interaction with lipid bilayers as well as the mechanism of channel gating are still not completely clear. As opposed to studies of the crystalline state, the polypeptide structures in the environment of fluid bilayers are difficult to probe. We have investigated the conformation of alamethicin in highly aligned stacks of model lipid membranes by synchrotron-based x-ray scattering. The (wide-angle) scattering distribution has been measured by reciprocal space mappings. A pronounced scattering signal is observed in samples of high molar peptide/lipid ratio which is distinctly different from the scattering distribution of an ideal helix in the transmembrane state. Beyond simple models of ideal helices, the data is analyzed in terms of models based on atomic coordinates from the Brookhaven Protein Data Bank, as well as from published molecular dynamics simulations. The results can be explained by assuming a wide distribution of helix tilt angles with respect to the membrane normal and a partial insertion of the N-terminus into the membrane.