End Effects Influence Short Model Peptide Conformation

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.

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.

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.

Bilayer thickness determines the alignment of model polyproline helices in lipid membranes

Our understanding of protein folds relies fundamentally on the set of secondary structures found in the proteomes. Yet, there also exist intriguing structures and motifs that are underrepresented in natural biopolymeric systems. One example is the polyproline II helix, which is usually considered to have a polar character and therefore does not form membrane spanning sections of membrane proteins. In our work, we have introduced specially designed polyproline II helices into the hydrophobic membrane milieu and used 19F NMR to monitor the helix alignment in oriented lipid bilayers. Our results show that these artificial hydrophobic peptides can adopt several different alignment states. If the helix is shorter than the thickness of the hydrophobic core of the membrane, it is submerged into the bilayer with its long axis parallel to the membrane plane. The polyproline helix adopts a transmembrane alignment when its length exceeds the bilayer thickness. If the peptide length roughly matches the lipid thickness, a coexistence of both states is observed. We thus show that the lipid thickness plays a determining role in the occurrence of a transmembrane polyproline II helix. We also found that the adaptation of polyproline II helices to hydrophobic mismatch is in some notable aspects different from α-helices. Finally, our results prove that the polyproline II helix is a competent structure for the construction of transmembrane peptide segments, despite the fact that no such motif has ever been reported in natural systems.

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.

Transmembrane Polyproline Helix

The third most abundant polypeptide conformation in nature, the polyproline-II helix, is a polar, extended secondary structure with a local organization stabilized by intercarbonyl interactions within the peptide chain. Here we design a hydrophobic polyproline-II helical peptide based on an oligomeric octahydroindole-2-carboxylic acid scaffold and demonstrate its transmembrane alignment in model lipid bilayers by means of solid-state ¹⁹F NMR. As result, we provide a first example of a purely artificial transmembrane peptide with a structural organization that is not based on hydrogen-bonding.

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.

Hydrolysis, polarity, and conformational impact of C-terminal partially fluorinated ethyl esters in peptide models

Fluorinated moieties are highly valuable to chemists due to the sensitive NMR detectability of the 19F nucleus. Fluorination of molecular scaffolds can also selectively influence a molecule's polarity, conformational preferences and chemical reactivity, properties that can be exploited for various chemical applications. A powerful route for incorporating fluorine atoms in biomolecules is last-stage fluorination of peptide scaffolds. One of these methods involves esterification of the C-terminus of peptides using a diazomethane species. Here, we provide an investigation of the physicochemical consequences of peptide esterification with partially fluorinated ethyl groups. Derivatives of N-acetylproline are used to model the effects of fluorination on the lipophilicity, hydrolytic stability and on conformational properties. The conformational impact of the 2,2-difluoromethyl ester on several neutral and charged oligopeptides was also investigated. Our results demonstrate that partially fluorinated esters undergo variable hydrolysis in biologically relevant buffers. The hydrolytic stability can be tailored over a broad pH range by varying the number of fluorine atoms in the ester moiety or by introducing adjacent charges in the peptide sequence.

Conformational Plasticity of the Cell-Penetrating Peptide SAP As Revealed by Solid-State 19F-NMR and Circular Dichroism Spectroscopies

The cell-penetrating peptide SAP, which was designed as an amphipathic poly-l-proline helix II (PPII), was suggested to self-assemble into regular fibrils that are relevant for its internalization. Herein we have analyzed the structure of SAP in the membrane-bound state by solid-state ¹⁹F-NMR, which revealed other structural states, in addition to the expected surface-aligned PPII. Trifluoromethyl-bicyclopentyl-glycine (CF₃-Bpg) and two rigid isomers of trifluoromethyl-4,5-methanoprolines (CF₃-MePro) were used as labels for ¹⁹F-NMR analysis. The equilibria between different conformations of SAP were studied and were found to be shifted by the substituents at Pro-11. Synchrotron-CD results suggested that substituting Pro-11 by CF₃-MePro governed the coil-to-PPII equilibrium in solution and in the presence of a lipid bilayer. Using CD and ¹⁹F-NMR, we examined the slow kinetics of the association of SAP with membranes and the dependence of the SAP conformational dynamics on the lipid composition. The peptide did not bind to lipids in the solid ordered phase and aggregated only in the liquid ordered "raft"-like bilayers. Self-association could not be detected in solution or in the presence of liquid disordered membranes. Surface-bound amphipathic SAP in a nonaggregated state was structured as a mixture of nonideal extended conformations reflecting the equilibrium already present in solution, i.e., before binding to the membrane.

Populations of the Minor α-Conformation in AcGXGNH2 and the α-Helical Nucleation Propensities

Intrinsic backbone conformational preferences of different amino acids are important for understanding the local structure of unfolded protein chains. Recent evidence suggests α-structure is relatively minor among three major backbone conformations for unfolded proteins. The α-helices are the dominant structures in many proteins. For these proteins, how could the α-structures occur from the least in unfolded to the most in folded states? Populations of the minor α-conformation in model peptides provide vital information. Reliable determination of populations of the α-conformers in these peptides that exist in multiple equilibriums of different conformations remains a challenge. Combined analyses on data from AcGXPNH₂ and AcGXGNH₂ peptides allow us to derive the populations of PII, β and α in AcGXGNH₂. Our results show that on average residue X in AcGXGNH₂ adopt PII, β, and α 44.7%, 44.5% and 10.8% of time, respectively. The contents of α-conformations for different amino acids define an α-helix nucleation propensity scale. With derived PII, β and α-contents, we can construct a free energy-conformation diagram on each AcGXGNH₂ in aqueous solution for the three major backbone conformations. Our results would have broad implications on early-stage events of protein folding.

Effects of phosphorylation on the intrinsic propensity of backbone conformations of serine/threonine

Each amino acid has its intrinsic propensity for certain local backbone conformations, which can be further modulated by the physicochemical environment and post-translational modifications. In this work, we study the effects of phosphorylation on the intrinsic propensity for different local backbone conformations of serine/threonine by molecular dynamics simulations. We showed that phosphorylation has very different effects on the intrinsic propensity for certain local backbone conformations for the serine and threonine. The phosphorylation of serine increases the propensity of forming polyproline II, whereas that of threonine has the opposite effect. Detailed analysis showed that such different responses to phosphorylation mainly arise from their different perturbations to the backbone hydration and the geometrical constraints by forming side-chain-backbone hydrogen bonds due to phosphorylation. Such an effect of phosphorylation on backbone conformations can be crucial for understanding the molecular mechanism of phosphorylation-regulated protein structures/dynamics and functions.

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.

Peptide Conformation Analysis Using an Integrated Bayesian Approach

Unlike native proteins that are amenable to structural analysis at atomic resolution, unfolded proteins occupy a manifold of dynamically interconverting structures. Defining the conformations of unfolded proteins is of significant interest and importance, for folding studies and for understanding the properties of intrinsically disordered proteins. Short chain protein fragments, i.e., oligopeptides, provide an excellent test-bed in efforts to define the conformational ensemble of unfolded chains. Oligomers of alanine in particular have been extensively studied as minimalist models of the intrinsic conformational preferences of the peptide backbone. Even short alanine peptides occupy an ensemble of substates that are distinguished by small free energy differences, so that the problem of quantifying the conformational preferences of the backbone remains a fundamental challenge in protein biophysics. Here, we demonstrate an integrated computational-experimental-Bayesian approach to quantify the conformational ensembles of the model trialanine peptide in water. In this approach, peptide conformational substates are first determined objectively by clustering molecular dynamics snapshots based on both structural and dynamic information. Next, a set of spectroscopic data for each conformational substate is computed. Finally, a Bayesian statistical analysis of both experimentally measured spectroscopic data and computational results is carried out to provide a current best estimate of the substate population ensemble together with corresponding confidence intervals. This distribution of substates can be further systematically refined with additional high-quality experimental data and more accurate computational modeling. Using an experimental data set of NMR coupling constants, we have also applied this approach to characterize the conformation ensemble of trivaline in water.

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.

Evaluating protocols and analytical methods for peptide adsorption experiments

This paper evaluates analytical techniques that are relevant for performing reliable quantitative analysis of peptide adsorption on surfaces. Two salient problems are addressed: determining the solution concentrations of model GG-X-GG, X5, and X10 oligopeptides (G = glycine, X = a natural amino acid), and quantitative analysis of these peptides following adsorption on surfaces. To establish a uniform methodology for measuring peptide concentrations in water across the entire GG-X-GG and X n series, three methods were assessed: UV spectroscopy of peptides having a C-terminal tyrosine, the bicinchoninic acid (BCA) protein assay, and amino acid (AA) analysis. Due to shortcomings or caveats associated with each of the different methods, none were effective at measuring concentrations across the entire range of representative model peptides. In general, reliable measurements were within 30% of the nominal concentration based on the weight of as-received lyophilized peptide. In quantitative analysis of model peptides adsorbed on surfaces, X-ray photoelectron spectroscopy (XPS) data for a series of lysine-based peptides (GGKGG, K5, and K10) on Au substrates, and for controls incubated in buffer in the absence of peptides, suggested a significant presence of aliphatic carbon species. Detailed analysis indicated that this carbonaceous contamination adsorbed from the atmosphere after the peptide deposition. The inferred adventitious nature of the observed aliphatic carbon was supported by control experiments in which substrates were sputter-cleaned by Ar(+) ions under ultra-high vacuum (UHV) then re-exposed to ambient air. In contrast to carbon contamination, no adventitious nitrogen species were detected on the controls; therefore, the relative surface densities of irreversibly-adsorbed peptides were calculated by normalizing the N/Au ratios by the average number of nitrogen atoms per residue.

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.

Disorder and order in unfolded and disordered peptides and proteins: a view derived from tripeptide conformational analysis. I. Tripeptides with long and predominantly hydrophobic side chains

We performed a conformational analysis of the central residues of three tripeptides glycyl-L-isoleucyl-glycine (GIG), glycyl-L-tyrosyl-glycine (GYG) and glycyl-L-arginyl-glycine (GRG) in aqueous solution, based on a global analysis of amide I' band profiles and NMR J-coupling constants. The results are compared with recently reported distributions of GVG, GFG and GEG. For GIG and GYG, we found that even though the polyproline II (pPII) fraction is below 0.5, it is still the most populated conformation, whereas GVG and GFG show both a larger β-strand fraction. For GRG, we observed a clear dominance of pPII over β-strand, reminiscent of observations for GEG and GKG. This finding indicates that terminal charges on otherwise hydrophobic residue side chains stabilize pPII over β-strand conformations. For all peptides investigated we found that a variety of compact and turn-like conformations constitute nearly 20 percent of their conformational distributions. Attempts to analyze our data with a simple two-state pPII-->/<--β model therefore do not yield any satisfactory reproduction of experimental results. A comparison of the obtained GxG ensembles with conformational distributions of GxG segments in truncated coil libraries (helices and sheets omitted) revealed a much larger fraction of type II β(i+2) and type III β like conformations for the latter. Thus, a comparison of conformational distributions of unfolded peptide segments in solution and in coil libraries reveal interesting information on how the interplay between intrinsic propensities of amino acid residues and non-local interactions in polypeptide chains determine the conformations of loop segments in proteins.