Unraveling motion in proteins by combining NMR relaxometry and molecular dynamics simulations: A case study on ubiquitin

Candide Champion,Marc Lehner,Albert A. Smith,Fabien Ferrage,Nicolas Bolik-Coulon,Sereina Riniker
DOI: https://doi.org/10.1063/5.0188416
IF: 4.304
2024-03-11
The Journal of Chemical Physics
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.
chemistry, physical,physics, atomic, molecular & chemical
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to analyze the motions in proteins by combining nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations, especially on the time scale from picoseconds to nanoseconds. Specifically, the researchers focus on how to use these combined analysis methods to obtain descriptions of the internal motion mechanisms of proteins and the influence of these motions on protein functions. The paper takes ubiquitin as a test system and focuses on the motions of its aliphatic side chains because the backbone of this protein is relatively rigid. In the paper, five computational methods are compared. These methods allow for the combined analysis of data from MD simulations and NMR relaxation experiments, aiming to explore the relative advantages and application scopes of these methods in interpreting MD simulation dynamics and show how they can be used to analyze the dynamics of ubiquitin. In addition, through simulations using eight different combinations of force fields and water models, the study also explores how different metrics can be used to identify the deficiencies of force fields. This not only helps to verify the ability of the current state - of - the - art biomolecular force fields to predict dynamics, but also provides a perspective for the parameterization and verification of future force fields.