Effect of base stacking on the relative thermodynamic stability of oligonucleotide complexes: a spectroscopic study

Three-strand oligonucleotide complexes are employed to assess the effect of base stacking and base pair mismatch on the relative thermodynamic stabilities of oligonucleotide duplexes. The melting behavior of three-strand oligonucleotide complexes incorporating nicks and gaps as well as internal single base mismatches is monitored using temperature-dependent optical absorption spectroscopy. A sequential three-state equilibrium model is used to analyze the measured melting profiles and evaluate thermodynamic parameters associated with dissociation of the complexes. The free-energy of stabilization of a nick complex compared to a gap complex due to base stacking is determined to be -1.9 kcal/mol. The influence of a mispaired base in these systems is shown to destabilize a nick complex by 3.1 kcal/mol and a gap complex by 2.8 kcal/mol, respectively.     

Thermodynamic analysis of stacking hybridization of oligonucleotides with DNA template.

Contiguous stacking hybridization of oligodeoxyribonucleotides with DNA as template was investigated using three types of complexes: oligonucleotide contiguously stacked with the stem of the preformed minihairpin (complexes I), oligonucleotide tandems containing two (complexes II) or three (complexes III) short oligomers with a common DNA template. Enthalpy Delta H degrees and entropy Delta S degrees of the coaxial stacking of adjacent duplexes were determined for GC/G*pC, GT/A*pC, AC/G*pT, AT/A*pT, CT/A*pG, AG/C*pT, AA/T*pT and TT/A*pA nicked (*) dinucleotide base pairs. The maximal efficiency of co-operative interaction was found for the GC/G*pC interface (Delta G degrees(NN/N*pN)=-2.7 kcal/mol) and the minimal one for the AA/T*pT interface (Delta G degrees(NN/N*pN)=-1.2 kcal/mol) at 37 degrees C. As a whole, the efficiency of the base pairs interaction Delta G degrees(NN/N*pN) in the nick is not lower than that within the intact DNA helix (Delta G degrees(NN/NN)).These observed Delta G degrees(NN/N*pN) values are proposed may include the effect of the partial removal of fraying at the adjacent helix ends additionally to the effect of the direct stacking of the terminal base pairs in the duplex junction (Delta G degrees(NN/NN). The thermodynamic parameters have been found to describe adequately the formation of all tandem complexes of the II and III types with oligonucleotides of various length and hybridization properties. The performed thermodynamic analysis reveals features of stacking oligonucleotide hybridization which allow one to predict the temperature dependence of association of oligonucleotides and the DNA template within tandem complexes as well as to determine optimal concentration for formation of these complexes characterized by high co-operativity level.     

The thermodynamic advantage of DNA oligonucleotide 'stacking hybridization' reactions: energetics of a DNA nick.

'Stacking hybridization reactions' wherein two or more short DNA oligomers hybridize in a contiguous tandem orientation onto a longer complementary DNA single strand have been employed to enhance a variety of analytical oligonucleotide hybridization schemes. If the short oligomers anneal in perfect head-to-tail register the resulting duplex contains a nick at every boundary between hybridized oligomers. Alternatively, if the short oligomers do not hybridize precisely in register, i.e. single strand regions on the longer strand are left unbound, gaps are formed between regions where short oligomers bind. The resulting gapped DNA duplexes are considerably less stable than their nicked duplex analogs. Formation of base pair stacking interactions between neighboring oligomers at the nicks that do not occur in gapped duplexes has been proposed as the source of the observed added stability. However, quantitative evidence supporting this hypothesis for DNA has not been reported. Until now, a direct comparison of the thermodynamics of DNA nicks versus DNA gaps has not been performed. In this communication we report such a comparison. Analysis of optical melting experiments in a well defined molecular context enabled quantitative evaluations of the relative thermodynamic difference between nicked and gapped DNA duplexes. Results of the analysis reveal that a nick may be energetically favored over a gap by at least 1.4 kcal/mol and perhaps as much as 2.4 kcal/mol. The presence of a 5'phosphate at a nick or gap fails to significantly affect their stabilities.     

Synthesis and hybridization properties of the conjugates of oligonucleotides and stabilization agents. Part 3

New compounds having tri- or pentamethylenamine linker functions were synthesized. These derivatives were covalently attached through the 5'-phosphoramide linkage to heptanucleotide pd(CCAAACA). Complementary complexes of the octanucleotide pd(TGTTTGGC) and above oligonucleotide conjugates were tested for their thermodynamic response. The T(m) data and thermodynamic parameters for complex formation confirmed the ability of chromone (gamma-pyrone) derivatives to stabilize strongly the 7-mer/8-mer complementary complex. Moreover, benzochromone (naphthopyrane) and, surprisingly, tetrahydropyrimidinethanone derivatives showed the capacity of stabilizing this 7-mer/8-mer complementary complex. The effect of all these compounds on the stability of the oligonucleotide complexes (DeltaDeltaG at 37 degrees C ranged from -1.2 to -2.0 kcal/mol) was shown to be comparable to the effect of one nucleotide base pair and similar to the effect (DeltaDeltaG at 37 degrees C ranged from -1.5 to -2.0 kcal/mol) found for acridine-oligonucleotide conjugates, which served as a reference in this study.     

The Influence of Nearest Neighbours on the Efficiency of Coaxial Stacking at Contiguous Stacking Hybridization of Oligodeoxyribonucleotides

Contiguous stacking hybridization of oligodeoxyribonucleotides with a stem of preformed minihairpin structure of a DNA template was studied with the use of UV‐melting technique. It was shown that the free‐energy of the coaxial stacking interaction (ΔG°_ST at 37°C, 1 M NaCl, pH 7.4) at the complementary interface XA^*pTY/ZATV (an asterisk stands for a nick) strongly depends on the type of nearest neighbor bases X and Y flanking the nicked dinucleotide step. The maximum efficiency of the coaxial stacking was observed for the PuA^*pTPy/PuATPy interface, whereas the minimum efficiency was obtained for the PyA^*pTPu/PyATPu interface. A 5′‐phosphate residue in the nick enhances the coaxial stacking. In dependence on duplex structure the observed efficiency of A^*T/AT coaxial stacking varied from (? 0.97 kcal/mol) for unphosphorylated TA^*TA/TATA interface to three‐fold higher value (? 2.78 kcal/mol) for GA^*pTT/AATC interface.     

Cooperative nonenzymic base recognition--a study of interaction of oligoguanylic acid with poly C

The interaction between poly C and (Gp)_nG(n = 1,2) in dilute solution was investigated spectrophotometrically in 0.1M phosphate buffer pH 7.2 under conditions unfavorable for the formation of self-associated complexes of oligoguanylic acids. Two isosbestic points were observed when poly C was titrated gradually with GpGpG, one at 232–233 mμ(in the range of 0–33% poly C) and one around 238 mμ (in the range of 50–100% poly C). The melting temperature (T_m) of the 1:1 poly C: (Gp)nG complexes (n = 1,2) of varying concentration were determined. The equilibrium properties of the 1:1 complexes can be described by two interaction parameters, namely, (i) cooperative stacking interaction between the first nearest neighbor of the adsorbed oligomer, and (ii) intrinsic association constant of the adsorbed oligomer with its polymeric site, since the cooperative helix–coil transition particularly in the smaller oligonucleotide can be described by an “all or none” model. Based on such a model the enthalpy of stacking inteaction-dependent T_m values yielded directly the sum of the enthalpy of stacking interaction and of basepairing (which is dependent on the chain length of the oligomer) and the value of S, the stability constant of a G–C pair within a helix. The enthalpy of formation of G–C pair is then calculated as ?6.3 kcal/base pair either from the chain length dependent enthalpy term or from the temperature coefficient of S values. From the S value and the association constant of 1:1 GpGpGpC:GpCpCpC complex, other thermodynamic parameters such as nucleation parameter (β) and free energy of stacking interaction can be obtained.     

Synthesis and high stability of complementary complexes of N-(2-hydroxyethyl)phenazinium derivatives of oligonucleotides.

Two simple methods for the synthesis of oligonucleotides bearing a N-(2-hydroxyethyl)phenazinium (Phn) residue at the 5'- and/or 3'-terminal phosphate groups are proposed. By forming complexes between a dodecanucleotide d(pApApCpCpTpGpTpTpTpGpGpC), a heptanucleotide d(pCpCpApApApCpA), and Phn derivatives of the latter, it is shown that the introduction of a dye at the end of an oligonucleotide chain strongly stabilizes its complementary complexes. The Tmax and the thermodynamic parameters (delta H, delta S, delta G) of complex formation were determined. According to these data, coupling of a dye with the 5'-terminal phosphate group is the most advantageous: delta G(37 degrees C) is increased by 3.59 +/- 0.04 kcal/mol compared to 2.06 +/- 0.04 kcal/mol for 3'-Phn derivatives. The elongation of the linker, which connects the dye to the oligonucleotide, from a dimethylene up to a heptamethylene usually leads to destabilization of the oligonucleotide complex. The complementary complex formed by the 3',5'-di-Phn derivative of the heptanucleotide was found to be the most stable among all duplexes investigated. Relative to the unmodified complex the increase in free energy was 4.96 +/- 0.04 kcal/mol. The association constant of this modified complex at 37 degrees C is 9.5 x 10(6) M-1, whereas the analogous value for the unmodified complex is only 3 x 10(3) M-1.     

Calorimetric unfolding of the bimolecular and i-motif complexes of the human telomere complementary strand, d(C(3)TA(2))(4)

A combination of spectroscopic and calorimetric techniques is used to determine the unfolding thermodynamics of the complexes formed by the complementary sequence of the human telomere, d(C(3)TA(2))(4), in the pH range of 4.2 to 6. Calorimetric melting curves show biphasic transitions; both transitions are shifted to higher temperatures as the pH is decreased, indicative of cytosine protonation, which favors the formation of C*C(+) base pairs. Furthermore, the transition temperature, T(M), of the lower transition depends on strand concentration, while the T(M) of the higher transition is independent of strand concentration, indicating the following sequential melting: bimolecular complex(s)-->intramolecular complex-->random coil. The thermodynamic profiles for the formation of each complex, bimolecular and i-motif reveals small favorable free energy terms resulting from favorable enthalpy-unfavorable entropy compensations, uptake of protons, marginal uptake of counterions (i-motif) and marginal release of water molecules (i-motif). Furthermore, an enthalpy of 3.2 kcal/mol (bimolecular complex) and 5.0 kcal/mol (i-motif) is estimated for a single C*C(+)/C*C(+) base-pair stack.     

Comparison of the base pairing properties of a series of nitroazole nucleobase analogs in the oligodeoxyribonucleotide sequence 5'-d(CGCXAATTYGCG)-3'.

The nucleoside analogs 1-(2'-deoxy-beta-D-ribofuranosyl)- 3-nitropyrrole (9), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitropyrazole (10), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitroimidazole (11) and 1-(2'-deoxy-beta-D-ribofuranosyl)-5-nitroindole (21) were incorporated into the oligonucleotide 5'-d(CGCXAATTYGCG)-3'in the fourth position from the 5'-end. Procedures for synthesis of two of the nitroazole nucleosides, 10 and 11, were developed for this study. Each of the nitroazoles was converted into a 3'-phosphoramidite for oligonucleotide synthesis by conventional automated protocols. Four oligonucleotides were synthesized for each modified nucleoside in order to obtain duplexes in which each of the four natural bases was placed opposite (position 9) the nitroazole. In order to assess the role of the nitro group on base stacking interaction, sequences were also synthesized in which the fourth base was 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole. Corresponding sequences containing an abasic site, as well as sequences containing inosine, were synthesized for comparison. Thermal melting studies yielded T m values and thermodynamic parameters. Each nucleoside analog displayed a unique pattern of base pairing preferences. The least discriminating analog was 3-nitropyrrole, for which T m values differed by 5 degrees C and Delta G 25 degrees C ranged from -6.1 to -6.5 kcal/mol. 5-Nitroindole gave duplexes with significantly higher thermal stability, with Tm values varying from 35.0 to 46.5 degrees C and -Delta G 25 degrees C ranging from 7.7 to 8.5 kcal/mol. Deoxyinosine (22), a natural analog which has found extensive use as a universal nucleoside, is far less non-discriminating than any of the nitroazole derivatives. Tm values ranged from 35.4 degrees C when paired with G to 62.3 degrees C when paired with C. The significance of the nitro substituent was determined by comparison of the base pairing properties of a simple azole nucleoside, 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole (12). The pyrazole-containing sequences melt at 10-20 degrees C lower than the corresponding nitropyrazole-containing sequences. On average, the pyrazole-containing sequences were equivalent in stability (average Delta G = -4.8 kcal/mol) to the sequences containing an abasic site (average Delta G = -4.7 kcal/mol).     

Calorimetric study of 2U:1A three-stranded complexes formed between poly(U) and adenine dinucleotides: ApA and diastereoisomers of nonionic adenine dideoxyribonucleoside methyl phosphonate

Melting parameters of 2U:1A complexes formed by polyuridylic acid [poly(U)] and three adenine dinucleotides, diribonucleoside monophosphonate ApA and diastereoisomers of dideoxyribonucleoside methyl phosphonate [(dApA)₁ and (dApA)₂], in 1M NaCl and at a number of dinucleotide concentrations were obtained from differential scanning microcalorimetric data and interpreted in terms of the theory of helix–coil equilibrium in oligonucleotide–polynucleotide systems. The apparent binding constant, 1/c_m, at 39°C and melting temperatures, T_m, at 1 × 10^?3 M dinucleotide concentration indicate the following order of thermodynamic stability of the complexes: 2 poly(U) · (dApA)₂ (2.27 × 10³M^?1, 44.2°C) > 2 poly(U) · (dApA)₁ (9.9 × 10²M¹, 39.2°C) > 2 poly(U) · (ApA) (5.9 × 10²M^?1, 35.8°C). Corresponding calorimetric enthalpies of melting, ΔH_m: 13.5, 12.7, and 12.8 kcal/mol (UUA base triplets) were found to be considerably lower than the van't Hoff enthalpies, ΔH_app: 29.4, 16.2, and 16.2 kcal/mol, respectively, evaluated from the dependence of the melting temperatures on dinucleotide concentration. Self-association of dinucleotides and their simultaneous binding as monomers, dimers, and higher-order associated species is suggested as the most probable cause of the differences between ΔH_m and ΔH_app values. The differences in thermodynamic properties of the complexes formed by (dApA)₁ and (dApA)₂ diastereoisomers are discussed in connection with their known conformational properties. The higher and essentially enthalpic stability of the 2 poly(U) · (dApA)₂ complex correlates with a lower degree of intramolecular stacking of the (dApA)₂ isomer. The hydrophobically enhanced strong self-association of the latter greatly influences the thermodynamics of its complex formation with poly(U) and results in ΔH_app/ΔH_m = 2.3.     

Cooperative nonenzymic base recognition II. thermodynamics of the helix-coil transition of oligoadenylic + oligouridylic acids

The properties of oligonucleotide helices of adeuylic- and uridylic acid oligomers have been investigated by measurements of hypo-and hyperchromieity. High ionic strengths favor the formation of triple helices. Thus, the double helix-coil transition can be studied (without interference by triple helices) only at low ionic-strength. A “phase diagram” is given representing the T_m-values of the various transitions at different ionic strengths for the system A(pA)₁₇ + U(pU)₁₇. Oligonucleolides of chain lengths <8 always form both double and triple helices at the nucleotide concentrations required for base pairing. For this reason the double helix-coil transition without coupling of the triple helix equilibrium can only be measured for chain lengths higher than 7. Melting curves corresponding to this transition have been determined for chain lengths 8, 9, 10, 11, 14 and 18 at different concentrations. An increase in nucleotide concentration leads to an increase in melting temperature. The shorter the chain length the lower the T_m-value and the broader the helix-coil transition. The experimental transition curves have been analysed according to a staggering zipper model with consideration of the stacking of the adeuylic acid single strands and the electrostatic repulsion of tlip phosphate charges on opposite strands. The temperature dependence of the nucleation parameter has been accounted for by a slacking factor x. The stacking factor expresses the magnitude of the stacking enthalpy. By curve fitting xwas computed to be 0.7, corresponding to a stacking enthalpy of about S kcal/mole. The model described allows the reproduction of the experimental transition curves with relatively high accuracy. In an appendix the thermodynamic parameters of the stacking equilibrium of poly A and of the helix-coil equilibria of poly A + poly U at neutral pH are calculated (ΔH_A = ?7.9 kcal/mole for the poly A stacking and ΔH₁₂ = ?10.9 kcal/mole for the formation of the double helix from the randomly coiled single strands). A formula for the configurational entropy of polymers derived by Flory on the basis of a liquid lattice model is adapted to calculate the stacking entropies of adenylic oligomers.     

Thermodynamic contribution and nearest-neighbor parameters of pseudouridine-adenosine base pairs in oligoribonucleotides

Pseudouridine (Ψ) is the most common noncanonical nucleotide present in naturally occurring RNA and serves a variety of roles in the cell, typically appearing where structural stability is crucial to function. Ψ residues are isomerized from native uridine residues by a class of highly conserved enzymes known as pseudouridine synthases. In order to quantify the thermodynamic impact of pseudouridylation on U-A base pairs, 24 oligoribonucleotides, 16 internal and eight terminal Ψ-A oligoribonucleotides, were thermodynamically characterized via optical melting experiments. The thermodynamic parameters derived from two-state fits were used to generate linearly independent parameters for use in secondary structure prediction algorithms using the nearest-neighbor model. On average, internally pseudouridylated duplexes were 1.7 kcal/mol more stable than their U-A counterparts, and terminally pseudouridylated duplexes were 1.0 kcal/mol more stable than their U-A equivalents. Due to the fact that Ψ-A pairs maintain the same Watson-Crick hydrogen bonding capabilities as the parent U-A pair in A-form RNA, the difference in stability due to pseudouridylation was attributed to two possible sources: the novel hydrogen bonding capabilities of the newly relocated imino group as well as the novel stacking interactions afforded by the electronic configuration of the Ψ residue. The newly derived nearest-neighbor parameters for Ψ-A base pairs may be used in conjunction with other nearest-neighbor parameters for accurately predicting the most likely secondary structure of A-form RNA containing Ψ-A base pairs.     

Free Energy Calculations Predict Sequence Specificity in DNA-Drug Complexes

Abstract

Molecular dynamics simulations of DNA-netropsin complexes in water were performed using the thermodynamic cycle-perturbation method to calculate the free energy difference between complexes with an adenine-containing binding site and corresponding complexes where adenines are replaced by 2,6-diaminopurines (dap). The calculations predict a free energy difference of 3.7±0.9 kcal/mol (at 300K) in favour of netropsin binding to an (AATT)₂ DNA sequence compared to a (dapdapTT)₂ sequence.     

Complex of human apolipoprotein C-1 with phospholipid: thermodynamic or kinetic stability?

Thermal unfolding of discoidal complexes of apolipoprotein (apo) C-1 with dimyristoyl phosphatidylcholine (DMPC) reveals a novel mechanism of lipoprotein stabilization that is based on kinetics rather than thermodynamics. Far-UV CD melting curves recorded at several heating/cooling rates from 0.047 to 1.34 K/min show hysteresis and scan rate dependence characteristic of slow nonequilibrium transitions. At slow heating rates, the apoC-1 unfolding in the complexes starts just above 25 degrees C and has an apparent melting temperature T(m) approximately 48 +/- 1.5 degrees C, close to T(m) = 51 +/- 1.5 degrees C of free protein. Thus, DMPC binding may not substantially increase the low apparent thermodynamic stability of apoC-1, DeltaG(25 degrees C) < 2 kcal/mol. The scan rate dependence of T(m) and Arrhenius analysis of the kinetic data suggest an activation enthalpy E(a) = 25 +/- 5 kcal/mol that provides the major contribution to the free energy barrier for the protein unfolding on the disk, DeltaG > or = 17 kcal/mol. Consequently, apoC-1/DMPC disks are kinetically but not thermodynamically stable. To explore the origins of this kinetic stability, we utilized dynode voltage measured in CD experiments that shows temperature-dependent contribution from UV light scattering of apoC-1/DMPC complexes (d approximately 20 nm). Correlation of CD and dynode voltage melting curves recorded at 222 nm indicates close coupling between protein unfolding and an increase in the complex size and/or lamellar structure, suggesting that the enthalpic barrier arises from transient disruption of lipid packing interactions upon disk-to-vesicle fusion. We hypothesize that a kinetic mechanism may provide a general strategy for lipoprotein stabilization that facilitates complex stability and compositional variability in the absence of high packing specificity.     

Interaction of poly-5-bromouridylic acid. II. Thermodynamic analysis of its interaction with adenosine and 9-methyladenine

The interaction of poly-5-bromouridylic acid [poly(BU)] with adenosine and 9-methyladenine was studied by equilibrium dialysis, optical melting, and microcalorimetry. The stacking free energy, ω, was estimated as ?17.6 kJ/mol for adenosine·2poly(BU) and ?18.8 kJ/mol for 9-methyladenine·2poly(BU) from the binding isotherms constructed from equilibrium dialysis results. The binding isotherms constructed from a series of melting curves also gave ω values for adenosine·2poly(BU). The thermal stability of the complex depends on monomer concentration, and the partial molar enthalpies of the complex formation at the midpoint of the transition were evaluated from the T_m coefficients as a function of free monomer concentration. The values of ?92.0 and ?90.4 kJ/mol were obtained for adenosine·2poly(BU) and 9-methyladenine·2poly(BU) in 0.4M NaCl–0.02M Na-cacodylate–5 × 10^?4M EDTA (pH 7.0), respectively. Microcalorimetric measurements provided lower integral heats of reaction values for these complexes, i.e., ?73.2 kJ/mol for adenosine·2poly(BU) and ?71.5 kJ/mol for 9-methyladenine·2poly(BU). A comparison with a polyribouridylic acid system provided a quantitative understanding of a stabilization by bromination in terms of thermodynamic parameters.     

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.     

Hydrogen-bonded Trimers of DNA Bases and their Interaction with Metal Cations: Ab initio Quantum-chemical and Empirical Potential Study

Abstract

Neutral (G.GC, A. AT, G.AT, T. AT, and C (imino).GC) and protonated (CH⁺.GC and AH⁺.GC) hydrogen-bonded trimers of nucleic acid bases were characterized by ab initio methods with the inclusion of electron correlation. In addition, the influence of metal cations on the third-strand binding in Purine-Purine-Pyrimidine (Pu.PuPy) reverse-Hoogsteen triplets has been studied. The ab initio calculations were compared with those from recently introduced force fields (AMBER4.1, CHARMM23, and CFF95). The three-body term in neutral trimers is mostly negligible, and the use of empirical potentials is justified. The only exception is the neutral G.GC Hoogsteen trimer with a three-body term of -4 kcal/mol. Protonated trimers are stabilized by molecular ion—;molecular dipole attraction and the interaction within the complex is nonadditive, with the three-body term on the order of -3 kcal/mol. There is a significant induction interaction between the third-strand protonated base and guanine. The calculations indicate an enhancement of the third-strand binding in the G.GC reverse-Hoogsteen trimer due to metal cation coordination to the N7/06 position of the third-strand guanine. Interactions between metal cations and complexes of DNA bases are in general highly nonadditive; the three-body term is above -10 kcal/mol in a complex of a divalent cation (Ca²⁺) with the GG reverse-Hoogsteen pair. The pairwise additive empirical potentials qualitatively underestimate the binding energy between cation and base.     

A remarkably stable kissing-loop interaction defines substrate recognition by the Neurospora Varkud Satellite ribozyme

Kissing loops are tertiary structure elements that often play key roles in functional RNAs. In the Neurospora VS ribozyme, a kissing-loop interaction between the stem–loop I (SLI) substrate and stem–loop V (SLV) of the catalytic domain is known to play an important role in substrate recognition. In addition, this I/V kissing-loop interaction is associated with a helix shift in SLI that activates the substrate for catalysis. To better understand the role of this kissing-loop interaction in substrate recognition and activation by the VS ribozyme, we performed a thermodynamic characterization by isothermal titration calorimetry using isolated SLI and SLV stem–loops. We demonstrate that preshifted SLI variants have higher affinity for SLV than shiftable SLI variants, with an energetic cost of 1.8–3 kcal/mol for the helix shift in SLI. The affinity of the preshifted SLI for SLV is remarkably high, the interaction being more stable by 7–8 kcal/mol than predicted for a comparable duplex containing three Watson–Crick base pairs. The structural basis of this remarkable stability is discussed in light of previous NMR studies. Comparative thermodynamic studies reveal that kissing-loop complexes containing 6–7 Watson–Crick base pairs are as stable as predicted from comparable RNA duplexes; however, those with 2–3 Watson–Crick base pairs are more stable than predicted. Interestingly, the stability of SLI/ribozyme complexes is similar to that of SLI/SLV complexes. Thus, the I/V kissing loop interaction represents the predominant energetic contribution to substrate recognition by the trans-cleaving VS ribozyme.     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.