Chemical exchange in biomacromolecules: past, present, and future

An optimal control approach to design entire relaxation dispersion experiments

A general approach is introduced to optimize experiments for the analysis of spin systems in the presence of chemical exchange. Rather than optimizing individual pulse sequence elements, such as refocusing pulses, entire relaxation dispersion sequences are optimized in the form of a single shaped pulse. This is achieved by defining a performance index that is only based on the remaining signal after the relaxation dispersion sequence for a range of exchange, relaxation, offset, and rf inhomogeneity parameters. The approach is demonstrated by optimizing energy-limited broadband relaxation dispersion sequences that closely approach the overall effect of ideal CPMG sequences. As illustrated both theoretically and experimentally, significant improvements are found compared to standard amplitude or energy-limited CPMG sequences.

An excited state underlies gene regulation of a transcriptional riboswitch

Riboswitches control gene expression through ligand-dependent structural rearrangements of the sensing aptamer domain. However, we found that the Bacillus cereus fluoride riboswitch aptamer adopts identical tertiary structures in solution with and without ligand. Using chemical-exchange saturation transfer (CEST) NMR spectroscopy, we revealed that the structured ligand-free aptamer transiently accesses a low-populated (∼1%) and short-lived (∼3 ms) excited conformational state that unravels a conserved 'linchpin' base pair to signal transcription termination. Upon fluoride binding, this highly localized, fleeting process is allosterically suppressed, which activates transcription. We demonstrated that this mechanism confers effective fluoride-dependent gene activation over a wide range of transcription rates, which is essential for robust toxicity responses across diverse cellular conditions. These results unveil a novel switching mechanism that employs ligand-dependent suppression of an aptamer excited state to coordinate regulatory conformational transitions rather than adopting distinct aptamer ground-state tertiary architectures, exemplifying a new mode of ligand-dependent RNA regulation.

Second-order quadrupolar line shapes under molecular dynamics: An additional transition in the extremely fast regime

NMR spectroscopy is a powerful tool for probing molecular dynamics. For the classic case of two-site exchange, NMR spectra go through the transition from exchange broadening through coalescence and then motional narrowing as the exchange rate increases passing through the difference between the resonance frequencies of the two sites. For central-transition spectra of half-integer quadrupolar nuclei in solids, line shape change due to molecular dynamics occurs in two stages. The first stage occurs when the exchange rate is comparable to the second-order quadrupolar interaction. The second spectral transition comes at a faster exchange rate which approaches the Larmor frequency and generally reduces the isotropic quadrupolar shift. Such a two-stage transition phenomenon is unique to half-integer quadrupolar nuclei. A quantum mechanical formalism in full Liouville space is presented to explain the physical origin of the two-stage phenomenon and for use in spectral simulations. Variable-temperature ¹⁷O NMR of solid NaNO₃ in which the NO₃^- ion undergoes 3-fold jumps confirms the two-stage transition process. The spectra of NaNO₃ acquired in the temperature range of 173-413K agree well with simulations using the quantum mechanical formalism. The rate constants for the 3-fold NO₃^- ion jumps span eight orders of magnitude (10²-10¹⁰s^-1) covering both transitions of the dynamic ¹⁷O line shape.

Recent advances in measuring the kinetics of biomolecules by NMR relaxation dispersion spectroscopy

Protein function can be modulated or dictated by the amplitude and timescale of biomolecular motion, therefore it is imperative to study protein dynamics. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful technique capable of studying timescales of motion that range from those faster than molecular reorientation on the picosecond timescale to those that occur in real-time. Across this entire regime, NMR observables can report on the amplitude of atomic motion, and the kinetics of atomic motion can be ascertained with a wide variety of experimental techniques from real-time to milliseconds and several nanoseconds to picoseconds. Still a four orders of magnitude window between several nanoseconds and tens of microseconds has remained elusive. Here, we highlight new relaxation dispersion NMR techniques that serve to cover this "hidden-time" window up to hundreds of nanoseconds that achieve atomic resolution while studying the molecule under physiological conditions.

Exploration of the Conformational Dynamics of Major Histocompatibility Complex Molecules

Major histocompatibility complex (MHC) molecules are loaded with a wide variety of self- and non-self-peptides in their binding grooves and present these to T cell receptors (TCRs) in order to activate the adaptive immune system. A large number of crystal structures of different MHC alleles with different bound peptides have been determined, and they have been found to be quite similar to one another regardless of the bound peptide sequence. The structures do not change markedly even when forming complexes with TCRs. Nonetheless, the degree of TCR activation does differ markedly depending on the peptide presented by the MHC. Recent structural studies in solution rather than as crystals have suggested that the conformational dynamics of MHC molecules may be responsible for the MHC stability differences. Furthermore, it was shown that the conformational dynamics of MHC molecules is important for peptide loading and presentation to TCR. Here, we describe the static and dynamic structures of MHC molecules and appropriate methods to analyze them. We focus particularly on nuclear magnetic resonance (NMR), one of the most powerful tools to study dynamic properties of proteins. The number of such studies in the literature is limited, but in this review, we show that NMR is valuable for elucidating the structural dynamics of MHC molecules.

Insights into Watson-Crick/Hoogsteen breathing dynamics and damage repair from the solution structure and dynamic ensemble of DNA duplexes containing m1A

In the canonical DNA double helix, Watson-Crick (WC) base pairs (bps) exist in dynamic equilibrium with sparsely populated (∼0.02-0.4%) and short-lived (lifetimes ∼0.2-2.5 ms) Hoogsteen (HG) bps. To gain insights into transient HG bps, we used solution-state nuclear magnetic resonance spectroscopy, including measurements of residual dipolar couplings and molecular dynamics simulations, to examine how a single HG bp trapped using the N1-methylated adenine (m1A) lesion affects the structural and dynamic properties of two duplexes. The solution structure and dynamic ensembles of the duplexes reveals that in both cases, m1A forms a m1A•T HG bp, which is accompanied by local and global structural and dynamic perturbations in the double helix. These include a bias toward the BI backbone conformation; sugar repuckering, major-groove directed kinking (∼9°); and local melting of neighboring WC bps. These results provide atomic insights into WC/HG breathing dynamics in unmodified DNA duplexes as well as identify structural and dynamic signatures that could play roles in m1A recognition and repair.

Applications of NMR and computational methodologies to study protein dynamics

Overwhelming evidence now illustrates the defining role of atomic-scale protein flexibility in biological events such as allostery, cell signaling, and enzyme catalysis. Over the years, spin relaxation nuclear magnetic resonance (NMR) has provided significant insights on the structural motions occurring on multiple time frames over the course of a protein life span. The present review article aims to illustrate to the broader community how this technique continues to shape many areas of protein science and engineering, in addition to being an indispensable tool for studying atomic-scale motions and functional characterization. Continuing developments in underlying NMR technology alongside software and hardware developments for complementary computational approaches now enable methodologies to routinely provide spatial directionality and structural representations traditionally harder to achieve solely using NMR spectroscopy. In addition to its well-established role in structural elucidation, we present recent examples that illustrate the combined power of selective isotope labeling, relaxation dispersion experiments, chemical shift analyses, and computational approaches for the characterization of conformational sub-states in proteins and enzymes.

Protein folding by NMR

Protein folding is a highly complex process proceeding through a number of disordered and partially folded nonnative states with various degrees of structural organization. These transiently and sparsely populated species on the protein folding energy landscape play crucial roles in driving folding toward the native conformation, yet some of these nonnative states may also serve as precursors for protein misfolding and aggregation associated with a range of devastating diseases, including neuro-degeneration, diabetes and cancer. Therefore, in vivo protein folding is often reshaped co- and post-translationally through interactions with the ribosome, molecular chaperones and/or other cellular components. Owing to developments in instrumentation and methodology, solution NMR spectroscopy has emerged as the central experimental approach for the detailed characterization of the complex protein folding processes in vitro and in vivo. NMR relaxation dispersion and saturation transfer methods provide the means for a detailed characterization of protein folding kinetics and thermodynamics under native-like conditions, as well as modeling high-resolution structures of weakly populated short-lived conformational states on the protein folding energy landscape. Continuing development of isotope labeling strategies and NMR methods to probe high molecular weight protein assemblies, along with advances of in-cell NMR, have recently allowed protein folding to be studied in the context of ribosome-nascent chain complexes and molecular chaperones, and even inside living cells. Here we review solution NMR approaches to investigate the protein folding energy landscape, and discuss selected applications of NMR methodology to studying protein folding in vitro and in vivo. Together, these examples highlight a vast potential of solution NMR in providing atomistic insights into molecular mechanisms of protein folding and homeostasis in health and disease.

Protein NMR Studies of Substrate Binding to Human Blood Group A and B Glycosyltransferases

Donor and acceptor substrate binding to human blood group A and B glycosyltransferases (GTA, GTB) has been studied by a variety of protein NMR experiments. Prior crystallographic studies had shown these enzymes to adopt an open conformation in the absence of substrates. Binding either of the donor substrate UDP-Gal or of UDP induces a semiclosed conformation. In the presence of both donor and acceptor substrates, the enzymes shift towards a closed conformation with ordering of an internal loop and the C-terminal residues, which then completely cover the donor-binding pocket. Chemical-shift titrations of uniformly ² H,¹⁵ N-labeled GTA or GTB with UDP affected about 20 % of all crosspeaks in ¹ H,¹⁵ N TROSY-HSQC spectra, reflecting substantial plasticity of the enzymes. On the other hand, it is this conformational flexibility that impedes NH backbone assignments. Chemical-shift-perturbation experiments with δ1-[¹³ C]methyl-Ile-labeled samples revealed two Ile residues-Ile123 at the bottom of the UDP binding pocket, and Ile192 as part of the internal loop-that were significantly disturbed upon stepwise addition of UDP and H-disaccharide, also revealing long-range perturbations. Finally, methyl TROSY-based relaxation dispersion experiments do not reveal micro- to millisecond timescale motions. Although this study reveals substantial conformational plasticity of GTA and GTB, the matter of how binding of substrates shifts the enzymes into catalytically competent states remains enigmatic.

Molecular Cross-Talk between Nonribosomal Peptide Synthetase Carrier Proteins and Unstructured Linker Regions

Nonribosomal peptide synthetases (NRPSs) employ multiple domains separated by linker regions to incorporate substrates into natural products. During synthesis, substrates are covalently tethered to carrier proteins that translocate between catalytic partner domains. The molecular parameters that govern translocation and associated linker remodeling remain unknown. Here, we used NMR to characterize the structure, dynamics, and invisible states of a peptidyl carrier protein flanked by its linkers. We showed that the N-terminal linker stabilizes and interacts with the protein core while modulating dynamics at specific sites involved in post-translational modifications and/or domain interactions. The results detail the molecular communication between peptidyl carrier proteins and their linkers and could guide efforts in engineering NRPSs to obtain new pharmaceuticals.

Probing conformational dynamics in biomolecules via chemical exchange saturation transfer: a primer

Although Chemical Exchange Saturation Transfer (CEST) type NMR experiments have been used to study chemical exchange processes in molecules since the early 1960s, there has been renewed interest in the past several years in using this approach to study biomolecular conformational dynamics. The methodology is particularly powerful for the study of sparsely populated, transiently formed conformers that are recalcitrant to investigation using traditional biophysical tools, and it is complementary to relaxation dispersion and magnetization transfer experiments that have traditionally been used to study chemical exchange processes. Here we discuss the concepts behind the CEST experiment, focusing on practical aspects as well, we review some of the pulse sequences that have been developed to characterize protein and RNA conformational dynamics, and we discuss a number of examples where the CEST methodology has provided important insights into the role of dynamics in biomolecular function.

Confinement and Stabilization of Fyn SH3 Folding Intermediate Mimetics within the Cavity of the Chaperonin GroEL Demonstrated by Relaxation-Based NMR

The interaction of two folding intermediate mimetics of the model protein substrate Fyn SH3 with the chaperonin GroEL, a supramolecular foldase/unfoldase machine, has been investigated by ¹⁵N relaxation-based nuclear magnetic resonance spectroscopy (lifetime line broadening, dark state exchange saturation transfer, and relaxation dispersion). The two mimetics comprise C-terminal truncations of wild-type and triple-mutant (A39V/N53P/V55L) Fyn SH3 in which the C-terminal strand of the SH3 domain is unfolded, while preserving the remaining domain structure. Quantitative analysis of the data reveals that a mobile state of the SH3 domain confined and tethered within the cavity of GroEL, possibly through interactions with the disordered, methionine-rich C-terminal tail(s), can be detected, and that the native state of the folding intermediate mimetics is stabilized by both confinement within and binding to apo GroEL. These data provide a basis for understanding the passive activity of GroEL as a foldase/unfoldase: the unfolded state, in the absence of hydrophobic GroEL-binding consensus sequences, is destabilized within the cavity because of its larger radius of gyration compared to that of the folding intermediate, while the folding intermediate is stabilized relative to the native state because of exposure of a hydrophobic patch that favors GroEL binding.

Evaluation of the Coupled Two-Dimensional Main Chain Torsional Potential in Modeling Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) carry out crucial biological functions in essential biological processes of life. Because of the highly dynamic and conformationally heterogeneous nature of the disordered states of IDPs, molecular dynamics simulations are becoming an indispensable tool for the investigation of the conformational ensembles and dynamic properties of IDPs. Nevertheless, there is still no consensus on the most reliable force field in molecular dynamics simulations for IDPs hitherto. In this work, the recently proposed AMBER99SB^2D force field is evaluated in modeling some disordered polypeptides and proteins by checking its ability to reproduce experimental NMR data. The results highlight that when the ildn side-chain corrections are included, AMBER99SB^2D-ildn exhibits reliable results that agree with experiments compared with its predecessors, the AMBER14SB, AMBER99SB, AMBER99SB-ildn, and AMBER99SB^2D force fields, and that decreasing the overall magnitude of protein-protein interactions in favor of protein-water interactions is a key ingredient behind the improvement.

General expressions for R1ρ relaxation for N-site chemical exchange and the special case of linear chains

Exploration of dynamic processes in proteins and nucleic acids by spin-locking NMR experiments has been facilitated by the development of theoretical expressions for the R_1ρ relaxation rate constant covering a variety of kinetic situations. Herein, we present a generalized approximation to the chemical exchange, R_ex, component of R_1ρ for arbitrary kinetic schemes, assuming the presence of a dominant major site population, derived from the negative reciprocal trace of the inverse Bloch-McConnell evolution matrix. This approximation is equivalent to first-order truncation of the characteristic polynomial derived from the Bloch-McConnell evolution matrix. For three- and four-site chemical exchange, the first-order approximations are sufficient to distinguish different kinetic schemes. We also introduce an approach to calculate R_1ρ for linear N-site schemes, using the matrix determinant lemma to reduce the corresponding 3N×3N Bloch-McConnell evolution matrix to a 3×3 matrix. The first- and second order-expansions of the determinant of this 3×3 matrix are closely related to previously derived equations for two-site exchange. The second-order approximations for linear N-site schemes can be used to obtain more accurate approximations for non-linear N-site schemes, such as triangular three-site or star four-site topologies. The expressions presented herein provide powerful means for the estimation of R_ex contributions for both low (CEST-limit) and high (R_1ρ-limit) radiofrequency field strengths, provided that the population of one state is dominant. The general nature of the new expressions allows for consideration of complex kinetic situations in the analysis of NMR spin relaxation data.

Excited States of Nucleic Acids Probed by Proton Relaxation Dispersion NMR Spectroscopy

In this work an improved stable isotope labeling protocol for nucleic acids is introduced. The novel building blocks eliminate/minimize homonuclear (13) C and (1) H scalar couplings thus allowing proton relaxation dispersion (RD) experiments to report accurately on the chemical exchange of nucleic acids. Using site-specific (2) H and (13) C labeling, spin topologies are introduced into DNA and RNA that make (1) H relaxation dispersion experiments applicable in a straightforward manner. The novel RNA/DNA building blocks were successfully incorporated into two nucleic acids. The A-site RNA was previously shown to undergo a two site exchange process in the micro- to millisecond time regime. Using proton relaxation dispersion experiments the exchange parameters determined earlier could be recapitulated, thus validating the proposed approach. We further investigated the dynamics of the cTAR DNA, a DNA transcript that is involved in the viral replication cycle of HIV-1. Again, an exchange process could be characterized and quantified. This shows the general applicablility of the novel labeling scheme for (1) H RD experiments of nucleic acids.

The Dynamic Basis for Signal Propagation in Human Pin1-WW

Allostery is the structural manifestation of information transduction in biomolecules. Its hallmark is conformational change induced by perturbations at a distal site. An increasing body of evidence demonstrates the presence of allostery in very flexible and even disordered proteins, encouraging a thermodynamic description of this phenomenon. Still, resolving such processes at atomic resolution is difficult. Here we establish a protocol to determine atomistic thermodynamic models of such systems using high-resolution solution state nuclear magnetic resonance data and extensive molecular simulations. Using this methodology, we study information transduction in the WW domain of a key cell-cycle regulator Pin1. Pin1 binds promiscuously to phospho-Ser/Thr-Pro motifs, however, disparate structural and dynamic responses have been reported upon binding different ligands. Our model consists of two topologically distinct states whose relative population may be specifically skewed by an incoming ligand. This model provides a canonical basis for the understanding of multi-functionality in Pin1.

Studying Dynamics by Magic-Angle Spinning Solid-State NMR Spectroscopy: Principles and Applications to Biomolecules

Magic-angle spinning solid-state NMR spectroscopy is an important technique to study molecular structure, dynamics and interactions, and is rapidly gaining importance in biomolecular sciences. Here we provide an overview of experimental approaches to study molecular dynamics by MAS solid-state NMR, with an emphasis on the underlying theoretical concepts and differences of MAS solid-state NMR compared to solution-state NMR. The theoretical foundations of nuclear spin relaxation are revisited, focusing on the particularities of spin relaxation in solid samples under magic-angle spinning. We discuss the range of validity of Redfield theory, as well as the inherent multi-exponential behavior of relaxation in solids. Experimental challenges for measuring relaxation parameters in MAS solid-state NMR and a few recently proposed relaxation approaches are discussed, which provide information about time scales and amplitudes of motions ranging from picoseconds to milliseconds. We also discuss the theoretical basis and experimental measurements of anisotropic interactions (chemical-shift anisotropies, dipolar and quadrupolar couplings), which give direct information about the amplitude of motions. The potential of combining relaxation data with such measurements of dynamically-averaged anisotropic interactions is discussed. Although the focus of this review is on the theoretical foundations of dynamics studies rather than their application, we close by discussing a small number of recent dynamics studies, where the dynamic properties of proteins in crystals are compared to those in solution.

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

Evaluating the influence of initial magnetization conditions on extracted exchange parameters in NMR relaxation experiments: applications to CPMG and CEST

Transient excursions of native protein states to functionally relevant higher energy conformations often occur on the μs-ms timescale. NMR spectroscopy has emerged as an important tool to probe such processes using techniques such as Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion and Chemical Exchange Saturation Transfer (CEST). The extraction of kinetic and structural parameters from these measurements is predicated upon mathematical modeling of the resulting relaxation profiles, which in turn relies on knowledge of the initial magnetization conditions at the start of the CPMG/CEST relaxation elements in these experiments. Most fitting programs simply assume initial magnetization conditions that are given by equilibrium populations, which may be incorrect in certain implementations of experiments. In this study we have quantified the systematic errors in extracted parameters that are generated from analyses of CPMG and CEST experiments using incorrect initial boundary conditions. We find that the errors in exchange rates (k ex ) and populations (p E ) are typically small (<10 %) and thus can be safely ignored in most cases. However, errors become larger and cannot be fully neglected (20-40 %) as k ex falls near the lower limit of each method or when short CPMG/CEST relaxation elements are used in these experiments. The source of the errors can be rationalized and their magnitude given by a simple functional form. Despite the fact that errors tend to be small, it is recommended that the correct boundary conditions be implemented in fitting programs so as to obtain as robust estimates of exchange parameters as possible.

The molecular choreography of protein synthesis: translational control, regulation, and pathways

Translation of proteins by the ribosome regulates gene expression, with recent results underscoring the importance of translational control. Misregulation of translation underlies many diseases, including cancer and many genetic diseases. Decades of biochemical and structural studies have delineated many of the mechanistic details in prokaryotic translation, and sketched the outlines of eukaryotic translation. However, translation may not proceed linearly through a single mechanistic pathway, but likely involves multiple pathways and branchpoints. The stochastic nature of biological processes would allow different pathways to occur during translation that are biased by the interaction of the ribosome with other translation factors, with many of the steps kinetically controlled. These multiple pathways and branchpoints are potential regulatory nexus, allowing gene expression to be tuned at the translational level. As research focus shifts toward eukaryotic translation, certain themes will be echoed from studies on prokaryotic translation. This review provides a general overview of the dynamic data related to prokaryotic and eukaryotic translation, in particular recent findings with single-molecule methods, complemented by biochemical, kinetic, and structural findings. We will underscore the importance of viewing the process through the viewpoints of regulation, translational control, and heterogeneous pathways.

STARD6 on steroids: solution structure, multiple timescale backbone dynamics and ligand binding mechanism

START domain proteins are conserved α/β helix-grip fold that play a role in the non-vesicular and intracellular transport of lipids and sterols. The mechanism and conformational changes permitting the entry of the ligand into their buried binding sites is not well understood. Moreover, their functions and the identification of cognate ligands is still an active area of research. Here, we report the solution structure of STARD6 and the characterization of its backbone dynamics on multiple time-scales through ¹⁵N spin-relaxation and amide exchange studies. We reveal for the first time the presence of concerted fluctuations in the Ω₁ loop and the C-terminal helix on the microsecond-millisecond time-scale that allows for the opening of the binding site and ligand entry. We also report that STARD6 binds specifically testosterone. Our work represents a milestone for the study of ligand binding mechanism by other START domains and the elucidation of the biological function of STARD6.

Biophysical and computational methods to analyze amino acid interaction networks in proteins

Globular proteins are held together by interacting networks of amino acid residues. A number of different structural and computational methods have been developed to interrogate these amino acid networks. In this review, we describe some of these methods, including analyses of X-ray crystallographic data and structures, computer simulations, NMR data, and covariation among protein sequences, and indicate the critical insights that such methods provide into protein function. This information can be leveraged towards the design of new allosteric drugs, and the engineering of new protein function and protein regulation strategies.

Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch

Intermediates play important roles in RNA folding but can be difficult to characterize when short-lived or not significantly populated. By combining (15)N relaxation dispersion NMR with chemical probing, we visualized a fast (kex=k1+k-1≈423 s(-1)) secondary structural switch directed towards a low-populated (∼3%) partially folded intermediate in tertiary folding of the P5abc subdomain of the 'Tetrahymena' group I intron ribozyme. The secondary structure switch changes the base-pairing register across the P5c hairpin, creating a native-like structure, and occurs at rates of more than two orders of magnitude faster than tertiary folding. The switch occurs robustly in the absence of tertiary interactions, Mg(2+) or even when the hairpin is excised from the three-way junction. Fast, highly modular secondary structural switches may be quite common during RNA tertiary folding where they may help smoothen the folding landscape by allowing folding to proceed efficiently via additional pathways.

Distribution of Pico- and Nanosecond Motions in Disordered Proteins from Nuclear Spin Relaxation

Intrinsically disordered proteins and intrinsically disordered regions (IDRs) are ubiquitous in the eukaryotic proteome. The description and understanding of their conformational properties require the development of new experimental, computational, and theoretical approaches. Here, we use nuclear spin relaxation to investigate the distribution of timescales of motions in an IDR from picoseconds to nanoseconds. Nitrogen-15 relaxation rates have been measured at five magnetic fields, ranging from 9.4 to 23.5 T (400-1000 MHz for protons). This exceptional wealth of data allowed us to map the spectral density function for the motions of backbone NH pairs in the partially disordered transcription factor Engrailed at 11 different frequencies. We introduce an approach called interpretation of motions by a projection onto an array of correlation times (IMPACT), which focuses on an array of six correlation times with intervals that are equidistant on a logarithmic scale between 21 ps and 21 ns. The distribution of motions in Engrailed varies smoothly along the protein sequence and is multimodal for most residues, with a prevalence of motions around 1 ns in the IDR. We show that IMPACT often provides better quantitative agreement with experimental data than conventional model-free or extended model-free analyses with two or three correlation times. We introduce a graphical representation that offers a convenient platform for a qualitative discussion of dynamics. Even when relaxation data are only acquired at three magnetic fields that are readily accessible, the IMPACT analysis gives a satisfactory characterization of spectral density functions, thus opening the way to a broad use of this approach.

Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies

It is now common knowledge that enzymes are mobile entities relying on complex atomic-scale dynamics and coordinated conformational events for proper ligand recognition and catalysis. However, the exact role of protein dynamics in enzyme function remains either poorly understood or difficult to interpret. This mini-review intends to reconcile biophysical observations and biological significance by first describing a number of common experimental and computational methodologies employed to characterize atomic-scale residue motions on various timescales in enzymes, and second by illustrating how the knowledge of these motions can be used to describe the functional behavior of enzymes and even act upon it. Two biologically relevant examples will be highlighted, namely the HIV-1 protease and DNA polymerase β enzyme systems.

The energetics of a three-state protein folding system probed by high-pressure relaxation dispersion NMR spectroscopy

The energetic and volumetric properties of a three-state protein folding system, comprising a metastable triple mutant of the Fyn SH3 domain, have been investigated using pressure-dependent (15) N-relaxation dispersion NMR from 1 to 2500 bar. Changes in partial molar volumes (ΔV) and isothermal compressibilities (ΔκT ) between all the states along the folding pathway have been determined to reasonable accuracy. The partial volume and isothermal compressibility of the folded state are 100 mL mol(-1) and 40 μL mol(-1)  bar(-1) , respectively, higher than those of the unfolded ensemble. Of particular interest are the findings related to the energetic and volumetric properties of the on-pathway folding intermediate. While the latter is energetically close to the unfolded state, its volumetric properties are similar to those of the folded protein. The compressibility of the intermediate is larger than that of the folded state reflecting the less rigid nature of the former relative to the latter.

Bloch Equations for Proton Exchange Reactions in an Aqueous Solution

The extension of the Bloch equations for acid-base reactions in an aqueous solution is revisited. The acid-base reactions are second-order, and several reactions catalyzed by distinct catalysts may happen simultaneously. By constructing pseudo first-order reactions and assuming fast dissemination of protons from catalysts to solvent water, this extension converges to the well-known Bloch-McConnell equations for a two-site first-order exchange. Thus, explicit relationships between the parameters appearing in the reactions and the Bloch-McConnell equations are established. The dependencies of exchange rates and chemical exchange saturation transfer effects on pH were numerically and experimentally investigated for representative examples.

A dynamic look backward and forward

The 2015 Gunther Laukien Prize recognized solution NMR studies of protein dynamics and thermodynamics. This Perspective surveys aspects of the development and application of NMR spin relaxation for investigations of protein flexibility and function over multiple time scales in solution. Methods highlighted include analysis of overall rotational diffusion, theoretical descriptions of R1ρ relaxation, and molecular dynamics simulations to interpret NMR spin relaxation. Applications are illustrated for the zinc-finger domain Xfin-31, the calcium-binding proteins calbindin D9k and calmodulin, and the bZip transcription factor of GCN4.

NMR Methods to Study Dynamic Allostery

Nuclear magnetic resonance (NMR) spectroscopy provides a unique toolbox of experimental probes for studying dynamic processes on a wide range of timescales, ranging from picoseconds to milliseconds and beyond. Along with NMR hardware developments, recent methodological advancements have enabled the characterization of allosteric proteins at unprecedented detail, revealing intriguing aspects of allosteric mechanisms and increasing the proportion of the conformational ensemble that can be observed by experiment. Here, we present an overview of NMR spectroscopic methods for characterizing equilibrium fluctuations in free and bound states of allosteric proteins that have been most influential in the field. By combining NMR experimental approaches with molecular simulations, atomistic-level descriptions of the mechanisms by which allosteric phenomena take place are now within reach.

Residence Times of Molecular Complexes in Solution from NMR Data of Intermolecular Hydrogen-Bond Scalar Coupling

The residence times of molecular complexes in solution are important for understanding biomolecular functions and drug actions. We show that NMR data of intermolecular hydrogen-bond scalar couplings can yield information on the residence times of molecular complexes in solution. The molecular exchange of binding partners via the breakage and reformation of a complex causes self-decoupling of intermolecular hydrogen-bond scalar couplings, and this self-decoupling effect depends on the residence time of the complex. For protein-DNA complexes, we investigated the salt concentration dependence of intermolecular hydrogen-bond scalar couplings between the protein side-chain (15)N and DNA phosphate (31)P nuclei, from which the residence times were analyzed. The results were consistent with those obtained by (15)Nz-exchange spectroscopy. This self-decoupling-based kinetic analysis is unique in that it does not require any different signatures for the states involved in the exchange, whereas such conditions are crucial for kinetic analyses by typical NMR and other methods.

Protein stabilizer, NDSB-195, enhances the dynamics of the β4 -α2 loop of ubiquitin

Non-detergent sulfobetaines (NDSBs) are a new group of small, synthetic protein stabilizers, which have advantages over classical compatible osmolytes, such as polyol, amines, and amino acids: they do not increase solution viscosity, unlike polyols, and they are zwitterionic at all pH ranges, unlike amines and amino acids. NDSBs also facilitate the crystallization and refolding of proteins. The mechanism whereby NDSBs exhibit such activities, however, remains elusive. To gain insight into this mechanism, we studied, using nuclear magnetic resonance (NMR), the effects of dimethylethylammonium propane sulfonate (NDSB-195) on the dynamics of ubiquitin, on which a wealth of information has been accumulated. By analyzing the line width of amide proton resonances and the transverse relaxation rates of nitrogen atoms, we found that NDSB-195 enhances the microsecond-millisecond dynamics of a β4 -α2 loop of ubiquitin. Although those compounds that enhance protein dynamics are generally considered to destabilize protein molecules, NDSB-195 enhanced the stability of ubiquitin against guanidium chloride denaturation. Thus, the simultaneous enhancement of stability and flexibility by a single compound can be attained.

Conformational Selection in a Protein-Protein Interaction Revealed by Dynamic Pathway Analysis

Molecular recognition plays a central role in biology, and protein dynamics has been acknowledged to be important in this process. However, it is highly debated whether conformational changes happen before ligand binding to produce a binding-competent state (conformational selection) or are caused in response to ligand binding (induced fit). Proposals for both mechanisms in protein/protein recognition have been primarily based on structural arguments. However, the distinction between them is a question of the probabilities of going via these two opposing pathways. Here, we present a direct demonstration of exclusive conformational selection in protein/protein recognition by measuring the flux for rhodopsin kinase binding to its regulator recoverin, an important molecular recognition in the vision system. Using nuclear magnetic resonance (NMR) spectroscopy, stopped-flow kinetics, and isothermal titration calorimetry, we show that recoverin populates a minor conformation in solution that exposes a hydrophobic binding pocket responsible for binding rhodopsin kinase. Protein dynamics in free recoverin limits the overall rate of binding.

Speeding-up exchange-mediated saturation transfer experiments by Fourier transform

Protein motions over various time scales are crucial for protein function. NMR relaxation dispersion experiments play a key role in explaining these motions. However, the study of slow conformational changes with lowly populated states remained elusive. The recently developed exchange-mediated saturation transfer experiments allow the detection and characterization of such motions, but require extensive measurement time. Here we show that, by making use of Fourier transform, the total acquisition time required to measure an exchange-mediated saturation transfer profile can be reduced by twofold in case that one applies linear prediction. In addition, we demonstrate that the analytical solution for R1ρ experiments can be used for fitting the exchange-mediated saturation transfer profile. Furthermore, we show that simultaneous analysis of exchange-mediated saturation transfer profiles with two different radio-frequency field strengths is required for accurate and precise characterization of the exchange process and the exchanging states.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

Theoretical tools for the design of NMR relaxation dispersion pulse sequences

Recent decades have witnessed tremendous progress in the development of new experimental methods for studying biomolecules, particularly in the field of NMR relaxation dispersion. Here we review the theoretical frameworks that provided the insights necessary for such progress. The effect of radio-frequency manipulations on spin systems is discussed using Average Hamiltonian Theory (AHT), Average Liouvillian Theory (ALT), and Bloch-Wangsness-Redfield (BWR) relaxation theory. We illustrate these concepts using the case of Heteronuclear Double Resonance (HDR) methods.

Cofactor-Mediated Conformational Dynamics Promote Product Release From Escherichia coli Dihydrofolate Reductase via an Allosteric Pathway

The enzyme dihydrofolate reductase (DHFR, E) from Escherichia coli is a paradigm for the role of protein dynamics in enzyme catalysis. Previous studies have shown that the enzyme progresses through the kinetic cycle by modulating the dynamic conformational landscape in the presence of substrate dihydrofolate (DHF), product tetrahydrofolate (THF), and cofactor (NADPH or NADP(+)). This study focuses on the quantitative description of the relationship between protein fluctuations and product release, the rate-limiting step of DHFR catalysis. NMR relaxation dispersion measurements of millisecond time scale motions for the E:THF:NADP(+) and E:THF:NADPH complexes of wild-type and the Leu28Phe (L28F) point mutant reveal conformational exchange between an occluded ground state and a low population of a closed state. The backbone structures of the occluded ground states of the wild-type and mutant proteins are very similar, but the rates of exchange with the closed excited states are very different. Integrated analysis of relaxation dispersion data and THF dissociation rates measured by stopped-flow spectroscopy shows that product release can occur by two pathways. The intrinsic pathway consists of spontaneous product dissociation and occurs for all THF-bound complexes of DHFR. The allosteric pathway features cofactor-assisted product release from the closed excited state and is utilized only in the E:THF:NADPH complexes. The L28F mutation alters the partitioning between the pathways and results in increased flux through the intrinsic pathway relative to the wild-type enzyme. This repartitioning could represent a general mechanism to explain changes in product release rates in other E. coli DHFR mutants.

Dynamic Equilibria of Short-Range Electrostatic Interactions at Molecular Interfaces of Protein-DNA Complexes

Intermolecular ion pairs (salt bridges) are crucial for protein-DNA association. For two protein-DNA complexes, we demonstrate that the ion pairs of protein side-chain NH3+ and DNA phosphate groups undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. While the crystal structures of the complexes show only the solvent-separated ion pair (SIP) state for some interfacial lysine side chains, our NMR hydrogen-bond scalar coupling data clearly indicate the presence of the contact ion pair (CIP) state for the same residues. The 0.6-μs molecular dynamics (MD) simulations confirm dynamic transitions between the CIP and SIP states. This behavior is consistent with our NMR order parameters and scalar coupling data for the lysine side chains. Using the MD trajectories, we also analyze the free energies of the CIP-SIP equilibria. This work illustrates the dynamic nature of short-range electrostatic interactions in DNA recognition by proteins.

Intrinsic unfoldase/foldase activity of the chaperonin GroEL directly demonstrated using multinuclear relaxation-based NMR

The prototypical chaperonin GroEL assists protein folding through an ATP-dependent encapsulation mechanism. The details of how GroEL folds proteins remain elusive, particularly because encapsulation is not an absolute requirement for successful re/folding. Here we make use of a metastable model protein substrate, comprising a triple mutant of Fyn SH3, to directly demonstrate, by simultaneous analysis of three complementary NMR-based relaxation experiments (lifetime line broadening, dark state exchange saturation transfer, and Carr-Purcell-Meinboom-Gill relaxation dispersion), that apo GroEL accelerates the overall interconversion rate between the native state and a well-defined folding intermediate by about 20-fold, under conditions where the "invisible" GroEL-bound states have occupancies below 1%. This is largely achieved through a 500-fold acceleration in the folded-to-intermediate transition of the protein substrate. Catalysis is modulated by a kinetic deuterium isotope effect that reduces the overall interconversion rate between the GroEL-bound species by about 3-fold, indicative of a significant hydrophobic contribution. The location of the GroEL binding site on the folding intermediate, mapped from (15)N, (1)HN, and (13)Cmethyl relaxation dispersion experiments, is composed of a prominent, surface-exposed hydrophobic patch.

Coupling of Conformational Transitions in the N-terminal Domain of the 51-kDa FK506-binding Protein (FKBP51) Near Its Site of Interaction with the Steroid Receptor Proteins

Interchanging Leu-119 for Pro-119 at the tip of the β4-β5 loop in the first FK506 binding domain (FK1) of the FKBP51 and FKBP52 proteins, respectively, has been reported to largely reverse the inhibitory (FKBP51) or stimulatory (FKBP52) effects of these co-chaperones on the transcriptional activity of glucocorticoid and androgen receptor-protein complexes. Previous NMR relaxation studies have identified exchange line broadening, indicative of submillisecond conformational motion, throughout the β4-β5 loop in the FK1 domain of FKBP51, which are suppressed by the FKBP52-like L119P substitution. This substitution also attenuates exchange line broadening in the underlying β2 and β3a strands that is centered near a bifurcated main chain hydrogen bond interaction between these two strands. The present study demonstrates that these exchange line broadening effects arise from two distinct coupled conformational transitions, and the transition within the β2 and β3a strands samples a transient conformation that resembles the crystal structures of the selectively inhibited FK1 domain of FKBP51 recently reported. Although the crystal structures for their series of inhibitors were interpreted as evidence for an induced fit mechanism of association, the presence of a similar conformation being significantly populated in the unliganded FKBP51 domain is more consistent with a conformational selection binding process. The contrastingly reduced conformational plasticity of the corresponding FK1 domain of FKBP52 is consistent with the current model in which FKBP51 binds to both the apo- and hormone-bound forms of the steroid receptor to modulate its affinity for ligand, whereas FKBP52 binds selectively to the latter state.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

Characterizing RNA Excited States Using NMR Relaxation Dispersion

Changes in RNA secondary structure play fundamental roles in the cellular functions of a growing number of noncoding RNAs. This chapter describes NMR-based approaches for characterizing microsecond-to-millisecond changes in RNA secondary structure that are directed toward short-lived and low-populated species often referred to as "excited states." Compared to larger scale changes in RNA secondary structure, transitions toward excited states do not require assistance from chaperones, are often orders of magnitude faster, and are localized to a small number of nearby base pairs in and around noncanonical motifs. Here, we describe a procedure for characterizing RNA excited states using off-resonance R1ρ NMR relaxation dispersion utilizing low-to-high spin-lock fields (25-3000 Hz). R1ρ NMR relaxation dispersion experiments are used to measure carbon and nitrogen chemical shifts in base and sugar moieties of the excited state. The chemical shift data are then interpreted with the aid of secondary structure prediction to infer potential excited states that feature alternative secondary structures. Candidate structures are then tested by using mutations, single-atom substitutions, or by changing physiochemical conditions, such as pH and temperature, to either stabilize or destabilize the candidate excited state. The resulting chemical shifts of the mutants or under different physiochemical conditions are then compared to those of the ground and excited states. Application is illustrated with a focus on the transactivation response element from the human immune deficiency virus type 1, which exists in dynamic equilibrium with at least two distinct excited states.

Visualizing transient Watson-Crick-like mispairs in DNA and RNA duplexes

Rare tautomeric and anionic nucleobases are believed to have fundamental biological roles, but their prevalence and functional importance has remained elusive because they exist transiently, in low abundance, and involve subtle movements of protons that are difficult to visualize. Using NMR relaxation dispersion, we show here that wobble dG•dT and rG•rU mispairs in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-populated Watson-Crick-like mispairs that are stabilized by rare enolic or anionic bases. These mispairs can evade Watson-Crick fidelity checkpoints and form with probabilities (10(-3) to 10(-5)) that strongly imply a universal role in replication and translation errors. Our results indicate that rare tautomeric and anionic bases are widespread in nucleic acids, expanding their structural and functional complexity beyond that attainable with canonical bases.

Characterizing excited conformational states of RNA by NMR spectroscopy

Conformational dynamics is a hallmark of diverse non-coding RNA functions. During these functional processes, RNA molecules almost ubiquitously undergo conformational transitions that are tuned to meet distinct structural and kinetic requirements for proper function. A complete mechanistic understanding of RNA function requires comprehensive structural and dynamic knowledge of these complex transitions, which often involve alternative higher-energy conformational states that pose a major challenge for high-resolution structural study by conventional methods. In this review, we describe recent progress in RNA NMR that has started to unveil detailed structural, thermodynamic and kinetic insights into some of these excited conformational states of RNA and their functional roles in biology.