Stable conformations of tripeptides in aqueous solution studied by UV circular dichroism spectroscopy

Conformational Manifold Sampled by Two Short Linear Motif Segments Probed by Circular Dichroism, Vibrational, and Nuclear Magnetic Resonance Spectroscopy

Disordered protein segments called short linear motifs (SLiM) serve as recognition sites for a variety of biological processes and act as targeting signals, modification, and ligand binding sites. While SLiMs do not adopt one of the known regular secondary structures, the conformational distribution might still reflect the structural propensities of their amino acid residues and possible interactions between them. In the past, conformational analyses of short peptides provided compelling evidence for the notion that individual residues are less conformationally flexible than locally expected for a random coil. Here, we combined various spectroscopies (NMR, IR, vibrational, and UV circular dichroism) to determine the Ramachandran plots of two SLiM motifs, i.e., GRRDSG and GRRTSG. They are two representatives of RxxS motifs that are capable of being phosphorylated by protein kinase A, an enzyme that plays a fundamental role in a variety of biological processes. Our results reveal that the nearest and non-nearest interactions between residues cause redistributions between polyproline II and β-strand basins while concomitantly stabilizing extended relative to turn-forming and helical structures. They also cause shifts in basin positions. With increasing temperature, β-strand populations become more populated at the expense of polyproline II. While molecular dynamics simulations with Amber ff14SB and CHARMM 36m force fields indicate residue-residue interactions, they do not account for the observed structural changes.

The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins

Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.

Introducing Aliphatic Fluoropeptides: Perspectives on Folding Properties, Membrane Partition and Proteolytic Stability

A de novo designed class of peptide-based fluoropolymers composed of fluorinated aliphatic amino acids as main components is reported. Structural characterization provided insights into fluorine-induced alterations on β-strand to α-helix transition upon an increase in SDS content and revealed the unique formation of PPII structures for trifluorinated fluoropeptides. A combination of circular dichroism, fluorescence-based leaking assays and surface enhanced infrared absorption spectroscopy served to examine the insertion and folding processes into unilamellar vesicles. While partitioning into lipid bilayers, the degree of fluorination conducts a decrease in α-helical content. Furthermore, this study comprises a report on the proteolytic stability of peptides exclusively built up by fluorinated amino acids and proved all sequences to be enzymatically degradable despite the degree of fluorination. Herein presented fluoropeptides as well as the distinctive properties of these artificial and polyfluorinated foldamers with enzyme-degradable features will play a crucial role in the future development of fluorinated peptide-based biomaterials.

Exploring Nearest Neighbor Interactions and Their Influence on the Gibbs Energy Landscape of Unfolded Proteins and Peptides

The Flory isolated pair hypothesis (IPH) is one of the corner stones of the random coil model, which is generally invoked to describe the conformational dynamics of unfolded and intrinsically disordered proteins (IDPs). It stipulates, that individual residues sample the entire sterically allowed space of the Ramachandran plot without exhibiting any correlations with the conformational dynamics of its neighbors. However, multiple lines of computational, bioinformatic and experimental evidence suggest that nearest neighbors have a significant influence on the conformational sampling of amino acid residues. This implies that the conformational entropy of unfolded polypeptides and proteins is much less than one would expect based on the Ramachandran plots of individual residues. A further implication is that the Gibbs energies of residues in unfolded proteins or polypeptides are not additive. This review provides an overview of what is currently known and what has yet to be explored regarding nearest neighbor interactions in unfolded proteins.

Stereostructure Dependence Phenomenon on the Self-Assembly of Ala-Ala-Ala Lipotripeptides

A series of lipotripeptide stereoisomers based on alanine were synthesized, and their self-assembling behaviors were studied by means of circular dichroism spectra, ATR-IR, temperature-dependent ¹H NMR, and X-ray diffraction patterns. In the mixed solvent of hexafluoroisopropanol/H₂O (1/9, v/v), eight lipotripeptides were able to self-assembled into nanoflakes or nanoribbons driven by the hydrophobic association of alkyl chains, intermolecular hydrogen bonding among carboxyl groups at C-terminal and amide groups of alanine moieties in the peptide segment. It was found that the stacking chirality of carbonyl groups was determined by the chirality of alanine residue at C-terminal (i.e., "C-terminal determination" rule). Moreover, our research also highlighted the intermolecular hydrogen bonding on amide groups of each alanine residue, terminal carboxyl as well as the molecular packing structures can be subtly manipulated by changing the stereochemical sequence of peptide segment.

Conformation and Metal Cation Binding of Zwitterionic Alanine Tripeptide in Saline Solutions by Infrared Vibrational Spectroscopy and Molecular Dynamics Simulations

In this work, linear infrared (IR) spectroscopy and molecular dynamics (MD) simulations were used to examine the interaction of different metal cations (Na⁺, Ca²⁺, Mg²⁺, and Zn²⁺) with backbone (amide C═O) and C-terminal carboxylate (COO^-) groups in zwitterionic alanine tripeptide (Ala3) in aqueous solutions with varying saline concentrations. Circular dichroism spectra and MD results suggest that Ala3 is predominantly in polyproline-II (PP_II) conformation, whose amide-I and asymmetric carboxylate stretching IR vibration signatures are also supported by quantum-chemistry calculations. The zwitterionic form of Ala3 separates the two amide-I modes in frequency, which are weakly coupled modes, as revealed by two-dimensional IR measurement, and can be used to probe backbone-cation interactions at different scenarios (near charged or neutral chemical groups respectively). Cation concentration-dependent IR frequency red shifts in the amide-I mode are seen for both amide-I modes, whereas blue shifts are also seen in the amide-I mode far from the NH₃⁺ group. The observed spectral changes are discussed from the perspective of the salting-in and salting-out abilities of the cations. In addition, all the metal cations studied here (except Zn²⁺) can specifically coordinate to the COO^- group in bidentate and pseudo-bridging forms simultaneously. For Zn²⁺, only the pseudo-bridging form exists. Our results shed light on the macroscopic protein salting-in and salting-out phenomena from the perspective of key chemical bonds in peptides.

Reversible electrochemical triggering and optical interrogation of polylysine α-helix formation

Reversible electrochemical triggering of the random coil to α-helix conformational transition of polylysine (Lys₁₀, Lys₂₀, Lys₅₀) was accomplished at a Pt electrode at potentials < |1| V vs. Ag/AgCl. Direct electroreduction of the N-terminus vs ε-amino groups in Lys sidechains, as well as hydronium reduction and electrolysis, could be easily distinguished and deconvolved using differential pulse voltammetry. Electrochemistry was coupled with in situ UV absorbance and circular dichroism spectroscopies to dynamically follow the evolution of α-helix formation at different potentials. Isotope experiments in H₂O vs. D₂O unequivocally confirm that direct electroreduction of ε-NH₃⁺/ND₃⁺ groups in Lys sidechains, rather than electrochemically generated pH gradient-induced deprotonation, leads to subsequent α-helix formation. The site-selective electrochemistry and optical methodologies presented herein can be generalized and extended to interrogate other protonation-sensitive biomolecular systems, and potentially provide access to early intermediates and control over the dynamic structural evolution of peptides and proteins.

Novel chiral naphthalimide-cycloalkanediamine conjugates: Design, synthesis and antitumor activity

A novel series of enantiopure naphthalimide-cycloalkanediamine conjugates were designed, synthetized and evaluated for in vitro cytotoxicity against human colon adenocarcinoma (LoVo), human lung adenocarcinoma (A549), human cervical carcinoma (Hela) and human promyelocytic leukemia cell lines (HL-60). The cytotoxicity of the compounds was highly dependent on size and relative stereochemistry of the cycloalkyl ring as well as length of the spacer. By contrast, any kind of enantioselection was observed for each pair of enantiomers. Flow cytometric analysis indicated that compounds 22 and 23 could effectively induce G₂/M arrest in the four previous cell lines despite a mild apoptotic effect.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

Multiconformational analysis of tripeptides upon consideration of implicit and explicit hydration effects

During the last two decades, numerous observed data obtained by various physical techniques, also supported by molecular modeling approaches, have highlighted the structuring features of tripeptides, as well as their aggregation properties. Herein, we focus on the structural dynamics of four trimers, i.e., Gly-Gly-Gly, Gly-Ala-Gly, Ala-Ala-Ala and Ala-Phe-Ala, in an aqueous environment. Density functional theory calculations (DFT) were carried out to assess the stability of four types of secondary structures, i.e., β-strand, polyproline-II (pP-II), α-helix and γ-turn, of which the formation had been described in these tripeptides. Both implicit and explicit hydration effects were analyzed on the conformational and energetic features of trimers. It has been shown that the use of M062X functional (versus B3LYP) improve the stability of intramolecular H-bonds, especially in inverse γ-turn structures, as well as the energetic and conformational equilibrium in all tripeptides. Explicit hydration reflected by the presence of five water molecules around the backbone polar sites (NH₃⁺, N-H, CO and NH₂) considerably changes the conformational landscapes of the trimers. Characteristic intramolecular and intermolecular interactions evidenced by the calculations, were emphasized.

Water-Mediated Electronic Structure of Oligopeptides Probed by Their UV Circular Dichroism, Absorption Spectra, and Time-Dependent DFT Calculations

We investigate the UV absorption spectra of a series of cationic GxG peptides (where x denotes a guest residue) in aqueous solution and find that only a subset of these spectra show a strong dependence with temperature. To explore whether or not this observation reflects conformational dependencies, we carry out time-dependent density functional calculations for the polyproline II (pPII) and β-strand conformations in implicit and explicit water. We find that the calculated CD spectra for pPII can qualitatively account for the experimental spectra irrespective of the water model. The β-strand UV-CD spectra, however, require the explicit consideration of water. Contrary to conventional wisdom, we find that both the NV₁ and NV₂ band are the envelopes of contributions from multiple transitions that involve more than just the HOMOs and LUMOs of the peptide groups. A natural transition orbital analysis reveals that some of the transitions have a charge-transfer character. The overall manifold of transitions depends on the peptide's backbone conformation, peptide hydration, and side chain of the guest residue. Our results reveal that peptide groups, side chains, and hydration shells must be considered as an entity for a physically valid characterization of UV absorbance and circular dichroism.

Anticooperative Nearest-Neighbor Interactions between Residues in Unfolded Peptides and Proteins

Growing evidence suggests that the conformational distributions of amino acid residues in unfolded peptides and proteins depend on the nature of the nearest neighbors. To explore whether the underlying interactions would lead to a breakdown of the isolated pair hypothesis of the classical random coil model, we further analyzed the conformational propensities that were recently obtained for the two guest residues (x,y) of GxyG tetrapeptides. We constructed a statistical thermodynamics model that allows for cooperative as well as for anticooperative interactions between adjacent residues adopting either a polyproline II or a β-strand conformation. Our analysis reveals that the nearest-neighbor interactions between most of the central residues in the investigated GxyG peptides are anticooperative. Interaction Gibbs energies are rather large at high temperatures (350 K), at which point many proteins undergo thermal unfolding. At room temperature, these interaction energies are less pronounced. We used the obtained interaction parameter in a Zimm-Bragg/Ising-type approach to calculate the temperature dependence of the ultraviolet circular dichroism (CD) of the MAX3 peptide, which is predominantly built by KV repeats. The agreement between simulation and experimental data was found to be satisfactory. Finally, we analyzed the temperature dependence of the CD and ³J(H^NH^α) parameters of the amyloid β_1-9 fragment. The results of this analysis and a more qualitative consideration of the temperature dependence of denatured proteins probed by CD spectroscopy further corroborate the dominance of anticooperative nearest-neighbor interactions. Generally, our results show that unfolded peptides-and most likely also proteins-exhibit some similarity with antiferromagnetic systems.

A novel thermo-responsive multiblock architecture composed of a sequential peptide and an amino acid-derived vinyl polymer: toward protein-mimicking single-chain folding

A novel multiblock hybrid copolymer composed of a β-sheet peptide and a glycine-derived vinyl polymer was developed to achieve single-chain folding into well-defined nanoparticles.

Understanding the Interaction Between Oligopeptide and Water in Aqueous Solution Using Temperature-Dependent Near-Infrared Spectroscopy

Investigating the interaction between oligopeptide and water is essential for understanding the structure, dynamics and function of proteins. Temperature-dependent near-infrared (NIR) spectroscopy and independent component analysis (ICA) were employed to study the interaction between oligopeptide and water in aqueous solution. The NIR spectra of two homo-oligopeptides, penta-aspartic acid (D₅) and penta-lysine (K₅), in aqueous solution of different concentration were measured at different temperature (30-90 ℃). Independent component analysis was performed to extract the spectral information that changes with temperature. The independent components (ICs) representing the spectral information of NH and CH₂ groups were obtained. Compared with D₅, the two groups in K₅ change significantly at higher temperature. The result may suggest that K₅ has stronger interaction with water than D₅. Moreover, three ICs that contain the spectral information of the water species with no (S₀), one (S₁), and two (S₂) hydrogen-bonds were obtained. It was shown that the spectral intensity of S₀ and S₁ increases while that of S₂ decreases with the temperature, and the changes of oligopeptide solutions are weaker than those of pure water. The results indicate that water structure is sensitive to temperature and the oligopeptide in aqueous solution improves the thermal stability of the water species. When oligopeptide is added, the spectral intensity of S₀ and S₂ decreases and that of S₁ increases for D₅ solution, but the intensity of all the three species decreases for K₅ solution. Furthermore, the concentration effect of K₅ was found to be stronger than D₅. The result may reveal that D₅ combines with water molecule through forming one hydrogen bond but K₅ interacts with water through a different way.

Molecular dynamics-derived rotamer libraries for d-amino acids within homochiral and heterochiral polypeptides

Computational resources have contributed to the design and engineering of novel proteins by integrating genomic, structural and dynamic aspects of proteins. Non-canonical amino acids, such as d-amino acids, expand the available sequence space for designing and engineering proteins; however, the rotamer libraries for d-amino acids are usually constructed as the mirror images of l-amino acid rotamer libraries, an assumption that has not been tested. To this end, we have performed molecular dynamics (MD) simulations of model host-guest peptide systems containing d-amino acids. Our simulations systematically address the applicability of the mirror image convention as well as the effects of neighboring residue chirality. Rotamer libraries derived from these systems provide realistic rotamer distributions suitable for use in both rational and computational design workflows. Our simulations also address the impact of chirality on the intrinsic conformational preferences of amino acids, providing fundamental insights into the relationship between chirality and biomolecular dynamics. While d-amino acids are rare in naturally occurring proteins, they are used in designed proteins to stabilize a desired conformation, increase bioavailability or confer favorable biochemical and physical attributes. Here, we present d-amino acid rotamer libraries derived from MD simulations of alanine-based host-guest pentapeptides and show how certain residues can deviate from mirror image symmetry. Our simulations directly model d-amino acids as guest residues within the chiral l-Ala and d-Ala pentapeptide series to explicitly incorporate any contributions resulting from the chiralities of neighboring residues.

Quantum Mechanics Study on Hydrophilic and Hydrophobic Interactions in the Trivaline-Water System

With the aim to elucidate hydrophobic effects in the unfolded state of peptides, DFT-M062X computations on the Val₃H⁺· nH₂O ( n up to 22) clusters have been accomplished. As far as the main chain is concerned, four conformers with β-strand and/or polyproline type II conformations, PPII (indicated as β-β, β-PPII, PPII-β, and PPII-PPII), have been found by changing the ϕ and ψ angles. For bare peptide, the side chain (isopropyl) of each residue can independently take on three different orientations with negligible effects on energetics. The great isopropyl spatial separations in β-β and β-PPII conformers allow for the construction of synergic and extensive water-water and water-peptide H-bonding in the minimal hydration Val₃H⁺·22H₂O models without significant steric encumbrance. Conversely, due to the proximity of the isopropyl of the central residue with the other two, some restrictions in the water shell construction around the peptide become evident for the PPII-PPII conformer and the number of energetically accessible structures decreases. This is indicative of correlated motion involving isopropyls and backbone mediated by water molecules, the origin of the nearest neighbor effects. Comparing the thermodynamic data of Ala₃H⁺·22H₂O and Val₃H⁺·22H₂O, what emerges is that both hydration enthalpy and entropy drive the β-strand stability of the latter.

EXAT: EXcitonic analysis tool

We introduce EXcitonic Analysis Tool (EXAT), a program able to compute optical spectra of large excitonic systems directly from the output of quantum mechanical calculations performed with the popular Gaussian 16 package. The software is able to combine in an excitonic scheme the single-chromophore properties and exciton couplings to simulate energies, coefficients, and excitonic spectra (UV-vis, CD, and LD). The effect of the environment can also be included using a Polarizable Continuum Model. EXAT also presents a simple graphical user interface, which shows on-screen both site and exciton properties. To show the potential of the method, we report two applications on a a chiral perturbed BODIPY system and DNA G-quadruplexes, respectively. The program is available online at http://molecolab.dcci.unipi.it/tools/. © 2017 Wiley Periodicals, Inc.

N-terminal diproline and charge group effects on the stabilization of helical conformation in alanine-based short peptides: CD studies with water and methanol as solvent

Protein folding problem remains a formidable challenge as main chain, side chain and solvent interactions remain entangled and have been difficult to resolve. Alanine-based short peptides are promising models to dissect protein folding initiation and propagation structurally as well as energetically. The effect of N-terminal diproline and charged side chains is assessed on the stabilization of helical conformation in alanine-based short peptides using circular dichroism (CD) with water and methanol as solvent. A1 (Ac-Pro-Pro-Ala-Lys-Ala-Lys-Ala-Lys-Ala-NH₂ ) is designed to assess the effect of N-terminal homochiral diproline and lysine side chains to induce helical conformation. A2 (Ac-Pro-Pro-Glu-Glu-Ala-Ala-Lys-Lys-Ala-NH₂ ) and A3 (Ac-dPro-Pro-Glu-Glu-Ala-Ala-Lys-Lys-Ala-NH₂ ) with N-terminal homochiral and heterochiral diproline, respectively, are designed to assess the effect of Glu...Lys (i, i + 4) salt bridge interactions on the stabilization of helical conformation. The CD spectra of A1, A2 and A3 in water manifest different amplitudes of the observed polyproline II (PPII) signals, which indicate different conformational distributions of the polypeptide structure. The strong effect of solvent substitution from water to methanol is observed for the peptides, and CD spectra in methanol evidence A2 and A3 as helical folds. Temperature-dependent CD spectra of A1 and A2 in water depict an isodichroic point reflecting coexistence of two conformations, PPII and β-strand conformation, which is consistent with the previous studies. The results illuminate the effect of N-terminal diproline and charged side chains in dictating the preferences for extended-β, semi-extended PPII and helical conformation in alanine-based short peptides. The results of the present study will enhance our understanding on stabilization of helical conformation in short peptides and hence aid in the design of novel peptides with helical structures. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Conformational sampling of metastable states: Tq-REM as a novel replica exchange method

Although the replica exchange methods (REMs) were developed as efficient conformational sampling methods for bio-molecular simulations, their application to very large bio-systems is somewhat limited. We propose a new replica exchange scheme (Tq-REM) created by combining the conventional temperature-REM (T-REM) and one of the Hamiltonian-REMs, q-REM, using the effective potential with reduced barriers. In the proposed Tq-REM scheme, high temperature replicas in T-REM are substituted with q-replicas. This combined scheme is expected to exploit advantages of the T-REM and q-REM resulting in improved sampling efficiency while minimizing the drawbacks of both approaches. We investigated the performance of Tq-REM compared with T-REM by performing all-atom MD simulations on Met-enkephalin, (AAQAA)₃, and Trpzip2. It was found that convergence of the free energy surfaces was improved by Tq-REM over the conventional T-REM. In particular, the trajectories of Tq-REM were able to sample the relevant conformations for all of the metastable folding intermediates, while some of the local minimum structures are poorly represented by T-REM. The results of the present study suggest that Tq-REM can provide useful tools to investigate systems where metastable states play important roles.

Phosphorylation Increases Persistence Length and End-to-End Distance of a Segment of Tau Protein

Intrinsically disordered regions of proteins, which lack unique tertiary structure under physiological conditions, are enriched in phosphorylation sites and in significant local bias toward the polyproline II conformation. The overrepresented coincidence of this posttranslational regulatory signal and local conformational bias within unstructured regions raises a question: can phosphorylation serve to manipulate the conformational preferences of a disordered protein? In this study, we use time-resolved fluorescence resonance energy transfer and a, to our knowledge, novel data analysis method to directly measure the end-to-end distance distribution of a phosphorylatable peptide derived from the human microtubule associated protein tau. Our results show that phosphorylation at threonine or serine extends the end-to-end distance and increases the effective persistence length of the tested model peptides. Unexpectedly, the extension is independent of salt concentration, suggestive of a nonelectrostatic origin. The phosphorylation extension and stiffening effect provides a peptide-scale physical interpretation for the posttranslational regulation of the highly abundant protein-protein interactions found in disordered proteins, as well as a potential insight into the regulatory mechanism of the tau protein's microtubule binding activity.

The interplay of aggregation, fibrillization and gelation of an unexpected low molecular weight gelator: glycylalanylglycine in ethanol/water

Hydrogels formed by polypeptides could be much-favored tools for drug delivery because their main ingredients are generally biodegradable. However, the gelation of peptides in aqueous solution generally requires a minimal length of the peptide as well as distinct sequences of hydrophilic and hydrophobic residues. The aggregation of short peptides like tripeptides, which are relatively cheap and offer a high degree of biodegradability, are generally thought to require a high hydrophobicity of their residues. We found that contrary to this expectation cationic glycylalanylglycine in 55 mol% ethanol/45 mol% water forms a gel below a melting temperature of ca. 36 °C. A pure hydrogel state can be obtained after allowing the ethanol component to evaporate. The gel phase consists of crystalline fibrils of several 100 μm, which form a sample-spanning network. Rheological data reveal a soft elastic solid gel. We investigated the kinetics of the various processes that lead to the final gel state of the ternary mixture by a unique combination of UV circular dichroism, infrared, vibrational circular dichroism (VCD) and rheological measurements. A mathematical analysis of our data show that gelation is preceded by the formation of peptide β-sheet like tapes or ribbons, which give rise to a significant enhancement of the amide I' VCD signal, and the subsequent formation of rather thick and long fibrils. The VCD signals indicate that the tapes exhibit a right-handed helicity at temperatures above 16 °C and a left-handed helicity below. The tapes'/ribbons' helicity change occurs at a temperature where the UVCD data reflect a relatively long nucleation process. The kinetics of gel formation probed by the storage and loss moduli are composed of a fast process that follows tape/ribbon/fibril formation and is clearly identifiable in a movie that shows the gelation process and a slow process that causes an additional gel stabilization. The rheological data indicate that left-handed fibrils observed at low temperatures form a more solid-like structure than their right-handed counterparts formed at higher temperatures. Taken together our data reveal GAG as an unexpected gelator, the formation of which is underlied by a set of distinguishable kinetic processes.

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.

Structure of Penta-Alanine Investigated by Two-Dimensional Infrared Spectroscopy and Molecular Dynamics Simulation

We have studied the structure of (Ala)5, a model unfolded peptide, using a combination of 2D IR spectroscopy and molecular dynamics (MD) simulation. Two different isotopomers, each bis-labeled with (13)C═O and (13)C═(18)O, were strategically designed to shift individual site frequencies and uncouple neighboring amide-I' modes. 2D IR spectra taken under the double-crossed ⟨π/4, -π/4, Y, Z⟩ polarization show that the labeled four-oscillator systems can be approximated by three two-oscillator systems. By utilizing the different polarization dependence of diagonal and cross peaks, we extracted the coupling constants and angles between three pairs of amide-I' transition dipoles through spectral fitting. These parameters were related to the peptide backbone dihedral angles through DFT calculated maps. The derived dihedral angles are all located in the polyproline-II (ppII) region of the Ramachandran plot. These results were compared to the conformations sampled by Hamiltonian replica-exchange MD simulations with three different CHARMM force fields. The C36 force field predicted that ppII is the dominant conformation, consistent with the experimental findings, whereas C22/CMAP predicted similar population for α+, β, and ppII, and the polarizable Drude-2013 predicted dominating β structure. Spectral simulation based on MD representative conformations and structure ensembles demonstrated the need to include multiple 2D spectral features, especially the cross-peak intensity ratio and shape, in structure determination. Using 2D reference spectra defined by the C36 structure ensemble, the best spectral simulation is achieved with nearly 100% ppII population, although the agreement with the experimental cross-peak intensity ratio is still insufficient. The dependence of population determination on the choice of reference structures/spectra and the current limitations on theoretical modeling relating peptide structures to spectral parameters are discussed. Compared with the previous results on alanine based oligopeptides, the dihedral angles of our fitted structure, and the most populated ppII structure from the C36 simulation are in good agreement with those suggesting a major ppII population. Our results provide further support for the importance of ppII conformation in the ensemble of unfolded peptides.

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.

Computational scrutiny of the effect of N-terminal proline and residue stereochemistry in the nucleation of α-helix fold

N-Terminal l- to d-residue mutation nucleate helical fold in Ac–^DAla–^LAla₃–NHMe (Ib, m2), Ac–^DPro–^LAla₃–NHMe (IIb, m1), and Ac–^DPro–^LPro–^LAla₂–NHMe (IIIb, m2) peptides.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Engineering of a peptide probe for β-amyloid aggregates

A non-self-aggregating peptide ligand for β-amyloid aggregates created by simple point mutation of an β-amyloid-derived segment.

Intrinsic α-helical and β-sheet conformational preferences: a computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.

Ionic strength of electrospray droplets affects charging of DNA oligonucleotides

The fundamental aspects of charging in electrospray ionization (ESI) are hotly debated. In the present study, ESI charging of DNA oligonucleotides was explored in both positive (ESI+) and negative (ESI-) polarity using mass spectrometry detection. Single-stranded 12-mer CCCCAATTCCCC in buffer solution (aqueous NH4Ac, 100 mM) produced similar charge state distribution (CSD) in either ESI+ or ESI-. Similarity of CSD in ESI+ and ESI- was also observed for the double-stranded 12-mer CGCGAATTCGCG. By adding typical low-vapor reagents (e.g. m-nitro benzyl alcohol, m-NBA; sulfolane) into the same buffer solution (<0.5% w/v), both CCCCAATTCCCC and CGCGAATTCGCG revealed strong supercharging (SC) effect in ESI-, while very little or no SC effect was observed in ESI+. With either sulfolane or m-NBA, the CGCGAATTCGCG duplex dissociated into single strands in ESI-. No SC was observed in both ESI+ and ESI- for thermally denatured CGCGAATTCGCG duplex in NH4 Ac buffer without the reagents. These findings are difficult to reconcile with the earlier model, which attributes SC in aqueous buffer solution to the conformational changes of analytes. Our observations suggest that the ionic strength of ESI droplets strongly affects the CSD of biopolymers such as DNA oligonucleotides and that SC effect is related to the depletion of ionic strength during the ESI process.

Predictive atomic resolution descriptions of intrinsically disordered hTau40 and α-synuclein in solution from NMR and small angle scattering

The development of molecular descriptions of intrinsically disordered proteins (IDPs) is essential for elucidating conformational transitions that characterize common neurodegenerative disorders. We use nuclear magnetic resonance, small angle scattering, and molecular ensemble approaches to characterize the IDPs Tau and α-synuclein. Ensemble descriptions of IDPs are highly underdetermined due to the inherently large number of degrees of conformational freedom compared with available experimental measurements. Using extensive cross-validation we show that five different types of independent experimental parameters are predicted more accurately by selected ensembles than by statistical coil descriptions. The improvement increases in regions whose local sampling deviates from statistical coil, validating the derived conformational description. Using these approaches we identify enhanced polyproline II sampling in aggregation-nucleation sites, supporting suggestions that this region of conformational space is important for aggregation.

The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development

The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.

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.

Solvent interactions stabilize the polyproline II conformation of glycosylated oligoprolines

In nature, proline residues carry several post-translational modifications (PTMs), including 4R hydroxylation and glycosylation. A recent study synthesized contiguously hydroxylated and glycosylated nonaproline peptides and revealed that both PTMs lead to a significant increase in the thermal stability of PPII relative to the unmodified oligoproline. The increased stability of the hydroxylated peptide can be explained by increased stability of the trans isomer due to stereoelectronic effects. However, the effects of glycosylation cannot be completely explained by stereoelectronics since previous experimental results indicate that 4R-glycosylation does not produce observable changes in the trans preference compared to 4R-hydroxylation. We therefore used sophisticated molecular modeling techniques to determine the reason for the further increase in thermal stability upon glycosylation. Free energy estimates obtained from adaptively biased molecular dynamics calculations in implicit (explicit) solvent are -9 kcal mol(-1) (-20 kcal mol(-1)) for the hydroxylated compound and -9 kcal mol(-1) (-46 kcal mol(-1)) for the glycosylated compound, indicating that direct solvent-peptide interactions are vital for explaining the glycosylation effects on PPII stability. Our data reveals for the first time that interactions between the hydroxyl groups in the glycosylated compound and water act in a complementary fashion with stereoelectronic effects to stabilize the PPII conformation in these substituted oligoproline peptides.

Effects of intramolecular distance between amyloidogenic domains on amyloid aggregation

Peptide/protein aggregation is implicated in many amyloid diseases. Some amyloidogenic peptides/proteins, such as those implicated in Alzheimer’s and Parkinson’s diseases, contain multiple amyloidogenic domains connected by “linker” sequences displaying high propensities to form turn structures. Recent studies have demonstrated the importance of physicochemical properties of each amino acid contained in the polypeptide sequences in amyloid aggregation. However, effects on aggregation related to the intramolecular distance between amyloidogenic domains, which may be determined by a linker length, have yet to be examined. In the study presented here, we created peptides containing two copies of KFFE, a simple four-residue amyloidogenic domain, connected by GS-rich linker sequences with different lengths yet similar physicochemical properties. Our experimental results indicate that aggregation occurred most rapidly when KFFE domains were connected by a linker of an intermediate length. Our experimental findings were consistent with estimated entropic contribution of a linker length toward formation of (partially) structured intermediates on the aggregation pathway. Moreover, inclusion of a relatively short linker was found to inhibit formation of aggregates with mature fibril morphology. When the results are assimilated, our study demonstrates that intramolecular distance between amyloidogenic domains is an important yet overlooked factor affecting amyloid aggregation.

Antifreeze glycoprotein agents: structural requirements for activity

Antifreeze glycoproteins (AFGPs) are considered to be the most efficient means to reduce ice damage to cell tissues since they are able to inhibit growth and crystallization of ice. The key element of antifreeze proteins is to act in a non-colligative manner which allows them to function at concentrations 300-500 times lowers than other dissolved solutes. During the past decade, AFGPs have demonstrated tremendous potential for many pharmaceutical and food applications. Presently, the only route to obtain AFGPs involves the time consuming and expensive process of isolation and purification from deep-sea polar fishes. Unfortunately, it is not amenable to mass production and commercial applications. The lack of understanding of the mechanism through which the AFGPs inhibit ice growth has also hampered the realization of industrial and biotechnological applications. Here we report the structural motifs that are essential for antifreeze activity of AFGPs, and propose a unified mechanism based on both recent studies of short alanine peptides and structure activity relationship of synthesized AFGPs.

Effects of solvents on the intrinsic propensity of peptide backbone conformations

We investigated the effects of solvents on the intrinsic propensity of peptide backbone conformations based on molecular dynamics simulations. The results show that compared with pure water, aqueous urea decreases the helix propensity. In comparison, methanol decreases the polyproline II (PPII) propensity. Such a solvent dependence of the intrinsic propensity of the backbone conformation is correlated with the solvent dependence of the hydration of the backbone groups and the formation probability of the local intrapeptide hydrogen bonds. Aqueous urea which has low ability to stabilize the local intrapeptide hydrogen bonds disfavors the helical conformation. Whereas, methanol which has low ability to hydrate the backbone groups disfavors the polyproline II conformation. In addition, the solvent effects can be further modulated by the side chains of the peptides. The solvent effects of the intrinsic propensity of peptide backbone conformations observed in this work suggest that changing the intrinsic propensity of the protein backbone conformations can partly contribute to the solvent-induced protein structure and dynamics variations. These results will be useful in understanding the solvent dependence of the conformational distributions of the unfolded proteins or peptides (or intrinsically disordered proteins) in which the global tertiary interactions are less important than that in the well-folded proteins.

Phosphorylation alters backbone conformational preferences of serine and threonine peptides

Despite the notion that a control of protein function by phosphorylation works mainly by inducing its conformational changes, the phosphorylation effects on even small peptide conformation have not been fully understood yet. To study its possible effects on serine and threonine peptide conformations, we recently carried out pH- and temperature-dependent circular dichroism (CD) as well as (1)H NMR studies of the phosphorylated serine and threonine peptides and compared them with their unphosphorylated analogs. In the present article, by performing the self-consistent singular value decomposition analysis of the temperature-dependent CD spectra and by analyzing the (3)J(H(N),H(α)) coupling constants extracted from the NMR spectra, the populations of the polyproline II (PPII) and β-strand conformers of the phosphorylated Ser and Thr peptides are determined. As temperature is increased, the β-strand populations of both phosphorylated serine and threonine peptides increase. However, the dependences of PPII/β-strand population ratio on pH are different for these two cases. The phosphorylation of the serine peptide enhances the PPII propensity, whereas that of the threonine peptide has the opposite effect. This suggests that the serine and threonine phosphorylations can alter the backbone conformational propensity via direct but selective intramolecular hydrogen-bonding interactions with the peptide N--H groups. This clearly indicates that the phosphoryl group actively participates in modulating the peptide backbone conformations.