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.     

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

Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson–Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson–Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson–Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3′ or 5′ of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.     

Stability of single-nucleotide bulge loops embedded in a GAAA RNA hairpin stem

Forty-six RNA hairpins containing combinations of 3' or 5' bulge loops and a 3' or 5' fluorescein label were optically melted in 1 M NaCl, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each hairpin were determined. The bulge loops were of the group I variety, in which the identity of the bulge is known, and the group II variety, in which the bulged nucleotide is identical to one of its nearest neighbors, leading to ambiguity as to the exact position of the bulge. The fluorescein label at either the 3' end or 5' end of the hairpin did not significantly influence the stability of the hairpin. As observed with bulge loops inserted into a duplex motif, the insertion of a bulge loop into the stem of a hairpin loop was destabilizing. The model developed to predict the influence of bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability. Specifically, the influence of the bulge is related to the stability of the hairpin stem distal from the hairpin loop.     

Analysis of the stability and flexibility of RNA complexes containing bulge loops of different sizes

Molecular dynamics simulations of RNA molecules consisting of an antisense oligonucleotide forming a complex with a target strand thereby creating an internal bulge-loop with 3, 4, or 5 nucleotides have been performed with and without O2' methylation of the antisense strand. The methylation influcences minor groove hydration, in particular near guanines but also around the methylated O2', and it also reduces the flexibility of both RNA strands. A G.U wobble pair adjacent to the bulge-loop is also found to increase the flexibility of the bulge nucleotides, compared to the situation with a standard Watson-Crick G-C base-pair in the same position.     

Sequence dependence of the stability of RNA hairpin molecules with six nucleotide loops

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequence of the types GCGXUAAUYCGC and GGUXUAAUYACC with Watson-Crick loop closure, where XY is the set of 10 possible mismatch base pairs. A nearest-neighbor analysis of the data indicates the free energy of loop formation at 37 degrees C varies from 3.1 to 5.1 kcal/mol. These results agree with the model previously developed [Vecenie, C. J., and Serra, M. J. (2004) Biochemistry 43, 11813] to predict the stability of RNA hairpin loops: DeltaG degrees (37L(n) = DeltaG degrees (37i(n) + DeltaG degrees (37MM) - 0.8 (if first mismatch is GA or UU) - 0.8 (if first mismatch is GG and loop is closed on the 5' side by a purine). Here, DeltaG degrees (37i(n) is the free energy for initiating a loop of n nucleotides, and DeltaG degrees (37MM) is the free energy for the interaction of the first mismatch with the closing base pair. Thermodynamic parameters are also reported for hairpin formation in 1 M NaCl by RNA sequence of the types GACGXUAAUYUGUC and GGUXUAAUYGCC with GU base pair closure, where XY is the set of 10 possible mismatch base pairs. A nearest-neighbor analysis of the data indicates the free energy of loop formation at 37 degrees C varies from 3.6 to 5.3 kcal/mol. These results allow the development of a model for predicting the stability of hairpin loops closed by GU base pairs. DeltaG degrees (37L(n) (kcal/mol) = DeltaG degrees (37i(n) - 0.8 (if the first mismatch is GA) - 0.8 (if the first mismatch is GG and the loop is closed on the 5' side by a purine). Note that for these hairpins, the stability of the loops does not depend on DeltaG degrees (37MM). For hairpin loops closed by GU base pairs, the DeltaG degrees (37i(n) values, when n = 4, 5, 6, 7, and 8, are 4.9, 5.0, 4.6, 5.0, and 4.8 kcal/mol, respectively. The model gives good agreement when tested against six naturally occurring hairpin sequences. Thermodynamic values for terminal mismatches adjacent to GC, GU, and UG base pairs are also reported.     

Determination of the secondary structure of group II bulge loops using the fluorescent probe 2-aminopurine

Eleven RNA hairpins containing 2-aminopurine (2-AP) in either base-paired or single nucleotide bulge loop positions were optically melted in 1 M NaCl; and, the thermodynamic parameters ΔH°, ΔS°, ΔG°₃₇, and T_M for each hairpin were determined. Substitution of 2-AP for an A (adenosine) at a bulge position (where either the 2-AP or A is the bulge) in the stem of a hairpin, does not affect the stability of the hairpin. For group II bulge loops such as AA/U, where there is ambiguity as to which of the A residues is paired with the U, hairpins with 2-AP substituted for either the 5′ or 3′ position in the hairpin stem have similar stability. Fluorescent melts were performed to monitor the environment of the 2-AP. When the 2-AP was located distal to the hairpin loop on either the 5′ or 3′ side of the hairpin stem, the change in fluorescent intensity upon heating was indicative of an unpaired nucleotide. A database of phylogenetically determined RNA secondary structures was examined to explore the presence of naturally occurring bulge loops embedded within a hairpin stem. The distribution of bulge loops is discussed and related to the stability of hairpin structures.     

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.     

The stability of RNA hairpin loops containing A-U-G: An-U-G-Um

RNA oligomers with the sequence A_n-U-G-U_m, n = 7–9, m = 5–10, have been synthesized and found to form hairpin loops in 21 mM or 101 mM Na⁺. The hairpin loops displayed melting temperatures 13°–29°C greater than that of the hairpin loop A₆-C₈-U₆ in the same solvent. The increased stability of these hairpin loops is attributed to the presence of the trinucleotide A-U-G in the loop. Circular dichroism (CD) spectra were taken of the hairpin loops A₇-U-G-U_6,7,8, A₈-U-G-U_6,7,8, and A₉-U-G-U₆ in 21 mM Na⁺, and compared with circular dichroism spectra of A₆-U₆ in 1 M Na⁺. Difference spectra were calculated between each A_n-U-G-U_m and 11.5 mM (nucleotide) A₆-U₆ at similar temperatures and identical singlestrand fractions to give the “experimental” CD spectrum of the unbase-paired nucleotides in the loop, assuming, five, six, or seven base pairs. CD spectra were calculated for each of the assumed unbase-paired sequences using the measured CD spectra of ApA, A-U-G-, A, and U, and compared with the experimental spectra. The best agreement was found for hairpin-loop models containing five base pairs and five to eight unbase-paired nucleotides in the loop.     

Continuum solvent studies of the stability of RNA hairpin loops and helices

We apply continuum solvent models to investigate the relative stability of various conformational forms for two RNA sequences, GGAC(UUCG)GUCC and GGUG(UGAA)CACC. In the first part, we compare alternate hairpin conformations to explore the reliability of these models to discriminate between different local conformations. A second part looks at the hairpin-duplex conversion for the UUCG sequence, identifying major contributors to the thermodynamics of a much large scale transition. Structures were taken as snapshots from multi-nanosecond molecular dynamics simulations computed in a consistent fashion using explicit solvent and with long-range electrostatics accounted for using the Particle-Mesh Ewald procedure. The electrostatic contribution to solvation energies were computed using both a finite-difference Poisson-Boltzmann (PB) model and a pairwise Generalized Born model; non-electrostatic contributions were estimated with a surface-area dependent term. To these solvation free energies were added the mean solute internal energies (determined from a molecular mechanics potential) and estimates of the solute entropy (from a harmonic analysis). Consistent with experiment and with earlier solvated molecular dynamics simulations, the UUCG hairpin was found to prefer conformers close to a recent NMR structure determination in preference to those from an earlier NMR study. Similarly, results for the UGAA hairpin favored an NMR-derived structure over that to be expected for a generic GNRA hairpin loop. Experimental free energies are not known for the hairpin/duplex conversion, but must be close to zero since hairpins are seen in solution and duplexes in crystals; out calculations find a value near zero and illustrate the expected interplay of solvation, salt effects and entropy in affecting this equilibrium.     

Factors affecting the thermodynamic stability of small asymmetric internal loops in RNA

Optical melting experiments were used to determine the thermodynamic parameters for oligoribonucleotides containing small asymmetric internal loops. The results show a broad range of thermodynamic stabilities, which depend on loop size, asymmetry, sequence, closing base pairs, and length of helix stems. Imino proton NMR experiments provide evidence for possible hydrogen bonding in GA and UU mismatches in some asymmetric loops. The stabilizing effects of GA, GG, and UU mismatches on the thermodynamic stability of internal loops vary depending on the size and asymmetry of the loop. The dependence of loop stability on Watson-Crick closing base pairs may be explained by an account of hydrogen bonds. Models are presented for approximating the free energy increments of 2 x 3 and 1 x 3 internal loops.     

Thermal stability of RNA hairpins containing a four-membered loop and a bulge nucleotide

Fourteen RNA hairpins containing a four-membered loop and a bulge nucleotide were synthesized and their thermal stabilities determined. The combined contribution of a four-membered loop and bulge A to the free energy of a hairpin is calculated to be 9.3 kcal/mol at 37 degrees C and successfully predicts the stability of an independent RNA hairpin. The introduction of a bulge nucleotide to the helical stem of an RNA hairpin destabilizes the molecule in a sequence-dependent manner. The individual thermodynamic contributions of a four-membered loop and bulge A, G, and U residues to the stability of an RNA hairpin loop are presented.     

Experimentally derived nearest-neighbor parameters for the stability of RNA three- and four-way multibranch loops.

Algorithms for predicting RNA secondary structure require approximations for the free energies of multibranch loops, also called junctions. The stabilities of 62 RNA duplexes with three- and four-way multibranch loops were determined by optical melting. To account for the observed sequence dependence, a revised loop free-energy approximation is proposed that accounts for the strain in three-way junctions with fewer than two unpaired nucleotides, penalizes asymmetry in the distribution of unpaired nucleotides, and gives a bonus for four-way loops relative to three-way loops. Parameters for this equation were determined by linear regression.     

The complete nucleotide sequence of the group II RNA coliphage GA

The complete nucleotide sequence of the RNA coliphage GA, a group II phage, is presented. The entire genome comprises 3466 bases. Three large open reading frames were identified, which correspond to the maturation protein gene (390 amino acids), the coat protein gene (129 amino acids) and the replicase beta-subunit protein gene (531 amino acids). In addition, untranslated regions occur at the 5' (135 bases) and 3' (122 bases) ends of the molecule. Two intercistronic untranslated regions occur between the cistrons for the maturation and coat proteins, and between the coat and beta-subunit proteins. We have compared the nucleotide sequence of GA RNA with the published sequence of MS2 RNA, and show that they are related. The comparative structures of two important regulatory regions are presented; the coat protein binding site which is involved in translational repression of the replicase beta-subunit protein gene, and a hairpin in a region proximal to the lysis protein gene.     

Sample size determination for classifiers based on single-nucleotide polymorphisms

Single-nucleotide polymorphisms (SNPs), believed to determine human differences, are widely used to predict risk of diseases. Typically, clinical samples are limited and/or the sampling cost is high. Thus, it is essential to determine an adequate sample size needed to build a classifier based on SNPs. Such a classifier would facilitate correct classifications, while keeping the sample size to a minimum, thereby making the studies cost-effective. For coded SNP data from 2 classes, an optimal classifier and an approximation to its probability of correct classification (PCC) are derived. A linear classifier is constructed and an approximation to its PCC is also derived. These approximations are validated through a variety of Monte Carlo simulations. A sample size determination algorithm based on the criterion, which ensures that the difference between the 2 approximate PCCs is below a threshold, is given and its effectiveness is illustrated via simulations. For the HapMap data on Chinese and Japanese populations, a linear classifier is built using 51 independent SNPs, and the required total sample sizes are determined using our algorithm, as the threshold varies. For example, when the threshold value is 0.05, our algorithm determines a total sample size of 166 (83 for Chinese and 83 for Japanese) that satisfies the criterion.     

Analysis of the complete nucleotide sequence of the group IV RNA coliphage SP.

We report the nucleotide sequence of the Group IV RNA bacteriophage SP. The entire sequence is 4276 nucleotides long. Four cistrons have been identified by comparison with the related Group III phage Q beta. The maturation protein contains 449 amino acids, the coat protein contains 131 amino acids, the read-through protein contains 330 amino acids and the replicase beta-subunit contains 575 amino acids. SP is 59 nucleotides longer than Q beta. We have analyzed both sequence and structural conservation between SP and Q beta and shown that the sequences for the coat and central region of the replicase are strongly conserved between the two genomes. We also show that the S and M replicase binding sites of Q beta are strongly conserved in SP. Interestingly, the base composition of SP and Q beta differ significantly from one another, and most of the differences can be accounted for by a strong preponderance of U in the third position of each codon of Q beta relative to SP. We also compare conserved hairpins associated with potential coat protein and replicase binding sites.     

Complete nucleotide sequence of the group I RNA bacteriophage fr

We report the complete nucleotide sequence of the group I RNA bacteriophage fr. The entire genome consists of 3575 nucleotides, six nucleotides more than the only other sequenced group I representative, MS2. The greatest divergence between these phages occurs in the 5' terminal region of the A gene, while the lysis-replicase gene overlap, the coat gene and the central region of the replicase gene are highly conserved. Overall sequence homology between fr and MS2 is 77%. Here, we present a general comparison between the two phages. In the accompanying paper we use phylogenetic sequence comparison between MS2 and fr to deduce the secondary structure at the 3' untranslated region.     

Influence of the 2'-hydroxyl group conformation on the stability of A-form helices in RNA

The 2'-hydroxyl group plays fundamental roles in both the structure and the function of RNA, and is the major determinant of the conformational and thermodynamic differences between RNA and DNA. Here, we report a conformational analysis of 2'-OH groups of the HIV-2 TAR RNA by means of NMR scalar coupling measurements in solution. Our analysis supports the existence of a network of water molecules spanning the minor groove of an RNA A-form helix, as has been suggested on the basis of a high-resolution X-ray study of an RNA duplex. The 2'-OH protons of the lower stem nucleotides of the TAR RNA project either towards the O3' or towards the base, where the 2'-OH group can favorably participate in H-bonding interactions with a water molecule situated in the nucleotide base plane. We observe that the k(ex) rate of the 2'-OH proton with the bulk solvent anti-correlates with the base-pair stability, confirming the involvement of the 2'-OH group in a collective network of H-bonds, which requires the presence of canonical helical secondary structure. The methodology and conformational analysis presented here are broadly applicable and facilitate future studies aimed to correlate the conformation of the 2'-OH group with both the structure and the function of RNA and RNA-ligand complexes.     

A new set of peptide-based group heat capacities for use in protein stability calculations.

The partial molar heat capacities of the tripeptides of the sequence glycyl-X-glycine, where X is one of the amino acids leucine, threonine, glutamine, phenylalanine, histidine, cysteine, proline, glutamic acid or arginine, and of the two tetrapeptides tetraglycine and glycyltryptophanylglycylglycine in aqueous solution over the temperature range 10-100 degrees C have been determined using high sensitivity scanning microcalorimetry. These results were used to derive the partial molar heat capacities of the various amino acid side-chains. This report completes our programme to derive reliable side-chain heat capacities for all 20 amino acids of proteins over a wide temperature range using the tripeptides Gly-X-Gly as realistic model compounds. Included in the study is a summary of the partial molar heat capacities of all 20 amino acid side-chains. These results, along with the heat capacity of the peptide backbone group, were used to calculate the partial molar heat capacities of some oligopeptides and of the random coil form of some unfolded proteins in water. The calculated heat capacities of the proteins obtained using this new set of heat capacities for the constituent groups are consistent with the heat capacities of the denatured state determined experimentally.     

Escherichia coli mutant temperature sensitive for group I RNA bacteriophages.

The temperature-sensitive conjugational transfer-deficient mutant Escherichia coli JCFL39, carrying a traD(Ts) mutation, is herein described as also being temperature sensitive for group I RNA phages (MS2, f2, and R17) but not for Q beta. Temperature shift experiments showed that the growth of group I phage MS2 in the mutant could be inhibited by a post-penetration event at high temperature. A possible role for the traD cistron of sex factor F in the intracellular development of MS2 is suggested.