1
|
Zumpfe K, Berbon M, Habenstein B, Loquet A, Smith AA. Analytical Framework to Understand the Origins of Methyl Side-Chain Dynamics in Protein Assemblies. J Am Chem Soc 2024; 146:8164-8178. [PMID: 38476076 PMCID: PMC10979401 DOI: 10.1021/jacs.3c12620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Revised: 02/23/2024] [Accepted: 02/23/2024] [Indexed: 03/14/2024]
Abstract
Side-chain motions play an important role in understanding protein structure, dynamics, protein-protein, and protein-ligand interactions. However, our understanding of protein side-chain dynamics is currently limited by the lack of analytical tools. Here, we present a novel analytical framework employing experimental nuclear magnetic resonance (NMR) relaxation measurements at atomic resolution combined with molecular dynamics (MD) simulation to characterize with a high level of detail the methyl side-chain dynamics in insoluble protein assemblies, using amyloid fibrils formed by the prion HET-s. We use MD simulation to interpret experimental results, where rotameric hops, including methyl group rotation and χ1/χ2 rotations, cannot be completely described with a single correlation time but rather sample a broad distribution of correlation times, resulting from continuously changing local structure in the fibril. Backbone motion similarly samples a broad range of correlation times, from ∼100 ps to μs, although resulting from mostly different dynamic processes; nonetheless, we find that the backbone is not fully decoupled from the side-chain motion, where changes in side-chain dynamics influence backbone motion and vice versa. While the complexity of side-chain motion in protein assemblies makes it very challenging to obtain perfect agreement between experiment and simulation, our analytical framework improves the interpretation of experimental dynamics measurements for complex protein assemblies.
Collapse
Affiliation(s)
- Kai Zumpfe
- Institute
for Medical Physics and Biophysics, Leipzig
University, Härtelstraße
16-18, 04107 Leipzig, Germany
| | - Mélanie Berbon
- University
of Bordeaux, CNRS, Bordeaux INP, CBMN, UMR 5248, IECB, 33600 Pessac, France
| | - Birgit Habenstein
- University
of Bordeaux, CNRS, Bordeaux INP, CBMN, UMR 5248, IECB, 33600 Pessac, France
| | - Antoine Loquet
- University
of Bordeaux, CNRS, Bordeaux INP, CBMN, UMR 5248, IECB, 33600 Pessac, France
| | - Albert A. Smith
- Institute
for Medical Physics and Biophysics, Leipzig
University, Härtelstraße
16-18, 04107 Leipzig, Germany
| |
Collapse
|
2
|
Champion C, Lehner M, Smith AA, Ferrage F, Bolik-Coulon N, Riniker S. Unraveling motion in proteins by combining NMR relaxometry and molecular dynamics simulations: A case study on ubiquitin. J Chem Phys 2024; 160:104105. [PMID: 38465679 DOI: 10.1063/5.0188416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Accepted: 02/20/2024] [Indexed: 03/12/2024] Open
Abstract
Nuclear magnetic resonance (NMR) relaxation experiments shine light onto the dynamics of molecular systems in the picosecond to millisecond timescales. As these methods cannot provide an atomically resolved view of the motion of atoms, functional groups, or domains giving rise to such signals, relaxation techniques have been combined with molecular dynamics (MD) simulations to obtain mechanistic descriptions and gain insights into the functional role of side chain or domain motion. In this work, we present a comparison of five computational methods that permit the joint analysis of MD simulations and NMR relaxation experiments. We discuss their relative strengths and areas of applicability and demonstrate how they may be utilized to interpret the dynamics in MD simulations with the small protein ubiquitin as a test system. We focus on the aliphatic side chains given the rigidity of the backbone of this protein. We find encouraging agreement between experiment, Markov state models built in the χ1/χ2 rotamer space of isoleucine residues, explicit rotamer jump models, and a decomposition of the motion using ROMANCE. These methods allow us to ascribe the dynamics to specific rotamer jumps. Simulations with eight different combinations of force field and water model highlight how the different metrics may be employed to pinpoint force field deficiencies. Furthermore, the presented comparison offers a perspective on the utility of NMR relaxation to serve as validation data for the prediction of kinetics by state-of-the-art biomolecular force fields.
Collapse
Affiliation(s)
- Candide Champion
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Marc Lehner
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Albert A Smith
- Institute for Medical Physics and Biophysics, Leipzig University, Härtelstrasse 16-18, 04107 Leipzig, Germany
| | - Fabien Ferrage
- Laboratoire des Biomolécules, LBM, Département de Chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Nicolas Bolik-Coulon
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada
- Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Department of Biochemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Sereina Riniker
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
3
|
Agback T, Lesovoy D, Han X, Lomzov A, Sun R, Sandalova T, Orekhov VY, Achour A, Agback P. Combined NMR and molecular dynamics conformational filter identifies unambiguously dynamic ensembles of Dengue protease NS2B/NS3pro. Commun Biol 2023; 6:1193. [PMID: 38001280 PMCID: PMC10673835 DOI: 10.1038/s42003-023-05584-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Accepted: 11/14/2023] [Indexed: 11/26/2023] Open
Abstract
The dengue protease NS2B/NS3pro has been reported to adopt either an 'open' or a 'closed' conformation. We have developed a conformational filter that combines NMR with MD simulations to identify conformational ensembles that dominate in solution. Experimental values derived from relaxation parameters for the backbone and methyl side chains were compared with the corresponding back-calculated relaxation parameters of different conformational ensembles obtained from free MD simulations. Our results demonstrate a high prevalence for the 'closed' conformational ensemble while the 'open' conformation is absent, indicating that the latter conformation is most probably due to crystal contacts. Conversely, conformational ensembles in which the positioning of the co-factor NS2B results in a 'partially' open conformation, previously described in both MD simulations and X-ray studies, were identified by our conformational filter. Altogether, we believe that our approach allows for unambiguous identification of true conformational ensembles, an essential step for reliable drug discovery.
Collapse
Affiliation(s)
- Tatiana Agback
- Department of Molecular Sciences, Swedish University of Agricultural Sciences, PO Box 7015, SE-750 07, Uppsala, Sweden
| | - Dmitry Lesovoy
- Department of Structural Biology, Shemyakin-Ovchinnikov, Institute of Bioorganic Chemistry RAS, 117997, Moscow, Russia
- Swedish NMR Centre, University of Gothenburg, Box 465, 40530, Gothenburg, Sweden
| | - Xiao Han
- Science for Life Laboratory, Department of Medicine, Karolinska Institute, and Division of Infectious Diseases, Karolinska University Hospital, SE-171 76, Stockholm, Sweden
| | - Alexander Lomzov
- Laboratory of Structural Biology, Institute of Chemical Biology and Fundamental Medicine SB RAS, 630090, Novosibirsk, Russia
| | - Renhua Sun
- Science for Life Laboratory, Department of Medicine, Karolinska Institute, and Division of Infectious Diseases, Karolinska University Hospital, SE-171 76, Stockholm, Sweden
| | - Tatyana Sandalova
- Science for Life Laboratory, Department of Medicine, Karolinska Institute, and Division of Infectious Diseases, Karolinska University Hospital, SE-171 76, Stockholm, Sweden
| | - Vladislav Yu Orekhov
- Swedish NMR Centre, University of Gothenburg, Box 465, 40530, Gothenburg, Sweden
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 465, 40530, Gothenburg, Sweden
| | - Adnane Achour
- Science for Life Laboratory, Department of Medicine, Karolinska Institute, and Division of Infectious Diseases, Karolinska University Hospital, SE-171 76, Stockholm, Sweden.
| | - Peter Agback
- Department of Molecular Sciences, Swedish University of Agricultural Sciences, PO Box 7015, SE-750 07, Uppsala, Sweden.
| |
Collapse
|
4
|
Kümmerer F, Orioli S, Lindorff-Larsen K. Fitting Force Field Parameters to NMR Relaxation Data. J Chem Theory Comput 2023. [PMID: 37276045 DOI: 10.1021/acs.jctc.3c00174] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We present an approach to optimize force field parameters using time-dependent data from NMR relaxation experiments. To do so, we scan parameters in the dihedral angle potential energy terms describing the rotation of the methyl groups in proteins and compare NMR relaxation rates calculated from molecular dynamics simulations with the modified force fields to deuterium relaxation measurements of T4 lysozyme. We find that a small modification of Cγ methyl groups improves the agreement with experiments both for the protein used to optimize the force field and when validating using simulations of CI2 and ubiquitin. We also show that these improvements enable a more effective a posteriori reweighting of the MD trajectories. The resulting force field thus enables more direct comparison between simulations and side-chain NMR relaxation data and makes it possible to construct ensembles that better represent the dynamics of proteins in solution.
Collapse
Affiliation(s)
- Felix Kümmerer
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Simone Orioli
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
- Structural Biophysics, Niels Bohr Institute, Faculty of Science, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
5
|
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
|
6
|
How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation. Curr Opin Struct Biol 2022; 77:102459. [PMID: 36148743 DOI: 10.1016/j.sbi.2022.102459] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2022] [Revised: 07/29/2022] [Accepted: 08/05/2022] [Indexed: 12/14/2022]
Abstract
Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from 15N relaxation and side-chain dynamics from 13C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.
Collapse
|
7
|
Bolik-Coulon N, Languin-Cattoën O, Carnevale D, Zachrdla M, Laage D, Sterpone F, Stirnemann G, Ferrage F. Explicit Models of Motion to Understand Protein Side-Chain Dynamics. PHYSICAL REVIEW LETTERS 2022; 129:203001. [PMID: 36462011 DOI: 10.1103/physrevlett.129.203001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Accepted: 09/04/2022] [Indexed: 06/17/2023]
Abstract
Nuclear magnetic relaxation is widely used to probe protein dynamics. For decades, most analyses of relaxation in proteins have relied successfully on the model-free approach, forgoing mechanistic descriptions of motion. Model-free types of correlation functions cannot describe a large carbon-13 relaxation dataset in protein side chains. Here, we use molecular dynamics simulations to design explicit models of motion and solve Fokker-Planck diffusion equations. These models of motion provide better agreement with relaxation data, mechanistic insight, and a direct link to configuration entropy.
Collapse
Affiliation(s)
- Nicolas Bolik-Coulon
- Laboratoire des Biomolécules, LBM, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
| | - Olivier Languin-Cattoën
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologique Physico-Chimique, Université Paris Cité, PSL University, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Diego Carnevale
- Laboratoire des Biomolécules, LBM, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
| | - Milan Zachrdla
- Laboratoire des Biomolécules, LBM, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologique Physico-Chimique, Université Paris Cité, PSL University, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Damien Laage
- PASTEUR, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
| | - Fabio Sterpone
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologique Physico-Chimique, Université Paris Cité, PSL University, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Guillaume Stirnemann
- CNRS Laboratoire de Biochimie Théorique, Institut de Biologique Physico-Chimique, Université Paris Cité, PSL University, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabien Ferrage
- Laboratoire des Biomolécules, LBM, Département de chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
8
|
Hussain A, Paukovich N, Henen MA, Vögeli B. Advances in the exact nuclear Overhauser effect 2018-2022. Methods 2022; 206:87-98. [PMID: 35985641 PMCID: PMC9596134 DOI: 10.1016/j.ymeth.2022.08.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2022] [Revised: 08/05/2022] [Accepted: 08/12/2022] [Indexed: 11/26/2022] Open
Abstract
The introduction of the exact nuclear Overhauser enhancement (eNOE) methodology to solution-state nuclear magnetic resonance (NMR) spectroscopy results in tighter distance restraints from NOEs than in convention analysis. These improved restraints allow for higher resolution in structure calculation and even the disentanglement of different conformations of macromolecules. While initial work primarily focused on technical development of the eNOE, structural studies aimed at the elucidation of spatial sampling in proteins and nucleic acids were published in parallel prior to 2018. The period of 2018-2022 saw a continued series of technical innovation, but also major applications addressing biological questions. Here, we review both aspects, covering topics from the implementation of non-uniform sampling of NOESY buildups, novel pulse sequences, adaption of the eNOE to solid-state NMR, advances in eNOE data analysis, and innovations in structural ensemble calculation, to applications to protein, RNA, and DNA structure elucidation.
Collapse
Affiliation(s)
- Alya Hussain
- Department of Biochemistry & Molecular Genetics, School of Medicine, University of Colorado, 12801 E. 17(th) Avenue, Aurora, CO 80045, USA
| | - Natasia Paukovich
- Department of Biochemistry & Molecular Genetics, School of Medicine, University of Colorado, 12801 E. 17(th) Avenue, Aurora, CO 80045, USA
| | - Morkos A Henen
- Department of Biochemistry & Molecular Genetics, School of Medicine, University of Colorado, 12801 E. 17(th) Avenue, Aurora, CO 80045, USA; Department of Pharmaceutical Organic Chemistry, Faculty of Pharmacy, Mansoura University, Mansoura 35516, Egypt
| | - Beat Vögeli
- Department of Biochemistry & Molecular Genetics, School of Medicine, University of Colorado, 12801 E. 17(th) Avenue, Aurora, CO 80045, USA.
| |
Collapse
|
9
|
Bolik-Coulon N, Ferrage F. Explicit models of motions to analyze NMR relaxation data in proteins. J Chem Phys 2022; 157:125102. [DOI: 10.1063/5.0095910] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Nuclear Magnetic Resonance (NMR) is a tool of choice to characterize molecular motions. In biological macromolecules, pico- to nano-second motions, in particular, can be probed by nuclear spin relaxation rates which depend on the time fluctuations of the orientations of spin interaction frames. For the past 40 years, relaxation rates have been successfully analyzed using the Model Free (MF) approach which makes no assumption on the nature of motions and reports on the effective amplitude and time-scale of the motions. However, obtaining a mechanistic picture of motions from this type of analysis is difficult at best, unless complemented with molecular dynamics (MD) simulations. In spite of their limited accuracy, such simulations can be used to obtain the information necessary to build explicit models of motions designed to analyze NMR relaxation data. Here, we present how to build such models, suited in particular to describe motions of methyl-bearing protein side-chains and compare them with the MF approach. We show on synthetic data that explicit models of motions are more robust in the presence of rotamer jumps which dominate the relaxation in methyl groups of protein side-chains. We expect this work to motivate the use of explicit models of motion to analyze MD and NMR data.
Collapse
Affiliation(s)
| | - Fabien Ferrage
- Departement de chimie, Ecole Normale Superieure Departement de Chimie, France
| |
Collapse
|
10
|
Love O, Pacheco Lima MC, Clark C, Cornillie S, Roalstad S, Cheatham TE. Evaluating the accuracy of the AMBER protein force fields in modeling dihydrofolate reductase structures: misbalance in the conformational arrangements of the flexible loop domains. J Biomol Struct Dyn 2022:1-15. [PMID: 35838167 PMCID: PMC9840716 DOI: 10.1080/07391102.2022.2098823] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Protein flexible loop regions were once thought to be simple linkers between other more functional secondary structural elements. However, as it becomes clearer that these loop domains are critical players in a plethora of biological processes, accurate conformational sampling of 3D loop structures is vital to the advancement of drug design techniques and the overall growth of knowledge surrounding molecular systems. While experimental techniques provide a wealth of structural information, the resolution of flexible loop domains is sometimes low or entirely absent due to their complex and dynamic nature. This highlights an opportunity for de novo structure prediction using in silico methods with molecular dynamics (MDs). This study evaluates some of the AMBER protein force field's (ffs) ability to accurately model dihydrofolate reductase (DHFR) conformations, a protein complex characterized by specific arrangements and interactions of multiple flexible loops whose conformations are determined by the presence or absence of bound ligands and cofactors. Although the AMBER ffs, including ff19SB, studied well model most protein structures with rich secondary structure, results obtained here suggest the inability to significantly sample the expected DHFR loop-loop conformations - of the six distinct protein-ligand systems simulated, a majority lacked consistent stabilization of experimentally derived metrics definitive the three enzyme conformations. Although under-sampling and the chosen ff parameter combinations could be the cause, given past successes with these MD approaches for many protein systems, this suggests a potential misbalance in available ff parameters required to accurately predict the structure of multiple flexible loop regions present in proteins.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Olivia Love
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT, USA
| | | | - Casey Clark
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT, USA
| | - Sean Cornillie
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT, USA
| | - Shelly Roalstad
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT, USA
| | - Thomas E. Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, Salt Lake City, UT, USA
| |
Collapse
|
11
|
Yang D, Gronenborn AM, Chong LT. Development and Validation of Fluorinated, Aromatic Amino Acid Parameters for Use with the AMBER ff15ipq Protein Force Field. J Phys Chem A 2022; 126:2286-2297. [PMID: 35352936 PMCID: PMC9014858 DOI: 10.1021/acs.jpca.2c00255] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 03/15/2022] [Indexed: 12/27/2022]
Abstract
We developed force field parameters for fluorinated, aromatic amino acids enabling molecular dynamics (MD) simulations of fluorinated proteins. These parameters are tailored to the AMBER ff15ipq protein force field and enable the modeling of 4, 5, 6, and 7F-tryptophan, 3F- and 3,5F-tyrosine, and 4F- or 4-CF3-phenylalanine. The parameters include 181 unique atomic charges derived using the implicitly polarized charge (IPolQ) scheme in the presence of SPC/Eb explicit water molecules and 9 unique bond, angle, or torsion terms. Our simulations of benchmark peptides and proteins maintain expected conformational propensities on the μs time scale. In addition, we have developed an open-source Python program to calculate fluorine relaxation rates from MD simulations. The extracted relaxation rates from protein simulations are in good agreement with experimental values determined by 19F NMR. Collectively, our results illustrate the power and robustness of the IPolQ lineage of force fields for modeling the structure and dynamics of fluorine-containing proteins at the atomic level.
Collapse
Affiliation(s)
- Darian
T. Yang
- Molecular
Biophysics and Structural Biology Graduate Program, University of Pittsburgh and Carnegie Mellon University, Pittsburgh, Pennsylvania 15260, 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
| | - Angela M. Gronenborn
- 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
| | - Lillian T. Chong
- Department
of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
12
|
Hoffmann F, Mulder FAA, Schäfer LV. How Much Entropy Is Contained in NMR Relaxation Parameters? J Phys Chem B 2021; 126:54-68. [PMID: 34936366 DOI: 10.1021/acs.jpcb.1c07786] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Solution-state NMR relaxation experiments are the cornerstone to study internal protein dynamics at an atomic resolution on time scales that are faster than the overall rotational tumbling time τR. Since the motions described by NMR relaxation parameters are connected to thermodynamic quantities like conformational entropies, the question arises how much of the total entropy is contained within this tumbling time. Using all-atom molecular dynamics simulations of the T4 lysozyme, we found that entropy buildup is rather fast for the backbone, such that the majority of the entropy is indeed contained in the short-time dynamics. In contrast, the contribution of the slow dynamics of side chains on time scales beyond τR on the side-chain conformational entropy is significant and should be taken into account for the extraction of accurate thermodynamic properties.
Collapse
Affiliation(s)
- Falk Hoffmann
- Center for Theoretical Chemistry, Ruhr University Bochum, D-44 780 Bochum, Germany
| | - Frans A A Mulder
- Interdisciplinary Nanoscience Center and Department of Chemistry, University of Aarhus, DK-8000 Aarhus, Denmark
| | - Lars V Schäfer
- Center for Theoretical Chemistry, Ruhr University Bochum, D-44 780 Bochum, Germany
| |
Collapse
|
13
|
Abstract
Relaxation in nuclear magnetic resonance is a powerful method for obtaining spatially resolved, timescale-specific dynamics information about molecular systems. However, dynamics in biomolecular systems are generally too complex to be fully characterized based on NMR data alone. This is a familiar problem, addressed by the Lipari-Szabo model-free analysis, a method that captures the full information content of NMR relaxation data in case all internal motion of a molecule in solution is sufficiently fast. We investigate model-free analysis, as well as several other approaches, and find that model-free, spectral density mapping, LeMaster's approach, and our detector analysis form a class of analysis methods, for which behavior of the fitted parameters has a well-defined relationship to the distribution of correlation times of motion, independent of the specific form of that distribution. In a sense, they are all "model-free." Of these methods, only detectors are generally applicable to solid-state NMR relaxation data. We further discuss how detectors may be used for comparison of experimental data to data extracted from molecular dynamics simulation, and how simulation may be used to extract details of the dynamics that are not accessible via NMR, where detector analysis can be used to connect those details to experiments. We expect that combined methodology can eventually provide enough insight into complex dynamics to provide highly accurate models of motion, thus lending deeper insight into the nature of biomolecular dynamics.
Collapse
Affiliation(s)
- Kai Zumpfe
- Institute for Medical Physics and Biophysics, Medical Faculty, Leipzig University, Leipzig, Germany
| | - Albert A Smith
- Institute for Medical Physics and Biophysics, Medical Faculty, Leipzig University, Leipzig, Germany
| |
Collapse
|
14
|
Hicks A, MacAinsh M, Zhou HX. Removing Thermostat Distortions of Protein Dynamics in Constant-Temperature Molecular Dynamics Simulations. J Chem Theory Comput 2021; 17:5920-5932. [PMID: 34464112 DOI: 10.1021/acs.jctc.1c00448] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Molecular dynamics simulations are widely used to determine equilibrium and dynamic properties of proteins. Nearly all simulations, currently, are carried out at constant temperature, with a Langevin thermostat among the most widely used. Thermostats distort protein dynamics, but whether or how such distortions can be corrected has long been an open question. Here, we show that constant-temperature simulations with a Langevin thermostat dilate protein dynamics and present a correction scheme to remove the dynamic distortions. Specifically, ns-scale time constants for overall rotation are dilated significantly but sub-ns time constants for internal motions are dilated modestly, while all motional amplitudes are unaffected. The correction scheme involves contraction of the time constants, with the contraction factor a linear function of the time constant to be corrected. The corrected dynamics of eight proteins are validated by NMR data for rotational diffusion and for backbone amide and side-chain methyl relaxation. The present work demonstrates that even for complex systems like proteins with dynamics spanning multiple timescales, one can predict how thermostats distort protein dynamics and remove such distortions. The correction scheme will have wide applications, facilitating force-field parameterization and propelling simulations to be on par with NMR and other experimental techniques in determining dynamic properties of proteins.
Collapse
|
15
|
Kümmerer F, Orioli S, Harding-Larsen D, Hoffmann F, Gavrilov Y, Teilum K, Lindorff-Larsen K. Fitting Side-Chain NMR Relaxation Data Using Molecular Simulations. J Chem Theory Comput 2021; 17:5262-5275. [PMID: 34291646 DOI: 10.1021/acs.jctc.0c01338] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Proteins display a wealth of dynamical motions that can be probed using both experiments and simulations. We present an approach to integrate side-chain NMR relaxation measurements with molecular dynamics simulations to study the structure and dynamics of these motions. The approach, which we term ABSURDer (average block selection using relaxation data with entropy restraints), can be used to find a set of trajectories that are in agreement with relaxation measurements. We apply the method to deuterium relaxation measurements in T4 lysozyme and show how it can be used to integrate the accuracy of the NMR measurements with the molecular models of protein dynamics afforded by the simulations. We show how fitting of dynamic quantities leads to improved agreement with static properties and highlight areas needed for further improvements of the approach.
Collapse
Affiliation(s)
- Felix Kümmerer
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Simone Orioli
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark.,Structural Biophysics, Niels Bohr Institute, Faculty of Science, University of Copenhagen, DK-2100 Copenhagen, Denmark
| | - David Harding-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Falk Hoffmann
- Theoretical Chemistry, Ruhr University Bochum, D-44780 Bochum, Germany
| | - Yulian Gavrilov
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Kaare Teilum
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
16
|
Karunanithy G, Shukla VK, Hansen DF. Methodological advancements for characterising protein side chains by NMR spectroscopy. Curr Opin Struct Biol 2021; 70:61-69. [PMID: 33989947 DOI: 10.1016/j.sbi.2021.04.002] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 03/31/2021] [Accepted: 04/01/2021] [Indexed: 11/18/2022]
Abstract
The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.
Collapse
Affiliation(s)
- Gogulan Karunanithy
- Department of Structural and Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, United Kingdom
| | - Vaibhav Kumar Shukla
- Department of Structural and Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, United Kingdom
| | - D Flemming Hansen
- Department of Structural and Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, United Kingdom.
| |
Collapse
|
17
|
Grohe K, Patel S, Hebrank C, Medina S, Klein A, Rovó P, Vasa SK, Singh H, Vögeli B, Schäfer LV, Linser R. Protein Motional Details Revealed by Complementary Structural Biology Techniques. Structure 2020; 28:1024-1034.e3. [PMID: 32579946 DOI: 10.1016/j.str.2020.06.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Revised: 05/05/2020] [Accepted: 06/03/2020] [Indexed: 01/16/2023]
Abstract
Proteins depend on defined molecular plasticity for their functionality. How to comprehensively capture dynamics correctly is of ubiquitous biological importance. Approaches commonly used to probe protein dynamics include model-free elucidation of site-specific motion by NMR relaxation, molecular dynamics (MD)-based approaches, and capturing the substates within a dynamic ensemble by recent eNOE-based multiple-structure approaches. Even though MD is sometimes combined with ensemble-averaged NMR restraints, these approaches have largely been developed and used individually. Owing to the different underlying concepts and practical requirements, it has remained unclear how they compare, and how they cross-validate and complement each other. Here, we extract and compare the differential information contents of MD simulations, NMR relaxation measurements, and eNOE-based multi-state structures for the SH3 domain of chicken α-spectrin. The data show that a validated, consistent, and detailed picture is feasible both for timescales and actual conformational states sampled in the dynamic ensemble. This includes the biologically important side-chain plasticity, for which experimentally cross-validated assessment is a significant challenge.
Collapse
Affiliation(s)
- Kristof Grohe
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany; Faculty of Chemistry and Chemical Biology, Technical University Dortmund, 44227 Dortmund, Germany
| | - Snehal Patel
- Theoretical Chemistry, Ruhr University Bochum, 44801 Bochum, Germany
| | - Cornelia Hebrank
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany
| | - Sara Medina
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany; Faculty of Chemistry and Chemical Biology, Technical University Dortmund, 44227 Dortmund, Germany
| | - Alexander Klein
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany; Faculty of Chemistry and Chemical Biology, Technical University Dortmund, 44227 Dortmund, Germany
| | - Petra Rovó
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany
| | - Suresh K Vasa
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany; Faculty of Chemistry and Chemical Biology, Technical University Dortmund, 44227 Dortmund, Germany
| | - Himanshu Singh
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany; Faculty of Chemistry and Chemical Biology, Technical University Dortmund, 44227 Dortmund, Germany
| | - Beat Vögeli
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver, Aurora, CO 80045, USA
| | - Lars V Schäfer
- Theoretical Chemistry, Ruhr University Bochum, 44801 Bochum, Germany.
| | - Rasmus Linser
- Faculty for Chemistry and Pharmacy, Ludwig-Maximilians-University Munich, 81377 Munich, Germany; Faculty of Chemistry and Chemical Biology, Technical University Dortmund, 44227 Dortmund, Germany.
| |
Collapse
|
18
|
Marcellini M, Nguyen MH, Martin M, Hologne M, Walker O. Accurate Prediction of Protein NMR Spin Relaxation by Means of Polarizable Force Fields. Application to Strongly Anisotropic Rotational Diffusion. J Phys Chem B 2020; 124:5103-5112. [DOI: 10.1021/acs.jpcb.0c01922] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Moreno Marcellini
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Minh-Ha Nguyen
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Marie Martin
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Maggy Hologne
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Olivier Walker
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| |
Collapse
|
19
|
Anderson JS, Hernández G, LeMaster DM. 13C NMR Relaxation Analysis of Protein GB3 for the Assessment of Side Chain Dynamics Predictions by Current AMBER and CHARMM Force Fields. J Chem Theory Comput 2020; 16:2896-2913. [PMID: 32268062 DOI: 10.1021/acs.jctc.0c00050] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Molecular simulations with seven current AMBER- and CHARMM-based force fields yield markedly differing internal bond vector autocorrelation function predictions for many of the 223 methine and methylene H-C bonds of the 56-residue protein GB3. To enable quantification of accuracy, 13C R1, R2, and heteronuclear NOE relaxation rates have been determined for the methine and stereochemically assigned methylene Cα and Cβ positions. With only three experimental relaxation values for each bond vector, central to this analysis is the accuracy with which MD-derived autocorrelation curves can be represented by a 3-parameter equation which, in turn, maps onto the NMR relaxation values. In contrast to the more widely used extended Lipari-Szabo order parameter representation, 95% of these MD-derived internal autocorrelation curves for GB3 can be fitted to within 1.0% rmsd over the time frame from 30 ps to 4 ns by a biexponential Larmor frequency-selective representation (LF-S2). Applying the LF-S2 representation to the experimental relaxation rates and uncertainties serves to determine the boundary range for the autocorrelation function of each bond vector consistent with the experimental data. Not surprisingly, all seven force fields predict the autocorrelation functions for the more motionally restricted 1Hα-13Cα and 1Hβ-13Cβ bond vectors with reasonable accuracy. However, for the 1Hβ-13Cβ bond vectors exhibiting aggregate order parameter S2 values less than 0.85, only 1% of the MD-derived predictions lie with 1 σ of the experimentally determined autocorrelation functions and only 7% within 2 σ. On the other hand, substantial residue type-specific improvements in predictive performance were observed among the recent AMBER force fields. This analysis indicates considerable potential for the use of 13C relaxation measurements in guiding the optimization of the side chain dynamics characteristics of protein molecular simulations.
Collapse
Affiliation(s)
- Janet S Anderson
- Department of Chemistry, Union College, Schenectady, New York 12308, United States
| | - Griselda Hernández
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, New York 12201, United States
| | - David M LeMaster
- Wadsworth Center, New York State Department of Health, Empire State Plaza, Albany, New York 12201, United States
| |
Collapse
|
20
|
Rashid S, Lee BL, Wajda B, Spyracopoulos L. Nucleotide Binding and Active Site Gate Dynamics for the Hsp90 Chaperone ATPase Domain from Benchtop and High Field 19F NMR Spectroscopy. J Phys Chem B 2020; 124:2984-2993. [PMID: 32212608 DOI: 10.1021/acs.jpcb.0c00626] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein turnover in cells is regulated by the ATP dependent activity of the Hsp90 chaperone. In concert with accessory proteins, ATP hydrolysis drives the obligate Hsp90 dimer through a cycle between open and closed states that is critical for assisting the folding and stability of hundreds of proteins. Cycling is initiated by ATP binding to the ATPase domain, with the chaperone and the active site gates in the dimer in open states. The chaperone then adopts a short-lived, ATP bound closed state with a closed active site gate. The structural and dynamic changes induced in the ATPase domain and active site gate upon nucleotide binding, and their impact on dimer closing are not well understood. We site-specifically 19F-labeled the ATPase domain at the active site gate to enable benchtop and high field 19F NMR spectroscopic studies. Combined with MD simulations, this allowed accurate characterization of pico- to nanosecond time scale motions of the active site gate, as well as slower micro- to millisecond time scale processes resulting from nucleotide binding. ATP binding induces increased flexibility at one of the hinges of the active site gate, a necessary prelude to release of the second hinge and eventual gate closure in the intact chaperone.
Collapse
Affiliation(s)
- Suad Rashid
- Department of Biochemistry, University of Alberta, Edmonton, AB T6G 2H7, Canada
| | - Brian L Lee
- Department of Biochemistry, University of Alberta, Edmonton, AB T6G 2H7, Canada
| | - Benjamin Wajda
- Department of Biochemistry, University of Alberta, Edmonton, AB T6G 2H7, Canada
| | - Leo Spyracopoulos
- Department of Biochemistry, University of Alberta, Edmonton, AB T6G 2H7, Canada
| |
Collapse
|
21
|
Hoffmann F, Mulder FAA, Schäfer LV. Predicting NMR relaxation of proteins from molecular dynamics simulations with accurate methyl rotation barriers. J Chem Phys 2020; 152:084102. [DOI: 10.1063/1.5135379] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Falk Hoffmann
- Theoretical Chemistry, Ruhr University Bochum, D-44780 Bochum, Germany
| | - Frans A. A. Mulder
- Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry, Aarhus University, DK-8000 Aarhus, Denmark
| | - Lars V. Schäfer
- Theoretical Chemistry, Ruhr University Bochum, D-44780 Bochum, Germany
| |
Collapse
|
22
|
Emani PS, Yimer YY, Davidowski SK, Gebhart RN, Ferreira HE, Kuprov I, Pfaendtner J, Drobny GP. Combining Molecular and Spin Dynamics Simulations with Solid-State NMR: A Case Study of Amphiphilic Lysine-Leucine Repeat Peptide Aggregates. J Phys Chem B 2019; 123:10915-10929. [PMID: 31769684 DOI: 10.1021/acs.jpcb.9b09245] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Interpreting dynamics in solid-state molecular systems requires characterization of the potentially heterogeneous environmental contexts of molecules. In particular, the analysis of solid-state nuclear magnetic resonance (ssNMR) data to elucidate molecular dynamics (MD) involves modeling the restriction to overall tumbling by neighbors, as well as the concentrations of water and buffer. In this exploration of the factors that influence motion, we utilize atomistic MD trajectories of peptide aggregates with varying hydration to mimic an amorphous solid-state environment and predict ssNMR relaxation rates. We also account for spin diffusion in multiply spin-labeled (up to 19 nuclei) residues, with several models of dipolar-coupling networks. The framework serves as a general approach to determine essential spin couplings affecting relaxation, benchmark MD force fields, and reveal the hydration dependence of dynamics in a crowded environment. We demonstrate the methodology on a previously characterized amphiphilic 14-residue lysine-leucine repeat peptide, LKα14 (Ac-LKKLLKLLKKLLKL-c), which has an α-helical secondary structure and putatively forms leucine-burying tetramers in the solid state. We measure the R1 relaxation rates of uniformly 13C-labeled and site-specific 2H-labeled leucines in the hydrophobic core of LKα14 at multiple hydration levels. Studies of 9 and 18 tetramer bundles reveal the following: (a) for the incoherent component of 13C relaxation, the nearest-neighbor spin interactions dominate, while the 1H-1H interactions have minimal impact; (b) the AMBER ff14SB dihedral barriers for the leucine Cγ-Cδ bond ("methyl rotation barriers") must be lowered by a factor of 0.7 to better match the 2H data; (c) proton-driven spin diffusion explains some of the discrepancy between experimental and simulated rates for the Cβ and Cα nuclei; and (d) 13C relaxation rates are mostly underestimated in the MD simulations at all hydrations, and the discrepancies identify likely motions missing in the 50 ns MD trajectories.
Collapse
Affiliation(s)
- Prashant S Emani
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Yeneneh Y Yimer
- Department of Chemical Engineering , University of Washington , 105 Benson Hall, Box 351750 , Seattle , Washington 98195-1750 , United States
| | - Stephen K Davidowski
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Rachel N Gebhart
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Helen E Ferreira
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| | - Ilya Kuprov
- Department of Chemistry , University of Southampton , Highfield, Southampton SO17 1BJ , U.K
| | - Jim Pfaendtner
- Department of Chemical Engineering , University of Washington , 105 Benson Hall, Box 351750 , Seattle , Washington 98195-1750 , United States
| | - Gary P Drobny
- Department of Chemistry , University of Washington , Box 351700 , Seattle , Washington 98195-1700 , United States
| |
Collapse
|
23
|
Hata H, Nishiyama M, Kitao A. Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics. Biochim Biophys Acta Gen Subj 2019; 1864:129395. [PMID: 31302180 DOI: 10.1016/j.bbagen.2019.07.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Revised: 07/03/2019] [Accepted: 07/08/2019] [Indexed: 12/17/2022]
Abstract
BACKGROUND Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations. SCOPE OF REVIEW First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure. MAJOR CONCLUSIONS Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening. GENERAL SIGNIFICANCE MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.
Collapse
Affiliation(s)
- Hiroaki Hata
- School of Life Science and Technology, Tokyo Institute of Technology, Ookayama, 2-12-1 Meguro-ku, Tokyo 152-8550, Japan
| | - Masayoshi Nishiyama
- Department of Physics, Kindai University, 3-4-1 Kowakae, Higashiosaka, Osaka 577-8502, Japan
| | - Akio Kitao
- School of Life Science and Technology, Tokyo Institute of Technology, Ookayama, 2-12-1 Meguro-ku, Tokyo 152-8550, Japan.
| |
Collapse
|
24
|
Polimeno A, Zerbetto M, Abergel D. Stochastic modeling of macromolecules in solution. I. Relaxation processes. J Chem Phys 2019; 150:184107. [PMID: 31091939 DOI: 10.1063/1.5077065] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
A framework for the stochastic description of relaxation processes in flexible macromolecules, including dissipative effects, is introduced from an atomistic point of view. Projection-operator techniques are employed to obtain multidimensional Fokker-Planck operators governing the relaxation of internal coordinates and global degrees of freedom and depending upon parameters fully recoverable from classic force fields (energetics) and continuum models (friction tensors). A hierarchy of approaches of different complexity is proposed in this unified context, aimed primarily at the interpretation of magnetic resonance relaxation experiments. In particular, a model based on a harmonic internal Hamiltonian is discussed as a test case.
Collapse
Affiliation(s)
- Antonino Polimeno
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, I-35131 Padova, Italy
| | - Mirco Zerbetto
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, I-35131 Padova, Italy
| | - Daniel Abergel
- Laboratoire des Biomolécules, LBM, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
25
|
Rashid S, Lee BL, Wajda B, Spyracopoulos L. Side-Chain Dynamics of the Trifluoroacetone Cysteine Derivative Characterized by 19F NMR Relaxation and Molecular Dynamics Simulations. J Phys Chem B 2019; 123:3665-3671. [PMID: 30973726 DOI: 10.1021/acs.jpcb.9b01741] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
19F NMR spectroscopy is a powerful tool for the study of the structures, dynamics, and interactions of proteins bearing cysteine residues chemically modified with a trifluoroacetone group (CYF residue). 19F NMR relaxation rates for the fluoromethyl group of CYF residues are sensitive to overall rotational tumbling of proteins, fast rotation about the CF3 methyl axis, and the internal motion of the CYF side-chain. To develop a quantitative understanding of these various motional contributions, we used the model-free approach to extend expressions for 19F- T2 NMR relaxation to include side-chain motions for the CYF residue. We complemented the NMR studies with atomic views of methyl rotation and side-chain motions using molecular dynamics simulations. This combined methodology allows for quantitative separation of the contributions of fast pico- to nanosecond dynamics from micro- to millisecond exchange processes to the 19F line width and highlights the utility of the CYF residue as a sensitive reporter of side-chain environment and dynamics in proteins.
Collapse
Affiliation(s)
- Suad Rashid
- Department of Biochemistry , University of Alberta , Edmonton , Alberta T6G 2H7 , Canada
| | - Brian L Lee
- Department of Biochemistry , University of Alberta , Edmonton , Alberta T6G 2H7 , Canada
| | - Benjamin Wajda
- Department of Biochemistry , University of Alberta , Edmonton , Alberta T6G 2H7 , Canada
| | - Leo Spyracopoulos
- Department of Biochemistry , University of Alberta , Edmonton , Alberta T6G 2H7 , Canada
| |
Collapse
|
26
|
Tchaicheeyan O, Mendelman N, Zerbetto M, Meirovitch E. Local Ordering at Mobile Sites in Proteins: Combining Perspectives from NMR Relaxation and Molecular Dynamics. J Phys Chem B 2019; 123:2745-2755. [DOI: 10.1021/acs.jpcb.8b10801] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Oren Tchaicheeyan
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel
| | - Netanel Mendelman
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel
| | - Mirco Zerbetto
- Department of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Eva Meirovitch
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel
| |
Collapse
|
27
|
Abstract
NMR-based studies of protein dynamics and molecular simulations have a synergistic relationship. Molecular simulations, in combination with interpretative theoretical models, leverage the dynamical information obtained from NMR. They provide the concrete physical schema underlying the quantities measured by NMR, and help extend the range of applications beyond the strictly dynamic properties. NMR data in turn provide concrete data to test and improve the potential functions used for simulation of dynamics of proteins. The concept of time correlation functions is central to the understanding of many dynamical processes. Their evaluation through atomistic simulations is discussed, with application to different dynamical quantities measured by NMR. While advances in computers have made such atomistic simulations almost routine, the companion use of simple interpretive models is stressed, to provide not just numbers but physical insight.
Collapse
Affiliation(s)
- Kim A Sharp
- Department of Biochemistry and Biophysics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA, United States.
| |
Collapse
|
28
|
Cousin SF, Kadeřávek P, Bolik-Coulon N, Gu Y, Charlier C, Carlier L, Bruschweiler-Li L, Marquardsen T, Tyburn JM, Brüschweiler R, Ferrage F. Time-Resolved Protein Side-Chain Motions Unraveled by High-Resolution Relaxometry and Molecular Dynamics Simulations. J Am Chem Soc 2018; 140:13456-13465. [DOI: 10.1021/jacs.8b09107] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Samuel F. Cousin
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Pavel Kadeřávek
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Nicolas Bolik-Coulon
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Yina Gu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Cyril Charlier
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Ludovic Carlier
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Lei Bruschweiler-Li
- Campus Chemical Instrument Center, The Ohio State University, Columbus, Ohio 43210, United States
| | | | - Jean-Max Tyburn
- Bruker BioSpin, 34 rue de l’Industrie BP 10002, 67166 Wissembourg Cedex, France
| | - Rafael Brüschweiler
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Campus Chemical Instrument Center, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Biological Chemistry and Pharmacology, The Ohio State University, Columbus, Ohio 43210, United States
| | - Fabien Ferrage
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
29
|
Ollila OHS, Heikkinen HA, Iwaï H. Rotational Dynamics of Proteins from Spin Relaxation Times and Molecular Dynamics Simulations. J Phys Chem B 2018; 122:6559-6569. [PMID: 29812937 PMCID: PMC6150695 DOI: 10.1021/acs.jpcb.8b02250] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
![]()
Conformational
fluctuations and rotational tumbling of proteins
can be experimentally accessed with nuclear spin relaxation experiments.
However, interpretation of molecular dynamics from the experimental
data is often complicated, especially for molecules with anisotropic
shape. Here, we apply classical molecular dynamics simulations to
interpret the conformational fluctuations and rotational tumbling
of proteins with arbitrarily anisotropic shape. The direct calculation
of spin relaxation times from simulation data did not reproduce the
experimental data. This was successfully corrected by scaling the
overall rotational diffusion coefficients around the protein inertia
axes with a constant factor. The achieved good agreement with experiments
allowed the interpretation of the internal and overall dynamics of
proteins with significantly anisotropic shape. The overall rotational
diffusion was found to be Brownian, having only a short subdiffusive
region below 0.12 ns. The presented methodology can be applied to
interpret rotational dynamics and conformation fluctuations of proteins
with arbitrary anisotropic shape. However, a water model with more
realistic dynamical properties is probably required for intrinsically
disordered proteins.
Collapse
Affiliation(s)
- O H Samuli Ollila
- Research Program in Structural Biology and Biophysics, Institute of Biotechnology , University of Helsinki , 00014 Helsinki , Finland.,Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 117 20 Prague 6 , Czech Republic
| | - Harri A Heikkinen
- Research Program in Structural Biology and Biophysics, Institute of Biotechnology , University of Helsinki , 00014 Helsinki , Finland
| | - Hideo Iwaï
- Research Program in Structural Biology and Biophysics, Institute of Biotechnology , University of Helsinki , 00014 Helsinki , Finland
| |
Collapse
|
30
|
Methyl NMR spectroscopy: Measurement of dynamics in viral RNA-directed RNA polymerases. Methods 2018; 148:100-114. [PMID: 29857193 DOI: 10.1016/j.ymeth.2018.05.021] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Revised: 05/22/2018] [Accepted: 05/23/2018] [Indexed: 11/23/2022] Open
Abstract
Measurement of nuclear spin relaxation provides a powerful approach to access information about biomolecular conformational dynamics over several orders of magnitude in timescale. In several cases this knowledge in combination with spatial information from three-dimensional structures yields unique insight into protein stability and the kinetics and thermodynamics of their interactions and function. However, due to intrinsic difficulties in studying large systems using solution state nuclear magnetic resonance (NMR) approaches, until recently these measurements were limited to small-to-medium-sized systems. However, the development of a wide range of novel strategies that allow the selective isotope labeling of methyl groups in proteins have allowed the exploitation of the unique relaxation properties of this spin-system. This has in turn enabled the extension of NMR approaches to high molecular weight proteins including a variety of enzymes and their complexes. Here, we recount our experiences in obtaining assignments of the methyl resonances for two representative members of a class of RNA-directed RNA polymerases (RdRps) encoded by bacteriophages of the Cystoviridae family. We demonstrate the utility of these methyl probes, limited in number for one case and more numerous for the other, to investigate the conformational dynamics of RdRps on the fast (ps-ns) and slow (μs-ms) timescales.
Collapse
|
31
|
Hoffmann F, Xue M, Schäfer LV, Mulder FAA. Narrowing the gap between experimental and computational determination of methyl group dynamics in proteins. Phys Chem Chem Phys 2018; 20:24577-24590. [DOI: 10.1039/c8cp03915a] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
A synergistic analysis of methyl NMR relaxation data and MD simulations identifies ways to improve studies of protein dynamics.
Collapse
Affiliation(s)
- Falk Hoffmann
- Theoretical Chemistry
- Ruhr-University Bochum
- D-44780 Bochum
- Germany
| | - Mengjun Xue
- Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry
- University of Aarhus
- DK-8000 Aarhus
- Denmark
| | - Lars V. Schäfer
- Theoretical Chemistry
- Ruhr-University Bochum
- D-44780 Bochum
- Germany
| | - Frans A. A. Mulder
- Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry
- University of Aarhus
- DK-8000 Aarhus
- Denmark
| |
Collapse
|