Free energy of imperfect nucleic acid helices. II. Small hairpin loops

Physical studies of enzymically synthesized oligonucleotides of defined sequence are used to evaluate quantitatively the stability of small RNA hairpin loops and helices. The series (Ap)₄G(pC) _N(pU)₄, N = 4, 5 or 6, exists as monomolecular hairpin helices when N ≥ 5, and as imperfect dimer helices when N ≤ 4. In this size range, hairpin loops become more favorable (less destabilizing thermodynamically) as they increase in size from 3 to 4 to 5 unbonded nucleotides. Very small hairpin loops are particularly destabilizing; molecules whose base sequence would imply a hairpin loop of three nucleotides will generally exist with a loop of five, including a broken terminal base pair.Thermodynamic parameters for base pair and loop formation are calculated by a method which makes unnecessary the use of measured enthalpies of polynucleotide melting. Literature data on oligonucleotide double helices yield estimates of the free energy contribution from each of the six types of stacking interactions between three possible neighboring base pairs. The advantage of this approach is that the properties of oligonucleotides are used in predicting the stability of small RNA helices, avoiding the long extrapolation from the properties of high polymers.We provide Tables of temperature-dependent free energies that allow one to predict the stability and thermal transition temperature of many simple RNA secondary structures (applicable to ~1 m-Na⁺ concentration). As an example, we apply the rules to an isolated fragment of tRNA^Ser (yeast) (Coutts, 1971), whose properties were not used in calculating the free-energy parameters. The experimental melting temperature of 88 °C is predicted with an error margin of 5 deg. C.     

Topography of nucleic acid helices in solutions. 3. Interactions of spermine and spermidine derivatives with polyadenylic-polyuridylic and polyinosinic-polycytidylic acid helices

The synthesis of spermine derivatives (II), \documentclass{article}\pagestyle{empty}\begin{document}$ {\rm R}_1 {\rm R}_{\rm 2} {\rm R}_{\rm 3} \mathop {\rm N}\limits^ + \left( {{\rm CH}_2 } \right)_3 \mathop {\rm N}\limits^ + {\rm R}_{\rm 1} {\rm R}_{\rm 2} \left( {{\rm CH}_2 } \right)_2 ]_2 \cdot 4{\rm X}^ - $\end{document}, and spermidine derivatives (III), \documentclass{article}\pagestyle{empty}\begin{document}$ {\rm R}_1 {\rm R}_{\rm 2} {\rm R}_{\rm 3} \mathop {\rm N}\limits^ + \left( {{\rm CH}_2 } \right)_4 \mathop {\rm N}\limits^ + {\rm R}_{\rm 1} {\rm R}_{\rm 2} \left( {{\rm CH}_2 } \right)_3 \mathop {\rm N}\limits^ + {\rm R}_{\rm 1} {\rm R}_{\rm 2} {\rm R}_3 \cdot 3{\rm X}^ - $\end{document}, are reported. The effects of these salts on the helix–coil transition of rA–rU and rI–rC helices were examined. Increasing the size of the hydrophobic substituents, R¹, R², and R³ lowers the degree of stabilization of the helical structure. The disproportionation reaction, 2rA–rU→rA–rU₂ + rA occurs readily with salts II and III, especially when the substituents, R₁, R₂, and R₃ are small, i.e., H or Me. Spermine is found to stabilize the rA–rU² and rI–rC helices to approximately the same extent; however, large differences between the degree of stabilization of rA–rU₂ and rI-rC helices are observed when the substituents R₁, R₂, and R₃ are large hydrophobic groups. Similar results are also obtained for the spermidine series. Finally, differences in the interactions of the salts II and III with rA–rU₂ and rI–rC helices suggest that the latter helix is denser.     

Ring current shielding effects in nucleic acid double helices.

Values of ring current shielding in parts-per-million have been calculated for double helical nucleic acids in the A-RNA (RNA-11), A' -RNA (RNA-12) and B-DNA geometries. Atomic coordinates determined previously from x-ray diffraction of fibers were used to calculate the positions of protons relative to both nearest and second nearest neighboring bases, including those on the complementary strand. The magnitude of the diamagnetic shielding was then calculated for each aromatic ring. From these calculations tables were constructed for use in determining the shielding expected for carbon-bound and ring nitrogen-bound protons of any double helical nucleic acid sequence. The results are compared with available experimental data for several oligonucleotides and with previous ring current shielding calculations where differences of up to 0.4 ppm are found.     

Effects of GA mismatches on the structure and thermodynamics of RNA internal loops

UV melting, CD and NMR studies indicate rGCGAGCG and rGCAGGCG from unusually stable duplexes of type a and b. The observed delta G degree 37's in 1 M NaCl are -6.7 and -6.3 kcal/mol, respectively. For the related duplex, c, delta G degree 37 is -4.2 kcal/mol. The predicted delta G degree 37 from nearest-neighbor parameters (formula; see text) for all three duplexes is -4.7 kcal/mol (Freier, S.M., Kierzek, R., Jaeger, J.A., Sugimoto, N., Caruthers, M.H., Neilson, T., & Turner, D.H. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9373-9377). The results suggest a special interaction in the duplexes containing GA mismatches. Presumably, this is hydrogen bonding between G and A. While the thermodynamics for (rGCGAGCG)2 and (rGCAGGCG)2 are similar, CD and the imino region of the proton NMR spectra indicate their structures are different. In particular, (rGCAGGCG)2 exhibits a CD spectrum typical of A-form geometry with a weak negative band at 280 nm. In contrast, the CD spectrum for (rGCGAGCG)2 has an intense positive band at 285 nm. The NMR spectrum of (rGCAGGCG)2 has a resonance corresponding to a hydrogen-bonded GA mismatch, while for (rGCGAGCG)2 no hydrogen-bonded imino proton is observed for the mismatch. The glycosidic torsion angles of the bases in the GA mismatches of (rGCAGGCG)2 and (rCGCAGGCG)2 are anti. Duplexes of type d, where X is A, G, or U, are more stable than e, and the stability differences are similar to those (formula; see text) observed for f versus g. Thus, 3'-dangling ends in this system make contributions to duplex stability that are similar to contributions observed with fully paired duplexes.     

Free energy contributions of G.U and other terminal mismatches to helix stability

Thermodynamic parameters of helix formation were measured spectroscopically for seven hexaribonucleotides containing a GC tetramer core and G.U or other terminal mismatches. The free energies of helix formation are compared with those for the tetramer core alone and with those for the hexamer with six Watson-Crick base pairs. In 1 M NaCl, at 37 degrees C, the free energy of a terminal G.U mismatch is about equal to that of the corresponding A.U pair. Although other terminal mismatches studied add between -1.0 and -1.6 kcal/mol to delta G0 37 for helix formation, all are less stable than the corresponding Watson-Crick pairs. Comparisons of the stability increments for terminal G.U mismatches and G.C pairs suggest when stacking is weak the additional hydrogen bond in the G.C pair adds roughly -1 kcal/mol to the favorable free energy of duplex formation.     

Interactions of nucleic acid double helices induced by electric field pulses

Electric field pulses induce a substantial increase of the light scattering intensity of double-helical DNA. The relative change of light scattering and also the reciprocal relaxation time constants under electric field pulses increase with increasing nucleotide concentration. These observations, together with a large difference between dichroism orientation time constants and light scattering time constants under electric field pulses, demonstrate that the main part of the light scattering effect is due not to field-induced orientation but to interactions between DNA helices. From the concentration dependence of the light scattering time constants we obtain, according to an isodesmic reaction model, association rate constants in the range 3 × 10¹⁰ M^?1 helices s^?1 for DNA with approx. 300 base-pairs. These values are at the limit of a diffusion-controlled DNA association and do not show any dependence upon the field strength. The dissociation rate constants k_d decrease strongly with increasing field strength E and thus demonstrate that the interactions between the helices are induced by the electric field. This conclusion is consistent with independent measurements which do not reveal any DNA association at zero field strength. The observed linear relation between log(k_d) and E² suggests a field-induced reaction driven by dipole changes. According to this interpretation the change of dipole moment should be in the range of approx. 1400 debye. The dissociation rates for DNA helices with approx. 300 to approx. 800 base-pairs strongly increase with increasing sail concentration (measured in the range 1–5 mM ionic strength), whereas the association rate constants remain virtually unchanged. Measurements of the linear dichroism in the same range of DNA chain length demonstrate that for long field pulses of e.g., 40 μs, the amplitude approaches a maximum value and then decreases. The dichroism relaxation curves observed after long field pulses exhibit a component with a positive dichroism and an increased decay time. These observations suggest the formation of a DNA aggregate with an unusual arrangement of the bases.     

Topography of nucleic acid helices in solutions. The interaction specificities of optically active amino acid derivatives

The synthesis and interactions of the d- and l-enantiomers of the amino acid amide derivatives [Formula: see text] (I) and lysyl dipeptides [Formula: see text] (II) with poly rI.poly rC, poly rA.poly rU and calf thymus DNA is reported. The following results were found. (1) The degree of stabilization of the helices as measured by the T(m) (;melting' temperature) of the helix-coil transition was dependent on the nature of the amino acid. (2) For the poly rI.poly rC helix, the l-enantiomers of salts (I) and (II) stabilized more than the d-enantiomers. The same was true for calf thymus DNA in the presence of salts (II) and for poly rA.poly rU in the presence of salts (II) and the proline derivatives of salts (I). (3) As R increased in size and became more apolar, the amount of stabilization of the poly rI.poly rC helix in the presence of salts (I) decreased. On the other hand, the amount of stabilization increased with more polar substituents. An attempt was then made to determine whether the difference in stabilization of the double-stranded helices at the T(m) in the presence of the l- and d-enantiomers of salts (I) is due to the interaction with the helix, the random coil or both. A new method was developed for determining the binding of small ions to polyions that involves a competition between an insoluble polystyrene ion-exchange resin and the soluble polyion for the counterion. Dissociation constants are obtained for the complexes of single- and double-stranded helices with the salts (I). The results are illuminating and indicate that with certain helices, i.e. poly rA.poly rU, the interactions of salts (I) with the single strands may not be ignored. It is concluded that the high optical specificity found in Nature, i.e. d-ribose in nucleic acids and l-amino acids in proteins, cannot be attributed solely to monomer-polymer interactions described by Gabbay (1968).     

Interaction of the ruthenium red cation with nucleic acid double helices

The hexapositive complex cation ruthenium red very effectively stabilizes DNA and RNA double helices against thermal denaturation. In the presence of nucleic acid helices, this symmetric cation acquires an extrinsic CD spectrum near the wavelength of the dye's maximum absorbance. Competition experiments with single-stranded polyd(T) show this induced CD to be the result of selective binding to helical sites. The preferential affinity of ruthenium red for double helical binding sites is so great that it brings about biphasic absorbance- temperature profiles of polyd(A-T) at low [cation]: [polynucleotide phosphate]. The visible CD signal and fraction of helix melting at the upper transition increases with ruthenium red concentration until approximate charge neutrality is reached. These interactions, which have been studied in detail with the poly(U-U) helix as well as polyd(A-T), are likely largely electrostatic, since sufficient [NaCl] eliminates the bipliasic melting of polyd(A-T), renders the ultraviolet absorbance of poly (U) insensitive to ruthenium red, and abolishes the induced CD effects. The bipliasic melting of polyd(A-T) at intermediate [dye] is attributed to saturation of remaining double helical segments by cation migration from newly melted regions- Furthermore, virtually no change was observed in the induced CD upon melting through the first transition, whereas the effect is destroyed upon inciting through the second transition. A quantitative treatment of the data is used to obtain binding site size and association constant for the complex. The induced effect may prove useful in the exploration of exposed nucleic acid helical structure in such complex particles as nucleosomes or ribosomes.     

Recognition of nucleic acid double helices by homopyrimidine 2', 5'-linked RNA.

We have studied the effect of a 2',5'-RNA third strand backbone on the stability of triple helices with a 'pyrimidine motif' targeting the polypurine strand of duplex DNA, duplex RNA and DNA/RNA hybrids. Comparative experiments were run in parallel with DNA and the regioisomeric RNA as third strands adopting the experimental design of Roberts and Crothers. The results reveal that 2',5'-RNA is indeed able to recognize double helical DNA (DD) and DNA (purine):RNA (pyrimidine) hybrids (DR). However, when the duplex purine strand is RNA and the duplex pyrimidine strand is DNA or RNA (i.e. RD or RR), triplex formation is not observed. These results exactly parallel what is observed for DNA third strands. Based on T m data, the affinities of 2',5'-RNA and DNA third strands towards DD and DR duplexes were similar. The RNA third strand formed triplexes with all four hairpins, as previously demonstrated. In analogy to the arabinose and 2'-deoxyribose third strands, the possible C2'- endo pucker of 2',5'-linked riboses together with the lack of an alpha-2'-OH group are believed to be responsible for the selective binding of 2',5'-RNA to DD and DR duplexes, over RR and RD duplexes. These studies indicate that the use of other oligonucleotide analogues will prove extremely useful in dissecting the contributions of backbone and/or sugar puckering to the recognition of nucleic acid duplexes.     

Analysis of birefringence decay profiles for nucleic acid helices possessing bends: the tau-ratio approach.

For nucleic acid helices in the 100-200-bp range, a central bend or point of flexibility increases the rate of rotational diffusion. In a transient electric birefringence (TEB) experiment, this increase is manifest as a reduction in the terminal (slowest) birefringence decay time. Previous experimental and theoretical work has demonstrated that the ratio of the decay times for a bent/flexible molecule and its fully duplex (linear) counterpart represents a sensitive, quantifiable measure of the apparent bend angle (tau-ratio approach). In the current work, we have examined the influence of helix parameters (e.g., persistence length, helix rise, diameter) on the tau-ratio for a given bend. The tau-ratio is found to be remarkably insensitive to variations and/or uncertainties in the helix parameters, provided that one employs bent and control molecules with the same sequence and length (apart from the bend itself). Although a single tau-ratio determination normally does not enable one to distinguish between fixed and flexible bends, such a distinction can be made from a set of tau-ratios for molecules possessing two variably phased bends. A number of additional uncertainties are examined, including errors in the estimation of the dimensions of nonhelix elements that are responsible for bends; such errors can, in principle, be estimated by performing a series of measurements for molecules of varying length.     

Thermodynamics of internal C.T mismatches in DNA.

Thermodynamics of 23 oligonucleotides with internal single C.T mismatches were obtained by measuring UV absorbance as a function of temperature. Results from these 23 duplexes were combined with three measurements from the literature to derive nearest-neighbor thermodynamic parameters for seven linearly independent trimer sequences with internal C.T mismatches. The data show that the nearest-neighbor model is adequate for predicting thermodynamics of oligonucleotides with internal C.T with average deviations for Delta G degrees37, Delta H degrees, Delta S degrees and T m of 6.4%, 9.9%, 10.6%, and 1.9 degreesC respectively. C.T mismatches destabilize the duplex in all sequence contexts. The thermodynamic contribution of C. T mismatches to duplex stability varies weakly depending on the orientation of the mismatch and its context and ranges from +1.02 kcal/mol for GCG/CTC and CCG/GTC to +1.95 kcal/mol for TCC/ATG.     

Minimum energy conformations of proline-containing helices.

Proline occurs frequently in transmembrane alpha-helices of transport and receptor proteins even though statistical surveys demonstrate the overwhelming preference of this residue for a non-alpha-helical, hydrophilic environment. As a result, membrane-buried proline has been proposed to be functionally important, with function arising from structural discontinuity or destabilization of the helix. Destabilization may occur by Pro-mediated conformational transitions between discrete states, and may be manifested in membrane protein systems through reversible processes such as channel opening and closing or signal transduction. In this study, computer modeling of a model transmembrane alpha-helix, (Ala)8-Leu-Pro-Phe-(Ala)8, in a medium of low polarity (dielectric = 2), is used to examine the occurrence and energetic accessibility of Pro-mediated conformational interconversions. Leu psi and chi 1, Pro psi, and Phe phi and chi 1 torsion angles were assigned random values so that a data base of 200 conformations for each of the cis and trans states was generated. The conformations were minimized and low-energy structures organized into families. This analysis demonstrated that the most populated lowest energy family is the Trans-I conformation, corresponding to proline in a kinked alpha-helix. Two additional trans structures, Trans-II and Trans-III, as well as a cis conformation, Cis-I, are also energetically competitive. Interconversions between the trans states could thus be mediated by changes at a single torsion angle, accompanied by minor local hydrogen-bonding rearrangements. This work substantiates that membrane-buried proline can provide the basis for conformational transitions between discrete alpha-helix-based structures in a nonpolar environment.     

The energy of formation of internal loops in triple-helical collagen polypeptides

The sequence-dependent local destabilization in the interior of the collagen triple helix has been evaluated by means of conformational energy computations. Using a model poly(Gly-Pro-Pro) triple helix as the reference state, a method was developed for generating local loops, i.e., internal deformations, and analyzing their conformations. A seven-residue Gly-Pro-Pro-Gly-Pro-Pro-Gly fragment was replaced by the Gly-Pro-Ala-Gly-Ala-Ala-Gly sequence in one, two, or all three of the strands of the loop region. A set of loop conformations was generated in which the ends of the loop were initially fixed in the triple-helical structure. The potential energy of the entire deformed triple helix was then minimized, resulting in a variety of structures that contained deformed loops. The conformations of the triple helices at the two ends of the loops remained essentially unchanged in many of the low-energy conformations. In numerous high-energy conformations, however, the triple-helical segments were also partially or totally disrupted. The minimum-energy conformations of the whole structures were compared in terms of rms deviations of atomic coordinates with respect to the original triple helix, and of the shapes of the loops (using a distance function derived from differential geometry). Three new geometrical parameters—stretch S, kink K, and unwinding U—were defined to describe the changes in the overall orientation of the triple helices at the two ends of the loop. It is shown that, when the number of Pro residues in a short fragment is reduced, the triple helical structure can accomodate internal loops (i.e., distortions) within a 5 kcal/mol cutoff from the essentially unperturbed triple helical structure. For structures with a Gly-Pro-Ala-Gly-Ala-Ala-Gly sequence in all three strands, the probability of finding conformations with internal loops is small, i.e., 0.06. Internal loops affect the overall orientation of these structures, as measured by the helix-distortion parameters S, K, and U. © 1995 John Wiley & Sons, Inc.     

Thermodynamics of RNA hairpins containing single internal mismatches.

Thermodynamic parameters and circular dichroism spectra are presented for RNA hairpins containing single internal mismatches in the stem regions. Three different sequence contexts for the G*U mismatch and two contexts for C*A, G*A, U*U, A*C and U*G mismatches were examined and compared with Watson-Crick base-pair stabilities. The RNA hairpins employed were a microhelix and tetraloop representing the Escherichia coli tRNAAlaacceptor stem and sequence variants that have been altered at the naturally occurring G*U mismatch site. UV melting studies were carried out under different conditions to evaluate the effects of sodium ion concentration and pH on the stability of mismatch-containing hairpins. Our main findings are that single internal mismatches exhibit a range of effects on hairpin stability. In these studies, the size and sequence of the loop and stem are shown to influence the overall stability of the RNA, and have a minor effect on the relative mismatch stabilities. The relationship of these results to RNA-ligand interactions involving mismatch base-pairs is discussed.     

Topography of nucleic acid helices in solutions. II. Structure of the double-stranded rA-rU, rI-rC, acid rA, and the triple-stranded rA-rU2 and rA-rI2 helices

Information concerning the structures of rA–rU, rA–rU₂ rI–rC, rA–rI₂, and acid rA helices in solutions is reported. Through the use of diquaternary ammonium salts of the general structure, \documentclass{article}\pagestyle{empty}\begin{document}$ {\rm R}_1 {\rm R}_2 {\rm R}_3 \mathop {\rm N}\limits^ + ({\rm CH}_2 )n\mathop {\rm N}\limits^ + {\rm R}_1 {\rm R}_2 {\rm R}_3 \cdot 2{\rm Br}^ - $\end{document} (I), it is shown that (1) the distances between adjacent negatively charged oxygen atoms on the helix increases in the following order rA–rI₂ < rI–rC < rA–rU ? rA–rU₂; (2) the density of the helices increases in the order. rA–rI₂ < rA–rU < rA–rU₂ < rI–rC; (3) there is a large hydrophobia site in rA–rI₂ and possibly also in rA–rU, rA–rU₂, and rI–rC helices; (4) the results of the interactions between the salts of type I and the helices may be formulated in semi-quantitative terms by the use of two parameters, α, and β which are shown to be related to the charge separation and the density of the helices, respectively; (5) the studies in solutions compare favorably with the x-ray studies on the fibers; and (6) the acid rA helix differs significantly from the other helices by the fact that the electrostatic interstrand interactions between the negatively charged oxygen atom of a phosphate group and the positively charged 10-amino group of adenine contribute significantly to the stabilization of the helix, and thus it is found that the presence of the salts, I, leads to a significant destabilization of the acid rA helix.     

Topography of nucleic acid helices in solutions. 28. Evidence for a dynamic structure of DNA in solution

Proton magnetic resonance (pmr), ultraviolet absorption, induced circular dichroism (CD), and viscometric evidence is presented which show that reporter molecules $\text{1}$ and $\text{2}$ bind to DNA via an intercalation process. Preliminary kinetic studies show that the DNA·$\text{1}$ complex forms rapidly (i.e., <1 msec), whereas the DNA·$\text{2}$ complex forms at a considerably slower rate ($\text{t}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{> 100 msec}$). The kinetic results, and the steric requirements for intercalation of $\text{2}$ can be explained on the basis of a dynamic structure of DNA.