Aggregation of polyalanine in a hydrophobic environment

How Efficient Is Replica Exchange Molecular Dynamics? An Analytic Approach

Replica exchange molecular dynamics (REMD) has become a standard technique for accelerating relaxation in biosimulations. Despite its widespread use, questions remain about its efficiency compared with conventional, constant temperature molecular dynamics (MD). An analytic approach is taken to describe the relative efficiency of REMD with respect to MD. This is applied to several simple two-state models and to several real proteins-protein L and the B domain of protein A-to predict the relative efficiency of REMD with respect to MD in actual applications. In agreement with others, we find the following: as long as there is a positive activation energy for folding, REMD is more efficient than MD; the effectiveness of REMD is strongly dependent on the activation enthalpy; and the efficiency of REMD for actual proteins is a strong function of the maximum temperature. Choosing the maximum temperature too high can result in REMD becoming significantly less efficient than conventional MD. A good rule of thumb appears to be to choose the maximum temperature of the REMD simulation slightly above the temperature at which the enthalpy for folding vanishes. Additionally, we find that the number of replicas in REMD, while important for simulations shorter than one or two relaxation times, has a minimal effect on the asymptotic efficiency of the method.

Conformational behavior of polyalanine peptides with and without protecting groups of varying chain lengths: population of PP-II structure!

Oculopharyngeal muscular dystrophy (OPMD), a polyalanine myopathy, occurs due to expansion of homo-polyalanine stretch in normal polyadenylating binding protein nuclear 1 (PABPN1) protein from Ala10 to Ala11-17. Therefore, the conformational behavior of polyalanine peptides with n = 10-17, with and without terminal protecting groups, have been investigated with different starting geometries in water by molecular dynamics simulation studies. Alanine peptides are shown to give rise to unordered structure irrespective of starting geometry and not more than two residues at a stretch have the same/similar set of φ, ψ values. However, the final structure with terminal protecting groups look like β-strand. Unprotected poly-Ala peptides adopt twisted β-hairpin/multi hairpin like structure with increasing chain length. The number of residues having φ, ψ values in collagen region is found to be less in peptides with unprotected termini as compared to peptides with protected termini of same chain length. The results have been supported by recent synchrotron radiation circular dichroism spectroscopy of polyproline II and unordered secondary structures. Opening of the helical structure in poly-Ala peptides with protecting groups has been shown to take place from C-terminal and in peptides without protecting groups opening of helix starts from both terminals. Further, opening of helix takes more time in poly-Ala peptides without terminal protecting groups. The deviations in amide bond planarity have been discussed and compared with available experimental and computational results.

Involvement of non-polyalanine (polyA) residues in aggregation of polyA proteins: Clue for inhibition of aggregation

Presence of polyalanine (polyA) stretches in some proteins is found to be associated with their aggregation, which causes disorders in various developmental processes. In this work, inherent propensities towards aggregation of some residues, which are not part of the polyA stretches, have been identified by using the primary sequences of seven polyA proteins with the help of Betascan, PASTA and Tango programs and explored unambiguously. This provides a basis for proposing molecular mechanism of this type of aggregation. Reported suppression of aggregation of polyA proteins by chaperones like HSP40 and HSP70 is substantiated through molecular docking. The hydrophobic residues of identified aggregating region are found to be interacting with hydrophobic surface of chaperones. This suggests a crucial clue for possible way to inhibit the aggregation of such proteins.

Gay-Berne and electrostatic multipole based coarse-grain potential in implicit solvent

A general, transferable coarse-grain (CG) framework based on the Gay-Berne potential and electrostatic point multipole expansion is presented for polypeptide simulations. The solvent effect is described by the Generalized Kirkwood theory. The CG model is calibrated using the results of all-atom simulations of model compounds in solution. Instead of matching the overall effective forces produced by atomic models, the fundamental intermolecular forces such as electrostatic, repulsion-dispersion, and solvation are represented explicitly at a CG level. We demonstrate that the CG alanine dipeptide model is able to reproduce quantitatively the conformational energy of all-atom force fields in both gas and solution phases, including the electrostatic and solvation components. Replica exchange molecular dynamics and microsecond dynamic simulations of polyalanine of 5 and 12 residues reveal that the CG polyalanines fold into "alpha helix" and "beta sheet" structures. The 5-residue polyalanine displays a substantial increase in the "beta strand" fraction relative to the 12-residue polyalanine. The detailed conformational distribution is compared with those reported from recent all-atom simulations and experiments. The results suggest that the new coarse-graining approach presented in this study has the potential to offer both accuracy and efficiency for biomolecular modeling.

Molecular dynamics simulations of the folding of poly(alanine) peptides

The secondary structures and the shapes of long-chain polyalanine (PA) molecules were investigated by all-atom molecular dynamics simulations using a modified Amber force field. Homopolymers of polyaminoacids such as PA are convenient models to study the mechanism of protein folding. It was found that the conformational structures of PA peptides are highly sensitive to the chain length. In the absence of solvent, straight α-helices dominate in short (n ∼ 20) peptides at room temperature. A shape transition occurs at a chain length n of 40-45; the compact helix-turn-helix structure (the double-leg hairpin) becomes favored over a straight α-helix. For n=60, double-leg and the triple-leg hairpins are the only structures present in PA molecules. An exploration of a chain organization in a cubic cavity revealed a clear predisposition of PA molecules for additional breaks in α-helices and the formation of multifolded hairpins. Furthermore, under confinement the hairpin structure becomes much looser, the antiparallel positions of helical stems are disturbed, and a sizeable proportion of the helical stems are transformed from α-helices into 3(10)-helices.

Self-assembly of polypeptides into left-handedly twisted fibril-like structures

In this paper, we investigated the spontaneous formation of aggregation structures of amyloid-forming peptide (GGVVIA) using a coarse-grained model and Monte Carlo simulations. The effects of concentration and temperature on the formation of different aggregation structures were studied. Three types of aggregation structures, single-layer beta sheet, amorphous beta-sheet aggregate, and fibril-like structures, were observed in our simulations. The fibril-like structures obtained in simulations have a common cross-beta spine structure in which beta sheets twist in a left-handed fashion. The averaged twisting angle of the beta sheet in the fibril-like structures is 12 degrees +/-2 degrees. Moreover, it was found that the peptides in the same beta sheets prefer to arrange in a parallel way, which is consistent with the corresponding GGVVIA crystalline structure. On the other hand, it was found that there is a rich family of beta-sheet stacking patterns in the fibril-like structures suggesting that the fibril structures are more complex than the corresponding crystalline structure and there exist many local free-energy minima rather than a distinct global minimum.

Association of helical beta-peptides and their aggregation behavior from the potential of mean force in explicit solvent

Helical beta-peptides have been shown to fold into well-defined structures. In aqueous solution, some beta-peptides self-assemble into nanoscale fibers, aggregates, and liquid crystalline phases. Molecular simulations, at an atomistic level, are used to examine, in a systematic manner, the interactions between distinct beta-peptide molecules. The relationship between side-chain chemistry (and position along the backbone) and, in particular, aggregation behaviors, is assessed by calculating the potential of mean force or dimerization free energy of two peptides in explicit water. The free energy profiles as a function of separation for helical, amphiphilic beta-peptides are consistent with experimental observations, and help explain the origins of aggregate or fiber formation in solution. Close examination of the energetic and entropic contributions to the free energy reveals that, depending on the position of certain side groups along the molecule, the tendency of two peptides to aggregate can be driven by entropy or by energy, respectively. In contrast to findings from previous works that employed a coarse representation of the solvent, it is shown that water-peptide interactions play key roles in the association behavior of beta-peptides.

Computational approaches for the design of peptides with anti-breast cancer properties

Background: Breast cancer is the most common cancer among women. Tamoxifen is the preferred drug for estrogen receptor-positive breast cancer treatment, yet many of these cancers are intrinsically resistant to tamoxifen or acquire resistance during treatment. Therefore, scientists are searching for breast cancer drugs that have different molecular targets. Methodology: Recently, a computational approach was used to successfully design peptides that are new lead compounds against breast cancer. We used replica exchange molecular dynamics to predict the structure and dynamics of active peptides, leading to the discovery of smaller bioactive peptides. Conclusions: These analogs inhibit estrogen-dependent cell growth in a mouse uterine growth assay, a test showing reliable correlation with human breast cancer inhibition. We outline the computational methods that were tried and used along with the experimental information that led to the successful completion of this research.

Helix forming tendency of valine substituted poly-alanine: a molecular dynamics investigation

In this study, classical molecular dynamics simulations have been carried out on the valine (guest) substituted poly alanine (host) using the host-guest peptide approach to understand the role of valine in the formation and stabilization of helix. Valine has been substituted in the host peptide starting from N terminal to C terminal. Various structural parameters have been obtained from the molecular dynamics simulation to understand the tolerance of helical motif to valine. Depending on the position of valine in the host peptide, it stabilizes (or destabilizes) the formation of the helical structure. The substitution of valine in the poly alanine at some positions has no effect on the helix formation (deformation). It is interesting to observe the coexistence of 3 10 and alpha-helix in the peptides due to the dynamical nature of the hydrogen bonding interaction and sterical interactions.

Modulation of self-association and subsequent fibril formation in an alanine-rich helical polypeptide

Thermal unfolding, reversible self-association, and irreversible aggregation were investigated for an alanine-rich helical polypeptide, 17-H-6, with sequence [AAAQEAAAAQAAAQAEAAQAAQ] 6. Dynamic light scattering, transmission electron microscopy, and thermal unfolding measurements indicate that 17-H-6 spontaneously and reversibly self-associates at acidic pH and low temperature. The resulting multimers have a compact, globular morphology with an average hydrodynamic radius approximately 10-20 nm and reversibly dissociate to monomers upon an increase to pH 7.4. Both free monomer and 17-H-6 chains within the multimers are alpha-helical and folded at low temperature. Reversible unfolding of the monomer occurs upon heating of solutions at pH 7.4. At pH 2.3, heating first causes incomplete dissociation and unfolding of the constituent chains. Further incubation at elevated temperature induces additional structural and morphological changes and results in fibrils with a beta-sheet 2 degrees structure and a characteristic diameter of 5-10 nm (7 nm mean). The ability to modulate association and aggregation suggests opportunities for this class of polypeptides in nanotechnology and biomedical applications.

Computational design and experimental discovery of an antiestrogenic peptide derived from alpha-fetoprotein

Breast cancer is the most common cancer among women, and tamoxifen is the preferred drug for estrogen receptor-positive breast cancer treatment. Many of these cancers are intrinsically resistant to tamoxifen or acquire resistance during treatment. Consequently, there is an ongoing need for breast cancer drugs that have different molecular targets. Previous work has shown that 8-mer and cyclic 9-mer peptides inhibit breast cancer in mouse and rat models, interacting with an unsolved receptor, while peptides smaller than eight amino acids did not. We show that the use of replica exchange molecular dynamics predicts the structure and dynamics of active peptides, leading to the discovery of smaller peptides with full biological activity. Simulations identified smaller peptide analogues with the same conserved reverse turn demonstrated in the larger peptides. These analogues were synthesized and shown to inhibit estrogen-dependent cell growth in a mouse uterine growth assay, a test showing reliable correlation with human breast cancer inhibition.

Two-Component Polypeptides Modeled with Effective Pair Potentials

We present Monte Carlo simulations performed within a model based on a set of distance-dependent effective potentials which are used to describe the interactions between a pair of distinct amino acids. These effective potentials are extracted from experimental correlation functions through the Ornstein-Zernike equations and adequate closure approximations. We focus our attention on the sequences of two specific residues, namely, alanine and glycine. The studied sequences are (a) (Ala)(12)-(Gly)(4)-(Ala)(12) and (b) three interacting chains of alternating alanines and glycines (with five residues per chain). The resulting structures are combinations of known secondary structures. More importantly, we verify that our simulated structures are in thermodynamic equilibrium by means of an estimation of the density of states.