1
|
Ali AAAI, Hoffmann F, Schäfer LV, Mulder FAA. Probing Methyl Group Dynamics in Proteins by NMR Cross-Correlated Dipolar Relaxation and Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:7722-7732. [PMID: 36326619 DOI: 10.1021/acs.jctc.2c00568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Nuclear magnetic resonance (NMR) spin relaxation is the most informative approach to experimentally probe the internal dynamics of proteins on the picosecond to nanosecond time scale. At the same time, molecular dynamics (MD) simulations of biological macromolecules are steadily improving through better physical models, enhanced sampling methods, and increased computational power, and they provide exquisite information about flexibility and its role in protein stability and molecular interactions. Many examples have shown that MD is now adept in probing protein backbone motion, but improvements are still required toward a quantitative description of the dynamics of side chains, for example, probed by the dynamics of methyl groups. Thus far, the comparison of computation with experiment for side chain dynamics has primarily focused on the relaxation of 13C and 2H nuclei induced by autocorrelated variation of spin interactions. However, the cross-correlation of 13C-1H dipolar interactions in methyl groups offers an attractive alternative. Here, we establish a computational framework to extract cross-correlation relaxation parameters of methyl groups in proteins from all-atom MD simulations. To demonstrate the utility of the approach, cross-correlation relaxation rates of ubiquitin are computed from MD simulations performed with the AMBER99SB*-ILDN and CHARMM36 force fields. Simulation results were found to agree well with those obtained by experiment. Moreover, the data obtained with the two force fields are highly consistent.
Collapse
Affiliation(s)
- Ahmed A A I Ali
- Theoretical Chemistry, Ruhr University Bochum, D-44780Bochum, Germany
| | - Falk Hoffmann
- Theoretical Chemistry, Ruhr University Bochum, D-44780Bochum, Germany
| | - Lars V Schäfer
- Theoretical Chemistry, Ruhr University Bochum, D-44780Bochum, Germany
| | - Frans A A Mulder
- Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry, Aarhus University, DK-8000Aarhus, Denmark
| |
Collapse
|
2
|
Wang J, Miao Y. Peptide Gaussian accelerated molecular dynamics (Pep-GaMD): Enhanced sampling and free energy and kinetics calculations of peptide binding. J Chem Phys 2021; 153:154109. [PMID: 33092378 DOI: 10.1063/5.0021399] [Citation(s) in RCA: 61] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play an important role in cellular signaling. However, it is challenging to simulate both binding and unbinding of peptides and calculate peptide binding free energies through conventional molecular dynamics, due to long biological timescales and extremely high flexibility of the peptides. Based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique, we have developed a new computational method "Pep-GaMD," which selectively boosts essential potential energy of the peptide in order to effectively model its high flexibility. In addition, another boost potential is applied to the remaining potential energy of the entire system in a dual-boost algorithm. Pep-GaMD has been demonstrated on binding of three model peptides to the SH3 domains. Independent 1 µs dual-boost Pep-GaMD simulations have captured repetitive peptide dissociation and binding events, which enable us to calculate peptide binding thermodynamics and kinetics. The calculated binding free energies and kinetic rate constants agreed very well with available experimental data. Furthermore, the all-atom Pep-GaMD simulations have provided important insights into the mechanism of peptide binding to proteins that involves long-range electrostatic interactions and mainly conformational selection. In summary, Pep-GaMD provides a highly efficient, easy-to-use approach for unconstrained enhanced sampling and calculations of peptide binding free energies and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| |
Collapse
|
3
|
Sheedlo MJ, Kenny S, Podkorytov IS, Brown K, Ma J, Iyer S, Hewitt CS, Arbough T, Mikhailovskii O, Flaherty DP, Wilson MA, Skrynnikov NR, Das C. Insights into Ubiquitin Product Release in Hydrolysis Catalyzed by the Bacterial Deubiquitinase SdeA. Biochemistry 2021; 60:584-596. [PMID: 33583181 DOI: 10.1021/acs.biochem.0c00760] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report the co-crystal structure of the (catalytic Cys)-to-Ala mutant of the deubiquitinase domain of the Legionella pneumophila effector SdeA (SdeADUB) with its ubiquitin (Ub) product. Most of the intermolecular interactions are preserved in this product-bound structure compared to that of the previously characterized complex of SdeADUB with the suicide inhibitor ubiquitin vinylmethyl ester (Ub-VME), whose structure models the acyl-enzyme thioester intermediate. Nuclear magnetic resonance (NMR) titration studies show a chemical shift perturbation pattern that suggests that the same interactions also exist in solution. Isothermal titration calorimetry and NMR titration data reveal that the affinity of wild-type (WT) SdeADUB for Ub is significantly lower than that of the Cys-to-Ala mutant. This is potentially due to repulsive interaction between the thiolate ion of the catalytic Cys residue in WT SdeADUB and the carboxylate group of the C-terminal Gly76 residue in Ub. In the context of SdeADUB catalysis, this electrostatic repulsion arises after the hydrolysis of the scissile isopeptide bond in the acyl-enzyme intermediate and the consequent formation of the C-terminal carboxylic group in the Ub fragment. We hypothesize that this electrostatic repulsion may expedite the release of the Ub product by SdeADUB. We note that similar repulsive interactions may also occur in other deubiquitinases and hydrolases of ubiquitin-like protein modifiers and may constitute a fairly general mechanism of product release within this family. This is a potentially important feature for a family of enzymes that form extensive protein-protein interactions during enzyme-substrate engagement.
Collapse
Affiliation(s)
- Michael J Sheedlo
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Sebastian Kenny
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ivan S Podkorytov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg 199034, Russia
| | - Kwame Brown
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Jia Ma
- Bindley Bioscience Center, West Lafayette, Indiana 47906, United States
| | - Shalini Iyer
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Chad S Hewitt
- Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, Indiana 47906, United States
| | - Trent Arbough
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Oleg Mikhailovskii
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States.,Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg 199034, Russia
| | - Daniel P Flaherty
- Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, Indiana 47906, United States
| | - Mark A Wilson
- Department of Biochemistry and Redox Biology Center, University of Nebraska, Lincoln, Nebraska 68588, United States
| | - Nikolai R Skrynnikov
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States.,Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg 199034, Russia
| | - Chittaranjan Das
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
4
|
Girodat D, Pati AK, Terry DS, Blanchard SC, Sanbonmatsu KY. Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein. PLoS Comput Biol 2020; 16:e1008293. [PMID: 33151943 PMCID: PMC7643941 DOI: 10.1371/journal.pcbi.1008293] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Accepted: 08/27/2020] [Indexed: 12/15/2022] Open
Abstract
Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BPSS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BPSS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BPSS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BPSS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R0) and fluorophore orientation factor (κ2) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET. Förster (fluorescence) resonance energy transfer (FRET) has been used extensively by biophysicists as a molecular-scale ruler that yields fundamental structural and kinetic insights into transient processes including complex formation and conformational rearrangements required for biological function. FRET techniques require the identification of informative fluorophore labeling sites, spaced at defined distances to inform on a reaction coordinate of interest and consideration of noise sources that have the potential to obscure quantitative interpretations. Here, we describe an approach to leverage advancements in computationally efficient all-atom structure-based molecular dynamics simulations in which structural dynamics observed via FRET can be interpreted in full atomistic detail on commensurate time scales. We demonstrate the potential of this approach using a model FRET system, the amino acid binding protein LIV-BPSS conjugated to self-healing organic fluorophores. LIV-BPSS exhibits large scale, sub-millisecond clamshell-like conformational changes between open and closed conformations associated with ligand unbinding and binding, respectively. Our findings inform on the molecular basis of the dynamics observed by smFRET and on strategies to optimize fluorophore labeling sites, the manner of fluorophore attachment, and fluorophore composition.
Collapse
Affiliation(s)
- Dylan Girodat
- Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | - Avik K Pati
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee, United States of America
| | - Daniel S Terry
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee, United States of America
| | - Scott C Blanchard
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee, United States of America
| | - Karissa Y Sanbonmatsu
- Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America.,New Mexico Consortium, Los Alamos, New Mexico, United States of America
| |
Collapse
|
5
|
Gerlach GJ, Carrock R, Stix R, Stollar EJ, Ball KA. A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway. PLoS Comput Biol 2020; 16:e1007815. [PMID: 32925900 PMCID: PMC7514057 DOI: 10.1371/journal.pcbi.1007815] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Revised: 09/24/2020] [Accepted: 08/15/2020] [Indexed: 12/20/2022] Open
Abstract
Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.
Collapse
Affiliation(s)
- Gabriella J. Gerlach
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| | - Rachel Carrock
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| | - Robyn Stix
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| | - Elliott J. Stollar
- School of Life Sciences, University of Liverpool, Liverpool, United Kingdom
| | - K. Aurelia Ball
- Department of Chemistry, Skidmore College, Saratoga Springs, New York, United States
| |
Collapse
|
6
|
Ding W, Tan HY, Zhang JX, Wilczek LA, Hsieh KR, Mulkin JA, Bianco PR. The mechanism of Single strand binding protein-RecG binding: Implications for SSB interactome function. Protein Sci 2020; 29:1211-1227. [PMID: 32196797 PMCID: PMC7184773 DOI: 10.1002/pro.3855] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2020] [Revised: 03/11/2020] [Accepted: 03/13/2020] [Indexed: 01/10/2023]
Abstract
The Escherichia coli single-strand DNA binding protein (SSB) is essential to viability where it functions to regulate SSB interactome function. Here it binds to single-stranded DNA and to target proteins that comprise the interactome. The region of SSB that links these two essential protein functions is the intrinsically disordered linker. Key to linker function is the presence of three, conserved PXXP motifs that mediate binding to oligosaccharide-oligonucleotide binding folds (OB-fold) present in SSB and its interactome partners. Not surprisingly, partner OB-fold deletions eliminate SSB binding. Furthermore, single point mutations in either the PXXP motifs or, in the RecG OB-fold, obliterate SSB binding. The data also demonstrate that, and in contrast to the view currently held in the field, the C-terminal acidic tip of SSB is not required for interactome partner binding. Instead, we propose the tip has two roles. First, and consistent with the proposal of Dixon, to regulate the structure of the C-terminal domain in a biologically active conformation that prevents linkers from binding to SSB OB-folds until this interaction is required. Second, as a secondary binding domain. Finally, as OB-folds are present in SSB and many of its partners, we present the SSB interactome as the first family of OB-fold genome guardians identified in prokaryotes.
Collapse
Affiliation(s)
- Wenfei Ding
- Center for Single Molecule BiophysicsUniversity at BuffaloBuffaloNew YorkUnited States
- Department of BiochemistryUniversity at BuffaloBuffaloNew YorkUnited States
| | - Hui Yin Tan
- Center for Single Molecule BiophysicsUniversity at BuffaloBuffaloNew YorkUnited States
- Present address:
Department of Chemistry and BiochemistryUniversity of Notre DameSouth BendIndianaUnited States
| | - Jia Xiang Zhang
- Department of BiochemistryUniversity at BuffaloBuffaloNew YorkUnited States
| | - Luke A. Wilczek
- Center for Single Molecule BiophysicsUniversity at BuffaloBuffaloNew YorkUnited States
- Department of BiochemistryUniversity at BuffaloBuffaloNew YorkUnited States
- Present address:
Department of ChemistryBrown UniversityProvidenceRhode IslandUnited States
| | - Karin R. Hsieh
- Center for Single Molecule BiophysicsUniversity at BuffaloBuffaloNew YorkUnited States
| | - Jeffrey A. Mulkin
- Center for Single Molecule BiophysicsUniversity at BuffaloBuffaloNew YorkUnited States
| | - Piero R. Bianco
- Center for Single Molecule BiophysicsUniversity at BuffaloBuffaloNew YorkUnited States
- Department of BiochemistryUniversity at BuffaloBuffaloNew YorkUnited States
| |
Collapse
|
7
|
Stadmiller SS, Aguilar JS, Waudby CA, Pielak GJ. Rapid Quantification of Protein-Ligand Binding via 19F NMR Lineshape Analysis. Biophys J 2020; 118:2537-2548. [PMID: 32348722 DOI: 10.1016/j.bpj.2020.03.031] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2020] [Accepted: 03/19/2020] [Indexed: 12/14/2022] Open
Abstract
Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional 19F NMR lineshape analysis to quantify the kinetics and equilibrium thermodynamics for the binding of a fluorine-labeled Src homology 3 (SH3) protein domain to four proline-rich peptides. SH3 domains are one of the largest and most well-characterized families of protein recognition domains and have a multitude of functions in eukaryotic cell signaling. First, we showe that fluorine incorporation into SH3 causes only minor structural changes to both the free and bound states using amide proton temperature coefficients. We then compare the results from lineshape analysis of one-dimensional 19F spectra to those from two-dimensional 1H-15N heteronuclear single quantum coherence spectra. Their agreement demonstrates that one-dimensional 19F lineshape analysis is a robust, low-cost, and fast alternative to traditional heteronuclear single quantum coherence-based experiments. The data show that binding is diffusion limited and indicate that the transition state is highly similar to the free state. We also measured binding as a function of temperature. At equilibrium, binding is enthalpically driven and arises from a highly positive activation enthalpy for association with small entropic contributions. Our results agree with those from studies using different techniques, providing additional evidence for the utility of 19F NMR lineshape analysis, and we anticipate that this analysis will be an effective tool for rapidly characterizing the energetics of protein interactions.
Collapse
Affiliation(s)
| | - Jhoan S Aguilar
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina
| | - Christopher A Waudby
- Department of Structural and Molecular Biology, University College London, London, United Kingdom
| | - Gary J Pielak
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina; Department of Biochemistry and Biophysics, University of North Carolina, Chapel Hill, North Carolina; Lineberger Comprehensive Cancer Center, University of North Carolina, Chapel Hill, North Carolina; Integrative Program for Biological and Genome Sciences, University of North Carolina, Chapel Hill, North Carolina.
| |
Collapse
|
8
|
Ozmaian M, Makarov DE. Transition path dynamics in the binding of intrinsically disordered proteins: A simulation study. J Chem Phys 2019; 151:235101. [PMID: 31864244 DOI: 10.1063/1.5129150] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Association of proteins and other biopolymers is a ubiquitous process in living systems. Recent single-molecule measurements probe the dynamics of association in unprecedented detail by measuring the properties of association transition paths, i.e., short segments of molecular trajectories between the time the proteins are close enough to interact and the formation of the final complex. Interpretation of such measurements requires adequate models for describing the dynamics of experimental observables. In an effort to develop such models, here we report a simulation study of the association dynamics of two oppositely charged, disordered polymers. We mimic experimental measurements by monitoring intermonomer distances, which we treat as "experimental reaction coordinates." While the dynamics of the distance between the centers of mass of the molecules is found to be memoryless and diffusive, the dynamics of the experimental reaction coordinates displays significant memory and can be described by a generalized Langevin equation with a memory kernel. We compute the most commonly measured property of transition paths, the distribution of the transition path time, and show that, despite the non-Markovianity of the underlying dynamics, it is well approximated as one-dimensional diffusion in the potential of mean force provided that an apparent value of the diffusion coefficient is used. This apparent value is intermediate between the slow (low frequency) and fast (high frequency) limits of the memory kernel. We have further studied how the mean transition path time depends on the ionic strength and found only weak dependence despite strong electrostatic attraction between the polymers.
Collapse
Affiliation(s)
- Masoumeh Ozmaian
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| | - Dmitrii E Makarov
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
9
|
Kämpf K, Izmailov SA, Rabdano SO, Groves AT, Podkorytov IS, Skrynnikov NR. What Drives 15N Spin Relaxation in Disordered Proteins? Combined NMR/MD Study of the H4 Histone Tail. Biophys J 2018; 115:2348-2367. [PMID: 30527335 DOI: 10.1016/j.bpj.2018.11.017] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 11/07/2018] [Accepted: 11/12/2018] [Indexed: 12/26/2022] Open
Abstract
Backbone (15N) NMR relaxation is one of the main sources of information on dynamics of disordered proteins. Yet, we do not know very well what drives 15N relaxation in such systems, i.e., how different forms of motion contribute to the measurable relaxation rates. To address this problem, we have investigated, both experimentally and via molecular dynamics simulations, the dynamics of a 26-residue peptide imitating the N-terminal portion of the histone protein H4. One part of the peptide was found to be fully flexible, whereas the other part features some transient structure (a hairpin stabilized by hydrogen bonds). The following motional modes proved relevant for 15N relaxation. 1) Sub-picosecond librations attenuate relaxation rates according to S2 ∼0.85-0.90. 2) Axial peptide-plane fluctuations along a stretch of the peptide chain contribute to relaxation-active dynamics on a fast timescale (from tens to hundreds of picoseconds). 3) φ/ψ backbone jumps contribute to relaxation-active dynamics on both fast (from tens to hundreds of picoseconds) and slow (from hundreds of picoseconds to a nanosecond) timescales. The major contribution is from polyproline II (PPII) ↔ β transitions in the Ramachandran space; in the case of glycine residues, the major contribution is from PPII ↔ (β) ↔ rPPII transitions, in which rPPII is the mirror-image (right-handed) version of the PPII geometry, whereas β geometry plays the role of an intermediate state. 4) Reorientational motion of certain (sufficiently long-lived) elements of transient structure, i.e., rotational tumbling, contributes to slow relaxation-active dynamics on ∼1-ns timescale (however, it is difficult to isolate this contribution). In conclusion, recent advances in the area of force-field development have made it possible to obtain viable Molecular Dynamics models of protein disorder. After careful validation against the experimental relaxation data, these models can provide a valuable insight into mechanistic origins of spin relaxation in disordered peptides and proteins.
Collapse
Affiliation(s)
- Kerstin Kämpf
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Sergei A Izmailov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Sevastyan O Rabdano
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Adam T Groves
- Department of Chemistry, Purdue University, West Lafayette, Indiana
| | - Ivan S Podkorytov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia
| | - Nikolai R Skrynnikov
- Laboratory of Biomolecular NMR, St. Petersburg State University, St. Petersburg, Russia; Department of Chemistry, Purdue University, West Lafayette, Indiana.
| |
Collapse
|
10
|
Sturzenegger F, Zosel F, Holmstrom ED, Buholzer KJ, Makarov DE, Nettels D, Schuler B. Transition path times of coupled folding and binding reveal the formation of an encounter complex. Nat Commun 2018; 9:4708. [PMID: 30413694 PMCID: PMC6226497 DOI: 10.1038/s41467-018-07043-x] [Citation(s) in RCA: 63] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2018] [Accepted: 10/15/2018] [Indexed: 11/09/2022] Open
Abstract
The association of biomolecules is the elementary event of communication in biology. Most mechanistic information of how the interactions between binding partners form or break is, however, hidden in the transition paths, the very short parts of the molecular trajectories from the encounter of the two molecules to the formation of a stable complex. Here we use single-molecule spectroscopy to measure the transition path times for the association of two intrinsically disordered proteins that form a folded dimer upon binding. The results reveal the formation of a metastable encounter complex that is electrostatically favored and transits to the final bound state within tens of microseconds. Such measurements thus open a new window into the microscopic events governing biomolecular interactions.
Collapse
Affiliation(s)
| | - Franziska Zosel
- Department of Biochemistry, University of Zurich, 8057, Zurich, Switzerland.,Novo Nordisk A/S, Novo Nordisk Park 1, 2760, Måløv, Denmark
| | - Erik D Holmstrom
- Department of Biochemistry, University of Zurich, 8057, Zurich, Switzerland
| | - Karin J Buholzer
- Department of Biochemistry, University of Zurich, 8057, Zurich, Switzerland
| | - Dmitrii E Makarov
- Department of Chemistry and Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, TX, 78712, USA
| | - Daniel Nettels
- Department of Biochemistry, University of Zurich, 8057, Zurich, Switzerland
| | - Benjamin Schuler
- Department of Biochemistry, University of Zurich, 8057, Zurich, Switzerland. .,Department of Physics, University of Zurich, 8057, Zurich, Switzerland.
| |
Collapse
|
11
|
Rodríguez Nassif A, de la Arada I, Arrondo JL, Pastrana-Rios B. 2D IR Correlation Spectroscopy in the Determination of Aggregation and Stability of KH Domain GXXG Loop Peptide in the Presence and Absence of Trifluoroacetate. Anal Chem 2017; 89:5765-5775. [PMID: 28459550 DOI: 10.1021/acs.analchem.6b04800] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Trifluoroacetate (TFA) is a strong anion byproduct of solid-phase peptide synthesis. Fourier transform infrared (FT-IR) spectroscopy can be used to ascertain the presence of this excipient in peptide samples for quality assessment. TFA absorbs as a strong sharp peak (1675 cm-1) within the amide I' band of the spectral region. A peptide sample and the TFA excipient can be studied simultaneously by FT-IR and 2D IR correlation spectroscopies. In addition, these techniques are able to determine the effect of TFA on the stability of the peptide. Herein, we describe the spectroscopic characterization of the GXXG loop peptide (GXXGlp), which is present in KH domain containing proteins. The sequence of the Homo sapiens Krr1 GXXGlp is evolutionarily conserved (165KRRQRLIGPKGSTLKALELLTNCY189) and has been associated with ssDNA interaction and ribosome biogenesis. Our goal was to determine the structural elements present in this peptide and evaluate whether TFA affects the stability of GXXGlp during thermal stress. We observed differences in the molecular behavior of the synthetic peptide in the presence and absence of TFA at various peptide concentrations. Finally, 2D IR correlation spectroscopy was used for the determination of the unfolding process, mechanism and extent of peptide aggregation, and the effect of TFA on the stability of the peptide. This spectroscopic method can be applied to the characterization of any synthetic peptide.
Collapse
Affiliation(s)
- Aslin Rodríguez Nassif
- Department of Chemistry, University of Puerto Rico , Mayagüez Campus, Mayagüez, Puerto Rico 00681-9019, United States
| | - Igor de la Arada
- Biofisika Institute and Biochemistry and Molecular Biology Department, CSIC and University of Basque Country , Bilbao, 48080, Spain
| | - José Luis Arrondo
- Biofisika Institute and Biochemistry and Molecular Biology Department, CSIC and University of Basque Country , Bilbao, 48080, Spain
| | - Belinda Pastrana-Rios
- Department of Chemistry, University of Puerto Rico , Mayagüez Campus, Mayagüez, Puerto Rico 00681-9019, United States.,Protein Research Center, University of Puerto Rico , Mayagüez Campus, Mayagüez, Puerto Rico 00681-9019, United States
| |
Collapse
|
12
|
Bessonov K, Vassall KA, Harauz G. Docking and molecular dynamics simulations of the Fyn-SH3 domain with free and phospholipid bilayer-associated 18.5-kDa myelin basic protein (MBP)-Insights into a noncanonical and fuzzy interaction. Proteins 2017; 85:1336-1350. [PMID: 28380689 DOI: 10.1002/prot.25295] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2016] [Revised: 03/03/2017] [Accepted: 03/27/2017] [Indexed: 01/06/2023]
Abstract
The molecular details of the association between the human Fyn-SH3 domain, and the fragment of 18.5-kDa myelin basic protein (MBP) spanning residues S38-S107 (denoted as xα2-peptide, murine sequence numbering), were studied in silico via docking and molecular dynamics over 50-ns trajectories. The results show that interaction between the two proteins is energetically favorable and heavily dependent on the MBP proline-rich region (P93-P98) in both aqueous and membrane environments. In aqueous conditions, the xα2-peptide/Fyn-SH3 complex adopts a "sandwich""-like structure. In the membrane context, the xα2-peptide interacts with the Fyn-SH3 domain via the proline-rich region and the β-sheets of Fyn-SH3, with the latter wrapping around the proline-rich region in a form of a clip. Moreover, the simulations corroborate prior experimental evidence of the importance of upstream segments beyond the canonical SH3-ligand. This study thus provides a more-detailed glimpse into the context-dependent interaction dynamics and importance of the β-sheets in Fyn-SH3 and proline-rich region of MBP. Proteins 2017; 85:1336-1350. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Kyrylo Bessonov
- Systems and Modeling Unit, Montefiore Institute, Université de Liège, Quartier Polytech 1, Allée de la Découverte 10, Liège, 4000, Belgium
| | - Kenrick A Vassall
- Department of Molecular and Cellular Biology, Biophysics Interdepartmental Group, and Collaborative Program in Neuroscience, University of Guelph, 50 Stone Road East, Guelph, Ontario, N1G 2W1, Canada
| | - George Harauz
- Department of Molecular and Cellular Biology, Biophysics Interdepartmental Group, and Collaborative Program in Neuroscience, University of Guelph, 50 Stone Road East, Guelph, Ontario, N1G 2W1, Canada
| |
Collapse
|
13
|
Olsson S, Noé F. Mechanistic Models of Chemical Exchange Induced Relaxation in Protein NMR. J Am Chem Soc 2016; 139:200-210. [DOI: 10.1021/jacs.6b09460] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Simon Olsson
- Computational Molecular
Biology,
FB Mathematik und Informatik, Freie Universität Berlin, Berlin 14195, Germany
| | - Frank Noé
- Computational Molecular
Biology,
FB Mathematik und Informatik, Freie Universität Berlin, Berlin 14195, Germany
| |
Collapse
|
14
|
Debiec KT, Cerutti DS, Baker LR, Gronenborn AM, Case DA, Chong LT. Further along the Road Less Traveled: AMBER ff15ipq, an Original Protein Force Field Built on a Self-Consistent Physical Model. J Chem Theory Comput 2016; 12:3926-47. [PMID: 27399642 PMCID: PMC4980686 DOI: 10.1021/acs.jctc.6b00567] [Citation(s) in RCA: 134] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present the AMBER
ff15ipq force field for proteins, the second-generation
force field developed using the Implicitly Polarized Q (IPolQ) scheme
for deriving implicitly polarized atomic charges in the presence of
explicit solvent. The ff15ipq force field is a complete rederivation
including more than 300 unique atomic charges, 900 unique torsion
terms, 60 new angle parameters, and new atomic radii for polar hydrogens.
The atomic charges were derived in the context of the SPC/Eb water model, which yields more-accurate rotational diffusion of
proteins and enables direct calculation of nuclear magnetic resonance
(NMR) relaxation parameters from molecular dynamics simulations. The
atomic radii improve the accuracy of modeling salt bridge interactions
relative to contemporary fixed-charge force fields, rectifying a limitation
of ff14ipq that resulted from its use of pair-specific Lennard-Jones
radii. In addition, ff15ipq reproduces penta-alanine J-coupling constants
exceptionally well, gives reasonable agreement with NMR relaxation
rates, and maintains the expected conformational propensities of structured
proteins/peptides, as well as disordered peptides—all on the
microsecond (μs) time scale, which is a critical regime for
drug design applications. These encouraging results demonstrate the
power and robustness of our automated methods for deriving new force
fields. All parameters described here and the mdgx program used to
fit them are included in the AmberTools16 distribution.
Collapse
Affiliation(s)
- Karl T Debiec
- Molecular Biophysics and Structural Biology Graduate Program, University of Pittsburgh and Carnegie Mellon University , Pittsburgh, Pennsylvania, United States.,Department of Structural Biology, University of Pittsburgh School of Medicine , Pittsburgh, Pennsylvania 15260, United States.,Department of Chemistry, University of Pittsburgh , Pittsburgh, Pennsylvania 15260, United States
| | - David S Cerutti
- Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| | - Lewis R Baker
- Department of Chemistry, University of Pittsburgh , Pittsburgh, Pennsylvania 15260, United States
| | - Angela M Gronenborn
- Department of Structural Biology, University of Pittsburgh School of Medicine , Pittsburgh, Pennsylvania 15260, United States
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University , New Brunswick, New Jersey 08854, United States
| | - Lillian T Chong
- Department of Chemistry, University of Pittsburgh , Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|