Coupling of sequential transitions in a DNA double hairpin: energetics, ion binding, and hydration

Effect of Internal and Bulge Loops on the Thermal Stability of Small DNA Duplexes

The thermal stabilities of DNA duplexes analogous to the let-7 microRNA:lin-41 mRNA complex from Caenorhabditis elegans have been measured by free solution capillary electrophoresis. DNA duplexes with the same stems but different types of internal or bulge loops and a control with no loop have also been studied. The melting temperatures of the DNA derivatives increased linearly with the logarithm of the Na+ or K+ ion concentration in the solution. Peaks in the electropherograms corresponding to duplexes with internal or bulge loops exhibited extensive tailing at high temperatures, suggesting that denaturation occurred by slow exchange between the duplexes and their component single strands. The single strands did not separate completely from the duplexes in aqueous solutions; instead, they appeared as small subpeaks on the tails of the duplex peaks. However, complete separation of the duplexes from their component single strands was observed at 20 °C in solutions containing 300 mM tetrapropylammonium ions. In addition, counterion condensation appears to be significantly reduced in DNA duplexes containing internal or bulge loops.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets

We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C₅ loop closing the triplex stem was replaced with ^5'-CT_nC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a ^5'-GCAA or ^5'-AACG tetraloop. For the triplexes with a ^5'-CT_nC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.

Investigation of the Melting Behavior of DNA Three-Way Junctions in the Closed and Open States

Intramolecular three-way junctions are commonly found in both DNA and RNA. These structures are functionally relevant in ribozymes, riboswitches, rRNA, and during replication. In this work, we present a thermodynamic description of the unfolding of DNA intramolecular three-way junctions. We used a combination of spectroscopic and calorimetric techniques to investigate the folding/unfolding thermodynamics of two three-way junctions with a closed (Closed-J) or open (Open-J) junction and their appropriate control stem-loop motifs (GAAATT-Hp, CTATC-Hp, and Dumbbell). The overall results show that both junctions are stable over a wide range of salt concentrations. However, Open-J is more stable due to a higher enthalpy contribution from the formation of a higher number of basepair stacks whereas Closed-J has a defined structure and retains the basepair stacking of all three stems. The comparison of the experimental results of Closed-J and Open-J with those of their component stem-loop motifs allowed us to be more specific about their cooperative unfolding. For instance, Closed-J sacrifices thermal stability of the Dumbbell structure to maintain an overall folded state. At higher salt concentration, the simultaneous unfolding of the above domains is lost, resulting in the unfolding of the three separate stems. In contrast, the junction of Open-J in low salt retains the thermal and enthalpic stability of the Dumbbell structure although sacrificing stability of the CTATC stem. The relative stability of Dumbbell is the primary reason for the higher ΔG°₍₅₎, or free energy, value seen for Open-J at low salt. Higher salt not only maintains thermal stability of the Dumbbell structure in Open-J but causes the CTATC stem to fully fold.

Thermodynamic profiles and nuclear magnetic resonance studies of oligonucleotide duplexes containing single diastereomeric spiroiminodihydantoin lesions

The spiroiminodihydantoins (Sp) are highly mutagenic oxidation products of guanine and 8-oxo-7,8-dihydroguanine in DNA. The Sp lesions have recently been detected in the liver and colon of mice infected with Helicobacter hepaticus that induces inflammation and the development of liver and colon cancers in murine model systems [Mangerich, A., et al. (2012) Proc. Natl. Acad. Sci. U.S.A. 109, E1820-E1829]. The impact of Sp lesions on the thermodynamic characteristics and the effects of the diastereomeric Sp-R and Sp-S lesions on the conformational features of double-stranded 11-mer oligonucleotide duplexes have been studied by a combination of microcalorimetric methods, analysis of DNA melting curves, and two-dimensional nuclear magnetic resonance methods. The nonplanar, propeller-like shapes of the Sp residues strongly diminish the extent of local base stacking interactions that destabilize the DNA duplexes characterized by unfavorable enthalpy contributions. Relative to that of an unmodified duplex, the thermally induced unfolding of the duplexes with centrally positioned Sp-R and Sp-S lesions into single strands is accompanied by a smaller release of cationic counterions (Δn(Na⁺) = 0.6 mol of Na⁺/mol of duplex) and water molecules (Δn(w) = 17 mol of H₂O/mol of duplex). The unfolding parameters are similar for the Sp-R and Sp-S lesions, although their orientations in the duplexes are different. The structural disturbances radiate one base pair beyond the flanking C:G pair, although Watson-Crick hydrogen bonding is maintained at all flanking base pairs. The observed relatively strong destabilization of B-form DNA by the physically small Sp lesions is expected to have a significant impact on the processing of these lesions in biological environments.

Tandem GA residues on opposite sides of the loop in molecular beacon-like DNA hairpins compact the loop and increase hairpin stability

The free solution electrophoretic mobilities and thermal stabilities of hairpins formed by two complementary 26-nucleotide oligomers have been measured by capillary electrophoresis. The oligomers are predicted to form molecular beacon-like hairpins with 5 bp stems and 16 nucleotides in the loop. One hairpin, called hairpin2 (hp2), migrates with a relatively fast free solution mobility and exhibits melting temperatures that are reasonably well predicted by the popular structure-prediction program Mfold. Its complement, called hairpin1 (hp1), migrates with a slower free solution mobility and forms a stable hairpin only in solutions containing ≥200 mM Na(+). The melting temperatures observed for hp1 are ~18 °C lower than those observed for hp2 and ~20 °C lower than those predicted by Mfold. The greater thermal stability of hp2 is due to the presence of tandem GA residues on opposite sides of the loop. If the corresponding TC residues in the hp1 loop are replaced by tandem GA residues, the melting temperatures of the modified hairpin are close to those observed for hp2. Eliminating the tandem GA residues in the hp2 loop significantly decreases the thermal stability of hp2. If the loops are replaced by a loop of 16 thymine residues, the free solution mobilities and thermal stabilities of the T-loop hairpin are equal to those observed for hp1. Hence, the loop of hp1 appears to be relatively unstructured, with few base-base stacking interactions. Interactions between tandem GA residues on opposite sides of the hp2 loop appear to compact the loop and increase hairpin stability.

Recognition of thymine in DNA bulges by a Zn(II) macrocyclic complex

A Zn(II) macrocyclic complex with appended quinoline is a bifunctional recognition agent that uses both the Zn(II) center and the pendent aromatic group to bind to thymine in bulges with good selectivity over DNA containing G, C or A bulges. Spectroscopic studies show that the stem containing the bulge stays largely intact in a DNA hairpin with the Zn(II) complex bound to the thymine bulge.

DNA interstrand cross-links of an antitumor trinuclear platinum(II) complex: thermodynamic analysis and chemical probing

The trinuclear platinum compound [{trans-PtCl(NH(3))(2)}(2)(μ-trans-Pt(NH(3))(2){NH(2)(CH(2))(6)NH(2)}(2))](4+) (BBR3464) belongs to the polynuclear class of platinum-based anticancer agents. These agents form in DNA long-range (Pt,Pt) interstrand cross-links, whose role in the antitumor effects of BBR3464 predominates. Our results show for the first time that the interstrand cross-links formed by BBR3464 between two guanine bases in opposite strands separated by two base pairs (1,4-interstrand cross-links) exist as two distinct conformers, which are not interconvertible, not only if these cross-links are formed in the 5'-5', but also in the less-usual 3'-3' direction. Analysis of the conformers by differential scanning calorimetry, chemical probes of DNA conformation, and minor groove binder Hoechst 33258 demonstrate that each of the four conformers affects DNA in a distinctly different way and adopts a different conformation. The results also support the thesis that the molecule of antitumor BBR3464 when forming DNA interstrand cross-links may adopt different global structures, including different configurations of the linker chain of BBR3464 in the minor groove of DNA. Our findings suggest that the multiple DNA interstrand cross-links available to BBR3464 may all contribute substantially to its cytotoxicity.

Unfolding thermodynamics of DNA pyrimidine triplexes with different molecularities

Nucleic acid oligonucleotides (ODNs), as drugs, present an exquisite selectivity and affinity that can be used in antigene and antisense strategies for the control of gene expression. In this work we try to answer the following question: How does the molecularity of a DNA triplex affect its overall stability and melting behavior? To this end, we used a combination of temperature-dependent UV spectroscopy and calorimetric (differential scanning calorimetry) techniques to investigate the melting behavior of DNA triplexes with a similar helical stem, TC+TC+TC+T/AGAGAGA/TCTCTCT, but formed with different strand molecularity. We determined standard thermodynamic profiles and the differential binding of protons and counterions accompanying their unfolding. The formation of a triplex is accompanied by a favorable free energy term, resulting from the typical compensation of favorable enthalpy-unfavorable entropy contributions, i.e., the folding of a particular triplex is enthalpy driven. The magnitude of the favorable enthalpy contributions corresponds to the number and strength of the base-triplet stacks formed, which are helped by stacking contributions due to the incorporation of dangling ends or loops. Triplex stability is in the following order: monomolecular > bimolecular > trimolecular; this is explained in terms of additional stacking contributions due to the inclusion of loops. As expected, acidic pH stabilized all triplexes by allowing protonation of the cytosines in the third strand; however, the percentage of protonation increases as the molecularity decreases. The results help to choose adequate solution conditions for the study of triplexes containing different ratios of CGC+ and TAT base triplets and to aid in the design of oligonucleotide sequences as targeting reagents that could effectively react with mRNA sequences involved in human diseases, thereby increasing the feasibility of using the antisense strategy for therapeutic purposes.

Monovalent ion dependence of neomycin B binding to an RNA aptamer characterized by spectroscopic methods

Understanding the interactions of small molecules like antibiotics with RNA is a prerequisite for the development of novel drugs. In this study we address structural and thermodynamic features of such interactions by using a simple model system: the binding of the highly charged antibiotic neomycin B to a short hairpin RNA molecule. Nucleotide A16, which acts as a flap over the neomycin B binding pocket, was substituted by the fluorescent adenine analogue 2-aminopurine (2-AP). Steady-state and time-resolved fluorescence measurements were complemented by UV-melting and circular dichroism studies. The binding of neomycin B at three sites was found to have a strong inverse correlation with Na(+) concentration. For the highest-affinity site, both fluorescence and UV absorption experiments were consistent with a model assuming at least three neomycin NH(3) (+) groups participating in addition to hydrogen bonds in electrostatic interactions with the RNA. The variation of fluorescence intensity and lifetime upon neomycin B binding indicated unstacking of 2-AP16 from neighbouring bases as it flipped over the binding pocket. RNA conformational changes upon binding of the antibiotic were confirmed by circular dichroism. The two weaker binding sites were characterized as unspecific binding to the aptamer, while the high-affinity binding event was shown to be highly specific even at high ionic concentration. In addition, 2-AP was confirmed to be a noninvasive fluorescent probe; it serves as a sensitive spectroscopic tool to investigate details of the interactions between small molecules and RNA.

Effects of monofunctional adducts of platinum(II) complexes on thermodynamic stability and energetics of DNA duplexes

Effects of adducts of [PtCl(NH3)3]Cl or chlorodiethylenetriamineplatinum(II) on DNA stability were studied with emphasis on thermodynamic origins of that stability. Oligodeoxyribonucleotide duplexes (15-bp) containing the single, site-specific monofunctional adduct at G-residues of the central sequences TGT/ACA or 5'-AGT/5'-ACT were prepared and analyzed by differential scanning calorimetry, temperature-dependent ultraviolet absorption and circular dichroism. The unfolding of the platinated duplexes was accompanied by relatively small unfavorable free energy terms. This destabilization was enthalpic in origin. On the other hand, a relatively large reduction of melting temperature (T(m)) was observed as a consequence of the monofunctional adduct in the TGT sequence, whereas T(m) due to the adduct in the AGT sequence was reduced only slightly. We also examined the efficiency of the mammalian nucleotide excision repair system to remove from DNA the monofunctional adducts and found that these lesions were not recognized by this repair system. Thus, rather thermodynamic than thermal characterization of DNA adducts of monofunctional platinum compounds is a property implicated in the modulation of downstream effects such as protein recognition and repair.

Thermodynamics of i-tetraplex formation in the nuclease hypersensitive element of human c-myc promoter

More than 85% of c-myc transcription is controlled by the nuclease hypersensitive element III(1) upstream of the P1 promoter of this oncogene. The purine-rich sequence in the anti-sense strand forms a G-quadruplex, which has been recently implicated in colorectal cancer, and is proposed as a silencer element [Proc. Natl. Acad. Sci. USA 101 (2004) 6140]. This prompted us to characterize the thermodynamics and proton/counterion effect of the complementary pyrimidine-rich sequence, which forms a C-tetraplex. We report the thermodynamic parameters for folding of the pyrimidine-rich DNA fragment from this region into a C-tetraplex. At 20 degrees C, we observed a DeltaG of -10.36+/-0.13kcalmol(-1) with favorable enthalpy (DeltaH=75.99+/-0.99kcalmol(-1)) and unfavorable entropy (TDeltaS=65.63+/-0.88 kcalmol(-1)) at pH 5.3 in 20mM NaCl for tetraplex folding. Similar characteristic stabilizing enthalpy and destabilizing entropy were observed at other pH and ionic strengths. Folding was induced by uptake of about two to three protons per mole of tetraplex while a marginal (0.5-1mol/mol) counterion uptake was observed. In the context of current understanding of c-myc transcription we envisage a role of the i-motif in remodeling the G-quadruplex silencer.

Effects of small neutral osmolytes on the supercoiling free energy and intrinsic twist of p30? DNA

Both theory and experiments are employed to investigate the effects of small neutral osmolytes on the average intrinsic twist (l0), the torsion and bending elastic constants, and the twist energy parameter (ET) that governs the supercoiling free energy. The experimental data for ethylene glycol and acetamide at 37 degrees C suggest, and are interpreted in terms of, a model wherein the DNA exhibits an equilibrium between two distinct conformational states that possess different numbers of bound water molecules and exhibit different intrinsic twists and torsion and bending elastic constants. Expressions are derived to relate the effective ET and l0 to the equilibrium constant, water activity (aw), and number (n) of bound water molecules released per cooperative domain undergoing the two-state transition. The variations of l0 and ET with -ln(aw) are similar for acetamide and ethylene glycol at 37 degrees C. Fitting the theory to those data yields the range n = 103-125 for ethylene glycol and n = 71-113 for acetamide, depending on the assumed value of ET for the dehydrated state. The cooperative domain size of the two-state transition is estimated to exceed about 25-30 base pairs (bp). Between 0 and 19.4 w/v % ethylene glycol, the torsion elastic constant, measured by time-resolved fluorescence polarization anisotropy (FPA), increases by 1.37-fold, whereas the measured ET decreases by 1.15-fold over that same range. The implied decrease in bending rigidity over that range is by a factor of about 0.7. The variations of l0 and ET with increasing -ln(aw) due to added ethylene glycol at 37 degrees C are far smaller than the corresponding variations observed previously at 14 and 15 degrees C. However, at 21 degrees C, upon adding either ethylene glycol or acetamide, l0 and ET initially decline steeply with increasing -ln(aw), with slopes possibly comparable to those seen at 14 and 15 degrees C, but then flatten out and follow curves similar to those at 37 degrees C. Possible origins of such mixed behavior are discussed. The effects of betaine at both 37 and 21 degrees C differ qualitatively and quantitatively in various respects from those of ethylene glycol and acetamide. Upon adding sucrose, l0 initially jumps to higher plateaus at both 37 and 21 degrees C, but its effects on ET cannot be reliably assessed, due to the limited range of -ln(aw).

Energetic contributions for the formation of TAT/TAT, TAT/CGC(+), and CGC(+)/CGC(+) base triplet stacks

We used a combination of spectroscopic and calorimetric techniques to determine complete thermodynamic profiles accompanying the folding of a set of triple helices and control duplexes. Specifically, we studied the sequences: d(A(7)C(5)T(7)C(5)T(7)), d(A(6)C(5)T(6)C(5)T(6)), d(A(6)C(5)T(6)), d(AGAGAGAC(5)TCTCTCTC(5)TCTCTCT), d(AGAGAC(5)TCTCTC(5)TCTCT), d(AGAGAC(5)TCTCTC(2)), d(AAGGAC(5)TCCTTC(5)TTCCT), d(AGGAAC(5)TTCCTC(5)TCCTT), and d(GAAAGC(5)CTTTCC(5)CTTTC). Circular dichroism spectroscopy indicated that all triplexes and duplexes are in the "B" conformation. DSC melting experiments revealed that the formation of triplexes is accompanied by a favorable free energy change, which arises from the compensation of a large and favorable enthalpic contribution with an unfavorable entropic contribution. Comparison of the thermodynamic profiles of these triplexes yielded enthalpic contributions of -24 kcal/mol, -23 kcal/mol, and -22 kcal/mol for the formation of TAT/TAT, TAT/CGC(+), and CGC(+)/CGC(+) base triplet stacks, respectively. UV melts as a function of sodium concentration show sodium ions stabilize the triplexes that contain only TAT triplets but destabilize the triplexes that contain CGC(+) triplets. UV melts as a function of pH indicate that the protonation of the third strand and loop cytosines stabilizes the triplexes that contain CGC(+) and TAT triplets, respectively. Our overall results suggest that the triplex to duplex transition of triplexes that contain CGC(+) triplets is accompanied by a release of protons and an uptake of sodium, while their duplex to random coil transition is accompanied by a release of sodium ions. A consequence of this opposite sodium dependence is that their coupled transitions are nearly independent of sodium concentration but are dependent on the experimental pH.

DNA oligonucleotide duplexes containing intramolecular platinated cross-links: energetics, hydration, sequence, and ionic effects

The anticancer activity of cisplatin arises from its ability to bind covalently to DNA, forming primarily intrastrand cross-links to adjacent purine residues; the most common adducts involve d(GpG) (65%) and d(ApG) (25%) intrastrand cross-links. The incorporation of these platinum adducts in a B-DNA helix induces local distortions, causing bending and unwinding of the DNA. In this work, we used temperature-dependent UV spectroscopy to investigate the unfolding thermodynamics, and associated ionic effects, of two sets of DNA decamer duplexes containing either cis-[Pt(NH(3))(2)[d(GpG]] or cis-[Pt(NH(3))(2) [d(ApG]] cross-links, and their corresponding unmodified duplexes. The platinated duplexes are less stable and unfold with lower T(M)s (and Delta G degrees s) in enthalpy-driven reactions, which indicates a loss of favorable base-pair stacking interactions. The folding thermodynamics and hydration effects for the first set of decamers containing the d(GpG) cross-link was investigated by a combination of titration calorimetry, density, and ultrasound techniques. The hydration parameters showed an uptake of structural water by the platinated duplex and a release of electrostricted water by the control duplex. Relative to the unmodified duplex, the folding of the platinated duplex at 20 degrees C yielded a positive Delta Delta G degrees term [and positive Delta Delta H-Delta(T Delta S) compensation] and a negative differential volume change. The opposite signs of the Delta Delta G degrees and Delta Delta V terms confirmed its uptake of structural water. Further, solvent-accessible surface areas calculations for a similar pair of dodecamer duplexes indicated that the modified duplex has a 503 oeA(2) higher polar and nonpolar surface area that is exposed to the solvent. Therefore, the incorporation of a platinum adduct in duplex DNA disrupts favorable base-pair stacking interactions, yielding a greater exposure of aromatic bases to the solvent, which in turn immobilizes structural water. The overall results correlate nicely with the results reported in the available structural data of nuclear magnetic resonance solution studies.

Thermodynamic and hydration effects for the incorporation of a cationic 3-aminopropyl chain into DNA

The introduction of cationic 5-(omega-aminoalkyl)-2'-deoxypyrimidines into duplex DNA has been shown to induce DNA bending. In order to understand the energetic and hydration contributions for the incorporation of a cationic side chain in DNA a combination of spectroscopy, calorimetry and density techniques were used. Specifically, the temperature unfolding and isothermal formation was studied for a pair of duplexes with sequence d(CGTAGUCG TGC)/d(GCACGACTACG), where U represents 2'-deoxyuridine ('control') or 5-(3-aminopropyl)-2'-deoxyuridine ('modified'). Continuous variation experiments confirmed 1:1 stoichiometries for each duplex and the circular dichroism spectra show that both duplexes adopted the B conformation. UV and differential scanning calorimetry melting experiments reveal that each duplex unfolds in two-state transitions. In low salt buffer, the 'modified' duplex is more stable and unfolds with a lower endothermic heat and lower release of counterion and water. This electrostatic stabilization is entropy driven and disappears at higher salt concentrations. Complete thermodynamic profiles at 15 degrees C show that the favorable formation of each duplex results from the compensation of a favorable exothermic heat with an unfavorable entropy contribution. However, the isothermal profiles yielded a differential enthalpy of 8.8 kcal/mol, which is 4.3 kcal/mol higher than the differential enthalpy observed in the unfolding profiles. This indicates that the presence of the aminopropyl chain induces an increase in base stacking interactions in the modified single strand and a decrease in base stacking interactions in the modified duplex. Furthermore, the formation of the 'control' duplex releases water while the 'modified' duplex takes up water. Relative to the control duplex, formation of the modified duplex at 15 degrees C yielded a marginal differential DeltaG degrees term, positive DeltaDeltaH(ITC)-Delta(TDeltaS) compensation, negative DeltaDeltaV and a net release of counterions. The opposite signs of the differential enthalpy-entropy compensation and differential volume change terms show a net uptake of structural water around polar and non-polar groups. This indicates that incorporation of the aminopropyl chain induces a higher exposure of aromatic bases to the solvent, which may be consistent with a small and local bend in the 'modified' duplex.

Incorporation of a cationic aminopropyl chain in DNA hairpins: thermodynamics and hydration

We report on the physicochemical effects resulting from incorporating a 5-(3-aminopropyl) side chain onto a 2'-deoxyuridine (dU) residue in a short DNA hairpin. A combination of spectroscopy, calorimetry, density and ultrasound techniques were used to investigate both the helix-coil transition of a set of hairpins with the following sequence: d(GCGACTTTTTGNCGC) [N = dU, deoxythymidine (dT) or 5-(3-aminopropyl)-2'-deoxyuridine (dU*)], and the interaction of each hairpin with Mg(2+). All three molecules undergo two-state transitions with melting temperatures (T(M)) independent of strand concentration that indicates their intramolecular hairpin formation. The unfolding of each hairpin takes place with similar T(M) values of 64-66 degrees C and similar thermodynamic profiles. The unfavorable unfolding free energies of 6.4-6.9 kcal/mol result from the typical compensation of unfavorable enthalpies, 36-39 kcal/mol, and favorable entropies of approximately 110 cal/mol. Furthermore, the stability of each hairpin increases as the salt concentration increases, the T(M)-dependence on salt yielded slopes of 2.3-2.9 degrees C, which correspond to counterion releases of 0.53 (dU and dT) and 0.44 (dU*) moles of Na(+) per mole of hairpin. Absolute volumetric and compressibility measurements reveal that all three hairpins have similar hydration levels. The electrostatic interaction of Mg(2+) with each hairpin yielded binding affinities in the order: dU > dT > dU*, and a similar release of 2-4 electrostricted water molecules. The main result is that the incorporation of the cationic 3-aminopropyl side chain in the major groove of the hairpin stem neutralizes some local negative charges yielding a hairpin molecule with lower charge density.

A thermodynamic study on the formation and stability of DNA duplex at transcription site for DNA binding proteins GCN4

Using isothermal titration calorimetry (ITC), we studied the thermodynamic parameters of the 15-mer duplex dsDNA [d(GAGATGACTCATCTC)].[d(GAGATGAGTCATCTC)] formation from its two complementary single strands (S1 and S2) over a range of temperatures. The two complementary single strands d(GAGATGACTCATCTC) (herein called S1) and d(GAGATGAGTCATCTC) (herein called S2) containing palindromic sequences may assume ordered structures at low temperatures, which made the duplex dsDNA formation rather complicated. The thermodynamic parameters for the duplex formation, such as the binding constants (Kb), the enthalpies (delta H0), the free energies (delta G0), the entropies (delta S0) are strongly temperature-dependent. The thermally-induced disruptions of the duplex and its two complementary single strands, S1 and S2, were measured using differential scanning calorimetry (DSC) and CD spectroscopy, the results demonstrate that the DNA duplex is very stable, and its component single strands have an ordered structure at low temperature. This 15-mer specific sequence DNA may act as recognition site for DNA binding proteins GCN4 and plays a key role in transcription regulation of gene expression. Our analyses of the thermodynamic data suggest that the duplex formation is a coupled process between conformational transitions in the two single strands and their binding to form duplex dsDNA.

Thermodynamic investigation of the association of ethidium, propidium and bis-ethidium to DNA hairpins

We have used a combination of calorimetric and spectroscopic techniques to investigate the association of the bis-intercalator ethidium homodimer (bis-ethidium) to short DNA hairpins with sequences: d(GCGCT5GCGC) and d(CGCGT5CGCG). The helix-coil transition of each hairpin, investigated by UV and calorimetric melting protocol, takes place in monomolecular two-state transitions with characteristic enthalpies of approximately 37 kcal mol-1 for disrupting the four dG-dC base pairs of the hairpin stems. Deconvolution of the bis-ethidium-hairpin calorimetric titration curves indicate that each hairpin contains two distinct binding sites for the ligand: a high affinity site in the stem (Kb approximately 10(7)) that accommodates one bis-ethidium molecule and a lower affinity site (Kb approximately 10(6)) located probably at the loop that accommodates two bis-ethidium molecules. The overall stoichiometries of three ligands per hairpin are in agreement with those obtained in continuous variation experiments using visible spectroscopy. The interaction of bis-ethidium for each type of sites results in enthalpy driven reactions, with average binding enthalpies, delta Hb, of -13.1 and -12.1 kcal mol-1 for the stem and loop sites, respectively. Comparison to the thermodynamic profiles of ethidium and propidium binding reveals that the bis-ethidium binding to the stem site of each hairpin has a more favorable free energy term of -1.4 kcal mol-1 and more favorable enthalpy of -4.2 kcal mol-1. These suggest that only one phenanthridine ring of bis-ethidium intercalates in the stem, while the second planar ring is exposed to solvent or weakly associated to the surface of DNA.

Thermodynamics of DNA hairpins: contribution of loop size to hairpin stability and ethidium binding

A combination of calorimetric and spectroscopic techniques was used to evaluate the thermodynamic behavior of a set of DNA hairpins with the sequence d(GCGCTnGCGC), where n = 3, 5 and 7, and the interaction of each hairpin with ethidium. All three hairpins melt in two-state monomolecular transitions, with tm's ranging from 79.1 degrees C (T3) to 57.5 degrees C (T7), and transition enthalpies of approximately 38.5 kcal mol-1. Standard thermodynamic profiles at 20 degrees C reveal that the lower stability of the T5 and T7 hairpins corresponds to a delta G degree term of +0.5 kcal mol-1 per thymine residue, due to the entropic ordering of the thymine loops and uptake of counterions. Deconvolution of the ethidium-hairpin calorimetric titration curves indicate two sets of binding sites that correspond to one ligand in the stem with binding affinity, Kb, of approximately 1.8 x 10(6) M-1, and two ligands in the loops with Kb of approximately 4.3 x 10(4) M-1. However, the binding enthalpy, delta Hb, ranges from -8.6 (T3) to -11.6 kcal mol-1 (T7) for the stem site, and -6.6 (T3) to -12.7 kcal mol-1 (T7) for the loop site. Relative to the T3 hairpin, we obtained an overall thermodynamic contribution (per dT residue) of delta delta Hb = delta(T delta Sb) = -0.7(5) kcal mol-1 for the stem sites and delta delta Hb = delta(T delta Sb) = -1.5 kcal mol-1 for the loop sites. Therefore, the induced structural perturbations of ethidium binding results in a differential compensation of favorable stacking interactions with the unfavorable ordering of the ligands.

Sampling of the conformations of the d(CGCTGCGGC) hairpin in solution by two-dimensional nuclear magnetic resonance and theoretical methods

Most NMR studies of DNA oligomers have focused on rigid structures that show a strong preference for one or a small set of ground-state conformations. There is an increasing interest in extending NMR methods to investigate DNA systems in which this preference does not exist. A DNA hairpin is one such system where a large number of low-energy structures coexist in solution. In this article we show how 1D/2D NMR data of the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin are used to map the conformational space of this molecule. First, we characterize the gross morphology of the hairpin by monitoring the exchangeable imino signals in the molecule. Second, we extract a set of inter-proton distances (i.e., the average values and the associated dispersions) for various pairwise interactions by performing full-matrix NOESY simulation with respect to the observed NOESY data for mixing times of 250 and 100 ms. Third, we use these distances as structural constraints to perform a 300-ps molecular dynamics simulation at 500 K. Fourth, we extract 600 snapshots (one after every 0.5 ps) from the MD trajectory and perform constrained energy minimization to map local minima on the sampled energy surface (we call this the rapid temperature quenching step). Fifth, we assign 600 structures to 14 disjoint clusters such that conformationally similar hairpins belong to the same cluster while conformationally distinct hairpins belong to different clusters. Finally, we interpret the NOESY data in terms of conformationally distinct structures by recalculating NOESY contributions taken from representative structures of different clusters. Our analyses clearly demonstrate that the NMR data correspond to an ensemble of distinct structures, i.e., a set of energetically stable but conformationally distinct structures that satisfies the constraints of loop folding in the d(C1-G2-C3-T4-G5-C6-G7-G8-C9) hairpin. Two types of loop folding consistent with NMR data are obtained: (i) a hairpin with two G.C pairs in the stem and four residues in the loop and (ii) a hairpin with two G.C pairs and a reverse wobble G.T pair in the stem plus two residues in the loop.

Theoretical predictions of DNA hairpin loop conformations: correlations with thermodynamic and spectroscopic data

A computational procedure for generating conformations of DNA hairpin loop structures from a broad range of low-energy starting states is described. The starting point of the modeling is the distribution of oligonucleotide chain conformations obtained from Monte Carlo simulations of feasible dinucleotide steps. Structures which meet the spatial criteria for hairpin loop formation are selected from the distributions and subsequently minimized using all-atom molecular mechanics. Both d(CTnG) and d(CAnG) oligomers, where n = 3, 4, or 5, are modeled. These sequences are chosen because of the large number of published NMR and thermodynamic studies on DNA hairpins containing thymine or adenine residues. The minimized three-dimensional hairpin loop structures are compared with one another as well as analyzed in terms of available experimental data. The computational approach provides the first detailed analysis of DNA hairpin loop structure in terms of a multistate conformational model. Investigation of the minimized conformations reveals several interesting structural features. First, hairpin loops of the same sequence adopt several distinctly different conformations, as opposed to minor variants of the same equilibrium structure, that could potentially interconvert in solution. Second, in contrast to double-helical nucleic acids, the hairpin loop models exhibit hydrophobic and hydrophilic surfaces. The different disposition of hydrophobic groups in loops versus duplexes could modulate both protein-nucleic acid interactions and nucleic acid self-associations. Third, perpendicular aromatic interactions of loop residues are observed in many of the computed hairpins. This sort of interaction might be important in the stabilization of non-hydrogen-bonded nucleic acid secondary and tertiary structures. The predicted structural features in the models help, in addition, to account for the unusual thermodynamic properties of DNA hairpin loops. Comparison of the theoretically-generated NOEs in different structures further reveals that very different molecular structures and interactions can, in principle, produce the same NOEs. The multistate description suggested by this observation differs from the conventional interpretation of DNA solution structure in terms of the fluctuations about a single preferred chain conformation. There is not necessarily only one set of closely related structures consistent with the observed data.

Conformational studies of hairpin sequences from the ColE1 cruciform

Inverted repeat sequences derived from the ColE1 cruciform were investigated by nuclear magnetic resonance (NMR) and UV spectroscopy. It was shown that 15 different sequences exist as stable hairpin structures over a range of buffer conditions and DNA concentrations. Experiments with six oligomers (1-6) containing the native stem sequence and five base loops, found that the two hairpins with the wild-type loops (1-2) served as upper and lower bounds for the thermodynamic stability of all the other sequences. NMR experiments, including rotational correlation time measurements and NOESY spectra, were then performed on 1, the most stable hairpin sequence to begin to uncover a structural basis of its stability.