Monte Carlo Studies of Folding, Dynamics, and Stability in α-Helices

Folding simulations of polyalanine peptides were carried out using an off-lattice Monte Carlo simulation technique. The peptide was represented as a chain of residues, each of which contains two interaction sites: one corresponding to the C(α) atom and the other to the side chain. A statistical potential was used to describe the interaction between these sites. The preferred conformations of the peptide chain on the energy surface, starting from several initial conditions, were searched by perturbations on its generalized coordinates with the Metropolis criterion. We observed that, at low temperatures, the effective energy was low and the helix content high. The calculated helix propagation (s) and nucleation (σ) parameters of the Zimm-Bragg model were in reasonable agreement with the empirical data. Exploration of the energy surface of the alanine-based peptides (AAQAA)(3) and AAAAA(AAARA)(3)A demonstrated that their behavior is similar to that of polyalanine, in regard to their effective energy, helix content, and the temperature-dependence of their helicity. In contrast, stable secondary structures were not observed for (Gly)(20) at similar temperatures, which is consistent with the nonfolder nature of this peptide. The fluctuations in the slowest dynamics mode, which describe the elastic behavior of the chain, showed that as the temperature decreases, the polyalanine peptides become stiffer and retain conformations with higher helix content. Clustering of conformations during the folding phase implied that polyalanine folds into a helix through fewer numbers of intermediate conformations as the temperature decreases

Interactions of Hydrophobic Peptides with Lipid Bilayers: Monte Carlo Simulations with M2δ

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone α-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2δ and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2δ were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2δ is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14° between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 Å compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides

Glycosylation May Reduce Protein Thermodynamic Stability by Inducing a Conformational Distortion

Glycosylation plays not only a functional role but can also modify the biophysical properties of the modified protein. Usually, natural glycosylation results in protein stabilization; however, in vitro and in silico studies showed that sometimes glycosylation results in thermodynamic destabilization. Here, we applied coarse-grained and all-atom molecular dynamics simulations to understand the mechanism underlying the loss of stability of the MM1 protein by glycosylation. We show that the origin of the destabilization is a conformational distortion of the protein caused by the interaction of the monosaccharide with the protein surface. Though glycosylation creates new short-range glycan–protein interactions that stabilize the conjugated protein, it breaks long-range protein–protein interactions. This has a destabilizing effect because the probability of long- and short-range interactions forming differs between the folded and unfolded states. The destabilization originates not from simple loss of interactions but due to a trade-off between the short- and long-range interactions

Interactions of Cationic-Hydrophobic Peptides with Lipid Bilayers: A Monte Carlo Simulation Method

We present a computational model of the interaction between hydrophobic cations, such as the antimicrobial peptide, Magainin2, and membranes that include anionic lipids. The peptide's amino acids were represented as two interaction sites: one corresponds to the backbone α-carbon and the other to the side chain. The membrane was represented as a hydrophobic profile, and its anionic nature was represented by a surface of smeared charges. Thus, the Coulombic interactions between the peptide and the membrane were calculated using the Gouy-Chapman theory that describes the electrostatic potential in the aqueous phase near the membrane. Peptide conformations and locations near the membrane, and changes in the membrane width, were sampled at random, using the Metropolis criterion, taking into account the underlying energetics. Simulations of the interactions of heptalysine and the hydrophobic-cationic peptide, Magainin2, with acidic membranes were used to calibrate the model. The calibrated model reproduced structural data and the membrane-association free energies that were measured also for other basic and hydrophobic-cationic peptides. Interestingly, amphipathic peptides, such as Magainin2, were found to adopt two main membrane-associated states. In the first, the peptide resided mostly outside the polar headgroups region. In the second, which was energetically more favorable, the peptide assumed an amphipathic-helix conformation, where its hydrophobic face was immersed in the hydrocarbon region of the membrane and the charged residues were in contact with the surface of smeared charges. This dual behavior provides a molecular interpretation of the available experimental data

Context-Dependent Effects of Asparagine Glycosylation on Pin WW Folding Kinetics and Thermodynamics

Context-Dependent Effects of Asparagine Glycosylation on Pin WW Folding Kinetics and Thermodynamic

Criteria for Selecting PEGylation Sites on Proteins for Higher Thermodynamic and Proteolytic Stability

PEGylation of protein side chains has been used for more than 30 years to enhance the pharmacokinetic properties of protein drugs. However, there are no structure- or sequence-based guidelines for selecting sites that provide optimal PEG-based pharmacokinetic enhancement with minimal losses to biological activity. We hypothesize that globally optimal PEGylation sites are characterized by the ability of the PEG oligomer to increase protein conformational stability; however, the current understanding of how PEG influences the conformational stability of proteins is incomplete. Here we use the WW domain of the human protein Pin 1 (WW) as a model system to probe the impact of PEG on protein conformational stability. Using a combination of experimental and theoretical approaches, we develop a structure-based method for predicting which sites within WW are most likely to experience PEG-based stabilization, and we show that this method correctly predicts the location of a stabilizing PEGylation site within the chicken Src SH3 domain. PEG-based stabilization in WW is associated with enhanced resistance to proteolysis, is entropic in origin, and likely involves disruption by PEG of the network of hydrogen-bound solvent molecules that surround the protein. Our results highlight the possibility of using modern site-specific PEGylation techniques to install PEG oligomers at predetermined locations where PEG will provide optimal increases in conformational and proteolytic stability