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...

Descripción completa

Guardado en:
Detalles Bibliográficos
Autores principales: Juan Alcantara, Robyn Stix, Katherine Huang, Acadia Connor, Ray East, Valeria Jaramillo-Martinez, Elliott J. Stollar, K. Aurelia Ball
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