Direct identification of NH...N hydrogen bonds in non-canonical base pairs of RNA by NMR spectroscopy

Direct Detection of Hydrogen Bonds in Supramolecular Systems Using 1H-15N Heteronuclear Multiple Quantum Coherence Spectroscopy

Hydrogen-bonded supramolecular systems are usually characterized in solution through analysis of NMR data such as complexation-induced shifts and nuclear Overhauser effects (nOe). Routine direct detection of hydrogen bonding particularly in multicomponent mixtures, even with the aid of 2D NMR experiments for full assignment, is more challenging. We describe an elementary rapid ¹H-¹⁵N HMQC NMR experiment which addresses these challenges without the need for complex pulse sequences. Under readily accessible conditions (243/263 K, 50 mM solutions) and natural ¹⁵N abundance, unambiguous assignment of ¹⁵N resonances facilitates direct detection of intra- and intermolecular hydrogen bonds in mechanically interlocked structures and quadruply hydrogen-bonded dimers─of dialkylaminoureidopyrimidinones, ureidopyrimidinones, and diamidonaphthyridines─in single or multicomponent mixtures to establish tautomeric configuration, conformation, and, to resolve self-sorted speciation.

RNA modifications stabilize the tertiary structure of tRNAfMet by locally increasing conformational dynamics

A plethora of modified nucleotides extends the chemical and conformational space for natural occurring RNAs. tRNAs constitute the class of RNAs with the highest modification rate. The extensive modification modulates their overall stability, the fidelity and efficiency of translation. However, the impact of nucleotide modifications on the local structural dynamics is not well characterized. Here we show that the incorporation of the modified nucleotides in tRNAfMet from Escherichia coli leads to an increase in the local conformational dynamics, ultimately resulting in the stabilization of the overall tertiary structure. Through analysis of the local dynamics by NMR spectroscopic methods we find that, although the overall thermal stability of the tRNA is higher for the modified molecule, the conformational fluctuations on the local level are increased in comparison to an unmodified tRNA. In consequence, the melting of individual base pairs in the unmodified tRNA is determined by high entropic penalties compared to the modified. Further, we find that the modifications lead to a stabilization of long-range interactions harmonizing the stability of the tRNA's secondary and tertiary structure. Our results demonstrate that the increase in chemical space through introduction of modifications enables the population of otherwise inaccessible conformational substates.

1H, 13C and 15N chemical shift assignment of the stem-loop 5a from the 5'-UTR of SARS-CoV-2

The SARS-CoV-2 (SCoV-2) virus is the causative agent of the ongoing COVID-19 pandemic. It contains a positive sense single-stranded RNA genome and belongs to the genus of Betacoronaviruses. The 5'- and 3'-genomic ends of the 30 kb SCoV-2 genome are potential antiviral drug targets. Major parts of these sequences are highly conserved among Betacoronaviruses and contain cis-acting RNA elements that affect RNA translation and replication. The 31 nucleotide (nt) long highly conserved stem-loop 5a (SL5a) is located within the 5'-untranslated region (5'-UTR) important for viral replication. SL5a features a U-rich asymmetric bulge and is capped with a 5'-UUUCGU-3' hexaloop, which is also found in stem-loop 5b (SL5b). We herein report the extensive ¹H, ¹³C and ¹⁵N resonance assignment of SL5a as basis for in-depth structural studies by solution NMR spectroscopy.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.

The structure of the SAM/SAH-binding riboswitch

S-adenosylmethionine (SAM) is a central metabolite since it is used as a methyl group donor in many different biochemical reactions. Many bacteria control intracellular SAM concentrations using riboswitch-based mechanisms. A number of structurally different riboswitch families specifically bind to SAM and mainly regulate the transcription or the translation of SAM-biosynthetic enzymes. In addition, a highly specific riboswitch class recognizes S-adenosylhomocysteine (SAH)—the product of SAM-dependent methyl group transfer reactions—and regulates enzymes responsible for SAH hydrolysis. High-resolution structures are available for many of these riboswitch classes and illustrate how they discriminate between the two structurally similar ligands SAM and SAH. The so-called SAM/SAH riboswitch class binds both ligands with similar affinities and is structurally not yet characterized. Here, we present a high-resolution nuclear magnetic resonance structure of a member of the SAM/SAH-riboswitch class in complex with SAH. Ligand binding induces pseudoknot formation and sequestration of the ribosome binding site. Thus, the SAM/SAH-riboswitches are translational ‘OFF’-switches. Our results establish a structural basis for the unusual bispecificity of this riboswitch class. In conjunction with genomic data our structure suggests that the SAM/SAH-riboswitches might be an evolutionary late invention and not a remnant of a primordial RNA-world as suggested for other riboswitches.

Viroid research and its significance for RNA technology and basic biochemistry

Viroids were described 47 years ago as the smallest RNA molecules capable of infecting plants and autonomously self-replicating without an encoded protein. Work on viroids initiated the development of a number of innovative methods. Novel chromatographic and gelelectrophoretic methods were developed for the purification and characterization of viroids; these methods were later used in molecular biology, gene technology and in prion research. Theoretical and experimental studies of RNA folding demonstrated the general biological importance of metastable structures, and nuclear magnetic resonance spectroscopy of viroid RNA showed the partially covalent nature of hydrogen bonds in biological macromolecules. RNA biochemistry and molecular biology profited from viroid research, such as in the detection of RNA as template of DNA-dependent polymerases and in mechanisms of gene silencing. Viroids, the first circular RNA detected in nature, are important for studies on the much wider spectrum of circular RNAs and other non-coding RNAs.

NMR resonance assignments for the SAM/SAH-binding riboswitch RNA bound to S-adenosylhomocysteine

Riboswitches are structured RNA elements in the 5'-untranslated regions of bacterial mRNAs that are able to control the transcription or translation of these mRNAs in response to the specific binding of small molecules such as certain metabolites. Riboswitches that bind with high specificity to either S-adenosylmethionine (SAM) or S-adenosylhomocysteine (SAH) are widespread in bacteria. Based on differences in secondary structure and sequence these riboswitches can be grouped into a number of distinct classes. X-ray structures for riboswitch RNAs in complex with SAM or SAH established a structural basis for understanding ligand recognition and discrimination in many of these riboswitch classes. One class of riboswitches-the so-called SAM/SAH riboswitch class-binds SAM and SAH with similar affinity. However, this class of riboswitches is structurally not yet characterized and the structural basis for its unusual bispecificity is not established. In order to understand the ligand recognition mode that enables this riboswitch to bind both SAM and SAH with similar affinities, we are currently determining its structure in complex with SAH using NMR spectroscopy. Here, we present the NMR resonance assignment of the SAM/SAH binding riboswitch (env9b) in complex with SAH as a prerequisite for a solution NMR-based high-resolution structure determination.

Adenine protonation enables cyclic-di-GMP binding to cyclic-GAMP sensing riboswitches

In certain structural or functional contexts, RNA structures can contain protonated nucleotides. However, a direct role for stably protonated nucleotides in ligand binding and ligand recognition has not yet been demonstrated unambiguously. Previous X-ray structures of c-GAMP binding riboswitch aptamer domains in complex with their near-cognate ligand c-di-GMP suggest that an adenine of the riboswitch either forms two hydrogen bonds to a G nucleotide of the ligand in the unusual enol tautomeric form or that the adenine in its N1 protonated form binds the G nucleotide of the ligand in its canonical keto tautomeric state. By using NMR spectroscopy we demonstrate that the c-GAMP riboswitches bind c-di-GMP using a stably protonated adenine in the ligand binding pocket. Thereby, we provide novel insights into the putative biological functions of protonated nucleotides in RNA, which in this case influence the ligand selectivity in a riboswitch.

Structure of the 30 kDa HIV-1 RNA Dimerization Signal by a Hybrid Cryo-EM, NMR, and Molecular Dynamics Approach

Cryoelectron microscopy (cryo-EM) and nuclear magnetic resonance (NMR) spectroscopy are routinely used to determine structures of macromolecules with molecular weights over 65 and under 25 kDa, respectively. We combined these techniques to study a 30 kDa HIV-1 dimer initiation site RNA ([DIS]2; 47 nt/strand). A 9 Å cryo-EM map clearly shows major groove features of the double helix and a right-handed superhelical twist. Simulated cryo-EM maps generated from time-averaged molecular dynamics trajectories (10 ns) exhibited levels of detail similar to those in the experimental maps, suggesting internal structural flexibility limits the cryo-EM resolution. Simultaneous inclusion of the cryo-EM map and 2H-edited NMR-derived distance restraints during structure refinement generates a structure consistent with both datasets and supporting a flipped-out base within a conserved purine-rich bulge. Our findings demonstrate the power of combining global and local structural information from these techniques for structure determination of modest-sized RNAs.

Structure modeling of RNA using sparse NMR constraints

RNAs fold into distinct molecular conformations that are often essential for their functions. Accurate structure modeling of complex RNA motifs, including ubiquitous non-canonical base pairs and pseudoknots, remains a challenge. Here, we present an NMR-guided all-atom discrete molecular dynamics (DMD) platform, iFoldNMR, for rapid and accurate structure modeling of complex RNAs. We show that sparse distance constraints from imino resonances, which can be readily obtained from routine NMR experiments and easier to compile than laborious assignments of non-solvent-exchangeable protons, are sufficient to direct a DMD search for low-energy RNA conformers. Benchmarking on a set of RNAs with complex folds spanning up to 56 nucleotides in length yields structural models that recapitulate experimentally determined structures with all-heavy-atom RMSDs ranging from 2.4 to 6.5 Å. This platform represents an efficient approach for high-throughput RNA structure modeling and will facilitate analysis of diverse, newly discovered functional RNAs.

NMR experiments for the rapid identification of P=O···H-X type hydrogen bonds in nucleic acids

Hydrogen bonds involving the backbone phosphate groups occur with high frequency in functional RNA molecules. They are often found in well-characterized tertiary structural motifs presenting powerful probes for the rapid identification of these motifs for structure elucidation purposes. We have shown recently that stable hydrogen bonds to the phosphate backbone can in principle be detected by relatively simple NMR-experiments, providing the identity of both the donor hydrogen and the acceptor phosphorous within the same experiment (Duchardt-Ferner et al., Angew Chem Int Ed Engl 50:7927-7930, 2011). However, for imino and hydroxyl hydrogen bond donor groups rapidly exchanging with the solvent as well as amino groups broadened by conformational exchange experimental sensitivity is severely hampered by extensive line broadening. Here, we present improved methods for the rapid identification of hydrogen bonds to phosphate groups in nucleic acids by NMR. The introduction of the SOFAST technique into ¹H,³¹P-correlation experiments as well as a BEST-HNP experiment exploiting ^3hJ_N,P rather than ^2hJ_H,P coupling constants enables the rapid and sensitive identification of these hydrogen bonds in RNA. The experiments are applicable for larger RNAs (up to ~ 100-nt), for donor groups influenced by conformational exchange processes such as amino groups and for hydrogen bonds with rather labile hydrogens such as 2'-OH groups as well as for moderate sample concentrations. Interestingly, the size of the through-hydrogen bond scalar coupling constants depends not only on the type of the donor group but also on the structural context. The largest coupling constants were measured for hydrogen bonds involving the imino groups of protonated cytosine nucleotides as donors.

NMR resonance assignments for the tetramethylrhodamine binding RNA aptamer 3 in complex with the ligand 5-carboxy-tetramethylrhodamine

RNA aptamers are used in a wide range of biotechnological or biomedical applications. In many cases the high resolution structures of these aptamers in their ligand-complexes have revealed fundamental aspects of RNA folding and RNA small molecule interactions. Fluorescent RNA-ligand complexes in particular find applications as optical sensors or as endogenous fluorescent tags for RNA tracking in vivo. Structures of RNA aptamers and aptamer ligand complexes constitute the starting point for rational function directed optimization approaches. Here, we present the NMR resonance assignment of an RNA aptamer binding to the fluorescent ligand tetramethylrhodamine (TMR) in complex with the ligand 5-carboxy-tetramethylrhodamine (5-TAMRA) as a starting point for a high-resolution structure determination using NMR spectroscopy in solution.

The synthesis of 15N(7)-Hoogsteen face-labeled adenosine phosphoramidite for solid-phase RNA synthesis

ABSTRACT

We have developed an efficient route for the synthesis of 15N(7)-labeled adenosine as phosphoramidite building block for site- and atom-specific incorporation into RNA by automated solid-phase synthesis. Such labeled RNA is required for the evaluation of selected non-canonical base pair interactions in folded RNA using NMR spectroscopic methods.

GRAPHICAL ABSTRACT

An intermolecular G-quadruplex as the basis for GTP recognition in the class V-GTP aptamer

Many naturally occurring or artificially created RNAs are capable of binding to guanine or guanine derivatives with high affinity and selectivity. They bind their ligands using very different recognition modes involving a diverse set of hydrogen bonding and stacking interactions. Apparently, the potential structural diversity for guanine, guanosine, and guanine nucleotide binding motifs is far from being fully explored. Szostak and coworkers have derived a large set of different GTP-binding aptamer families differing widely in sequence, secondary structure, and ligand specificity. The so-called class V-GTP aptamer from this set binds GTP with very high affinity and has a complex secondary structure. Here we use solution NMR spectroscopy to demonstrate that the class V aptamer binds GTP through the formation of an intermolecular two-layered G-quadruplex structure that directly incorporates the ligand and folds only upon ligand addition. Ligand binding and G-quadruplex formation depend strongly on the identity of monovalent cations present with a clear preference for potassium ions. GTP binding through direct insertion into an intermolecular G-quadruplex is a previously unobserved structural variation for ligand-binding RNA motifs and rationalizes the previously observed specificity pattern of the class V aptamer for GTP analogs.

NMR resonance assignments for the class II GTP binding RNA aptamer in complex with GTP

The structures of RNA-aptamer-ligand complexes solved in the last two decades were instrumental in realizing the amazing potential of RNA for forming complex tertiary structures and for molecular recognition of small molecules. For GTP as ligand the sequences and secondary structures for multiple families of aptamers were reported which differ widely in their structural complexity, ligand affinity and ligand functional groups involved in RNA-binding. However, for only one of these families the structure of the GTP-RNA complex was solved. In order to gain further insights into the variability of ligand recognition modes we are currently determining the structure of another GTP-aptamer--the so-called class II aptamer--bound to GTP using NMR-spectroscopy in solution. As a prerequisite for a full structure determination, we report here (1)H, (13)C, (15)N and partial (31)P-NMR resonance assignments for the class II GTP-aptamer bound to GTP.

NMR Methods for Characterization of RNA Secondary Structure

Knowledge of RNA secondary structure is often sufficient to identify relationships between the structure of RNA and processing pathways, and the design of therapeutics. Nuclear magnetic resonance (NMR) can identify types of nucleotide base pairs and the sequence, thus limiting possible secondary structures. Because NMR experiments, like chemical mapping, are performed in solution, not in single crystals, experiments can be initiated as soon as the biomolecule is expressed and purified. This chapter summarizes NMR methods that permit rapid identification of RNA secondary structure, information that can be used as supplements to chemical mapping, and/or as preliminary steps required for 3D structure determination. The primary aim is to provide guidelines to enable a researcher with minimal knowledge of NMR to quickly extract secondary structure information from basic datasets. Instrumental and sample considerations that can maximize data quality are discussed along with some details for optimal data acquisition and processing parameters. Approaches for identifying base pair types in both unlabeled and isotopically labeled RNA are covered. Common problems, such as missing signals and overlaps, and approaches to address them are considered. Programs under development for merging NMR data with structure prediction algorithms are briefly discussed.

A conditional random fields method for RNA sequence-structure relationship modeling and conformation sampling

Accurate tertiary structures are very important for the functional study of non-coding RNA molecules. However, predicting RNA tertiary structures is extremely challenging, because of a large conformation space to be explored and lack of an accurate scoring function differentiating the native structure from decoys. The fragment-based conformation sampling method (e.g. FARNA) bears shortcomings that the limited size of a fragment library makes it infeasible to represent all possible conformations well. A recent dynamic Bayesian network method, BARNACLE, overcomes the issue of fragment assembly. In addition, neither of these methods makes use of sequence information in sampling conformations. Here, we present a new probabilistic graphical model, conditional random fields (CRFs), to model RNA sequence–structure relationship, which enables us to accurately estimate the probability of an RNA conformation from sequence. Coupled with a novel tree-guided sampling scheme, our CRF model is then applied to RNA conformation sampling. Experimental results show that our CRF method can model RNA sequence–structure relationship well and sequence information is important for conformation sampling. Our method, named as TreeFolder, generates a much higher percentage of native-like decoys than FARNA and BARNACLE, although we use the same simple energy function as BARNACLE.

Contact:zywang@ttic.edu; j3xu@ttic.edu

Supplementary Information:Supplementary data are available at Bioinformatics online.

Rapid global structure determination of large RNA and RNA complexes using NMR and small-angle X-ray scattering

Among the greatest advances in biology today are the discoveries of various roles played by RNA in biological processes. However, despite significant advances in RNA structure determination using X-ray crystallography [1] and solution NMR [2-4], the number of bona fide RNA structures is very limited, in comparison with the growing number of known functional RNAs. This is because of great difficulty in growing crystals or/and obtaining phase information, and severe size constraints on structure determination by solution NMR spectroscopy. Clearly, there is an acute need for new methodologies for RNA structure determination. The prevailing approach for structure determination of RNA in solution is a "bottom-up" approach that was basically transplanted from the approach used for determining protein structures, despite vast differences in both structural features and chemical compositions between these two types of biomacromolecules. In this chapter, we describe a new method, which has been reported recently, for rapid global structure determination of RNAs using solution-based NMR spectroscopy and small-angle X-ray scattering. The method treats duplexes as major building blocks of RNA structures. By determining the global orientations of the duplexes and the overall shape, the global structure of an RNA can be constructed and further regularized using Xplor-NIH. The utility of the method was demonstrated in global structure determination of two RNAs, a 71-nt and 102-nt RNAs with an estimated backbone RMSD ∼3.0Å. The global structure opens door to high-resolution structure determination in solution.

Structures of HIV TAR RNA-ligand complexes reveal higher binding stoichiometries

Target TAR by NMR: Tripeptides containing arginines as terminal residues and non-natural amino acids as central residues are good leads for drug design to target the HIV trans-activation response element (TAR). The structural characterization of the RNA-ligand complex by NMR spectroscopy reveals two specific binding sites that are located at bulge residue U23 and around the pyrimidine-stretch U40-C41-U42 directly adjacent to the bulge.

A method for helical RNA global structure determination in solution using small-angle x-ray scattering and NMR measurements

We report a "top-down" method that uses mainly duplexes' global orientations and overall molecular dimension and shape restraints, which were extracted from experimental NMR and small-angle X-ray scattering data, respectively, to determine global architectures of RNA molecules consisting of mostly A-form-like duplexes. The method is implemented in the G2G (from global measurement to global structure) toolkit of programs. We demonstrate the efficiency and accuracy of the method by determining the global structure of a 71-nt RNA using experimental data. The backbone root-mean-square deviation of the ensemble of the calculated global structures relative to the X-ray crystal structure is 3.0+/-0.3 A using the experimental data and is only 2.5+/-0.2 A for the three duplexes that were orientation restrained during the calculation. The global structure simplifies interpretation of multidimensional nuclear Overhauser spectra for high-resolution structure determination. The potential general application of the method for RNA structure determination is discussed.

Longitudinal-relaxation-enhanced NMR experiments for the study of nucleic acids in solution

Atomic-resolution information on the structure and dynamics of nucleic acids is essential for a better understanding of the mechanistic basis of many cellular processes. NMR spectroscopy is a powerful method for studying the structure and dynamics of nucleic acids; however, solution NMR studies are currently limited to relatively small nucleic acids at high concentrations. Thus, technological and methodological improvements that increase the experimental sensitivity and spectral resolution of NMR spectroscopy are required for studies of larger nucleic acids or protein-nucleic acid complexes. Here we introduce a series of imino-proton-detected NMR experiments that yield an over 2-fold increase in sensitivity compared to conventional pulse schemes. These methods can be applied to the detection of base pair interactions, RNA-ligand titration experiments, measurement of residual dipolar (15)N-(1)H couplings, and direct measurements of conformational transitions. These NMR experiments employ longitudinal spin relaxation enhancement techniques that have proven useful in protein NMR spectroscopy. The performance of these new experiments is demonstrated for a 10 kDa TAR-TAR*(GA) RNA kissing complex and a 26 kDa tRNA.

A NMR strategy to unambiguously distinguish nucleic acid hairpin and duplex conformations applied to a Xist RNA A-repeat

All RNA sequences that fold into hairpins possess the intrinsic potential to form intermolecular duplexes because of their high self-complementarity. The thermodynamically more stable duplex conformation is favored under high salt conditions and at high RNA concentrations, posing a challenging problem for structural studies of small RNA hairpin conformations. We developed and applied a novel approach to unambiguously distinguish RNA hairpin and duplex conformations for the structural analysis of a Xist RNA A-repeat. Using a combination of a quantitative HNN-COSY experiment and an optimized double isotope-filtered NOESY experiment we could define the conformation of the 26-mer A-repeat RNA. In contrast to a previous secondary structure prediction of a double hairpin structure, the NMR data show that only the first predicted hairpin is formed, while the second predicted hairpin mediates dimerization of the A-repeat by duplex formation with a second A-repeat. The strategy employed here will be generally applicable to identify and quantify populations of hairpin and duplex conformations and to define RNA folding topology from inter- and intra-molecular base-pairing patterns.

Stability and kinetics of G-quadruplex structures

In this review, we give an overview of recent literature on the structure and stability of unimolecular G-rich quadruplex structures that are relevant to drug design and for in vivo function. The unifying theme in this review is energetics. The thermodynamic stability of quadruplexes has not been studied in the same detail as DNA and RNA duplexes, and there are important differences in the balance of forces between these classes of folded oligonucleotides. We provide an overview of the principles of stability and where available the experimental data that report on these principles. Significant gaps in the literature have been identified, that should be filled by a systematic study of well-defined quadruplexes not only to provide the basic understanding of stability both for design purposes, but also as it relates to in vivo occurrence of quadruplexes. Techniques that are commonly applied to the determination of the structure, stability and folding are discussed in terms of information content and limitations. Quadruplex structures fold and unfold comparatively slowly, and DNA unwinding events associated with transcription and replication may be operating far from equilibrium. The kinetics of formation and resolution of quadruplexes, and methodologies are discussed in the context of stability and their possible biological occurrence.

NMR-assisted prediction of RNA secondary structure: identification of a probable pseudoknot in the coding region of an R2 retrotransposon

As the rate of functional RNA sequence discovery escalates, high-throughput techniques for reliable structural determination are becoming crucial for revealing the essential features of these RNAs in a timely fashion. Computational predictions of RNA secondary structure quickly generate reasonable models but suffer from several approximations, including overly simplified models and incomplete knowledge of significant interactions. Similar problems limit the accuracy of predictions for other self-folding polymers, including DNA and peptide nucleic acid (PNA). The work presented here demonstrates that incorporating unassigned data from simple nuclear magnetic resonance (NMR) experiments into a dynamic folding algorithm greatly reduces the potential folding space of a given RNA and therefore increases the confidence and accuracy of modeling. This procedure has been packaged into an NMR-assisted prediction of secondary structure (NAPSS) algorithm that can produce pseudoknotted as well as non-pseudoknotted secondary structures. The method reveals a probable pseudoknot in the part of the coding region of the R2 retrotransposon from Bombyx mori that orchestrates second-strand DNA cleavage during insertion into the genome.

Remarkable metal counterion effect on the internucleotide J-couplings and chemical shifts of the N-H...N hydrogen bonds in the W-C base pairs

The effects of metal ion binding on the (2h) J(NN)-coupling and delta( (1)H)/Deltadelta( (15)N) chemical shifts of N-H...N H-bond units in internucleotide base pairs were explored by a combination of density functional theory calculations and molecular dynamics (MD) simulations. Results indicate that the NMR parameters vary considerably upon cation binding to the natural GC or AT base pairs, and thus can be used to identify the status of the base pairs, if cation-perturbed. The basic trend is that cation perturbation causes (2h) J(NN) to increase, Deltadelta( (15)N) to decrease, and delta( (1)H) to shift upfield for GC, and in the opposite directions for AT. The magnitudes of variation are closely related to the Lewis acidity of the metal ions. For both base pair series (M(z+)GC and M(z+)AT), these NMR parameters are linearly correlated among themselves. Their values depend strongly on the energy gaps (Delta(ELP-->sigma*)) and the second-order interaction energies ( E(2)) between the donor N lone pair (LP(N)) and the acceptor sigma* N-H localized NBO orbitals. In addition, the (2h) J NN changes are also sensitive to the amount of sigma charge transfer from LP(N) to sigma*(N-H) NBOs or from the purine to the pyrimidine moieties. The different trends are a consequence of the different H-bond patterns combined with the polarization effect of the metal ions in the cationized M(z+)AT series, M(z+) <-- A --> T, and the cationized GC series, M(z+) <-- G <-- C. The predicted cation-induced systematic trends of (2h) J(NN) and delta( (15)N, (1)H) in N-H...N H-bond units may provide a new approach to the determination of H-bond structure and strength in Watson-Crick base pairs, and provide an alternative probe of the heterogeneity of DNA sequences.

A review with comprehensive data on experimental indirect scalar NMR spin-spin coupling constants across hydrogen bonds

Scalar NMR spin-spin coupling constants across hydrogen bonds are fundamental in structural studies and as test grounds for theoretical calculations. Since they are scattered among many articles of different kinds, it seems useful to collect them in the most comprehensive way.

NMR methods for studying quadruplex nucleic acids

Solution NMR spectroscopy has traditionally played a central role in examining quadruplex structure, dynamics, and interactions. Here, an overview is given of the methods currently applied to structural, dynamics, thermodynamics, and kinetics studies of nucleic acid quadruplexes and associated cations.

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.

Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain

Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg2+-ions. Furthermore, a role for Mg2+-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA-adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg2+, but becomes partially pre-organized for ligand binding in the presence of Mg2+. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg2+ is more pronounced.

Characterizing the Cooperativity in H-Bonded Amino Structures

Density functional theory calculations were used to examine the effect of H-bond cooperativity on the magnitude of the NMR chemical shifts and spin-spin coupling constants in a C4h-symmetric G-quartet and in structures consisting of six cyanamide monomers. These included two ring structures (a planar C6h-symmetric structure and a nonplanar S6-symmetric structure) and two linear chain structures (a fully optimized planar Cs-symmetric chain and a planar chain structure where all intra- and intermolecular parameters were constrained to be identical). The NMR parameters were computed for the G-quartet and cyanamide structures, as well as for shorter fragments derived from these assemblies without reoptimization. In the ring structures and the chain with identical monomers, the intra- and intermolecular geometries of the cyanamides were identical, thereby allowing the study of cooperative effects in the absence of geometry changes. The magnitude of the |1JNH| coupling, 1H and 15N chemical shifts of the H-bonding amino N-H group, and the |h2JNN| H-bond coupling increased, whereas the size of the |1JNH| coupling of the non-H-bonded amino N-H bonds of the first amino group in the chain, which are roughly perpendicular to the H-bonding network, decreased in magnitude when H-bonding monomers were progressively added to extending ring or chain structures. These effects are attributed to electron redistribution induced by the presence of the nearby H-bonding guanine or cyanamide molecules.

Synthesis and characterization of a [3-15N]-labeled cis-syn thymine dimer-containing DNA duplex

Cis-syn thymine dimers are the major photoproducts of DNA and are the principal cause of mutations induced by sunlight. To better understand the nature of base pairing with cis-syn thymine dimers, we have synthesized a decamer oligodeoxynucleotide (ODN) containing a cis-syn thymine dimer labeled at the N3 of both T's with 15N by two efficient routes from [3-15N]-thymidine phosphoramidite. In the postsynthetic irradiation route, an ODN containing an adjacent pair of [3-15N]-labeled T's was irradiated and the cis-syn dimer-containing ODN isolated by HPLC. In the mixed building block route, a mixture of cis-syn and trans-syn dimer-containing ODNs was synthesized from a mixture of [3-15N]-labeled thymine dimer phosphoramidites after which the cis-syn dimer-containing ODN was isolated by HPLC. The N3-nitrogen and imino proton signals of an (15)N-labeled thymine dimer-containing decamer duplex were assigned by 2D 1H-15N heterocorrelated HSQC NMR spectroscopy, and the 15N-1H coupling constant was found to be 1.8 Hz greater for the 5'-T than for the 3'-T. The larger coupling constant is indicative of weaker H-bonding that is consistent with the more distorted nature of the 5'-base pair found in solution state NMR and crystallographic structures.

Solution structure of a purine rich hexaloop hairpin belonging to PGY/MDR1 mRNA and targeted by antisense oligonucleotides

A preferential target of antisense oligonucleotides directed against human PGY/MDR1 mRNA is a hairpin containing a stem with a G*U wobble pair, capped by the purine-rich 5'r(GGGAUG)3' hexaloop. This hairpin is studied by multidimensional NMR and restrained molecular dynamics, with special emphasis on the conformation of south sugars and non-standard phosphate linkages evidenced in both the stem and the loop. The hairpin is found to be highly structured. The G*U wobble pair, a strong counterion binding site, displays structural particularities that are characteristic of this type of mismatch. The upper part of the stem undergoes distortions that optimize its interactions with the beginning of the loop. The loop adopts a new fold in which the single-stranded GGGA purine tract is structured in A-like conformation stacked in continuity of the stem and displays an extensive hydrogen bonding surface for recognition. The remarkable hairpin stability results from classical inter- and intra-strand interactions reinforced by numerous hydrogen bonds involving unusual backbone conformations and ribose 2'-hydroxyl groups. Overall, this work emphasizes numerous features that account for the well-ordered structure of the whole hairpin and highlights the loop properties that facilitate interaction with antisense oligonucleotides.

RNA unrestrained molecular dynamics ensemble improves agreement with experimental NMR data compared to single static structure: a test case

Nuclear magnetic resonance (NMR) provides structural and dynamic information reflecting an average, often non-linear, of multiple solution-state conformations. Therefore, a single optimized structure derived from NMR refinement may be misleading if the NMR data actually result from averaging of distinct conformers. It is hypothesized that a conformational ensemble generated by a valid molecular dynamics (MD) simulation should be able to improve agreement with the NMR data set compared with the single optimized starting structure. Using a model system consisting of two sequence-related self-complementary ribonucleotide octamers for which NMR data was available, 0.3 ns particle mesh Ewald MD simulations were performed in the AMBER force field in the presence of explicit water and counterions. Agreement of the averaged properties of the molecular dynamics ensembles with NMR data such as homonuclear proton nuclear Overhauser effect (NOE)-based distance constraints, homonuclear proton and heteronuclear (1)H-(31)P coupling constant (J) data, and qualitative NMR information on hydrogen bond occupancy, was systematically assessed. Despite the short length of the simulation, the ensemble generated from it agreed with the NMR experimental constraints more completely than the single optimized NMR structure. This suggests that short unrestrained MD simulations may be of utility in interpreting NMR results. As expected, a 0.5 ns simulation utilizing a distance dependent dielectric did not improve agreement with the NMR data, consistent with its inferior exploration of conformational space as assessed by 2-D RMSD plots. Thus, ability to rapidly improve agreement with NMR constraints may be a sensitive diagnostic of the MD methods themselves.

Characterization of the monovalent ion position and hydrogen-bond network in guanine quartets by DFT calculations of NMR parameters

Conformational stability of G-quartets found in telomeric DNA quadruplex structures requires the coordination of monovalent ions. Here, an extensive Hartree-Fock and density functional theory analysis of the energetically favored position of Li+, Na+, and K+ ions is presented. The calculations show that at quartet-quartet distances observed in DNA quadruplex structures (3.3 A), the Li+ and Na+ ions favor positions of 0.55 and 0.95 A outside the plane of the G-quartet, respectively. The larger K+ ion prefers a central position between successive G-quartets. The energy barrier separating the minima in the quartet-ion-quartet model are much smaller for the Li+ and Na+ ions compared with the K+ ion; this suggests that K+ ions will not move as freely through the central channel of the DNA quadruplex. Spin-spin coupling constants and isotropic chemical shifts in G-quartets extracted from crystal structures of K+- and Na+-coordinated DNA quadruplexes were calculated with B3LYP/6-311G(d). The results show that the sizes of the trans-hydrogen-bond couplings are influenced primarily by the hydrogen bond geometry and only slightly by the presence of the ion. The calculations show that the R(N2N7) distance of the N2-H2...N7 hydrogen bond is characterized by strong correlations to both the chemical shifts of the donor group atoms and the (h2)J(N2N7) couplings. In contrast, weaker correlations between the (h3)J(N1C6') couplings and single geometric factors related to the N1-H1...O6=C6 hydrogen bond are observed. As such, deriving geometric information on the hydrogen bond through the use of trans-hydrogen-bond couplings and chemical shifts is more complex for the N1-H1...O6=C6 hydrogen bond than for the N2-H2...N7 moiety. The computed trans-hydrogen-bond couplings are shown to correlate with the experimentally determined couplings. However, the experimental values do not show such strong geometric dependencies.

NMR techniques for very large proteins and rnas in solution

Three-dimensional structure determination of small proteins and oligonucleotides by solution NMR is established. With the development of novel NMR and labeling techniques, structure determination is now feasible for proteins with a molecular mass of up to approximately 100 kDa and RNAs of up to 35 kDa. Beyond these molecular masses special techniques and approaches are required for applying NMR as a multiprobe method for structural investigations of proteins and RNAs. It is the aim of this review to summarize the NMR techniques and approaches available to advance the molecular mass limit of NMR both for proteins (up to 1 MDa) and RNAs (up to 100 kDa). Physical pictures of the novel techniques, their experimental applications, as well as labeling and assignment strategies are discussed and accompanied by future perspectives.

Dissecting non-canonical interactions in frameshift-stimulating mRNA pseudoknots

A variety of powerful NMR experiments have been introduced over the last few years that allow for the direct identification of different combinations of donor and acceptor atoms involved in hydrogen bonds in biomolecules. This ability to directly observe tertiary structural hydrogen bonds in solution tremendously facilitates structural studies of nucleic acids. We show here that an adiabatic HNN-COSY pulse scheme permits observation and measurement of J(N,N) couplings for nitrogen sites that are separated by up to 140 ppm in a single experiment at a proton resonance frequency of 500 MHz. Crucial hydrogen bond acceptor sites in nucleic acids, such as cytidine N3 nitrogens, can be unambiguously identified even in the absence of detectable H41 and H42 amino protons using a novel triple-resonance two-dimensional experiment, denoted H5(C5C4)N3. The unambiguous identification of amino nitrogen donor and aromatic nitrogen acceptor sites associated with both major groove as well as minor groove triple base pairs reveal the details of hydrogen bonding networks that stabilize the complex architecture of frameshift-stimulating mRNA pseudoknots. Another key tertiary interaction involving a 2'-OH hydroxyl proton that donates a hydrogen bond to an aromatic nitrogen acceptor in a cis Watson-Crick/sugar edge interaction can also be directly detected using a quantitative J(H,N) 1H,15N-HSQC experiment.

NMR methods for studying the structure and dynamics of RNA

Proper functioning of RNAs requires the formation of complex three-dimensional structures combined with the ability to rapidly interconvert between multiple functional states. This review covers recent advances in isotope-labeling strategies and NMR experimental approaches that have promise for facilitating solution structure determinations and dynamics studies of biologically active RNAs. Improved methods for the production of isotopically labeled RNAs combined with new multidimensional heteronuclear NMR experiments make it possible to dramatically reduce spectral crowding and simplify resonance assignments for RNAs. Several novel applications of experiments that directly detect hydrogen-bonding interactions are discussed. These studies demonstrate how NMR spectroscopy can be used to distinguish between possible secondary structures and identify mechanisms of ligand binding in RNAs. A variety of recently developed methods for measuring base and sugar residual dipolar couplings are described. NMR residual dipolar coupling techniques provide valuable data for determining the long-range structure and orientation of helical regions in RNAs. A number of studies are also presented where residual dipolar coupling constraints are used to determine the global structure and dynamics of RNAs. NMR relaxation data can be used to probe the dynamics of macromolecules in solution. The power dependence of transverse rotating-frame relaxation rates was used here to study dynamics in the minimal hammerhead ribozyme. Improved methods for isotopically labeling RNAs combined with new types of structural data obtained from a growing repertoire of NMR experiments are facilitating structural and dynamic studies of larger RNAs.

NMR structure and Mg2+ binding of an RNA segment that underlies the L7/L12 stalk in the E.coli 50S ribosomal subunit

Helix 42 of Domain II of Escherichia coli 23S ribosomal RNA underlies the L7/L12 stalk in the ribosome and may be significant in positioning this feature relative to the rest of the 50S ribosomal subunit. Unlike the Haloarcula marismortui and Deinococcus radiodurans examples, the lower portion of helix 42 in E.coli contains two consecutive G*A oppositions with both adenines on the same side of the stem. Herein, the structure of an analog of positions 1037-1043 and 1112-1118 in the helix 42 region is reported. NMR spectra and structure calculations support a cis Watson-Crick/Watson-Crick (cis W.C.) G*A conformation for the tandem (G*A)2 in the analog and a minimally perturbed helical duplex stem. Mg2+ titration studies imply that the cis W.C. geometry of the tandem (G*A)2 probably allows O6 of G20 and N1 of A4 to coordinate with a Mg2+ ion as indicated by the largest chemical shift changes associated with the imino group of G20 and the H8 of G20 and A4. A cross-strand bridging Mg2+ coordination has also been found in a different sequence context in the crystal structure of H.marismortui 23S rRNA, and therefore it may be a rare but general motif in Mg2+ coordination.

NMR spectroscopy of RNA

NMR spectroscopy is a powerful tool for studying proteins and nucleic acids in solution. This is illustrated by the fact that nearly half of all current RNA structures were determined by using NMR techniques. Information about the structure, dynamics, and interactions with other RNA molecules, proteins, ions, and small ligands can be obtained for RNA molecules up to 100 nucleotides. This review provides insight into the resonance assignment methods that are the first and crucial step of all NMR studies, into the determination of base-pair geometry, into the examination of local and global RNA conformation, and into the detection of interaction sites of RNA. Examples of NMR investigations of RNA are given by using several different RNA molecules to illustrate the information content obtainable by NMR spectroscopy and the applicability of NMR techniques to a wide range of biologically interesting RNA molecules.

Cofactor-Apoprotein Hydrogen Bonding in Oxidized and Fully Reduced Flavodoxin Monitored by Trans-Hydrogen-Bond Scalar Couplings

Hydrogen bonding plays a key role in the tight binding of the FMN cofactor and the regulation of its redox properties in flavodoxins. Hydrogen bonding interactions can be directly observed in solution by multidimensional heteronuclear NMR spectroscopy through the scalar couplings between donor and acceptor nuclei. Here we report on the detection of intermolecular trans-hydrogen-bond couplings ((h)J) between the flavin ring system and the backbone of Desulfovibrio vulgaris flavodoxin in the oxidized and the two-electron reduced states. For this purpose, experiments are adapted from pulse sequences previously applied to determining (h)J coupling constants in nucleic acid-base pairs and proteins. The resulting (h2)J(N,N), (h4)J(N,N), (h3)J(C,N), and (h1)J(H,N) couplings involve the (15)N(1), (13)C(2), and (15)N(3) nuclei of the pyrimidine moiety of FMN, whereas no such interactions are detectable for (13)C(4) and (15)N(5). Several long-range (15)N-(15)N, (13)C-(15)N, and (1)H-(15)N J-coupling constants within the flavin are obtained as "by-products". The magnitudes of both (h)J and regular J couplings are found to be dependent on the redox state. In general, good correlations between (h)J coupling constants and donor-group (1)H chemical shifts and also crystallographic donor-acceptor distances are observed.

Method for direct discrimination of intra- and intermolecular hydrogen bonds, and characterization of the G(:A):G(:A):G(:A):G heptad, with scalar couplings across hydrogen bonds

Discrimination of intra- and intermolecular hydrogen bonds in a symmetric multimer has not been accomplished yet, although such discrimination would provide a crucial basis for construction of the multimeric architecture of nucleic acids by NMR. We have developed a direct and unambiguous method for such discrimination involving the use of scalar couplings across hydrogen bonds. The method has been validated with a symmetric dimer of d(GGGCTTTTGGGC), for which the structure including both intra- and intermolecular hydrogen bonds was already reported. This has demonstrated that our method can clearly discriminate these two kinds of hydrogen bonds. Then, the method was applied to a symmetric dimer of d(GGAGGAGGAGGA) and has provided decisive information on its multimeric architecture. Additionally, the values for scalar couplings across hydrogen bonds for G:G and G:A base pairs in the G(:A):G(:A):G(:A):G heptad formed by d(GGAGGAGGAGGA) were determined for the first time. This determination has provided an insight into the nature of the heptad.

Identification of NH...N hydrogen bonds by magic angle spinning solid state NMR in a double-stranded RNA associated with myotonic dystrophy

RNA plays a central role in biological processes and exhibits a variety of secondary and tertiary structural features that are often stabilized via hydrogen bonds. The distance between the donor and acceptor nitrogen nuclei involved in NH...N hydrogen bonds in nucleic acid base pairs is typically in the range of 2.6-2.9 A. Here, we show for the first time that such spatial proximity between 15N nitrogen nuclei can be conveniently monitored via magic angle spinning solid state NMR on a uniformly 15N-labelled RNA. The presence of NH.N hydrogen bonds is reflected as cross-peaks between the donor and acceptor nitrogen nuclei in 2D 15N dipolar chemical shift correlation spectra. The RNA selected for this experimental study was a CUG repeat expansion implicated in the neuromuscular disease myotonic dystrophy. The results presented provide direct evidence that the CUG repeat expansion adopts a double-stranded conformation.

The structure of the stemloop D subdomain of coxsackievirus B3 cloverleaf RNA and its interaction with the proteinase 3C

Stemloop D (SLD) of the 5' cloverleaf RNA is the cognate ligand of the coxsackievirus B3 (CVB3) 3C proteinase (3Cpro). Both are indispensable components of the viral replication initiation complex. SLD is a structurally autonomous subunit of the 5' cloverleaf. The SLD structure was solved by NMR spectroscopy to an rms deviation of 0.66 A (all heavy atoms). SLD contains a novel triple pyrimidine mismatch motif with a central Watson-Crick type C:U pair. SLD is capped by an apical uCACGg tetraloop adopting a structure highly similar to stable cUNCGg tetraloops. Binding of CVB3 3Cpro induces changes in NMR spectra for nucleotides adjacent to the triple pyrimidine mismatch and of the tetraloop implying them as sites of specific SLD:3Cpro interaction. The binding of 3Cpro to SLD requires the integrity of those structural elements, strongly suggesting that 3Cpro recognizes a structural motif instead of a specific sequence.

Structural basis of the highly efficient trapping of the HIV Tat protein by an RNA aptamer

An RNA aptamer containing two binding sites exhibits extremely high affinity to the HIV Tat protein. We have determined the structure of the aptamer complexed with two argininamide molecules. Two adjacent U:A:U base triples were formed, which widens the major groove to make space for the two argininamide molecules. The argininamide molecules bind to the G bases through hydrogen bonds. The binding is stabilized through stacking interactions. The structure of the aptamer complexed with a Tat-derived arginine-rich peptide was also characterized. It was suggested that the aptamer structure is similar for both complexes and that the aptamer interacts with two different arginine residues of the peptide simultaneously at the two binding sites, which could explain the high affinity to Tat.