Effects of single-base bulges on intercalator binding to small RNA and DNA hairpins and a ribosomal RNA fragment

An asymmetric mesoscopic model for single bulges in RNA

Simple one-dimensional DNA or RNA mesoscopic models are of interest for their computational efficiency while retaining the key elements of the molecular interactions. However, they only deal with perfectly formed DNA or RNA double helices and consider the intra-strand interactions to be the same on both strands. This makes it difficult to describe highly asymmetric structures such as bulges and loops and, for instance, prevents the application of mesoscopic models to determine RNA secondary structures. Here we derived the conditions for the Peyrard-Bishop mesoscopic model to overcome these limitations and applied it to the calculation of single bulges, the smallest and simplest of these asymmetric structures. We found that these theoretical conditions can indeed be applied to any situation where stacking asymmetry needs to be considered. The full set of parameters for group I RNA bulges was determined from experimental melting temperatures using an optimization procedure, and we also calculated average opening profiles for several RNA sequences. We found that guanosine bulges show the strongest perturbation on their neighboring base pairs, considerably reducing the on-site interactions of their neighboring base pairs.

The HIV-2 Rev-response element: determining secondary structure and defining folding intermediates

Interaction between the viral protein Rev and the RNA motifs known as Rev response elements (RREs) is required for transport of unspliced and partially spliced human immunodeficiency virus (HIV)-1 and HIV-2 RNAs from the nucleus to the cytoplasm during the later stages of virus replication. A more detailed understanding of these nucleoprotein complexes and the host factors with which they interact should accelerate the development of new antiviral drugs targeting cis-acting RNA regulatory signals. In this communication, the secondary structures of the HIV-2 RRE and two RNA folding precursors have been identified using the SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) chemical probing methodology together with a novel mathematical approach for determining the secondary structures of RNA conformers present in a mixture. A complementary chemical probing technique was also used to support these secondary structure models, to confirm that the RRE2 RNA undergoes a folding transition and to obtain information about the relative positioning of RRE2 substructures in three dimensions. Our analysis collectively suggests that the HIV-2 RRE undergoes two conformational transitions before assuming the energetically most favorable conformer. The 3D models for the HIV-2 RRE and folding intermediates are also presented, wherein the Rev-binding stem–loops (IIB and I) are located coaxially in the former, which is in agreement with previous models for HIV-1 Rev-RRE binding.

Structural basis for bifunctional zinc(II) macrocyclic complex recognition of thymine bulges in DNA

The zinc(II) complex of 1-(4-quinoylyl)methyl-1,4,7,10-tetraazacyclododecane (cy4q) binds selectively to thymine bulges in DNA and to a uracil bulge in RNA. Binding constants are in the low-micromolar range for thymine bulges in the stems of hairpins, for a thymine bulge in a DNA duplex, and for a uracil bulge in an RNA hairpin. Binding studies of Zn(cy4q) to a series of hairpins containing thymine bulges with different flanking bases showed that the complex had a moderate selectivity for thymine bulges with neighboring purines. The dissociation constants of the most strongly bound Zn(cy4q)-DNA thymine bulge adducts were 100-fold tighter than similar sequences with fully complementary stems or than bulges containing cytosine, guanine, or adenine. In order to probe the role of the pendent group, three additional zinc(II) complexes containing 1,4,7,10-tetraazacyclododecane (cyclen) with aromatic pendent groups were studied for binding to DNA including 1-(2-quinolyl)methyl-1,4,7,10-tetraazacyclododecane (cy2q), 1-(4-biphenyl)methyl-1,4,7,10-tetraazacyclododecane (cybp), and 5-(1,4,7,10-tetraazacyclododecan-1-ylsulfonyl)-N,N-dimethylnaphthalen-1-amine (dsc). The Zn(cybp) complex binds with moderate affinity but little selectivity to DNA hairpins with thymine bulges and to DNA lacking bulges. Similarly, Zn(dsc) binds weakly both to thymine bulges and hairpins with fully complementary stems. The zinc(II) complex of cy2q has the 2-quinolyl moiety bound to the Zn(II) center, as shown by (1)H NMR spectroscopy and pH-potentiometric titrations. As a consequence, only weak (500 μM) binding is observed to DNA with no appreciable selectivity. An NMR structure of a thymine-bulge-containing hairpin shows that the thymine is extrahelical but rotated toward the major groove. NMR data for Zn(cy4q) bound to DNA containing a thymine bulge is consistent with binding of the zinc(II) complex to the thymine N3(-) and stacking of the quinoline on top of the thymine. The thymine-bulge bound zinc(II) complex is pointed into the major groove, and there are interactions with the guanine positioned 5' to the thymine bulge.

Native-like RNA tertiary structures using a sequence-encoded cleavage agent and refinement by discrete molecular dynamics

The difficulty of analyzing higher order RNA structure, especially for folding intermediates and for RNAs whose functions require domains that are conformationally flexible, emphasizes the need for new approaches for modeling RNA tertiary structure accurately. Here, we report a concise approach that makes use of facile RNA structure probing experiments that are then interpreted using a computational algorithm, carefully tailored to optimize both the resolution and refinement speed for the resulting structures, without requiring user intervention. The RNA secondary structure is first established using SHAPE chemistry. We then use a sequence-directed cleavage agent, which can be placed arbitrarily in many helical motifs, to obtain high quality inter-residue distances. We interpret this in-solution chemical information using a fast, coarse grained, discrete molecular dynamics engine in which each RNA nucleotide is represented by pseudoatoms for the phosphate, ribose, and nucleobase groups. By this approach, we refine base paired positions in yeast tRNA(Asp) to 4 A rmsd without any preexisting information or assumptions about secondary or tertiary structures. This blended experimental and computational approach has the potential to yield native-like models for the diverse universe of functionally important RNAs whose structures cannot be characterized by conventional structural methods.

Binding of helix-threading peptides to E. coli 16S ribosomal RNA and inhibition of the S15-16S complex

Helix-threading peptides (HTPs) constitute a new class of small molecules that bind selectively to duplex RNA structures adjacent to helix defects and project peptide functionality into the dissimilar duplex grooves. To further explore and develop the capabilities of the HTP design for binding RNA selectively, we identified helix 22 of the prokaryotic ribosomal RNA 16S as a target. This helix is a component of the binding site for the ribosomal protein S15. In addition, the S15-16S RNA interaction is important for the ordered assembly of the bacterial ribosome. Here we present the synthesis and characterization of helix-threading peptides that bind selectively to helix 22 of E. coli 16S RNA. These compounds bind helix 22 by threading intercalation placing the N termini in the minor groove and the C termini in the major groove. Binding is dependent on the presence of a highly conserved purine-rich internal loop in the RNA, whereas removal of the loop minimally affects binding of the classical intercalators ethidium bromide and methidiumpropyl-EDTAFe (MPEFe). Moreover, binding selectivity translates into selective inhibition of formation of the S15-16S complex.

RNA binding and thiolytic stability of a quinoline-containing helix-threading peptide

Helix-threading peptides (HTPs) bind selectively to sites predisposed to intercalation in folded RNA molecules placing peptide functional groups into the dissimilar grooves of the duplex. Here we report the design and synthesis of new HTPs with quinoline as the intercalation domain. A quinoline-containing HTP is shown to bind selectively to duplex RNA binding sites. Furthermore, the affinity cleavage pattern generated using an EDTA.Fe modified derivative is consistent with minor groove localization of its N-terminus. This compound binds base-pair steps flanked by single nucleotide bulges on the 3' side on both strands, whereas bulges on the 5' side of the intercalation site do not support binding. Furthermore, unlike acridine HTPs, the quinoline compound is resistant to thiolytic degradation that leads to loss of RNA-binding activity. The RNA-binding selectivity and stability observed for quinoline-containing HTPs make them excellent candidates for further development as regulators of intracellular RNA function.

Inhibitory properties of nucleic acid-binding ligands on protein synthesis

The use of small molecule inhibitors in the study of cellular processes is a powerful approach to understanding gene function. During the course of a high throughput screen for novel inhibitors of eukaryotic translation, we identified a number of nucleic acid binding ligands that showed activity in our assay. When tested on a panel of mRNA transcripts displaying different modes of translation initiation, these ligands showed a range of biological activities--with some inhibiting both cap-dependent and internal initiation and others preferentially blocking internal initiation. We used this information to identify a novel threading intercalator that inhibits Hepatitis C virus internal initiation.

The Thermodynamics of DNA Structural Motifs

DNA secondary structure plays an important role in biology, genotyping diagnostics, a variety of molecular biology techniques, in vitro-selected DNA catalysts, nanotechnology, and DNA-based computing. Accurate prediction of DNA secondary structure and hybridization using dynamic programming algorithms requires a database of thermodynamic parameters for several motifs including Watson-Crick base pairs, internal mismatches, terminal mismatches, terminal dangling ends, hairpins, bulges, internal loops, and multibranched loops. To make the database useful for predictions under a variety of salt conditions, empirical equations for monovalent and magnesium dependence of thermodynamics have been developed. Bimolecular hybridization is often inhibited by competing unimolecular folding of a target or probe DNA. Powerful numerical methods have been developed to solve multistate-coupled equilibria in bimolecular and higher-order complexes. This review presents the current parameter set available for making accurate DNA structure predictions and also points to future directions for improvement.

Two-dimensional conformation-dependent electrophoresis (2D-CDE) to separate DNA fragments containing unmatched bulge from complex DNA samples

DNA fragments containing mispaired and modified bases, bulges, lesions and specific sequences have altered conformation. Methods for separating complex samples of DNA fragments based on conformation but independent of length have many applications, including (i) separation of mismatched or unmatched DNA fragments from those perfectly matched; (ii) simultaneous, diagnostic, mismatch scanning of multiple fragments; (iii) isolation of damaged DNA fragments from undamaged fragments; and (iv) estimation of reannealing efficiency of complex DNA samples. We developed a two-dimensional conformation-dependent electrophoresis (2D-CDE) method for separating DNA fragments based on length and conformation in the first dimension and only on length in the second dimension. Differences in migration velocity due to conformation were minimized during second dimension electrophoresis by introducing an intercalator. To test the method, we constructed 298 bp DNA fragments containing cytosine bulges ranging from 1 to 5 nt. Bulge-containing DNA fragments had reduced migration velocity in the first dimension due to altered conformation. After 2D-CDE, bulge-containing DNA fragments had migrated in front of an arc comprising heterogeneous fragments with regular conformation. This simple and robust method could be used in both analytical and preparative applications involving complex DNA samples.

Binding of Hoechst 33258 to the TAR RNA of HIV-1. Recognition of a pyrimidine bulge-dependent structure

The transactivation response region (TAR) RNA is an essential component in transcriptional regulation of the human immunodeficiency virus type-1 (HIV-1) genome. We have examined the interaction between TAR RNA and the bisbenzimidazole derivative Hoechst 33258. Previous studies have shown that this drug, which is well known as an AT-selective DNA minor groove binder, can also interact with GC-rich sequences in DNA as well as with RNA, possibly by intercalation. Absorption spectroscopy, circular dichroism and electric linear dichroism, as well as RNase A footprinting, were employed to compare binding of Hoechst 33258 to wild-type RNA and its analogue lacking the pyrimidine bulge. The uridine bulge, which is an important contributor to the structural stability of TAR, plays an essential role in drug binding. Deletion of the bulge destabilizes both free and drug-bound forms of TAR and markedly affects the orientation of the drug chromophore complexed with the RNA. According to the linear dichroism data, the bisbenzimidazole is oriented more or less perpendicular to the RNA helix axis. The data are compatible with a model in which the bisbenzimidazole chromophore is inserted into the existing cavity created by the pyrimidine bulge. The footprinting experiments, showing that the drug binds to a unique site opposite the unpaired uridine residues, also support this model. The binding of Hoechst 33258 to the sequence 5'-GCUCU, which delimits the cavity, reflects the greater accessibility of that region compared with other sites in the RNA molecule. The identification of a binding site for small molecules in TAR offers promising perspectives for developing drugs that would block the access of TAR RNA to proteins and therefore for the design of anti-HIV agents.

Quantification of total RNA by ethidium bromide fluorescence may not accurately reflect the RNA mass

The fluorescent signal from purified large and small ribosomal RNA subunits stained with ethidium bromide was quantified using agarose gels and image analysis. A significantly smaller fluorescent signal was observed for 23 S rRNA compared to 16 S rRNA when normalized to mass. Therefore, small changes in amounts of 23 S or 16 S rRNA within a sample can have a profound impact on observed bulk fluorescence. A comparison of 23 S/16 S rRNA fluorescence ratios for 4 marine bacteria grown in batch culture indicates that the lower fluorescence of 23 S rRNA is widespread. The 23 S/16 S rRNA fluorescence ratios also suggest that intracellular concentrations of ribosomal subunits do not always adhere to a stoichiometry of 2.0.

Crystal structures of an A-form duplex with single-adenosine bulges and a conformational basis for site-specific RNA self-cleavage

BACKGROUND

Bulged nucleotides are common secondary structural motifs in RNA molecules and are often involved in RNA-RNA and RNA-protein interactions. RNA is selectively cleaved at bulge sites (when compared to other sites within stems) in the presence of divalent metal cations. The effects of bulge nucleotides on duplex stability and topology have been extensively investigated, but no detailed X-ray structures of bulge-containing RNA fragments have been available.

RESULTS

We have crystallized a self-complementary RNA-DNA chimeric 11-nucleotide sequence containing single-adenosine bulges under two different conditions, giving two distinct crystal forms. In both lattices the adenosines are looped out, leaving the stacking interactions in the duplex virtually unaffected. The bulges cause the duplex to kink in both cases. In one of the structures, the conformation of the bulged nucleotide places its modeled 2'-oxygen in line with the adjacent phosphate on the 3' side, where it is poised for nucleophilic attack.

CONCLUSIONS

Single adenosine bulges cause a marked opening of the normally narrow RNA major groove in both crystal structures, rendering the bases more accessible to interacting molecules compared with an intact stem. The geometries around the looped-out adenosines are different in the two crystal forms, indicating that bulges can confer considerable local plasticity on the usually rigid RNA double helix. The results provide a conformational basis for the preferential, metal-assisted self-cleavage of RNA at bulged sites.

Mapping RNA regions in eukaryotic ribosomes that are accessible to methidiumpropyl-EDTA.Fe(II) and EDTA.Fe(II)

Methidiumpropyl-EDTA.Fe(II) [MPE.Fe(II)] and EDTA.Fe(II) were used to investigate the structure of Drosophila melanogaster ribosomes. Cleavage reactions were performed on intact ribosomes in cell lysates in vitro and analyzed by primer extension with reverse transcriptase using oligodeoxynucleotide primers. Regions of 18S and 28S ribosomal RNAs (rRNAs) which are accessible to MPE.Fe(II) and EDTA.Fe(II) are located almost exclusively within expansion segments. The accessibility of these regions to cleavage indicates that they are likely exposed on the surface of eukaryotic ribosomes. These results provide information about the overall tertiary structure of rRNA in ribosomes.

Encephalomyocarditis virus internal ribosomal entry site RNA-protein interactions

Translational initiation of encephalomyocarditis virus (EMCV) mRNA occurs by ribosomal entry into the 5' nontranslated region of the EMCV mRNA, rather than by ribosomal scanning. Internal ribosomal binding requires a cis-acting element termed the internal ribosomal entry site (IRES). IRES elements have been proposed to be involved in the translation of picornavirus mRNAs and some cellular mRNAs. Internal ribosome binding likely requires the interaction of trans-acting factors that recognize both the mRNA and the ribosomal complex. Five cellular proteins (p52, p57, p70, p72, and p100) cross-link the EMCV IRES or fragments of the IRES. For one of these proteins, p57, binding to the IRES correlates with translation. Recently, p57 was identified to be very similar, if not identical, to polypyrimidine tract-binding protein. On the basis of cross-linking results with 21 different EMCV IRES fragments and cytoplasmic HeLa extract or rabbit reticulocyte lysate as the source of polypeptides, consensus binding sites for p52, p57, p70, and p100 are proposed. It is suggested that each of these proteins recognizes primarily a structural feature of the RNA rather than a specific sequence.

Implication of RNA structure on antisense oligonucleotide hybridization kinetics

A 47-nucleotide transcript of the activated Ha-ras gene was prepared and determined, by enzymatic structure mapping, to form a stable hairpin structure. Six antisense decaribonucleotides were designed, and association constants (Ka) for the hairpin- and length-matched complements were measured. Two of the antisense oligonucleotides targeted to the loop had nearly equal affinity for the transcript compared to the complement. The others, including one oligonucleotide complementary to the 3' side of the single-stranded loop, bound 10(5)-10(6)-fold less tightly to the transcript than to the short complement. We propose the difference in affinity is due to the target structure, both the secondary structure of the stem and the structure in the loop. Measurement of the bimolecular association rate constant, k1, and the dissociation rate constant, k-1, for these oligonucleotides indicates the observed relationship between affinity and structure is primarily due to k1.

A molecular mechanics investigation of RNA complexes. I. Ethidium intercalation in an HIV-1 TAR RNA sequence with an unpaired adenosine

Nucleic acid complexes with ethidium intercalated into different sites in a segment of HIV-1 TAR RNA with an unpaired A base, along with corresponding complexes with a normal RNA sequence without an unpaired base were studied by molecular mechanics energy minimization methods. Different intercalation geometries as well as different orientations of the ethidium molecule in the intercalation sites were tested. A general binding affinity enhancement for the ethidium binding to the bulge sequence compared with the normal RNA segment was obtained. With the unpaired adenosine base stacked in the duplex, the binding site adjacent to the 3' side of the bulge was found to be the most energetically favorable binding site, and the intercalation site 5' to the bulge in the same sequence is much less favorable. Unique correlated backbone conformational changes on binding of ethidium to the intercalation site 3' to the bulge were found to relieve backbone strains caused by the stacking of the unpaired base into the helix. These backbone conformational changes present a plausible molecular basis for the experimentally observed ethidium binding preference in this bulge RNA segment (L.S. Ratmeyer, R. Vinayak, G. Zon and W.D. Wilson, J. Med. Chem. 35, 966, 1992).

Modeling of nucleic acid complexes with cationic ligands: a specialized molecular mechanics force field and its application

A potential energy force field designed for modeling nucleic acids and particularly their complexes with cationic ligands is presented. The force field is a modified version of that developed by Weiner, S.J., Kollman, P.A., Nguyen, D.T. and Case, D.A.,J. Comp. Chem. 7,230-252 (1986) and is based upon the use of a distance dependent dielectric constant, epsilon = 4rij, and partially neutralized phosphates to represent solvent and counterion. Changes from the Weiner et al. force field include additional atom types and modifications to van der Waals, electrostatic, hydrogen bonding and torsional parameters. Molecular modeling test cases of the force field are presented for a number of simple small molecules, as well as uracil and benzene dimerization, thymine-adenine and cytosine-guanine base pair formation, and adenosine/deoxyadenosine pseudorotation. Several DNA and RNA oligomers and DNA/RNA intercalation complexes with ethidium are also modeled with the force field. In all cases, the modeling results compare favorably with available experimental results. Additionally, conformational trends observed experimentally for nucleic acids by NMR and X-ray crystallographic techniques are reproduced. The modeling results for ethidium intercalation indicate a complex in which the favorable interactions are primarily van der Waals contacts, and in which electrostatic interactions are a relatively minor component. We feel the force field is particularly useful for molecular mechanics aided drug design, and an analysis of modeling results with respect to design of drugs which bind selectively to RNA is presented.

Effects of mutations of the bulged nucleotide in the conserved P7 pairing element of the phage T4 td intron on ribozyme function

The P7 element of group I introns contains a semiconserved "bulged" nucleotide, a C in group IA introns (nt 870 in the td intron) and an A in group IB introns [Cech, T.R. (1988) Gene 73, 259-271]. Variants U870, G870, and A870, isolated by a combination of in vitro and in vivo genetic strategies, indicate that C and A at position 870 are consistent with splicing whereas U and G are not. Although mutants G870 and U870 could be activated in vitro by increasing the Mg2+ concentration, their Km for GTP at pH 7 was 20-100-fold elevated, and they were unable to undergo site-specific hydrolysis. The dependence of the mutants on high guanosine concentrations could be substantially overcome by an increase in pH, suggesting that a tautomeric change, which makes U and G mimic C and A, is responsible for restoring function. In contrast to the striking Km effect, Vmax for the mutants differed by less than a factor of 2 from the wild type. Furthermore, streptomycin, an aminoglycoside antibiotic that competes with guanosine for its binding site, inhibited splicing of the U870 and G870 constructs at least as well as of the C870 and A870 variants, indicating that the guanosine-binding site of the mutants is proficient at interacting with a guanidino group. While our experiments argue against a hydrogen-bonding interaction between the C6-O of the cofactor and C4-NH2 of the bulged nucleotide, they are consistent with other models in which the C4-NH2 and/or N3 groups of the bulged C are involved in establishing an active ribozyme.

Effects of proflavin and photoactivated proflavin on the template function of single-stranded DNA

DNA context-specific effects of the association of proflavin, single-stranded DNA and DNA polymerase on DNA polymerization reactions were examined. Frameshift mutations induced by the presence of proflavin during in vitro DNA replication of a single-stranded DNA template by the Klenow fragment of Escherichia coli DNA polymerase I were sequenced. More than 80% of the frameshifts were one base-pair deletions opposite purine bases that were immediately 3' to pyrimidines. Purines (Pu) that were not adjacent to pyrimidines (Py) were not deletion sites. The remaining deletions were opposite template pyrimidines that were also immediately 3' to another pyrimidine. All pyrimidine site deletions occurred in the context 5' PyPyPu 3'. In additional experiments, the site-specific inhibition of processive DNA polymerization by proflavin was examined. A novel inhibition of polymerization was found opposite all pyrimidines in the template when proflavin-template complexes were exposed to ten seconds of white light. This inhibition of polymerization is reversible. Longer photoactivation led to an altered pattern of DNA sequence-specific inhibition that was not reversible. The role of DNA sequence-specific interactions of proflavin with DNA in proflavin mutagenesis is discussed.

Thermodynamics of protein-RNA recognition in a highly conserved region of the large-subunit ribosomal RNA

Ribosomal protein L11 from Escherichia coli specifically binds to a highly conserved region of 23S ribosomal RNA. The thermodynamics of forming a complex between this protein and several different rRNA fragments have been investigated, by use of a nitrocellulose filter binding assay. A 57-nucleotide region of the RNA (C1052-U1108) contains all the protein recognition features, and an RNA fragment containing this region binds L11 10(3)-10(4)-fold more tightly than tRNA. Binding constants are on the order of 10 microM-1 and are only weakly dependent on K+ concentration (delta log K/delta log [K+] = -1.4) or temperature. Binding requires multivalent cations; Mg2+ is taken up into the complex with an affinity of approximately 3 mM-1. Other multivalent cations tested, Ca2+ and Co(NH3)63+, promote binding nearly as well. The pH dependence of binding is a bell-shaped curve with a maximum near neutral pH, but the entire curve is shifted to higher pH for the smaller of two RNA fragments tested. This result suggests that the smaller fragment favors a conformation stabilizing protonated forms of the RNA recognition site and is potentially relevant to a hypothesis that this rRNA region undergoes an ordered series of conformational changes during the ribosome cycle.

Scission of RNA by the chemical nuclease of 1,10-phenanthroline-copper ion: preference for single-stranded loops

The scission of RNA by the chemical nuclease activity of 1,10-phenanthroline-copper (OP-Cu) has been studied using a lac mRNA fragment and tRNAphe as substrates. Since the chemical mechanism of scission involves oxidative attack on the ribose, scission is observed at all nucleotides including dihydrouridine and Y-bases. Specificity for single-stranded loop regions is apparent from the similarity of the reactivity of OP-Cu to the single-strand specific reagents dimethyl sulfate and diethyl pyrocarbonate using the fragment of lac mRNA as a substrate. Similar preference is observed in the reaction with tRNA although scission in the helical acceptor stem is also observed.