Chemical Exchange

Conformational penalties: New insights into nucleic acid recognition

The energy cost accompanying changes in the structures of nucleic acids when they bind partner molecules is a significant but underappreciated thermodynamic contribution to binding affinity and specificity. This review highlights recent advances in measuring conformational penalties and determining their contribution to the recognition, folding, and regulatory activities of nucleic acids. Notable progress includes methods for measuring and structurally characterizing lowly populated conformational states, obtaining ensemble information in high throughput, for large macromolecular assemblies, and in complex cellular environments. Additionally, quantitative and predictive thermodynamic models have been developed that relate conformational penalties to nucleic acid-protein association and cellular activity. These studies underscore the crucial role of conformational penalties in nucleic acid recognition.

Insights into Ligand-Mediated Activation of an Oligomeric Ring-Shaped Gene-Regulatory Protein from Solution- and Solid-State NMR

The 91 kDa oligomeric ring-shaped ligand binding protein TRAP (trp RNA binding attenuation protein) regulates the expression of a series of genes involved in tryptophan (Trp) biosynthesis in bacilli. When cellular Trp levels rise, the free amino acid binds to sites buried in the interfaces between each of the 11 (or 12, depending on the species) protomers in the ring. Crystal structures of Trp-bound TRAP show the Trp ligands are sequestered from solvent by a pair of loops from adjacent protomers that bury the bound ligand via polar contacts to several threonine residues. Binding of the Trp ligands occurs cooperatively, such that successive binding events occur with higher apparent affinity but the structural basis for this cooperativity is poorly understood. We used solution methyl-TROSY NMR relaxation experiments focused on threonine and isoleucine sidechains, as well as magic angle spinning solid-state NMR 13C-13C and 15N-13C chemical shift correlation spectra on uniformly labeled samples recorded at 800 and 1200 MHz, to characterize the structure and dynamics of the protein. Methyl 13C relaxation dispersion experiments on ligand-free apo TRAP revealed concerted exchange dynamics on the µs-ms time scale, consistent with transient sampling of conformations that could allow ligand binding. Cross-correlated relaxation experiments revealed widespread disorder on fast timescales. Chemical shifts for methyl-bearing side chains in apo- and Trp-bound TRAP revealed subtle changes in the distribution of sampled sidechain rotameric states. These observations reveal a pathway and mechanism for induced conformational changes to generate homotropic Trp-Trp binding cooperativity.

Gating residues govern ligand unbinding kinetics from the buried cavity in HIF-2α PAS-B

While transcription factors have been generally perceived as "undruggable," an exception is the HIF-2 hypoxia-inducible transcription factor, which contains an internal cavity that is sufficiently large to accommodate a range of small-molecules, including the therapeutically used inhibitor belzutifan. Given the relatively long ligand residence times of these small molecules and the lack of any experimentally observed pathway connecting the cavity to solvent, there has been great interest in understanding how these drug ligands exit the buried receptor cavity. Here, we focus on the relevant PAS-B domain of hypoxia-inducible factor 2α (HIF-2α) and examine how one such small molecule (THS-017) exits from the buried cavity within this domain on the seconds-timescale using atomistic simulations and ZZ-exchange NMR. To enable the simulations, we applied the weighted ensemble path sampling strategy, which generates continuous pathways for a rare-event process [e.g., ligand (un)binding] with rigorous kinetics in orders of magnitude less computing time compared to conventional simulations. Results reveal the formation of an encounter complex intermediate and two distinct classes of pathways for ligand exit. Based on these pathways, we identified two pairs of conformational gating residues in the receptor: one for the major class (N288 and S304) and another for the minor class (L272 and M309). ZZ-exchange NMR validated the kinetic importance of N288 for ligand unbinding. Our results provide an ideal simulation dataset for rational manipulation of ligand unbinding kinetics.

Slow Conformational Exchange between Partially Folded and Near-Native States of Ubiquitin: Evidence for a Multistate Folding Model

The mechanism by which small proteins fold, i.e., via intermediates or via a two-state mechanism, is a subject of intense investigation. Intermediate states in the folding pathways of these proteins are sparsely populated due to transient lifetimes under normal conditions rendering them transparent to a majority of the biophysical methods employed for structural, thermodynamic, and kinetic characterization, which attributes are essential for understanding the cooperative folding/unfolding of such proteins. Dynamic NMR spectroscopy has enabled the characterization of folding intermediates of ubiquitin that exist in equilibrium under conditions of low pH and denaturants. At low pH, an unlocked state defined as N' is in fast exchange with an invisible state, U″, as observed by CEST NMR. Addition of urea to ubiquitin at pH 2 creates two new states F' and U', which are in slow exchange (kF'→U' = 0.14 and kU'→F' = 0.28 s-1) as indicated by longitudinal ZZ-magnetization exchange spectroscopy. High-resolution solution NMR structures of F' show it to be in an "unlocked" conformation with measurable changes in rotational diffusion, translational diffusion, and rotational correlational times. U' is characterized by the presence of just the highly conserved N-terminal β1-β2 hairpin. The folding of ubiquitin is cooperative and is nucleated by the formation of an N-terminal β-hairpin followed by significant hydrophobic collapse of the protein core resulting in the formation of bulk of the secondary structural elements stabilized by extensive tertiary contacts. U' and F' may thus be described as early and late folding intermediates in the ubiquitin folding pathway.

NMR and Single-Molecule FRET Insights into Fast Protein Motions and Their Relation to Function

Proteins often undergo large-scale conformational transitions, in which secondary and tertiary structure elements (loops, helices, and domains) change their structures or their positions with respect to each other. Simple considerations suggest that such dynamics should be relatively fast, but the functional cycles of many proteins are often relatively slow. Sophisticated experimental methods are starting to tackle this dichotomy and shed light on the contribution of large-scale conformational dynamics to protein function. In this review, we focus on the contribution of single-molecule Förster resonance energy transfer and nuclear magnetic resonance (NMR) spectroscopies to the study of conformational dynamics. We briefly describe the state of the art in each of these techniques and then point out their similarities and differences, as well as the relative strengths and weaknesses of each. Several case studies, in which the connection between fast conformational dynamics and slower function has been demonstrated, are then introduced and discussed. These examples include both enzymes and large protein machines, some of which have been studied by both NMR and fluorescence spectroscopies.

Increasing the accuracy of exchange parameters reporting on slow dynamics by performing CEST experiments with 'high' B1 fields

Over the last decade chemical exchange saturation transfer (CEST) NMR methods have emerged as powerful tools to characterize biomolecular conformational dynamics occurring between a visible major state and 'invisible' minor states. The ability of the CEST experiment to detect these minor states, and provide precise exchange parameters, hinges on using appropriate B1 field strengths during the saturation period. Typically, a pair of B1 fields with ω1 (=2πB1) values around the exchange rate kex are chosen. Here we show that the transverse relaxation rate of the minor state resonance (R2,B) also plays a crucial role in determining the B1 fields that lead to the most informative datasets. Using [Formula: see text]  ≥ kex, to guide the choice of B1, instead of kex, leads to data wherefrom substantially more accurate exchange parameters can be derived. The need for higher B1 fields, guided by K, is demonstrated by studying the conformational exchange in two mutants of the 71 residue FF domain with kex ∼ 11 s-1 and ∼ 72 s-1, respectively. In both cases analysis of CEST datasets recorded using B1 field values guided by kex lead to imprecise exchange parameters, whereas using B1 values guided by K resulted in precise site-specific exchange parameters. The conclusions presented here will be valuable while using CEST to study slow processes at sites with large intrinsic relaxation rates, including carbonyl sites in small to medium sized proteins, amide 15N sites in large proteins and when the minor state dips are broadened due to exchange among the minor states.

Insights into Ligand-Mediated Activation of an Oligomeric Ring-Shaped Gene-Regulatory Protein from Solution- and Solid-State NMR

The 91 kDa oligomeric ring-shaped ligand binding protein TRAP (trp RNA binding attenuation protein) regulates the expression of a series of genes involved in tryptophan (Trp) biosynthesis in bacilli. When cellular Trp levels rise, the free amino acid binds to sites buried in the interfaces between each of the 11 (or 12, depending on the species) protomers in the ring. Crystal structures of Trp-bound TRAP show the Trp ligands are sequestered from solvent by a pair of loops from adjacent protomers that bury the bound ligand via polar contacts to several threonine residues. Binding of the Trp ligands occurs cooperatively, such that successive binding events occur with higher apparent affinity but the structural basis for this cooperativity is poorly understood. We used solution methyl-TROSY NMR relaxation experiments focused on threonine and isoleucine sidechains, as well as magic angle spinning solid-state NMR 13C-13C and 15N-13C chemical shift correlation spectra on uniformly labeled samples recorded at 800 and 1200 MHz, to characterize the structure and dynamics of the protein. Methyl 13C relaxation dispersion experiments on ligand-free apo TRAP revealed concerted exchange dynamics on the μs-ms time scale, consistent with transient sampling of conformations that could allow ligand binding. Cross-correlated relaxation experiments revealed widespread disorder on fast timescales. Chemical shifts for methyl-bearing side chains in apo- and Trp-bound TRAP revealed subtle changes in the distribution of sampled sidechain rotameric states. These observations reveal a pathway and mechanism for induced conformational changes to generate homotropic Trp-Trp binding cooperativity.

Comparison of the Backbone Dynamics of Dehaloperoxidase-Hemoglobin Isoenzymes

Dehaloperoxidase (DHP) is a multifunctional hemeprotein with a functional switch generally regulated by the chemical class of the substrate. Its two isoforms, DHP-A and DHP-B, differ by only five amino acids and have an almost identical protein fold. However, the catalytic efficiency of DHP-B for oxidation by a peroxidase mechanism ranges from 2- to 6-fold greater than that of DHP-A depending on the conditions. X-ray crystallography has shown that many substrates and ligands have nearly identical binding in the two isoenzymes, suggesting that the difference in catalytic efficiency could be due to differences in the conformational dynamics. We compared the backbone dynamics of the DHP isoenzymes at pH 7 through heteronuclear relaxation dynamics at 11.75, 16.45, and 19.97 T in combination with four 300 ns MD simulations. While the overall dynamics of the isoenzymes are similar, there are specific local differences in functional regions of each protein. In DHP-A, Phe35 undergoes a slow chemical exchange between two conformational states likely coupled to a swinging motion of Tyr34. Moreover, Asn37 undergoes fast chemical exchange in DHP-A. Given that Phe35 and Asn37 are adjacent to Tyr34 and Tyr38, it is possible that their dynamics modulate the formation and migration of the active tyrosyl radicals in DHP-A at pH 7. Another significant difference is that both distal and proximal histidines have a 15-18% smaller S2 value in DHP-B, thus their greater flexibility could account for the higher catalytic activity. The distal histidine grants substrate access to the distal pocket. The greater flexibility of the proximal histidine could also accelerate H2O2 activation at the heme Fe by increased coupling of an amino acid charge relay to stabilize the ferryl Fe(IV) oxidation state in a Poulos-Kraut "push-pull"-type peroxidase mechanism.

Beyond slow two-state protein conformational exchange using CEST: applications to three-state protein interconversion on the millisecond timescale

Although NMR spectroscopy is routinely used to study the conformational dynamics of biomolecules, robust analyses of the data are challenged in cases where exchange is more complex than two-state, such as when a 'visible' major conformer exchanges with two 'invisible' minor states on the millisecond timescale. It is becoming increasingly clear that chemical exchange saturation transfer (CEST) NMR experiments that were initially developed to study systems undergoing slow interconversion are also sensitive to intermediate-fast timescale biomolecular conformational exchange. Here we investigate the utility of the amide 15N CEST experiment to characterise protein three-state exchange occurring on the millisecond timescale by studying the interconversion between the folded (F) state of the FF domain from human HYPA/FBP11 (WT FF) and two of its folding intermediates I1 and I2. Although 15N CPMG experiments are consistent with the F state interconverting with a single minor state on the millisecond timescale, 15N CEST data clearly establish an exchange process between F and a pair of minor states. A unique three-state exchange model cannot be obtained by analysis of 15N CEST data recorded at a single temperature. However, including the relative sign of the difference in the chemical shifts of the two minor states based on a simple two-state analysis of CEST data recorded at multiple temperatures, results in a robust three-state model in which the F, I1 and I2 states interconvert with each other on the millisecond timescale ( k e x , F I 1 ~ 550 s-1, k e x , F I 2 ~ 1200 s-1, k e x , I 1 I 2 ~ 5000 s-1), with I1 and I2 sparsely populated at ~ 0.15% and ~ 0.35%, respectively, at 15 °C. A computationally demanding grid-search of exchange parameter space is not required to extract the best-fit exchange parameters from the CEST data. The utility of the CEST experiment, thus, extends well beyond studies of conformers in slow exchange on the NMR chemical shift timescale, to include systems with interconversion rates on the order of thousands/second.

Measuring Hydroxyl Exchange Rates in Glycans Using a Synergistic Combination of Saturation Transfer and Relaxation Dispersion NMR

Molecular recognition events mediated by glycans play pivotal roles in controlling the fate of diverse biological processes such as cellular communication and the immune response. The affinity of glycans for their target receptors is governed primarily by the hydrogen bonds formed by hydroxyl groups decorating the glycan surface. Hydroxyl exchange rate constants are therefore vital parameters that report on glycan structure and dynamics. Here we present a strategy for characterizing hydroxyl hydrogen/deuterium (H/D) exchange in glycans that employs a synergistic combination of 13C chemical exchange saturation transfer (CEST) and Carr-Purcell-Meiboom-Gill relaxation dispersion (CPMG) NMR methods. We show that the combination of CEST and CPMG experiments facilitates the sensitive detection of the small (∼0.1 ppm) two-bond deuterium isotope shift on a 13C nucleus when the attached hydroxyl group fluctuates between protonated and deuterated states. This shift is leveraged for measuring site-specific kinetic H/D exchange rate constants as well as thermodynamic free energies of isotope fractionation. The CEST and CPMG modules are integrated with a selective J-cross-polarization scheme that provides the flexibility for rapid characterization of H/D exchange at a specific hydroxyl site. Moreover, our approach enables the precise isothermal measurement of hydroxyl exchange rate constants without the need for cumbersome isotope labeling. The H/D exchange rate constants of three different glycans assessed using this method highlight its potential for detecting transient intra- and intermolecular hydrogen bonds. In addition, the trends in H/D exchange rate constants establish site-specific steric accessibility as a key determinant of solvent exchange dynamics in glycans.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

Relaxation and single site multiple mutations to identify and control allosteric networks

Advances in Nuclear Magnetic Resonance (NMR) spectroscopy have allowed for the identification and characterization of movements in enzymes over the last 20 years that has also revealed the complexities of allosteric coupling. For example, many of the inherent movements of enzymes, and proteins in general, have been shown to be highly localized but nonetheless still coupled over long distances. Such partial couplings provide challenges to both identifying allosteric networks of dynamic communication and determining their roles in catalytic function. We have developed an approach to help identify and engineer enzyme function, called Relaxation And Single Site Multiple Mutations (RASSMM). This approach is a powerful extension of mutagenesis and NMR that is based on the observation that multiple mutations to a single site distal to the active site allosterically induces different effects to networks. Such an approach generates a panel of mutations that can also be subjected to functional studies in order to match catalytic effects with changes to coupled networks. In this review, the RASSMM approach is briefly outlined together with two applications that include cyclophilin-A and Biliverdin Reductase B.

Dynamic basis for dA•dGTP and dA•d8OGTP misincorporation via Hoogsteen base pairs

Replicative errors contribute to the genetic diversity needed for evolution but in high frequency can lead to genomic instability. Here, we show that DNA dynamics determine the frequency of misincorporating the A•G mismatch, and altered dynamics explain the high frequency of 8-oxoguanine (8OG) A•8OG misincorporation. NMR measurements revealed that Aanti•Ganti (population (pop.) of >91%) transiently forms sparsely populated and short-lived Aanti+•Gsyn (pop. of ~2% and kex = kforward + kreverse of ~137 s-1) and Asyn•Ganti (pop. of ~6% and kex of ~2,200 s-1) Hoogsteen conformations. 8OG redistributed the ensemble, rendering Aanti•8OGsyn the dominant state. A kinetic model in which Aanti+•Gsyn is misincorporated quantitatively predicted the dA•dGTP misincorporation kinetics by human polymerase β, the pH dependence of misincorporation and the impact of the 8OG lesion. Thus, 8OG increases replicative errors relative to G because oxidation of guanine redistributes the ensemble in favor of the mutagenic Aanti•8OGsyn Hoogsteen state, which exists transiently and in low abundance in the A•G mismatch.

Synergism between x-ray crystallography and NMR residual dipolar couplings in characterizing protein dynamics

The important role of structural dynamics in protein function is widely recognized. Thermal or B-factors and their anisotropy, seen in x-ray analysis of protein structures, report on the presence of atomic coordinate heterogeneity that can be attributed to motion. However, their quantitative evaluation in terms of protein dynamics by x-ray ensemble refinement remains challenging. NMR spectroscopy provides quantitative information on the amplitudes and time scales of motional processes. Unfortunately, with a few exceptions, the NMR data do not provide direct insights into the atomic details of dynamic trajectories. Residual dipolar couplings, measured by solution NMR, are very precise parameters reporting on the time-averaged bond-vector orientations and may offer the opportunity to derive correctly weighted dynamic ensembles of structures for cases where multiple high-resolution x-ray structures are available. Applications to the SARS-CoV-2 main protease, Mpro, and ubiquitin highlight this complementarity of NMR and crystallography for quantitative assessment of internal motions.

Multi-Site Conformational Exchange in the Synthetic Neomycin-Sensing Riboswitch Studied by 19 F NMR

The synthetic neomycin-sensing riboswitch interacts with its cognate ligand neomycin as well as with the related antibiotics ribostamycin and paromomycin. Binding of these aminoglycosides induces a very similar ground state structure in the RNA, however, only neomycin can efficiently repress translation initiation. The molecular origin of these differences has been traced back to differences in the dynamics of the ligand:riboswitch complexes. Here, we combine five complementary fluorine based NMR methods to accurately quantify seconds to microseconds dynamics in the three riboswitch complexes. Our data reveal complex exchange processes with up to four structurally different states. We interpret our findings in a model that shows an interplay between different chemical groups in the antibiotics and specific bases in the riboswitch. More generally, our data underscore the potential of 19 F NMR methods to characterize complex exchange processes with multiple excited states.

An NMR portrait of functional and dysfunctional allosteric cooperativity in cAMP-dependent protein kinase A

The cAMP-dependent protein kinase A (PKA) is the archetypical eukaryotic kinase. The catalytic subunit (PKA-C) structure is highly conserved among the AGC-kinase family. PKA-C is a bilobal enzyme with a dynamic N-lobe, harbouring the Adenosine-5'-triphosphate (ATP) binding site and a more rigid helical C-lobe. The substrate-binding groove resides at the interface of the two lobes. A distinct feature of PKA-C is the positive binding cooperativity between nucleotide and substrate. Several PKA-C mutations lead to the development of adenocarcinomas, myxomas, and other rare forms of liver tumours. Nuclear magnetic resonance (NMR) spectroscopy shows that these mutations disrupt the allosteric communication between the two lobes, causing a drastic decrease in binding cooperativity. The loss of cooperativity correlates with changes in substrate fidelity and reduced kinase affinity for the endogenous protein kinase inhibitor (PKI). The similarity between PKI and the inhibitory sequence of the kinase regulatory subunits suggests that the overall mechanism of regulation of the kinase may be disrupted. We surmise that a reduced or obliterated cooperativity may constitute a common trait for both orthosteric and allosteric mutations of PKA-C that may lead to dysregulation and disease.

Slow conformational changes in the rigid and highly stable chymotrypsin inhibitor 2

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used ¹⁵ N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1 °C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes. This article is protected by copyright. All rights reserved.

Site-based description of R1ρ relaxation in local reference frames

Nuclear magnetic spin relaxation in the presence of an applied radiofrequency field depends critically on chemical exchange processes that transfer nuclear spins between chemical or conformational environments with distinct resonance frequencies. Characterization of chemical exchange processes in R1ρ relaxation dispersion, CEST, and DEST experiments provides powerful insights into chemical and conformational kinetics of biological macromolecules. The present work reformulates expressions for magnetization evolution and the R1ρ relaxation rate constant by focussing on the orientations of the tilted rotating frames of reference for magnetization components in different sites, by treating the spin-locking field strength as a perturbation to free-precession evolution, and by applying the Homotopy Analysis and Laplace transform methods to approximate solutions to the Bloch-McConnell equations. The results provide an expression for R1ρ that is invariant to the topology of the kinetic scheme, an approximate equation for evolution of spin-locked z-magnetization, and an approach for effective simplification of chemical exchange topologies for 2- and N-site chemical exchange processes. The theoretical approach also provides an accurate approximation for relaxation during a constant-amplitude radiofrequency field in the absence of exchange.

Real-time Exchange of the Lipid-bound Intermediate and Post-fusion States of the HIV-1 gp41 Ectodomain

The envelope glycoprotein gp41 of the HIV-1 virus mediates its entry into the host cell. During this process, gp41 undergoes large conformational changes and the energy released in the remodeling events is utilized to overcome the barrier associated with fusing the viral and host membranes. Although the structural intermediates of this fusion process are attractive targets for drug development, no detailed high-resolution structural information or quantitative thermodynamic characterization are available. By measuring the dynamic equilibrium between the lipid-bound intermediate and the post-fusion six-helical bundle (6HB) states of the gp41 ectodomain in the presence of bilayer membrane mimetics, we derived both the reaction kinetics and energies associated with these two states by solution NMR spectroscopy. At equilibrium, an exchange time constant of about 12 seconds at 38 °C is observed, and the post-fusion conformation is energetically more stable than the lipid-bound state by 3.4 kcal mol-1. The temperature dependence of the kinetics indicates that the folding occurs through a high-energy transition state which may resemble a 5HB structure. The energetics and kinetics of gp41 folding in the context of membrane bilayers provide a molecular basis for an improved understanding of viral membrane fusion.

Global protein dynamics as communication sensors in peptide synthetase domains

Biological activity is governed by the timely redistribution of molecular interactions, and static structural snapshots often appear insufficient to provide the molecular determinants that choreograph communication. This conundrum applies to multidomain enzymatic systems called nonribosomal peptide synthetases (NRPSs), which assemble simple substrates into complex metabolites, where a dynamic domain organization challenges rational design to produce new pharmaceuticals. Using a nuclear magnetic resonance (NMR) atomic-level readout of biochemical transformations, we demonstrate that global structural fluctuations help promote substrate-dependent communication and allosteric responses, and impeding these global dynamics by a point-site mutation hampers allostery and molecular recognition. Our results establish global structural dynamics as sensors of molecular events that can remodel domain interactions, and they provide new perspectives on mechanisms of allostery, protein communication, and NRPS synthesis.

Dynamical Models of Chemical Exchange in Nuclear Magnetic Resonance Spectroscopy

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 inter-conversion between chemical and conformational states of proteins, nucleic acids, and other biological macromolecules. A simple and often used theoretical description of chemical exchange in NMR spectroscopy is based on an idealized two-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 using chemical reaction rate theory. The time scale 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. Theoretical modeling and numerical 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 numerical simulations is provided in the form of Python and Fortran source codes.

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Conformational Dynamics of the Hepatitis B Virus Pre-genomic RNA on Multiple Time Scales: Implications for Viral Replication

Human hepatitis B virus (HBV) replication is initiated by the binding of the viral polymerase (P) to epsilon (ε), an ≈85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 5'-end of the pre-genomic RNA (pgRNA). This interaction triggers P and pgRNA packaging and protein-primed reverse transcription and is therefore an attractive therapeutic target. Our recent nuclear magnetic resonance (NMR) structure of ε provides a useful starting point toward a detailed understanding of HBV replication, and hints at the functional importance of ε dynamics. Here, we present a detailed description of ε motions on the ps to ns and μs to ms time scales by NMR spin relaxation and relaxation dispersion, respectively. We also carried out molecular dynamics simulations to provide additional insight into ε conformational dynamics. These data outline a series of complex motions on multiple time scales within ε. Moreover, these motions occur in mostly conserved nucleotides from structural regions (i.e., priming loop, pseudo-triloop, and U43 bulge) that biochemical and mutational studies have shown to be essential for P binding, P-pgRNA packaging, protein-priming, and DNA synthesis. Taken together, our work implicates RNA dynamics as an integral feature that governs HBV replication.

Conformational Dynamics of Intrinsically Disordered Proteins Regulate Biomolecular Condensate Chemistry

Motions in biomolecules are critical for biochemical reactions. In cells, many biochemical reactions are executed inside of biomolecular condensates formed by ultradynamic intrinsically disordered proteins. A deep understanding of the conformational dynamics of intrinsically disordered proteins in biomolecular condensates is therefore of utmost importance but is complicated by diverse obstacles. Here we review emerging data on the motions of intrinsically disordered proteins inside of liquidlike condensates. We discuss how liquid–liquid phase separation modulates internal motions across a wide range of time and length scales. We further highlight the importance of intermolecular interactions that not only drive liquid–liquid phase separation but appear as key determinants for changes in biomolecular motions and the aging of condensates in human diseases. The review provides a framework for future studies to reveal the conformational dynamics of intrinsically disordered proteins in the regulation of biomolecular condensate chemistry.

1H R1ρ relaxation dispersion experiments in aromatic side chains

Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on microsecond to millisecond timescales can be studied by relaxation dispersion experiments that attenuate conformational exchange contributions to the transverse relaxation rate by varying the refocusing frequency of applied radio-frequency fields implemented as either CPMG pulse trains or continuous spin-lock periods. Here we present an aromatic ¹H R_1ρ relaxation dispersion experiment enabling studies of two to three times faster exchange processes than achievable by existing experiments for aromatic side chains. We show that site-specific isotope labeling schemes generating isolated ¹H-¹³C spin pairs with vicinal ²H-¹²C moieties are necessary to avoid anomalous relaxation dispersion profiles caused by Hartmann-Hahn matching due to the ³J_HH couplings and limited chemical shift differences among ¹H spins in phenylalanine, tyrosine and the six-ring moiety of tryptophan. This labeling pattern is sufficient in that remote protons do not cause additional complications. We validated the approach by measuring ring-flip kinetics in the small protein GB1. The determined rate constants, k_flip, agree well with previous results from ¹³C R_1ρ relaxation dispersion experiments, and yield ¹H chemical shift differences between the two sides of the ring in good agreement with values measured under slow-exchange conditions. The aromatic¹H R_1ρ relaxation dispersion experiment in combination with the site-selective ¹H-¹³C/²H-¹²C labeling scheme enable measurement of exchange rates up to k_ex = 2k_flip = 80,000 s^-1, and serve as a useful complement to previously developed ¹³C-based methods.

The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways

Conformational dynamics play critical roles in protein folding, misfolding, function, misfunction, and aggregation. While detecting and studying the different conformational states populated by protein molecules on their free energy surfaces (FESs) remain a challenge, NMR spectroscopy has emerged as an invaluable experimental tool to explore the FES of a protein, as conformational dynamics can be probed at atomic resolution over a wide range of timescales. Here, we use chemical exchange saturation transfer (CEST) to detect "invisible" minor states on the energy landscape of the A39G mutant FF domain that exhibited "two-state" folding kinetics in traditional experiments. Although CEST has mostly been limited to studies of processes with rates between ∼5 to 300 s^-1 involving sparse states with populations as low as ∼1%, we show that the line broadening that is often associated with minor state dips in CEST profiles can be exploited to inform on additional conformers, with lifetimes an order of magnitude shorter and populations close to 10-fold smaller than what typically is characterized. Our analysis of CEST profiles that exploits the minor state linewidths of the 71-residue A39G FF domain establishes a folding mechanism that can be described in terms of a four-state exchange process between interconverting states spanning over two orders of magnitude in timescale from ∼100 to ∼15,000 μs. A similar folding scheme is established for the wild-type domain as well. The study shows that the folding of this small domain proceeds through a pair of sparse, partially structured intermediates via two discrete pathways on a volcano-shaped FES.

Advances in NMR Spectroscopy of Weakly Aligned Biomolecular Systems

The measurement and application of residual dipolar couplings (RDCs) in solution NMR studies of biological macromolecules has become well established over the past quarter of a century. Numerous methods for generating the requisite anisotropic orientational molecular distribution have been demonstrated, each with its specific strengths and weaknesses. In parallel, an enormous number of pulse schemes have been introduced to measure the many different types of RDCs, ranging from the most widely measured backbone amide ¹⁵N-¹H RDCs, to ¹H-¹H RDCs and couplings between low-γ nuclei. Applications of RDCs range from structure validation and refinement to the determination of relative domain orientations, the measurement of backbone and domain motions, and de novo structure determination. Nevertheless, it appears that the power of the RDC methodology remains underutilized. This review aims to highlight the practical aspects of sample preparation and RDC measurement while describing some of the most straightforward applications that take advantage of the exceptionally precise information contained in such data. Some emphasis will be placed on more recent developments that enable the accurate measurement of RDCs in larger systems, which is key to the ongoing shift in focus of biological NMR spectroscopy from structure determination toward gaining improved understanding of how molecular flexibility drives protein function.

Analysis of conformational exchange processes using methyl-TROSY-based Hahn echo measurements of quadruple-quantum relaxation

Transverse nuclear spin relaxation is a sensitive probe of chemical exchange on timescales on the order of microseconds to milliseconds. Here we present an experiment for the simultaneous measurement of the relaxation rates of two quadruple-quantum transitions in 13 CH3 -labelled methyl groups. These coherences are protected against relaxation by intra-methyl dipolar interactions and so have unexpectedly long lifetimes within perdeuterated biomacromolecules. However, these coherences also have an order of magnitude higher sensitivity to chemical exchange broadening than lower order coherences and therefore provide ideal probes of dynamic processes. We show that analysis of the static magnetic field dependence of zero-, double- and quadruple-quantum Hahn echo relaxation rates provides a robust indication of chemical exchange and can determine the signed relative magnitudes of proton and carbon chemical shift differences between ground and excited states. We also demonstrate that this analysis can be combined with established Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements, providing improved precision in parameter estimates, particularly in the determination of 1 H chemical shift differences.

Reverse-dialysis can be misleading for drug release studies in emulsions as demonstrated by NMR dilution experiments

Emulsions are an important class of carriers for the delivery of hydrophobic drugs. While knowledge of drug release kinetics is critical to optimizing drug carrying emulsions, there remain many open questions about the validity of standard characterization methods such as the commonly used reverse-dialysis. In this paper, the kinetic parameters of isoflurane release in perfluorotributylamine emulsions determined from both reverse-dialysis and nuclear magnetic resonance (NMR) dilution experiments are compared. The NMR-determined kinetic parameters of isoflurane release were found to be approximately seven orders of magnitude larger than those determined from conventional reverse-dialysis and were also shown to be consistent with prior in vivo observations of the anesthetization of rats.

Nuclear singlet relaxation by chemical exchange

The population imbalance between nuclear singlet states and triplet states of strongly coupled spin-1/2 pairs, also known as nuclear singlet order, is well protected against several common relaxation mechanisms. We study the nuclear singlet relaxation of ¹³C pairs in aqueous solutions of 1,2-¹³C₂ squarate over a range of pH values. The ¹³C singlet order is accessed by introducing ¹⁸O nuclei in order to break the chemical equivalence. The squarate dianion is in chemical equilibrium with hydrogen-squarate (SqH^-) and squaric acid (SqH₂) characterized by the dissociation constants pK₁ = 1.5 and pK₂ = 3.4. Surprisingly, we observe a striking increase in the singlet decay time constants T_S when the pH of the solution exceeds ∼10, which is far above the acid-base equilibrium points. We derive general rate expressions for chemical-exchange-induced nuclear singlet relaxation and provide a qualitative explanation of the T_S behavior of the squarate dianion. We identify a kinetic contribution to the singlet relaxation rate constant, which explicitly depends on kinetic rate constants. Qualitative agreement is achieved between the theory and the experimental data. This study shows that infrequent chemical events may have a strong effect on the relaxation of nuclear singlet order.

Determining chemical exchange rate constants in nanoemulsions using nuclear magnetic resonance

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.

Deuteron Chemical Exchange Saturation Transfer for the Detection of Slow Motions in Rotating Solids

We utilized the ²H Chemical Exchange Saturation Transfer (CEST) technique under magic angle spinning (MAS) conditions to demonstrate the feasibility of the method for studies of slow motions in the solid state. For the quadrupolar anisotropic interaction, the essence of CEST is to scan the saturation pattern over a range of offsets corresponding to the entire spectral region(s) for all conformational states involved, which translates into a range of −60–+ 60 kHz for methyl groups. Rotary resonances occur when the offsets are at half-and full-integer of the MAS rates. The choice of the optimal MAS rate is governed by the condition to reduce the number of rotary resonances in the CEST profile patterns and retain a sufficiently large quadrupolar interaction active under MAS to maintain sensitivity to motions. As examples, we applied this technique to a well-known model compound dimethyl-sulfone (DMS) as well as amyloid-β fibrils selectively deuterated at a single methyl group of A2 belonging to the disordered domain. It is demonstrated that the obtained exchange rate between the two rotameric states of DMS at elevated temperatures fell within known ranges and the fitted model parameters for the fibrils agree well with the previously obtained value using static ²H NMR techniques. Additionally, for the fibrils we have observed characteristic broadening of rotary resonances in the presence of conformational exchange, which provides implications for model selection and refinement. This work sets the stage for future potential extensions of the ²H CEST under MAS technique to multiple-labeled samples in small molecules and proteins.

From Selection to Instruction and Back: Competing Conformational Selection and Induced Fit Pathways in Abiotic Hosts

Two limiting cases of molecular recognition, induced fit (IF) and conformational selection (CS), play a central role in allosteric regulation of natural systems. The IF paradigm states that a substrate "instructs" the host to change its shape after complexation, while CS asserts that a guest "selects" the optimal fit from an ensemble of preexisting host conformations. With no studies that quantitatively address the interplay of two limiting pathways in abiotic systems, we herein and for the first time describe the way by which twisted capsule M-1, encompassing two conformers M-1(+) and M-1(-), trap CX₄ (X=Cl, Br) to give CX₄ ⊂M-1(+) and CX₄ ⊂M-1(-), with all four states being in thermal equilibrium. With the assistance of 2D EXSY, we found that CBr₄ would, at its lower concentrations, bind M-1 via a M-1(+)→M-1(-)→CBr₄ ⊂M-1(-) pathway corresponding to conformational selection. For M-1 complexing CCl₄ though, data from 2D EXSY measurements and 1D NMR line-shape analysis suggested that lower CCl₄ concentrations would favor CS while the IF pathway prevailed at higher proportions of the guest. Since CS and IF are not mutually exclusive, we reason that our work sets the stage for characterizing the dynamics of a wide range of already existing hosts to broaden our fundamental understanding of their action. The objective is to master the way in which encapsulation takes place for designing novel and allosteric sequestering agents, catalysts and chemosensors akin to those found in nature.

Allosteric and dynamic control of RNA-dependent RNA polymerase function and fidelity

All RNA viruses encode an RNA-dependent RNA polymerase (RdRp) responsible for genome replication. It is now recognized that enzymes in general, and RdRps specifically, are dynamic macromolecular machines such that their moving parts, including active site loops, play direct functional roles. While X-ray crystallography has provided deep insight into structural elements important for RdRp function, this methodology generally provides only static snapshots, and so is limited in its ability to report on dynamic fluctuations away from the lowest energy conformation. Nuclear magnetic resonance (NMR), molecular dynamics (MD) simulations and other biophysical techniques have brought new insight into RdRp function by their ability to characterize the trajectories, kinetics and thermodynamics of conformational motions. In particular, these methodologies have identified coordinated motions among conserved structural motifs necessary for nucleotide selection and incorporation. Disruption of these motions through amino acid substitutions or inhibitor binding impairs RdRp function. Understanding and re-engineering these motions thus provides exciting new avenues for anti-viral strategies. This chapter outlines the basics of these methodologies, summarizes the dynamic motions observed in different RdRps important for nucleotide selection and incorporation, and illustrates how this information can be leveraged towards rational vaccine strain development and anti-viral drug design.

Hydrogen/Deuterium Exchange and Nuclear Magnetic Resonance Spectroscopy Reveal Dynamic Allostery on Multiple Time Scales in the Serine Protease Thrombin

A deeper understanding of how hydrogen/deuterium exchange mass spectrometry (HDX-MS) reveals allostery is important because HDX-MS can reveal allostery in systems that are not amenable to nuclear magnetic resonance (NMR) spectroscopy. We were able to study thrombin and its complex with thrombomodulin, an allosteric regulator, by both HDX-MS and NMR. In this Perspective, we compare and contrast the results from both experiments and from molecular dynamics simulations. NMR detects changes in the chemical environment around the protein backbone N-H bond vectors, providing residue-level information about the conformational exchange between distinct states. HDX-MS detects changes in amide proton solvent accessibility and H-bonding. Taking advantage of NMR relaxation dispersion measurements of the time scale of motions, we draw conclusions about the motions reflected in HDX-MS experiments. Both experiments detect allostery, but they reveal different components of the allosteric transition. The insights gained from integrating NMR and HDX-MS into thrombin dynamics enable a clearer interpretation of the evidence for allostery revealed by HDX-MS in larger protein complexes and assemblies that are not amenable to NMR.

Dimerization of Cadherin-11 involves multi-site coupled unfolding and strand swapping

Cadherin extracellular domain 1 (EC1) mediates homophilic dimerization in adherens junctions. Conserved Trp2 and Trp4 residues in type II cadherins anchor the EC1 A strand intermolecularly in strand-swapped dimers. Herein, NMR spectroscopy is used to elucidate the roles of Trp2 and Trp4 in Cadherin-11 dimerization. The monomeric state, with the A strand and Trp side chains packed intramolecularly, is in equilibrium with sparsely populated partially and fully A-strand-exposed states, in which Trp2 (and Trp4, respectively) side-chain packing is disrupted. Exchange kinetics between the major state and the partially (fully) A-strand-exposed state is slow-intermediate (intermediate-fast). A separate very fast process exchanges ordered and random-coil BC-loop conformations with populations dependent on A-strand exposure and dimerization status. In addition, very slow processes connect the folded A-strand-exposed conformation to partially unfolded states, which may represent additional domain-swapping intermediates. The dimerization mechanism of type II cadherins is revealed as coupled folding and strand swapping.

Recent developments in deuterium solid-state NMR for the detection of slow motions in proteins

Slow timescale dynamics in proteins are essential for a variety of biological functions spanning ligand binding, enzymatic catalysis, protein folding and misfolding regulations, as well as protein-protein and protein-nucleic acid interactions. In this review, we focus on the experimental and theoretical developments of 2H static NMR methods applicable for studies of microsecond to millisecond motional modes in proteins, particularly rotating frame relaxation dispersion (R1ρ), quadrupolar Carr-Purcell-Meiboom-Gill (QCPMG) relaxation dispersion, and quadrupolar chemical exchange saturation transfer NMR experiments (Q-CEST). With applications chosen from amyloid-β fibrils, we show the complementarity of these approaches for elucidating the complexities of conformational ensembles in disordered domains in the non-crystalline solid state, with the employment of selective deuterium labels. Combined with recent advances in relaxation dispersion backbone measurements for 15N/13C/1H nuclei, these techniques provide powerful tools for studies of biologically relevant timescale dynamics in disordered domains in the solid state.

NMR spectroscopy captures the essential role of dynamics in regulating biomolecular function

Biomolecules are in constant motion. To understand how they function, and why malfunctions can cause disease, it is necessary to describe their three-dimensional structures in terms of dynamic conformational ensembles. Here, we demonstrate how nuclear magnetic resonance (NMR) spectroscopy provides an essential, dynamic view of structural biology that captures biomolecular motions at atomic resolution. We focus on examples that emphasize the diversity of biomolecules and biochemical applications that are amenable to NMR, such as elucidating functional dynamics in large molecular machines, characterizing transient conformations implicated in the onset of disease, and obtaining atomic-level descriptions of intrinsically disordered regions that make weak interactions involved in liquid-liquid phase separation. Finally, we discuss the pivotal role that NMR has played in driving forward our understanding of the biomolecular dynamics-function paradigm.

The Role of Conformational Dynamics in the Recognition and Regulation of Ubiquitination

The ubiquitination pathway is central to many cell signaling and regulatory events. One of the intriguing aspects of the pathway is the combinatorial sophistication of substrate recognition and ubiquitin chain building determinations. The abundant structural and biological data portray several characteristic protein folds among E2 and E3 proteins, and the understanding of the combinatorial complexity that enables interaction with much of the human proteome is a major goal to developing targeted and selective manipulation of the pathway. With the commonality of some folds, there are likely other aspects that can provide differentiation and recognition. These aspects involve allosteric effects and conformational dynamics that can direct recognition and chain building processes. In this review, we will describe the current state of the knowledge for conformational dynamics across a wide timescale, address the limitations of present approaches, and illustrate the potential to make new advances in connecting dynamics with ubiquitination regulation.

Algebraic expressions for Carr-Purcell-Meiboom-Gill relaxation dispersion for N-site chemical exchange

The Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiment measures the effective relaxation rate constant during a train of spin-echo pulse sequence elements as a function of the echo time. The CPMG experiment is a powerful method for characterizing chemical and conformational dynamic processes, termed chemical and conformational exchange, on μs-ms time scales, comparable to the experimentally accessible echo times. Approximate theoretical expressions for the effective relaxation rate constant for N-site chemical exchange have been reported (H. Koss, M. Rance, and A. G. Palmer, Biochemistry 57, 4753-4763 (2018)). Expressions for the effective relaxation rate constant have been improved by using the Cayley-Hamilton theorem to obtain simple and accurate approximations of the average Liouvillian for the CPMG experiment. The improved accuracy of the results allows efficient analyses of experimental data. In addition, the relationship is clarified between the approach of Koss and coworkers and that of Jen (1978).

Impact of 5-formylcytosine on the melting kinetics of DNA by 1H NMR chemical exchange

5-Formylcytosine (5fC) is a chemically edited, naturally occurring nucleobase which appears in the context of modified DNA strands. The understanding of the impact of 5fC on dsDNA physical properties is to date limited. In this work, we applied temperature-dependent 1H Chemical Exchange Saturation Transfer (CEST) NMR experiments to non-invasively and site-specifically measure the thermodynamic and kinetic influence of formylated cytosine nucleobase on the melting process involving dsDNA. Incorporation of 5fC within symmetrically positioned CpG sites destabilizes the whole dsDNA structure-as witnessed from the ∼2°C decrease in the melting temperature and 5-10 kJ mol-1 decrease in ΔG°-and affects the kinetic rates of association and dissociation. We observed an up to ∼5-fold enhancement of the dsDNA dissociation and an up to ∼3-fold reduction in ssDNA association rate constants, over multiple temperatures and for several proton reporters. Eyring and van't Hoff analysis proved that the destabilization is not localized, instead all base-pairs are affected and the transition states resembles the single-stranded conformation. These results advance our knowledge about the role of 5fC as a semi-permanent epigenetic modification and assist in the understanding of its interactions with reader proteins.

A CEST NMR experiment to obtain glycine 1Hα chemical shifts in 'invisible' minor states of proteins

Chemical exchange saturation transfer (CEST) experiments are routinely used to study protein conformational exchange between a 'visible' major state and 'invisible' minor states because they can detect minor states with lifetimes varying from ~ 3 to ~ 100 ms populated to just ~ 0.5%. Consequently several ¹H, ¹⁵N and ¹³C CEST experiments have been developed to study exchange and obtain minor state chemical shifts at almost all backbone and sidechain sites in proteins. Conspicuously missing from this extensive set of CEST experiments is a ¹H CEST experiment to study exchange at glycine (Gly) ¹H^α sites as the existing ¹H CEST experiments that have been designed to study dynamics in amide ¹H-¹⁵N spin systems and methyl ¹³CH₃ groups with three equivalent protons while suppressing ¹H-¹H NOE induced dips are not suitable for studying exchange in methylene ¹³CH₂ groups with inequivalent protons. Here a Gly ¹H^α CEST experiment to obtain the minor state Gly ¹H^α chemical shifts is presented. The utility of this experiment is demonstrated on the L99A cavity mutant of T4 Lysozyme (T4L L99A) that undergoes conformational exchange between two compact conformers. The CEST derived minor state Gly ¹H^α chemical shifts of T4L L99A are in agreement with those obtained previously using CPMG techniques. The experimental strategy presented here can also be used to obtain methylene proton minor state chemical shifts from protein sidechain and nucleic acid backbone sites.

Achieving pure spin effects by artifact suppression in methyl adiabatic relaxation experiments

Recent methyl adiabatic relaxation dispersion experiments provide examination of conformational dynamics across a very wide timescale (102-105 s-1) and, particularly, provide insight into the hydrophobic core of proteins and allosteric effects associated with modulators. The experiments require efficient decoupling of 1H and 13C spin interactions, and some artifacts have been discovered, which are associated with the design of the proton decoupling scheme. The experimental data suggest that the original design is valid; however, pulse sequences with either no proton decoupling or proton decoupling with imperfect pulses can potentially exhibit complications in the experiments. Here, we demonstrate that pulse imperfections in the proton decoupling scheme can be dramatically alleviated by using a single composite π pulse and provide pure single-exponential relaxation data. It allows the opportunity to access high-quality methyl adiabatic relaxation dispersion data by removing the cross-correlation between dipole-dipole interaction and chemical shift anisotropy. The resulting high-quality data is illustrated with the binding of an allosteric modulator (G2BR) to the ubiquitin conjugating enzyme Ube2g2.

Compact expressions for R1ρ relaxation for N-site chemical exchange using Schur decomposition

The rotating-frame spin relaxation rate constant, R_1ρ, is a powerful probe of macromolecular chemical and conformational dynamics in relaxation dispersion, CEST, and DEST NMR experiments. The R_1ρ relaxation rate constant is given by the absolute value of the largest (least negative) eigenvalue of the Bloch-McConnell evolution matrix; however, estimation of this eigenvalue require inversion of 3 N × 3 N dimensional matrices, in which N is the number of interconverting sites or states for a given nuclear spin in a molecule. The Schur complement is used to reduce the problem of calculating the characteristic polynomial of a 3 N × 3 N matrix to that of calculating the characteristic polynomial of a 3 × 3 matrix. The resulting expressions for N-site chemical exchange are more numerically tractable, because the largest matrix inversion also is of dimension 3 × 3. In addition, the simplifications offered by the Schurr complement conveniently illustrate the effects of fast or slow kinetic steps within an N-site kinetic topology.

Deuteron Quadrupolar Chemical Exchange Saturation Transfer (Q-CEST) Solid-State NMR for Static Powder Samples: Approach and Applications to Amyloid-β Fibrils

We provide an experimental and computational framework for 2 H quadrupolar chemical exchange saturation transfer NMR experiments (Q-CEST) under static solid-state conditions for the quantification of dynamics on μs-ms timescales. Simulations using simple 2-site exchange models provide insights into the relation between spin dynamics and motions. Biological applications focus on two sites of amyloid-β fibrils in the 3-fold symmetric polymorph. The first site, the methyl group of A2 of the disordered N-terminal domain, undergoes diffusive motions and conformational exchange due to transient interactions. Earlier 2 H rotating frame relaxation and quadrupolar CPMG measurements are combined with the Q-CEST approach to characterize the multiple conformational states of the domain. The second site, the methyl group of M35, spans the water-accessible cavity inside the fibrils' core and undergoes extensive rotameric exchange. Q-CEST permits us to refine the rotameric exchange model for this site and allows the more precise determination of populations and rotameric exchange rate constants than line shape analysis.

Methyl TROSY spectroscopy: A versatile NMR approach to study challenging biological systems

A major goal in structural biology is to unravel how molecular machines function in detail. To that end, solution-state NMR spectroscopy is ideally suited as it is able to study biological assemblies in a near natural environment. Based on methyl TROSY methods, it is now possible to record high-quality data on complexes that are far over 100 kDa in molecular weight. In this review, we discuss the theoretical background of methyl TROSY spectroscopy, the information that can be extracted from methyl TROSY spectra and approaches that can be used to assign methyl resonances in large complexes. In addition, we touch upon insights that have been obtained for a number of challenging biological systems, including the 20S proteasome, the RNA exosome, molecular chaperones and G-protein-coupled receptors. We anticipate that methyl TROSY methods will be increasingly important in modern structural biology approaches, where information regarding static structures is complemented with insights into conformational changes and dynamic intermolecular interactions.