Isotropic Reorientational Eigenmode Dynamics Complements NMR Relaxation Measurements for RNA

NMR studies of nucleic acid dynamics

Nucleic acid structures have to satisfy two diametrically opposite requirements; on one hand they have to adopt well-defined 3D structures that can be specifically recognized by proteins; on the other hand, their structures must be sufficiently flexible to undergo very large conformational changes that are required during key biochemical processes, including replication, transcription, and translation. How do nucleic acids introduce flexibility into their 3D structure without losing biological specificity? Here, I describe the development and application of NMR spectroscopic techniques in my laboratory for characterizing the dynamic properties of nucleic acids that tightly integrate a broad set of NMR measurements, including residual dipolar couplings, spin relaxation, and relaxation dispersion with sample engineering and computational approaches. This approach allowed us to obtain fundamental new insights into directional flexibility in nucleic acids that enable their structures to change in a very specific functional manner.

Domain-elongation NMR spectroscopy yields new insights into RNA dynamics and adaptive recognition

By simplifying the interpretation of nuclear magnetic resonance spin relaxation and residual dipolar couplings data, recent developments involving the elongation of RNA helices are providing new atomic insights into the dynamical properties that allow RNA structures to change functionally and adaptively. Domain elongation, in concert with spin relaxation measurements, has allowed the detailed characterization of a hierarchical network of local and collective motional modes occurring at nanosecond timescale that mirror the structural rearrangements that take place following adaptive recognition. The combination of domain elongation with residual dipolar coupling measurements has allowed the experimental three-dimensional visualization of very large amplitude rigid-body helix motions in HIV-1 transactivation response element (TAR) that trace out a highly choreographed trajectory in which the helices twist and bend in a correlated manner. The dynamic trajectory allows unbound TAR to sample many of its ligand bound conformations, indicating that adaptive recognition occurs by "conformational selection" rather than "induced fit." These studies suggest that intrinsic flexibility plays essential roles directing RNA conformational changes along specific pathways.

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.

Pressure effects on the ensemble dynamics of ubiquitin inspected with molecular dynamics simulations and isotropic reorientational eigenmode dynamics

According to NMR chemical shift data, the ensemble of ubiquitin is a mixture of "open" and "closed" conformations at rapid equilibrium. Pressure perturbations provide the means to study the transition between the two conformers by imposing an additional constraint on the system's partial molar volume. Here we use nanosecond-timescale molecular dynamics simulations to characterize the network of correlated motions accessible to the conformers at low- and high-pressure conditions. Using the isotropic reorientational eigenmode dynamics formalism to analyze our simulation trajectories, we reproduce NMR relaxation data without fitting any parameters of our model. Comparative analysis of our results suggests that the two conformations behave very differently. The dynamics of the "closed" conformation are almost unaffected by pressure and are dominated by large-amplitude correlated motions of residues 23-34 in the extended alpha-helix. The "open" conformation under conditions of normal pressure displays increased mobility, focused on the loop residues 17-20, 46-55, and 58-59 at the bottom of the core of the structure, as well as the C-terminal residues 69-76, that directly participate in key protein-protein interactions. For the same conformation, a pressure increase induces a loss of separability between molecular tumbling and internal dynamics, while motions between different backbone sites become uncorrelated.

Probing chemical shifts of invisible states of proteins with relaxation dispersion NMR spectroscopy: how well can we do?

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy has evolved into a powerful approach for the study of low populated, invisible conformations of biological molecules. One of the powerful features of the experiment is that chemical shift differences between the exchanging conformers can be obtained, providing structural information about invisible excited states. Through the development of new labeling approaches and NMR experiments it is now possible to measure backbone 13C(alpha) and 13CO relaxation dispersion profiles in proteins without complications from 13C-13C couplings. Such measurements are presented here, along with those that probe exchange using 15N and 1HN nuclei. A key experimental design has been the choice of an exchanging system where excited-state chemical shifts were known from independent measurement. Thus it is possible to evaluate quantitatively the accuracy of chemical shift differences obtained in dispersion experiments and to establish that in general very accurate values can be obtained. The experimental work is supplemented by computations that suggest that similarly accurate shifts can be measured in many cases for systems with exchange rates and populations that fall within the range of those that can be quantified by relaxation dispersion. The accuracy of the extracted chemical shifts opens up the possibility of obtaining quantitative structural information of invisible states of the sort that is now available from chemical shifts recorded on ground states of proteins.

An improved 15N relaxation dispersion experiment for the measurement of millisecond time-scale dynamics in proteins

A new (15)N constant-time relaxation dispersion pulse scheme for the quantification of millisecond time-scale exchange dynamics in proteins is presented. The experiment differs from previously developed sequences in that it includes (1)H continuous-wave decoupling during the (15)N Carr-Purcell-Meiboom-Gill (CPMG) pulse train that significantly improves the relaxation properties of (15)N magnetization, leading to sensitivity gains in experiments. Moreover, it is shown that inclusion of an additional (15)N 180 degrees refocusing pulse (phase cycled +/- x) in the center of the CPMG pulse train, consisting of 1(5)N 180 degrees (y) pulses, provides compensation for pulse imperfections beyond the normal CPMG scheme. Relative to existing relaxation-compensated constant-time relaxation dispersion pulse schemes, nu(CPMG) values that are only half as large can be employed, offering increased sensitivity to slow time-scale exchange processes. The robustness of the methodology is illustrated with applications involving a pair of proteins: an SH3 domain that does not show millisecond time-scale exchange and an FF domain with significant chemical exchange contributions.

An exchange-free measure of 15N transverse relaxation: an NMR spectroscopy application to the study of a folding intermediate with pervasive chemical exchange

A series of experiments are presented that provide an exchange-free measure of dipole-dipole (15)N transverse relaxation, R(dd), that can then be substituted for (15)N R(1rho) or R(2) rates in the study of internal protein dynamics. The method is predicated on the measurement of a series of relaxation rates involving (1)H-(15)N longitudinal order, anti-phase (1)H and (15)N single-quantum coherences, and (1)H-(15)N multiple quantum coherences; the relaxation rates of all coherences are measured under conditions of spin-locking. Results from detailed simulations and experiments on a number of protein systems establish that R(dd) values are independent of exchange and systematic errors from dipolar interactions with proximal protons are calculated to be less than 1-2%, on average, for applications to perdeuterated proteins. Simulations further indicate that the methodology is rather insensitive to the exact level of deuteration so long as proteins are reasonably highly deuterated (>50%). The utility of the methodology is demonstrated with applications involving protein L, ubiquitin, and a stabilized folding intermediate of apocytochrome b(562) that shows large contributions to (15)N R(1rho) relaxation from chemical exchange.

Separating degenerate (1)H transitions in methyl group probes for single-quantum (1)H-CPMG relaxation dispersion NMR spectroscopy

A relaxation dispersion pulse scheme is presented for quantifying chemical exchange processes in proteins that exploits 1H chemical shifts as probes of changes in conformation. The experiment selects 1H single-quantum magnetization from the I = 1/2 manifolds of the methyl group, which behave like AX spin systems, while suppressing coherences that derive from the 3/2 manifold that are extremely sensitive to pulse imperfections and that would otherwise severely compromise the accuracy of the experiment. The utility of the sequence is first demonstrated with an application to a protein system that is known not to undergo chemical exchange and flat dispersion profiles are obtained. Subsequently, the methodology is applied to study the folding of a G48M mutant of the Fyn SH3 domain that has been shown previously to undergo exchange between folded and unfolded states on the millisecond time scale.

Structural plasticity and Mg2+ binding properties of RNase P P4 from combined analysis of NMR residual dipolar couplings and motionally decoupled spin relaxation

The P4 helix is an essential element of ribonuclease P (RNase P) that is believed to bind catalytically important metals. Here, we applied a combination of NMR residual dipolar couplings (RDCs) and a recently introduced domain-elongation strategy for measuring "motionally decoupled" relaxation data to characterize the structural dynamics of the P4 helix from Bacillus subtilis RNase P. In the absence of divalent ions, the two P4 helical domains undergo small amplitude (approximately 13 degrees) collective motions about an average interhelical angle of 10 degrees. The highly conserved U7 bulge and helical residue C8, which are proposed to be important for substrate recognition and metal binding, are locally mobile at pico- to nanosecond timescales and together form the pivot point for the collective domain motions. Chemical shift mapping reveals significant association of Mg2+ ions at the P4 major groove near the flexible pivot point at residues (A5, G22, G23) previously identified to bind catalytically important metals. The Mg2+ ions do not, however, significantly alter the structure or dynamics of P4. Analysis of results in the context of available X-ray structures of the RNA component of RNase P and structural models that include the pre-tRNA substrate suggest that the internal motions observed in P4 likely facilitate adaptive changes in conformation that take place during folding and substrate recognition, possibly aided by interactions with Mg2+ ions. Our results add to a growing view supporting the existence of functionally important internal motions in RNA occurring at nanosecond timescales.

Off-Resonance TROSY-Selected R1ρ Experiment with Improved Sensitivity for Medium- and High-Molecular-Weight Proteins

NMR spin relaxation techniques that utilize relaxation interference phenomena (TROSY) enable chemical exchange processes to be characterized in high-molecular-weight proteins. A TROSY-selected (TS) approach for measuring off-resonance R1rho relaxation in the spin-locked rotating reference frame is developed using three principles: (i) deuteration of nonexchangeable 1H sites to minimize remote dipole-dipole interactions, (ii) selective excitation of the slowly relaxing 15N doublet component to obtain optimal initial conditions, and (iii) selective inversion of one of the 15N doublet components to suppress cross-relaxation during the spin-lock period. The method is validated using [90%-15N, 70%-2H] ubiquitin at 280 K. The TROSY-selected R1rho experiment enables characterization of backbone dynamics on the microsecond time scale in large proteins.

CEESY: characterizing the conformation of unobservable protein states

Protein conformations that are only marginally populated often play important roles as intermediate states in many processes such as ligand binding, enzyme catalysis, allostery, and protein folding. An NMR method is presented that can give valuable information about the structure of these "excited states" by measuring the relative position of exchanging excited- and ground-state resonances using a single 2D spectrum. This new approach can be applied to any nucleus, which will facilitate a complete structural characterization of these states.

Multiple-Site Exchange in Proteins Studied with a Suite of Six NMR Relaxation Dispersion Experiments: An Application to the Folding of a Fyn SH3 Domain Mutant

The three-site exchange folding reaction of an (15)N-labeled, highly deuterated Gly48Met mutant of the Fyn SH3 domain has been characterized at 25 degrees C using a suite of six CPMG-type relaxation dispersion experiments that measure exchange contributions to backbone (1)H and (15)N transverse relaxation rates in proteins. It is shown that this suite of experiments allows the extraction of all the parameters of this multisite exchange process in a robust manner, including chemical shift differences between exchanging states, from a data set recorded at only a single temperature. The populations of the exchanging folded, intermediate, and unfolded states that are fit are 94, 0.7, and 5%, respectively. Despite the small fraction of the intermediate, structural information is obtained for this state that is consistent with the picture of SH3 domain folding that has emerged from other studies. Taken together, the six dispersion experiments facilitate the complete reconstruction of (1)H-(15)N correlation spectra for the unfolded and intermediate states that are "invisible" in even the most sensitive of NMR experiments.

Iron responsive element RNA flexibility described by NMR and isotropic reorientational eigenmode dynamics

The first example of the application of reorientational eigenmode dynamics (RED) to RNA is shown here for the small and floppy Iron Responsive Element (IRE) RNA hairpin. Order parameters calculated for bases and riboses from a 12 ns molecular dynamics trajectory are compared to experimentally determined order parameters from 13C-1H NMR relaxation experiments, and shown to be in qualitative agreement. Given the small size of the IRE hairpin and its very flexible loop, isotropic RED (iRED) was also used to analyze the trajectory in order to describe its dynamic motions. iRED analysis shows that the global and internal dynamics of the IRE are not rigorously separable, which will result in inaccurate experimental order parameters. In addition, the iRED analysis described the many correlated motions that comprise the dynamics of the IRE RNA. The combined use of NMR relaxation, RED, and iRED provide a uniquely detailed description of IRE RNA dynamics.