Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model.

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A -form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.     

Conformational transitions in thymidine bulge-containing deoxytridecanucleotide duplexes. Role of flanking sequence and temperature in modulating the equilibrium between looped out and stacked thymidine bulge states

Structural features at extra thymidine bulge sites in DNA duplexes have been elucidated from a two-dimensional NMR analysis of through-bond and through-space connectivities in the otherwise self-complementary d(C-C-G-T-G-A-A-T-T-C-C-G-G) (GTG 13-mer) and d(C-C-G-G-A-A-T-T-C-T-C-G-G) (CTC 13-mer) duplexes in aqueous solution. These studies establish that the extra thymidine flanked by guanosines in the GTG 13-mer duplex is in a conformational equilibrium between looped out and stacked states. The looped-out state is favored at low temperature (0 degrees C), whereas the equilibrium shifts in favor of the stacked state at elevated temperatures (35 degrees C) prior to the onset of the duplex-strand transition. By contrast, the extra thymidine flanked by cytidines in the CTC 13-mer duplex is looped out independent of temperature in the duplex state. Our results demonstrate that temperature and flanking sequence modulate the equilibrium between looped-out and stacked conformations of single base thymidine bulges in DNA oligomer duplexes.     

The structure of helix III in Xenopus oocyte 5 S rRNA: an RNA stem containing a two-nucleotide bulge

The solution structure of an oligonucleotide containing the helix III sequence from Xenopus oocyte 5 S rRNA has been determined by NMR spectroscopy. Helix III includes two unpaired adenosine residues, flanked on either side by G:C base-pairs, that are required for binding of ribosomal protein L5. The consensus conformation of helix III in the context provided by this oligonucleotide has the two adenosine residues located in the minor groove and stacked upon the 3' flanking guanosine residue, consistent with biochemical studies of free 5 S rRNA in solution. A distinct break in stacking that occurs between the first adenosine residue of the bulge and the flanking 5' guanosine residue exposes the base of the adenosine residue in the minor groove and the base of the guanosine residue in the major groove. The major groove of the helix is widened at the site of the unpaired nucleotides and the helix is substantially bent; nonetheless, the G:C base-pairs flanking the bulge are intact. The data indicate that there may be conformational heterogeneity centered in the bulge region. The corresponding adenosine residues in the Haloarcula marismortui 50 S ribosomal subunit form a dinucleotide platform, which is quite different from the motif seen in solution. Thus, the conformation of helix III probably changes when 5 S rRNA is incorporated into the ribosome.     

Investigation of potential RNA bulge stabilizing elements

As a part of our interest in recognition and cleavage of RNA we carried out thermal melting studies with the aim of screening a number of simple oligonucleotide modifications for their potential as modifying elements for RNA bulge stabilizing oligonucleotides. A specific model system from our studies on oligonucleotide-based artificial nuclease (OBAN) systems was chosen and the bulge size was varied from one to five nucleotides. Introduction of single 2'-modified nucleoside moieties (2'-O-methyl, 2'-deoxy and 2'-deoxy-2'-amino) with different conformational preferences adjacent to the bulge revealed that a higher preference for the north conformers gave more stable bulges across the whole range of bulge sizes. Changing a bulge closing a G-U wobble base pair to a G-C pair resulted in the interesting observation that, although the fully complementary complex and small bulges were highly stabilized, there was little difference in the stability of the larger bulges. The wobble base pair even gave a slight stabilization of the 5 nt bulge system. Introduction of a uridine C-5 linker with a single ammonium group was clearly bulge stabilizing (DeltaT(m) + 4.6 to + 5.4 degrees C for the three most stabilized bulges), although with limited selectivity for different bulge sizes since the fully complementary duplex was also stabilized. Introduction of a naphthoyl group on a 2'-aminolinker mostly gave a destabilizing effect, while introduction of a 5-aminoneocuproine moiety on the same linker resulted in stabilization of all bulges, in particular those with two or four unpaired nucleotides (DeltaT(m) + 3.6 and + 2.9 degrees C respectively). The aromatic groups destabilize the fully complementary duplex, resulting in higher selectivity towards stabilization of bulges. A combination of the studied partial element should be suitable for future designs of modified oligonucleotides that, apart from standard base pairing, can also provide additional non-Watson-Crick recognition of RNA.     

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations.

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.     

Kinetics of Bulge Bases in Small RNAs and the Effect of Pressure on It

Due to their self-catalytic properties, small RNAs with bulge bases are hypothesized to be primordial molecules which could form elementary translation systems. Using molecular dynamics simulations, we study the binding propensity of small RNAs by calculating the free energy barrier corresponding to the looped out conformations of bulge bases, which presumably act as the binding sites for ligands in these small RNAs. We find that base flipping kinetics can proceed at atmospheric pressure but with a very small propensity. Furthermore, the free energy barrier associated with base flipping depends on the stacking with neighboring bases. Next, we studied the base flipping kinetics with pressure. We find that the free energy associated with base looping out increases monotonically as the pressure is increased. Furthermore, we calculate the mean first-passage time of conformational looping out of the bulge base using the diffusion of reaction coordinate associated with the base flipping on the underlying free energy surface. We find that the mean first-passage time associated with bulge looping out increases slowly upon increasing pressures [Formula: see text] up to [Formula: see text] atm but changes dramatically for [Formula: see text] atm. Finally, we discuss our results in the light of the role of hydration shell of water around RNA. Our results are relevant for the RNA world hypothesis.     

Conserved nucleotides in the TAR RNA stem of human immunodeficiency virus type 1 are critical for Tat binding and trans activation: model for TAR RNA tertiary structure.

Interaction between the human immunodeficiency virus type 1 (HIV-1) trans-activator Tat and its cis-acting responsive RNA element TAR is necessary for activation of HIV-1 gene expression. We investigated the hypothesis that the essential uridine residue at position 23 in the bulge of TAR RNA is involved in intramolecular hydrogen bonding to stabilize an unique RNA structure required for recognition by Tat. Nucleotide substitutions in the two base pairs of the TAR stem directly above the essential trinucleotide bulge that maintain base pairing but change sequence prevent complex formation with Tat in vitro. Corresponding mutations tested in a trans-activation assay strongly affect the biological activity of TAR in vivo, suggesting an important role for these nucleotides in the Tat-TAR interaction. On the basis of these data, a model is proposed which implicates uridine 23 in a stable tertiary interaction with the GC pair directly above the bulge. This interaction would cause widening of the major groove of the RNA, thereby exposing its hydrogen-bonding surfaces for possible interaction with Tat. The model also predicts a gap between uridine 23 and the first base pair in the stem above, which would require one or more unpaired nucleotides to close, but does not predict any other role for such nucleotides. In accordance with this prediction, synthetic propyl phosphate linkers of equivalent length to 1 or 2 nucleotides, were found to be fully acceptable substitutes in the bulge above uridine 23, demonstrating that neither the bases nor the ribose moieties at these positions are implicated in the recognition of TAR RNA by Tat.     

Crystal structure of an RNA duplex r(gugucgcac)(2) with uridine bulges.

The crystal structure of a nonamer RNA duplex with a uridine bulge in each strand, r(gugucgcac)(2), was determined at 1.4 A resolution. The structure was solved by multiple anomalous diffraction phasing method using a three-wavelength data set collected at the Advanced Protein Source and refined to a final R(work)/R(free) of 21.2 %/23.4 % with 33,271 independent reflections (Friedel pairs unmerged). The RNA duplex crystallized in the tetragonal space group P4(1)22 with two independent molecules in the asymmetric unit. The unit cell dimensions are a=b=47.18 A and c=80.04 A. The helical region of the nonamer adopts the A-form conformation. The uridine bulges assume similar conformations, with uracils flipping out and protruding into the minor groove. The presence of the bulge induces very large twist angles (approximately +50 degrees) between the base-pairs flanking the bulges while causing profound kinks in the helix axis at the bulges. This severe twist and the large kink in turn produces a very narrow major groove at the middle of the molecule. The ribose sugars of the guanosines before the bulges adopt the C2'-endo conformation while the rest, including the bulges, are in the C3'-endo conformation. The intrastrand phosphate-phosphate (P-P) distance of the phosphate groups flanking the bulges (approximately 4.4 A) are significantly shorter than the average P-P distance in the duplex (6.0 A). This short distance between the two phosphate groups brings the non-bridging oxygen atoms close to each other where a calcium ion is bound to each strand. The calcium ions in molecule 1 are well defined while the calcium ions in molecule 2 are disordered.     

Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies

Asymmetric bulge loop motifs are widely dispersed in all types of functional RNAs. They are frequently occurring structural motifs in folded RNA structures and appear commonly in pre-microRNA and ribosomes, where they are involved in specific RNA–RNA and RNA–protein interactions. It is therefore necessary to understand such motifs from a structural point of view. We analyzed all available RNA structures and identified quite a few fragments of double helices that contain bulges. We found that these discontinuities often introduce kinks into the double helices, which also affects the stacking overlap between the base pairs across the irregularity. In order to understand the influence of these bulges on stability and flexibility, we carried out molecular dynamics simulations of three different single-residue bulge-containing RNA helices using the CHARMM36 force field. The structural variability at the junctions of RNA bulges is expected to differ from that in continuous double-helical stretches. The structural features of the junction region were observed to vary noticeably depending on the orientation of the bulge residue. When the base of the bulge residue is looped out, the RNA stretch behaves like a standard long A-form RNA double helix, whereas the entire RNA behaves differently when the base of the bulge residue is intercalated between base pairs inside the RNA stem. Such single-base intercalation was found to introduce a permanent kink into the composite double helix, which could be a recognition element for Dicer during the maturation of miRNA.     

New complex of post-activated neocarzinostatin chromophore with DNA: bulge DNA binding from the minor groove

Neocarzinostatin (NCS-chrom), a natural enediyne antitumor antibiotic, undergoes either thiol-dependent or thiol-independent activation, resulting in distinctly different DNA cleavage patterns. Structures of two different post-activated NCS-chrom complexes with DNA have been reported, revealing strikingly different binding modes that can be directly related to the specificity of DNA chain cleavage caused by NCS-chrom. The third structure described herein is based on recent studies demonstrating that glutathione (GSH) activated NCS-chrom efficiently cleaves DNA at specific single-base sites in sequences containing a putative single-base bulge. In this structure, the GSH post-activated NCS-chrom (NCSi-glu) binds to a decamer DNA, d(GCCAGAGAGC), from the minor groove. This binding triggers a conformational switch in DNA from a loose duplex in the free form to a single-strand, tightly folded hairpin containing a bulge adenosine embedded between a three base pair stem. The naphthoate aromatic moiety of NCSi-glu intercalates into a GG step flanked by the bulge site, and its substituent groups, the 2-N-methylfucosamine carbohydrate ring and the tetrahydroindacene, form a complementary minor groove binding surface, mostly interacting with the GCC strand in the duplex stem of DNA. The bulge site is stabilized by the interactions involving NCSi-glu naphthoate and GSH tripeptide. The positioning of NCSi-glu is such that only single-chain cleavage via hydrogen abstraction at the 5'-position of the third base C (which is opposite to the putative bulge base) in GCC is possible, explaining the observed single-base cleavage specificity. The reported structure of the NCSi-glu-bulge DNA complex reveals a third binding mode of the antibiotic and represents a new family of minor groove bulge DNA recognition structures. We predict analogue structures of NCSi-R (R = glu or other substituent groups) may be versatile probes for detecting the existence of various structures of nucleic acids. The NMR structure of this complex, in combination with the previously reported NCSi-gb-bulge DNA complex, offers models for specific recognition of DNA bulges of various sizes through binding to either the minor or the major groove and for single-chain cleavage of bulge DNA sequences.     

Characterization of the solution conformations of unbound and Tat peptide-bound forms of HIV-1 TAR RNA.

Basic peptides from the carboxy terminus of the HIV-1 Tat protein bind to the apical stem-loop region of TAR RNA with high affinity and moderate specificity. The conformations of the unbound and 24 residue Tat peptide (Tfr24)-bound forms of TAR RNA have been characterized by NMR spectroscopy. The unbound form of TAR exists in major and minor forms having different trinucleotide bulge conformations. A specific TAR RNA conformational change is observed upon complex formation with Tfr24, consisting of coaxial stacking of helical stems and base triple formation. A U23-A27-U38 base triple is proposed based on exchangeable proton NMR data, where U23 forms a base pair with A27 in the major groove. No evidence for base triple formation was found for Tat peptides in which lysine residues are extensively substituted for arginine.     

Thermodynamic and spectroscopic study of bulge loops in oligoribonucleotides

Thermodynamic parameters for bulge loops of one to three nucleotides in oligoribonucleotide duplexes have been measured by optical melting. The results indicate bulges Bn of An and Un have similar stabilities in the duplexes, GCGBmGCG + CGCCGC. The stability increment for a bulge depends on more than its adjacent base pairs. For example, the stability increment for a bulge is affected more than 1 kcal/mol by changing two nonadjacent base pairs or by adding terminal unpaired nucleotides (dangling ends) three base pairs away. Thus a nearest-neighbor approximation for helixes with bulges is oversimplified. Many of the non-self-complementary strands used in this study were observed to form homoduplexes. Such duplexes with GA mismatches were particularly stable.     

The global structures of a wild-type and poorly functional plant luteoviral mRNA pseudoknot are essentially identical

The helical junction region of a -1 frameshift stimulating hairpin-type mRNA pseudoknot from sugarcane yellow leaf virus (ScYLV) is characterized by a novel C27.(G7-C14) loop 2-stem 1 minor groove base triple, which is stacked on a C8+.(G12-C28) loop 1-stem 2 major groove base triple. Substitution of C27 with adenosine reduces frameshifting efficiency to a level just twofold above the slip-site alone. Here, we show that the global structure of the C27A ScYLV RNA is nearly indistinguishable from the wild-type counterpart, despite the fact that the helical junction region is altered and incorporates the anticipated isostructural A27.(G7-C14) minor groove base triple. This interaction mediates a 2.3-A displacement of C8+ driven by an A27 N6-C8+ O2 hydrogen bond as part of an A(n-1).C+.G-Cn base quadruple. The helical junction regions of the C27A ScYLV and the beet western yellows virus (BWYV) pseudoknots are essentially superimposable, the latter of which contains an analogous A25.(G7-C14) minor groove base triple. These results reveal that the global ground-state structure is not strongly correlated with frameshift stimulation and point to a reduced thermodynamic stability and/or enhanced kinetic lability that derives from an altered helical junction architecture in the C27A ScYLV RNA as a significant determinant for setting frameshifting efficiencies in plant luteoviral mRNA pseudoknots.     

The crystal structure of the octamer [r(guauaca)dC]2 with six Watson-Crick base-pairs and two 3' overhang residues

The crystal structure of an alternating RNA octamer, r(guauaca)dC (RNA bases are in lower case while the only DNA base is in upper case), with two 3' overhang residues one of them a terminal deoxycytosine and the other a ribose adenine, has been determined at 2.2 A resolution. The refined structure has an Rwork 18.6% and Rfree 26.8%. There are two independent duplexes (molecules I and II) in the asymmetric unit cell, a = 24.95, b = 45.25 and c = 73.67 A, with space group P2(1)2(1)2(1). Instead of forming a blunt end duplex with two a+.c mispairs and six Watson-Crick base-pairs, the strands in the duplex slide towards the 3' direction forming a two-base overhang (radC) and a six Watson-Crick base-paired duplex. The duplexes are bent (molecule I, 20 degrees; molecule II, 25 degrees) and stack head-to-head to form a right-handed superhelix. The overhang residues are looped out and the penultimate adenines of the two residues at the top end (A15) are anti and at the bottom (A7) end are syn. The syn adenine bases form minor groove A*(G.C) base triples with C8-H...N2 hydrogen bonds. The anti adenine in molecule II also forms a triple and a different C2-H...N3 hydrogen bond, while the other anti adenine in molecule I does not, it stacks on the looped out overhang base dC. The 3' terminal deoxycytosines form two stacked hemiprotonated trans d(C.C)+ base-pairs and the pseudo dyad related molecules form four consecutive deoxyribose and ribose zipper hydrogen bonds in the minor groove.     

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).     

Thermodynamic parameters for an expanded nearest-neighbor model for the formation of RNA duplexes with single nucleotide bulges

Thirty-four RNA duplexes containing single nucleotide bulges were optically melted, and the thermodynamic parameters deltaH degrees, deltaS degrees, deltaG degrees (37), and T(M) for each sequence were determined. Data from this study were combined with data from previous thermodynamic data [Longfellow, C. E., Kierzek, R., and Turner, D. H. (1990) Biochemistry 29, 278-85] to develop a model that will more accurately predict the free energy of an RNA duplex containing a single nucleotide bulge. Differences between purine and pyrimidine bulges as well as differences between Group I duplexes, those in which the bulge is not identical to either neighboring nucleotide, and Group II duplexes, those in which the bulge is identical to at least one neighboring nucleotide, were considered. The length of the duplex, non-nearest-neighbor effects, and bulge location were also examined. A model was developed which divides sequences into two groups: those with pyrimidine bulges and those with purine bulges. The proposed model for pyrimidine bulges predicts deltaG degrees (37,bulge) = 3.9 kcal/mol + 0.10deltaG degrees (37,nn) + beta, while the model for purine bulges predicts deltaG degrees (37,bulge) = 3.3 kcal/mol - 0.30deltaG degrees (37,nn) + beta, where beta has a value of 0.0 and -0.8 kcal/mol for Group I and Group II sequences, respectively, and deltaG degrees (37,nn) is the nearest-neighbor free energy of the base pairs surrounding the bulge. The conformation of bulge loops present in rRNA was examined. Three distinct families of structures were identified. The bulge loop was either extrahelical, intercalated, or in a "side-step" conformation.     

The effect of base sequence on the stability of RNA and DNA single base bulges

Forty-eight RNA duplexes were constructed that contained all common single base bulges at six different locations. The stabilities of the RNAs were determined by temperature gradient gel electrophoresis (TGGE). The relative stability of a single base bulge was dependent on both base identity and the nearest neighbor context. The single base bulges were placed into two categories. A bulged base with no identical neighboring base was defined as a Group I base bulge. Group II-bulged bases had at least one neighboring base identical to it. Group II bulges were generally more stable than Group I bulges in the same nearest neighbor environments. This indicates that position degeneracy of an unpaired base enhances stability. Differences in the mobility transition temperatures between the RNA fragments with bulges and the completely base-paired reference RNAs were related to free energy differences. Simple models for estimating the free energy contribution of single base bulges were evaluated from the free energy difference data. The contribution of a Group I bulge 5'-(XNZ)-3'.5'-(Z'-X')-3' where N is the unpaired base and X.X' and Z.Z' the neighboring base pairs, could be well-represented (+/-0.34 kcal/mol) by the equation, DeltaG((X)(N)()(Z))(.)((Z)(')(-)(X)(')()) = 3.11 + 0. 40DeltaG(s)()((XZ))(.)((Z)(')(X)(')()). DeltaG(s)()((XZ))(. )((Z)(')(X)(')()) is the stacking energy of the closing base pair doublet. By adding a constant term, delta = -0.3 kcal/mol, to the right side of the above equation, free energies of Group II bulges could also be predicted with the same accuracy. The term delta represents the stabilizing effect due to position degeneracy. A similar equation/model was applied to previous data from 32 DNA fragments with single base bulges. It predicted the free energy differences with a similar standard deviation.     

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

The way in which a single-base bulge might affect the structure of an RNA helix has been examined by preparing a series of six RNA hairpins, all with seven base pairs and a four-nucleotide loop. Five of the hairpins have single-base bulges at different positions. The intercalating cleavage reagent (methidiumpropyl)-EDTA-Fe(II) [MPE-Fe(II)] binds preferentially at a CpG sequence in the helix lacking a bulge and in four of the five hairpins with bulges. Hairpins with a bulge one or two bases to the 3' side of the CpG sequence bind ethidium 4-5-fold more strongly than the others. V1 RNase, which is sensitive to RNA backbone conformation in helices, detects a conformational change in all of the helices when ethidium binds; the most dramatic changes, involving the entire hairpin stem, are in one of the two hairpins with enhanced ethidium affinity. Only a slight conformational change is detected in the hairpin lacking a bulge. A bulge adjacent to a CpG sequence in a 100-nucleotide ribosomal RNA fragment enhances MPE-Fe(II) binding by an order of magnitude. These results extend our previous observations of bulges at a single position in an RNA hairpin [White, S. A., & Draper, D.E. (1987) Nucleic Acids Res. 15, 4049] and show that (1) a structural change in an RNA helix may be propagated for several base pairs, (2) bulges tend to increase the number of conformations available to a helix, and (3) the effects observed in small RNA hairpins are relevant to larger RNAs with more extensive structure. A bulge in a DNA hairpin identical in sequence with the RNA hairpins does not enhance MPE-Fe(II) binding affinity, relative to a control DNA hairpin. The effects of bulges on ethidium intercalation are evidently modulated by helix structure.     

Dynamic ensemble view of the conformational landscape of HIV-1 TAR RNA and allosteric recognition

RNA conformational dynamics and the resulting structural heterogeneity play an important role in RNA functions, e.g., recognition. Recognition of HIV-1 TAR RNA has been proposed to occur via a conformational capture mechanism. Here, using ultrafast time-resolved fluorescence spectroscopy, we have probed the complexity of the conformational landscape of HIV-1 TAR RNA and monitored the position-dependent changes in the landscape upon binding of a Tat protein-derived peptide and neomycin B. In the ligand-free state, the TAR RNA samples multiple families of conformations with various degrees of base stacking around the three-nucleotide bulge region. Some subpopulations partially resemble those ligand-bound states, but the coaxially stacked state is below the detection limit. When Tat or neomycin B binds, the bulge region as an ensemble undergoes a conformational transition in a position-dependent manner. Tat and neomycin B induce mutually exclusive changes in the TAR RNA underlying the mechanism of allosteric inhibition at an ensemble level with residue-specific details. Time-resolved anisotropy decay measurements revealed picosecond motions of bases in both ligand-free and ligand-bound states. Mutation of a base pair at the bulge--stem junction has differential effects on the conformational distributions of the bulge bases. A dynamic model of the ensemble view of the conformational landscape for HIV-1 TAR RNA is proposed, and the implication of the general mechanism of RNA recognition and its impact on RNA-based therapeutics are discussed.