NMR solution structure of the lead-dependent ribozyme: evidence for dynamics in RNA catalysis

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Single-round deoxyribozyme discovery

Artificial evolution experiments typically use libraries of ∼1015 sequences and require multiple rounds of selection to identify rare variants with a desired activity. Based on the simple structures of some aptamers and nucleic acid enzymes, we hypothesized that functional motifs could be isolated from significantly smaller libraries in a single round of selection followed by high-throughput sequencing. To test this idea, we investigated the catalytic potential of DNA architectures in which twelve or fifteen randomized positions were embedded in a scaffold present in all library members. After incubating in either the presence or absence of lead (which promotes the nonenzymatic cleavage of RNA), library members that cleaved themselves at an RNA linkage were purified by PAGE and characterized by high-throughput sequencing. These selections yielded deoxyribozymes with activities 8- to 30-fold lower than those previously isolated under similar conditions from libraries containing 1014 different sequences, indicating that the disadvantage of using a less diverse pool can be surprisingly small. It was also possible to elucidate the sequence requirements and secondary structures of deoxyribozymes without performing additional experiments. Due to its relative simplicity, we anticipate that this approach will accelerate the discovery of new catalytic DNA and RNA motifs.

Dynamics-Function Analysis in Catalytic RNA Using NMR Spin Relaxation and Conformationally Restricted Nucleotides

A full understanding of biomolecular function requires an analysis of both the dynamic properties of the system of interest and the identification of those dynamics that are required for function. We describe NMR methods based on metabolically directed specific isotope labeling for the identification of molecular disorder and/or conformational transitions on the RNA backbone ribose groups. These analyses are complemented by the use of synthetic covalently modified nucleotides constrained to a single sugar pucker, which allow functional assessment of dynamics by selectively removing a minor conformer identified by NMR from the structural ensemble.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme

In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.

Crystal structure of an RNA-cleaving DNAzyme

In addition to storage of genetic information, DNA can also catalyze various reactions. RNA-cleaving DNAzymes are the catalytic DNAs discovered the earliest, and they can cleave RNAs in a sequence-specific manner. Owing to their great potential in medical therapeutics, virus control, and gene silencing for disease treatments, RNA-cleaving DNAzymes have been extensively studied; however, the mechanistic understandings of their substrate recognition and catalysis remain elusive. Here, we report three catalytic form 8-17 DNAzyme crystal structures. 8-17 DNAzyme adopts a V-shape fold, and the Pb²⁺ cofactor is bound at the pre-organized pocket. The structures with Pb²⁺ and the modification at the cleavage site captured the pre-catalytic state of the RNA cleavage reaction, illustrating the unexpected Pb²⁺-accelerated catalysis, intrinsic tertiary interactions, and molecular kink at the active site. Our studies reveal that DNA is capable of forming a compacted structure and that the functionality-limited bio-polymer can have a novel solution for a functional need in catalysis.

Finite-size effect on the charging free energy of protein in explicit solvent

The finite-size effect in periodic system is examined for the charging free energy of protein in explicit solvent over a variety of charged states. The key to the finite-size correction is the self-energy, which is defined as the interaction energy of the solute with its own periodic images and the neutralizing background. By employing the thermodynamic-integration method with systematically varied sizes of the unit cell of molecular dynamics (MD) simulations, we show for ubiquitin that the self-energy corrects the finite-size effect on the charging free energy within 1 kcal/mol at total charges of -5e, -1e, neutral, and +1e and within 5 kcal/mol even for a highly charged state with +8e. We then sought the additional correction from the solvation effect using the numerical solution to the Poisson equation of the protein with implicit solvent. This correction reduces the cell-size dependence of the charging free energy at +8e to 3 kcal/mol and is well expressed as the self-energy divided by the dielectric constant of solvent water.

The Conservation and Function of RNA Secondary Structure in Plants

RNA transcripts fold into secondary structures via intricate patterns of base pairing. These secondary structures impart catalytic, ligand binding, and scaffolding functions to a wide array of RNAs, forming a critical node of biological regulation. Among their many functions, RNA structural elements modulate epigenetic marks, alter mRNA stability and translation, regulate alternative splicing, transduce signals, and scaffold large macromolecular complexes. Thus, the study of RNA secondary structure is critical to understanding the function and regulation of RNA transcripts. Here, we review the origins, form, and function of RNA secondary structure, focusing on plants. We then provide an overview of methods for probing secondary structure, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods. Combining these latter methods with high-throughput sequencing has enabled them to scale across whole transcriptomes, yielding tremendous new insights into the form and function of RNA secondary structure.

Pf1 filamentous phage as an alignment tool for generating local and global structural information in nucleic acids

Abstract Pf1 filamentous phage represent a simple versatile method for generating partially ordered macromolecules in solution. The phage allow tunable degrees of alignment of macromolecules under a wide range of temperature and solvent conditions. The negatively charged phage are ideal for aligning negatively charged nucleic acids and these phage-nucleic acid solutions are stable indefinitely. We have used Pf1 phage to align various DNA and RNA molecules in solution for measurement of dipolar coupling interactions. These dipolar couplings can be used to improve the local structure of nucleic acids. More importantly they also contain information on the global structure, such as DNA bending, which presently cannot be obtained by standard NMR methods. The principles involved in using Pf1 phage to generate solutions of partially order macromolecules will be discussed. The use of (1)H-(1)H, (1)H-(13)C and (1)H-(15)N dipolar couplings for generating angle constraints for structure refinement of nucleic acids will also be discussed.

Searching for a DNAzyme Version of the Leadzyme

The leadzyme refers to a small ribozyme that cleaves a RNA substrate in the presence of Pb(2+). In an optimized form, the enzyme strand contains only two unpaired nucleotides. Most RNA-cleaving DNAzymes are much longer. Two classical Pb(2+)-dependent DNAzymes, 8-17 and GR5, both contain around 15 nucleotides in the enzyme loop. This is also the size of most RNA-cleaving DNAzymes that use other metal ions for their activity. Such large enzyme loops make spectroscopic characterization difficult and so far no high-resolution structural information is available for active DNAzymes. The goal of this work is to search for DNAzymes with smaller enzyme loops. A simple replacement of the ribonucleotides in the leadzyme by deoxyribonucleotides failed to produce an active enzyme. A Pb(2+)-dependent in vitro selection combined with deep sequencing was then performed. After sequence alignment and DNA folding, a new DNAzyme named PbE22 was identified, which contains only 5 nucleotides in the enzyme catalytic loop. The biochemical characteristics of PbE22 were compared with those of the leadzyme and the two classical Pb(2+)-dependent DNAzymes. The rate of PbE22 rises with increase in Pb(2+) concentration, being 1.7 h(-1) in the presence of 100 μM Pb(2+) and reaching 3.5 h(-1) at 500 µM Pb(2+). The log of PbE22 rate rises linearly in a pH-dependent fashion (20 µM Pb(2+)) with a slope of 0.74. In addition, many other abundant sequences in the final library were studied. These sequences are quite varied in length and nucleotide composition, but some contain a few conserved nucleotides consistent with the GR5 structure. Interestingly, some sequences are active with Pb(2+) but none of them were active with even 50 mM Mg(2+), which is reminiscent of the difference between the GR5 and 8-17 DNAzymes.

Structure, dynamics, and mechanism of the lead-dependent ribozyme

Leadzyme is a small catalytic RNA that was identified by in vitro selection for Pb2+-dependent cleavage from a tRNA library. Leadzyme employs a unique two-step Pb2+-specific mechanism to cleave within its active site. NMR and crystal structures of the active site revealed different folding patterns, but neither features the in-line alignment for attack by the 2'-OH nucleophilic group. These experimentally determined structures most likely represent ground states and are catalytically inactive. There are significant dynamics of the active site and the motif samples multiple conformations at the ground states. Various metal ion binding sites have been identified, including one that may be occupied by a catalytic Pb2+. Based on functional group analysis, a computational model of the transition state has been proposed. This model features a unique base triple that is consistent with sequence and functional group requirements for catalysis. This structure is likely only populated transiently, but imposing appropriate conformational constraints may significantly stabilize this state thereby promoting catalysis. Other ions may inhibit the cleavage by competing for the Pb2+ binding site, or by stabilizing the ground state thereby suppressing its transition to the catalytically active conformation. Some rare earth ions can enhance the reaction via an unknown mechanism. Because of its unique chemistry and dynamic behavior, leadzyme can continue to serve as an excellent model system for teaching us RNA biology and chemistry.

Role of helical constraints of the EBS1-IBS1 duplex of a group II intron on demarcation of the 5' splice site

Recognition of the 5' splice site by group II introns involves pairing between an exon binding sequence (EBS) 1 within the ID3 stem-loop of domain 1 and a complementary sequence at the 3' end of exon 1 (IBS1). To identify the molecular basis for splice site definition of a group IIB ai5γ intron, we probed the solution structure of the ID3 stem-loop alone and upon binding of its IBS1 target by solution NMR. The ID3 stem was structured. The base of the ID3 loop was stacked but displayed a highly flexible EBS1 region. The flexibility of EBS1 appears to be a general feature of the ai5γ and the smaller Oceanobacillus iheyensis (O.i.) intron and may help in effective search of conformational space and prevent errors in splicing as a result of fortuitous base-pairing. Binding of IBS1 results in formation of a structured seven base pair duplex that terminates at the 5' splice site in spite of the potential for additional A-U and G•U pairs. Comparison of these data with conformational features of EBS1-IBS1 duplexes extracted from published structures suggests that termination of the duplex and definition of the splice site are governed by constraints of the helical geometry within the ID3 loop. This feature and flexibility of the uncomplexed ID3 loop appear to be common for both the ai5γ and O.i. introns and may help to fine-tune elements of recognition in group II introns.

Solution structure of a DNA mimicking motif of an RNA aptamer against transcription factor AML1 Runt domain

AML1/RUNX1 is an essential transcription factor involved in the differentiation of hematopoietic cells. AML1 binds to the Runt-binding double-stranded DNA element (RDE) of target genes through its N-terminal Runt domain. In a previous study, we obtained RNA aptamers against the AML1 Runt domain by systematic evolution of ligands by exponential enrichment and revealed that RNA aptamers exhibit higher affinity for the Runt domain than that for RDE and possess the 5'-GCGMGNN-3' and 5'-N'N'CCAC-3' conserved motif (M: A or C; N and N' form Watson-Crick base pairs) that is important for Runt domain binding. In this study, to understand the structural basis of recognition of the Runt domain by the aptamer motif, the solution structure of a 22-mer RNA was determined using nuclear magnetic resonance. The motif contains the AH(+)-C mismatch and base triple and adopts an unusual backbone structure. Structural analysis of the aptamer motif indicated that the aptamer binds to the Runt domain by mimicking the RDE sequence and structure. Our data should enhance the understanding of the structural basis of DNA mimicry by RNA molecules.

Towards Accurate Prediction of Protonation Equilibrium of Nucleic Acids

The role of protonated nucleotides in modulating the pH-dependent properties of nucleic acids is one of the emerging frontiers in the field of nucleic acid biology. The recent development of a constant pH molecular dynamics simulation (CPHMD^MSλD) framework for simulating nucleic acids has provided a tool for realistic simulations of pH-dependent dynamics. We enhanced the CPHMD^MSλD framework with pH-based replica exchange (pH-REX), which significantly improves the sampling of both titration and spatial coordinates. The results from our pK_a calculations for the GAAA tetraloop, which was predicted with lower accuracy previously due to sampling challenges, demonstrates that pH-REX reduces the average unsigned error (AUE) to 0.7 pK_a units, and the error of the most poorly predicted residue A17 was drastically reduced from 2.9 to 1.2 pK_a unit. Lastly, we show that pH-REX CPHMD^MSλD simulations can be used to identify the dominant conformation of nucleic acid structures in alternate pH environments. This work suggests that pH-REX CPHMD^MSλD simulations provide a practical tool for predicting nucleic acid protonation equilibrium from first-principles, and offering structural and mechanistic insight into the study of pH-dependent properties of nucleic acids.

pH-dependent dynamics of complex RNA macromolecules

The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMSλD) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMSλD to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMSλD accurately predicts the direction of the pKa shifts and reproduces experimentally-measured microscopic pKa values with an average unsigned error of 1.3 pKa units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMSλD simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provides a powerful tool to investigate pH-dependent processes in nucleic acids.

RNA-PAIRS: RNA probabilistic assignment of imino resonance shifts

The significant biological role of RNA has further highlighted the need for improving the accuracy, efficiency and the reach of methods for investigating RNA structure and function. Nuclear magnetic resonance (NMR) spectroscopy is vital to furthering the goals of RNA structural biology because of its distinctive capabilities. However, the dispersion pattern in the NMR spectra of RNA makes automated resonance assignment, a key step in NMR investigation of biomolecules, remarkably challenging. Herein we present RNA Probabilistic Assignment of Imino Resonance Shifts (RNA-PAIRS), a method for the automated assignment of RNA imino resonances with synchronized verification and correction of predicted secondary structure. RNA-PAIRS represents an advance in modeling the assignment paradigm because it seeds the probabilistic network for assignment with experimental NMR data, and predicted RNA secondary structure, simultaneously and from the start. Subsequently, RNA-PAIRS sets in motion a dynamic network that reverberates between predictions and experimental evidence in order to reconcile and rectify resonance assignments and secondary structure information. The procedure is halted when assignments and base-parings are deemed to be most consistent with observed crosspeaks. The current implementation of RNA-PAIRS uses an initial peak list derived from proton-nitrogen heteronuclear multiple quantum correlation ((1)H-(15)N 2D HMQC) and proton-proton nuclear Overhauser enhancement spectroscopy ((1)H-(1)H 2D NOESY) experiments. We have evaluated the performance of RNA-PAIRS by using it to analyze NMR datasets from 26 previously studied RNAs, including a 111-nucleotide complex. For moderately sized RNA molecules, and over a range of comparatively complex structural motifs, the average assignment accuracy exceeds 90%, while the average base pair prediction accuracy exceeded 93%. RNA-PAIRS yielded accurate assignments and base pairings consistent with imino resonances for a majority of the NMR resonances, even when the initial predictions are only modestly accurate. RNA-PAIRS is available as a public web-server at http://pine.nmrfam.wisc.edu/RNA/.

Prevalence of syn nucleobases in the active sites of functional RNAs

Biological RNAs, like their DNA counterparts, contain helical stretches, which have standard Watson-Crick base pairs in the anti conformation. Most functional RNAs also adopt geometries with far greater complexity such as bulges, loops, and multihelical junctions. Occasionally, nucleobases in these regions populate the syn conformation wherein the base resides close to or over the ribose sugar, which leads to a more compact state. The importance of the syn conformation to RNA function is largely unknown. In this study, we analyze 51 RNAs with tertiary structure, including aptamers, riboswitches, ribozymes, and ribosomal RNAs, for number, location, and properties of syn nucleobases. These RNAs represent the set of nonoverlapping, moderate- to high-resolution structures available at present. We find that syn nucleobases are much more common among purines than pyrimidines, and that they favor C2'-endo-like conformations especially among those nucleobases in the intermediate syn conformation. Strikingly, most syn nucleobases participate in tertiary stacking and base-pairing interactions: Inspection of RNA structures revealed that the majority of the syn nucleobases are in regions assigned to function, with many syn nucleobases interacting directly with a ligand or ribozyme active site. These observations suggest that judicious placement of conformationally restricted nucleotides biased into the syn conformation could enhance RNA folding and catalysis. Such changes could also be useful for locking RNAs into functionally competent folds for use in X-ray crystallography and NMR.

A CA(+) pair adjacent to a sheared GA or AA pair stabilizes size-symmetric RNA internal loops

RNA internal loops are often important sites for folding and function. Residues in internal loops can have pK_a values shifted close to neutral pH because of the local structural environment. A series of RNA internal loops were studied at different pH by UV absorbance versus temperature melting experiments and imino proton nuclear magnetic resonance (NMR). A stabilizing CA pair forms at pH 7 in the and nearest neighbors when the CA pair is the first noncanonical pair (loop-terminal pair) in 3 × 3 nucleotide and larger size-symmetric internal loops. These and nearest neighbors, with CA adjacent to a closing Watson−Crick pair, are further stabilized when the pH is lowered from 7 to 5.5. The results are consistent with a significantly larger fraction (from ∼20% at pH 7 to ∼90% at pH 5.5) of adenines being protonated at the N1 position to form stabilizing wobble CA⁺ pairs adjacent to a sheared GA or AA pair. The noncanonical pair adjacent to the GA pair in can either stabilize or destabilize the loop, consistent with the sequence-dependent thermodynamics of GA pairs. No significant pH-dependent stabilization is found for most of the other nearest neighbor combinations involving CA pairs (e.g., and ), which is consistent with the formation of various nonwobble pairs observed in different local sequence contexts in crystal and NMR structures. A revised free-energy model, including stabilization by wobble CA⁺ pairs, is derived for predicting stabilities of medium-size RNA internal loops.

Conformational distribution and ultrafast base dynamics of leadzyme

The dynamic nature of ribozymes represents a significant challenge in elucidating their structure-dynamics-function relationship. Here, using femtosecond time-resolved spectroscopy and other biophysical tools, we demonstrate that the active site of leadzyme does not have a unique structure, but rather samples an ensemble of conformations that undergo picosecond structural changes. Various base modifications have a profound context-dependent impact on the catalysis.

A catalytic metal ion interacts with the cleavage Site G.U wobble in the HDV ribozyme

The HDV ribozyme self-cleaves by a chemical mechanism involving general acid-base catalysis to generate 2',3'-cyclic phosphate and 5'-hydroxyl termini. Biochemical studies from several laboratories have implicated C75 as the general acid and hydrated magnesium as the general base. We have previously shown that C75 has a pK(a) shifted >2 pH units toward neutrality [Gong, B., Chen, J. H., Chase, E., Chadalavada, D. M., Yajima, R., Golden, B. L., Bevilacqua, P. C., and Carey, P. R. (2007) J. Am. Chem. Soc. 129, 13335-13342], while in crystal structures, it is well-positioned for proton transfer. However, no evidence for a hydrated magnesium poised to serve as a general base in the reaction has been observed in high-resolution crystal structures of various reaction states and mutants. Herein, we use solution kinetic experiments and parallel Raman crystallographic studies to examine the effects of pH on the rate and Mg(2+) binding properties of wild-type and 7-deazaguanosine mutants of the HDV ribozyme. These data suggest that a previously unobserved hydrated magnesium ion interacts with N7 of the cleavage site G.U wobble base pair. Integrating this metal ion binding site with the available crystal structures provides a new three-dimensional model for the active site of the ribozyme that accommodates all available biochemical data and appears competent for catalysis. The position of this metal is consistent with a role of a magnesium-bound hydroxide as a general base as dictated by biochemical data.

Computational strategies for the automated design of RNA nanoscale structures from building blocks using NanoTiler

One approach to designing RNA nanoscale structures is to use known RNA structural motifs such as junctions, kissing loops or bulges and to construct a molecular model by connecting these building blocks with helical struts. We previously developed an algorithm for detecting internal loops, junctions and kissing loops in RNA structures. Here we present algorithms for automating or assisting many of the steps that are involved in creating RNA structures from building blocks: (1) assembling building blocks into nanostructures using either a combinatorial search or constraint satisfaction; (2) optimizing RNA 3D ring structures to improve ring closure; (3) sequence optimisation; (4) creating a unique non-degenerate RNA topology descriptor. This effectively creates a computational pipeline for generating molecular models of RNA nanostructures and more specifically RNA ring structures with optimized sequences from RNA building blocks. We show several examples of how the algorithms can be utilized to generate RNA tecto-shapes.

Conformationally restricted nucleotides as a probe of structure-function relationships in RNA

We introduce the use of commercially available locked nucleic acids (LNAs) as a functional probe in RNA. LNA nucleotides contain a covalent linkage that restricts the pseudorotation phase of the ribose to C3'-endo (A-form). Introduction of an LNA at a single site thus allows the role of ribose structure and dynamics in RNA function to be assessed. We apply LNA probing at multiple sites to analyze self-cleavage in the lead-dependent ribozyme (leadzyme), thermodynamic stability in the UUCG tetraloop, and the kinetics of recognition of U1A protein by U1 snRNA hairpin II. In the leadzyme, locking a single guanosine residue into the C3'-endo pucker increases the catalytic rate by a factor of 20, despite the fact that X-ray crystallographic and NMR structures of the leadzyme ground state reported a C2'-endo conformation at this site. These results strongly suggest that a conformational change at this position is critical for catalytic function. Functional insights obtained in all three systems demonstrate the highly general applicability of LNA probing in analysis of the role of ribose orientation in RNA structure, dynamics, and function.

A conformationally restricted guanosine analog reveals the catalytic relevance of three structures of an RNA enzyme

Recent studies indicate that RNA function can be enhanced by the incorporation of conformationally restricted nucleotides. Herein, we use 8-bromoguanosine, a nucleotide analog with an enforced syn conformation, to elucidate the catalytic relevance of ribozyme structures. We chose to study the lead-dependent ribozyme (leadzyme) because structural models derived from NMR, crystal, and computational (MC-Sym) studies differ in which of the three active site guanosines (G7, G9, or G24) have a syn glycosidic torsion angle. Kinetic assays were carried out on 8BrG variants at these three guanosine positions. These data indicate that an 8BrG24 leadzyme is hyperactive, while 8BrG7 and 8BrG9 leadzymes have reduced activity. These findings support the computational model of the leadzyme, rather than the NMR and crystal structures, as being the most relevant to phosphodiester bond cleavage.

NMR (13)C-relaxation study of base and sugar dynamics in GCAA RNA hairpin tetraloop

Intramolecular dynamics of a 14-mer RNA hairpin including GCAA tetraloop was investigated by (13)C NMR relaxation. R(1) and R(1rho) relaxation rates were measured for all protonated base carbons as well as for C1' carbons of ribose sugars at several magnetic field strengths. The data has been interpreted in the framework of modelfree analysis [G. Lipari and A. Szabo. J Am Chem Soc 104, 4546-4559 (1982); G. Lipari and A. Szabo. J Am Chem Soc 104, 4559-4570 (1982)] characterizing the internal dynamics of the molecule by order parameters and correlation times for fast motions on picosecond to nanosecond time scale and by contributions of the chemical exchange. The fast dynamics reveals a rather rigid stem and a significantly more flexible loop. The cytosine and the last adenine bases in the loop as well as all the loop sugars exhibit a significant contribution of conformational equilibrium on microsecond to millisecond time scale. The high R(1rho) values detected on both base and sugar moieties of the loop indicate coordinated motions in this region. A semiquantitative analysis of the conformational equilibrium suggests the exchange rates on the order of 10(4) s(-1). The results are in general agreement with dynamics studies of GAAA loops by NMR relaxation and fluorescent spectroscopy and support the data on the GCAA loop dynamics obtained by MD simulations.

NMR structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: probing the structural features that shape the thermodynamic stability of GA pairs

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.

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.

A suite of 2H NMR spin relaxation experiments for the measurement of RNA dynamics

A suite of (2)H-based spin relaxation NMR experiments is presented for the measurement of molecular dynamics in a site-specific manner in uniformly (13)C, randomly fractionally deuterated ( approximately 50%) RNA molecules. The experiments quantify (2)H R(1) and R(2) relaxation rates that can subsequently be analyzed to obtain information about dynamics on a pico- to nanosecond time scale. Sensitivity permitting, the consistency of the data can be evaluated by measuring all five rates that are accessible for a spin 1 particle and establishing that the rates obey relations that are predicted from theory. The utility of the methodology is demonstrated with studies of the dynamics of a 14-mer RNA containing the UUCG tetraloop at temperatures of 25 and 5 degrees C. The high quality of the data, even at 5 degrees C, suggests that the experiments will be of use for the study of RNA molecules that are as large as 30 nucleotides.

Linkage between proton binding and folding in RNA: a thermodynamic framework and its experimental application for investigating pKa shifting

Perturbation of pKa values can change the favored protonation states of the nucleobases at biological pH and thereby modulate the function of RNA and DNA molecules. In an effort to understand the driving forces for pKa shifting specific to nucleic acids, we developed a thermodynamic framework that relates proton binding to the nucleobases and the helix-coil transition. Key features that emerge from the treatment are a comprehensive description of all the actions of proton binding on RNA folding: acid and alkaline denaturation of the helix and pKa shifting in the folded state. Practical experimental approaches for measuring pKas from thermal denaturation experiments are developed. Microscopic pka values (where ka is the acid dissociation constant) for the unfolded state were determined directly by experiments on unstructured oligonucleotides, which led to a macroscopic pKa for the ensemble of unfolded states shifted toward neutrality. The formalism was then applied to pH-dependent UV melting data for model DNA oligonucleotides. Folded-state pka) values were in good agreement with the outcome of pH titrations, and the acid and alkaline denaturation regions were well described. The formalism developed here is similar to that of Draper and coworkers for Mg2+ binding to RNA, except that the unfolded state is described explicitly owing to the presence of specific proton-binding sites on the bases. A principal conclusion is that it should be possible to attain large pKa shifts by designing RNA molecules that fold cooperatively.

Dynamics in the U6 RNA intramolecular stem-loop: a base flipping conformational change

The U6 RNA intramolecular stem-loop (ISL) structure is an essential component of the spliceosome and binds a metal ion required for pre-messenger RNA splicing. The metal binding internal loop region of the stem contains a partially protonated C67-(+)A79 base pair (pK(a) = 6.5) and an unpaired U80 nucleotide that is stacked within the helix at pH 7.0. Here, we determine that protonation occurs with an exchange lifetime of approximately 20 micros and report the solution structures of the U6 ISL at pH 5.7. The differences between pH 5.7 and 7.0 structures reveal that the pH change significantly alters the RNA conformation. At lower pH, U80 is flipped out into the major groove. Base flipping involves a purine stacking interaction of flanking nucleotides, inversion of the sugar pucker 5' to the flipped base, and phosphodiester backbone rearrangement. Analysis of residual dipolar couplings as a function of pH indicates that base flipping is not restricted to a local conformational change. Rather, base flipping alters the alignment of the upper and lower helices. The alternative conformations of the U6 ISL reveal striking structural similarities with both the NMR and crystal structures of domain 5 of self-splicing group II introns. These structures suggest that base flipping at an essential metal binding site is a conserved feature of the splicing machinery for both the spliceosome and group II self-splicing introns.

The solution structure of the VS ribozyme active site loop reveals a dynamic "hot-spot"

The VS ribozyme is the largest ribozyme in its class and is also the least structurally characterized thus far. The current working model of the VS ribozyme locates the active site in stem-loop VI. The solution structure of this active site loop was determined using high resolution NMR spectroscopy. The structure reveals that the ground-state conformation of the active site differs significantly from that determined previously from chemical structure probing and mutational analysis of the ribozyme in its active conformation, which contains several looped out bases. In contrast, the base-pairing scheme found for the isolated loop contains three mismatched base-pairs: an A+-C, a G-U wobble, and a sheared G-A base-pair and no looped out bases. Dynamics observed within the active site loop provide insight into the mechanism by which the RNA can rearrange its secondary structure into an "activated" conformation prior to cleavage. These findings lend support to the idea that RNA secondary structure is more fluid than once believed and that a better understanding of structure and dynamic features of ribozymes is required to unravel the intricacies of their catalytic abilities.

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

The full length, positive-strand genome of the Moloney Murine Leukemia Virus contains a "core encapsidation signal" that is essential for efficient genome packaging during virus assembly. We have determined the structure of a 101-nucleotide RNA that contains this signal (called mPsi) using a novel isotope-edited NMR approach. The method is robust and should be generally applicable to larger RNAs. mPsi folds into three stem loops, two of which (SL-C and SL-D) co-stack to form an extended helix. The third stem loop (SL-B) is connected to SL-C by a flexible, four-nucleotide linker. The structure contains five mismatched base-pairs, an unusual C.CG base-triple platform, and a novel "A-minor K-turn," in which unpaired adenosine bases A340 and A341 of a GGAA bulge pack in the minor groove of a proximal stem, and a bulged distal uridine (U319) forms a hydrogen bond with the phosphodiester of A341. Phylogenetic analyses indicate that these essential structural elements are conserved among the murine C-type retroviruses.

Crystal structure of the leadzyme at 1.8 A resolution: metal ion binding and the implications for catalytic mechanism and allo site ion regulation

The leadzyme is a small ribozyme, derived from in vitro selection, which catalyzes site specific, Pb(2+)-dependent RNA cleavage. Pb(2+) is required for activity; Mg(2+) inhibits activity, while many divalent and trivalent ions enhance it. The leadzyme structure consists of an RNA duplex interrupted by a trinucleotide bulge. Here, crystal structures determined to 1.8 A resolution, both with Mg(2+) as the sole divalent counterion and with Mg(2+) and Sr(2+) (which mimics Pb(2+) with respect to binding but not catalysis), reveal the metal ion interactions with both the ground state and precatalytic conformations of the leadzyme. Mg(H(2)O)(6)(2+) ions bridge complementary strands of the duplex at multiple locations by binding tandem purines of one RNA strand in the major groove. At one site, Mg(H(2)O)(6)(2+) ligates the phosphodiester backbone of the trinucleotide bulge in the ground state conformation, but not in the precatalytic conformation, suggesting (a) Mg(2+) may inhibit leadzyme activity by stabilizing the ground state and (b) metal ions which displace Mg(2+) from this site may activate the leadzyme. Binding of Sr(2+) to the presumed catalytic Pb(2+) site in the precatalytic leadzyme induces local structural changes in a manner that would facilitate alignment of the catalytic ribose 2'-hydroxyl with the scissile bond for cleavage. These data support a model wherein binding of a catalytic ion to a precatalytic conformation of the leadzyme, in conjunction with the flexibility of the trinucleotide bulge, may facilitate structural rearrangements around the scissle phosphodiester bond favoring configurations that allow bond cleavage.

A lead-dependent DNAzyme with a two-step mechanism

A detailed biochemical and mechanistic study of in vitro selected variants of 8-17 DNAzymes is presented. Even though the 8-17 DNAzyme motif has been obtained through in vitro selection under three different conditions involving 10 mM Mg(2+) (called 8-17), 0.5 mM Mg(2+)/50 mM histidine (called Mg5), or 100 microM Zn(2+) (called 17E), all variants are shown to be the most active with Pb(2+) (8-17: k(obs) approximately 0.5 min(-1); Mg5: k(obs) approximately 2 min(-1); 17E: k(obs) approximately 1 min(-1) with 200 microM Pb(2+) at pH 5.0). For the 17E variant of the 8-17 DNAzyme, the single-turnover rate constants followed the order of Pb(2+) >> Zn(2+) >> Mn(2+) approximately Co(2+) > Ni(2+) > Mg(2+) approximately Ca(2+) > Sr(2+) approximately Ba(2+). The catalytic rate is half-maximal at 13.5 microM Pb(2+), 0.97 mM Zn(2+), or 10.5 mM Mg(2+), suggesting that the metal-binding affinity of the DNAzymes is in the order of Pb(2+) > Zn(2+) > Mg(2+). The Pb(2+)-dependent activity increases linearly with pH and the slope of the plot of log k(obs) versus pH is approximately 1, suggesting a single deprotonation in the rate-limiting step of the reaction. Sequence variations of the DNAzyme confirm the importance of the G*T wobble pair, the two loops and the intervening stem in maintaining the active conformation of the system. While Mg(2+) and Zn(2+) catalyze only a transesterification reaction with formation of a product containing a 2',3'-cyclic phosphate, Pb(2+) catalyzes a transesterification reaction followed by hydrolysis of the 2',3'-cyclic phosphate. Although this two-step mechanism has shown to be operative in protein ribonucleases and in the leadzyme RNAzyme, it is now demonstrated for the first time that this DNAzyme may also use the same mechanism. Therefore, the two-step mechanism is observed in metalloenzymes of all classes, and this 8-17 DNAzyme provides a simple, stable, and cost-effective model system for understanding the structure of Pb(2+)-binding sites and their roles in the two-step mechanism.

The global conformation of the hammerhead ribozyme determined using residual dipolar couplings

The global structure of the hammerhead ribozyme was determined in the absence of Mg(2+) by solution NMR experiments. The hammerhead ribozyme motif forms a branched structure consisting of three helical stems connected to a catalytic core. The (1)H-(15)N and (1)H-(13)C residual dipolar couplings were measured in a set of differentially (15)N/(13)C-labeled ribozymes complexed with an unlabeled noncleavable substrate. The residual dipolar couplings provide orientation information on both the local and the global structure of the molecule. Analysis of the residual dipolar couplings demonstrated that the local structure of the three helical stems in solution is well modeled by an A-form conformation. However, the global structure of the hammerhead in solution in the absence of Mg(2+) is not consistent with the Y-shaped conformation observed in crystal structures of the hammerhead. The residual dipolar couplings for the helical stems were combined with standard NOE and J coupling constant NMR data from the catalytic core. The NOE data show formation of sheared G-A base pairs in domain 2. These NMR data were used to determine the global orientation of the three helical stems in the hammerhead. The hammerhead forms a rather extended structure under these conditions with a large angle between stems I and II ( approximately 153 degrees ), a smaller angle between stems II and III ( approximately 100 degrees ), and the smallest angle between stems I and III ( approximately 77 degrees ). The residual dipolar coupling data also contain information on the dynamics of the molecule and were used here to provide qualitative information on the flexibility of the helical domains in the hammerhead ribozyme-substrate complex.

Parametrization of a force field for metals complexed to biomacromolecules: applications to Fe(II), Cu(II) and Pb(II)

Although techniques for the simulation of biomolecules, such as proteins and RNAs, have greatly advanced in the last decade, modeling complexes of biomolecules with metal ions remains problematic. Precise calculations can be done with quantum mechanical methods but these are prohibitive for systems the size of macromolecules. More qualitative modeling can be done with molecular mechanical potentials but the parametrization of force fields for metals is often difficult, particularly if the bonding between the metal and the groups in its coordination shell has significant covalent character. In this paper we present a method for deriving bond and bond-angle parameters for metal complexes from experimental bond and bond-angle distributions obtained from the Cambridge Structural Database. In conjunction with this method, we also introduce a non-standard energy term of gaussian form that allows us to obtain a stable description of the coordination about a metal center during a simulation. The method was evaluated on Fe(II)-porphyrin complexes, on simple Cu(II) ion complexes and a number of complexes of the Pb(II) ion.

Stabilities and isomeric equilibria in solutions of monomeric metal-ion complexes of guanosine 5'-triphosphate (GTP4-) and inosine 5'-triphosphate (ITP4-) in comparison with those of adenosine 5'-triphosphate (ATP4-)

Under experimental conditions in which the self-association of the purine-nucleoside 5'-triphosphates (PuNTPs) GTP and ITP is negligible, potentiometric pH titrations were carried out to determine the stabilities of the M(H;PuNTP) and M(PuNTP)2-complexes where M2+ = Mg2+, Ca2+, Sr2+. Ba2+, Mn2+, Co2+, Ni2+, Cu2+, Zn2+, or Cd2+ (I = 0.1 M, 25 degrees C). The stabilities of all M(GTP)2- and M(ITP)2- complexes are significantly larger than those of the corresponding complexes formed with pyrimidine-nucleoside 5'-triphosphates (PyNTPs), which had been determined previously under the same conditions. This increased complex stability is attributed, in agreement with previous 1H MNR shift studies, to the formation of macrochelates of the phosphate-coordinated metal ions with N7 of the purine residues. A similar enhanced stability (despite relatively large error limits) was observed for the M(H;PuNTP) complexes, in which H+ is bound to the terminal y-phosphate group, relative to the stability of the M(H;PyNTP)- species. The percentage of the macrochelated isomers in the M(GTP)2- and M(ITP)2- systems was quantified by employing the difference log KMM(PuNTP)-log KMM(PyNTP); the lowest and highest formation degrees of the macrochelates were observed for Mg(ITP)2- and Cu(GTP)2- with 17 +/- 11% and 97 +/- 1%, respectively. From previous studies of M(ATP)2- complexes, it is known that innersphere and outersphere macrochelates may form; that is, in the latter case a water molecule is between N7 and the phosphate-coordinated M2+. Similar conclusions are reached now by comparisons with earlier 1H MNR shift measurements, that is, that Mg(GTP)2- (21 +/- 11%), for example, exists largely in the form of an outersphere macrochelate and Zn(GTP)2- (68 +/- 4%) as an innersphere one. Generally, the overall percentage of macrochelate falls off for a given metal ion in the order M(GTP)2- > M(ITP)2- > M(ATP)2-; this is in accord with the decreasing basicity of N7 and the steric inhibition of the (C6)NH2 group in the adenine residue. Furthermore, although the absolute stability constants of the previously studied M(GMP), M(IMP), and M(AMP) complexes differ by about two to three log units from the present M(PuNTP)2- results, the formation degrees of the macrochelates are astonishingly similar for the two series of nucleotides for a given metal ion and purine-nucleobase residue. The conclusion that N7 of the guanine residue is an especially favored binding site for metal ions is also in accord with observations made for nucleic acids.

Structural characterization of a six-nucleotide RNA hairpin loop found in Escherichia coli, r(UUAAGU)

The binding region of the Escherichia coli S2 ribosomal protein contains a conserved UUAAGU hairpin loop. The structure of the hairpin formed by the oligomer r(GCGU4U5A6A7G8U9CGCA), which has an r(UUAAGU) hairpin loop, was determined by NMR and molecular modeling techniques as part of a study aimed at characterizing the structure and thermodynamics of RNA hairpin loops. Thermodynamic data obtained from melting curves for this RNA oligomer show that it forms a hairpin in solution with the following parameters: DeltaH degrees = -42.8 +/- 2.2 kcal/mol, DeltaS degrees = -127.6 +/- 6.5 eu, and DeltaG degrees (37) = -3.3 +/- 0.2 kcal/mol. Two-dimensional NOESY WATERGATE spectra show an NOE between U imino protons, which suggests that U4 and U9 form a hydrogen bonded U.U pair. The U5(H2') proton shows NOEs to both the A6(H8) proton and the A7(H8) proton, which is consistent with formation of a "U" turn between nucleotides U5 and A6. An NOE between the A7(H2) proton and the U9(H4') proton shows the proximity of the A7 base to the U9 sugar, which is consistent with the structure determined for the six-nucleotide loop. In addition to having a hydrogen-bonded U.U pair as the first mismatch and a U turn, the r(UUAAGU) loop has the G8 base protruding into the solvent. The solution structure of the r(UUAAGU) loop is essentially identical to the structure of an identical loop found in the crystal structure of the 30S ribosomal subunit where the guanine in the loop is involved in tertiary interactions with RNA bases from adjacent regions [Wimberly, B. T., Brodersen, D. E., Clemons, W. M., Morgan-Warren, R. J., Carter, A. P., Vonrhein, C., Hartsch, T., and Ramakrishnan, V. (2000) Nature 407, 327-339]. The similarity of the solution and solid-state structures of this hairpin loop suggests that formation of this hairpin may facilitate folding of 16S RNA.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

Intramolecular secondary structure rearrangement by the kissing interaction of the Neurospora VS ribozyme

Kissing interactions in RNA are formed when bases between two hairpin loops pair. Intra- and intermolecular kissing interactions are important in forming the tertiary or quaternary structure of many RNAs. Self-cleavage of the wild-type Varkud satellite (VS) ribozyme requires a kissing interaction between the hairpin loops of stem-loops I and V. In addition, self-cleavage requires a rearrangement of several base pairs at the base of stem I. We show that the kissing interaction is necessary for the secondary structure rearrangement of wild-type stem-loop I. Surprisingly, isolated stem-loop V in the absence of the rest of the ribozyme is sufficient to rearrange the secondary structure of isolated stem-loop I. In contrast to kissing interactions in other RNAs that are either confined to the loops or culminate in an extended intermolecular duplex, the VS kissing interaction causes changes in intramolecular base pairs within the target stem-loop.

Quantitative analysis of nucleic acid three-dimensional structures

A new computer program to annotate DNA and RNA three-dimensional structures, MC-Annotate, is introduced. The goals of annotation are to efficiently extract and manipulate structural information, to simplify further structural analyses and searches, and to objectively represent structural knowledge. The input of MC-Annotate is a PDB formatted DNA or RNA three-dimensional structure. The output of MC-Annotate is composed of a structural graph that contains the annotations, and a series of HTML documents, one for each nucleotide conformation and base-base interaction present in the input structure. The atomic coordinates of all nucleotides and the homogeneous transformation matrices of all base-base interactions are stored in the structural graph. Symbolic classifications of nucleotide conformations, using sugar puckering modes and nitrogen base orientations around the glycosyl bond, and base-base interactions, using stacking and hydrogen bonding information, are introduced. Peculiarity factors of nucleotide conformations and base-base interactions are defined to indicate their marginalities with all other examples. The peculiarity factors allow us to identify irregular regions and possible stereochemical errors in 3-D structures without interactive visualization. The annotations attached to each nucleotide conformation include its class, its torsion angles, a distribution of the root-mean-square deviations with examples of the same class, the list of examples of the same class, and its peculiarity value. The annotations attached to each base-base interaction include its class, a distribution of distances with examples of the same class, the list of examples of the same class, and its peculiarity value. The distance between two homogeneous transformation matrices is evaluated using a new metric that distinguishes between the rotation and the translation of a transformation matrix in the context of nitrogen bases. MC-Annotate was used to build databases of nucleotide conformations and base-base interactions. It was applied to the ribosomal RNA fragment that binds to protein L11, which annotations revealed peculiar nucleotide conformations and base-base interactions in the regions where the RNA contacts the protein. The question of whether the current database of RNA three-dimensional structures is complete is addressed.

A pH controlled conformational switch in the cleavage site of the VS ribozyme substrate RNA

The VS ribozyme is a 154 nucleotide sequence found in certain natural strains of Neurospora. The RNA can be divided into a substrate and a catalytic domain. Here we present the solution structure of the substrate RNA that is cleaved in a trans reaction by the catalytic domain in the presence of Mg2+. The 30 nucleotide substrate RNA forms a compact helix capped by a flexible loop. The cleavage site bulge contains three non-canonical base-pairs, including an A+.C pair with a protonated adenine. This adenine (A622) is a pH controlled conformational switch that opens up the internal loop at higher pH. The possible significance of this switch for substrate recognition and cleavage is discussed.

Solution structure of an A-tract DNA bend

The solution structure of a DNA dodecamer d(GGCAAAAAACGG)/d(CCGTTTTTTGCC) containing an A-tract has been determined by NMR spectroscopy with residual dipolar couplings. The structure shows an overall helix axis bend of 19 degrees in a geometry consistent with solution and gel electrophoresis experiments. Fourteen degrees of the bending occurs in the GC regions flanking the A-tract. The remaining 5 degrees is spread evenly over its six AT base-pairs. The A-tract is characterized by decreasing minor groove width from the 5' to the 3' direction along the A strand. This is a result of propeller twist in the AT pairs and the increasing negative inclination of the adenine bases at the 3' side of the run of adenine bases. The four central thymine bases all have negative inclination throughout the A-tract with an average value of -6.1 degrees. Although this negative inclination makes the geometry of the A-tract different from all X-ray structures, the proton on N6 of adenine and the O4 of thymine one step down the helix are within distance to form bifurcated hydrogen bonds. The 5' bend of 4 degrees occurs at the junction between the GC flank and the A-tract through a combination of tilt and roll. The larger 3' bend, 10 degrees, occurs in two base steps: the first composed of tilt, -4.1 degrees, and the second a combination of tilt, -4.2 degrees, and roll, 6.0 degrees. This second step is a direct consequence of the change in inclination between an adjacent cytosine base, which has an inclination of -12 degrees, and the next base, a guanine, which has 3 degrees inclination. This bend is a combination of tilt and roll. The large change in inclination allows the formation of a hydrogen bond between the protons of N4 of the 3' cytosine and the O6 of the next 3' base, a guanine, stabilizing the roll component in the bend. These structural features differ from existing models for A-tract bends.For comparison, we also determined the structure of the control sequence, d(GGCAAGAAACGG)/d(CCGTTTCTTGCC), with an AT to GC transition in the center of the A-tract. This structure has no negative inclination in most of the bases within the A-tract, resulting in a bend of only 9 degrees. When ligated in phase, the control sequence has nearly normal mobility in gel electrophoresis experiments.

Ribozyme structures and mechanisms

The past few years have seen exciting advances in understanding the structure and function of catalytic RNA. Crystal structures of several ribozymes have provided detailed insight into the folds of RNA molecules. Models of other biologically important RNAs have been constructed based on structural, phylogenetic, and biochemical data. However, many questions regarding the catalytic mechanisms of ribozymes remain. This review compares the structures and possible catalytic mechanisms of four small self-cleaving RNAs: the hammerhead, hairpin, hepatitis delta virus, and in vitro-selected lead-dependent ribozymes. The organization of these small catalysts is contrasted to that of larger ribozymes, such as the group I intron.

Lead(II)-binding properties of the 5'-monophosphates of adenosine (AMP2-), inosine (IMP2-), and guanosine (GMP2-) in aqueous solution. Evidence for nucleobase-lead(II) interactions

The stability constants of the 1:1 complexes formed between Pb2+ and the nucleosides (Ns), adenosine and guanosine, as well as between the nucleotides (NMP2-), AMP2-, IMP2-, and GMP2-, were determined by potentiometric pH titrations in aqueous solution (25 degrees C; I = 0.1 M, NaNO3). Based on previously established log KPb(R-PO3)Pb versus pKH(R-PO3)H straight-line plots (R-PO3(2-) = simple phosphate monoester or phosphonate ligands where R is a noninteracting site), it is shown that the Pb(IMP) and Pb(GMP) complexes are more stable than is expected on the basis of the basicity of the phosphate group of IMP2- and GMP2-. This means that macrochelates are formed, where the phosphate-coordinated Pb2+ also interacts with N7 of the nucleobase residue. In contrast, the stability of the Pb(AMP) complex is governed by the basicity of the AMP2- phosphate group. These results agree with the observations made for the Pb(Ns)2+ complexes: Pb(adenosine)2+ is very unstable in contrast to Pb(guanosine)2+, the stability of which is very similar to the one of Pb(cytidine)2+ studied previously. The stability constants of the Pb(Ns)2+ complexes also allowed an evaluation of the structure in solution of the monoprotonated Pb(H;NMP)+ complexes, the stabilities of which were also determined. We were able to show that the proton is located at the phosphate group and Pb2+ at the N7/(C6)O site of H(GMP)-; in the case of H(AMP)- Pb2+ is probably about equally distributed between the adenine residue and the monoprotonated phosphate group. On the basis of the stability constants of these complexes and their structures in solution, it is possible to provide a series which reflects the decreasing affinity for Pb2+ of nucleobase residues in single-stranded nucleic acids: guanine approximately equal to cytosine > (hypoxanthine) > adenine > uracil approximately equal to thymine. The Pb2+ affinity of the phosphodiester linkage, -PO3(-)-, is similar to the one of the adenine residue, but is expected to be more significant due to its larger abundance. The relevance of these results for lead-activated ribozymes is briefly discussed.

Multiple folding pathways for the P4-P6 RNA domain

We recently described site-specific pyrene labeling of RNA to monitor Mg(2+)-dependent equilibrium formation of tertiary structure. Here we extend these studies to follow the folding kinetics of the 160-nucleotide P4-P6 domain of the Tetrahymena group I intron RNA, using stopped-flow fluorescence with approximately 1 ms time resolution. Pyrene-labeled P4-P6 was prepared using a new phosphoramidite that allows high-yield automated synthesis of oligoribonucleotides with pyrene incorporated at a specific 2'-amino-2'-deoxyuridine residue. P4-P6 forms its higher-order tertiary structure rapidly, with k(obs) = 15-31 s(-1) (t(1/2) approximately 20-50 ms) at 35 degrees C and [Mg(2+)] approximately 10 mM in Tris-borate (TB) buffer. The folding rate increases strongly with temperature from 4 to 45 degrees C, demonstrating a large activation enthalpy DeltaH(double dagger) approximately 26 kcal/mol; the activation entropy DeltaS(double dagger) is large and positive. In low ionic strength 10 mM sodium cacodylate buffer at 35 degrees C, a slow (t(1/2) approximately 1 s) folding component is also observed. The folding kinetics are both ionic strength- and temperature-dependent; the slow phase vanishes upon increasing [Na(+)] in the cacodylate buffer, and the kinetics switch completely from fast at 30 degrees C to slow at 40 degrees C. Using synchrotron hydroxyl radical footprinting, we confirm that fluorescence monitors the same kinetic events as hydroxyl radical cleavage, and we show that the previously reported slow P4-P6 folding kinetics apply only to low ionic strength conditions. One model to explain the fast and slow folding kinetics postulates that some tertiary interactions are present even without Mg(2+) in the initial state. The fast kinetic phase reflects folding that is facilitated by these interactions, whereas the slow kinetics are observed when these interactions are disrupted at lower ionic strength and higher temperature.