Measurement of methyl axis orientations in invisible, excited states of proteins by relaxation dispersion NMR spectroscopy

Characterization of a transitionally occupied state and thermal unfolding of domain 1.1 of σ A factor of RNA polymerase from Bacillus subtilis

σ factors are essential parts of bacterial RNA polymerase (RNAP) as they allow to recognize promotor sequences and initiate transcription. Domain 1.1 of vegetative σ factors occupies the primary channel of RNAP and also prevents binding of the σ factor to promoter DNA alone. Here, we show that domain 1.1 of Bacillus subtilisσ A exists in more structurally distinct variants in dynamic equilibrium. The major conformation at room temperature is represented by a previously reported well-folded structure solved by nuclear magnetic resonance (NMR), but 4% of the protein molecules are present in a less thermodynamically favorable state. We show that this population increases with temperature and we predict its significant elevation at higher but still biologically relevant temperatures. We characterized the minor state of the domain 1.1 using specialized methods of NMR. We found that, in contrast to the major state, the detected minor state is partially unfolded. Its propensity to form secondary structure elements is especially decreased for the first and third α helices, while the second α helix and β strand close to the C-terminus are more stable. We also analyzed thermal unfolding of the domain 1.1 and performed functional experiments with full lengthσ A and its shortened version lacking domain 1.1 (σ A _ Δ 1.1 ). The results revealed that while full lengthσ A increases transcription activity of RNAP with increasing temperature, transcription withσ A _ Δ 1.1 remains constant. In summary, this study reveals conformational dynamics of domain 1.1 and provides a basis for studies of its interaction with RNAP and effects on transcription regulation.

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.

Systematic Approach to Find the Global Minimum of Relaxation Dispersion Data for Protein-Induced B-Z Transition of DNA

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion spectroscopy is commonly used for quantifying conformational changes of protein in μs-to-ms timescale transitions. To elucidate the dynamics and mechanism of protein binding, parameters implementing CPMG relaxation dispersion results must be appropriately determined. Building an analytical model for multi-state transitions is particularly complex. In this study, we developed a new global search algorithm that incorporates a random search approach combined with a field-dependent global parameterization method. The robust inter-dependence of the parameters carrying out the global search for individual residues (GSIR) or the global search for total residues (GSTR) provides information on the global minimum of the conformational transition process of the Zα domain of human ADAR1 (hZαADAR1)-DNA complex. The global search results indicated that a α-helical segment of hZαADAR1 provided the main contribution to the three-state conformational changes of a hZαADAR1-DNA complex with a slow B-Z exchange process. The two global exchange rate constants, kex and kZB, were found to be 844 and 9.8 s-1, respectively, in agreement with two regimes of residue-dependent chemical shift differences-the "dominant oscillatory regime" and "semi-oscillatory regime". We anticipate that our global search approach will lead to the development of quantification methods for conformational changes not only in Z-DNA binding protein (ZBP) binding interactions but also in various protein binding processes.

A Single-Point Mutation in d-Arginine Dehydrogenase Unlocks a Transient Conformational State Resulting in Altered Cofactor Reactivity

Proteins are inherently dynamic, and proper enzyme function relies on conformational flexibility. In this study, we demonstrated how an active site residue changes an enzyme's reactivity by modulating fluctuations between conformational states. Replacement of tyrosine 249 (Y249) with phenylalanine in the active site of the flavin-dependent d-arginine dehydrogenase yielded an enzyme with both an active yellow FAD (Y249F-y) and an inactive chemically modified green FAD, identified as 6-OH-FAD (Y249F-g) through various spectroscopic techniques. Structural investigation of Y249F-g and Y249F-y variants by comparison to the wild-type enzyme showed no differences in the overall protein structure and fold. A closer observation of the active site of the Y249F-y enzyme revealed an alternative conformation for some active site residues and the flavin cofactor. Molecular dynamics simulations probed the alternate conformations observed in the Y249F-y enzyme structure and showed that the enzyme variant with FAD samples a metastable conformational state, not available to the wild-type enzyme. Hybrid quantum/molecular mechanical calculations identified differences in flavin electronics between the wild type and the alternate conformation of the Y249F-y enzyme. The computational studies further indicated that the alternate conformation in the Y249F-y enzyme is responsible for the higher spin density at the C6 atom of flavin, which is consistent with the formation of 6-OH-FAD in the variant enzyme. The observations in this study are consistent with an alternate conformational space that results in fine-tuning the microenvironment around a versatile cofactor playing a critical role in enzyme function.

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.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics

NMR spectroscopy is a powerful tool for studying protein dynamics. Conventionally, NMR studies on protein dynamics have probed motions of protein backbone NH, side-chain aromatic, and CH₃ groups. Recently, there has been remarkable progress in NMR methodologies that can characterize motions of cationic groups in protein side chains. These NMR methods allow investigations of the dynamics of positively charged lysine (Lys) and arginine (Arg) side chains and their hydrogen bonds as well as their electrostatic interactions important for protein function. Here, describing various practical aspects, we provide an overview of the NMR methods for dynamics studies of Lys and Arg side chains. Some example data on protein-DNA complexes are shown. We will also explain how molecular dynamics (MD) simulations can facilitate the interpretation of the NMR data on these basic side chains. Studies combining NMR and MD have revealed the highly dynamic nature of short-range electrostatic interactions via ion pairs, especially those involving Lys side chains.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

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.

Self-Assembly of Human Profilin-1 Detected by Carr-Purcell-Meiboom-Gill Nuclear Magnetic Resonance (CPMG NMR) Spectroscopy

Protein oligomerization in the cell has important implications for both health and disease, and an understanding of the mechanisms by which proteins can self-associate is, therefore, of critical interest. Initial stages of the oligomerization process can be hard to detect, as they often involve the formation of sparsely populated and transient states that are difficult to characterize by standard biophysical approaches. Using relaxation dispersion nuclear magnetic resonance spectroscopy, we study the oligomerization of human profilin-1, a protein that regulates the polymerization of actin. We show that in solution and at millimolar concentrations profilin-1 is predominantly monomeric. However, fits of concentration-dependent relaxation data are consistent with the formation of a higher-order oligomer that is generated via a multistep process. Together with crystallographic data for profilin-2, a homologue of the protein studied here, our results suggest that profilin-1 forms a sparsely populated tetrameric conformer in solution.

NMR studies on the dynamics of hydrogen bonds and ion pairs involving lysine side chains of proteins

Hydrogen bonds and ion pairs involving side chains play vital roles in protein functions such as molecular recognition and catalysis. Despite the wealth of structural information about hydrogen bonds and ion pairs at functionally crucial sites on proteins, the dynamics of these fundamental chemical interactions are not well understood largely due to the lack of suitable experimental tools in the past. NMR spectroscopy is a powerful tool for investigations of protein dynamics, but the vast majority of NMR methods had been applicable only to the backbone or methyl groups. Recently, a substantial progress has been made in the research on the dynamics of hydrogen bonds and ion pairs involving lysine side-chain NH3+ groups. Together with computational/theoretical approaches, the new NMR methods provide unique insights into the dynamics of hydrogen bonds and ion pairs involving lysine side chains. Here, the methodology and its applications are reviewed.

Harnessing NMR relaxation interference effects to characterise supramolecular assemblies

Supramolecular assemblies can show unusual splitting patterns in their NMR spectra, with the magnitude of the effect increasing with molecular size. A simple NMR relaxation experiment reveals the origin of the effect and provides information on tumbling rates and local electronic structure.

An exact solution for R2,eff in CPMG experiments in the case of two site chemical exchange

The Carr-Purcell-Meiboom-Gill (CPMG) experiment is widely used to quantitatively analyse the effects of chemical exchange on NMR spectra. In a CPMG experiment, the effective transverse relaxation rate, R2,eff, is typically measured as a function of the pulse frequency, νCPMG. Here, an exact expression for how R2,eff varies with νCPMG is derived for the commonly encountered scenario of two-site chemical exchange of in-phase magnetisation. This result, summarised in Appendix A, generalises a frequently used equation derived by Carver and Richards, published in 1972. The expression enables more rapid analysis of CPMG data by both speeding up calculation of R2,eff over numerical methods by a factor of ca. 130, and yields exact derivatives for use in data analysis. Moreover, the derivation provides insight into the physical principles behind the experiment.

Analyzing protein-ligand interactions by dynamic NMR spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy can provide detailed information on protein-ligand interactions that is inaccessible using other biophysical techniques. This chapter focuses on NMR-based approaches for extracting affinity and rate constants for weakly binding transient protein complexes with lifetimes of less than about a second. Several pulse sequences and analytical techniques are discussed, including line-shape simulations, spin-echo relaxation dispersion methods (CPMG), and magnetization exchange (EXSY) experiments.

Solution NMR and computational methods for understanding protein allostery

Allosterism is an essential biological regulatory mechanism. In enzymes, allosteric regulation results in an activation or inhibition of catalytic turnover. The mechanisms by which this is accomplished are unclear and vary significantly depending on the enzyme. It is commonly the case that a metabolite binds to the enzyme at a site distant from the catalytic site, yet its binding is coupled to and sensed by the active site. This coupling can manifest in changes in structure, dynamics, or both at the active site. These interactions between the allosteric and active site, which are often quite distant from one another, involve numerous atoms as well as complex conformational rearrangements of the protein secondary and tertiary structure. Interrogation of this complex biological phenomenon necessitates multiple experimental approaches. In this article, we outline a combined solution NMR spectroscopic and computational approach using molecular dynamics and network models to uncover mechanistic aspects of allostery in the enzyme imidazole glycerol phosphate synthase.

Structural insights into the Trp-cage folding intermediate formation

The 20 residue long Trp-cage is the smallest protein known, and thus has been the subject of several in vitro and in silico folding studies. Here, we report the multistate folding scenario of the miniprotein in atomic detail. We detected and characterized different intermediate states by temperature dependent NMR measurements of the (15)N and (13)C/(15)N labeled protein, both at neutral and acidic pH values. We developed a deconvolution technique to characterize the invisible--fully folded, unfolded and intermediate--fast exchanging states. Using nonlinear fitting methods we can obtain both the thermodynamic parameters (ΔH(F-I), T(m)(F-I), ΔC(p)(F-I) and ΔH(I-U), T(m)(I-U), ΔC(p)(I-U)) and the NMR chemical shifts of the conformers of the multistate unfolding process. During the unfolding of Trp-cage distinct intermediates evolve: a fast-exchanging intermediate is present under neutral conditions, whereas a slow-exchanging intermediate-pair emerges at acidic pH. The fast-exchanging intermediate has a native-like structure with a short α-helix in the G(11)-G(15) segment, whereas the slow-exchanging intermediate-pair presents elevated dynamics, with no detectable native-like residue contacts in which the G(11)-P(12) peptide bond has either cis or trans conformation. Heteronuclear relaxation studies combined with MD simulations revealed the source of backbone mobility and the nature of structural rearrangements during these transitions. The ability to detect structural and dynamic information about folding intermediates in vitro provides an excellent opportunity to gain new insights into the energetic aspects of the energy landscape of protein folding. Our new experimental data offer exceptional testing ground for further computational simulations.

Visualizing transient low-populated structures of RNA

The visualization of RNA conformational changes has provided fundamental insights into how regulatory RNAs carry out their biological functions. The RNA structural transitions that have been characterized to date involve long-lived species that can be captured by structure characterization techniques. Here, we report the Nuclear Magnetic Resonance visualization of RNA transitions towards invisible ‘excited states’ (ES), which exist in too little abundance (2–13%) and for too short periods of time (45–250 μs) to allow structural characterization by conventional techniques. Transitions towards ESs result in localized rearrangements in base-pairing that alter building block elements of RNA architecture, including helix-junction-helix motifs and apical loops. The ES can inhibit function by sequestering residues involved in recognition and signaling or promote ATP-independent strand exchange. Thus, RNAs do not adopt a single conformation, but rather exist in rapid equilibrium with alternative ESs, which can be stabilized by cellular cues to affect functional outcomes.

Measurement of the signs of methyl 13C chemical shift differences between interconverting ground and excited protein states by R(1ρ): an application to αB-crystallin

Carr-Purcell-Meiboom-Gill relaxation dispersion (CPMG RD) NMR spectroscopy has emerged as a powerful tool for quantifying the kinetics and thermodynamics of millisecond time-scale exchange processes involving the interconversion between a visible ground state and one or more minor, sparsely populated invisible 'excited' conformational states. Recently it has also become possible to determine atomic resolution structural models of excited states using a wide array of CPMG RD approaches. Analysis of CPMG RD datasets provides the magnitudes of the chemical shift differences between the ground and excited states, Δϖ, but not the sign. In order to obtain detailed structural insights from, for example, excited state chemical shifts and residual dipolar coupling measurements, these signs are required. Here we present an NMR experiment for obtaining signs of (13)C chemical shift differences of (13)CH(3) methyl groups using weak field off-resonance R(1ρ) relaxation measurements. The accuracy of the method is established by using an exchanging system where the invisible, excited state can be converted to the visible, ground state by altering sample conditions so that the signs of Δϖ values obtained from the spin-lock approach can be validated against those measured directly. Further, the spin-lock experiments are compared with the established H(S/M)QC approach for measuring signs of chemical shift differences and the relative strengths of each method are discussed. In the case of the 650 kDa human αB-crystallin complex where there are large transverse relaxation differences between ground and excited state spins the R(1ρ) method is shown to be superior to more 'traditional' experiments for sign determination.

Dynamics of lysine side-chain amino groups in a protein studied by heteronuclear 1H−15N NMR spectroscopy

Despite their importance in macromolecular interactions and functions, the dynamics of lysine side-chain amino groups in proteins are not well understood. In this study, we have developed the methodology for the investigations of the dynamics of lysine NH3(+) groups by NMR spectroscopy and computation. By using 1H−15N heteronuclear correlation experiments optimized for 15NH3(+) moieties, we have analyzed the dynamic behavior of individual lysine NH3(+) groups in human ubiquitin at 2 °C and pH 5. We modified the theoretical framework developed previously for CH3 groups and used it to analyze 15N relaxation data for the NH3(+) groups. For six lysine NH3(+) groups out of seven in ubiquitin, we have determined model-free order parameters, correlation times for bond rotation, and reorientation of the symmetry axis occurring on a pico- to nanosecond time scale. From CPMG relaxation dispersion experiment for lysine NH3(+) groups, slower dynamics occurring on a millisecond time scale have also been detected for Lys27. The NH3(+) groups of Lys48, which plays a key role as the linkage site in ubiquitination for proteasomal degradation, was found to be highly mobile with the lowest order parameter among the six NH3(+) groups analyzed by NMR. We compared the experimental order parameters for the lysine NH3(+) groups with those from a 1 μs molecular dynamics simulation in explicit solvent and found good agreement between the two. Furthermore, both the computer simulation and the experimental correlation times for the bond rotations of NH3(+) groups suggest that their hydrogen bonding is highly dynamic with a subnanosecond lifetime. This study demonstrates the utility of combining NMR experiment and simulation for an in-depth characterization of the dynamics of these functionally most important side-chains of ubiquitin.

Paramagnetic-based NMR restraints lift residual dipolar coupling degeneracy in multidomain detergent-solubilized membrane proteins

Residual dipolar couplings (RDCs) are widely used as orientation-dependent NMR restraints to improve the resolution of the NMR conformational ensemble of biomacromolecules and define the relative orientation of multidomain proteins and protein complexes. However, the interpretation of RDCs is complicated by the intrinsic degeneracy of analytical solutions and protein dynamics that lead to ill-defined orientations of the structural domains (ghost orientations). Here, we illustrate how restraints from paramagnetic relaxation enhancement (PRE) experiments lift the orientational ambiguity of multidomain membrane proteins solubilized in detergent micelles. We tested this approach on monomeric phospholamban (PLN), a 52-residue membrane protein, which is composed of two helical domains connected by a flexible loop. We show that the combination of classical solution NMR restraints (NOEs and dihedral angles) with RDC and PRE constraints resolves topological ambiguities, improving the convergence of the PLN structural ensemble and giving the depth of insertion of the protein within the micelle. The combination of RDCs with PREs will be necessary for improving the accuracy and precision of membrane protein conformational ensembles, where three-dimensional structures are dictated by interactions with the membrane-mimicking environment rather than compact tertiary folds common in globular proteins.

13CHD2 methyl group probes of millisecond time scale exchange in proteins by 1H relaxation dispersion: an application to proteasome gating residue dynamics

A pulse scheme is presented for quantifying millisecond time scale chemical exchange processes in proteins by measuring (1)H CPMG relaxation dispersion profiles of (13)CHD(2) methyl groups. The use of (13)CHD(2) isotopomers for (1)H methyl dispersion experiments eliminates problems with interconversion between differentially relaxing proton transitions that complicate the extraction of accurate exchange parameters when (13)CH(3) probes are used. Good agreement is demonstrated between extracted chemical shift differences from fits of dispersion profiles and the corresponding differences measured independently on a model exchanging system, validating the experiment. The methodology is applied to the gating residues of the T. acidiphilium proteasome that are shown to undergo extensive motion on the millisecond time scale.

An introduction to NMR-based approaches for measuring protein dynamics

Proteins are inherently flexible at ambient temperature. At equilibrium, they are characterized by a set of conformations that undergo continuous exchange within a hierarchy of spatial and temporal scales ranging from nanometers to micrometers and femtoseconds to hours. Dynamic properties of proteins are essential for describing the structural bases of their biological functions including catalysis, binding, regulation and cellular structure. Nuclear magnetic resonance (NMR) spectroscopy represents a powerful technique for measuring these essential features of proteins. Here we provide an introduction to NMR-based approaches for studying protein dynamics, highlighting eight distinct methods with recent examples, contextualized within a common experimental and analytical framework. The selected methods are (1) Real-time NMR, (2) Exchange spectroscopy, (3) Lineshape analysis, (4) CPMG relaxation dispersion, (5) Rotating frame relaxation dispersion, (6) Nuclear spin relaxation, (7) Residual dipolar coupling, (8) Paramagnetic relaxation enhancement. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

NMR spectroscopy brings invisible protein states into focus

Molecular dynamics are essential for protein function. In some cases these dynamics involve the interconversion between ground state, highly populated conformers and less populated higher energy structures ('excited states') that play critical roles in biochemical processes. Here we describe recent advances in NMR spectroscopy methods that enable studies of these otherwise invisible excited states at an atomic level and that help elucidate their important relation to function. We discuss a range of examples from molecular recognition, ligand binding, enzyme catalysis and protein folding that illustrate the role that motion plays in 'funneling' conformers along preferred pathways that facilitate their biological function.