An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations
Disordered proline-rich motifs are common across the proteomes of many species and are often involved in protein-protein interactions. Proline is a unique amino acid due to the covalent bond between the backbone nitrogen and the proline side chain. The resulting five-membered ring allows proline to...
Guardado en:
Autores principales: | , , , , , , , |
---|---|
Formato: | article |
Lenguaje: | EN |
Publicado: |
Frontiers Media S.A.
2021
|
Materias: | |
Acceso en línea: | https://doaj.org/article/f1791d82aafe4d9dbee25ef021463a37 |
Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
id |
oai:doaj.org-article:f1791d82aafe4d9dbee25ef021463a37 |
---|---|
record_format |
dspace |
spelling |
oai:doaj.org-article:f1791d82aafe4d9dbee25ef021463a372021-11-15T06:44:24ZAn Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations2296-889X10.3389/fmolb.2021.734169https://doaj.org/article/f1791d82aafe4d9dbee25ef021463a372021-11-01T00:00:00Zhttps://www.frontiersin.org/articles/10.3389/fmolb.2021.734169/fullhttps://doaj.org/toc/2296-889XDisordered proline-rich motifs are common across the proteomes of many species and are often involved in protein-protein interactions. Proline is a unique amino acid due to the covalent bond between the backbone nitrogen and the proline side chain. The resulting five-membered ring allows proline to sample the cis state about its peptide bond, which other residues cannot do as readily. Because proline-rich disordered sequences exist as ensembles that likely include structures with the proline peptide bond in cis, a robust methodology to accurately account for these conformations in the overall ensemble is crucial. Observing the cis conformations of proline in a disordered sequence is challenging both experimentally and computationally. Nitrogen-hydrogen NMR spectroscopy cannot directly observe proline residues, which lack an amide bond, and computational methods struggle to overcome the large kinetic barrier between the cis and trans states, since isomerization usually occurs on the order of seconds. In the current work, Gaussian accelerated molecular dynamics was used to overcome this free energy barrier and simulate proline isomerization in a tetrapeptide (KPTP) and in the 12-residue proline-rich SH3 binding peptide, ArkA. We found that Gaussian accelerated molecular dynamics, when combined with a lowered peptide bond dihedral angle potential energy barrier (15 kcal/mol), allowed sufficient sampling of the proline cis and trans states on a microsecond timescale. All ArkA prolines spend a significant fraction of time in cis, leading to a more compact ensemble with less polyproline II helix structure than an ArkA ensemble with all peptide bonds in trans. The ensemble containing cis prolines also matches more closely to in vitro circular dichroism data than the all-trans ensemble. The ability of the ArkA prolines to isomerize likely affects the peptide’s ability to bind its partner SH3 domain, and should be studied further. This is the first molecular dynamics simulation study of proline isomerization in a biologically relevant proline-rich sequence that we know of, and a similar protocol could be applied to study multi-proline isomerization in other proline-containing proteins to improve conformational diversity and agreement with in vitro data.Juan AlcantaraRobyn StixKatherine HuangAcadia ConnorRay EastValeria Jaramillo-MartinezElliott J. StollarK. Aurelia BallFrontiers Media S.A.articleprolineisomerizationmolecular dynamicsintrinsically disordered proteinscircular dichroismBiology (General)QH301-705.5ENFrontiers in Molecular Biosciences, Vol 8 (2021) |
institution |
DOAJ |
collection |
DOAJ |
language |
EN |
topic |
proline isomerization molecular dynamics intrinsically disordered proteins circular dichroism Biology (General) QH301-705.5 |
spellingShingle |
proline isomerization molecular dynamics intrinsically disordered proteins circular dichroism Biology (General) QH301-705.5 Juan Alcantara Robyn Stix Katherine Huang Acadia Connor Ray East Valeria Jaramillo-Martinez Elliott J. Stollar K. Aurelia Ball An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations |
description |
Disordered proline-rich motifs are common across the proteomes of many species and are often involved in protein-protein interactions. Proline is a unique amino acid due to the covalent bond between the backbone nitrogen and the proline side chain. The resulting five-membered ring allows proline to sample the cis state about its peptide bond, which other residues cannot do as readily. Because proline-rich disordered sequences exist as ensembles that likely include structures with the proline peptide bond in cis, a robust methodology to accurately account for these conformations in the overall ensemble is crucial. Observing the cis conformations of proline in a disordered sequence is challenging both experimentally and computationally. Nitrogen-hydrogen NMR spectroscopy cannot directly observe proline residues, which lack an amide bond, and computational methods struggle to overcome the large kinetic barrier between the cis and trans states, since isomerization usually occurs on the order of seconds. In the current work, Gaussian accelerated molecular dynamics was used to overcome this free energy barrier and simulate proline isomerization in a tetrapeptide (KPTP) and in the 12-residue proline-rich SH3 binding peptide, ArkA. We found that Gaussian accelerated molecular dynamics, when combined with a lowered peptide bond dihedral angle potential energy barrier (15 kcal/mol), allowed sufficient sampling of the proline cis and trans states on a microsecond timescale. All ArkA prolines spend a significant fraction of time in cis, leading to a more compact ensemble with less polyproline II helix structure than an ArkA ensemble with all peptide bonds in trans. The ensemble containing cis prolines also matches more closely to in vitro circular dichroism data than the all-trans ensemble. The ability of the ArkA prolines to isomerize likely affects the peptide’s ability to bind its partner SH3 domain, and should be studied further. This is the first molecular dynamics simulation study of proline isomerization in a biologically relevant proline-rich sequence that we know of, and a similar protocol could be applied to study multi-proline isomerization in other proline-containing proteins to improve conformational diversity and agreement with in vitro data. |
format |
article |
author |
Juan Alcantara Robyn Stix Katherine Huang Acadia Connor Ray East Valeria Jaramillo-Martinez Elliott J. Stollar K. Aurelia Ball |
author_facet |
Juan Alcantara Robyn Stix Katherine Huang Acadia Connor Ray East Valeria Jaramillo-Martinez Elliott J. Stollar K. Aurelia Ball |
author_sort |
Juan Alcantara |
title |
An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations |
title_short |
An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations |
title_full |
An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations |
title_fullStr |
An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations |
title_full_unstemmed |
An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations |
title_sort |
unbound proline-rich signaling peptide frequently samples cis conformations in gaussian accelerated molecular dynamics simulations |
publisher |
Frontiers Media S.A. |
publishDate |
2021 |
url |
https://doaj.org/article/f1791d82aafe4d9dbee25ef021463a37 |
work_keys_str_mv |
AT juanalcantara anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT robynstix anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT katherinehuang anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT acadiaconnor anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT rayeast anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT valeriajaramillomartinez anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT elliottjstollar anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT kaureliaball anunboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT juanalcantara unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT robynstix unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT katherinehuang unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT acadiaconnor unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT rayeast unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT valeriajaramillomartinez unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT elliottjstollar unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations AT kaureliaball unboundprolinerichsignalingpeptidefrequentlysamplescisconformationsingaussianacceleratedmoleculardynamicssimulations |
_version_ |
1718428592387391488 |