1
|
Banayan NE, Hsu A, Hunt JF, Palmer AG, Friesner RA. Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations. J Chem Theory Comput 2024; 20:6316-6327. [PMID: 38957960 DOI: 10.1021/acs.jctc.4c00378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.
Collapse
Affiliation(s)
- Nooriel E Banayan
- Department of Biological Sciences, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Andrew Hsu
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John F Hunt
- Department of Biological Sciences, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Arthur G Palmer
- Department of Biochemistry and Molecular Biophysics, Columbia University, 701 West 168th Street, New York, New York 10032, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
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
|
Gilardoni I, Fröhlking T, Bussi G. Boosting Ensemble Refinement with Transferable Force-Field Corrections: Synergistic Optimization for Molecular Simulations. J Phys Chem Lett 2024; 15:1204-1210. [PMID: 38272001 DOI: 10.1021/acs.jpclett.3c03423] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2024]
Abstract
A novel method combining the force-field fitting approach and ensemble refinement by the maximum entropy principle is presented. Its formulation allows us to continuously interpolate between these two methods, which can thus be interpreted as two limiting cases. A cross-validation procedure enables us to correctly assess the relative weight of both of them, distinguishing scenarios in which the combined approach is meaningful from those in which either ensemble refinement or force-field fitting separately prevails. The efficacy of their combination is examined for a realistic case study of RNA oligomers. Within the new scheme, molecular dynamics simulations are integrated with experimental data provided by nuclear magnetic resonance measures. We show that force-field corrections are in general superior when applied to the appropriate force-field terms but are automatically discarded by the method when applied to inappropriate force-field terms.
Collapse
Affiliation(s)
- Ivan Gilardoni
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| | - Thorben Fröhlking
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
4
|
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
|