Search

A dynamic programming algorithm for prediction of RNA secondary structure has been revised to accommodate folding constraints determined by chemical modification and to include free energy increments for coaxial stacking of helices when they are either adjacent or separated by a single mismatch. Furthermore, free energy parameters are revised to account for recent experimental results for terminal mismatches and hairpin, bulge, internal, and multibranch loops. To demonstrate the applicability of this method, in vivo modification was performed on 5S rRNA in both Escherichia coli and Candida albicans with 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate, dimethyl sulfate, and kethoxal. The percentage of known base pairs in the predicted structure increased from 26.3% to 86.8% for the E. coli sequence by using modification constraints. For C. albicans, the accuracy remained 87.5% both with and without modification data. On average, for these sequences and a set of 14 sequences with known secondary structure and chemical modification data taken from the literature, accuracy improves from 67% to 76%. This enhancement primarily reflects improvement for three sequences that are predicted with <40% accuracy on the basis of energetics alone. For these sequences, inclusion of chemical modification constraints improves the average accuracy from 28% to 78%. For the 11 sequences with <6% pseudoknotted base pairs, structures predicted with constraints from chemical modification contain on average 84% of known canonical base pairs.

Candida albicans is a significant cause of disease in immunocompromised humans. Because the number of people infected by fungal pathogens is increasing, strategies are being developed to target RNAs in fungi. This work shows that oligonucleotides can serve as therapeutics against C. albicans. In particular, oligonucleotides are taken up from cell culture medium in an energy-dependent process. After uptake, oligonucleotides, including RNA, remain mostly intact after 12 h in culture. For culture conditions designed for mammalian cells, intracellular concentrations of oligonucleotides in C. albicans exceed those in COS-7 mammalian cells, suggesting that uptake can provide selective targeting of fungi over human cells. A 19-mer 2'OMe (oligonucleotide with a 2'-O-methyl backbone) hairpin is described that inhibits growth of a C. albicans strain at pH < 4.0. This pH is easily tolerated in some parts of the body subject to C. albicans infections. In vivo dimethyl sulfate modification of ribosomal RNA and the decreased rate of protein synthesis suggest that this hairpin's activity may be due to targeting the ribosome in a way that does not depend on base pairing. Addition of anti-C. albicans oligonucleotides to COS-7 mammalian cells has no effect on cell growth. Evidently, oligonucleotides can selectively serve as therapeutics toward C. albicans and, presumably, other pathogens. Information from genome sequencing and functional genomics studies on C. albicans and other pathogens should allow rapid design and testing of other approaches for oligonucleotide therapies.

RNA is becoming an important therapeutic target. Many potential RNA targets require secondary or tertiary structure for function. Examples include ribosomal RNAs, RNase P RNAs, mRNAs with untranslated regions that regulate translation, and group I and group II introns. Here, a method is described to inhibit RNA function by exploiting the propensity of RNA to adopt multiple folded states that are of similar free energy. This method, called oligonucleotide directed misfolding of RNA (ODMiR), uses short oligonucleotides to stabilize inactive structures. The ODMiR method is demonstrated with the group I intron from Candida albicans, a human pathogen. The oligonucleotides, (L)(TACCTTTC) and T(L)CT(L)AC(L)GA(L)CG(L)GC(L)C, with L denoting a locked nucleic acid residue, inhibit 50% of group I intron splicing in a transcription mixture at about 150 and 30 nM oligonucleotide concentration, respectively. Both oligonucleotides induce misfolds as determined by native gel electrophoresis and diethyl pyrocarbonate modification. The ODMiR approach provides a potential therapeutic strategy applicable to RNAs with secondary or tertiary structures required for function.

Binding enhancement by tertiary interactions is a strategy that takes advantage of the higher order folding of functionally important RNAs to bind short nucleic acid-based compounds tightly and more specifically than possible by simple base pairing. For example, tertiary interactions enhance binding of specific hexamers to a group I intron ribozyme from the opportunistic pathogen Pneumocystis carinii by 1,000- to 100,000-fold relative to binding by only base pairing. One such hexamer, d(AnTnGnAnCn)rU, contains an N3' --> P5' phosphoramidate deoxysugar-phosphate backbone (n) that is resistant to chemical and enzymatic decay. Here, it is shown that this hexamer is also a suicide inhibitor of the intron's self-splicing reaction in vitro. The hexamer is ligated in trans to the 3' exon of the precursor, producing dead-end products. At 4 mM Mg2+, the fraction of trans-spliced product is greater than normally spliced product at hexamer concentrations as low as 200 nM. This provides an additional level of specificity for compounds that can exploit the catalytic potential of complexes with RNA targets.

An RNA model system consisting of an oligomer binding to a 4-nt overhang at the 5' end of a hairpin stem provides thermodynamic parameters for helix-helix interfaces. In a sequence-dependent manner, oligomers bind up to 1000-fold more tightly adjacent to the hairpin stem than predicted for binding to a free tetramer at 37 degrees C. For the interface (/) in [formula: see text] additional free energy change, delta delta G 37 degrees, for binding is roughly the nearest-neighbor delta G 37 degrees for propagation of an uninterrupted helix of equivalent sequence, CGGC. When X and Z are omitted, the delta delta 37 degrees is even more favorable by approximately 1 kcal/mol (1 cal = 4.184J). On average, predictions of 11 RNA secondary structures improve from 67 to 74% accuracy by inclusion of similar stacking contributions.

The effects of guanosine 5'-monophosphate and 2'-deoxyguanosine 5'-monophosphate on the thermodynamics and kinetics of pyrene-labeled 5' exon mimic (pyCUCU) binding to the catalytic RNA (ribozyme) from Tetrahymena thermophila have been determined by fluorescence titration and kinetics experiments at 15 degrees C. pyCUCU binding to L-21 Sca I-truncated ribozyme is weaker by a factor of 5 in the presence of saturating guanosine 5'-monophosphate, whereas it is 4-fold stronger in the presence of saturating 2'-deoxyguanosine 5'-monophosphate. Results from kinetics experiments suggest that anticooperative effects in the presence of guanosine 5'-monophosphate arise primarily from slower formation of tertiary contacts between the catalytic core of the ribozyme and the P1 duplex formed by pyCUCU and GGAGGG of the ribozyme. Conversely, cooperative effects in the presence of 2'-deoxyguanosine 5'-monophosphate arise primarily from slower disruption of tertiary contacts between the catalytic core of the ribozyme and the P1 duplex. Additional experiments suggest that these cooperative and anticooperative effects are not a function of the pyrene label, are not caused by a salt effect, and are not specific to one renaturation procedure for the ribozyme.

The accuracy of computer predictions of RNA secondary structure from sequence data and free energy parameters has been increased to roughly 70%. Performance is judged by comparison with structures known from phylogenetic analysis. The algorithm also generates suboptimal structures. On average, the best structure within 10% of the lowest free energy contains roughly 90% of phylogenetically known helixes. The algorithm does not include tertiary interactions or pseudoknots and employs a crude model for single-stranded regions. The only favorable interactions are base pairing and stacking of terminal unpaired nucleotides at the ends of helixes. The excellent performance is consistent with these interactions being the primary interactions determining RNA secondary structure.

Thermodynamic parameters for prediction of RNA duplex stability are reported. One parameter for duplex initiation and 10 parameters for helix propagation are derived from enthalpy and free-energy changes for helix formation by 45 RNA oligonucleotide duplexes. The oligomer sequences were chosen to maximize reliability of secondary structure predictions. Each of the 10 nearest-neighbor sequences is well-represented among the 45 oligonucleotides, and the sequences were chosen to minimize experimental errors in delta GO at 37 degrees C. These parameters predict melting temperatures of most oligonucleotide duplexes within 5 degrees C. This is about as good as can be expected from the nearest-neighbor model. Free-energy changes for helix propagation at dangling ends, terminal mismatches, and internal G X U mismatches, and free-energy changes for helix initiation at hairpin loops, internal loops, or internal bulges are also tabulated.