NMR Relaxation Dispersion Experiments Research Articles

Existence of Singularities in NMR Relaxation Dispersion Profiles: Implications for Hidden Dynamics.

It is common for NMR relaxation dispersion experiments to suggest an absence of dynamics despite anecdotal indications of conformational dynamics. We explore the potential explanations and approaches to this conundrum. Some inconsistencies have been observed between two relaxation dispersion experiments, Carr-Purcell-Meiboom-Gill (CPMG) and adiabatic relaxation dispersion experiments, in recent dynamic studies of different biomolecules. Theoretical analyses show that such seemingly paradoxical results might come from a complex exchange topology that is concealed by the application of the simple two-site exchange model for interpretation of the relaxation dispersion data. Scenarios are explored and revealed in which the presence of complex millisecond conformational exchange could suppress the amplitude of CPMG relaxation dispersion profiles, even when the exchange rates are within the detectable range of the experiments. With experimental errors, the suppressed relaxation dispersion profiles could lead to the conclusion of "no millisecond conformational exchange". However, such hidden dynamics can potentially be detected by adiabatic relaxation dispersion experiments. Finally, we demonstrate some advantages of adiabatic relaxation dispersion experiments over conventional relaxation dispersion experiments and a simplified computational approach to analyze the adiabatic relaxation dispersion profiles.

Transient excited states of the metamorphic protein Mad2 and their implications for function.

The spindle checkpoint complex is a key surveillance mechanism in cell division that prevents premature separation of sister chromatids. Mad2 is an integral component of this spindle checkpoint complex that recognizes cognate substrates such as Mad1 and Cdc20 in its closed (C-Mad2) conformation by fastening a "seatbelt" around short peptide regions that bind to the substrate recognition site. Mad2 is also a metamorphic protein that adopts not only the fold found in C-Mad2, but also a structurally distinct open conformation (O-Mad2) which is incapable of binding substrates. Here, we show using chemical exchange saturation transfer (CEST) and relaxation dispersion (CPMG) NMR experiments that Mad2 transiently populates three other higher free energy states with millisecond lifetimes, two in equilibrium with C-Mad2 (E1 and E2) and one with O-Mad2 (E3). E1 is a mimic of substrate-bound C-Mad2 in which the N-terminus of one C-Mad2 molecule inserts into the seatbelt region of a second molecule of C-Mad2, providing a potential pathway for autoinhibition of C-Mad2. E2 is the "unbuckled" conformation of C-Mad2 that facilitates the triage of molecules along competing fold-switching and substrate binding pathways. The E3 conformation that coexists with O-Mad2 shows fluctuations at a hydrophobic lock that is required for stabilizing the O-Mad2 fold and we hypothesize that E3 represents an early intermediate on-pathway towards conversion to C-Mad2. Collectively, the NMR data highlight the rugged free energy landscape of Mad2 with multiple low-lying intermediates that interlink substrate-binding and fold-switching, and also emphasize the role of molecular dynamics in its function.

Metal Binding of Alzheimer's Amyloid-β and Its Effect on Peptide Self-Assembly.

ConspectusMetal ions have been identified as key factors modulating the aggregation of amyloid-β peptide (Aβ) implicated in Alzheimer's disease (AD). The presence of elevated levels of metal ions in the amyloid plaques in AD patients supports the notion that the dysfunction of metal homeostasis is connected to the development of AD pathology. Here, recent findings from high- and low-resolution biophysical techniques are put into perspective, providing detailed insights into the molecular structures and dynamics of metal-bound Aβ complexes and the effect of metal ions on the Aβ aggregation process. In particular, the development of theoretical kinetic models deducing different microscopic nucleation events from the macroscopic aggregation behavior has enabled deciphering of the effect of metal ions on specific nucleation processes. In addition to these macroscopic measurements of bulk aggregation to quantify microscopic rates, recent NMR studies have revealed details about the structures and dynamics of metal-Aβ complexes, thereby linking structural events to bulk aggregation. Interestingly, transition-metal ions, such as copper, zinc, and silver ions, form a compact complex with the N-terminal part of monomeric Aβ, respectively, where the metal-bound "folded" state is in dynamic equilibrium with an "unfolded" state. The rates and thermodynamic features of these exchange dynamics have been determined by using NMR relaxation dispersion experiments. Additionally, the application of specifically tailored paramagnetic NMR experiments on the Cu(II)-Aβ complex has been fruitful in obtaining structural constraints within the blind sphere of conventional NMR experiments. This enables the determination of molecular structures of the "folded" Cu(II)-coordinated N-terminal region of Aβ. Furthermore, the discussed transition-metal ions modulate Aβ self-assembly in a concentration-dependent manner, where low metal ion concentrations inhibit Aβ fibril formation, while at high metal ion concentrations other processes occur, resulting in amorphous aggregate formation. Remarkably, the metal-Aβ interaction predominately reduces one specific nucleation step, the fibril-end elongation, whereas primary and surface-catalyzed secondary nucleation mechanisms are less affected. Specific inhibition of fibril-end elongation theoretically predicts an enhanced generation of Aβ oligomers, which is an interesting contribution to understanding metal-Aβ-associated neurotoxic effects. Taken together, the metal binding process creates a metal-bound Aβ complex, which is seemingly inert to aggregation. This process hence efficiently reduces the aggregation-prone peptide pool, which on the macroscopic level is reflected as slower aggregation kinetics. Thus, the specific binding of metals to the Aβ monomer can be linked to the macroscopic inhibitory effect on Aβ bulk aggregation, providing a molecular understanding of the Aβ aggregation mechanism in the presence of metal ions, where the metal ion can be seen as a minimalist agent against Aβ self-assembly. These insights can help to target Aβ aggregation in vivo, where metal ions are key factors modulating the Aβ self-assembly and Aβ-associated neurotoxicity.

An integrated approach of NMR experiments and MD simulations visualizes structural dynamics of a cyclic multi-domain protein.

Cyclization can stabilize the structure of proteins, as previously demonstrated in single-domain proteins. Although Lys48-linked polyubiquitin, a multi-domain protein, is also known to be cyclized in human cells, the structural effects of cyclization remain unclear. Here, we examined the impact of cyclization on the structural stability and dynamics of cyclic Lys48-linked diubiquitin (Ub2 ). As expected, cyclization increased the thermal stability of Ub2 and its resistance to proteolytic digestion, indicating that cyclization stabilized the structure of Ub2 . Furthermore, cyclization repressed the interdomain motion in Ub2 , but cyclic Ub2 still exhibited microsecond conformational exchange in NMR relaxation dispersion experiments. A series of long coarse-grained (CG) MD simulations visualized how cyclization slowed down the intrinsic nanosecond open-closed domain motion of Ub2 to microseconds. Thus, CG-MD analysis helped to explain the unexpected NMR relaxation results, thereby facilitating characterization of the structural stabilization of cyclic Ub2 .

Insights into the Structure of Invisible Conformations of Large Methyl Group Labeled Molecular Machines from High Pressure NMR

Most proteins are highly flexible and can adopt conformations that deviate from the energetically most favorable ground state. Structural information on these lowly populated, alternative conformations is often lacking, despite the functional importance of these states. Here, we study the pathway by which the Dcp1:Dcp2 mRNA decapping complex exchanges between an autoinhibited closed and an open conformation. We make use of methyl Carr–Purcell–Meiboom–Gill (CPMG) NMR relaxation dispersion (RD) experiments that report on the population of the sparsely populated open conformation as well as on the exchange rate between the two conformations. To obtain volumetric information on the open conformation as well as on the transition state structure we made use of RD measurements at elevated pressures. We found that the open Dcp1:Dcp2 conformation has a lower molecular volume than the closed conformation and that the transition state is close in volume to the closed state. In the presence of ATP the volume change upon opening of the complex increases and the volume of the transition state lies in-between the volumes of the closed and open state. These findings show that ATP has an effect on the volume changes that are associated with the opening-closing pathway of the complex. Our results highlight the strength of pressure dependent NMR methods to obtain insights into structural features of protein conformations that are not directly observable. As our work makes use of methyl groups as NMR probes we conclude that the applied methodology is also applicable to high molecular weight complexes.

A "Steady-State" Relaxation Dispersion Nuclear Magnetic Resonance Experiment for Studies of Chemical Exchange in Degenerate 1H Transitions of Methyl Groups.

Degenerate spin-systems consisting of magnetically equivalent nuclear spins, such as a 1H3 spin-system in selectively 13CH3-labeled proteins, present considerable challenges for the design of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments to characterize chemical exchange on the micro-to-millisecond time-scale. Several approaches have been previously proposed for the elimination of deleterious artifacts observed in methyl 1H CPMG relaxation dispersion profiles obtained for (13C)1H3 groups. We describe an alternative, experimentally simple solution and design a "steady-state" methyl 1H CPMG scheme, where 90° or acute-angle (<90°) 1H radiofrequency pulses are applied after each CPMG echo in-phase with methyl 1H magnetization, resulting in the establishment of a "steady-state" for effective rates of magnetization decay. A simple computational procedure for quantitative analysis of the "steady-state" CPMG relaxation dispersion profiles is developed. The "steady-state" CPMG methodology is applied to two protein systems where exchange between major and minor species occurs in different regimes on the chemical shift time-scale.

Characterizing Fast Conformational Exchange of Aromatic Rings Using Residual Dipolar Couplings: Distinguishing Jumplike Flips from Other Exchange Mechanisms.

Aromatic ring flips are a hallmark of protein dynamics. They are experimentally studied by NMR spectroscopy, where recent advances have led to improved characterization across a wide range of time scales. Results on different proteins have been interpreted as continuous diffusive ring rotations or jumplike flips, leading to diverging views of the protein interior as being fluidlike or solidlike, respectively. It is challenging to distinguish between these mechanisms and other types of conformational exchange because chemical-shift-mediated line broadening provides only conclusive evidence for ring flips only if the system can be moved from the slow- to intermediate/fast-exchange regime. Moreover, whenever the chemical shift difference between the two symmetry-related sites is close to zero, it is not generally possible to determine the exchange time scale. Here we resolve these issues by measuring residual dipolar coupling (RDC)-mediated exchange contributions using NMR relaxation dispersion experiments on proteins dissolved in dilute liquid crystalline media. Excellent agreement is found between the experimental difference in RDC between the two symmetry-related sites and the value calculated from high-resolution X-ray structures, demonstrating that dynamics measured for F52 in the B1 domain of protein G reports on distinct, jumplike flips rather than other types of conformational exchange.

Inherent Minor Conformer of Bordetella Effector BteA Directs Chaperone-Mediated Unfolding.

The pathogen Bordetella pertussis uses a type-3 secretion system (T3SS) to inject its cytotoxic effector BteA into the host cell via a designated needle structure. Prior to injection BteA is bound to its cognate chaperone BtcA presumed to assist in effector unfolding en route to needle passage. We utilized NMR and EPR spectroscopy to uncover the molecular mechanism of BtcA-mediated unfolding of BteA. BtcA induces a global structural change in the effector, which adopts a more extended and partially unfolded conformation. EPR distance measurements further show that the structured helical-bundle form of free BteA exists in conformational equilibrium with a lowly populated minor species. The nature of this equilibrium was probed using NMR relaxation dispersion experiments. At 283 K structural effects are most pronounced for a contiguous surface spanning the A- and B-helices of BteA, extending at 303 K to a second surface including the D- and E-helices. Residues perturbed in the minor conformation coincide with those exhibiting a BtcA-induced increase in flexibility, identifying this conformation as the BtcA-bound form of the effector. Our findings hint at a conformational-selectivity mechanism for the chaperone interaction with the effector, a paradigm that may be common to effector-chaperones secretion complexes in this family of pathogens.

Solution structure and dynamics of the mitochondrial-targeted GTPase-activating protein (GAP) VopE by an integrated NMR/SAXS approach.

The bacterial pathogen Vibrio cholerae use a type III secretion system to inject effector proteins into a host cell. Recently, a putative Toxic GTPase Activating Protein (ToxGAP) called Vibrio outer protein E (VopE) was identified as a T3SS substrate and virulence factor that affected host mitochondrial dynamics and immune response. However, biophysical and structural characterization has been absent. Here, we describe solution NMR structure of the putative GTPase-activating protein (GAP) domain (73-204) of VopE. Using size exclusion chromatography coupled with small-angle x-ray scattering and residual dipolar coupling data, we restrained the MD process to efficiently determine the overall fold and improve the quality of the output calculated structures. Comparing the structure of VopE with other ToxGAP's revealed a similar overall fold with several features unique to VopE. Specifically, the "Bulge 1," α1 helix, and noteworthy "backside linker" elements on the N-terminus are dissimilar to the other ToxGAP's. By using NMR relaxation dispersion experiments, we demonstrate that these regions undergo motions on a > 6s-1 timescale. Based on the disposition of these mobile regions relative to the putative catalytic arginine residue, we hypothesize that the protein may undergo structural changes to bind cognate GTPases.

Conformational Transitions in Yeast Chorismate Mutase Important for Allosteric Regulation as Identified by Nuclear Magnetic Resonance Spectroscopy

Proteins fluctuate between different conformations in solution, and these conformational fluctuations can be important for protein function and allosteric regulation. The chorismate mutase from Saccharomyces cerevisiae (ScCM), a key enzyme in the biosynthesis of aromatic amino acids, is allosterically activated and inhibited by tryptophan and tyrosine, respectively. It was initially proposed that in the absence of effector, ScCM fluctuates between activated R and inhibited T conformations according to the Monod-Wyman-Changeux (MWC) model, although a more complex regulation pattern was later suggested by mutagenesis and kinetic data. Here we used NMR relaxation dispersion experiments to understand the conformational fluctuations on the microsecond-to-millisecond timescale that occur in ScCM. In the absence of allosteric effectors, ScCM did not exclusively exchange between T and R conformations, suggesting that the two-state MWC model is insufficient to explain conformational dynamics. Addition of tyrosine led to the quenching of much of the motion on this timescale, while new motions were identified in the presence of tryptophan. These new motions are consistent with conformational fluctuations into an alternative conformation that may be important for enzyme activity.

Studying protein folding landscape using NMR relaxation measurements with pressure perturbation

Repeat proteins greatly simplify the identification of the sequence determinants for protein folding cooperativity. We use the leucine rich repeat (LRR) domain of the tumor suppressor pp32 as a model. It is composed of five LRR with a capping domain on each of its termini. Previous pressure-dependent NMR experiments of cavity-containing mutants of pp32 revealed that the intermediates populated during pressure-induced unfolding depend upon the position of the cavity. In order to further characterize the nature of these intermediates, we carried out NMR relaxation experiments, relaxation dispersion (CPMG) and chemical exchange saturation transfer (CEST) experiments.

Theoretical classification of exchange geometries from the perspective of NMR relaxation dispersion

NMR relaxation dispersion experiments have been widely applied to probe important conformational exchange of macro-molecules in many biological systems. The current improvements in computational techniques as well as the theoretical breakthroughs make the quantitative data analysis of complex exchange models possible. However, the topology of a given exchange model is also one of the main factors affecting the solution of Bloch-McConnell equation. The lack of a theoretical analysis of the exchange topologies at n-site exchange hinders further progress of such data analysis. Here, using graph theory, we reveal the topological complexity of n-site exchange and present all exchange models when n is less than 6. Furthermore, we introduce an alternative way, using machine learning, to select an exchange model based on a set of relaxation dispersion data without fitting them with every individual exchange model.

Chemical shift prediction of RNA imino groups: application toward characterizing RNA excited states

NH groups in proteins or nucleic acids are the most challenging target for chemical shift prediction. Here we show that the RNA base pair triplet motif dictates imino chemical shifts in its central base pair. A lookup table is established that links each type of base pair triplet to experimental chemical shifts of the central base pair, and can be used to predict imino chemical shifts of RNAs to remarkable accuracy. Strikingly, the semiempirical method can well interpret the variations of chemical shifts for different base pair triplets, and is even applicable to non-canonical motifs. This finding opens an avenue for predicting chemical shifts of more complicated RNA motifs. Furthermore, we combine the imino chemical shift prediction with NMR relaxation dispersion experiments targeting both 15N and 1HN of the imino group, and verify a previously characterized excited state of P5abc subdomain including an earlier speculated non-native G•G mismatch.

Mapping the energy landscape of protein-ligand binding via linear free energy relationships determined by protein NMR relaxation dispersion.

Biochemical signaling is mediated by complexes between macromolecular receptors and their ligands, with the duration of the signal being directly related to the lifetime of the ligand–receptor complex. In the field of drug design, the recognition that drug efficacy in vivo depends on the lifetime of the drug–protein complex has spawned the concept of designing drugs with particular binding kinetics. To advance this field it is critical to investigate how the molecular details of designed ligands might affect the binding kinetics, as well as the equilibrium binding constant. Here we use protein NMR relaxation dispersion to determine linear free energy relationships involving the on- and off-rates and the affinity for a series of congeneric ligands targeting the carbohydrate recognition domain of galectin-3. Using this approach we determine the energy landscape and the position of the transition state along the reaction coordinate of protein–ligand binding. The results show that ligands exhibiting reduced off-rates achieve this by primarily stabilizing the bound state, but do not affect the transition state to any greater extent. The transition state forms early, that is, it is located significantly closer to the free state than to the bound state, suggesting a critical role of desolvation. Furthermore, the data suggest that different subclasses of ligands show different behavior with respect to these characteristics.

Structural dynamics of double-stranded DNA with epigenome modification.

Modification of cytosine plays an important role in epigenetic regulation of gene expression and genome stability. Cytosine is converted to 5-methylcytosine (5mC) by DNA methyltransferase; in turn, 5mC may be oxidized to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation enzyme. The structural flexibility of DNA is known to affect the binding of proteins to methylated DNA. Here, we have carried out a semi-quantitative analysis of the dynamics of double-stranded DNA (dsDNA) containing various epigenetic modifications by combining data from imino 1H exchange and imino 1H R1ρ relaxation dispersion NMR experiments in a complementary way. Using this approach, we characterized the base-opening (kopen) and base-closing (kclose) rates, facilitating a comparison of the base-opening and -closing process of dsDNA containing cytosine in different states of epigenetic modification. A particularly striking result is the increase in the kopen rate of hemi-methylated dsDNA 5mC/C relative to unmodified or fully methylated dsDNA, indicating that the Watson–Crick base pairs undergo selective destabilization in 5mC/C. Collectively, our findings imply that the epigenetic modulation of cytosine dynamics in dsDNA mediates destabilization of the GC Watson–Crick base pair to allow base-flipping in living cells.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe–S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol−1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

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

The rotating-frame spin relaxation rate constant, R1ρ, is a powerful probe of macromolecular chemical and conformational dynamics in relaxation dispersion, CEST, and DEST NMR experiments. The R1ρ 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.

Direct evidence for (G)O6···H2-N4(C)+ hydrogen bonding in transient G(syn)-C+ and G(syn)-m5C+ Hoogsteen base pairs in duplex DNA from cytosine amino nitrogen off-resonance R1ρ relaxation dispersion measurements.

NMR relaxation dispersion studies have shown that Watson-Crick G-C and A-T base pairs in duplex DNA exist in dynamic equilibrium with their Hoogsteen counterparts. Hoogsteen base pairs form through concurrent rotation of the purine base about the glycosidic bond from an anti to a syn conformation and constriction of the C1'-C1' distance across the base pair by ∼2 Å to allow Hoogsteen type hydrogen bonding. Owing to their unique structure, Hoogsteen base pairs can play important roles in DNA recognition, the accommodation, recognition, and repair of DNA damage, and in DNA replication. NMR relaxation dispersion experiments targeting imino nitrogen and protonated base and sugar carbons have provided insights into many structural features of transient Hoogsteen base pairs, including one of two predicted hydrogen bonds involving (G)N7···H-N3(C)+ and (A)N7···H-N3(T). Here, through measurement of cytosine amino (N4) R1ρ relaxation dispersion, we provide direct evidence for the second (G)O6···H2-N4(C)+ hydrogen bond in G(syn)-C+ transient Hoogsteen base pairs. The utility of cytosine N4 R1ρ relaxation dispersion as a new sensitive probe of transient Hoogsteen base pairs, and cytosine dynamics in general, is further demonstrated by measuring G(syn)-C+ Hoogsteen exchange near neutral pH and in the context of the naturally occurring DNA modification 5-methyl cytosine (m5C), in DNA samples prepared using chemical synthesis and a 15N labeled m5C phosphoramidite.

Accelerating 2D NMR relaxation dispersion experiments using iterated maps.

NMR relaxation dispersion experiments play a central role in exploring molecular motion over an important range of timescales, and are an example of a broader class of multidimensional NMR experiments that probe important biomolecules. However, resolving the spectral features of these experiments using the Fourier transform requires sampling the full Nyquist grid of data, making these experiments very costly in time. Practitioners often reduce the experiment time by omitting 1D experiments in the indirectly observed dimensions, and reconstructing the spectra using one of a variety of post-processing algorithms. In prior work, we described a fast, Fourier-based reconstruction method using iterated maps according to the Difference Map algorithm of Veit Elser (DiffMap). Here we describe coDiffMap, a new reconstruction method that is based on DiffMap, but which exploits the strong correlations between 2D data slices in a pseudo-3D experiment. We apply coDiffMap to reconstruct dispersion curves from an [Formula: see text] relaxation dispersion experiment, and demonstrate that the method provides fast reconstructions and accurate relaxation curves down to very low numbers of sparsely-sampled data points.

Equilibrium and Kinetic Unfolding of GB1: Stabilization of the Native State by Pressure.

NMR spectroscopy allows an all-atom view on pressure-induced protein folding, separate detection of different folding states, determination of their population, and the measurement of the folding kinetics at equilibrium. Here, we studied the folding of protein GB1 at pH 2 in a temperature and pressure dependent way. We find that the midpoints of temperature-induced unfolding increase with higher pressure. NMR relaxation dispersion experiments disclosed that the unfolding kinetics slow down at elevated pressure while the folding kinetics stay virtually the same. Therefore, pressure is stabilizing the native state of GB1. These findings extend the knowledge of the influence of pressure on protein folding kinetics, where so far typically a destabilization by increased activation volumes of folding was observed. Our findings thus point toward an exceptional section in the pressure-temperature phase diagram of protein unfolding. The stabilization of the native state could potentially be caused by a shift of p Ka values of glutamates and aspartates in favor of the negatively charged state as judged from pH sensitive chemical shifts.