Conformational Changes of Trialanine Induced by Direct Interactions between Alanine Residues and Alcohols in Binary Mixtures of Water with Glycerol and Ethanol

Temperature Dependence of Peptide Conformational Equilibria from Simulations at a Single Temperature

Understanding the structure of proteins is key to unraveling their function in biological processes. Thus, significant attention has been paid to the calculation of conformational free energies. In this paper, we demonstrate a simple extension of fluctuation theory that permits the calculation of the temperature derivative of the conformational free energy, and hence the internal energy and entropy, from single-temperature simulations. The method further enables the decomposition into the contribution of different interactions present in the system to the internal energy surface. We illustrate the method for the canonical test system of alanine dipeptide in aqueous solution, for which we examine the free energy as a function of two dihedral angles. This system, like many, is most effectively treated using accelerated sampling methods and we show how the present approach is compatible with an important class of these, those that introduce a bias potential, by implementing it within metadynamics.

Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis

Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.

Glycine in Water Favors the Polyproline II State

Conformational preferences of amino acid residues in water are determined by the backbone and side-chain properties. Alanine is known for its high polyproline II (pPII) propensity. The question of relative contributions of the backbone and side chain to the conformational preferences of alanine and other amino acid residues in water is not fully resolved. Because glycine lacks a heavy-atom side chain, glycine-based peptides can be used to examine to which extent the backbone properties affect the conformational space. Here, we use published spectroscopic data for the central glycine residue of cationic triglycine in water to demonstrate that its conformational space is dominated by the pPII state. We assess three commonly used molecular dynamics (MD) force fields with respect to their ability to capture the conformational preferences of the central glycine residue in triglycine. We show that pPII is the mesostate that enables the functional backbone groups of the central residue to form the most hydrogen bonds with water. Our results indicate that the pPII propensity of the central glycine in GGG is comparable to that of alanine in GAG, implying that the water-backbone hydrogen bonding is responsible for the high pPII content of these residues.

Peptide-based short single β-strand mimics without hydrogen bonding or aggregation

Generation of short peptides with a single β-strand structure in solution is difficult. Herein, we design a new class of single β-strand peptidic mimics that are stable without self-aggregation in protic and non-protic solvents. Introduction of one present β-strand mimic can induce and propagate the β-strand structure for at least a penta-peptide sequence.

Demixing of water and ethanol causes conformational redistribution and gelation of the cationic GAG tripeptide

The cationic tripeptide GAG undergoes three conformational changes in binary mixtures of water and ethanol. At 17 mol% of ethanol conformational sampling is shifted from pPII towards β-strands. A more pronounced shift in the same direction occurs at 40 mol%. At ca. 55 mol% of ethanol and above a peptide concentration of ca. 0.2 M the ternary peptide-water-ethanol mixture forms a hydrogel which is comprised of unusually large crystalline like non-β sheet fibrils forming a sample spanning matrix.

Assessing backbone solvation effects in the conformational propensities of amino acid residues in unfolded peptides

Conformational ensembles of individual amino acid residues within model GxG peptides (x representing different amino acid residues) are dominated by a mixture of polyproline II (pPII) and β-strand like conformations. We recently discovered rather substantial differences between the enthalpic and entropic contributions to this equilibrium for different amino acid residues. Isoleucine and valine exceed all other amino acid residues in terms of their rather large enthalpic stabilization and entropic destabilization of polyproline II. In order to shed light on these underlying physical mechanisms, we performed high-level DFT calculations to explore the energetics of four representative GxG peptides where x = alanine (A), leucine (L), valine (V), and isoleucine (I) in explicit water (10 H2O molecules with a polarizable continuum water model) and in vacuo. We found that the large energetic contributions to the stabilization of pPII result, to a major extent, from peptide-water, water-water interactions, and changes of the solvent self-energy. Differences between the peptide-solvent interaction energies of hydration in pPII and β-strand peptides are particularly important for the pPII ⇌ β equilibria of the more aliphatic peptides GIG and GLG. Furthermore, we performed a vibrational analysis of the four peptides in both conformations and discovered a rather substantial mixing between water motions and peptide vibrations below 700 cm(-1). We found that the respective vibrational entropies are substantially different for the considered conformations, and their contributions to the Gibbs/Helmholtz energy stabilize β-strand conformations. Taken together, our results underscore the notion of the solvent being the predominant determinant of peptide (and protein) conformations in the unfolded state.

Conformational Changes of Trialanine in Water Induced by Vibrational Relaxation of the Amide I Mode

Most of the protein-based diseases are caused by anomalies in the functionality and stability of these molecules. Experimental and theoretical studies of the conformational dynamics of proteins are becoming in this respect essential to understand the origin of these anomalies. However, a description of the conformational dynamics of proteins based on mechano-energetic principles still remains elusive because of the intrinsic high flexibility of the peptide chains, the participation of weak noncovalent interactions, and the role of the ubiquitous water solvent. In this work, the conformational dynamics of trialanine dissolved in water (D2O) is investigated through Molecular Dynamics (MD) simulations combined with instantaneous normal modes (INMs) analysis both at equilibrium and after the vibrational excitation of the C-terminal amide I mode. The conformational equilibrium between α and pPII conformers is found to be altered by the intramolecular relaxation of the amide I mode as a consequence of the different relaxation pathways of each conformer which modify the amount of vibrational energy stored in the torsional motions of the tripeptide, so the α → pPII and pPII → α conversion rates are increased differently. The selectivity of the process comes from the shifts of the vibrational frequencies with the conformational changes that modify the resonance conditions driving the intramolecular energy flows.

Water-Centered Interpretation of Intrinsic pPII Propensities of Amino Acid Residues: In Vitro-Driven Molecular Dynamics Study

Amino acid residues of unfolded peptides in water sample only a few basins in the Ramachandran plot, including prominent polyproline II-like (pPII) conformations. Dynamics of guest residues, X, in GXG peptides in water were recently reported to be dominated by pPII and β-strand-like (β) conformations, resulting in an enthalpy-entropy compensation at ∼300 K. Using molecular dynamics (MD) in explicit solvent, we here examine pPII and β conformational ensembles of 15 guest residues in GXG peptides, quantify local orientation of water around their side chains through novel water orientation plots, and study their hydration and hydrogen bonding properties. We show that pPII and β ensembles are characterized by distinct water orientations: pPII ensembles are associated with an increased population of water oriented in parallel to the side chain surface whereas β ensembles exhibit more heterogeneous water orientations. The backbone hydration is significantly higher in pPII than in β ensembles. Importantly, pPII to β hydration differences and the solvent accessible surface area of Cβ hydrogens both correlate with experimental pPII propensities. We propose that pPII conformations are stabilized by a local, hydrogen-bonded clathrate-like water structure and that residue-specific intrinsic pPII propensities reflect distinct abilities of side chains to template this water structure.

Randomizing the unfolded state of peptides (and proteins) by nearest neighbor interactions between unlike residues

To explore the influence of nearest neighbors on conformational biases in unfolded peptides, we combined vibrational and 2D NMR spectroscopy to obtain the conformational distributions of selected "GxyG" host-guest peptides in aqueous solution: GDyG, GSyG, GxLG, GxVG, where x/y=A, K, L, V. Large changes of conformational propensities were observed due to nearest-neighbor interactions, at variance with the isolated pair hypothesis. We found that protonated aspartic acid and serine lose their above-the-average preference for turn-like structures in favor of polyproline II (pPII) populations in the presence of neighbors with bulky side chains. Such residues also decrease the above-the-average pPII preference of alanine. These observations suggest that the underlying mechanism involves a disruption of the hydration shell. Thermodynamic analysis of (3) J(H(N) ,H(α) ) (T) data for each x,y residue reveals that modest changes in the conformational ensemble masks larger changes of enthalpy and entropy governing the pPII↔β equilibrium indicating a significant residue dependent temperature dependence of the peptides' conformational ensembles. These results suggest that nearest-neighbor interactions between unlike residues act as conformational randomizers close to the enthalpy-entropy compensation temperature, eliminating intrinsic biases in favor of largely balanced pPII/β dominated ensembles at physiological temperatures.

Local order in the unfolded state: conformational biases and nearest neighbor interactions

The discovery of Intrinsically Disordered Proteins, which contain significant levels of disorder yet perform complex biologically functions, as well as unwanted aggregation, has motivated numerous experimental and theoretical studies aimed at describing residue-level conformational ensembles. Multiple lines of evidence gathered over the last 15 years strongly suggest that amino acids residues display unique and restricted conformational preferences in the unfolded state of peptides and proteins, contrary to one of the basic assumptions of the canonical random coil model. To fully understand residue level order/disorder, however, one has to gain a quantitative, experimentally based picture of conformational distributions and to determine the physical basis underlying residue-level conformational biases. Here, we review the experimental, computational and bioinformatic evidence for conformational preferences of amino acid residues in (mostly short) peptides that can be utilized as suitable model systems for unfolded states of peptides and proteins. In this context particular attention is paid to the alleged high polyproline II preference of alanine. We discuss how these conformational propensities may be modulated by peptide solvent interactions and so called nearest-neighbor interactions. The relevance of conformational propensities for the protein folding problem and the understanding of IDPs is briefly discussed.

The TrpA protein of Trichodesmium erythraeum IMS101 is a non-fibril-forming collagen and a component of the outer sheath

Collagen molecules are structural in nature and primarily found in eukaryotic, multicellular organisms. Recently, a collagen-like protein, TrpA, was identified and characterized in the marine cyanobacterium Trichodesmium erythraeum IMS 101, and it was shown to be involved in maintaining the structural integrity of the trichomes. The TrpA protein contains one glycine interruption in the otherwise perfectly uninterrupted collagenous domain. In this study, we used phylogenetic analysis to determine that the TrpA protein sequence is most closely associated with non-fibril-forming collagen proteins. Structural modelling and circular dichroism data suggest that the glycine insertion decreases the stability of TrpA compared to uninterrupted collagen sequences. Additionally, scanning electron microscopy revealed that TrpA is expressed entirely on the surface of the trichomes, with no specific pattern of localization. These data indicate that the TrpA protein is part of the outer sheath of this organism. As such, this protein may function to promote adhesion between individual T. erythraeum trichomes, and between this organism and heterotrophic bacteria found in the same environment.

pH-Independence of trialanine and the effects of termini blocking in short peptides: a combined vibrational, NMR, UVCD, and molecular dynamics study

Several lines of evidence now well establish that unfolded peptides in general, and alanine in specific, have an intrinsic preference for the polyproline II (pPII) conformation. Investigation of local order in the unfolded state is, however, complicated by experimental limitations and the inherent dynamics of the system, which has in some cases yielded inconsistent results from different types of experiments. One method of studying these systems is the use of short model peptides, and specifically short alanine peptides, known for predominantly sampling pPII structure in aqueous solution. Recently, He et al. ( J. Am. Chem. Soc. 2012 , 134 , 1571 - 1576 ) proposed that unblocked tripeptides may not be suitable models for studying conformational propensities in unfolded peptides due to the presence of end effect, that is, electrostatic interactions between investigated amino acid residues and terminal charges. To determine whether changing the protonation states of the N- and C-termini influence the conformational manifold of the central amino acid residue in tripeptides, we have examined the pH-dependence of unblocked trialanine and the conformational preferences of alanine in the alanine dipeptide. To this end, we measured and globally analyzed amide I' band profiles and NMR J-coupling constants. We described conformational distributions as the superposition of two-dimensional Gaussian distributions assignable to specific subspaces of the Ramachandran plot. Results show that the conformational ensemble of trialanine as a whole, and the pPII content (χpPII = 0.84) in particular, remains practically unaffected by changing the protonation state. We found that compared to trialanine, the alanine dipeptide has slightly lower pPII content (χpPII = 0.74) and an ensemble more reminiscent of the unblocked Gly-Ala-Gly model peptide. In addition, a two-state thermodynamic analysis of the conformational sensitive Δε(T) and (3)J(H(N)H(α))(T) data obtained from electronic circular dichroism and H NMR spectra indicate that the free energy landscape of trialanine is similar in all protonation states. MD simulations for the investigated peptides corroborate this notion and show further that the hydration shell around unblocked trialanine is unaffected by the protonation/deprotonation of the C-terminal group. In contrast, the alanine dipeptide shows a reduced water density around the central residue as well as a less ordered hydration shell, which decreases the pPII propensity and reduces the lifetime of sampled conformations.

Disorder and order in unfolded and disordered peptides and proteins: a view derived from tripeptide conformational analysis. II. Tripeptides with short side chains populating asx and β-type like turn conformations

In the preceding paper, we found that ensembles of tripeptides with long or bulky chains can include up to 20% of various turns. Here, we determine the structural and thermodynamic characteristics of GxG peptides with short polar and/or ionizable central residues (D, N, C), whose conformational distributions exhibit higher than average percentage (>20%) of turn conformations. To probe the side-chain conformations of these peptides, we determined the (3)J(H(α),H(β)) coupling constants and derived the population of three rotamers with χ1 -angles of -60°, 180° and 60°, which were correlated with residue propensities by DFT-calculations. For protonated GDG, the rotamer distribution provides additional evidence for asx-turns. A comparison of vibrational spectra and NMR coupling constants of protonated GDG, ionized GDG, and the protonated aspartic acid dipeptide revealed that side chain protonation increases the pPII content at the expense of turn populations. The charged terminal groups, however, have negligible influence on the conformational properties of the central residue. Like protonated GDG, cationic GCG samples asx-turns to a significant extent. The temperature dependence of the UVCD spectra and (3)J(H(N)H(α)) constants suggest that the turn populations of GDG and GNG are practically temperature-independent, indicating enthalpic and entropic stabilization. The temperature-independent J-coupling and UVCD spectra of GNG require a three-state model. Our results indicate that short side chains with hydrogen bonding capability in GxG segments of proteins may serve as hinge regions for establishing compact structures of unfolded proteins and peptides.