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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Do molecular dynamics force fields accurately model Ramachandran distributions of amino acid residues in water?
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
Award ID(s):
1817650
PAR ID:
10335239
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
24
Issue:
5
ISSN:
1463-9076
Page Range / eLocation ID:
3259 to 3279
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Biomolecules)
    null (Ed.)
    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
  2. ; (, Chemistry – A European Journal)
    Abstract The type II polyproline helix (PPII) is a fundamental secondary structure of proteins, important in globular proteins, in intrinsically disordered proteins, and at protein‐protein interfaces. PPII is stabilized in part byn→π* interactions between consecutive carbonyls, via electron delocalization between an electron‐donor carbonyl lone pair (n) and an electron‐acceptor carbonyl (π*) on the subsequent residue. We previously demonstrated that changes to the electronic properties of the acyl donor can predictably modulate the strength ofn→π* interactions, with data from model compounds, in solution in chloroform, in the solid state, and computationally. Herein, we examined whether the electronic properties of acyl capping groups could modulate the stability of PPII in peptides in water. InX−PPGY‐NH2peptides (X=10 acyl capping groups), the effect of acyl group identity on PPII was quantified by circular dichroism and NMR spectroscopy. Electron‐rich acyl groups promoted PPII relative to the standard acetyl (Ac−) group, with the pivaloyl andiso‐butyryl groups most significantly increasing PPII. In contrast, acyl derivatives with electron‐withdrawing substituents and the formyl group relatively disfavored PPII. Similar results, though lesser in magnitude, were also observed inX−APPGY‐NH2peptides, indicating that the capping group can impact PPII conformation at both proline and non‐proline residues. The pivaloyl group was particularly favorable in promoting PPII. The effects of acyl capping groups were further analyzed inX–DfpPGY‐NH2andX−ADfpPGY‐NH2peptides, Dfp=4,4‐difluoroproline. Data on these peptides indicated that acyl groups induced order Piv‐ > Ac‐ > For‐. These results suggest that greater consideration should be given to the identity of acyl capping groups in inducing structure in peptides. 
    more » « less
  3. ; ; ; ; ; ; (, MRS Communications)
    Abstract Polyproline II (PPII) peptide sequences are recognized as promising biomaterials because of their attractive antifouling properties. However, the mechanisms behind their antifouling behavior have not been fully characterized. In this work we show that PPII peptide coverage, controlled by adsorption time, significantly reduces the fouling of bovine serum albumin (BSA, a model foulant). In addition, guest residues introduced into the PPII sequence are shown to significantly impact BSA adsorption as well as human mesenchymal stem cell (hMSC) spreading. This research will help guide future PPII peptide designs for incorporation into novel biomaterials. Graphical abstract 
    more » « less
  4. ; ; (, Proteins: Structure, Function, and Bioinformatics)
    Abstract Structures at serine‐proline sites in proteins were analyzed using a combination of peptide synthesis with structural methods and bioinformatics analysis of the PDB. Dipeptides were synthesized with the proline derivative (2S,4S)‐(4‐iodophenyl)hydroxyproline [hyp(4‐I‐Ph)]. The crystal structure of Boc‐Ser‐hyp(4‐I‐Ph)‐OMe had two molecules in the unit cell. One molecule exhibitedcis‐proline and a type VIa2 β‐turn (BcisD). Thecis‐proline conformation was stabilized by a C–H/O interaction between Pro C–Hαand the Ser side‐chain oxygen. NMR data were consistent with stabilization ofcis‐proline by a C–H/O interaction in solution. The other crystallographically observed molecule hadtrans‐Pro and both residues in the PPII conformation. Two conformations were observed in the crystal structure of Ac‐Ser‐hyp(4‐I‐Ph)‐OMe, with Ser adopting PPII in one and the β conformation in the other, each with Pro in the δ conformation andtrans‐Pro. Structures at Ser‐Pro sequences were further examined via bioinformatics analysis of the PDB and via DFT calculations. Ser‐Pro versus Ala–Pro sequences were compared to identify bases for Ser stabilization of local structures. C–H/O interactions between the Ser side‐chain Oγand Pro C–Hαwere observed in 45% of structures with Ser‐cis‐Pro in the PDB, with nearly all Ser‐cis‐Pro structures adopting a type VI β‐turn. 53% of Ser‐trans‐Pro sequences exhibited main‐chain COi•••HNi+3or COi•••HNi+4hydrogen bonds, with Ser as theiresidue and Pro as thei + 1 residue. These structures were overwhelmingly either type I β‐turns or N‐terminal capping motifs on α‐helices or 310‐helices. These results indicate that Ser‐Pro sequences are particularly potent in favoring these structures. In each, Ser is in either the PPII or β conformation, with the Ser Oγcapable of engaging in a hydrogen bond with the amide N–H of thei + 2 (type I β‐turn or 310‐helix; Serχ1t) ori + 3 (α‐helix; Serχ1g+) residue. Non‐prolinecisamide bonds can also be stabilized by C–H/O interactions. 
    more » « less
  5. ; (, International Journal of Molecular Sciences)
    The interaction between bone morphogenetic protein-2 (BMP-2) and the surface of biomaterials is essential for the restoration of bone and cartilage tissue, inducing cellular differentiation and proliferation. The properties of the surface, including topology features, regulate the conformation and bioactivity of the protein. In this research, we investigated the influence of nanopatterned surfaces on the interaction of a homodimer BMP-2 with graphite material by combining molecular dynamics (MD) and steered molecular dynamics (SMD) simulations. The graphite substrates were patterned as flat, linear grating, square, and circular profiles in combination with BMP-2 conformation in the side-on configuration. Ramachandran plots for the wrist and knuckle epitopes indicated no steric hindrances and provided binding sites to type I and type II receptors. Results showed two optimal patterns that increased protein adsorption of the lower monomer while preserving the secondary structure and leaving the upper monomer free to interact with the cells. Charged residues arginine and lysine and polar residues histidine and tyrosine were the main residues responsible for the strong interaction with the graphite surface. This research provides new molecular-level insights to further understand the mechanisms underlying protein adsorption on nanoscale patterned substrates. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email