Thermodynamic examination of trinucleotide bulged RNA in the context of HIV-1 TAR RNA

RNA structures contain many bulges and loops that are expected to be sites for inter- and intra-molecular interactions. Nucleotides in the bulge are expected to influence the structure and recognition of RNA. The same stability is assigned to all trinucleotide bulged RNA in the current secondary structure prediction models. In this study thermal denaturation experiments were performed on four trinucleotide bulged RNA, in the context of HIV-1 TAR RNA, to determine whether the bulge sequence affects RNA stability and its divalent ion interactions. Cytosine-rich bulged RNA were more stable than uracil-rich bulged RNA in 1 M KCl. Interactions of divalent ions were more favorable with uracil-rich bulged RNA by ~2 kcal/mol over cytosine-rich bulged RNA. The UCU-TAR RNA (wild type) is stabilized by 1.7 kcal/mol in 9.5 mM Ca²⁺ as compared with 1 M KCl, whereas no additional gain in stability is measured for CCC-TAR RNA. These results have implications for base substitution experiments traditionally employed to identify metal ion binding sites. To our knowledge, this is the first systematic study to quantify the effect of small sequence changes on RNA stability upon interactions with divalent ions.     

Structure and folding of a rare, natural kink turn in RNA with an A*A pair at the 2b*2n position

The kink turn (k-turn) is a frequently occurring motif, comprising a bulge followed by G•A and A•G pairs that introduces a sharp axial bend in duplex RNA. Natural k-turn sequences exhibit significant departures from the consensus, including the A•G pairs that form critical interactions stabilizing the core of the structure. Kt-23 found in the small ribosomal subunit differs from the consensus in many organisms, particularly in the second A•G pair distal to the bulge (2b•2n). Analysis of many Kt-23 sequences shows that the frequency of occurrence at the 2n position (i.e., on the nonbulged strand, normally G in standard k-turns) is U>C>G>A. Less than 1% of sequences have A at the 2n position, but one such example occurs in Thelohania solenopsae Kt-23. This sequence folds only weakly in the presence of Mg²⁺ ions but is induced to fold normally by the binding of L7Ae protein. Introduction of this sequence into the SAM-I riboswitch resulted in normal binding of SAM ligand, indicating that tertiary RNA contacts have resulted in k-turn folding. X-ray crystallography shows that the T. solenopsae Kt-23 adopts a standard k-turn geometry, making the key, conserved hydrogen bonds in the core and orienting the 1n (of the bulge-proximal A•G pair) and 2b adenine nucleobases in position facing the opposing minor groove. The 2b and 2n adenine nucleobases are not directly hydrogen bonded, but each makes hydrogen bonds to their opposing strands.     

Thermodynamic examination of 1- to 5-nt purine bulge loops in RNA and DNA constructs

Bulge loops are common features of RNA structures that are involved in the formation of RNA tertiary structures and are often sites for interactions with proteins and ions. Minimal thermodynamic data currently exist on the bulge size and sequence effects. Using thermal denaturation methods, thermodynamic properties of 1- to 5-nt adenine and guanine bulge loop constructs were examined in 10 mM MgCl₂ or 1 M KCl. The ΔG37∘ loop parameters for 1- to 5-nt purine bulge loops in RNA constructs were between 3.07 and 5.31 kcal/mol in 1 M KCl buffer. In 10 mM magnesium ions, the ΔΔG° values relative to 1 M KCl were 0.47–2.06 kcal/mol more favorable for the RNA bulge loops. The ΔG37∘ loop parameters for 1- to 5-nt purine bulge loops in DNA constructs were between 4.54 and 5.89 kcal/mol. Only 4- and 5-nt guanine constructs showed significant change in stability for the DNA constructs in magnesium ions. A linear correlation is seen between the size of the bulge loop and its stability. New prediction models are proposed for 1- to 5-nt purine bulge loops in RNA and DNA in 1 M KCl. We show that a significant stabilization is seen for small bulge loops in RNA in the presence of magnesium ions. A prediction model is also proposed for 1- to 5-nt purine bulge loop RNA constructs in 10 mM magnesium chloride.     

Enhancement of the rate of pyrophosphate hydrolysis by nonenzymatic catalysts and by inorganic pyrophosphatase

To estimate the proficiency of inorganic pyrophosphatase as a catalyst, ³¹P NMR was used to determine rate constants and thermodynamics of activation for the spontaneous hydrolysis of inorganic pyrophosphate (PP_i) in the presence and absence of Mg²⁺ at elevated temperatures. These values were compared with rate constants and activation parameters determined for the reaction catalyzed by Escherichia coli inorganic pyrophosphatase using isothermal titration calorimetry. At 25 °C and pH 8.5, the hydrolysis of MgPP_i²⁻ proceeds with a rate constant of 2.8 × 10⁻¹⁰ s⁻¹, whereas E. coli pyrophosphatase was found to have a turnover number of 570 s⁻¹ under the same conditions. The resulting rate enhancement (2 × 10¹²-fold) is achieved entirely by reducing the enthalpy of activation (ΔΔH^‡ = −16.6 kcal/mol). The presence of Mg²⁺ ions or the transfer of the substrate from bulk water to dimethyl sulfoxide was found to increase the rate of pyrophosphate hydrolysis by as much as ∼10⁶-fold. Transfer to dimethyl sulfoxide accelerated PP_i hydrolysis by reducing the enthalpy of activation. Mg²⁺ increased the rate of PP_i hydrolysis by both increasing the entropy of activation and reducing the enthalpy of activation.     

Specificity of RSG-1.2 peptide binding to RRE-IIB RNA element of HIV-1 over Rev peptide is mainly enthalpic in origin

Rev is an essential HIV-1 regulatory protein which binds to the Rev responsive element (RRE) present within the env gene of HIV-1 RNA genome. This binding facilitates the transport of the RNA to the cytoplasm, which in turn triggers the switch between viral latency and active viral replication. Essential components of this complex have been localized to a minimal arginine rich Rev peptide and stem IIB region of RRE. A synthetic peptide known as RSG-1.2 binds with high binding affinity and specificity to the RRE-IIB than the Rev peptide, however the thermodynamic basis of this specificity has not yet been addressed. The present study aims to probe the thermodynamic origin of this specificity of RSG-1.2 over Rev Peptide for RRE-IIB. The temperature dependent melting studies show that RSG-1.2 binding stabilizes the RRE structure significantly (ΔT _m = 4.3°C), in contrast to Rev binding. Interestingly the thermodynamic signatures of the binding have also been found to be different for both the peptides. At pH 7.5, RSG-1.2 binds RRE-IIB with a K_a = 16.2±0.6×10⁷ M⁻¹ where enthalpic change ΔH = −13.9±0.1 kcal/mol is the main driving force with limited unfavorable contribution from entropic change TΔS = −2.8±0.1 kcal/mol. A large part of ΔH may be due to specific stacking between U72 and Arg15. In contrast binding of Rev (K_a = 3.1±0.4×10⁷ M⁻¹) is driven mainly by entropy (ΔH = 0 kcal/mol and TΔS = 10.2±0.2 kcal/mol) which arises from major conformational changes in the RNA upon binding.     

Thermodynamic analysis of ligand binding and ligand binding-induced tertiary structure formation by the thiamine pyrophosphate riboswitch

The thi-box riboswitch regulates gene expression in response to the intracellular concentration of thiamine pyrophosphate (TPP) in archaea, bacteria, and eukarya. To complement previous biochemical, genetic, and structural studies of this phylogenetically widespread RNA domain, we have characterized its interaction with TPP by isothermal titration calorimetry. This shows that TPP binding is highly dependent on Mg²⁺ concentration. The dissociation constant decreases from ∼200 nM at 0.5 mM Mg²⁺ concentration to ∼9 nM at 2.5 mM Mg²⁺ concentration. Binding is enthalpically driven, but the unfavorable entropy of binding decreases as Mg²⁺ concentration rises, suggesting that divalent cations serve to pre-organize the RNA. Mutagenesis, biochemical analysis, and a new crystal structure of the riboswitch suggest that a critical element that participates in organizing the riboswitch structure is the tertiary interaction formed between the P3 and L5 regions. This tertiary contact is distant from the TPP binding site, but calorimetric analysis reveals that even subtle mutations in L5 can have readily detectable effects on TPP binding. The thermodynamic signatures of these mutations, namely decreased favorable enthalpy of binding and small effects on entropy of binding, are consistent with the P3–L5 association contributing allosterically to TPP-induced compaction of the RNA.     

Thermodynamic contributions of single internal rA·dA, rC·dC, rG·dG and rU·dT mismatches in RNA/DNA duplexes

The thermodynamic contributions of rA·dA, rC·dC, rG·dG and rU·dT single internal mismatches were measured for 54 RNA/DNA duplexes in a 1 M NaCl buffer using UV absorbance thermal denaturation. Thermodynamic parameters were obtained by fitting absorbance versus temperature profiles using the curve-fitting program Meltwin. The weighted average thermodynamic data were fit using singular value decomposition to determine the eight non-unique nearest-neighbor parameters for each internal mismatch. The new parameters predict the ΔG°₃₇, ΔH° and melting temperature (T_m) of duplexes containing these single mismatches within an average of 0.33 kcal/mol, 4.5 kcal/mol and 1.4°C, respectively. The general trend in decreasing stability for the single internal mismatches is rG·dG > rU·dT > rA·dA > rC·dC. The stability trend for the base pairs 5′ of the single internal mismatch is rG·dC > rC·dG > rA·dT > rU·dA. The stability trend for the base pairs 3′ of the single internal mismatch is rC·dG > rG·dC >> rA·dT > rU·dA. These nearest-neighbor values are now a part of a complete set of single internal mismatch thermodynamic parameters for RNA/DNA duplexes that are incorporated into the nucleic acid assay development software programs Visual oligonucleotide modeling platform (OMP) and ThermoBLAST.     

Ion-induced folding of a kink turn that departs from the conventional sequence

Kink turns (k-turns) are important structural motifs that create a sharp axial bend in RNA. Most conform to a consensus in which a three-nucleotide bulge is followed by consecutive G•A and A•G base pairs, and when these G•A pairs are modified in vitro this generally leads to a failure to adopt the k-turn conformation. Kt-23 in the 30S ribosomal subunit of Thermus thermophilus is a rare exception in which the bulge-distal A•G pair is replaced by a non-Watson–Crick A•U pair. In the context of the ribosome, Kt-23 adopts a completely conventional k-turn geometry. We show here that this sequence is induced to fold into a k-turn structure in an isolated RNA duplex by Mg²⁺ or Na⁺ ions. Therefore, the Kt-23 is intrinsically stable despite lacking the key A•G pair; its formation requires neither tertiary interactions nor protein binding. Moreover, the Kt-23 k-turn is stabilized by the same critical hydrogen-bonding interactions within the core of the structure that are found in more conventional sequences such as the near-consensus Kt-7. T. thermophilus Kt-23 has two further non-Watson–Crick base pairs within the non-canonical helix, three and four nucleotides from the bulge, and we find that the nature of these pairs influences the ability of the RNA to adopt k-turn conformation, although the base pair adjacent to the A•U pair is more important than the other.     

Effect of 6-thioguanine on the stability of duplex DNA

The incorporation of 6-thioguanine (S6G) into DNA is a prerequisite for its cytotoxic action, but duplex structure is not significantly perturbed by the presence of the lesion [J. Bohon and C. R. de los Santos (2003) Nucleic Acids Res., 31, 1331–1338]. It is therefore possible that the mechanism of cytotoxicity relies on a loss of stability rather than a pathway involving direct structural recognition. The research described here focuses on the changes in thermodynamic properties of duplex DNA owing to the introduction of S6G as well as the kinetic properties of base pairs involving S6G. Replacement of a guanine in a G•C pair by S6G results in ~1 kcal/mol less favorable Gibbs free energy of duplex formation at 37°C. S6G•T and G•T mismatch-containing duplexes have almost identical Gibbs free energy at 37°C, with values ~3 kcal/mol less favorable than that of the control. Base pair stability is affected by S6G. The lifetime of the normal G•C base pair is ~125 ms, whereas that of the G•T mismatch is below the detection limit. The lifetimes of S6G•C and S6G•T pairs are ~7 and 2 ms, respectively, demonstrating that, although S6G significantly decreases the stability of the pairing with cytosine, it slightly increases that of a mismatch.     

Twist-joints and double twist-joints in RNA structure

Analysis of available RNA crystal structures has allowed us to identify a new family of RNA arrangements that we call double twist-joints, or DTJs. Each DTJ is composed of a double helix that contains two bulges incorporated into different strands and separated from each other by 2 or 3 bp. At each bulge, the double helix is over-twisted, while the unpaired nucleotides of both bulges form a complex network of stacking and hydrogen-bonding with nucleotides of helical regions. In total, we identified 14 DTJ cases, which can be combined in three groups based on common structural characteristics. One DTJ is found in a functional center of the ribosome, another DTJ mediates binding of the pre-tRNA to the RNase P, and two more DTJs form the sensing domains in the glycine riboswitch.     

Mg2+-induced conformational changes in the btuB riboswitch from E. coli

Mg²⁺ has been shown to modulate the function of riboswitches by facilitating the ligand-riboswitch interactions. The btuB riboswitch from Escherichia coli undergoes a conformational change upon binding to its ligand, coenzyme B₁₂ (adenosyl-cobalamine, AdoCbl), and down-regulates the expression of the B₁₂ transporter protein BtuB in order to control the cellular levels of AdoCbl. Here, we discuss the structural folding attained by the btuB riboswitch from E. coli in response to Mg²⁺ and how it affects the ligand binding competent conformation of the RNA. The btuB riboswitch notably adopts different conformational states depending upon the concentration of Mg²⁺. With the help of in-line probing, we show the existence of at least two specific conformations, one being achieved in the complete absence of Mg²⁺ (or low Mg²⁺ concentration) and the other appearing above ∼0.5 mM Mg²⁺. Distinct regions of the riboswitch exhibit different dissociation constants toward Mg²⁺, indicating a stepwise folding of the btuB RNA. Increasing the Mg²⁺ concentration drives the transition from one conformation toward the other. The conformational state existing above 0.5 mM Mg²⁺ defines the binding competent conformation of the btuB riboswitch which can productively interact with the ligand, coenzyme B₁₂, and switch the RNA conformation. Moreover, raising the Mg²⁺ concentration enhances the ratio of switched RNA in the presence of AdoCbl. The lack of a AdoCbl-induced conformational switch experienced by the btuB riboswitch in the absence of Mg²⁺ indicates a crucial role played by Mg²⁺ for defining an active conformation of the riboswitch.     

RNA recognition by Tat-derived peptides: interaction in the major groove?

Replication of human immunodeficiency virus requires binding of the viral Tat protein to its RNA target sequence TAR; peptides derived from Tat bind to a TAR "contact site" spanning 5 bp and a trinucleotide pyrimidine bulge. We find that high affinity binding requires a U residue in the bulge loop and 2 specific adjacent base pairs. Other bulged RNAs bind in a lower affinity nonspecific manner; sequence-specific binding requires a bulge loop of more than 1 nucleotide. Reaction with diethyl pyrocarbonate indicates that one effect of the bulge is to make the otherwise deep and narrow RNA major groove accessible. A model consistent with these data involves local distortion of A-form geometry at the bulge, which bends the helix and permits protein binding and interactive access in the RNA major groove.     

A conserved RNA pseudoknot in a putative molecular switch domain of the 3'-untranslated region of coronaviruses is only marginally stable

The 3'-untranslated region (UTR) of the group 2 coronavirus mouse hepatitis virus (MHV) genome contains a predicted bulged stem-loop (designated P0ab), a conserved cis-acting pseudoknot (PK), and a more distal stem-loop (designated P2). Base-pairing to create the pseudoknot-forming stem (P1(pk)) is mutually exclusive with formation of stem P0a at the base of the bulged stem-loop; as a result, the two structures cannot be present simultaneously. Herein, we use thermodynamic methods to evaluate the ability of individual subdomains of the 3' UTR to adopt a pseudoknotted conformation. We find that an RNA capable of forming only the predicted PK (58 nt; 3' nucleotides 241-185) adopts the P2 stem-loop with little evidence for P1(pk) pairing in 0.1 M KCl and the absence of Mg(2+); as Mg(2+) or 1 M KCl is added, a new thermal unfolding transition is induced and assignable to P1(pk) pairing. The P1(pk) helix is only marginally stable, ΔG(25) ≈ 1.2 ± 0.3 kcal/mol (5.0 mM Mg(2+), 100 mM K(+)), and unfolded at 37°C. Similar findings characterize an RNA 5' extended through the P0b helix only (89 nt; 294-185). In contrast, an RNA capable of forming either the P0a helix or the pseudoknot (97 nt; 301-185) forms no P1(pk) helix. Thermal unfolding simulations are fully consistent with these experimental findings. These data reveal that the PK forms weakly and only when the competing double-hairpin structure cannot form; in the UTR RNA, the double hairpin is the predominant conformer under all solution conditions.     

Crystal structures of the thi-box riboswitch bound to thiamine pyrophosphate analogs reveal adaptive RNA-small molecule recognition

Riboswitches are noncoding mRNA elements that bind small-molecule metabolites with high affinity and specificity, and they regulate the expression of associated genes. The thi-box riboswitch can exhibit a 1000-fold higher affinity for thiamine pyrophosphate over closely related noncognate compounds such as thiamine monophosphate. To understand the chemical basis of thi-box pyrophosphate specificity, we have determined crystal structures of an E. coli thi-box bound to thiamine pyrophosphate, thiamine monophosphate, and the structural analogs benfotiamine and pyrithiamine. When bound to monophosphorylated compounds, the RNA elements that recognize the thiamine and phosphate moieties of the ligand move closer together. This allows the riboswitch to recognize the monophosphate in a manner similar to how it recognizes the beta-phosphate of thiamine pyrophosphate. In the pyrithiamine complex, the pyrophosphate binding site is largely unstructured. These results show how the riboswitch can bind to various metabolites, and why the thi-box preferentially binds thiamine pyrophosphate.     

A difference in the importance of bulged nucleotides and their parent base pairs in the binding of transcription factor IIIA to Xenopus 5S RNA and 5S RNA genes

Individual bulge loops present in Xenopus 5S RNA (positions 49A-A50 in helix III, C63 in helix II and A83 in helix IV), were deleted by site directed mutagenesis. The interaction of these mutant 5S RNA molecules with TFIIIA was measured by a direct binding assay and a competition assay. The results of these experiments show that none of the bulged nucleotides in Xenopus 5S RNA are required for the binding of TFIIIA. The affinity of the mutant 5S RNA genes for TFIIIA was also studied by a filter binding assay. In contrast to the effect that deleting bulged nucleotides had on the TFIIIA-RNA binding affinity, deletion of the corresponding A-T base pair at position +83 in 5S DNA was found to reduce the apparent association constant of TFIIIA by a factor of four-fold.     

Thermodynamic and kinetic characterization of ligand binding to the purine riboswitch aptamer domain

Riboswitches are cis-acting genetic regulatory elements found commonly in bacterial mRNAs that consist of a metabolite-responsive aptamer domain coupled to a regulatory switch. Purine riboswitches respond to intracellular concentrations of either adenine or guanine/hypoxanthine to control gene expression. The aptamer domain of the purine riboswitch contains a pyrimidine residue (Y74) that forms a Watson-Crick base-pairing interaction with the bound purine nucleobase ligand that discriminates between adenine and guanine. We sought to understand the structural basis of this specificity and the mechanism of ligand recognition by the purine riboswitch. Here, we present the 2,6-diaminopurine-bound structure of a C74U mutant of the xpt-pbuX guanine riboswitch, along with a detailed thermodynamic and kinetic analysis of nucleobase recognition by both the native and mutant riboswitches. These studies demonstrate clearly that the pyrimidine at position 74 is the sole determinant of purine riboswitch specificity. In addition, the mutant riboswitch binds adenine and adenine derivatives well compared with the guanine-responsive riboswitch. Under our experimental conditions, 2,6-diaminopurine binds the RNA with DeltaH=-40.3 kcal mol(-1), DeltaS=-97.6 cal mol(-1)K(-1), and DeltaG=-10.73 kcal mol(-1). A kinetic determination of the slow rate (0.15 x 10(5)M(-1)s(-1) and 2.1 x 10(5)mM(-1)s(-1) for 2-aminopurine binding the adenine-responsive mutant riboswitch and 7-deazaguanine-binding guanine riboswitch, respectively) of association under varying experimental conditions allowed us to propose a mechanism for ligand recognition by the purine riboswitch. A conformationally dynamic unliganded state for the binding pocket is stabilized first by the Watson-Crick base pairing between the ligand and Y74, and by the subsequent ordering of the J2/3 loop, enclosing the ligand within the three-way junction.     

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was 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. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.     

Day-night changes of energy-rich compounds in crassulacean acid metabolism (CAM) species utilizing hexose and starch

• Background and Aims Plants with crassulacean acid metabolism (CAM) can be divided into two groups according to the major carbohydrates used for malic acid synthesis, either polysaccharide (starch) or monosaccharide (hexose). This is related to the mechanism and affects energy metabolism in the two groups. In Kalanchoë pinnata and K. daigremontiana, which utilize starch, ATP-dependent phosphofructokinase (tonoplast inorganic pyrophosphatase) activity is greater than inorganic pyrophosphate-dependent phosphofructokinase (tonoplast adenosine triphosphatase) activity, but the reverse is the case in pineapple (Ananas comosus) utilizing hexose. To test the hypothesis that the energy metabolism of the two groups differs, day-night changes in the contents of ATP, ADP, AMP, inorganic phosphate (Pi), phosphoenolpyruvate (PEP) and inorganic pyrophosphate (PPi) in K. pinnata and K. daigremontiana leaves and in pineapple chlorenchyma were analysed.• Methods The contents of energy-rich compounds were measured spectrophotometrically in extracts of tissue sampled in the light and dark, using potted plants, kept for 15 d before the experiments in a growth chamber.• Key Results In the three species, ATP content and adenylate energy charge (AEC) increased in the dark and decreased in the light, in contrast to ADP and AMP. Changes in ATP and AEC were greater in Kalanchoë leaves than in pineapple chlorenchyma. PPi content in the three species increased in the dark, but on illumination it decreased rapidly and substantially, remaining little changed through the rest of the light period. Pi content of Kalanchoë leaves did not change between dark and light, whereas Pi in pineapple chlorenchyma increased in the dark and decreased in the light, and the changes were far greater than in Kalanchoë leaves. Light-dark changes in PEP content in the three species were similar.• Conclusions These results corroborate our hypothesis that day–night changes in the contents of energy-rich compounds differ between CAM species and are related to the carbohydrate used for malic acid synthesis.Key words: Ananas comosus, ATP, chlorenchyma, crassulacean acid metabolism, inorganic pyrophosphate, Kalanchoë daigremontiana, Kalanchoë pinnata, phosphoenolpyruvate     

Predicting loop–helix tertiary structural contacts in RNA pseudoknots

Tertiary interactions between loops and helical stems play critical roles in the biological function of many RNA pseudoknots. However, quantitative predictions for RNA tertiary interactions remain elusive. Here we report a statistical mechanical model for the prediction of noncanonical loop–stem base-pairing interactions in RNA pseudoknots. Central to the model is the evaluation of the conformational entropy for the pseudoknotted folds with defined loop–stem tertiary structural contacts. We develop an RNA virtual bond-based conformational model (Vfold model), which permits a rigorous computation of the conformational entropy for a given fold that contains loop–stem tertiary contacts. With the entropy parameters predicted from the Vfold model and the energy parameters for the tertiary contacts as inserted parameters, we can then predict the RNA folding thermodynamics, from which we can extract the tertiary contact thermodynamic parameters from theory–experimental comparisons. These comparisons reveal a contact enthalpy (ΔH) of −14 kcal/mol and a contact entropy (ΔS) of −38 cal/mol/K for a protonated C⁺•(G–C) base triple at pH 7.0, and (ΔH = −7 kcal/mol, ΔS = −19 cal/mol/K) for an unprotonated base triple. Tests of the model for a series of pseudoknots show good theory–experiment agreement. Based on the extracted energy parameters for the tertiary structural contacts, the model enables predictions for the structure, stability, and folding pathways for RNA pseudoknots with known or postulated loop–stem tertiary contacts from the nucleotide sequence alone.