Structures of synthetic polynucleotides in the A-RNA and A'-RNA conformations: x-ray diffraction analyses of the molecular conformations of polyadenylic acid--polyuridylic acid and polyinosinic acid--polycytidylic acid

Accurately Modeling RNA Stem-Loops in an Implicit Solvent Environment

Ribonucleic acid (RNA) molecules can adopt a variety of secondary and tertiary structures in solution, with stem-loops being one of the more common motifs. Here, we present a systematic analysis of 15 RNA stem-loop sequences simulated with molecular dynamics simulations in an implicit solvent environment. Analysis of RNA cluster ensembles showed that the stem-loop structures can generally adopt the A-form RNA in the stem region. Loop structures are more sensitive, and experimental structures could only be reproduced with modification of CH···O interactions in the force field, combined with an implicit solvent nonpolar correction to better model base stacking interactions. Accurately modeling RNA with current atomistic physics-based models remains challenging, but the RNA systems studied herein may provide a useful benchmark set for testing other RNA modeling methods in the future.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

Dual Time-Scale Proton Transfer and High-Energy, Long-Lived Excitons Unveiled by Broadband Ultrafast Time-Resolved Fluorescence in Adenine-Uracil RNA Duplexes

In contrast to the immense amount of research on electronically excited DNA, surprisingly little has been done about the excited states of RNA. Herein, we demonstrate an ultrafast broadband time-resolved fluorescence and fluorescence anisotropy study to probe directly the intrinsic fluorescence and overall dynamics of the fluorescence from a homopolymeric adenine·uracil RNA duplex adopting the A-form structure. The results unveiled complex deactivation through distinctive multichannels mediated by states of varied energy, a character of charge transfer, and a lifetime from sub-picosecond to nanoseconds. In particular, we observed an unprecedented kinetic isotopic effect and participation of unusual proton transfer from states in two discrete energies and time domains. We also identified a high-energy nanosecond emission that we attributed to its fluorescence anisotropy to long-lived weakly emissive excitons not reported in DNA. These distinguishing features originate from the stacking, pairing, and local hydration environment specific to the A-form conformation of the adenine·uracil double helix.

Molecular structure of a U•A-U-rich RNA triple helix with 11 consecutive base triples

Three-dimensional structures have been solved for several naturally occurring RNA triple helices, although all are limited to six or fewer consecutive base triples, hindering accurate estimation of global and local structural parameters. We present an X-ray crystal structure of a right-handed, U•A-U-rich RNA triple helix with 11 continuous base triples. Due to helical unwinding, the RNA triple helix spans an average of 12 base triples per turn. The double helix portion of the RNA triple helix is more similar to both the helical and base step structural parameters of A′-RNA rather than A-RNA. Its most striking features are its wide and deep major groove, a smaller inclination angle and all three strands favoring a C3′-endo sugar pucker. Despite the presence of a third strand, the diameter of an RNA triple helix remains nearly identical to those of DNA and RNA double helices. Contrary to our previous modeling predictions, this structure demonstrates that an RNA triple helix is not limited in length to six consecutive base triples and that longer RNA triple helices may exist in nature. Our structure provides a starting point to establish structural parameters of the so-called ‘ideal’ RNA triple helix, analogous to A-RNA and B-DNA double helices.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

Polarizable force field for RNA based on the classical drude oscillator

RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.

Coarse-Grained Double-Stranded RNA Model from Quantum-Mechanical Calculations

A coarse-grained model for simulating structural properties of double-stranded RNA is developed with parameters obtained from quantum-mechanical calculations. This model follows previous parametrization for double-stranded DNA, which is based on mapping the all-atom picture to a coarse-grained four-bead scheme. Chemical and structural differences between RNA and DNA have been taken into account for the model development. The parametrization is based on simulations using density functional theory (DFT) on separate units of the RNA molecule without implementing experimental data. The total energy is decomposed into four terms of physical significance: hydrogen bonding interaction, stacking interactions, backbone interactions, and electrostatic interactions. The first three interactions are treated within DFT, whereas the last one is included within a mean field approximation. Our double-stranded RNA coarse-grained model predicts stable helical structures for RNA. Other characteristics, such as structural or mechanical properties are reproduced with a very good accuracy. The development of the coarse-grained model for RNA allows extending the spatial and temporal length scales accessed by computer simulations and being able to model RNA-related biophysical processes, as well as novel RNA nanostructures.

Insights into the structural features and stability of peptide nucleic acid with a D-prolyl-2-aminocyclopentane carboxylic acid backbone that binds to DNA and RNA

Peptide nucleic acid (PNA) is a powerful biomolecule with a wide variety of important applications. In this work, the molecular structures and binding affinity of PNA with a D-prolyl-2-aminocyclopentane carboxylic acid backbone (acpcPNA) that binds to both DNA and RNA were studied using molecular dynamics simulations. The simulated structures of acpcPNA-DNA and acpcPNA-RNA duplexes more closely resembled the typical structures of B-DNA and A-RNA than the corresponding duplexes of aegPNA. The calculated binding free energies are in good agreement with the experimental results that the acpcPNA-DNA duplex is more stable than the acpcPNA-RNA duplex regardless of the base sequences. The results provide further insights in the relationship between structure and stability of this unique PNA system.

Triplet-Based Codon Organization Optimizes the Impact of Synonymous Mutation on Nucleic Acid Molecular Dynamics

Since the elucidation of the genetic code almost 50 years ago, many nonrandom aspects of its codon organization remain only partly resolved. Here, we investigate the recent hypothesis of ‘dual-use’ codons which proposes that in addition to allowing adjustment of codon optimization to tRNA abundance, the degeneracy in the triplet-based genetic code also multiplexes information regarding DNA’s helical shape and protein-binding dynamics while avoiding interference with other protein-level characteristics determined by amino acid properties. How such structural optimization of the code within eukaryotic chromatin could have arisen from an RNA world is a mystery, but would imply some preadaptation in an RNA context. We analyzed synonymous (protein-silent) and nonsynonymous (protein-altering) mutational impacts on molecular dynamics in 13823 identically degenerate alternative codon reorganizations, defined by codon transitions in 7680 GPU-accelerated molecular dynamic simulations of implicitly and explicitly solvated double-stranded aRNA and bDNA structures. When compared to all possible alternative codon assignments, the standard genetic code minimized the impact of synonymous mutations on the random atomic fluctuations and correlations of carbon backbone vector trajectories while facilitating the specific movements that contribute to DNA polymer flexibility. This trend was notably stronger in the context of RNA supporting the idea that dual-use codon optimization and informational multiplexing in DNA resulted from the preadaptation of the RNA duplex to resist changes to thermostability. The nonrandom and divergent molecular dynamics of synonymous mutations also imply that the triplet-based code may have resulted from adaptive functional expansion enabling a primordial doublet code to multiplex gene regulatory information via the shape and charge of the minor groove.

Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics

The backbone dihedral parameters of the Amber RNA force field were improved by fitting using multiple linear regression to potential energies determined by quantum chemistry calculations. Five backbone and four glycosidic dihedral parameters were fit simultaneously to reproduce the potential energies determined by a high-level density functional theory calculation (B97D3 functional with the AUG-CC-PVTZ basis set). Umbrella sampling was used to determine conformational free energies along the dihedral angles, and these better agree with the population of conformations observed in the protein data bank for the new parameters than for the conventional parameters. Molecular dynamics simulations performed on a set of hairpin loops, duplexes and tetramers with the new parameter set show improved modeling for the structures of tetramers CCCC, CAAU, and GACC, and an RNA internal loop of noncanonical pairs, as compared to the conventional parameters. For the tetramers, the new parameters largely avoid the incorrect intercalated structures that dominate the conformational samples from the conventional parameters. For the internal loop, the major conformation solved by NMR is stable with the new parameters, but not with the conventional parameters. The new force field performs similarly to the conventional parameters for the UUCG and GCAA hairpin loops and the [U(UA)6A]2 duplex.

High-resolution measurement of long-range distances in RNA: pulse EPR spectroscopy with TEMPO-labeled nucleotides

Distance measurements in RNAs by pulse EPR with TEMPO-labeled nucleotides allow for model free conversion of distances into base-pair separation.

Structural information at atomic resolution of biomolecular assemblies, such as RNA and RNA protein complexes, is fundamental to comprehend biological function. Modern spectroscopic methods offer exceptional opportunities in this direction. Here we present the capability of pulse EPR to report high-resolution long-range distances in RNAs by means of a recently developed spin labeled nucleotide, which carries the TEMPO group directly attached to the nucleobase and preserves Watson–Crick base-pairing. In a representative RNA duplex with spin-label separations up to 28 base pairs (≈8 nm) we demonstrate that the label allows for a model-free conversion of inter-spin distances into base-pair separation (Δbp) if broad-band pulse excitation at Q band frequencies (34 GHz) is applied. The observed distance distribution increases from ±0.2 nm for Δbp = 10 to only ±0.5 nm for Δbp = 28, consistent with only small deviations from the “ideal” A-form RNA structure. Molecular dynamics (MD) simulations conducted at 20 °C show restricted conformational freedom of the label. MD-generated structural deviations from an “ideal” A-RNA geometry help disentangle the contributions of local flexibility of the label and its neighboring nucleobases and global deformations of the RNA double helix to the experimental distance distributions. The study demonstrates that our simple but strategic spin labeling procedure can access detailed structural information on RNAs at atomic resolution over distances that match the size of macromolecular RNA complexes.

Computational Approaches to Nucleic Acid Origami

Recent advances in experimental DNA origami have dramatically expanded the horizon of DNA nanotechnology. Complex 3D suprastructures have been designed and developed using DNA origami with applications in biomaterial science, nanomedicine, nanorobotics, and molecular computation. Ribonucleic acid (RNA) origami has recently been realized as a new approach. Similar to DNA, RNA molecules can be designed to form complex 3D structures through complementary base pairings. RNA origami structures are, however, more compact and more thermodynamically stable due to RNA's non-canonical base pairing and tertiary interactions. With all these advantages, the development of RNA origami lags behind DNA origami by a large gap. Furthermore, although computational methods have proven to be effective in designing DNA and RNA origami structures and in their evaluation, advances in computational nucleic acid origami is even more limited. In this paper, we review major milestones in experimental and computational DNA and RNA origami and present current challenges in these fields. We believe collaboration between experimental nanotechnologists and computer scientists are critical for advancing these new research paradigms.

Sequence-controlled RNA self-processing: computational design, biochemical analysis, and visualization by AFM

Reversible chemistry allowing for assembly and disassembly of molecular entities is important for biological self-organization. Thus, ribozymes that support both cleavage and formation of phosphodiester bonds may have contributed to the emergence of functional diversity and increasing complexity of regulatory RNAs in early life. We have previously engineered a variant of the hairpin ribozyme that shows how ribozymes may have circularized or extended their own length by forming concatemers. Using the Vienna RNA package, we now optimized this hairpin ribozyme variant and selected four different RNA sequences that were expected to circularize more efficiently or form longer concatemers upon transcription. (Two-dimensional) PAGE analysis confirms that (i) all four selected ribozymes are catalytically active and (ii) high yields of cyclic species are obtained. AFM imaging in combination with RNA structure prediction enabled us to calculate the distributions of monomers and self-concatenated dimers and trimers. Our results show that computationally optimized molecules do form reasonable amounts of trimers, which has not been observed for the original system so far, and we demonstrate that the combination of theoretical prediction, biochemical and physical analysis is a promising approach toward accurate prediction of ribozyme behavior and design of ribozymes with predefined functions.

A single-stranded architecture for cotranscriptional folding of RNA nanostructures

Artificial DNA and RNA structures have been used as scaffolds for a variety of nanoscale devices. In comparison to DNA structures, RNA structures have been limited in size, but they also have advantages: RNA can fold during transcription and thus can be genetically encoded and expressed in cells. We introduce an architecture for designing artificial RNA structures that fold from a single strand, in which arrays of antiparallel RNA helices are precisely organized by RNA tertiary motifs and a new type of crossover pattern. We constructed RNA tiles that assemble into hexagonal lattices and demonstrated that lattices can be made by annealing and/or cotranscriptional folding. Tiles can be scaled up to 660 nucleotides in length, reaching a size comparable to that of large natural ribozymes.

NAFlex: a web server for the study of nucleic acid flexibility

We present NAFlex, a new web tool to study the flexibility of nucleic acids, either isolated or bound to other molecules. The server allows the user to incorporate structures from protein data banks, completing gaps and removing structural inconsistencies. It is also possible to define canonical (average or sequence-adapted) nucleic acid structures using a variety of predefined internal libraries, as well to create specific nucleic acid conformations from the sequence. The server offers a variety of methods to explore nucleic acid flexibility, such as a colorless wormlike-chain model, a base-pair resolution mesoscopic model and atomistic molecular dynamics simulations with a wide variety of protocols and force fields. The trajectories obtained by simulations, or imported externally, can be visualized and analyzed using a large number of tools, including standard Cartesian analysis, essential dynamics, helical analysis, local and global stiffness, energy decomposition, principal components and in silico NMR spectra. The server is accessible free of charge from the mmb.irbbarcelona.org/NAFlex webpage.

Hydrotalcite Intercalated siRNA: Computational Characterization of the Interlayer Environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.

The crystal structure of an oligo(U):pre-mRNA duplex from a trypanosome RNA editing substrate

Guide RNAs bind antiparallel to their target pre-mRNAs to form editing substrates in reaction cycles that insert or delete uridylates (Us) in most mitochondrial transcripts of trypanosomes. The 5' end of each guide RNA has an anchor sequence that binds to the pre-mRNA by base-pair complementarity. The template sequence in the middle of the guide RNA directs the editing reactions. The 3' ends of most guide RNAs have ∼15 contiguous Us that bind to the purine-rich unedited pre-mRNA upstream of the editing site. The resulting U-helix is rich in G·U wobble base pairs. To gain insights into the structure of the U-helix, we crystallized 8 bp of the U-helix in one editing substrate for the A6 mRNA of Trypanosoma brucei. The fragment provides three samples of the 5'-AGA-3'/5'-UUU-3' base-pair triple. The fusion of two identical U-helices head-to-head promoted crystallization. We obtained X-ray diffraction data with a resolution limit of 1.37 Å. The U-helix had low and high twist angles before and after each G·U wobble base pair; this variation was partly due to shearing of the wobble base pairs as revealed in comparisons with a crystal structure of a 16-nt RNA with all Watson-Crick base pairs. Both crystal structures had wider major grooves at the junction between the poly(U) and polypurine tracts. This junction mimics the junction between the template helix and the U-helix in RNA-editing substrates and may be a site of major groove invasion by RNA editing proteins.

Structure of the HIV-1 frameshift site RNA bound to a small molecule inhibitor of viral replication

Programmed -1 translational frameshifting is an essential event in the replication cycle of HIV. Frameshifting is required for expression of the viral Pol proteins, and drug-like molecules that target this process may inhibit HIV replication. A small molecule stimulator of HIV-1 frameshifting and inhibitor of viral replication, DB213 (RG501), was previously discovered from a high-throughput screen. However, the mechanistic basis for this compound's effects was unknown, and to date no structural information exists for small molecule effectors of frameshifting. Here, we investigate the binding of DB213 to the frameshift site RNA and have determined the structure of this complex by NMR. Binding of DB213 stabilizes the RNA and increases its melting temperature by 10 °C. The ligand binds to a primary site on the RNA stem-loop, although nonspecific interactions are also detected. The compound binds in the major groove and spans a distance of 9 base pairs. DB213 hydrogen bonds to phosphate groups on opposite sides of the major groove and alters the conformation of a conserved GGA bulge in the RNA. This study may provide a starting point for structure-based optimization of compounds targeting the HIV-1 frameshift site RNA.

RNA dimerization monitored by fluorescence correlation spectroscopy

Fluorescence correlation spectroscopy (FCS) provides a versatile tool to investigate molecular interaction under native conditions, approximating infinite dilution. One precondition for its application is a sufficient difference between the molecular weights of the fluorescence-labelled unbound and bound ligand. In previous studies, an 8-fold difference in molecular weights or correspondingly a 1.6-fold difference in diffusion coefficients was required to accurately distinguish between two diffusion species by FCS. In the presented work, the hybridization of two complementary equally sized RNA single strands was investigated at an excellent signal-to-noise ratio enabled by the highly photostable fluorophore Atto647N. The fractions of ssRNA and dsRNA were quantified by applying multicomponent model analysis of single autocorrelation functions and globally fitting several autocorrelation functions. By introducing a priori knowledge into the fitting procedure, 1.3- to 1.4-fold differences in diffusion coefficients of single- and double-stranded RNA of 26, 41, and 54 nucleotides could be accurately resolved. Global fits of autocorrelation functions of all titration steps enabled a highly accurate quantification of diffusion species fractions and mobilities. At a high signal-to-noise ratio, the median of individually fitted autocorrelation functions allowed a robust representation of heterogeneous data. These findings point out the possibility of studying molecular interaction of equally sized molecules based on their diffusional behavior, which significantly broadens the application spectrum of FCS.

Predicting translational diffusion of evolutionary conserved RNA structures by the nucleotide number

Ribonucleic acids are highly conserved essential parts of cellular life. RNA function is determined to a large extent by its hydrodynamic behaviour. The presented study proposes a strategy to predict the hydrodynamic behaviour of RNA single strands on the basis of the polymer size. By atom-level shell-modelling of high-resolution structures, hydrodynamic radius and diffusion coefficient of evolutionary conserved RNA single strands (ssRNA) were calculated. The diffusion coefficients D of 17–174 nucleotides (nt) containing ssRNA depended on the number of nucleotides N with D = 4.56 × 10⁻¹⁰ N^−0.39 m² s⁻¹. The hydrodynamic radius R_H depended on N with R_H = 5.00 × 10⁻¹⁰N^0.38 m. An average ratio of the radius of gyration and the hydrodynamic radius of 0.98 ± 0.08 was calculated in solution. The empirical law was tested by in solution measured hydrodynamic radii and radii of gyration and was found to be highly consistent with experimental data of evolutionary conserved ssRNA. Furthermore, the hydrodynamic behaviour of several evolutionary unevolved ribonucleic acids could be predicted. Based on atom-level shell-modelling of high-resolution structures and experimental hydrodynamic data, empirical models are proposed, which enable to predict the translational diffusion coefficient and molecular size of short RNA single strands solely on the basis of the polymer size.

Flexibility of Short-Strand RNA in Aqueous Solution as Revealed by Molecular Dynamics Simulation: Are A-RNA and A´-RNA Distinct Conformational Structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A′-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA – namely major groove width, inclination and the number of base pairs in a helical twist – over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A′-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors).

Crystallization and X-ray diffraction analysis of the Trp/amber editing site of hepatitis delta virus (+)RNA: a case of rational design

RNA editing by mammalian ADAR1 (Adenosine Deaminase Acting on RNA) is required for the life cycle of the hepatitis delta virus (HDV). Editing extends the single viral open reading frame to yield two protein products of alternate length. ADARs are believed to recognize double-stranded RNA substrates via a ;structure-based' readout mechanism. Crystals of 10-mer duplexes representing the HDV RNA-editing site diffracted to 1.35 A resolution, but suffered from merohedral twinning and averaging of the base registry. Expansion of the construct to include two flanking 3 x 1 internal loops yielded crystals in the primitive tetragonal space group P4(1)2(1)2 or P4(3)2(1)2. X-ray diffraction data were collected to 2.8 A resolution, revealing a unit cell with parameters a = 62.5, c = 63.5 A. The crystallization and X-ray analysis of multiple forms of the HDV RNA-editing substrate, encounters with common RNA crystal-growth defects and a strategy to overcome these problems are reported.

Single-molecule measurements of the persistence length of double-stranded RNA

Over the past few years, it has become increasingly apparent that double-stranded RNA (dsRNA) plays a far greater role in the life cycle of a cell than previously expected. Numerous proteins, including helicases, polymerases, and nucleases interact specifically with the double helix of dsRNA. To understand the detailed nature of these dsRNA-protein interactions, the (bio)chemical, electrostatic, and mechanical properties of dsRNA need to be fully characterized. We present measurements of the persistence length of dsRNA using two different single-molecule techniques: magnetic tweezers and atomic force microscopy. We deduce a mean persistence length for long dsRNA molecules of 63.8 +/- 0.7 nm from force-extension measurements with the magnetic tweezers. We present atomic force microscopy images of dsRNA and demonstrate a new method for analyzing these, which yields an independent, yet consistent value of 62 +/- 2 nm for the persistence length. The introduction of these single-molecule techniques for dsRNA analysis opens the way for real-time, quantitative analysis of dsRNA-protein interactions.

Atomic model of the Thermus thermophilus 70S ribosome developed in silico

The ribosome is a large molecular complex that consists of at least three ribonucleic acid molecules and a large number of proteins. It translates genetic information from messenger ribonucleic acid and makes protein accordingly. To better understand ribosomal function and provide information for designing biochemical experiments require knowledge of the complete structure of the ribosome. For expanding the structural information of the ribosome, we took on the challenge of developing a detailed Thermus thermophilus ribosomal structure computationally. By combining information derived from the low-resolution x-ray structure of the 70S ribosome (providing the overall fold), high-resolution structures of the ribosomal subunits (providing the local structure), sequences, and secondary structures, we have developed an atomic model of the T. thermophilus ribosome using a homology modeling approach. Our model is stereochemically sound with a consistent single-species sequence. The overall folds of the three ribosomal ribonucleic acids in our model are consistent with those in the low-resolution crystal structure (root mean-square differences are all <1.9 angstroms). The large overall interface area (approximately 2500 angstroms2) of intersubunit bridges B2a, B3, and B5, and the inherent flexibility in regions connecting the contact residues are consistent with these bridges serving as anchoring patches for the ratcheting and rolling motions between the two subunits during translocation.

Fluorescence energy transfer between fluorescein label and DNA intercalators to detect nucleic acids hybridization in homogeneous media

A general approach to detecting nucleic acid sequences in homogeneous media by means of steady-state fluorescence measurements is proposed. The methodology combines the use of a fluorescence-labeled single-strand DNA model probe, the complementary single-strand DNA target, and a DNA intercalator. The probe was fluorescein labeled to a spacer arm at the N4 position of the cytosine amino groups in polyribocytidylic acid (5'), poly(C), which acts as a model DNA probe. The complementary strand was polyriboinosinic acid (5'), poly(I), as a model of the target, and the energy transfer acceptor was an intercalator, either ethidium bromide or ethidium homodimer. In previous papers we have shown that the fluorescence intensity of the fluorescein label decreases when labeled poly(C) hybridizes with poly(I), and this fluorescence quenching can be used to detect DNA hybridization or renaturation in homogeneous media. In this paper we demonstrate that fluorescence resonance energy transfer (FRET) between fluorescein labeled to poly(C) and an intercalator agent takes place when single-stranded poly(C) hybridizes with poly(I), and we show how the fluorescence energy transfer further decreases the steady-state fluorescence intensity of the label, thus increasing the detection limit of the method. The main aim of this work was to develop a truly homogeneous detection system for specific nucleic acid hybridization in solution using steady-state fluorescence and FRET, but with the advantage of only having to label the probe with the energy donor since the energy acceptor is intercalated spontaneously. Moreover, the site label is not critical and can be labeled randomly in the DNA strand. Thus, the method is simpler than those published previously based on FRET. The experiments were carried out in both direct and competitive formats.

Structural basis for the unusual properties of 2',5' nucleic acids and their complexes with RNA and DNA

To provide insights into the unusual properties of 2',5' nucleic acids (iso nucleic acids), that includes their rejection by Nature as information molecules, modeling studies have been carried out to examine if they indeed possess the stereochemical ability to form helical duplexes and triplexes, just as their 3',5' linked constitutional isomers. The results show that the formation of helical duplexes with 2',5' linkages demands a mandatory displacement of the Watson and Crick base pairs from the helical axis, as a direct consequence of the lateral shift of the sugar-phosphate backbone from the periphery towards the interior of the helix. Thus, both duplexes and triplexes formed with a 2',5'-sugar-phosphate backbone possess this intrinsic trait, manifested normally only in A type duplexes of DNA and RNA. It was found that only a 10-fold symmetric parallel triplex with isomorphous T.AT triplets is stereochemically favorable for isoDNA with 'extended' nucleotide repeats, unlike the 12-fold symmetric triplex favored by DNA. The wider nature of a 12-fold triplex, concomitant with mandatory slide requirement for helix formation in isoDNA, demands even larger displacement, especially with 'extended' nucleotide structural repeats, thereby violating symmetry. However, a symmetric triplex possessing higher twist, can be naturally formed for isoDNA with a 'compact' nucleotide repeat. Two nanosecond molecular dynamics simulation of a 2',5'-B DNA duplex, formed with an intrinsic base pair displacement of -3.3 A, does not seem to favor a total transition to a typical A type duplex, although enhanced slide, X-displacement, decrease in helical rise and narrowing of the major groove during simulation seem to indicate a trend. Modeling of the interaction between the chimeric isoDNA.RNA duplex and E. coli RNase H has provided a structural basis for the inhibitory action of the enzyme. Interaction of residues Gln 80, Trp 81, Asn 16 and Lys 99, of E. coli RNase H with DNA of the DNA.RNA hybrid, are lost when the DNA backbone is replaced by isoDNA. Based on modeling and experimental observations, it is argued that 2',5' nucleic acids possess restricted conformational flexibility for helical polymorphism. The inability of isoDNA to favor the biologically relevant B form duplex and the associated topological inadequacies related to nucleic acid compaction and interactions with regulatory proteins may be some of the factors that might have led to the rejection of 2',5' links.

Visualization of unwinding activity of duplex RNA by DbpA, a DEAD box helicase, at single-molecule resolution by atomic force microscopy

The Escherichia coli protein DbpA is unique in its subclass of DEAD box RNA helicases, because it possesses ATPase-specific activity toward the peptidyl transferase center in 23S rRNA. Although its remarkable ATPase activity had been well defined toward various substrates, its RNA helicase activity remained to be characterized. Herein, we show by using biochemical assays and atomic force microscopy that DbpA exhibits ATP-stimulated unwinding activity of RNA duplex regardless of its primary sequence. This work presents an attempt to investigate the action of DEAD box proteins by a single-molecule visualization methodology. Our atomic force microscopy images enabled us to observe directly the unwinding reaction of a DEAD box helicase on long stretches of double-stranded RNA. Specifically, we could differentiate between the binding of DbpA to RNA in the absence of ATP and the formation of a Y-shaped intermediate after its progression through double-stranded RNA in the presence of ATP. Recent studies have questioned the designation of DbpA, in particular, and DEAD box proteins in general as RNA helicases. However, accumulated evidence and the results reported herein suggest that these proteins are indeed helicases that resemble in many aspects the DNA helicases.

Crystal structures of A-DNA duplexes

All crystal structures of A-DNA duplexes exhibit a typical crystal packing, with the termini of one molecule abutting the shallow grooves of symmetry related neighbors, while all other forms (B, Z, and RNA) tend to form infinitely stacked helices. The A-DNA arrangement leads to the formation of shallow groove base multiples that have implications for the structure of DNA in compacted states. The characteristic packing leaves big solvent channels, which can be sometimes occupied by B-DNA duplexes. Comparisons of the structures of the same oligomer crystallizing in two different space groups and of different sequences crystallizing in the same space group show that the lattice forces dominate the A-DNA conformation in the crystals, complicating the effort to elucidate the influence of the base sequence on the structures. Nevertheless, in both alternating and nonalternating fragments some sequence effects can still be uncovered. Furthermore, several studies have started to define the minimal sequence changes or chemical modifications that can interconvert the oligomers between different double-helical conformers (A-, B-, and Z-form). Overall, it is seen that the rigid nucleotide principle applies to the oligomeric fragments. Besides the structures of the naked DNAs, their interactions with water, polyamines, and metal ions have attracted considerable attention. There are conserved patterns in the hydration, involving both the grooves and the backbone, which are different from those of B-DNA or Z-DNA. Overall, A-DNA seems to be more economically hydrated than B-DNA, particularly around the sugar-phosphate backbone. Spermine was found to be able to bind exclusively to either of the grooves or to the phosphate groups of the backbone, or exhibit a mixed binding mode. The located metal cations prefer binding to guanine bases and phosphate groups. The only mispairs investigated in A-DNA are the wobble pairs, yielding structural insight into their effects on helix stabilities and hydration. G.T wobble pairs have been determined in various sequence contexts, where they differentially affect the conformations and stableness of the duplexes. The structure of a G.m5C base pair, which surprisingly also adopted the wobble conformation, suggests that a similar geometry may transiently exist for G.C pairs. These results from the crystalline state will be compared to the solution state and discussed in relation to their relevance in biology.

Specificity and mechanism analysis of hepatitis C virus RNA-dependent RNA polymerase

The RNA-dependent RNA polymerase encoded by the hepatitis C virus (HCV) NS5B gene has been expressed as a nonfusion protein in bacterial cells and purified to homogeneity using sequential chromatographic columns. The purified NS5B protein exhibited RNA-dependent RNA polymerase activity using poly(A) template and the K(m) and V(max) were determined as 8.4 microM and 1976 pmol/mg-min, respectively. This full-length NS5B protein exhibited much stronger binding affinity toward the 30-mer poly(G) than other homopolymeric RNAs of the same size. For the first time, we demonstrate that the HCV NS5B was able to bind various ribonucleotides. Using a panel of oligonucleotides varying in length, we studied the NS5B catalytic efficiency and proposed the size of the NS5B active site to be 8-10 nucleotides. The multifunctional nature of NS5B protein is also discussed and compared with other viral RNA polymerases.

Aminoglycoside binding in the major groove of duplex RNA: the thermodynamic and electrostatic forces that govern recognition

We use a combination of spectroscopic, calorimetric, viscometric and computer modeling techniques to characterize the binding of the aminoglycoside antibiotic, tobramycin, to the polymeric RNA duplex, poly(rI).poly(rC), which exhibits the characteristic A-type conformation that is conserved among natural and synthetic double-helical RNA sequences. Our results reveal the following significant features: (i) CD-detected binding of tobramycin to poly(rI).poly(rC) reveals an apparent site size of four base-pairs per bound drug molecule; (ii) tobramycin binding enhances the thermal stability of the host poly(rI).poly(rC) duplex, the extent of which decreases upon increasing in Na(+) concentration and/or pH conditions; (iii) the enthalpy of tobramycin- poly(rI).poly(rC) complexation increases with increasing pH conditions, an observation consistent with binding-induced protonation of one or more drug amino groups; (iv) the affinity of tobramycin for poly(rI).poly(rC) is sensitive to both pH and Na(+) concentration, with increases in pH and/or Na(+) concentration resulting in a concomitant reduction in binding affinity. The salt dependence of the tobramycin binding affinity reveals that the drug binds to the host RNA duplex as trication. (v) The thermodynamic driving force for tobramycin- poly(rI).poly(rC) complexation depends on pH conditions. Specifically, at pH< or =6.0, tobramycin binding is entropy driven, but is enthalpy driven at pH > 6.0. (vi) Viscometric data reveal non-intercalative binding properties when tobramycin complexes with poly(rI).poly(rC), consistent with a major groove-directed mode of binding. These data also are consistent with a binding-induced reduction in the apparent molecular length of the host RNA duplex. (vii) Computer modeling studies reveal a tobramycin-poly(rI). poly(rC) complex in which the drug fits snugly at the base of the RNA major groove and is stabilized, at least in part, by an array of hydrogen bonding interactions with both base and backbone atoms of the host RNA. These studies also demonstrate an inability of tobramycin to form a stable low-energy complex with the minor groove of the poly(rI).poly(rC) duplex. In the aggregate, our results suggest that tobramycin-RNA recognition is dictated and controlled by a broad range of factors that include electrostatic interactions, hydrogen bonding interactions, drug protonation reactions, and binding-induced alterations in the structure of the host RNA. These modulatory effects on tobramycin-RNA complexation are discussed in terms of their potential importance for the selective recognition of specific RNA structural motifs, such as asymmetric internal loops or hairpin loop-stem junctions, by aminoglycoside antibiotics and their derivatives.

Global structure and flexibility of hairpin ribozymes with extended terminal helices

Global structure and flexibility of three different hairpin ribozyme constructs have been analyzed by measuring their electric dichroism decay in various buffers at temperatures between 2 and 30 degrees C. The hairpin ribozyme is characterized by two independently folding domains A and B that are connected through a hinge and have to interact to enable catalysis. The analyzed constructs feature extended terminal helices 1 and 4 with 27 and 25 bp, respectively, to increase the sensitivity of the molecular rotational diffusion time constants with respect to the interdomain bending angle. Constructs HP1 and HP2 cannot cleave because of a G+1A change at the 3'-side of the cleavage site; in HP1 the helices 2 and 3 that flank the hinge form a continuous double helical segment; in HP2 and HP3, a six nucleotide bulge confers flexibility to the expected bending site; HP3 is a cleavable form of HP2 with a G+1-base. For comparison, a standard RNA double helix with 72 bp was included in our analysis. The dichroism decay curves of the hairpin constructs after pulses of low electric field strengths can be fitted to single exponentials taus, whereas the curves after pulses of high field strengths require two exponentials. In all cases, time constants increase with RNA concentration, indicating intermolecular interactions. Extrapolation of the tausvalues measured in standard buffer (50 mM Tris (pH 7.5) and 12 mM MgCl2) to zero RNA concentration provide values of 112, 93, and 73 ns for HP1, HP2 and HP3, respectively, at 30 degrees C, indicating increasingly compact structures. The 72 bp RNA reference under corresponding conditions did not show a dependence of its decay time constant on the RNA concentration nor on the field strength; its time constant is 175 ns (standard buffer, 30 degrees C). The observation of two relaxation processes for the hairpin constructs at high field strengths indicates stretching to a more elongated state; the fast process with a time constant of the order of 50 ns is assigned to reversion of stretching, the slow process to overall rotation. The overall rotational time of the stretched state at 20 degrees C is close to that for a completely stretched rigid state; at 30 degrees C the experimental values are around 70 % of that expected for a completely stretched rigid state, indicating flexibility and/or residual bending. Bead models were constructed to simulate dichroism decay curves. The time constants observed for the 72 bp RNA are as expected for a rigid rod with a rise of 2.8 A per base-pair. Based on this rise per base-pair for models of a V and a Y-shape, we estimate average bending angles of 80(+/-20) degrees and 105 (+/-25) degrees, respectively, for the catalytically active hairpin ribozyme HP3. The energy required for stretching is of the order of the thermal energy.

The role of RNA pseudoknot stem 1 length in the promotion of efficient -1 ribosomal frameshifting

The ribosomal frameshifting signal present in the genomic RNA of the coronavirus infectious bronchitis virus (IBV) contains a classic hairpin-type RNA pseudoknot that is believed to possess coaxially stacked stems of 11 bp (stem 1) and 6 bp (stem 2). We investigated the influence of stem 1 length on the frameshift process by measuring the frameshift efficiency in vitro of a series of IBV-based pseudoknots whose stem 1 length was varied from 4 to 13 bp in single base-pair increments. Efficient frameshifting depended upon the presence of a minimum of 11 bp; pseudoknots with a shorter stem 1 were either non-functional or had reduced frameshift efficiency, despite the fact that a number of them had a stem 1 with a predicted stability equal to or greater than that of the wild-type IBV pseudoknot. An upper limit for stem 1 length was not determined, but pseudoknots containing a 12 or 13 bp stem 1 were fully functional. Structure probing analysis was carried out on RNAs containing either a ten or 11 bp stem 1; these experiments confirmed that both RNAs formed pseudoknots and appeared to be indistinguishable in conformation. Thus the difference in frameshifting efficiency seen with the two structures was not simply due to an inability of the 10 bp stem 1 construct to fold into a pseudoknot. In an attempt to identify other parameters which could account for the poor functionality of the shorter stem 1-containing pseudoknots, we investigated, in the context of the 10 bp stem 1 construct, the influence on frameshifting of altering the slippery sequence-pseudoknot spacing distance, loop 2 length, and the number of G residues at the bottom of the 5'-arm of stem 1. For each parameter, it was possible to find a condition where a modest stimulation of frameshifting was observable (about twofold, from seven to a maximal 17 %), but we were unable to find a situation where frameshifting approached the levels seen with 11 bp stem 1 constructs (48-57 %). Furthermore, in the next smaller construct (9 bp stem 1), changing the bottom four base-pairs to G.C (the optimal base composition) only stimulated frameshifting from 3 to 6 %, an efficiency about tenfold lower than seen with the 11 bp construct. Thus stem 1 length is a major factor in determining the functionality of this class of pseudoknot and this has implications for models of the frameshift process.

Structure of a 16-mer RNA duplex r(GCAGACUUAAAUCUGC)2 with wobble C.A+ mismatches

The crystal structure of a 16-mer, the longest known RNA duplex, has been determined at 2.5 A resolution. The hexadecamer r(GCAGACUUAAAUCUGC) contains isolated C.A/A.C mismatches with two hydrogen bonds. The two hydrogen bonds in the mismatches suggests that N1 of A is protonated even though the crystallization was done at neutral pH. Therefore, the C.A mismatch is a C.A+ wobble similar to the G.U wobble. The two C.A+ pairs are isolated by four Watson-Crick pairs and flanked by five Watson-Crick base-pairs on either sides. Kinks/bends of 20 degrees are observed at the wobble sites. The Watson-Crick base-pair A5.U26 on the 5'-side of the first C6.A27(+) wobble has a twist angle of 27 degrees compared to the 3'-side U7.A28 pair of 36 degrees. The twist angles are reversed (37 degrees and 26 degrees) in the second A11(+).C22 wobble because of the approximate dyad in the molecule, the flanking base-pair sequences are A.U pairs. The wobbles expand the major groove to 7.1 A/7.3 A. The duplexes form helical columns and are tightly packed around the 31-screw axis. The minor grooves of adjacent columns in juxtaposition interact through the O2' atoms and the anionic phosphate oxygen atoms.

Nucleic acid vibrational circular dichroism, absorption, and linear dichroism spectra. II. A DeVoe theory approach

The DeVoe polarizability theory is used to calculate vibrational circular dichroism (VCD) and infrared (IR) absorption spectra of four polyribonucleotides: poly(rA) x poly(rU), poly(rU) x poly(rA) x poly(rU), poly(rG) x poly(rC), and poly(rC+) x poly(rI) x poly(rC). This is the first report on the use of the DeVoe theory to calculate VCD, oriented VCD, IR absorption, and IR linear dichroism (LD) spectra of double- and triple-stranded polyribonucleotides. Results are reported for DeVoe theory calculations--within the base-stretching 1750-1550 cm(-1) spectral region--on several proposed multistranded polyribonucleotide geometries. The calculated spectra obtained from these proposed geometries are compared with previously reported measured and calculated VCD and IR spectral results. Base-base hydrogen-bonding effects on the frequencies and magnitudes of the base carbonyl stretching modes are explicitly considered. The good agreements found between calculated and measured spectra are proposed to be further evidence of the usefulness of the DeVoe theory in drawing three-dimensional structural conclusions from measured polyribonucleotide VCD and IR spectra.

Interactions of the antiviral quinoxaline derivative 9-OH-B220 [2, 3-dimethyl-6-(dimethylaminoethyl)- 9-hydroxy-6H-indolo-[2, 3-b]quinoxaline] with duplex and triplex forms of synthetic DNA and RNA

The binding of an antiviral quinoxaline derivative, 2,3-dimethyl- 6 - (dimethylaminoethyl) - 9 - hydroxy - 6H - indolo - [2,3 - b]quinoxaline (9-OH-B220), to synthetic double and triple helical DNA (poly(dA).poly(dT) and poly(dA).2poly(dT)) and RNA (poly(rA). poly(rU) and poly (rA).2poly(rU)) has been characterized using flow linear dichroism (LD), circular dichroism (CD), fluorescence spectroscopy, and thermal denaturation. When either of the DNA structures or the RNA duplex serve as host polymers a strongly negative LD is displayed, consistent with intercalation of the chromophoric ring system between the base-pairs/triplets of the nucleic acid structures. Evidence for this geometry also includes weak induced CD signals and strong increments of the fluorescence emission intensities upon binding of the drug to each of these polymer structures. In agreement with intercalative binding, 9-OH-B220 is found to effectively enhance the thermal stability of both the double and triple helical states of DNA as well as the RNA duplex. In the case of poly(dA).2poly(dT), the drug provides an unusually large stabilization of its triple helical state; upon binding of 9-OH-B220 the triplex-to-duplex equilibrium is shifted towards higher temperature by 52.5 deg. C in a 10 mM sodium cacodylate buffer (pH 7.0) containing 100 mM NaCl and 1 mM EDTA. When triplex RNA serves as host structure, LD indicates that the average orientation angle between the drug chromophore plane and the helix axis of the triple helical RNA is only about 60 to 65 degrees. Moreover, the thermal stabilizing capability, as well as the fluorescence increment, CD inducing power and perturbations of the absorption envelope, of 9-OH-B220 in complex with the RNA triplex are all less pronounced than those observed for the complexes with DNA and duplex RNA. These features indicate binding of 9-OH-B220 in the wide and shallow minor groove of poly(rA).2poly(rU). Based on the present results, some implications for the applications of this low-toxic, antiviral and easily administered drug in an antigene strategy, as well as its potential use as an antiretroviral agent, are discussed.

Site size of cooperative single-stranded RNA binding by poliovirus RNA-dependent RNA polymerase

The poliovirus RNA-dependent RNA polymerase binds cooperatively to single-stranded RNA. We have determined the minimal RNA-binding site size of the poliovirus polymerase using binding titration with oligonucleotides of increasing length. A dramatic increase in affinity was observed when the length of the oligo(U) increased from 8 to 10 nucleotides (nt), arguing that the minimal size of RNA for polymerase binding is 10 nt. Another increase in affinity seen as the oligo(U) reached 24 nt suggests that a 24-nucleotide RNA can be occupied by two polymerase molecules. Direct binding of wild-type polymerase to oligo(U)12 and oligo(U)24 RNAs showed differences in affinity and cooperativity consistent with this model. The increase in binding affinity seen for oligo(U)10 suggests either that the RNA-binding determinants are widely spaced on the polymerase structure or that a substantial conformational change in the polymerase occurs upon the filling of its RNA-binding site.

Solution structure of RNA duplexes containing alternating CG base pairs: NMR study of r(CGCGCG)2 and 2'-O-Me(CGCGCG)2 under low salt conditions

Structures of r(CGCGCG)2 and 2'-O-Me(CGCGCG)2 have been determined by NMR spectroscopy under low salt conditions. All protons and phosphorus nuclei resonances have been assigned. Signals of H5'/5" have been assigned stereospecifically. All 3JH,H and 3JP,H coupling constants have been measured. The structures were determined and refined using an iterative relaxation matrix procedure (IRMA) and the restrained MD simulation. Both duplexes form half-turn, right-handed helices with several conformational features which deviate significantly from a canonical A-RNA structure. Duplexes are characterised as having C3'-endo sugar pucker, very low base-pair rise and high helical twist and inclination angles. Helices are overwound with <10 bp per turn. There is limited inter-strand guanine stacking for CG steps. Within CG steps of both duplexes, the planes of the inter-strand cytosines are not parallel while guanines are almost parallel. For the GC steps this pattern is reversed. The 2'-O-methyl groups are spatially close to the 5'-hydrogens of neighbouring residues from the 3'-side and are directed towards the minor groove of 2'-O-Me(CGCGCG)2 forming a hydrophobic layer. Solution structures of both duplexes are similar; the effect of 2'-O-methylation on the parent RNA structure is small. This suggests that intrinsic properties imposed by alternating CG base pairs govern the overall conformation of both duplexes.

Crystal structure of 2'-O-Me(CGCGCG)2, an RNA duplex at 1.30 A resolution. Hydration pattern of 2'-O-methylated RNA

The molecular and crystal structure of 2'-O-Me (CGCGCG)2 has been determined using synchrotron radiation at near-atomic resolution (1.30 A), the highest resolution to date in the RNA field. The crystal structure is a half-turn A-type helix with some helical parameters deviating from canonical A-RNA, such as low base pair rise, elevated helical twist and inclination angles. In CG steps, inter-strand guanines are parallel while cytosines are not parallel. In steps GC this motif is reversed. This type of regularity is not seen in other RNA crystal structures. The structure includes 44 water molecules and two hydrated Mg2+ions one of which lies exactly on the crystallographic 2-fold axis. There are distinct patterns of hydration in the major and the minor grooves. The major groove is stabilised by water clusters consisting of fused five- and six-membered rings. Minor groove contains only a single row of water molecules; each water bridges either two self-parallel cytosines or two self-parallel guanines by a pair of hydrogen bonds. The structure provides the first view of the hydration scheme of 2'-O-methylated RNA duplex.

Photochemical and photophysical studies of 3-amino-6-iodoacridine and the inactivation of lambda phage

The photochemistry and photophysics of 3-amino-6-iodoacridine (Acr-I) was studied. Photolysis (350 nm) of Acr-I (free base) generates products consistent with a free radical intermediate in methanol, benzene and carbon tetrachloride. The Acr-I hydrochloride is shown to bind to calf thymus DNA and to the self-complementary dinucleotide cytidylyl-(3'-5')-guanosine (CpG) miniduplex in a manner similar to that of proflavine (Acr-NH2), a known DNA intercalator. The Acr-I is shown to more efficiently nick supercoiled plasmid DNA pBR322 upon 350 nm or 420 nm photolysis than Acr-NH2. The efficiency of Acr-I-sensitized DNA nicking is not oxygen dependent. Photolysis of the Acr-I/(CpG)2 complex leads to cleavage of the dinucleotide and to cytidine base release by selective damage to a specific ribose moiety. Dinucleotide cleavage occurs equally well in the presence or absence of oxygen, thereby eliminating a singlet oxygen- or peroxyl radical-mediated process. Photolysis of Acr-I in the presence of a mononucleotide (GMP) or a non-self-complementary dinucleotide (uridylyl-[3'-5']-cytidine-UpC) does not lead to fragmentation and base release. Similarly, photolysis of the Acr-NH2/(CpG)2 complex does not lead to fragmentation and base release. The data indicate that photolysis of an iodinated intercalator bound to CpG or plasmid DNA generates an intercalated aryl radical and that the reactive intermediate initiates a sequence of reactions that efficiently nick nucleic acids. The inactivation of lambda phage sensitized by Acr-I with UV (350 nm) light is oxygen independent but with visible (420 nm) light is strongly oxygen dependent. The Acr-I fluoresces more intensely when excited at 446 than at 376 nm. Thus, UV photolysis may lead to C-I bond homolysis and free radical formation, a process that is not energetically feasible with visible light. The results demonstrate the difficulty of extrapolating model studies involving simple molecules and DNA to understanding the mechanism of viral inactivation with a particular sensitizer.

Understanding the thermodynamic stability of an RNA hairpin and its mutant

The thermodynamic stability of RNA hairpin loops has been a subject of considerable interest in the recent past (Wimberly et al., 1991). There have been experimental reports indicating that the hairpins with a C(UUCG)G loop sequence are thermodynamically very stable (Wimberly et al., 1991). We used the solution structure of GGAC(UUCG)GUCC (Cheong et al., 1990; Varani et al., 1991) as the starting conformation in our attempt to understand its thermodynamic stability. We carried out molecular dynamics/free energy simulations to understand the basis for the destabilization of the C(UUCG)G loop by mutating cytosine (C7)-->uracil. Because of the limited length of simulation and the presence of kinetic barriers (solvent intervention) to the uracil-->cytosine mutation, all of our computed free energy differences are based on multiple forward simulations. Based on these calculations we find that the cytosine-->uracil mutation in the loop destabilizes it by approximately 1.5kcal/mol relative to that of the reference state, an A-form RNA but with cytosine (C7) looped out. This is the same sign and magnitude as that observed in the thermodynamic studies carried out by Varani et al.(1991). We have carried out free energy component analysis to understand the effect of mutating the cytosine residue to uracil on the thermodynamic stability of the C(UUCG)G hairpin loops. Our calculations show that the most significant contribution to the stability is from the phosphate group linking U5 and U6, which favors the cytosine residue over uracil by about 6.0 kcal/mol. The residues U5, U6, and G8 in the loop region also contribute significantly to the stability. The contributions from the salt and solvent compensate each other, indicating the dynamic nature of interactions of the environment with the nucleic acid system and the coupling between these two components.

Prospects for the therapeutic use of antigene oligonucleotides

An outgrowth of classic nucleic acid interaction studies, oligonucleotide-directed triple helix formation is a unique method for creating highly specific chemical ligands that recognize and bind to particular sequences of duplex DNA. Under permissive conditions, these oligonucleotide-based compounds can approach or exceed the binding affinity and sequence specificity of natural DNA-binding proteins. Triple helix recognition has been found to be useful in certain cell-free applications including precise chromosome fragmentation. It has been proposed that such oligonucleotides could also form the basis for gene-targeted (antigene) drugs that might repress transcription from undesired genes in living cells. However, current strategies for oligonucleotide-directed triple helix formation suffer from important constraints involving requirements for stabilizing binding conditions, restrictions on permitted target sequences, and inefficient nuclear delivery of oligonucleotides. Implementation of oligonucleotide-directed triple helix formation as a viable approach to cancer therapy must therefore await clever solutions to a series of fascinating problems.

Berenil [1,3-bis(4'-amidinophenyl)triazene] binding to DNA duplexes and to a RNA duplex: evidence for both intercalative and minor groove binding properties

Berenil is an antitrypanosomal agent that binds to nucleic acid duplexes. The generally accepted mode of berenil binding is via complexation into the minor groove of AT-rich domains of DNA double helices. We find that berenil can bind to RNA as well as DNA duplexes, while exhibiting properties characteristic of both intercalation as well as minor groove binding. More specifically, we use spectroscopic, calorimetric, and hydrodynamic techniques to characterize berenil binding to four DNA duplexes and to one RNA duplex. Our results reveal the following features: (i) Berenil binding to the poly[d(A-T)]2, poly(dA).poly(dT), poly[d(I-C)]2, poly[d(G-C)]2, and poly(rA).poly(rU) duplexes exhibits intercalative as well as minor groove binding characteristics. (ii) The apparent "site sizes" associated with berenil binding to these five duplexes range from 1 to 13 base pairs per bound berenil and depend, in part, on the host duplex. One of the site sizes common to all five duplexes is consistent with berenil binding to the minor groove. (iii) The apparent berenil binding affinity follows the hierarchy: poly(dA).poly(dT) > poly-[d(A-T)]2 approximately poly[d(I-C)]2 >> poly(rA).poly(rU) > poly[d(G-C)]2. (iv) Viscometric data reveal properties characteristic of a significant contribution from an intercalative mode of binding when berenil interacts with the poly[d(A-T)]2, poly[d(I-C)]2, poly[d(G-C)]2, and poly(rA).poly(rU) duplexes, while revealing an apparent nonintercalative mode when the drug binds to the poly(dA).poly(dT) duplex. (v) Berenil binding unwinds negative supercoils in the pBR322 plasmid, an observation consistent with an intercalative mode of binding to duplex DNA. (vi) Salt-dependent melting data suggest that both positively charged amidino groups of berenil participate in the complexation of the drug to the poly[d(I-C)]2, poly[d(A-T)]2, poly(dA).poly(dT), and poly(rA).poly(rU) duplexes, while also suggesting that the binding event is site-specific. In the aggregate, our results suggest that, in contrast to the conventional wisdom, berenil can exhibit intercalative as well as minor groove binding properties when it binds to both DNA and RNA duplexes, while also exhibiting a preference for DNA duplexes with unobstructed minor grooves. We comment on the potential correlation between drugs, such as berenil, that exhibit "mixed" binding motifs and those that express anticancer activity via inhibition of topoisomerase I activity.