An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations



Alcantara, Juan, Stix, Robyn, Huang, Katherine, Connor, Acadia, East, Ray, Jaramillo-Martinez, Valeria, Stollar, Elliott J ORCID: 0000-0001-9895-0335 and Ball, K Aurelia
(2021) An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations. FRONTIERS IN MOLECULAR BIOSCIENCES, 8. 734169-.

Access the full-text of this item by clicking on the Open Access link.

Abstract

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 <i>cis</i> 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 <i>cis</i>, a robust methodology to accurately account for these conformations in the overall ensemble is crucial. Observing the <i>cis</i> 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 <i>cis</i> and <i>trans</i> 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 <i>cis</i> and <i>trans</i> states on a microsecond timescale. All ArkA prolines spend a significant fraction of time in <i>cis</i>, leading to a more compact ensemble with less polyproline II helix structure than an ArkA ensemble with all peptide bonds in <i>trans</i>. The ensemble containing <i>cis</i> prolines also matches more closely to <i>in vitro</i> circular dichroism data than the all-<i>trans</i> 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 <i>in vitro</i> data.

Item Type: Article
Uncontrolled Keywords: circular dichroism, intrinsically disordered proteins, isomerization, molecular dynamics, proline
Divisions: Faculty of Health and Life Sciences
Faculty of Health and Life Sciences > Institute of Systems, Molecular and Integrative Biology
Faculty of Health and Life Sciences > Institute of Systems, Molecular and Integrative Biology > School of Life Sciences
Depositing User: Symplectic Admin
Date Deposited: 07 Jun 2022 09:50
Last Modified: 18 Jan 2023 21:00
DOI: 10.3389/fmolb.2021.734169
Open Access URL: https://doi.org/10.3389/fmolb.2021.734169
Related URLs:
URI: https://livrepository.liverpool.ac.uk/id/eprint/3155995