skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Glycine in Water Favors the Polyproline II State
Conformational preferences of amino acid residues in water are determined by the backbone and side-chain properties. Alanine is known for its high polyproline II (pPII) propensity. The question of relative contributions of the backbone and side chain to the conformational preferences of alanine and other amino acid residues in water is not fully resolved. Because glycine lacks a heavy-atom side chain, glycine-based peptides can be used to examine to which extent the backbone properties affect the conformational space. Here, we use published spectroscopic data for the central glycine residue of cationic triglycine in water to demonstrate that its conformational space is dominated by the pPII state. We assess three commonly used molecular dynamics (MD) force fields with respect to their ability to capture the conformational preferences of the central glycine residue in triglycine. We show that pPII is the mesostate that enables the functional backbone groups of the central residue to form the most hydrogen bonds with water. Our results indicate that the pPII propensity of the central glycine in GGG is comparable to that of alanine in GAG, implying that the water-backbone hydrogen bonding is responsible for the high pPII content of these residues.  more » « less
Award ID(s):
1817650
PAR ID:
10251353
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Biomolecules
Volume:
10
Issue:
8
ISSN:
2218-273X
Page Range / eLocation ID:
1121
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Cysteine sulfonic acid (Cys-SO3H; cysteic acid) is an oxidative post-translational modification of cysteine, resulting from further oxidation from cysteine sulfinic acid (Cys- SO2H). Cysteine sulfonic acid is considered an irreversible post-translational modification, which serves as a biomarker of oxidative stress that has resulted in oxidative damage to proteins. Cysteine sulfonic acid is anionic, as a sulfonate (Cys-SO –; cysteate), in the ionization state that 3 is almost exclusively present at physiological pH (pKa ~ –2). In order to understand protein structural changes that can occur upon oxidation to cysteine sulfonic acid, we analyzed its conformational preferences, using experimental methods, bioinformatics, and DFT-based computational analysis. Cysteine sulfonic acid was incorporated into model peptides for α-helix and polyproline II helix (PPII). Within peptides, oxidation of cysteine to the sulfonic acid proceeds rapidly and efficiently at room temperature in solution with methyltrioxorhenium (MeReO3) and H2O2. Peptides containing cysteine sulfonic acid were also generated on solid phase using trityl-protected cysteine and oxidation with MeReO3 and H2O2. Using methoxybenzyl (Mob)-protected cysteine, solid-phase oxidation with MeReO3 and H2O2 generated the Mob sulfone precursor to Cys-SO – within fully synthesized peptides. These two solid-phase methods allow the synthesis of peptides containing either Cys-SO – or Cys-SO – in a 32 practical manner, with no solution-phase synthesis required. Cys-SO – had low PPII propensity 3 for PPII propagation, despite promoting a relatively compact conformation in φ. In contrast, in a PPII initiation model system, Cys-SO – promoted PPII relative to neutral Cys, with PPII initiation similar to Cys thiolate but less than Cys-SO – or Ala. In an α-helix model system, Cys- 2 SO – promoted α-helix near the N-terminus, due to favorable helix dipole interactions and 3 favorable α-helix capping via a sulfonate-amide side chain-main chain hydrogen bond. Across all peptides, the sulfonate side chain was significantly less ordered than that of the sulfinate. Analysis of Cys-SO – in the PDB revealed a very strong propensity for local (i/i or i/i+1) side 3 chain-main chain sulfonate-amide hydrogen bonds for Cys-SO –, with > 80% of Cys-SO – 33 residues exhibiting these interactions. DFT calculations conducted to explore these conformational preferences indicated that side chain-main chain hydrogen bonds of the sulfonate with the intraresidue amide and/or with the i+1 amide were favorable. However, hydrogen bonds to water or to amides, as well as interactions with oxophilic metals, were weaker for the sulfonate than the sulfinate, due to lower charge density on the oxygens in the sulfonate. 
    more » « less
  2. 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. 
    more » « less
  3. 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 3J(HNHCα) 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. 
    more » « less
  4. Molecular dynamics (MD) is a powerful tool for studying intrinsically disordered proteins, however, its reliability depends on the accuracy of the force field. We assess Amber ff19SB, Amber ff14SB, OPLS-AA/M, and CHARMM36m with respect to their capacity to capture intrinsic conformational dynamics of 14 guest residues x (=G, A, L, V, I, F, Y, D P , E P , R, C, N, S, T) in GxG peptides in water. The MD-derived Ramachandran distribution of each guest residue is used to calculate 5 J-coupling constants and amide I′ band profiles to facilitate a comparison to spectroscopic data through reduced χ 2 functions. We show that the Gaussian model, optimized to best fit the experimental data, outperforms all MD force fields by an order of magnitude. The weaknesses of the MD force fields are: (i) insufficient variability of the polyproline II (pPII) population among the guest residues; (ii) oversampling of antiparallel at the expense of transitional β-strand region; (iii) inadequate sampling of turn-forming conformations for ionizable and polar residues; and (iv) insufficient guest residue-specificity of the Ramachandran distributions. Whereas Amber ff19SB performs worse than the other three force fields with respect to χ 2 values, it accounts for residue-specific pPII content better than the other three force fields. Additional testing of residue-specific RSFF1 and Amber ff14SB combined with TIP4P/2005 on six guest residues x (=A, I, F, D P , R, S) reveals that residue specificity derived from protein coil libraries or an improved water model alone do not result in significantly lower χ 2 values. 
    more » « less
  5. Amphipathic peptides with amino acids arranged in alternating patterns of hydrophobic and hydrophilic residues efficiently self‐assemble intoβ‐sheet bilayer nanoribbons. Hydrophobic side chain functionality is effectively buried in the interior of the putative bilayer of these nanoribbons. This study investigates consequences on self‐assembly of increasing the surface area of aromatic side chain groups that reside in the hydrophobic core of nanoribbons derived from Ac‐(XKXE)2‐NH2peptides (X = hydrophobic residue). A series of Ac‐(XKXE)2‐NH2peptides incorporating aromatic amino acids of increasing molecular volume and steric profile (X = phenylalanine [Phe], homophenylalanine [Hph], tryptophan [Trp], 1‐naphthylalanine [1‐Nal], 2‐naphthylalanine [2‐Nal], or biphenylalanine [Bip]) were assessed to determine substitution effects on self‐assembly propensity and on morphology of the resulting nanoribbon structures. Additional studies were conducted to determine the effects of incorporating amino acids of differing steric profile in the hydrophobic core (Ac‐X1KFEFKFE‐NH2and Ac‐(X1,5KFE)‐NH2peptides, X = Trp or Bip). Spectroscopic analysis by circular dichroism (CD) and Fourier transform infrared (FT‐IR) spectroscopy indicatedβ‐sheet formation for all variants. Self‐assembly rate increased with peptide hydrophobicity; increased molecular volume of the hydrophobic side chain groups did not appear to induce kinetic penalties on self‐assembly rates. Transmission electron microscopy (TEM) imaging indicated variation in fibril morphology as a function of amino acid in the X positions. This study confirms that hydrophobicity of amphipathic Ac‐(XKXE)2‐NH2peptides correlates to self‐assembly propensity and that the hydrophobic core of the resulting nanoribbon bilayers has a significant capacity to accommodate sterically demanding functional groups. These findings provide insight that may be used to guide the exploitation of self‐assembled amphipathic peptides as functional biomaterials. 
    more » « less