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

NMR characterization and ligand binding site of the stem-loop 2 motif from the Delta variant of SARS-CoV-2

The stem-loop 2 motif (s2m) in SARS-CoV-2 (SCoV-2) is located in the 3'-UTR. Although s2m has been reported to display characteristics of a mobile genomic element that might lead to an evolutionary advantage, its function has remained unknown. The secondary structure of the original SCoV-2 RNA sequence (Wuhan-Hu-1) was determined by NMR in late 2020, delineating the base-pairing pattern and revealing substantial differences in secondary structure compared to SARS-CoV-1 (SCoV-1). The existence of a single G29742-A29756 mismatch in the upper stem of s2m leads to its destabilization and impedes a complete NMR analysis. With Delta, a variant of concern has evolved with one mutation compared to the original sequence that replaces G29742 by U29742. We show here that this mutation results in a more defined structure at ambient temperature accompanied by a rise in melting temperature. Consequently, we were able to identify >90% of the relevant NMR resonances using a combination of selective RNA labeling and filtered 2D NOESY as well as 4D NMR experiments. We present a comprehensive NMR analysis of the secondary structure, (sub)nanosecond dynamics, and ribose conformation of s2m Delta based on heteronuclear 13C NOE and T 1 measurements and ribose carbon chemical shift-derived canonical coordinates. We further show that the G29742U mutation in Delta has no influence on the druggability of s2m compared to the Wuhan-Hu-1 sequence. With the assignment at hand, we identify the flexible regions of s2m as the primary site for small molecule binding.

Dissecting Monomer-Dimer Equilibrium of an RNase P Protein Provides Insight Into the Synergistic Flexibility of 5' Leader Pre-tRNA Recognition

Ribonuclease P (RNase P) is a universal RNA-protein endonuclease that catalyzes 5' precursor-tRNA (ptRNA) processing. The RNase P RNA plays the catalytic role in ptRNA processing; however, the RNase P protein is required for catalysis in vivo and interacts with the 5' leader sequence. A single P RNA and a P protein form the functional RNase P holoenzyme yet dimeric forms of bacterial RNase P can interact with non-tRNA substrates and influence bacterial cell growth. Oligomeric forms of the P protein can also occur in vitro and occlude the 5' leader ptRNA binding interface, presenting a challenge in accurately defining the substrate recognition properties. To overcome this, concentration and temperature dependent NMR studies were performed on a thermostable RNase P protein from Thermatoga maritima. NMR relaxation (R₁, R₂), heteronuclear NOE, and diffusion ordered spectroscopy (DOSY) experiments were analyzed, identifying a monomeric species through the determination of the diffusion coefficients (D) and rotational correlation times (τ_c). Experimental diffusion coefficients and τ_c values for the predominant monomer (2.17 ± 0.36 * 10^-10 m²/s, τ _c = 5.3 ns) or dimer (1.87 ± 0.40* 10^-10 m²/s, τ _c = 9.7 ns) protein assemblies at 45°C correlate well with calculated diffusion coefficients derived from the crystallographic P protein structure (PDB 1NZ0). The identification of a monomeric P protein conformer from relaxation data and chemical shift information enabled us to gain novel insight into the structure of the P protein, highlighting a lack of structural convergence of the N-terminus (residues 1-14) in solution. We propose that the N-terminus of the bacterial P protein is partially disordered and adopts a stable conformation in the presence of RNA. In addition, we have determined the location of the 5' leader RNA in solution and measured the affinity of the 5' leader RNA-P protein interaction. We show that the monomer P protein interacts with RNA at the 5' leader binding cleft that was previously identified using X-ray crystallography. Data support a model where N-terminal protein flexibility is stabilized by holoenzyme formation and helps to accommodate the 5' leader region of ptRNA. Taken together, local structural changes of the P protein and the 5' leader RNA provide a means to obtain optimal substrate alignment and activation of the RNase P holoenzyme.

m6A minimally impacts the structure, dynamics, and Rev ARM binding properties of HIV-1 RRE stem IIB

N6-methyladenosine (m6A) is a ubiquitous RNA post-transcriptional modification found in coding as well as non-coding RNAs. m6A has also been found in viral RNAs where it is proposed to modulate host-pathogen interactions. Two m6A sites have been reported in the HIV-1 Rev response element (RRE) stem IIB, one of which was shown to enhance binding to the viral protein Rev and viral RNA export. However, because these m6A sites have not been observed in other studies mapping m6A in HIV-1 RNA, their significance remains to be firmly established. Here, using optical melting experiments, NMR spectroscopy, and in vitro binding assays, we show that m6A minimally impacts the stability, structure, and dynamics of RRE stem IIB as well as its binding affinity to the Rev arginine-rich-motif (ARM) in vitro. Our results indicate that if present in stem IIB, m6A is unlikely to substantially alter the conformational properties of the RNA. Our results add to a growing view that the impact of m6A on RNA depends on sequence context and Mg2+.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Structural insight into precursor tRNA processing by yeast ribonuclease P

Ribonuclease P (RNase P) is a universal ribozyme responsible for processing the 5'-leader of pre-transfer RNA (pre-tRNA). Here, we report the 3.5-angstrom cryo-electron microscopy structures of Saccharomyces cerevisiae RNase P alone and in complex with pre-tRNA^Phe The protein components form a hook-shaped architecture that wraps around the RNA and stabilizes RNase P into a "measuring device" with two fixed anchors that recognize the L-shaped pre-tRNA. A universally conserved uridine nucleobase and phosphate backbone in the catalytic center together with the scissile phosphate and the O3' leaving group of pre-tRNA jointly coordinate two catalytic magnesium ions. Binding of pre-tRNA induces a conformational change in the catalytic center that is required for catalysis. Moreover, simulation analysis suggests a two-metal-ion S_N2 reaction pathway of pre-tRNA cleavage. These results not only reveal the architecture of yeast RNase P but also provide a molecular basis of how the 5'-leader of pre-tRNA is processed by eukaryotic RNase P.

Inner-Sphere Coordination of Divalent Metal Ion with Nucleobase in Catalytic RNA

Identification of the function of metal ions and the RNA moieties, particularly nucleobases, that bind metal ions is important in RNA catalysis. Here we combine single-atom and abasic substitutions to probe functions of conserved nucleobases in ribonuclease P (RNase P). Structural and biophysical studies of bacterial RNase P propose direct coordination of metal ions by the nucleobases of conserved uridine and guanosine in helix P4 of the RNA subunit (P RNA). To biochemically probe the function of metal ion interactions, we substituted the universally conserved bulged uridine (U51) in the P4 helix of circularly permuted Bacillus subtilis P RNA with 4-thiouridine, 4-deoxyuridine, and abasic modifications and G378/379 with 2-aminopurine, N7-deazaguanosine, and 6-thioguanosine. The functional group modifications of U51 decrease RNase P-catalyzed phosphodiester bond cleavage 16- to 23-fold, as measured by the single-turnover cleavage rate constant. The activity of the 4-thiouridine RNase P is partially rescued by addition of Cd(II) or Mn(II) ions. This is the first time a metal-rescue experiment provides evidence for inner-sphere divalent metal ion coordination with a nucleobase. Modifications of G379 modestly decrease the cleavage activity of RNase P, suggesting outer-sphere coordination of O6 on G379 to a metal ion. These data provide biochemical evidence for catalytically important interactions of the P4 helix of P RNA with metal ions, demonstrating that the bulged uridine coordinates at least one catalytic metal ion through an inner-sphere interaction. The combination of single-atom and abasic nucleotide substitutions provides a powerful strategy to probe functions of conserved nucleobases in large RNAs.

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.

The Diversity of Ribonuclease P: Protein and RNA Catalysts with Analogous Biological Functions

Ribonuclease P (RNase P) is an essential endonuclease responsible for catalyzing 5' end maturation in precursor transfer RNAs. Since its discovery in the 1970s, RNase P enzymes have been identified and studied throughout the three domains of life. Interestingly, RNase P is either RNA-based, with a catalytic RNA subunit, or a protein-only (PRORP) enzyme with differential evolutionary distribution. The available structural data, including the active site data, provides insight into catalysis and substrate recognition. The hydrolytic and kinetic mechanisms of the two forms of RNase P enzymes are similar, yet features unique to the RNA-based and PRORP enzymes are consistent with different evolutionary origins. The various RNase P enzymes, in addition to their primary role in tRNA 5' maturation, catalyze cleavage of a variety of alternative substrates, indicating a diversification of RNase P function in vivo. The review concludes with a discussion of recent advances and interesting research directions in the field.

Dynamic motions of the HIV-1 frameshift site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDOint) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDOint = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.

Structural dynamics of a single-stranded RNA-helix junction using NMR

Many regulatory RNAs contain long single strands (ssRNA) that adjoin secondary structural elements. Here, we use NMR spectroscopy to study the dynamic properties of a 12-nucleotide (nt) ssRNA tail derived from the prequeuosine riboswitch linked to the 3' end of a 48-nt hairpin. Analysis of chemical shifts, NOE connectivity, (13)C spin relaxation, and residual dipolar coupling data suggests that the first two residues (A25 and U26) in the ssRNA tail stack onto the adjacent helix and assume an ordered conformation. The following U26-A27 step marks the beginning of an A6-tract and forms an acute pivot point for substantial motions within the tail, which increase toward the terminal end. Despite substantial internal motions, the ssRNA tail adopts, on average, an A-form helical conformation that is coaxial with the helix. Our results reveal a surprising degree of structural and dynamic complexity at the ssRNA-helix junction, which involves a fine balance between order and disorder that may facilitate efficient pseudoknot formation on ligand recognition.

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

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.

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields.

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.

Of proteins and RNA: the RNase P/MRP family

Nuclear ribonuclease (RNase) P is a ubiquitous essential ribonucleoprotein complex, one of only two known RNA-based enzymes found in all three domains of life. The RNA component is the catalytic moiety of RNases P across all phylogenetic domains; it contains a well-conserved core, whereas peripheral structural elements are diverse. RNA components of eukaryotic RNases P tend to be less complex than their bacterial counterparts, a simplification that is accompanied by a dramatic reduction of their catalytic ability in the absence of protein. The size and complexity of the protein moieties increase dramatically from bacterial to archaeal to eukaryotic enzymes, apparently reflecting the delegation of some structural functions from RNA to proteins and, perhaps, in response to the increased complexity of the cellular environment in the more evolutionarily advanced organisms; the reasons for the increased dependence on proteins are not clear. We review current information on RNase P and the closely related universal eukaryotic enzyme RNase MRP, focusing on their functions and structural organization.

A divalent cation stabilizes the active conformation of the B. subtilis RNase P x pre-tRNA complex: a role for an inner-sphere metal ion in RNase P

Metal ions interact with RNA to enhance folding, stabilize structure, and, in some cases, facilitate catalysis. Assigning functional roles to specifically bound metal ions presents a major challenge in analyzing the catalytic mechanisms of ribozymes. Bacillus subtilis ribonuclease P (RNase P), composed of a catalytically active RNA subunit (PRNA) and a small protein subunit (P protein), catalyzes the 5'-end maturation of precursor tRNAs (pre-tRNAs). Inner-sphere coordination of divalent metal ions to PRNA is essential for catalytic activity but not for the formation of the RNase P x pre-tRNA (enzyme-substrate, ES) complex. Previous studies have demonstrated that this ES complex undergoes an essential conformational change (to the ES* conformer) before the cleavage step. Here, we show that the ES* conformer is stabilized by a high-affinity divalent cation capable of inner-sphere coordination, such as Ca(II) or Mg(II). Additionally, a second, lower-affinity Mg(II) activates cleavage catalyzed by RNase P. Structural changes that occur upon binding Ca(II) to the ES complex were determined by time-resolved Förster resonance energy transfer measurements of the distances between donor-acceptor fluorophores introduced at specific locations on the P protein and pre-tRNA 5' leader. These data demonstrate that the 5' leader of pre-tRNA moves 4 to 6 A closer to the PRNA x P protein interface during the ES-to-ES* transition and suggest that the metal-dependent conformational change reorganizes the bound substrate in the active site to form a catalytically competent ES* complex.

NMR and XAS reveal an inner-sphere metal binding site in the P4 helix of the metallo-ribozyme ribonuclease P

Functionally critical metals interact with RNA through complex coordination schemes that are currently difficult to visualize at the atomic level under solution conditions. Here, we report a new approach that combines NMR and XAS to resolve and characterize metal binding in the most highly conserved P4 helix of ribonuclease P (RNase P), the ribonucleoprotein that catalyzes the divalent metal ion-dependent maturation of the 5' end of precursor tRNA. Extended X-ray absorption fine structure (EXAFS) spectroscopy reveals that the Zn(2+) bound to a P4 helix mimic is six-coordinate, with an average Zn-O/N bond distance of 2.08 A. The EXAFS data also show intense outer-shell scattering indicating that the zinc ion has inner-shell interactions with one or more RNA ligands. NMR Mn(2+) paramagnetic line broadening experiments reveal strong metal localization at residues corresponding to G378 and G379 in B. subtilis RNase P. A new "metal cocktail" chemical shift perturbation strategy involving titrations with , Zn(2+), and confirm an inner-sphere metal interaction with residues G378 and G379. These studies present a unique picture of how metals coordinate to the putative RNase P active site in solution, and shed light on the environment of an essential metal ion in RNase P. Our experimental approach presents a general method for identifying and characterizing inner-sphere metal ion binding sites in RNA in solution.

Solution NMR studies of Chlorella virus DNA ligase-adenylate

DNA ligases are essential guardians of genome integrity by virtue of their ability to recognize and seal 3'-OH/5'-phosphate nicks in duplex DNA. The substrate binding and three chemical steps of the ligation pathway are coupled to global and local changes in ligase structure, involving both massive protein domain movements and subtle remodeling of atomic contacts in the active site. Here we applied solution NMR spectroscopy to study the conformational dynamics of the Chlorella virus DNA ligase (ChVLig), a minimized eukaryal ATP-dependent ligase consisting of nucleotidyltransferase, OB, and latch domains. Our analysis of backbone (15)N spin relaxation and (15)N,(1)H residual dipolar couplings of the covalent ChVLig-AMP intermediate revealed conformational sampling on fast (picosecond to nanosecond) and slow timescales (microsecond to millisecond), indicative of interdomain and intradomain flexibility. We identified local and global changes in ChVLig-AMP structure and dynamics induced by phosphate. In particular, the chemical shift perturbations elicited by phosphate were clustered in the peptide motifs that comprise the active site. We hypothesize that phosphate anion mimics some of the conformational transitions that occur when ligase-adenylate interacts with the nick 5'-phosphate.

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.

Preparation, resonance assignment, and preliminary dynamics characterization of residue specific 13C/15N-labeled elongated DNA for the study of sequence-directed dynamics by NMR

DNA is a highly flexible molecule that undergoes functionally important structural transitions in response to external cellular stimuli. Atomic level spin relaxation NMR studies of DNA dynamics have been limited to short duplexes in which sensitivity to biologically relevant fluctuations occurring at nanosecond timescales is often inadequate. Here, we introduce a method for preparing residue-specific (13)C/(15)N-labeled elongated DNA along with a strategy for establishing resonance assignments and apply the approach to probe fast inter-helical bending motions induced by an adenine tract. Preliminary results suggest the presence of elevated A-tract independent end-fraying internal motions occurring at nanosecond timescales, which evade detection in short DNA constructs and that penetrate deep (7 bp) within the DNA helix and gradually fade away towards the helix interior.

Extensive backbone dynamics in the GCAA RNA tetraloop analyzed using 13C NMR spin relaxation and specific isotope labeling

Conformational dynamics play a key role in the properties and functions of proteins and nucleic acids. Heteronuclear NMR spin relaxation is a uniquely powerful site-specific probe of dynamics in proteins and has found increasing applications to nucleotide base side chains and anomeric sites in RNA. Applications to the nucleic acid ribose backbone, however, have been hampered by strong magnetic coupling among ring carbons in uniformly 13C-labeled samples. In this work, we apply a recently developed, metabolically directed isotope labeling scheme that places 13C with high efficiency and specificity at the nucleotide ribose C2' and C4' sites. We take advantage of this scheme to explore backbone dynamics in the well-studied GCAA RNA tetraloop. Using a combination of CPMG (Carr-Purcell-Meiboom-Gill) and R(1rho) relaxation dispersion spectroscopy to explore exchange processes on the microsecond to millisecond time scale, we find an extensive pattern of dynamic transitions connecting a set of relatively well-defined conformations. In many cases, the observed transitions appear to be linked to C3'-endo/C2'-endo sugar pucker transitions of the corresponding nucleotides, and may also be correlated across multiple nucleotides within the tetraloop. These results demonstrate the power of NMR spin relaxation based on alternate-site isotope labeling to open a new window into the dynamic properties of ribose backbone groups in RNA.

Constructing atomic-resolution RNA structural ensembles using MD and motionally decoupled NMR RDCs

A broad structural landscape often needs to be characterized in order to fully understand how regulatory RNAs perform their biological functions at the atomic level. We present a protocol for visualizing thermally accessible RNA conformations at atomic-resolution and with timescales extending up to milliseconds. The protocol combines molecular dynamics (MD) simulations with experimental residual dipolar couplings (RDCs) measured in partially aligned (13)C/(15)N isotopically enriched elongated RNA samples. The structural ensembles generated in this manner provide insights into RNA dynamics and its role in functionally important transitions.

Mapping metal-binding sites in the catalytic domain of bacterial RNase P RNA

Ribonuclease P (RNase P) is a ribonucleoprotein enzyme that contains a universally conserved, catalytically active RNA component. RNase P RNA requires divalent metal ions for folding, substrate binding, and catalysis. Despite recent advances in understanding the structure of RNase P RNA, no comprehensive analysis of metal-binding sites has been reported, in part due to the poor crystallization properties of this large RNA. We have developed an abbreviated yet still catalytic construct, Bst P7Delta RNA, which contains the catalytic domain of the bacterial RNase P RNA and has improved crystallization properties. We use this mutant RNA as well as the native RNA to map metal-binding sites in the catalytic core of the bacterial RNase P RNA, by anomalous scattering in diffraction analysis. The results provide insight into the interplay between RNA structure and focalization of metal ions, and a structural basis for some previous biochemical observations with RNase P. We use electrostatic calculations to extract the potential functional significance of these metal-binding sites with respect to binding Mg(2+). The results suggest that with at least one important exception of specific binding, these sites mainly map areas of diffuse association of magnesium ions.

Bacterial RNase P RNA is a drug target for aminoglycoside-arginine conjugates

The ribonuclease P (RNase P) holoenzymes are RNPs composed of RNase P RNA (PRNA) and a variable number of P protein subunits. Primary differences in structure and function between bacterial and eukaryotic RNase P and its indispensability for cell viability make the bacterial enzyme an attractive drug target. On the basis of our previous studies, aminoglycoside-arginine conjugates (AACs) bind to HIV-1 TAR and Rev responsive element (RRE) RNAs significantly more efficiently than neomycin B. Their specific inhibition of bacterial rRNA as well as the findings that the hexa-arginine neomycin derivative (NeoR6) is 500-fold more potent than neomycin B in inhibiting bacterial RNase P, led us to explore the structure-function relationships of AACs in comparison to a new set of aminoglycoside-polyarginine conjugates (APACs). We here present predicted binding modes of AACs and APACs to PRNA. We used a multistep docking approach comprising rigid docking full scans and final refinement of the obtained complexes. Our docking results suggest three possible mechanisms of RNase P inhibition by AACs and APACs: competition with the P protein and pre-tRNA on binding to P1-P4 multihelix junction and to J19/4 region (probably including displacement of Mg2+ ions from the P4 helix) of PRNA; competition with Mg2+ ions near the P15 loop; and competition with the P protein and/or pre-tRNA near the P15 helix and interfering with interactions between the P protein and pre-tRNA at this region. The APACs revealed about 10-fold lower intermolecular energy than AACs, indicating stronger interactions of APACs than AACs with PRNA.

Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis

The HIV-1 transactivation response element (TAR) RNA binds a variety of proteins and is a target for developing anti-HIV therapies. TAR has two primary binding sites: a UCU bulge and a CUGGGA apical loop. We used NMR residual dipolar couplings, carbon spin relaxation (R(1) and R(2)), and relaxation dispersion (R(1rho)) in conjunction with molecular dynamics and mutagenesis to characterize the dynamics of the TAR apical loop and investigate previously proposed long-range interactions with the distant bulge. Replacement of the wild-type apical loop with a UUCG loop did not significantly affect the structural dynamics at the bulge, indicating that the apical loop and the bulge act largely as independent dynamical recognition centers. The apical loop undergoes complex dynamics at multiple timescales that are likely important for adaptive recognition: U31 and G33 undergo limited motions, G32 is highly flexible at picosecond-nanosecond timescales, and G34 and C30 form a dynamic Watson-Crick basepair in which G34 and A35 undergo a slow (approximately 30 mus) likely concerted looping in and out motion, with A35 also undergoing large amplitude motions at picosecond-nanosecond timescales. Our study highlights the power of combining NMR, molecular dynamics, and mutagenesis in characterizing RNA dynamics.

Unwinding RNA's secrets: advances in the biology, physics, and modeling of complex RNAs

The rapid development of our understanding of the diverse biological roles fulfilled by non-coding RNA has motivated interest in the basic macromolecular behavior, structure, and function of RNA. We focus on two areas in the behavior of complex RNAs. First, we present advances in the understanding of how RNA folding is accomplished in vivo by presenting a mechanism for the action of DEAD-box proteins. Members of this family are intimately associated with almost all cellular processes involving RNA, mediating RNA structural rearrangements and chaperoning their folding. Next, we focus on advances in understanding, and characterizing the basic biophysical forces that govern the folding of complex RNAs. Ultimately we expect that a confluence and synergy between these approaches will lead to profound understanding of RNA and its biology.

Importance of RNA-protein interactions in bacterial ribonuclease P structure and catalysis

Ribonuclease P (RNase P) is a ribonucleoprotein (RNP) complex that catalyzes the metal-dependent maturation of the 5' end of precursor tRNAs (pre-tRNAs) in all organisms. RNase P is comprised of a catalytic RNA (P RNA), and at least one essential protein (P protein). Although P RNA is the catalytic subunit of the enzyme and is active in the absence of P protein under high salt concentrations in vitro, the protein is still required for enzyme activity in vivo. Therefore, the function of the P protein and how it interacts with both P RNA and pre-tRNA have been the focus of much ongoing research. RNA-protein interactions in RNase P serve a number of critical roles in the RNP including stabilizing the structure, and enhancing the affinity for substrates and metal ions. This review examines the role of RNA-protein interactions in bacterial RNase P from both structural and mechanistic perspectives.

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

Resolving fast and slow motions in the internal loop containing stem-loop 1 of HIV-1 that are modulated by Mg2+ binding: role in the kissing-duplex structural transition

Stem loop 1 (SL1) is a highly conserved hairpin in the 5'-leader of the human immunodeficiency virus type I that forms a metastable kissing dimer that is converted during viral maturation into a stable duplex with the aid of the nucleocapsid (NC) protein. SL1 contains a highly conserved internal loop that promotes the kissing-duplex transition by a mechanism that remains poorly understood. Using NMR, we characterized internal motions induced by the internal loop in an SL1 monomer that may promote the kissing-duplex transition. This includes micro-to-millisecond secondary structural transitions that cause partial melting of three base-pairs above the internal loop making them key nucleation sites for exchanging strands and nanosecond rigid-body stem motions that can help bring strands into spatial register. We show that while Mg2+ binds to the internal loop and arrests these internal motions, it preserves and/or activates local mobility at internal loop residues G272 and G273 which are implicated in NC binding. By stabilizing SL1 without compromising the accessibility of G272 and G273 for NC binding, Mg2+ may increase the dependence of the kissing-duplex transition on NC binding thus preventing spontaneous transitions from taking place and ensuring that viral RNA and protein maturation occur in concert.