A molecular dynamics simulation study of coaxial stacking in RNA

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5'-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na+ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39 degrees. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.     

Anti-Influenza Virus Activities of Nicked and Circular Dumbbell RNA/DNA Chimeric Oligonucleotides

Abstract

We have designed a new type of antisense oligonucleotide, containing two hairpin loop structures with RNA/DNA base pairs (sense (RNA) and antisense (DNA)) in the double helical stem (nicked and circular dumbbell DNA/RNA chimeric oligonucleotides). The reaction of the nicked and circular dumbbell DNA/RNA chimeric oligonucleotides with RNase H gave the corresponding anti-DNA together with the sense RNA cleavage products. These oligonucleotides were more resistant to exonuclease attack. We also describe the anti-Fluv activities of nicked and circular dumbbell DNMA chimeric oligonucleotides.     

Properties and Anti-HIV Activity of Nicked Dumbbell Oligonucleotides

Abstract

We have designed a new type of oligodeoxyribonucleotide. These oligodeoxyribonucleotides form two hairpin loop structures with base pairs (sense and antisense) in the double helical stem at the 3′ and 5′-ends (nicked dumbbell oligonucleotides). The nicked dumbbell oligonucleotides are molecules with free ends that are more resistant to exonuclease attack. Furthermore, the nicked dumbbell oligonucleotide containing phosphorothioate (P=S) bonds in the hairpin loops has increased nuclease resistance, as compared to the unmodified nicked oligonucleotide. The binding of the nicked dumbbell oligonucleotide to RNA is lower than that of a single-stranded DNA. We also describe the anti-HIV activity of nicked dumbbell oligonucleotides.

     

Packing Features of the all-AT Oligonucleotide d(AAATTT)

Abstract

We describe the packing features of the oligonucleotide duplex d(AAATTT)₂, as determined by X-ray diffraction. There is little information on sequences that only contain A and T bases. The present structure confirms that these sequences tend to pack as a helical arrangement of stacked oligonucleotides in a B conformation with Watson-Crick hydrogen bonding. Our results demonstrate that the virtual TA base step between stacked duplexes has a negative twist that improves base stacking. This observation is consistent with the low stability of TA base steps in B-form DNA.     

Interconversion between Parallel and Antiparallel Conformations of a 4H RNA junction in Domain 3 of Foot-and-Mouth Disease Virus IRES Captured by Dynamics Simulations

RNA junctions are common secondary structural elements present in a wide range of RNA species. They play crucial roles in directing the overall folding of RNA molecules as well as in a variety of biological functions. In particular, there has been great interest in the dynamics of RNA junctions, including conformational pathways of fully base-paired 4-way (4H) RNA junctions. In such constructs, all nucleotides participate in one of the four double-stranded stem regions, with no connecting loops. Dynamical aspects of these 4H RNAs are interesting because frequent interchanges between parallel and antiparallel conformations are thought to occur without binding of other factors. Gel electrophoresis and single-molecule fluorescence resonance energy transfer experiments have suggested two possible pathways: one involves a helical rearrangement via disruption of coaxial stacking, and the other occurs by a rotation between the helical axes of coaxially stacked conformers. Employing molecular dynamics simulations, we explore this conformational variability in a 4H junction derived from domain 3 of the foot-and-mouth disease virus internal ribosome entry site (IRES); this junction contains highly conserved motifs for RNA-RNA and RNA-protein interactions, important for IRES activity. Our simulations capture transitions of the 4H junction between parallel and antiparallel conformations. The interconversion is virtually barrier-free and occurs via a rotation between the axes of coaxially stacked helices with a transient perpendicular intermediate. We characterize this transition, with various interhelical orientations, by pseudodihedral angle and interhelical distance measures. The high flexibility of the junction, as also demonstrated experimentally, is suitable for IRES activity. Because foot-and-mouth disease virus IRES structure depends on long-range interactions involving domain 3, the perpendicular intermediate, which maintains coaxial stacking of helices and thereby consensus primary and secondary structure information, may be beneficial for guiding the overall organization of the RNA system in domain 3.     

A new principle of RNA folding based on pseudoknotting.

Tertiary interactions involving hairpin or interior loops of RNA can lead to extended quasi-continuous double helical stem regions, consisting of coaxially stacked segments of duplex RNA, bridged by single-stranded connections. This type of compact folding plays a role in various strategic regions of RNA molecules. Their role in ribosome functioning, RNA splicing and recognition of tRNA-like structures is discussed.     

RNA2D3D: A program for Generating,Viewing, and Comparing 3-Dimensional Models of RNA

Abstract

Using primary and secondary structure information of an RNA molecule, the program RNA2D3D automatically and rapidly produces a first-order approximation of a 3-dimensional conformation consistent with this information. Applicable to structures of arbitrary branching complexity and pseudoknot content, it features efficient interactive graphical editing for the removal of any overlaps introduced by the initial generating procedure and for making conformational changes favorable to targeted features and subsequent refinement. With emphasis on fast exploration of alternative 3D conformations, one may interactively add or delete base-pairs, adjacent stems can be coaxially stacked or unstacked, single strands can be shaped to accommodate special constraints, and arbitrary subsets can be defined and manipulated as rigid bodies. Compaction, whereby base stacking within stems is optimally extended into connecting single strands, is also available as a means of strategically making the structures more compact and revealing folding motifs. Subsequent refinement of the first-order approximation, of modifications, and for the imposing of tertiary constraints is assisted with standard energy refinement techniques. Previously determined coordinates for any part of the molecule are readily incorporated, and any part of the modeled structure can be output as a PDB or XYZ file. Illustrative applications in the areas of ribozymes, viral kissing loops, viral internal ribosome entry sites, and nanobiology are presented.     

Fast Folding of an RNA Tetraloop on a Rugged Energy Landscape Detected by a Stacking-Sensitive Probe

We investigate the microsecond-timescale kinetics of the RNA hairpin ga^∗cUUCGguc. The fluorescent nucleotide 2-aminopurine (a^∗) reports mainly on base stacking. Ten kinetic traces and the temperature denaturation curve are globally fitted to four-state models of the free-energy surface. In the best-fitting sequential model, the hairpin unfolds over successively larger barriers in at least three stages: stem fraying and increased base-stacking fluctuations; concerted loss of hydrogen bonding and partial unstacking; and additional unstacking of single strands at the highest temperatures. Parallel and trap models also provide adequate fits: such pathways probably also play a role in the complete free-energy surface of the hairpin. To interpret the model states structurally, 200 ns of molecular dynamics, including six temperature-jump simulations, were run. Although the sampling is by no means comprehensive, five different states were identified using hydrogen bonding and base stacking as reaction coordinates. The four to five states required to explain the experiments or simulations set a lower limit on the complexity of this small RNA hairpin's energy landscape.     

Structure of 2′,5′-Linked Tetraribonucleotide Loops: A Novel RNA Motif

Abstract

We report on the three dimensional structure of an RNA hairpin containing a 2′,5′-linked tetraribonucleotide loop, namely, 5′-rGGAC(UUCG)GUCC-3′ (where UUCG = U_2′p5′U_2′p5′C_2′p5′G_2′p5′). We show that the 2′,5′-linked RNA loop adopts a conformation that is quite different from that previously observed for the native 3′,5′-linked RNA loop. The 2′,5′- RNA loop is stabilized by (a) U:G wobble base pairing, with both bases in the anti conformation, (b) extensive base stacking, and (c) sugar–base contacts, all of which contribute to the extra stability of this hairpin structure.     

Dynamic behavior of the telomerase RNA hairpin structure and its relationship to dyskeratosis congenita

In this paper, we present the results from a comprehensive study of nanosecond-scale implicit and explicit solvent molecular dynamics simulations of the wild-type telomerase RNA hairpin. The effects of various mutations on telomerase RNA dynamics are also investigated. Overall, we found that the human telomerase hairpin is a very flexible molecule. In particular, periodically the molecule exhibits dramatic structural fluctuations represented by the opening and closing of a non-canonical base-pair region. These structural deviations correspond to significant disruptions of the direct hydrogen bonding network in the helix, widening of the major groove of the hairpin structure, and causing several U and C nucleotides to protrude into the major groove from the helix permitting them to hydrogen bond with, for example, the P3 domain of the telomerase RNA. We suggest that these structural fluctuations expose a nucleation point for pseudoknot formation. We also found that mutations in the pentaloop and non-canonical region stabilize the hairpin. Moreover, our results show that the hairpin with dyskeratosis congenita mutations is more stable and less flexible than the wild-type hairpin due to base stacking in the pentaloop. The results from our molecular dynamics simulations are in agreement with experimental observations. In addition, they suggest a possible mechanism for pseudoknot formation based on the dynamics of the hairpin structure and also may explain the mutational aspects of dyskeratosis congenita.     

The Impact of Dyskeratosis Congenita Mutations on the Structure and Dynamics of the Human Telomerase RNA Pseudoknot Domain

Abstract

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.     

A COMPARISON OF THE PROPERTIES AND THE SOLUTION STRUCTURE FOR RNA AND DNA QUADRUPLEXES WHICH CONTAIN TWO GGAGG SEQUENCES JOINED WITH A TETRANUCLEOTIDE LINKER

Abstract

We have determined solution structure of r(GGAGGUUUUGGAGG) (R14) by NMR; the RNA 14-mer forms an intra-strand parallel quadruplex with a G-tetrad and a hexad, in which a G-tetrad core is augmented by association of two A residues. The quadruplex further forms a dimer through stacking interaction between the hexads. In order to obtain insight into the difference between RNA and DNA quadruplexes, we synthesized the corresponding DNA 14-mer, d(GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. K⁺ ions increased the thermal stability of both R14 and D14 structures. The binding affinity of K⁺ ions to R14 was much higher than that to D14. The CD and gel electrophoretic studies suggest that D14 forms a quadruplex entirely different from that of R14 in the presence of K⁺ ions; two molecules of D14 form a quadruplex with both antiparallel and parallel strand alignments and with diagonal loops at both ends of the stacked G-tetrads. The NMR study also gave results that are consistent with such structure: alternate glycosidic conformation, 5′G(syn)-G(anti)3′, and characteristic chemical shift data observed for many quadruplexes containing diagonal TTTT loops.     

The Influence of Nearest Neighbours on the Efficiency of Coaxial Stacking at Contiguous Stacking Hybridization of Oligodeoxyribonucleotides

Contiguous stacking hybridization of oligodeoxyribonucleotides with a stem of preformed minihairpin structure of a DNA template was studied with the use of UV‐melting technique. It was shown that the free‐energy of the coaxial stacking interaction (ΔG°_ST at 37°C, 1 M NaCl, pH 7.4) at the complementary interface XA^*pTY/ZATV (an asterisk stands for a nick) strongly depends on the type of nearest neighbor bases X and Y flanking the nicked dinucleotide step. The maximum efficiency of the coaxial stacking was observed for the PuA^*pTPy/PuATPy interface, whereas the minimum efficiency was obtained for the PyA^*pTPu/PyATPu interface. A 5′‐phosphate residue in the nick enhances the coaxial stacking. In dependence on duplex structure the observed efficiency of A^*T/AT coaxial stacking varied from (? 0.97 kcal/mol) for unphosphorylated TA^*TA/TATA interface to three‐fold higher value (? 2.78 kcal/mol) for GA^*pTT/AATC interface.     

An examination of coaxial stacking of helical stems in a pseudoknot motif: the gene 32 messenger RNA pseudoknot of bacteriophage T2

The RNA pseudoknot located at the 5' end of the gene 32 messenger RNA of bacteriophage T2 contains two A-form helical stems connected by two loops, in an H-type pseudoknot topology. A combination of multidimensional NMR methods and isotope labeling were used to investigate the pseudoknot structure, resulting in a more detailed structural model than provided by earlier homonuclear NMR studies. Of particular significance, the interface between the stacked helical stems within the pseudoknot motif is described in detail. The two stems are stacked in a coaxial manner, with an approximately 18 degrees rotation of stem1 relative to stem2 about an axis that is parallel to the helical axis. This rotation serves to relieve what would otherwise be a relatively close phosphate-phosphate contact at the junction of the two stems, while preserving the stabilizing effects of base stacking. The ability of the NMR data to determine pseudoknot bending was critically assessed. The data were found to be a modestly precise indicator of pseudoknot bending, with the angle between the helical axes of stem1 and stem2 being in the range of 15+/-15 degrees. Pseudoknot models with bend angles within this range are equally consistent with the data, since they differ by only small amounts in the relatively short-range interproton distances from which the structure was derived. The gene 32 messenger RNA pseudoknot was compared with other RNA structures with coaxial or near-coaxial stacked helical stems.     

Index to Authors,Volume 3

Abstract

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

The Conformation of a Conserved Stem-Loop Structure in Ribosomal RNA

Abstract

The RNA of small ribosomal subunits contains a conserved stem-loop structure near the 3′ end. Characteristics for the hairpins are: (a) a nine-basepairs stem; (b) a conserved ^A-U _U-Gjunction in the stem; (c) a conserved sequence Gm⁶ ₂Am⁶ ₂A sequence in the loop (except yeast mitochondria and mutants from bacteria). We are using UV-optics, micro-calorimetry and 500 MHz-NMR to investigate fragments of about 50 nucleotides cleaved from the 3′ ends of small ribosomal subunit RNA's by bacteriocins. Our preliminary conclusions are: (1) Dimethylation of the adenines in the loop destabilizes the hairpin because of an increased stacking; (2) melting of the hairpin starts at the ends as well as in the middle at the 5-H junction; (3) basepair substitutions have an unexpectedly large effect on thermal stability.     

Long promoter sequences form higher-order G-quadruplexes: an integrative structural biology study of c-Myc,k-Ras and c-Kit promoter sequences

We report on higher-order G-quadruplex structures adopted by long promoter sequences obtained by an iterative integrated structural biology approach. Our approach uses quantitative biophysical tools (analytical ultracentrifugation, small-angle X-ray scattering, and circular dichroism spectroscopy) combined with modeling and molecular dynamics simulations, to derive self-consistent structural models. The formal resolution of our approach is 18 angstroms, but in some cases structural features of only a few nucleotides can be discerned. We report here five structures of long (34–70 nt) wild-type sequences selected from three cancer-related promoters: c-Myc, c-Kit and k-Ras. Each sequence studied has a unique structure. Three sequences form structures with two contiguous, stacked, G-quadruplex units. One longer sequence from c-Myc forms a structure with three contiguous stacked quadruplexes. A longer c-Kit sequence forms a quadruplex-hairpin structure. Each structure exhibits interfacial regions between stacked quadruplexes or novel loop geometries that are possible druggable targets. We also report methodological advances in our integrated structural biology approach, which now includes quantitative CD for counting stacked G-tetrads, DNaseI cleavage for hairpin detection and SAXS model refinement. Our results suggest that higher-order quadruplex assemblies may be a common feature within the genome, rather than simple single quadruplex structures.     

The Twist,Writhe, and Linking Number Distributions in Closed Circular DNA

Abstract

We discuss the predictions which follow from the assumption of statistically independent twist and writhe distributions of given variances in circular DNA with single-strand nicks. The nature of the topoisomer distribution produced upon covalent closure of the nicks is described, as well as the nature of the twist and writhe distributions in the fully-closed molecules. In particular, we show how the distributions depend on the magnitudes of the given variances, and how the relative magnitudes of the variances can be deduced from experiment. One additional consequence of the theory is the prediction of a necessary difference between the temperature coefficient of the twist in nicked versus fully-closed circular DNA. The ratio of the two twist coefficients turns out to depend only on the ratio of the twist and writhe variances in nicked DNA.