Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study

Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.     

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.     

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.     

Contribution of the trifluoroacetyl group in the thermodynamics of antigen–antibody binding

We analyzed the binding of the 7C8 antibody to the chloramphenicol phosphonate antigens—one containing a trifluoroacetyl group (CP‐F) and the other containing an acetyl group (CP‐H)—by using isothermal titration calorimetry (ITC). The thermodynamic difference due to the substitution of F by H was evaluated using free energy calculations based on molecular dynamics (MD) simulations. We have previously shown that another antibody, namely, 6D9, binds more weakly to CP‐H than to CP‐F, mainly due to the different hydration free energies of the dissociated state and not due to the unfavorable hydrophobic interactions with the antibody in the bound state. Unlike in the binding of the trifluoroacetyl group with 6D9, in its binding with 7C8, it is exposed to the solvent, as seen in the crystal structure of the complex of 7C8 with CP‐F. The thermodynamic analysis performed in this study showed that the binding affinity of 7C8 for CP‐H is similar to that for CP‐F, but this binding to CP‐H is accompanied with less favorable enthalpy and more favorable entropy changes. The free energy calculations indicated that, upon the substitution of F by H, enthalpy and entropy changes in the associated and dissociated states were decreased, but the magnitude of enthalpy and entropy changes in the dissociated state was larger than that in the associated state. The differences in binding free energy, enthalpy, and entropy changes determined by the free energy calculations for the substitution of F by H are in good agreement with the experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Calculation of cyclodextrin binding affinities: energy, entropy, and implications for drug design

The second generation Mining Minima method yields binding affinities accurate to within 0.8 kcal/mol for the associations of alpha-, beta-, and gamma-cyclodextrin with benzene, resorcinol, flurbiprofen, naproxen, and nabumetone. These calculations require hours to a day on a commodity computer. The calculations also indicate that the changes in configurational entropy upon binding oppose association by as much as 24 kcal/mol and result primarily from a narrowing of energy wells in the bound versus the free state, rather than from a drop in the number of distinct low-energy conformations on binding. Also, the configurational entropy is found to vary substantially among the bound conformations of a given cyclodextrin-guest complex. This result suggests that the configurational entropy must be accounted for to reliably rank docked conformations in both host-guest and ligand-protein complexes. In close analogy with the common experimental observation of entropy-enthalpy compensation, the computed entropy changes show a near-linear relationship with the changes in mean potential plus solvation energy.     

X-ray crystallographic observation of "in-line" and "adjacent" conformations in a bulged self-cleaving RNA/DNA hybrid

The RNA strand in an RNA/DNA duplex with unpaired ribonucleotides can undergo self-cleavage at bulge sites in the presence of a variety of divalent metal ions (Hüsken et al., Biochemistry, 1996, 35:16591-16600). Transesterification proceeds via an in-line mechanism, with the 2'-OH of the bulged nucleotide attacking the 3'-adjacent phosphate group. The site-specificity of the reaction is most likely a consequence of the greater local conformational freedom of the RNA backbone in the bulge region. A standard A-form backbone geometry prohibits formation of an in-line arrangement between 2'-oxygen and phosphate. However, the backbone in the region of an unpaired nucleotide appears to be conducive to an in-line approach. Therefore, the bulge-mediated phosphoryl transfer reaction represents one of the simplest RNA self-cleavage systems. Here we focus on the conformational features of the RNA that underlie site-specific cleavage. The structures of an RNA/DNA duplex with single ribo-adenosyl bulges were analyzed in two crystal forms, permitting observation of 10 individual conformations of the RNA bulge moiety. The bulge geometries cover a range of relative arrangements between the 2'-oxygen of the bulged nucleotide and the P-O5' bond (including adjacent and near in-line) and give a detailed picture of the conformational changes necessary to line up the 2'-OH nucleophile and scissile bond. Although metal ions are of crucial importance in the catalysis of analogous cleavage reactions by ribozymes, it is clear that local strain or conformational flexibility in the RNA also affect cleavage selectivity and rate (Soukup & Breaker, RNA, 1999, 5:1308-1325). The geometries of the RNA bulges frozen out in the crystals provide snapshots along the reaction pathway prior to the transition state of the phosphoryl transfer reaction.     

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.     

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.     

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.     

Equilibrium and kinetic analysis of nucleotide binding to the DEAD-box RNA helicase DbpA

The Escherichia coli DEAD-box protein A (DbpA) is an RNA helicase that utilizes the energy from ATP binding and hydrolysis to facilitate structural rearrangements of rRNA. We have used the fluorescent nucleotide analogues, mantADP and mantATP, to measure the equilibrium binding affinity and kinetic mechanism of nucleotide binding to DbpA in the absence of RNA. Binding generates an enhancement in mant-nucleotide fluorescence and a corresponding reduction in intrinsic DbpA fluorescence, consistent with fluorescence resonance energy transfer (FRET) from DbpA tryptophan(s) to bound nucleotides. Fluorescent modification does not significantly interfere with the affinities and kinetics of nucleotide binding. Different energy transfer efficiencies between DbpA-mantATP and DbpA-mantADP complexes suggest that DbpA adopts nucleotide-dependent conformations. ADP binds (K(d) approximately 50 microM at 22 degrees C) 4-7 times more tightly than ATP (K(d) approximately 400 microM at 22 degrees C). Both nucleotides bind with relatively temperature-independent association rate constants (approximately 1-3 microM(-1) s(-1)) that are much lower than predicted for a diffusion-limited reaction. Differences in the binding affinities are dictated primarily by the dissociation rate constants. ADP binding occurs with a positive change in the heat capacity, presumably reflecting a nucleotide-induced conformational rearrangement of DbpA. At low temperatures (<22 degrees C), the binding free energies are dominated by favorable enthalpic and unfavorable entropic contributions. At physiological temperatures (>22 degrees C), ADP binding occurs with positive entropy changes. We favor a mechanism in which ADP binding increases the conformational flexibility and dynamics of DbpA.     

Conformational Entropy of the RNA Phosphate Backbone and Its Contribution to the Folding Free Energy

While major contributors to the free energy of RNA tertiary structures such as basepairing, base-stacking, and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the main-chain torsion-angle conformation space between every pair of nucleobases along any sequence to compute the corresponding backbone entropy. To validate this method, we have compared the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several-thousand backbone conformations between nearest-neighbor nucleobases as well as against extensive computer simulations. Using this calculation, we analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized along their sequences. The total entropic penalty due to steric congestions in the backbone for the native fold can be as high as 2.4 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 0.72 kcal/mol per nucleotide. For these RNAs, we found that low-entropy high-strain residues are predominantly located at loops with tight turns and at tertiary interaction platforms with unusual structural motifs.     

Conformational Entropy of the RNA Phosphate Backbone and Its Contribution to the Folding Free Energy

While major contributors to the free energy of RNA tertiary structures such as basepairing, base-stacking, and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the main-chain torsion-angle conformation space between every pair of nucleobases along any sequence to compute the corresponding backbone entropy. To validate this method, we have compared the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several-thousand backbone conformations between nearest-neighbor nucleobases as well as against extensive computer simulations. Using this calculation, we analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized along their sequences. The total entropic penalty due to steric congestions in the backbone for the native fold can be as high as 2.4 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 0.72 kcal/mol per nucleotide. For these RNAs, we found that low-entropy high-strain residues are predominantly located at loops with tight turns and at tertiary interaction platforms with unusual structural motifs.     

Evaluating conformational changes in protein structures binding RNA

Many protein-RNA recognition events are known to exhibit conformational changes from qualitative observations of individual complexes. However, a quantitative estimation of conformational changes is required if protein-RNA docking and template-based methods for RNA binding site prediction are to be developed. This study presents the first quantitative evaluation of conformational changes that occur when proteins bind RNA. The analysis of twelve RNA-binding proteins in the bound and unbound states using error-scaled difference distance matrices is presented. The binding site residues are mapped to each structure, and the conformational changes that affect these residues are evaluated. Of the twelve proteins four exhibit greater movements in nonbinding site residues, and a further four show the greatest movements in binding site residues. The remaining four proteins display no significant conformational change. When interface residues are found to be in conformationally variable regions of the protein they are typically seen to move less than 2 A between the bound and unbound conformations. The current data indicate that conformational changes in the binding site residues of RNA binding proteins may not be as significant as previously suggested, but a larger data set is required before wider conclusions may be drawn. The implications of the observed conformational changes for protein function prediction are discussed.     

Modeling disordered regions in proteins using Rosetta

Protein structure prediction methods such as Rosetta search for the lowest energy conformation of the polypeptide chain. However, the experimentally observed native state is at a minimum of the free energy, rather than the energy. The neglect of the missing configurational entropy contribution to the free energy can be partially justified by the assumption that the entropies of alternative folded states, while very much less than unfolded states, are not too different from one another, and hence can be to a first approximation neglected when searching for the lowest free energy state. The shortcomings of current structure prediction methods may be due in part to the breakdown of this assumption. Particularly problematic are proteins with significant disordered regions which do not populate single low energy conformations even in the native state. We describe two approaches within the Rosetta structure modeling methodology for treating such regions. The first does not require advance knowledge of the regions likely to be disordered; instead these are identified by minimizing a simple free energy function used previously to model protein folding landscapes and transition states. In this model, residues can be either completely ordered or completely disordered; they are considered disordered if the gain in entropy outweighs the loss of favorable energetic interactions with the rest of the protein chain. The second approach requires identification in advance of the disordered regions either from sequence alone using for example the DISOPRED server or from experimental data such as NMR chemical shifts. During Rosetta structure prediction calculations the disordered regions make only unfavorable repulsive contributions to the total energy. We find that the second approach has greater practical utility and illustrate this with examples from de novo structure prediction, NMR structure calculation, and comparative modeling.     

A model on the origin of RNA

The theory of the liquid-glass transition is extended to describe the polymerization of RNA in a nucleotide-condensed state. In the glassy state the ribose subunits are joined by a 2'-5' or 3'-5' phosphodiester linkage to form the ribose-phosphate backbone similar to oligosaccharides and polysaccharides. The occurrence of the glass transition requires two conditions: (1) a supercooled state in a nucleotide-condensed state should exist below the temperature at which the whole RNA hydrolyzes; (2) the Gibbs free energy due to the Kauzmann entropy, which obeys a Curie law with a negative sign, must be larger than the height of the potential barrier for nucleotides to overcome to form the binding.     

Bulge loops used to measure the helical twist of RNA in solution.

Bulge loops are commonly found in helical segments of cellular RNAs. When incorporated into long double-stranded RNAs, they may introduce points of flexibility or permanent bend that can be detected by the altered electrophoretic gel mobility of the RNA. We find that a single An or Un bulge loop near the middle of a long RNA helix significantly retards the RNA during polyacrylamide gel electrophoresis if n greater than or equal to 2. The mobility of an RNA containing two A2 bulges various periodically with the number of base pairs between the bulges. We interpret this to mean that A2 bulges varies periodically with the number of base pairs between the bulges. We interpret this to mean that Z2 bulges form torsionally stiff bends in the helix; the gel mobility reaches a minimum when the total helical twist between the bulges rotates the arms of the molecule into a cis conformation. The gel mobilities are proportional to the predicted end-to-end distance of the RNA if the average RNA helical repeat is 11.8 +/- 0.2 bp/turn and there is no helical twist (3 +/- 9 degrees) associated with the bulge (data obtained in 0.15 M Na+). Other sizes and sequences of bulges have very different effects on RNA helix conformation and flexibility. U2 bulges bend the helix to a much smaller degree than A2 bulges, while longer A or U bulge sequences probably allow bends of 90 degrees or more; all of these may be fairly flexible joints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Non-nearest-neighbor dependence of the stability for RNA group II single-nucleotide bulge loops

Thirty-one RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each sequence were determined. The bulge loops were of the group II variety, where the bulged nucleotide is identical to one of its nearest neighbors, leading to ambiguity as to the exact position of the bulge. The data were used to develop a model to predict the free energy of an RNA duplex containing a single-nucleotide bulge. The destabilization of the duplex by the bulge was primarily related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. Since there is an ambiguity of the bulge position for group II bulges, several different stem combinations are possible. The destabilization of group II bulge loops is similar to the destabilization of group I bulge loops, if the second least stable stem is used to predict the influence of the group II bulge. In-line structure probing of the group II bulge loop embedded in a hairpin indicates that the bulged nucleotide is the one positioned farther from the hairpin loop.