Peptide folding simulations

Statistical mechanics of protein folding with separable energy functions

We have initiated an entirely new approach to statistical mechanical models of strongly interacting systems where the configurational parameters and the potential energy function are both constructed so that the canonical partition function can be evaluated analytically. For a simplified model of proteins consisting of a single, fairly short polypeptide chain without cross-links, we can adjust the energy parameters to favor the experimentally determined native state of seven proteins having diverse types of folds. Then 497 test proteins are predicted to have stable native folds, even though they are also structurally diverse, and 480 of them have no significant sequence similarity to any of the training proteins.

Atomically detailed simulations of helix formation with the stochastic difference equation

An algorithm is described to compute approximate classical trajectories as a boundary value problem with an integration step in the arc length. High-frequency motions are filtered out when a large integration step is used, maintaining the stability of the algorithm. At the limit of high filtering (large steps), but still offering an accurate description of the continuous path, the trajectory approaches the steepest descent path (SDP). The steepest descent path is widely used as a reaction coordinate in chemical systems. At intermediate step sizes, some inertial motions remain, interpolating between reaction coordinates and exact classical trajectories. Numerical studies of spatial and energetic properties of meta-trajectories are carried out. Two systems are considered here: valine dipeptide and the folding of a small helical protein. Although thermodynamic properties of meta-trajectories are affected by the filtering, the ordering of events remains similar for substantial differences in trajectory resolution.

Nature of structural inhomogeneities on folding a helix and their influence on spectral measurements

Extensive conformational sampling and calculations of vibrational coupling provide a quantitative basis for the structurally inhomogeneous spectra of the amide unit in aqueous solutions containing folded and unfolded state distributions of helices. Replica exchange molecular dynamics simulations of the capped helical peptide, AA(AAKAA)(3)AAY, is carried out over a range of temperatures, where the system populates the folded and unfolded states. This sampling defines a set of ensembles that characterizes the conformational variability for configurations identified by their fraction of helical content. The effects of hydrogen bonding, both internal and external (with water), and the coupling between amide-I modes are computed as a function of temperature and helical content. End-to-end distance and coupling distributions are also computed. The solvent H-bonding, which is present at all temperatures, shifts the amide-I band toward lower frequency compared with the unsolvated band. Upon thermal denaturation of the peptide, the amide-I band shifts to higher frequency because the increase in solvent H-bonding fails to compensate for the loss in internal (helical) H-bonds. The loss of uniformity of the mode coupling along the helix at higher temperatures accounts for the well-known thermal broadening of the amide IR spectrum. The calculated inhomogeneities of segments of the peptide predict experimental properties of isotope-edited helices.

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

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

Temperature dependence of the hydrophobic hydration and interaction of simple solutes: An examination of five popular water models

We examine five different popular rigid water models (SPC, SPCE, TIP3P, TIP4P, and TIP5P) using molecular dynamics simulations in order to investigate the hydrophobic hydration and interaction of apolar Lennard-Jones solutes as a function of temperature in the range between 275 and 375 K along the 0.1 MPa isobar. For all investigated models and state points we calculate the excess chemical potential for the noble gases and methane employing the Widom particle insertion technique. All water models exhibit too small hydration entropies, but show a clear hierarchy. TIP3P shows poorest agreement with experiment, whereas TIP5P is closest to the experimental data at lower temperatures and SPCE is closest at higher temperatures. As a first approximation, this behavior can be rationalized as a temperature shift with respect to the solvation behavior found in real water. A rescaling procedure inspired by the information theory model of Hummer et al. [Chem. Phys. 258, 349 (2000)] suggests that the different solubility curves for the different models and real water can be largely explained on the basis of the different density curves at constant pressure. In addition, the models that give a good representation of the water structure at ambient conditions (TIP5P, SPCE, and TIP4P) show considerably better agreement with the experimental data than the ones which exhibit less structured O-O correlation functions (SPC and TIP3P). In the second part of the paper we calculate the hydrophobic interaction between xenon particles directly from a series of 60 ns simulation runs. We find that the temperature dependence of the association is to a large extent related to the strength of the solvation entropy. Nevertheless, differences between the models seem to require a more detailed molecular picture. The TIP5P model shows by far the strongest temperature dependence. The suggested density rescaling is also applied to the chemical potential in the xenon-xenon contact-pair configuration, indicating the presence of a temperature where the hydrophobic interaction turns into purely repulsive. The predicted association for xenon in real water suggests the presence of a strong variation with temperature, comparable to the behavior found for TIP5P water. Comparing different water models and experimental data we conclude that a proper description of density effects is an important requirement for a water model to account correctly for the correct description of the hydrophobic effects. A water model exhibiting a density maximum at the correct temperature is desirable.

Temperature dependence of the thermodynamics of helix-coil transition

The temperature-induced helix to coil transition in a series of host peptides was monitored using circular dichroism spectroscopy (CD) and differential scanning calorimetry (DSC). Combination of these two techniques allowed direct determination of the enthalpy of helix-coil transition for the studied peptides. It was found that the enthalpy of the helix-coil transition differs for different peptides and this difference is related to the difference in the temperature for the midpoint of helix-coil transition. The enthalpy of the helix-coil transition decreases with the increase in temperature, thus providing the first experimental estimate for the heat capacity changes upon helix-coil transition, DeltaC(p). The values for DeltaC(p) of helix-coil transition are found to be negative, which is in contrast to the positive DeltaC(p) for protein unfolding. Analysis suggests that this negative DeltaC(p) of helix-coil transition is due to the exposure of the polar peptide backbone to solvent upon helix unfolding.