NMR structure of the 101-nucleotide core encapsidation signal of the Moloney murine leukemia virus

Mechanisms of human immunodeficiency virus type 2 RNA packaging: efficient trans packaging and selection of RNA copackaging partners

Human immunodeficiency virus type 2 (HIV-2) has been reported to have a distinct RNA packaging mechanism, referred to as cis packaging, in which Gag proteins package the RNA from which they were translated. We examined the progeny generated from dually infected cell lines that contain two HIV-2 proviruses, one with a wild-type gag/gag-pol and the other with a mutant gag that cannot express functional Gag/Gag-Pol. Viral titers and RNA analyses revealed that mutant viral RNAs can be packaged at efficiencies comparable to that of viral RNA from which wild-type Gag/Gag-Pol is translated. These results do not support the cis-packaging hypothesis but instead indicate that trans packaging is the major mechanism of HIV-2 RNA packaging. To further characterize the mechanisms of HIV-2 RNA packaging, we visualized HIV-2 RNA in individual particles by using fluorescent protein-tagged RNA-binding proteins that specifically recognize stem-loop motifs in the viral genomes, an assay termed single virion analysis. These studies revealed that >90% of the HIV-2 particles contained viral RNAs and that RNAs derived from different viruses were copackaged frequently. Furthermore, the frequencies of heterozygous particles in the viral population could be altered by changing a 6-nucleotide palindromic sequence at the 5'-untranslated region of the HIV-2 genome. This finding indicates that selection of copackaging RNA partners occurs prior to encapsidation and that HIV-2 Gag proteins primarily package one dimeric RNA rather than two monomeric RNAs. Additionally, single virion analyses demonstrated a similar RNA distribution in viral particles regardless of whether both viruses had a functional gag or one of the viruses had a nonfunctional gag, providing further support for the trans-packaging hypothesis. Together, these results revealed mechanisms of HIV-2 RNA packaging that are, contrary to previous studies, in many respects surprisingly similar to those of HIV-1.

Physics-based de novo prediction of RNA 3D structures

Current experiments on structural determination cannot keep up the pace with the steadily emerging RNA sequences and new functions. This underscores the request for an accurate model for RNA three-dimensional (3D) structural prediction. Although considerable progress has been made in mechanistic studies, accurate prediction for RNA tertiary folding from sequence remains an unsolved problem. The first and most important requirement for the prediction of RNA structure from physical principles is an accurate free energy model. A recently developed three-vector virtual bond-based RNA folding model ("Vfold") has allowed us to compute the chain entropy and predict folding free energies and structures for RNA secondary structures and simple pseudoknots. Here we develop a free energy-based method to predict larger more complex RNA tertiary folds. The approach is based on a multiscaling strategy: from the nucleotide sequence, we predict the two-dimensional (2D) structures (defined by the base pairs and tertiary contacts); based on the 2D structure, we construct a 3D scaffold; with the 3D scaffold as the initial state, we combine AMBER energy minimization and PDB-based fragment search to predict the all-atom structure. A key advantage of the approach is the statistical mechanical calculation for the conformational entropy of RNA structures, including those with cross-linked loops. Benchmark tests show that the model leads to significant improvements in RNA 3D structure prediction.

Structure of a conserved retroviral RNA packaging element by NMR spectroscopy and cryo-electron tomography

The 5'-untranslated regions of all gammaretroviruses contain a conserved "double-hairpin motif" (Ψ(CD)) that is required for genome packaging. Both hairpins (SL-C and SL-D) contain GACG tetraloops that, in isolated RNAs, are capable of forming "kissing" interactions stabilized by two intermolecular G-C base pairs. We have determined the three-dimensional structure of the double hairpin from the Moloney murine leukemia virus ([Ψ(CD)](2), 132 nt, 42.8 kDa) using a (2)H-edited NMR-spectroscopy-based approach. This approach enabled the detection of (1)H-(1)H dipolar interactions that were not observed in previous studies of isolated SL-C and SL-D hairpin RNAs using traditional (1)H-(1)H correlated and (1)H-(13)C-edited NMR methods. The hairpins participate in intermolecular cross-kissing interactions (SL-C to SL-D' and SLC' to SL-D) and stack in an end-to-end manner (SL-C to SL-D and SL-C' to SL-D') that gives rise to an elongated overall shape (ca 95 Å×45 Å×25 Å). The global structure was confirmed by cryo-electron tomography (cryo-ET), making [Ψ(CD)](2) simultaneously the smallest RNA to be structurally characterized to date by cryo-ET and among the largest to be determined by NMR. Our findings suggest that, in addition to promoting dimerization, [Ψ(CD)](2) functions as a scaffold that helps initiate virus assembly by exposing a cluster of conserved UCUG elements for binding to the cognate nucleocapsid domains of assembling viral Gag proteins.

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.

Comparisons between chemical mapping and binding to isoenergetic oligonucleotide microarrays reveal unexpected patterns of binding to the Bacillus subtilis RNase P RNA specificity domain

Microarrays with isoenergetic pentamer and hexamer 2'-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on the RNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1 M Na(+) and 10 mM Mg(2+).

Preformed protein-binding motifs in 7SK snRNA: structural and thermodynamic comparisons with retroviral TAR

The 7SK small nuclear RNA is a highly conserved non-coding RNA that regulates transcriptional elongation. 7SK utilizes the HEXIM proteins to sequester the transcription factor P-TEFb by a mechanism similar to that used by retroviral TAR RNA to engage Tat and P-TEFb. Tat has also recently been shown to bind 7SK directly and recruit P-TEFb to TAR. We report here the solution structures of the free and arginine-bound forms of stem loop 4 of 7SK (7SK-SL4). Comparison of the 7SK-SL4 and TAR structures demonstrates the presence of a common arginine sandwich motif. However, arginine binding to 7SK-SL4 is mechanistically distinct and occurs via docking into a pre-organized pocket resulting in a 1000-fold increased affinity. Furthermore, whereas formation of the binding pocket in TAR requires a critical base-triple, hydrogen-bond formation between the equivalent bases in 7SK-SL4 is not essential and the pocket is stabilized solely by a pseudo base-triple platform. In addition, this theme of preformed protein binding motifs also extends into the pentaloop. The configuration of the loop suggests that 7SK-SL4 is poised to make ternary contacts with P-TEFb and HEXIM or Tat. These key differences between 7SK-SL4 and TAR present an opportunity to understand RNA structural adaptation and have implications for understanding differential interactions with Tat.

Major groove width variations in RNA structures determined by NMR and impact of 13C residual chemical shift anisotropy and 1H-13C residual dipolar coupling on refinement

Ribonucleic acid structure determination by NMR spectroscopy relies primarily on local structural restraints provided by (1)H- (1)H NOEs and J-couplings. When employed loosely, these restraints are broadly compatible with A- and B-like helical geometries and give rise to calculated structures that are highly sensitive to the force fields employed during refinement. A survey of recently reported NMR structures reveals significant variations in helical parameters, particularly the major groove width. Although helical parameters observed in high-resolution X-ray crystal structures of isolated A-form RNA helices are sensitive to crystal packing effects, variations among the published X-ray structures are significantly smaller than those observed in NMR structures. Here we show that restraints derived from aromatic (1)H- (13)C residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) can overcome NMR restraint and force field deficiencies and afford structures with helical properties similar to those observed in high-resolution X-ray structures.

Site-specific labeling of nucleotides for making RNA for high resolution NMR studies using an E. coli strain disabled in the oxidative pentose phosphate pathway

Escherichia coli (E. coli) is a versatile organism for making nucleotides labeled with stable isotopes ((13)C, (15)N, and/or (2)H) for structural and molecular dynamics characterizations. Growth of a mutant E. coli strain deficient in the pentose phosphate pathway enzyme glucose-6-phosphate dehydrogenase (K10-1516) on 2-(13)C-glycerol and (15)N-ammonium sulfate in Studier minimal medium enables labeling at sites useful for NMR spectroscopy. However, (13)C-sodium formate combined with (13)C-2-glycerol in the growth media adds labels to new positions. In the absence of labeled formate, both C5 and C6 positions of the pyrimidine rings are labeled with minimal multiplet splitting due to (1)J(C5C6) scalar coupling. However, the C2/C8 sites within purine rings and the C1'/C3'/C5' positions within the ribose rings have reduced labeling. Addition of (13)C-labeled formate leads to increased labeling at the base C2/C8 and the ribose C1'/C3'/C5' positions; these new specific labels result in two- to three-fold increase in the number of resolved resonances. This use of formate and (15)N-ammonium sulfate promises to extend further the utility of these alternate site specific labels to make labeled RNA for downstream biophysical applications such as structural, dynamics and functional studies of interesting biologically relevant RNAs.

Isotope labeling strategies for NMR studies of RNA

The known biological functions of RNA have expanded in recent years and now include gene regulation, maintenance of sub-cellular structure, and catalysis, in addition to propagation of genetic information. As for proteins, RNA function is tightly correlated with structure. Unlike proteins, structural information for larger, biologically functional RNAs is relatively limited. NMR signal degeneracy, relaxation problems, and a paucity of long-range (1)H-(1)H dipolar contacts have limited the utility of traditional NMR approaches. Selective isotope labeling, including nucleotide-specific and segmental labeling strategies, may provide the best opportunities for obtaining structural information by NMR. Here we review methods that have been developed for preparing and purifying isotopically labeled RNAs, as well as NMR strategies that have been employed for signal assignment and structure determination.

Prediction of geometrically feasible three-dimensional structures of pseudoknotted RNA through free energy estimation

Accurate free energy estimation is essential for RNA structure prediction. The widely used Turner's energy model works well for nested structures. For pseudoknotted RNAs, however, there is no effective rule for estimation of loop entropy and free energy. In this work we present a new free energy estimation method, termed the pseudoknot predictor in three-dimensional space (pk3D), which goes beyond Turner's model. Our approach treats nested and pseudoknotted structures alike in one unifying physical framework, regardless of how complex the RNA structures are. We first test the ability of pk3D in selecting native structures from a large number of decoys for a set of 43 pseudoknotted RNA molecules, with lengths ranging from 23 to 113. We find that pk3D performs slightly better than the Dirks and Pierce extension of Turner's rule. We then test pk3D for blind secondary structure prediction, and find that pk3D gives the best sensitivity and comparable positive predictive value (related to specificity) in predicting pseudoknotted RNA secondary structures, when compared with other methods. A unique strength of pk3D is that it also generates spatial arrangement of structural elements of the RNA molecule. Comparison of three-dimensional structures predicted by pk3D with the native structure measured by nuclear magnetic resonance or X-ray experiments shows that the predicted spatial arrangement of stems and loops is often similar to that found in the native structure. These close-to-native structures can be used as starting points for further refinement to derive accurate three-dimensional structures of RNA molecules, including those with pseudoknots.

An RNA structural switch regulates diploid genome packaging by Moloney murine leukemia virus

Retroviruses selectively package two copies of their RNA genomes via mechanisms that have yet to be fully deciphered. Recent studies with small fragments of the Moloney murine leukemia virus (MoMuLV) genome suggested that selection may be mediated by an RNA switch mechanism, in which conserved UCUG elements that are sequestered by base-pairing in the monomeric RNA become exposed upon dimerization to allow binding to the cognate nucleocapsid (NC) domains of the viral Gag proteins. Here we show that a large fragment of the MoMuLV 5' untranslated region that contains all residues necessary for efficient RNA packaging (Psi(WT); residues 147-623) also exhibits a dimerization-dependent affinity for NC, with the native dimer ([Psi(WT)](2)) binding 12+/-2 NC molecules with high affinity (K(d)=17+/-7 nM) and with the monomer, stabilized by substitution of dimer-promoting loop residues with hairpin-stabilizing sequences (Psi(M)), binding 1-2 NC molecules. Identical dimer-inhibiting mutations in MoMuLV-based vectors significantly inhibit genome packaging in vivo (approximately 100-fold decrease), whereas a large deletion of nearly 200 nucleotides just upstream of the gag start codon has minimal effects. Our findings support the proposed RNA switch mechanism and further suggest that virus assembly may be initiated by a complex comprising as few as 12 Gag molecules bound to a dimeric packaging signal.

Secondary structure of the mature ex virio Moloney murine leukemia virus genomic RNA dimerization domain

Retroviral genomes are dimeric, comprised of two sense-strand RNAs linked at their 5' ends by noncovalent base pairing and tertiary interactions. Viral maturation involves large-scale morphological changes in viral proteins and in genomic RNA dimer structures to yield infectious virions. Structural studies have largely focused on simplified in vitro models of genomic RNA dimers even though the relationship between these models and authentic viral RNA is unknown. We evaluate the secondary structure of the minimal dimerization domain in genomes isolated from Moloney murine leukemia virions using a quantitative and single nucleotide resolution RNA structure analysis technology (selective 2'-hydroxyl acylation analyzed by primer extension, or SHAPE). Results are consistent with an architecture in which the RNA dimer is stabilized by four primary interactions involving two sets of intermolecular base pairs and two loop-loop interactions. The dimerization domain can independently direct its own folding since heating and refolding reproduce the same structure as visualized in genomic RNA isolated from virions. Authentic ex virio RNA has a SHAPE reactivity profile similar to that of a simplified transcript dimer generated in vitro, with the important exception of a region that appears to form a compact stem-loop only in the virion-isolated RNA. Finally, we analyze the conformational changes that accompany folding of monomers into dimers in vitro. These experiments support well-defined structural models for an authentic dimerization domain and also emphasize that many features of mature genomic RNA dimers can be reproduced in vitro using properly designed, simplified RNAs.

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.

Measurement of imino 1H-1H residual dipolar couplings in RNA

Imino (1)H-(15)N residual dipolar couplings (RDCs) provide additional structural information that complements standard (1)H-(1)H NOEs leading to improvements in both the local and global structure of RNAs. Here, we report measurement of imino (1)H-(1)H RDCs for the Iron Responsive Element (IRE) RNA and native E. coli tRNA(Val) using a BEST-Jcomp-HMQC2 experiment. (1)H-(1)H RDCs are observed between the imino protons in G-U wobble base pairs and between imino protons on neighboring base pairs in both RNAs. These imino (1)H-(1)H RDCs complement standard (1)H-(15)N RDCs because the (1)H-(1)H vectors generally point along the helical axis, roughly perpendicular to (1)H-(15)N RDCs. The use of longitudinal relaxation enhancement increased the signal-to-noise of the spectra by ~3.5-fold over the standard experiment. The ability to measure imino (1)H-(1)H RDCs offers a new restraint, which can be used in NMR domain orientation and structural studies of RNAs.

Solution structure of tRNAVal from refinement of homology model against residual dipolar coupling and SAXS data

A procedure is presented for refinement of a homology model of E. coli tRNA(Val), originally based on the X-ray structure of yeast tRNA(Phe), using experimental residual dipolar coupling (RDC) and small angle X-ray scattering (SAXS) data. A spherical sampling algorithm is described for refinement against SAXS data that does not require a globbic approximation, which is particularly important for nucleic acids where such approximations are less appropriate. Substantially higher speed of the algorithm also makes its application favorable for proteins. In addition to the SAXS data, the structure refinement employed a sparse set of NMR data consisting of 24 imino N-H(N) RDCs measured with Pf1 phage alignment, and 20 imino N-H(N) RDCs obtained from magnetic field dependent alignment of tRNA(Val). The refinement strategy aims to largely retain the local geometry of the 58% identical tRNA(Phe) by ensuring that the atomic coordinates for short, overlapping segments of the ribose-phosphate backbone and the conserved base pairs remain close to those of the starting model. Local coordinate restraints are enforced using the non-crystallographic symmetry (NCS) term in the XPLOR-NIH or CNS software package, while still permitting modest movements of adjacent segments. The RDCs mainly drive the relative orientation of the helical arms, whereas the SAXS restraints ensure an overall molecular shape compatible with experimental scattering data. The resulting structure exhibits good cross-validation statistics (jack-knifed Q (free) = 14% for the Pf1 RDCs, compared to 25% for the starting model) and exhibits a larger angle between the two helical arms than observed in the X-ray structure of tRNA(Phe), in agreement with previous NMR-based tRNA(Val) models.

Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA

The functions of most RNA molecules are critically dependent on the distinct local dynamics that characterize secondary structure and tertiary interactions and on structural changes that occur upon binding by proteins and small molecule ligands. Measurements of RNA dynamics at nucleotide resolution set the foundation for understanding the roles of individual residues in folding, catalysis, and ligand recognition. In favorable cases, local order in small RNAs can be quantitatively analyzed by NMR in terms of a generalized order parameter, S2. Alternatively, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry measures local nucleotide flexibility in RNAs of any size using structure-sensitive reagents that acylate the 2'-hydroxyl position. In this work, we compare per-residue RNA dynamics, analyzed by both S2 and SHAPE, for three RNAs: the HIV-1 TAR element, the U1A protein binding site, and the Tetrahymena telomerase stem loop 4. We find a very strong correlation between the two measurements: nucleotides with high SHAPE reactivities consistently have low S2 values. We conclude that SHAPE chemistry quantitatively reports local nucleotide dynamics and can be used with confidence to analyze dynamics in large RNAs, RNA-protein complexes, and RNAs in vivo.

Retroviral nucleocapsid proteins display nonequivalent levels of nucleic acid chaperone activity

Human immunodeficiency virus type 1 (HIV-1) nucleocapsid protein (NC) is a nucleic acid chaperone that facilitates the remodeling of nucleic acids during various steps of the viral life cycle. Two main features of NC's chaperone activity are its abilities to aggregate and to destabilize nucleic acids. These functions are associated with NC's highly basic character and with its zinc finger domains, respectively. While the chaperone activity of HIV-1 NC has been extensively studied, less is known about the chaperone activities of other retroviral NCs. In this work, complementary experimental approaches were used to characterize and compare the chaperone activities of NC proteins from four different retroviruses: HIV-1, Moloney murine leukemia virus (MLV), Rous sarcoma virus (RSV), and human T-cell lymphotropic virus type 1 (HTLV-1). The different NCs exhibited significant differences in their overall chaperone activities, as demonstrated by gel shift annealing assays, decreasing in the order HIV-1 approximately RSV > MLV >> HTLV-1. In addition, whereas HIV-1, RSV, and MLV NCs are effective aggregating agents, HTLV-1 NC, which exhibits poor overall chaperone activity, is unable to aggregate nucleic acids. Measurements of equilibrium binding to single- and double-stranded oligonucleotides suggested that all four NC proteins have moderate duplex destabilization capabilities. Single-molecule DNA-stretching studies revealed striking differences in the kinetics of nucleic acid dissociation between the NC proteins, showing excellent correlation between nucleic acid dissociation kinetics and overall chaperone activity.

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.

Key labeling technologies to tackle sizeable problems in RNA structural biology

The ability to adopt complex three-dimensional (3D) structures that can rapidly interconvert between multiple functional states (folding and dynamics) is vital for the proper functioning of RNAs. Consequently, RNA structure and dynamics necessarily determine their biological function. In the post-genomic era, it is clear that RNAs comprise a larger proportion (>50%) of the transcribed genome compared to proteins (< or =2%). Yet the determination of the 3D structures of RNAs lags considerably behind those of proteins and to date there are even fewer investigations of dynamics in RNAs compared to proteins. Site specific incorporation of various structural and dynamic probes into nucleic acids would likely transform RNA structural biology. Therefore, various methods for introducing probes for structural, functional, and biotechnological applications are critically assessed here. These probes include stable isotopes such as (2)H, (13)C, (15)N, and (19)F. Incorporation of these probes using improved RNA ligation strategies promises to change the landscape of structural biology of supramacromolecules probed by biophysical tools such as nuclear magnetic resonance (NMR) spectroscopy, X-ray crystallography and Raman spectroscopy. Finally, some of the structural and dynamic problems that can be addressed using these technological advances are outlined.

Multiple segmental and selective isotope labeling of large RNA for NMR structural studies

Multiple segmental and selective isotope labeling of RNA with three segments has been demonstrated by introducing an RNA segment, selectively labeled with ¹³C₉/¹⁵N₂/²H_{(1′, 3′, 4′, 5′, 5′′)}-labeled uridine residues, into the central position of the 20 kDa ε-RNA of Duck Hepatitis B Virus. The RNA molecules were produced via two efficient protocols: a two-step protocol, which uses T4 DNA ligase and T4 RNA ligase 1, and a one-pot protocol, which uses T4 RNA ligase 1 alone. With T4 RNA ligase 1 all not-to-be-ligated termini are usually protected to prevent formation of side products. We show that such labor-intensive protection of termini is not required, provided segmentation sites can be chosen such that the segments fold into the target structure or target-like structures and thus are not trapped into stable alternate structures. These sites can be reliably predicted via DINAMelt. The simplified NMR spectrum provided evidence for the presence of a U28 H₃-imino resonance, previously obscured in the fully labeled sample, and thus of the non-canonical base pair U28:C37. The demonstrated multiple segmental labeling protocols are generally applicable to large RNA molecules and can be extended to more than three segments.

6th international symposium on retroviral nucleocapsid

Retroviruses and LTR-retrotransposons are widespread in all living organisms and, in some instances such as for HIV, can be a serious threat to the human health. The retroviral nucleocapsid is the inner structure of the virus where several hundred nucleocapsid protein (NC) molecules coat the dimeric, genomic RNA. During the past twenty years, NC was found to play multiple roles in the viral life cycle (Fig. 1), notably during the copying of the genomic RNA into the proviral DNA by viral reverse transcriptase and integrase, and is therefore considered to be a prime target for anti-HIV therapy. The 6th NC symposium was held in the beautiful city of Amsterdam, the Netherlands, on the 20th and 21st of September 2007. All aspects of NC biology, from structure to function and to anti-HIV vaccination, were covered during this meeting.

Analyzing the flexibility of RNA structures by constraint counting

RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA(ASP) structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length >/=40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.

Preparation of large RNA oligonucleotides with complementary isotope-labeled segments for NMR structural studies

RNA structure determination by solution NMR spectroscopy is often restricted to small RNAs (<15 kDa) owing to the problem of chemical shift degeneracy. A fruitful coupling of novel NMR techniques with segmental RNA labeling methodologies could be a powerful tool to overcome the molecular mass limitation of RNA NMR spectroscopy. Herein, we describe a time- and cost-effective procedure to prepare and purify segmentally labeled large RNAs. Two sets of RNA fragments with complementary labeling schemes, such as one fragment (13)C- and the other (15)N-labeled, are prepared by in vitro transcription from a single plasmid DNA. The desired RNA fragments are excised from the primary transcript by two cis-acting hammerhead ribozymes, yielding the required engineered ends for subsequent, complementary ligation. The resulting RNA oligonucleotides display NMR spectra with greatly reduced resonance overlap and thus enable NMR studies of smaller labeled RNA segments within the native context of a large RNA. The procedure is expected to take 3-4 weeks to implement.

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.

Solution NMR of large molecules and assemblies

Solution NMR spectroscopy represents a powerful tool for examining the structure and function of biological macromolecules. The advent of multidimensional (2D-4D) NMR, together with the widespread use of uniform isotopic labeling of proteins and RNA with the NMR-active isotopes, 15N and 13C, opened the door to detailed analyses of macromolecular structure, dynamics, and interactions of smaller macromolecules (< approximately 25 kDa). Over the past 10 years, advances in NMR and isotope labeling methods have expanded the range of NMR-tractable targets by at least an order of magnitude. Here we briefly describe the methodological advances that allow NMR spectroscopy of large macromolecules and their complexes and provide a perspective on the wide range of applications of NMR to biochemical problems.

Solution structure of the Rous sarcoma virus nucleocapsid protein: muPsi RNA packaging signal complex

The 5'-untranslated region (5'-UTR) of retroviral genomes contains elements required for genome packaging during virus assembly. For many retroviruses, the packaging elements reside in non-contiguous segments that span most or all of the 5'-UTR. The Rous sarcoma virus (RSV) is an exception, in that its genome can be packaged efficiently by a relatively short, 82 nt segment of the 5'-UTR called muPsi. The RSV 5'-UTR also contains three translational start codons (AUG-1, AUG-2 and AUG-3) that have been controvertibly implicated in translation initiation and genome packaging, one of which (AUG-3) resides within the muPsi sequence. We demonstrated recently that muPsi is capable of binding to the cognate RSV nucleocapsid protein (NC) with high affinity (dissociation constant K(d) approximately 2 nM), and that residues of AUG-3 are essential for tight binding. We now report the solution structure of the NC:muPsi complex, determined using NMR data obtained for samples containing ((13)C,(15)N)-labeled NC and (2)H-enriched, nucleotide-specifically protonated RNAs. Upon NC binding, muPsi adopts a stable secondary structure that consists of three stem loops (SL-A, SL-B and SL-C) and an 8 bp stem (O3). Binding is mediated by the two zinc knuckle domains of NC. The N-terminal knuckle interacts with a conserved U(217)GCG tetraloop (a member of the UNCG family; N=A,U,G or C), and the C-terminal zinc knuckle binds to residues that flank SL-A, including residues of AUG-3. Mutations of critical nucleotides in these sequences compromise or abolish viral infectivity. Our studies reveal novel structural features important for NC:RNA binding, and support the hypothesis that AUG-3 is conserved for genome packaging rather than translational control.

The SL1-SL2 (stem-loop) domain is the primary determinant for stability of the gamma retroviral genomic RNA dimer

Retroviral genomes are assembled from two sense-strand RNAs by noncovalent interactions at their 5' ends, forming a dimer. The RNA dimerization domain is a potential target for antiretroviral therapy and represents a compelling RNA folding problem. The fundamental dimerization unit for the Moloney murine sarcoma gamma retrovirus spans a 170-nucleotide minimal dimerization active sequence. In the dimer, two self-complementary sequences, PAL1 and PAL2, form intermolecular duplexes, and an SL1-SL2 (stem-loop) domain forms loop-loop base pairs, mediated by GACG tetraloops, and extensive tertiary interactions. To develop a framework for assembly of the retroviral RNA dimer, we quantified the stability of and established nucleotide resolution secondary structure models for sequence variants in which each motif was compromised. Base pairing and tertiary interactions between SL1-SL2 domains contribute a large free energy increment of -10 kcal/mol. In contrast, even though the PAL1 and PAL2 intermolecular duplexes span 10 and 16 bp in the dimer, respectively, they contribute only -2.5 kcal/mol to stability, roughly equal to a single new base pair. First, these results emphasize that the energetic costs for disrupting interactions in the monomer state nearly balance the PAL1 and PAL2 base pairing interactions that form in the dimer. Second, intermolecular duplex formation plays a biological role distinct from simply stabilizing the structure of the retroviral genomic RNA dimer.

Interpreting oligonucleotide microarray data to determine RNA secondary structure: application to the 3' end of Bombyx mori R2 RNA

A method to deduce RNA secondary structure on the basis of data from microarrays of 2'-O-methyl RNA 9-mers immobilized in agarose film on glass slides is tested with a 249 nucleotide RNA from the 3' end of the R2 retrotransposon from Bombyx mori. Various algorithms incorporating binding data and free-energy minimization calculations were compared for interpreting the data to provide possible secondary structures. Two different methods give structures with 100 and 87% of the base pairs determined by sequence comparison. In contrast, structures predicted by free-energy minimization alone by Mfold and RNAstructure contain 52 and 72% of the known base pairs, respectively. This combination of high throughput microarray techniques with algorithms using free-energy calculations has potential to allow for fast determination of RNA secondary structure. It should also facilitate the design of antisense and siRNA oligonucleotides.

Fluorescence and solution NMR study of the active site of a 160-kDa group II intron ribozyme

We have reconstructed the group II intron from Pylaiella littoralis (PL) into a hydrolytic ribozyme, comprising domains 1-3 (D123) connected in cis plus domain 5 (D5) supplied in trans that efficiently cleaves spliced exon substrates. Using a novel gel-based fluorescence assay and nuclear magnetic resonance (NMR) spectroscopy, we monitored the direct binding of D5 to D123, characterized the kinetics of the spliced exon hydrolysis reaction (which is mechanistically analogous to the reverse of the second catalytic step of splicing), and identified the binding surface of D123 on D5. This PL ribozyme acts as an RNA endonuclease even at low monovalent (100 mM KCl) and divalent ion concentrations (1-10 mM MgCl(2)). This is in contrast to other group II intron ribozyme systems that require high levels of salt, making NMR analysis problematic. D5 binds tightly to D123 with a K(d) of 650 +/- 250 nM, a K(m) of approximately 300 nM, and a K(cat) of 0.02 min(-1) under single turnover conditions. Within the approximately 160-kDa D123-D5 binary complex, site-specific binding to D123 leads to dramatic chemical shift perturbation of residues localized to the tetraloop and internal bulge within D5, suggesting a structural switch model for D5-assisted splicing. This minimal ribozyme thus recapitulates the essential features of the reverse of the second catalytic step and represents a well-behaved system for ongoing high-resolution structural work to complement folding and catalytic functional studies.

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong 13C-13C magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to (13)C and (15)N resonances on the nucleotide base side chains. We report here the application of the alternate-site (13)C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing (13)C label in one of two patterns: Either C1' or C2' in addition to C4', termed (1'/2',4') labeling, or nearly complete labeling at the C2' and C4' sites only, termed (2',4') labeling. These patterns provide isolated 13C-1H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from 13C-13C scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to 13C-13C interaction in T (1) measurements of larger nucleic acids and in T (1rho) measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.

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.

Pseudo-CSA Restraints for NMR Refinement of Nucleic Acid Structure

Upon alignment of oligonucleotides in a magnetic field, the downfield TROSY component of the 13C-{1H} doublet changes its resonance frequency as a result of residual 13C-1H dipolar coupling (RDC) and residual 13C chemical shift anisotropy (RCSA), and the sum of these two second rank tensors is referred to as the pseudo-CSA. The experimentally measured difference in the resonance frequency of the 13C TROSY component in the aligned and isotropic samples is referred to as residual pseudo-CSA (RPCSA), and it can be used directly as a restraint during structure calculation. Because measurement of the RPCSA involves detection of the narrow TROSY 13C doublet component, it is applicable to systems with larger rotational correlation times than RDC measurement. The method is demonstrated for structure refinement of the helical region of a 24-nt stem-loop segment or ribosomal helix-35, uniformly enriched in 13C and 15N, with RPCSA values measured at 5 and 25 degrees C. Substantial cross-validated improvements in structural accuracy are obtained upon incorporation of RPCSA restraints.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

Induced fit and "lock and key" recognition of 5S RNA by zinc fingers of transcription factor IIIA

Transcription factor IIIA (TFIIIA) is a Cys2His2 zinc finger protein that regulates expression of the 5 S ribosomal RNA gene by binding specifically to the internal control element. TFIIIA also functions in transport and storage of 5 S RNA by binding directly to the RNA transcript. To obtain insights into the mechanism by which TFIIIA recognizes 5 S RNA, we determined the solution structure of the middle three zinc fingers bound to the central core of 5 S RNA. Finger 4 utilizes "lock and key" recognition to bind in the widened major groove of the pre-structured RNA loop E motif. This interaction is mediated by direct hydrogen bonding interactions with bases. In contrast, recognition of loop A, a flexible junction of three helices, occurs by an induced fit mechanism that involves reorganization of the conserved CAUA motif and structuring of the finger 5-finger 6 interface to form a complementary RNA binding surface.

The packaging signal of MLV is an integrated module that mediates intracellular transport of genomic RNAs

Packaging of MLV genomes requires four cis-acting stem-loops. Stem-loops A and B are self-complementary and bind Gag in their dimeric form, while the C and D elements mediate loop-loop interactions that facilitate RNA dimerization. Packaging also requires nuclear export of viral genomes, and their cytoplasmic transport toward the plasma membrane. For MLV, this is mediated by Gag and Env, and occurs on endosomal vesicles. Here, we report that MLV Psi acts at several steps during the transport of genomic RNAs. First, deletion of stem-loop B or C leads to the accumulation of genomic RNAs in the nucleus, suggesting that these elements are involved in export. Second, in chronically infected cells, mutation of the C and D loops impairs endosomal transport. This suggests that RNA dimerization is essential for vesicular transport, consistent with its proposed requirement for Gag binding. Surprisingly, deletion of stem-loop A blocks vesicular transport, whereas removal of stem-loop B has no effects. This suggests that stem-loop A has unique functions in packaging, not predicted from previous in vitro analyses. Finally, in packaging cells that do not express any Psi-containing RNA, endosomal RNA transport becomes sequence-independent. This non-specific activity of Gag likely promotes packaging of cellular mRNAs.

RNA structure: the long and the short of it

The database of RNA structure has grown tremendously since the crystal structure analyses of ribosomal subunits in 2000–2001. During the past year, the trend toward determining the structure of large, complex biological RNAs has accelerated, with the analysis of three intact group I introns, A- and B-type ribonuclease P RNAs, a riboswitch–substrate complex and other structures. The growing database of RNA structures, coupled with efforts directed at the standardization of nomenclature and classification of motifs, has resulted in the identification and characterization of numerous RNA secondary and tertiary structure motifs. Because a large proportion of RNA structure can now be shown to be composed of these recurring structural motifs, a view of RNA as a modular structure built from a combination of these building blocks and tertiary linkers is beginning to emerge. At the same time, however, more detailed analysis of water, metal, ligand and protein binding to RNA is revealing the effect of these moieties on folding and structure formation. The balance between the views of RNA structure either as strictly a construct of preformed building blocks linked in a limited number of ways or as a flexible polymer assuming a global fold influenced by its environment will be the focus of current and future RNA structural biology.

RNA helical packing in solution: NMR structure of a 30 kDa GAAA tetraloop-receptor complex

Tertiary interactions are critical for proper RNA folding and ribozyme catalysis. RNA tertiary structure is often condensed through long-range helical packing interactions mediated by loop-receptor motifs. RNA structures displaying helical packing by loop-receptor interactions have been solved by X-ray crystallography, but not by NMR. Here, we report the NMR structure of a 30 kDa GAAA tetraloop-receptor RNA complex. In order to stabilize the complex, we used a modular design in which the RNA was engineered to form a homodimer, with each subunit containing a GAAA tetraloop phased one helical turn apart from its cognate 11-nucleotide receptor domain. The structure determination utilized specific isotopic labeling patterns (2H, 13C and 15N) and refinement against residual dipolar couplings. We observe a unique and highly unusual chemical shift pattern for an adenosine platform interaction that reveals a spectroscopic fingerprint for this motif. The structure of the GAAA tetraloop-receptor interaction is well defined solely from experimental NMR data, shows minor deviations from previously solved crystal structures, and verifies the previously inferred hydrogen bonding patterns within this motif. This work demonstrates the feasibility of using engineered homodimers as modular systems for the determination of RNA tertiary interactions by NMR.

How retroviruses select their genomes

As retroviruses assemble in infected cells, two copies of their full-length, unspliced RNA genomes are selected for packaging from a cellular milieu that contains a substantial excess of non-viral and spliced viral RNAs. Understanding the molecular details of genome packaging is important for the development of new antiviral strategies and to enhance the efficacy of retroviral vectors used in human gene therapy. Recent studies of viral RNA structure in vitro and in vivo and high-resolution studies of RNA fragments and protein-RNA complexes are helping to unravel the mechanism of genome packaging and providing the first glimpses of the initial stages of retrovirus assembly.

Composition and sequence-dependent binding of RNA to the nucleocapsid protein of Moloney murine leukemia virus

All retroviruses package two copies of their genomes during virus assembly, both of which are required for strand transfer-mediated recombination during reverse transcription. Genome packaging is mediated by interactions between the nucleocapsid (NC) domains of assembling Gag polyproteins and an RNA packaging signal, located near the 5' end of the genome, called Psi. We recently discovered that the NC protein of the Moloney murine leukemia virus (MLV) can bind with high affinity to conserved UCUG elements within the MLV packaging signal [D'Souza, V., and Summers, M. F. (2004) Nature 431, 586-590]. Selective binding to dimeric RNA is regulated by a conformational RNA switch, in which the UCUG elements are sequestered by base pairing in the monomeric RNA and do not bind NC, but become exposed for NC binding upon dimerization. Dimerization-dependent structural changes occur in other regions of the Psi-site, exposing guanosine-containing segments that might also bind NC. Here we demonstrate that short RNAs containing three such sequences, ACAG, UUUG, and UCCG, can bind NC with significant affinity (K(d) = 94-315 nM). Titration experiments with oligoribonucleotides of varying lengths and compositions, combined with NMR-based structural studies, reveal that binding is strictly dependent on the presence of an unpaired guanosine, and that relative binding affinities can vary by more than 1 order of magnitude depending on the nature of the three upstream nucleotides. Binding is enhanced in short RNAs containing terminal phosphates, indicating that electrostatic interactions contribute significantly to binding. Our findings extend a previously published model for genome recognition, in which the NC domains of assembling Gag molecules interact with multiple X(i-3)-X(i-2)-X(i-1)-G(i) elements (X is a variable nucleotide) that appear to be preferentially exposed in the dimeric RNA.

RNA flexibility in the dimerization domain of a gamma retrovirus

Retroviruses are the causative agents of serious diseases, such as acquired immunodeficiency syndromes and several cancers, and are also useful gene therapy vectors. Retroviruses contain two sense-strand RNA genomes, which become linked at their 5' ends to form an RNA dimer. Understanding the molecular basis for dimerization may yield new approaches for controlling viral infectivity. Because this RNA domain is highly conserved within retrovirus groups, it has not been possible to define a consensus structure for the 5' dimerization domain by comparative sequence analysis. Here, we defined a 170-nucleotide minimal dimerization active sequence (MiDAS) for a representative gamma retrovirus, the Moloney murine sarcoma virus, by stringent competitive dimerization. We then analyzed the structure at every nucleotide in the MiDAS monomeric starting state with quantitative selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry. Notably, SHAPE analysis demonstrated that the RNA monomer contains an extensive flexible domain spanning 50 nucleotides. These findings support a structural model in which RNA flexibility directly facilitates retroviral genome dimerization by reducing the energetic cost of disrupting pre-existing base pairings in the monomer.

Accurate measurement of 15N-13C residual dipolar couplings in nucleic acids

New 3D HCN quantitative J (QJ) pulse schemes are presented for the precise and accurate measurement of one-bond 15N1/9-13C1', 15N1/9-13C6/8, and 15N1/9-13C2/4 residual dipolar couplings (RDCs) in weakly aligned nucleic acids. The methods employ 1H-13C multiple quantum (MQ) coherence or TROSY-type pulse sequences for optimal resolution and sensitivity. RDCs are obtained from the intensity ratio of H1'-C1'-N1/9 (MQ-HCN-QJ) or H6/8-C6/8-N1/9 (TROSY-HCN-QJ) correlations in two interleaved 3D NMR spectra, with dephasing intervals of zero (reference spectrum) and approximately 1/(2J(NC)) (attenuated spectrum). The different types of 15N-13C couplings can be obtained by using either the 3D MQ-HCN-QJ or TROSY-HCN-QJ pulse scheme, with the appropriate setting of the duration of the constant-time 15N evolution period and the offset of two frequency-selective 13C pulses. The methods are demonstrated for a uniformly 13C, 15N-enriched 24-nucleotide stem-loop RNA sequence, helix-35psi, aligned in the magnetic field using phage Pf1. For measurements of RDCs systematic errors are found to be negligible, and experiments performed on a 1.5 mM helix-35psi sample result in an estimated precision of ca. 0.07 Hz for 1D(NC), indicating the utility of the measured RDCs in structure validation and refinement. Indeed, for a complete set of 15N1/9-13C1', 15N1/9-13C6/8, and 15N1/9-13C2/4 dipolar couplings obtained for the stem nucleotides, the measured RDCs are in excellent agreement with those predicted for an NMR structure of helix-35psi, refined using independently measured observables, including 13C-1H, 13C-13C and 1H-1H dipolar couplings.

An intermolecular base triple as the basis of ligand specificity and affinity in the guanine- and adenine-sensing riboswitch RNAs

Riboswitches are highly structured RNA elements that control the expression of many bacterial genes by binding directly to small metabolite molecules with high specificity and affinity. In Bacillus subtilis, two classes of riboswitches have been described that discriminate between guanine and adenine despite an extremely high degree of homology both in their primary and secondary structure. We have identified intermolecular base triples between both purine ligands and their respective riboswitch RNAs by NMR spectroscopy. Here, specificity is mediated by the formation of a Watson-Crick base pair between the guanine ligand and a C residue or the adenine ligand and a U residue of the cognate riboswitch RNA, respectively. In addition, a second base-pairing interaction common to both riboswitch purine complexes involves a uridine residue of the RNA and the N3/N9 edge of the purine ligands. This base pairing is mediated by a previously undescribed hydrogen-bonding scheme that contributes to the affinity of the RNA-ligand interaction. The observed intermolecular hydrogen bonds between the purine ligands and the RNA rationalize the previously observed change in specificity upon a C to U mutation in the core of the purine riboswitch RNAs and the differences in the binding affinities for a number of purine analogs.

Structure Determination of Protein⧸RNA Complexes by NMR

Structure determination of protein?RNA complexes in solution provides unique insights into factors that are involved in protein/RNA recognition. Here, we review the methodology used in our laboratory to overcome the challenges of protein?RNA structure determination by nuclear magnetic resonance (NMR). We use as two examples complexes recently solved in our laboratory, the nucleolin RBD12/b2NRE and Rnt1p dsRBD/snR47h complexes. Topics covered are protein and RNA preparation, complex formation, identification of the protein/RNA interface, protein and RNA resonance assignment, intermolecular NOE assignment, and structure calculation and analysis.

Application of NMR and EPR methods to the study of RNA

The application of techniques based on magnetic resonance, specifically electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR), has provided a wealth of new information on RNA structures, as well as insights into the dynamics and function of these important biomolecules. NMR spectroscopy is very successful for determining the solution structures of small RNA domains, aptamers and ribozymes, and exploring their intramolecular dynamics and interactions with ligands. EPR-based methods have been used to map local dynamic and structural features of RNA, to explore different modes of RNA-ligand interaction, to obtain long-range structural restraints and to probe metal-ion-binding sites.

Changes in DNA bending induced by restricting nucleotide ring pucker studied by weak alignment NMR spectroscopy

Changes in bending of the DNA helix axis caused by the introduction of conformationally locked nucleotide analogs into the center region of the palindromic Dickerson dodecamer, d(CGCGAATTCGCG)(2), have been studied by NMR measurement of residual one-bond (13)C-(1)H dipolar couplings. Thymidine analogs, in which the deoxyribose was substituted by bicyclo[3.1.0]hexane, were incorporated in the T7, T8, and T7T8 positions. These nucleotide analogs restrict the ring pucker to the C2'-exo or "north" conformation, instead of C2'-endo or "south," which dominates in regular B-form DNA. For all three oligomers, bending toward the major groove is found relative to the native molecule. The effects are additive with bending of 5 +/- 1 degrees per locked nucleotide. Measurement of the change in bending is more accurate than measurement of the bending angle itself and requires far fewer experimental data.