Naturally occurring cross-links in yeast chromosomal DNA.

Chromosome-size yeast DNA molecules with a number average molecular weight (Mn) of 3-4 X 10(8) were isolated from sucrose gradients after sedimentation of lysed yeast spheroplasts. Resedimentation showed that the molecules were isolated without introducing appreciable single-strand or double-strand breaks. The presence of cross-links in these molecules was suggested by the observation that the apparent Mn in alkali was greater than expected for separated single strands. Since cross-linked molecules would have strands which fail to separate upon denaturation, this was tested more directly. Neutralization of alkaline denaturing conditions resulted in up to 70% of the intact molecules rapidly reforming duplex structures, as shown by equilibrium banding in CsCI. Experiments with larger E. coli DNA molecules (Mn = 5.2 X 10(8)) indicated that the conditions used were sufficient to denature completely molecules of this size. Results of enzyme treatments suggest that the cross-links are not RNA or protein. Experiments with density-labeled yeast DNA molecules showed that the rapid reformation of duplex DNA is not the consequence either of a bimolecular reaction between separated DNA strands or of intrastrand renaturation. The data indicate that when the yeast DNA molecules are completely denatured, the strands fail to separate. Hence they must be cross-linked. Experiments with sheared DNA show that there are small number of cross-links, one to four, permolecule.     

Phosphates, DNA, and the search for nonterrean life: a second generation model for genetic molecules

Phosphate groups are found and used widely in biological chemistry. We have asked whether phosphate groups are likely to be important to the functioning of genetic molecules, including DNA and RNA. From observations made on synthetic analogs of DNA and RNA where the phosphates are replaced by nonanionic linking groups, we infer a set of rules that highlight the importance of the phosphodiester backbone for the proper functioning of DNA as a genetic molecule. The polyanionic backbone appears to give DNA the capability of replication following simple rules, and evolving. The polyanionic nature of the backbone appears to be critical to prevent the single strands from folding, permitting them to act as templates, guiding the interaction between two strands to form a duplex in a way that permits simple rules to guide the molecular recognition event, and buffering the sensitivity of its physicochemical properties to changes in sequence. We argue that the feature of a polyelectrolyte (polyanion or polycation) may be required for a "self-sustaining chemical system capable of Darwinian evolution." The polyelectrolyte structure therefore may be a universal signature of life, regardless of its genesis, and unique to living forms as well.     

RNA facilitates RecA-mediated DNA pairing and strand transfer between molecules bearing limited regions of homology

The RecA protein ofEscherichia coli catalyzes homologous pairing and strand exchange between a wide range of molecules showing nucleotide sequence complementarity, including a linear duplex and a single-stranded DNA molecule. We demonstrate that RecA can promote formation of joint molecules when the duplex contains an RNA/DNA hairpin and a single-stranded circle serves as the pairing partner. A chimeric RNA/DNA hairpin can be used to form stable joint molecules with as little as 15 bases of shared homology as long as the RNA stretch contains complementarity to the circle. The joint molecule bears some resemblance to a triple helical structure composed of RNA residues surrounded by two DNA strands which are in a parallel orientation. Evidence is presented that supports the notion that short stretches of RNA can be used in homologous pairing reactions at lengths below that required for DNA-DNA heteroduplex formation.     

A potent 2'-O-methylated RNA-based microRNA inhibitor with unique secondary structures

MicroRNAs (miRNAs) are involved in various biological processes and human diseases. The development of strong low-molecular weight inhibitors of specific miRNAs is thus expected to be useful in providing tools for basic research or in generating promising new therapeutic drugs. We have previously described the development of 'Tough Decoy (TuD) RNA' molecules, which achieve the long-term suppression of specific miRNA activity in mammalian cells when expressed from a lentivirus vector. In our current study, we describe new synthetic miRNA inhibitors, designated as S-TuD (Synthetic TuD), which are composed of two fully 2'-O-methylated RNA strands. Each of these strands includes a miRNA-binding site. Following the hybridization of paired strands, the resultant S-TuD forms a secondary structure with two stems, which resembles the corresponding TuD RNA molecule. By analyzing the effects of S-TuD against miR-21, miR-200c, miR-16 and miR-106b, we have elucidated the critical design features of S-TuD molecules that will provide optimum inhibitory effects following transfection into human cell lines. We further show that the inhibitory effects of a single transfection of S-TuD-miR200c are quite long-lasting (>7 days) and induce partial EMT, the full establishment of which requires 11 days when using a lentivirus vector that expresses TuD-miR200c continuously.     

A 7.1 kb linear DNA molecule of Theileria parva has scrambled rDNA sequences and open reading frames for mitochondrially encoded proteins.

Theileria parva, an intralymphocytic protozoan parasite of cattle, contains a linear 7.1 kb DNA element with terminal inverted repeat sequences. The molecule is transcribed into low molecular weight RNA, and both DNA strands encode short stretches of unique sequences, usually < 100 nucleotides, which are similar to large (LSU) or small (SSU) ribosomal subunit RNA. Phylogenetically conserved conformational rRNA domains were assembled from the discontinuous rDNA sequences using comparative secondary structure modelling. For example, a minimum of four predicted sequences, two derived from each DNA strand, is required to assemble domain V of LSU rRNA which participates in peptidyl transferase activity. The discontinuities in the identified rRNA domains fall within regions of no known functional significance. Hence, it is likely that the element encodes fragmented rDNA genes and the mature rRNA is unconventional, consisting of several fragments of RNA, primarily held together by intermolecular and intramolecular base pairing. The element also has ORFs for components of the last two mitochondrial electron transport enzyme complexes. The structure of the parasite DNA element, its protein coding capacity and scrambled rDNA gene sequences, are reminiscent of the mitochondrial genome of Chlamydomonas reinhardtii. We propose that the 7.1 kb element is equivalent to the mitochondrial DNA of T. parva, although a number of its features are unusual for this family of extrachromosomal DNA molecules.     

The effects of deoxyribonucleic acid secondary structure on tertiary structure.

The secondary structure of supercoiled DNA was varied by changes in ionic strength. For I = 0.075-0.4 the structure remained in the previously established branched form with only minor alterations in molecular dimensions. In 4M-NaCl, which induces linear DNA to change its secondary structure to the C structure and brings about an increase in the superhelix density of the molecule, no extra branches were observed on the molecules. The limiting factors that dictate supercoil structure seem to be the number and position of potential branch points and the proximity with which the two intertwining DNA strands can approach each other on the arms of the branches. This value is close to 10nm under the conditions described, and is 14-15nm at I = 0.2. It is suggested that such values should be borne in mind when models of chromosome structure are being constructed.     

MINT: software to identify motifs and short-range interactions in trajectories of nucleic acids

Structural biology experiments and structure prediction tools have provided many high-resolution three-dimensional structures of nucleic acids. Also, molecular dynamics force field parameters have been adapted to simulating charged and flexible nucleic acid structures on microsecond time scales. Therefore, we can generate the dynamics of DNA or RNA molecules, but we still lack adequate tools for the analysis of the resulting huge amounts of data. We present MINT (Motif Identifier for Nucleic acids Trajectory) — an automatic tool for analyzing three-dimensional structures of RNA and DNA, and their full-atom molecular dynamics trajectories or other conformation sets (e.g. X-ray or nuclear magnetic resonance-derived structures). For each RNA or DNA conformation MINT determines the hydrogen bonding network resolving the base pairing patterns, identifies secondary structure motifs (helices, junctions, loops, etc.) and pseudoknots. MINT also estimates the energy of stacking and phosphate anion-base interactions. For many conformations, as in a molecular dynamics trajectory, MINT provides averages of the above structural and energetic features and their evolution. We show MINT functionality based on all-atom explicit solvent molecular dynamics trajectory of the 30S ribosomal subunit.     

Secondary structure prediction of interacting RNA molecules

Computational tools for prediction of the secondary structure of two or more interacting nucleic acid molecules are useful for understanding mechanisms for ribozyme function, determining the affinity of an oligonucleotide primer to its target, and designing good antisense oligonucleotides, novel ribozymes, DNA code words, or nanostructures. Here, we introduce new algorithms for prediction of the minimum free energy pseudoknot-free secondary structure of two or more nucleic acid molecules, and for prediction of alternative low-energy (sub-optimal) secondary structures for two nucleic acid molecules. We provide a comprehensive analysis of our predictions against secondary structures of interacting RNA molecules drawn from the literature. Analysis of our tools on 17 sequences of up to 200 nucleotides that do not form pseudoknots shows that they have 79% accuracy, on average, for the minimum free energy predictions. When the best of 100 sub-optimal foldings is taken, the average accuracy increases to 91%. The accuracy decreases as the sequences increase in length and as the number of pseudoknots and tertiary interactions increases. Our algorithms extend the free energy minimization algorithm of Zuker and Stiegler for secondary structure prediction, and the sub-optimal folding algorithm by Wuchty et al. Implementations of our algorithms are freely available in the package MultiRNAFold.     

A large inverted repeat sequence overlaps two acceptor splice sites in adenovirus

The distribution of nucleotide sequences resembling functional sites for mRNA splicing was examined by computer-directed searches in order to determine what factors may influence splice site selection in nuclear precursors. In particular, the distribution of large potentially stable hairpin structures or regions of extensive dyad symmetry was studied in adenovirus sequences. One region, spanning 106 nucleotides, was found at 66.4 map units, overlapping back-to-back acceptor sites for two mRNA molecules, those coding for the 100K protein and the 72K DNA binding protein, which are transcribed from opposite strands. This region displays exceptional dyad symmetry and is potentially capable of forming a single, highly stable hairpin when transcribed. It seems likely that the secondary structure as well as the primary structure of RNA plays a role in determining the correct splicing of these mRNA molecules.     

Secondary structure maps of ribosomal RNA and DNA: I. Processing of Xenopus laevis ribosomal RNA and structure of single-stranded ribosomal DNA

Precursor and mature ribosomal RNA molecules from Xenopus laevis were examined by electron microscopy. A reproducible arrangement of hairpin loops was observed in these molecules. Maps based on this secondary structure were used to determine the arrangement of sequences in precursor RNA molecules and to identify the position of mature rRNAs within the precursors. A processing scheme was derived in which the 40 S rRNA is cleaved to 38 S RNA, which then yields 34 S plus 18 S RNA. The 34 S RNA is processed to 30 S, and finally to 28 S rRNA. The pathway is analogous to that of L-cell rRNA but differs from HeLa rRNA in that no 20 S rRNA intermediate was found. X. laevis 40 S rRNA (M_r = 2.7 × 10⁶) is much smaller than HeLa or L-cell 45 8 rRNA (M_r = 4.7 × 10⁶), but the arrangement of mature rRNA sequences in all precursors is very similar. Experiments with ascites cell 3′-exonuclease show that the 28 S region is located at or close to the 5′-end of the 40 S rRNA.Secondary structure maps were obtained also for single-stranded molecules of ribosomal DNA. The region in the DNA coding for the 40 S rRNA could be identified by its regular structure, which closely resembles that of the RNA. Regions corresponding to the 40 S RNA gene alternate with non-transcribed spacer regions along strands of rDNA. The latter have a large amount of irregular secondary structure and vary in length between different repeating units. A detailed map of the rDNA repeating unit was derived from these experiments.Optical melting studies are presented, showing that rRNAs with a high (G + C) content exhibit significant hypochromicity in the formamide/urea-containing solution that was used for spreading.     

Vaccinia virus RNA helicase: nucleic acid specificity in duplex unwinding.

Vaccinia virus RNA helicase (NPH-II) catalyzes nucleoside triphosphate-dependent unwinding of duplex RNAs containing a single-stranded 3' RNA tail. In this study, we examine the structural features of the nucleic acid substrate that are important for helicase activity. Strand displacement was affected by the length of the 3' tail. Whereas NPH-II efficiently unwound double-stranded RNA substrates with 19- or 11-nucleotide (nt) 3' tails, shortening the 3' tail to 4 nt reduced unwinding by an order of magnitude. Processivity of the helicase was inferred from its ability to unwind a tailed RNA substrate containing a 96-bp duplex region. NPH-II exhibited profound asymmetry in displacing hybrid duplexes composed of DNA and RNA strands. A 34-bp RNA-DNA hybrid with a 19-nt 3' RNA tail was unwound catalytically, whereas a 34-bp DNA-RNA hybrid containing a 19-nt 3' DNA tail was 2 orders of magnitude less effective as a helicase substrate. NPH-II was incapable of displacing a 34-bp double-stranded DNA substrate of identical sequence. 3'-Tailed DNA molecules with 24- or 19-bp duplex regions were also inert as helicase substrates. On the basis of current models for RNA-DNA hybrid structures, we suggest the following explanation for these findings. (i) Unwinding of duplex nucleic acids by NPH-II is optimal when the polynucleotide strand of the duplex along which the enzyme translocates has adopted an A-form secondary structure, and (ii) a B-form secondary structure impedes protein translocation through DNA duplexes.     

RNA pseudoknot prediction in energy-based models.

RNA molecules are sequences of nucleotides that serve as more than mere intermediaries between DNA and proteins, e.g., as catalytic molecules. Computational prediction of RNA secondary structure is among the few structure prediction problems that can be solved satisfactorily in polynomial time. Most work has been done to predict structures that do not contain pseudoknots. Allowing pseudoknots introduces modeling and computational problems. In this paper we consider the problem of predicting RNA secondary structures with pseudoknots based on free energy minimization. We first give a brief comparison of energy-based methods for predicting RNA secondary structures with pseudoknots. We then prove that the general problem of predicting RNA secondary structures containing pseudoknots is NP complete for a large class of reasonable models of pseudoknots.     

Distinguishable populations report on the interactions of single DNA molecules with solid-state nanopores

Solid-state nanopores have received increasing interest over recent years because of their potential for genomic screening and sequencing. In particular, small nanopores (2-5 nm in diameter) allow the detection of local structure along biological molecules, such as proteins bound to DNA or possibly the secondary structure of RNA molecules. In a typical experiment, individual molecules are translocated through a single nanopore, thereby causing a small deviation in the ionic conductance. A correct interpretation of these conductance changes is essential for our understanding of the process of translocation, and for further sophistication of this technique. Here, we present translocation measurements of double-stranded DNA through nanopores down to the diameter of the DNA itself (1.8-7 nm at the narrowest constriction). In contrast to previous findings on such small nanopores, we find that single molecules interacting with these pores can cause three distinct levels of conductance blockades. We attribute the smallest conductance blockades to molecules that briefly skim the nanopore entrance without translocating, the intermediate level of conductance blockade to regular head-to-tail translocations, and the largest conductance blockades to obstruction of the nanopore entrance by one or multiple (duplex) DNA strands. Our measurements are an important step toward understanding the conductance blockade of biomolecules in such small nanopores, which will be essential for future applications involving solid-state nanopores.     

The bacterial porin superfamily: sequence alignment and structure prediction

The porins of Gram-negative bacteria are responsible for the 'molecular sieve' properties of the outer membrane. They form large water-filled channels which allow the diffusion of hydrophilic molecules into the periplasmic space. Owing to the strong hydrophilicity of their amino acid sequence and the nature of their secondary structure (beta strands), conventional hydropathy methods for predicting membrane topology are useless for this class of protein. The large number of available porin amino acid sequences was exploited to improve the accuracy of the prediction in combination with tools detecting amphipathicity of secondary structure. Using the constraints of beta-sheet structure these porins are predicted to contain 16 membrane-spanning strands, 14 of which are common to the two (enteric and the neisserial) porin subfamilies.     

Secondary structure of the r(CUUCGG) tetraloop.

The oligonucleotide r(GGACUUCGGUCC) has been observed to adopt a hairpin conformation by solution NMR and a double helical conformation by X-ray diffraction. In order to understand this apparent conflict, we used time-resolved fluorescence depolarization and 19fluorine NMR to follow the secondary structure of this dodecamer as the solution composition was changed stepwise from the NMR experimental conditions to those used for crystallization. Calculation of the dodecamer concentration in the crystal (180 mM strands) and the cation concentration needed for neutrality (>2 M) prompted investigation of a tethered species, in which two dodecamers are connected by a string of 4 nt, geometrically equivalent to approximately 100 mM strands, in 2.5 M NaCl. The RNA tetraloop and its DNA analog maintain a single-strand hairpin conformation in solution, even under the conditions used to grow the crystal. Under high salt conditions, the tethered RNA and DNA analogs of this sequence yield secondary components which could be the double helical conformation. Crystal contacts in addition to solvent changes and high RNA concentrations are needed to obtain the double helix as the predominant species.