skip to main content


Title: Dynamical Models of Chemical Exchange in Nuclear Magnetic Resonance Spectroscopy
ABSTRACT Chemical exchange line broadening is an important phenomenon in nuclear magnetic resonance (NMR) spectroscopy, in which a nuclear spin experiences more than one magnetic environment as a result of chemical or conformational changes of a molecule. The dynamic process of chemical exchange strongly affects the sensitivity and resolution of NMR experiments and increasingly provides a powerful probe of the interconversion between chemical and conformational states of proteins, nucleic acids, and other biologic macromolecules. A simple and often used theoretic description of chemical exchange in NMR spectroscopy is based on an idealized 2-state jump model (the random phase or telegraph signal). However, chemical exchange can also be represented as a barrier crossing event that can be modeled by using chemical reaction rate theory. The timescale of crossing is determined by the barrier height, the temperature, and the dissipation modeled as collisional or frictional damping. This tutorial explores the connection between the NMR theory of chemical exchange line broadening and strong collision models for chemical kinetics in statistical mechanics. Theoretic modeling and numeric simulation are used to map the rate of barrier crossing dynamics of a particle on a potential energy surface to the chemical exchange relaxation rate constant. By developing explicit models for the exchange dynamics, the tutorial aims to elucidate the underlying dynamical processes that give rise to the rich phenomenology of chemical exchange observed in NMR spectroscopy. Software for generating and analyzing the numeric simulations is provided in the form of Python and Fortran source codes.  more » « less
Award ID(s):
1900416
NSF-PAR ID:
10339187
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
The Biophysicist
ISSN:
2578-6970
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Mechanistic understanding of DNA recombination in the Cre-loxPsystem has largely been guided by crystallographic structures of tetrameric synaptic complexes. Those studies have suggested a role for protein conformational dynamics that has not been well characterized at the atomic level. We used solution nuclear magnetic resonance (NMR) spectroscopy to discover the link between intrinsic flexibility and function in Cre recombinase. Transverse relaxation-optimized spectroscopy (TROSY) NMR spectra show the N-terminal and C-terminal catalytic domains (CreNTDand CreCat) to be structurally independent. Amide15N relaxation measurements of the CreCatdomain reveal fast-timescale dynamics in most regions that exhibit conformational differences in active and inactive Cre protomers in crystallographic tetramers. However, the C-terminal helix αN, implicated in assembly of synaptic complexes and regulation of DNA cleavage activity viatransprotein–protein interactions, is unexpectedly rigid in free Cre. Chemical shift perturbations and intra- and intermolecular paramagnetic relaxation enhancement (PRE) NMR data reveal an alternative autoinhibitory conformation for the αN region of free Cre, wherein it packsin cisover the protein DNA binding surface and active site. Moreover, binding toloxPDNA induces a conformational change that dislodges the C terminus, resulting in acis-to-transswitch that is likely to enable protein–protein interactions required for assembly of recombinogenic Cre intasomes. These findings necessitate a reexamination of the mechanisms by which this widely utilized gene-editing tool selects target sites, avoids spurious DNA cleavage activity, and controls DNA recombination efficiency.

     
    more » « less
  2. Proteins and nucleic acids participate in essentially every biochemical process in living organisms, and the elucidation of their structure and motions is essential for our understanding how these molecular machines perform their function. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful versatile technique that provides critical information on the molecular structure and dynamics. Spin-relaxation data are used to determine the overall rotational diffusion and local motions of biological macromolecules, while residual dipolar couplings (RDCs) reveal local and long-range structural architecture of these molecules and their complexes. This information allows researchers to refine structures of proteins and nucleic acids and provides restraints for molecular docking. Several software packages have been developed by NMR researchers in order to tackle the complicated experimental data analysis and structure modeling. However, many of them are offline packages or command-line applications that require users to set up the run time environment and also to possess certain programming skills, which inevitably limits accessibility of this software to a broad scientific community. Here we present new science gateways designed for NMR/structural biology community that address these current limitations in NMR data analysis. Using the GenApp technology for scientific gateways (https://genapp.rocks), we successfully transformed ROTDIF and ALTENS, two offline packages for bio-NMR data analysis, into science gateways that provide advanced computational functionalities, cloud-based data management, and interactive 2D and 3D plotting and visualizations. Furthermore, these gateways are integrated with molecular structure visualization tools (Jmol) and with gateways/engines (SASSIE-web) capable of generating huge computer-simulated structural ensembles of proteins and nucleic acids. This enables researchers to seamlessly incorporate conformational ensembles into the analysis in order to adequately take into account structural heterogeneity and dynamic nature of biological macromolecules. ROTDIF-web offers a versatile set of integrated modules/tools for determining and predicting molecular rotational diffusion tensors and model-free characterization of bond dynamics in biomacromolecules and for docking of molecular complexes driven by the information extracted from NMR relaxation data. ALTENS allows characterization of the molecular alignment under anisotropic conditions, which enables researchers to obtain accurate local and long-range bond-vector restraints for refining 3-D structures of macromolecules and their complexes. We will describe our experience bringing our programs into GenApp and illustrate the use of these gateways for specific examples of protein systems of high biological significance. We expect these gateways to be useful to structural biologists and biophysicists as well as NMR community and to stimulate other researchers to share their scientific software in a similar way. 
    more » « less
  3. ABSTRACT

    The 70 kDa heat shock proteins (Hsp70) are a family of molecular chaperones involved in protein folding, aggregate prevention, and protein disaggregation. They consist of the substrate‐binding domain (SBD) that binds client substrates, and the nucleotide‐binding domain (NBD), whose cycles of nucleotide hydrolysis and exchange underpin the activity of the chaperone. To characterize the structure–function relationships that link the binding state of the NBD to its conformational behavior, we analyzed the dynamics of the NBD of the Hsp70 chaperone fromBos taurus(PDB 3C7N:B) by all‐atom canonical molecular dynamics simulations. It was found that essential motions within the NBD fall into three major classes: the mutual class, reflecting tendencies common to all binding states, and the ADP‐ and ATP‐unique classes, which reflect conformational trends that are unique to either the ADP‐ or ATP‐bound states, respectively. “Mutual” class motions generally describe “in‐plane” and/or “out‐of‐plane” (scissor‐like) rotation of the subdomains within the NBD. This result is consistent with experimental nuclear magnetic resonance data on the NBD. The “unique” class motions target specific regions on the NBD, usually surface loops or sites involved in nucleotide binding and are, therefore, expected to be involved in allostery and signal transmission. For all classes, and especially for those of the “unique” type, regions of enhanced mobility can be identified; these are termed “hot spots,” and their locations generally parallel those found by NMR spectroscopy. The presence of magnesium and potassium cations in the nucleotide‐binding pocket was also found to influence the dynamics of the NBD significantly. Proteins 2015; 83:282–299. © 2014 Wiley Periodicals, Inc.

     
    more » « less
  4. Abstract

    Many-body dynamical models in which Boltzmann statistics can be derived directly from the underlying dynamical laws without invoking the fundamental postulates of statistical mechanics are scarce. Interestingly, one such model is found in econophysics and in chemistry classrooms: the money game, in which players exchange money randomly in a process that resembles elastic intermolecular collisions in a gas, giving rise to the Boltzmann distribution of money owned by each player. Although this model offers a pedagogical example that demonstrates the origins of Boltzmann statistics, such demonstrations usually rely on computer simulations. In fact, a proof of the exponential steady-state distribution in this model has only become available in recent years. Here, we study this random money/energy exchange model and its extensions using a simple mean-field-type approach that examines the properties of the one-dimensional random walk performed by one of its participants. We give a simple derivation of the Boltzmann steady-state distribution in this model. Breaking the time-reversal symmetry of the game by modifying its rules results in non-Boltzmann steady-state statistics. In particular, introducing ‘unfair’ exchange rules in which a poorer player is more likely to give money to a richer player than to receive money from that richer player, results in an analytically provable Pareto-type power-law distribution of the money in the limit where the number of players is infinite, with a finite fraction of players in the ‘ground state’ (i.e. with zero money). For a finite number of players, however, the game may give rise to a bimodal distribution of money and to bistable dynamics, in which a participant’s wealth jumps between poor and rich states. The latter corresponds to a scenario where the player accumulates nearly all the available money in the game. The time evolution of a player’s wealth in this case can be thought of as a ‘chemical reaction’, where a transition between ‘reactants’ (rich state) and ‘products’ (poor state) involves crossing a large free energy barrier. We thus analyze the trajectories generated from the game using ideas from the theory of transition paths and highlight non-Markovian effects in the barrier crossing dynamics.

     
    more » « less
  5. In this work, the second-order kinetics of molecules exchanging between chemically distinct microenvironments, such as those found in nanoemulsions, is studied using nuclear magnetic resonance (NMR). A unique aspect of NMR exchange studies in nanoemulsions is that the difference in molecular resonance frequencies between the two phases, which determines whether the exchange is fast, intermediate, or slow on the NMR timescale, can depend upon the emulsion droplet composition, which is also determined by the kinetic exchange constants themselves. Within the fast-exchange regime, changes in resonance frequencies and line widths with dilution were used to extract the exchange rate constants from the NMR spectra in a manner analogous to determining the kinetic parameters in NMR ligand binding experiments. As a demonstration, the kinetic exchange parameters of isoflurane release from an emulsification of isoflurane and perflurotributylamine (FC43) were determined using NMR dilution and diffusion studies. 
    more » « less