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

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.

Statistical proofs of the interdependence between nearest neighbor effects on polypeptide backbone conformations

Backbone dihedral angles ϕ and ψ are the main structural descriptors of proteins and peptides. The distribution of these angles has been investigated over decades as they are essential for the validation and refinement of experimental measurements, as well as for structure prediction and design methods. The dependence of these distributions, not only on the nature of each amino acid but also on that of the closest neighbors, has been the subject of numerous studies. Although neighbor-dependent distributions are nowadays generally accepted as a good model, there is still some controversy about the combined effects of left and right neighbors. We have investigated this question using rigorous methods based on recently-developed statistical techniques. Our results unambiguously demonstrate that the influence of left and right neighbors cannot be considered independently. Consequently, three-residue fragments should be considered as the minimal building blocks to investigate polypeptide sequence-structure relationships.

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.

Randomizing of Oligopeptide Conformations by Nearest Neighbor Interactions between Amino Acid Residues

Flory’s random coil model assumes that conformational fluctuations of amino acid residues in unfolded poly(oligo)peptides and proteins are uncorrelated (isolated pair hypothesis, IPH). This implies that conformational energies, entropies and solvation free energies are all additive. Nearly 25 years ago, analyses of coil libraries cast some doubt on this notion, in that they revealed that aromatic, but also β-branched side chains, could change the ³J(H^NH^Cα) coupling of their neighbors. Since then, multiple bioinformatical, computational and experimental studies have revealed that conformational propensities of amino acids in unfolded peptides and proteins depend on their nearest neighbors. We used recently reported and newly obtained Ramachandran plots of tetra- and pentapeptides with non-terminal homo- and heterosequences of amino acid residues to quantitatively determine nearest neighbor coupling between them with a Ising type model. Results reveal that, depending on the choice of amino acid residue pairs, nearest neighbor interactions either stabilize or destabilize pairs of polyproline II and β-strand conformations. This leads to a redistribution of population between these conformations and a reduction in conformational entropy. Interactions between residues in polyproline II and turn(helix)-forming conformations seem to be cooperative in most cases, but the respective interaction parameters are subject to large statistical errors.

Repeating Aspartic Acid Residues Prefer Turn-like Conformations in the Unfolded State: Implications for Early Protein Folding

Protein folding can be described as a motion of the polypeptide chain in a potential energy funnel, where the conformational manifold is narrowed as the chain traverses from a completely unfolded state until it reaches the folded (native) state. The initial folding stages set the tone for this process by substantially narrowing the manifold of accessible conformations. In an ideally unfolded state with no long-range stabilizing forces, local conformations (i.e., residual structures) are likely to drive the folding process. While most amino acid residues tend to predominantly adopt extended structures in unfolded proteins and peptides, aspartic acid exhibits a relatively high intrinsic preference for turn-forming conformations. Regions in an unfolded polypeptide or protein that are rich in aspartic acid residues may therefore be crucial sites for protein folding steps. By combining NMR and vibrational spectroscopies, we observed that the conformational sampling of multiple sequentially neighbored aspartic acid residues in the model peptides GDDG and GDDDG even show an on average higher propensity for turn-forming structures than the intrinsic reference system D in GDG, which suggests that nearest neighbor interactions between adjacent aspartic acid residues stabilize local turn-forming structures. In the presence of the unlike neighbor phenylalanine, nearest neighbor interactions are of a totally different nature in that it they decrease the turn-forming propensities and mutually increase the sampling of polyproline II (pPII) conformations. We hypothesize the structural role of aspartic residues in intrinsically disordered proteins in general, and particularly in small linear motifs, that are very much determined by their respective neighbors.

Alternative Causal Link between Peptide Fibrillization and β-Strand Conformation

In the prevailing phenomenon of peptide fibrillization, β-strand conformation has long been believed to be an important structural basis for peptide assembly. According to a widely accepted theory, in most peptide fibrillization processes, peptide monomers need to intrinsically take or transform to β-strand conformation before they can undergo ordered packing to form nanofibers. In this study, we reported our findings on an alternative peptide fibrillization pathway starting from a disordered secondary structure, which could then transform to β-strand after fibrillization. By using circular dichroism, thioflavin-T binding test, and transmission electron microscopy, we studied the secondary structure and assembly behavior of Ac-RADARADARADARADA-NH2 (RADA16-I) in a low concentration range. The effects of peptide concentration, solvent polarity, pH, and temperature were investigated in detail. Our results showed that at very low concentrations, even though the peptide was in a disordered secondary structure, it could still form nanofibers through intermolecular assembly, and under higher peptide concentrations, the transformation from the disordered structure to β-strand could happen with the growth of nanofibers. Our results indicated that even without ordered β-strand conformation, driving forces such as hydrophobic interaction and electrostatic interaction could still play a determinative role in the self-assembly of peptides. At least in some cases, the formation of β-strand might be the consequence rather than the cause of peptide fibrillization.

Short peptides as predictors for the structure of polyarginine sequences in disordered proteins

Intrinsically disordered proteins and intrinsically disordered regions are frequently enriched in charged amino acids. Intrinsically disordered regions are regularly involved in important biological processes in which one or more charged residues is the driving force behind a protein-biomolecule interaction. Several lines of experimental and computational evidence suggest that polypeptides and proteins that carry high net charges have a high preference for extended conformations with average end-to-end distances exceeding expectations for self-avoiding random coils. Here, we show that charged arginine residues even in short glycine-capped model peptides (GRRG and GRRRG) significantly affect the conformational propensities of each other when compared with the intrinsic propensities of a mostly unperturbed arginine in the tripeptide GRG. A conformational analysis based on experimentally determined J-coupling constants from heteronuclear NMR spectroscopy and amide I' band profiles from vibrational spectroscopy reveals that nearest-neighbor interactions stabilize extended β-strand conformations at the expense of polyproline II and turn conformations. The results from molecular dynamics simulations with a CHARMM36m force field and TIP3P water reproduce our results only to a limited extent. The use of the Ramachandran distribution of the central residue of GRRRG in a calculation of end-to-end distances of polyarginines of different length yielded the expected power law behavior. The scaling coefficient of 0.66 suggests that such peptides would be more extended than predicted by a self-avoiding random walk. Our findings thus support in principle theoretical predictions.

On the Role of Entropy in the Stabilization of α-Helices

Protein folding evolves by exploring the conformational space with a subtle balance between enthalpy and entropy changes which eventually leads to a decrease of free energy upon reaching the folded structure. A complete understanding of this process requires, therefore, a deep insight into both contributions to free energy. In this work, we clarify the role of entropy in favoring the stabilization of folded structures in polyalanine peptides with up to 12 residues. We use a novel method referred to as K2V that allows us to obtain the potential-energy landscapes in terms of residue conformations extracted from molecular dynamics simulations at conformational equilibrium and yields folding thermodynamic magnitudes, which are in agreement with the experimental data available. Our results demonstrate that the folded structures of the larger polyalanine chains are stabilized with respect to the folded structures of the shorter chains by both an energetic contribution coming from the formation of the intramolecular hydrogen bonds and an entropic contribution coming from an increase of the entropy of the solvent with approximate weights of 60 and 40%, respectively, thus unveiling a key piece in the puzzle of protein folding. In addition, the ability of the K2V method to provide the enthalpic and entropic contributions for individual residues along the peptide chain makes it clear that the energetic and entropic stabilizations are basically governed by the nearest neighbor residue conformations, with the folding propensity being rationalized in terms of triads of residues.

Intrinsic Conformational Dynamics of Alanine in Water/Ethanol Mixtures: An Experiment-Driven Molecular Dynamics Study

In vitro, cationic glycylalanylglycine (GAG) forms a hydrogel in binary mixtures of water and ethanol. In water, alanine residue is known for its high polyproline II (pPII) content. Spectroscopic data, including three J-coupling constants and amide I' profiles, indicate that addition of 42% ethanol to water significantly reduces the pPII content of alanine residue in GAG. Here, experiment-based Gaussian Ramachandran distributions of alanine in GAG at different ethanol fractions are examined and three MD force fields are evaluated with respect to their ability to capture these ethanol-induced conformational changes. MD simulations on monomeric GAG in eight different water/ethanol mixtures within Amber ff14SB, OPLS-AA/M, and CHARMM36m reveal that only Amber ff14SB partially captures the ethanol-induced conformational changes of alanine residue in monomeric GAG when 42% ethanol is added to water. MD simulations of 200 mM GAG ensembles in pure water and in the aqueous solution with 42% ethanol showcase the ability of CHARMM36m to capture the effect of ethanol on the average pPII content of alanine in GAG and provide a plausible explanation for this effect, which may stem from an increased propensity of GAG to form oligomers in the presence of ethanol.

Formation of peptide-based oligomers in dimethylsulfoxide: identifying the precursor of fibril formation

The well-studied dipeptide fluorenylmethyloxycarbonyl-di-phenylalanine (FmocFF) forms a rigid hydrogel upon dissolving in dimethylsulfoxide (DMSO) and dilution in H2O. Here, we explored the pre-aggregation of the peptide in pure DMSO by vibrational spectroscopies, X-ray powder diffraction and dynamic light scattering. Our results show an equilibrium between a dominant population of amorphous oligomers (on a length scale of 2 nm) and a small number of protofibrils/fibrils (on a length scale of 30 nm in the centimolar and of 200 nm in the sub-molar region). To probe the mechanism underlying the formation of these protofilaments, we measured the 1H-NMR, IR and visible Raman spectra of DMSO containing different FmocFF concentrations, ranging between 10 and 300 mM. Our data reveal that interpeptide hydrogen bonding leads to the self-assembly of FmocFF in the centimolar region, while π-π stacking between Fmoc-groups is observed above 100 mM. The high 3J(HNHCα) coupling constant of the N-terminal amide proton indicates that the Fmoc end-cap of the peptide locks the N-terminal residue into a conformational ensemble centered at a φ-value of ca. -120°, which corresponds to a parallel β-sheet type conformation. The 3J(HNHCα) coupling constant of the C-terminal residue is indicative of a polyproline II (pPII)/βt mixture. Our results suggest that the gelation of FmocFF caused by the addition of a small amount of water to DMSO mixtures is facilitated by the formation of disordered protofibrils in pure DMSO.

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.

Do Molecular Dynamics Force Fields Capture Conformational Dynamics of Alanine in Water?

We examine the ability of six molecular dynamics (MD) force fields (Amber ff14SB, Amber ff99SBnmr1, Amber ff03ws, OPLS-AA/L, OPLS-AA/M, and CHARMM36) to reproduce conformational ensembles of the central alanine in GAG and AAA in a way that is consistent with five (GAG) or six (AAA) J coupling constants and amide I' profiles. MD-derived Ramachandran plots for all six force fields under study differ from those obtained by the Gaussian fit to experimental data in three major ways: (i) the polyproline II (pPII) basin in the Ramachandran plot is too concentrated, (ii) the antiparallel β (aβ) basin is overpopulated, and (iii) the transitional β (βt) basin is underpopulated. Amber ff14SB outperforms the other five MD force fields and yields the highest pPII populations of the central alanine residue in GAG (55%) and AAA (63%), in good agreement with the predictions of the Gaussian model (59 and 76%). The analysis of the hydration layer around the central alanine residue reveals considerable reorientation of water molecules and reduction in both the average number of water molecules and the average number of water-water hydrogen bonds when glycines (in GAG) are replaced by alanines (in AAA), elucidating water-mediated nearest neighbor effects on alanine's conformational dynamics.

A new interpretation of the structure and solvent dependence of the far UV circular dichroism spectrum of short oligopeptides

UV circular dichroism (UVCD) spectroscopy is a prominent tool for exploring secondary structures of polypeptides and proteins. In the unfolded state of these biomolecules, most of the individual residues primarily sample a conformation called polyproline II. Its CD spectrum contains a negatively biased positive couplet with a pronounced negative maximum below and a weak positive maximum above 200 nm. It is traditionally rationalized in terms of an excitonic coupling mechanism augmented by polarization effects. In this work, we carry out new time-dependent density functional theory calculations on the cationic tripeptide GAG in implicit and explicit water to determine the transitions that give rise to the observed CD signals of polyproline II and β-strand conformations. Our results reveal a plethora of electronic transitions that are governed by configurational interactions between multiple molecular orbital transitions of comparable energy. We also show that reproducing the CD spectra of polyproline II and β-strand conformations requires the explicit consideration of water molecules. The structure dependence of delocalized occupied orbitals contributes to the experimentally-observed invalidation of Flory's isolated pair hypothesis.

Amino acid conformations control the morphological and chiral features of the self-assembled peptide nanostructures: Young investigators perspective

HYPOTHESIS

A variety of nanostructures with different chiral features can be self-assembled from short peptides with highly similar sequences. We hypothesize that these supramolecular nanostructures are ruled by the constituent amino acid residues which adopt their conformations under the influence of intra-/inter-molecular interactions during peptide self-assembly.

APPROACH

Through reviewing recent advances in the self-assembly of short peptides and focusing on the relationship between amino acid conformations, peptide secondary structures and intra-/inter-molecular interactions within the supramolecular architectures, we aim to rationalize the complex interactive processes involved in the self-assembly of short, designed peptides.

RESULTS

Given the highly complexing interactive processes, the adoption of amino acid conformations and their control over peptide self-assembly consist of 4 main steps: (1) Each amino acid residue adopts its unique conformation in a specific sequence; (2) The sequence exhibits its own main chain geometry and determines the propensity of the intermolecular alignment within the building block; (3) The structural propensity of the building block and the packing mode between them determine the self-assembled structural features such as twisting, growth and chirality; (4) In addition to intra-/inter-molecular interactions, inter-sheet and inter-building block interactions could also affect the residue conformations and nanostructures, causing structural readjustment.

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.

Orientation of Oligopeptides in Self-Assembled Monolayers Inferred from Infrared Reflection-Absorption Spectroscopy

A series of a single tryptophan containing oligo-alanine peptides were recently characterized as conductive molecules that enable electron transport between electrodes. IR reflection-absorption of self-assembled monolayers of such peptides on gold surfaces revealed that the relative intensities of amide I and II bands in the respective spectra depend on the tryptophan residue position in the oligopeptide sequence. This indicates different average peptide orientations with respect to the normal onto the carrying gold surface. We developed a model which calculates the polarized reflectivities of the amide I and II bands as function of the angle of the incident light, the average peptide orientation and the relative orientations of peptide group at the N-terminal. The orientation and strength of vibrational transition dipole moments were calculated by employing an excitonic coupling approach which considers probable conformational distributions of the disordered peptides. Our results revealed that the position of the tryptophan can affect the effective tilt angle of the peptide as well as the orientation of transition dipole moments with respect to the reflection plane. We have also calculated the average end to end distances of the examined peptides and found them to be in reasonable agreement with experimental values determined by ellipsometry. Some evidence is obtained for the notion that increasing the tilt angle of the investigated peptides reduces their conductivity.

Hydration effects on Leu's polyproline II population in AcLXPNH2

Hydration is important in many fundamental processes. To investigate hydration effects on peptide conformations, we examined neighboring-residue and side-chain blocking effects in AcLXPNH2. A correlation between two effects suggests that hydration stabilizes PII more than β-structures. Our results are important for understanding the hydration effects on peptide conformations and hydration-forces in general.

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.

Prediction of nearest neighbor effects on backbone torsion angles and NMR scalar coupling constants in disordered proteins

Using fine-tuned hydrogen bonding criteria, a library of coiled peptide fragments has been generated from a large set of high-resolution protein X-ray structures. This library is shown to be an improved representation of ϕ/ψ torsion angles seen in intrinsically disordered proteins (IDPs). The ϕ/ψ torsion angle distribution of the library, on average, provides good agreement with experimentally observed chemical shifts and ³ J_HN-Hα coupling constants for a set of five disordered proteins. Inspection of the coil library confirms that nearest-neighbor effects significantly impact the ϕ/ψ distribution of residues in the coil state. Importantly, ³ J_HN-Hα coupling constants derived from the nearest-neighbor modulated backbone ϕ distribution in the coil library show improved agreement to experimental values, thereby providing a better way to predict ³ J_HN-Hα coupling constants for IDPs, and for identifying locations that deviate from fully random behavior.

Construction and comparison of the statistical coil states of unfolded and intrinsically disordered proteins from nearest-neighbor corrected conformational propensities of short peptides

Assessing the influence of nearest neighbors on the conformational ensemble of amino acid residues in unfolded and intrinsically disordered proteins and peptides is pivotal for a thorough understanding of the statistical coil state of unfolded proteins as well as of the energetics of the folding process. Research aimed at exploring nearest neighbor interactions has mostly focused on the analysis of restricted coil libraries that reflect conformational distributions in loops connecting more regular secondary structure segments. Recently, however, Toal et al. reported an experimentally based structural analysis of selected xy-pairs in GxyG tetrapeptides, which revealed quantitative information about conformational changes induced by nearest-neighbor interactions (Eur. J. Chem., 2015, 21, 5173-5192). Here, we perform analyses of Ramachandran plots of xy-pairs in GxyG and in coil libraries (Ting et al., PLOS CompBiol, 2010, 6, e1000763) using Hellinger distances as a quantitative measure of dissimilarities between Ramachandran distributions. Our analysis reveals that nearest-neighbor effects inferred from the above coil library are much less pronounced than corresponding structural changes observed for GxyG peptides. To determine whether nearest-neighbor induced conformational changes observed for GxyG can be utilized for the analysis of unfolded proteins, we analyzed sets of ³J(H^HH^α) coupling constants of three different unfolded proteins, namely the 130-residue fragment of the Staphylococcus aureus fibronectin-binding protein (FnBPc), denatured hen lysozyme, and the htau40 protein. For the first two proteins we found statistically meaningful correlations between predicted nearest-neighbor induced changes of ³J(H^HH^α) and experimentally observed deviations from corresponding coupling constants of GxG peptides in water, which we used as reference system with minimal nearest-neighbor interactions. This observation is in line with the NMR based understanding of these proteins being predominantly statistical coils. For htau40, however, which is known to exhibit residual structure and large deviations form statistical coil expectations, these correlations are weak or absent. Our results thus underscore the importance of nearest-neighbor interactions for a complete physical description of an ideal statistical coil state of a protein.

Dimension conversion and scaling of disordered protein chains

To extract protein dimension and energetics information from single-molecule fluorescence resonance energy transfer spectroscopy (smFRET) data, it is essential to establish the relationship between the distributions of the radius of gyration (Rg) and the end-to-end (donor-to-acceptor) distance (Ree). Here, we performed a coarse-grained molecular dynamics simulation to obtain a conformational ensemble of denatured proteins and intrinsically disordered proteins. For any disordered chain with fixed length, there is an excellent linear correlation between the average values of Rg and Ree under various solvent conditions, but the relationship deviates from the prediction of a Gaussian chain. A modified conversion formula was proposed to analyze smFRET data. The formula reduces the discrepancy between the results obtained from FRET and small-angle X-ray scattering (SAXS). The scaling law in a coil-globule transition process was examined where a significant finite-size effect was revealed, i.e., the scaling exponent may exceed the theoretical critical boundary [1/3, 3/5] and the prefactor changes notably during the transition. The Sanchez chain model was also tested and it was shown that the mean-field approximation works well for expanded chains.

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.

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.