Translational Entropy and DNA Duplex Stability

Investigation of folding/unfolding DNA duplexes of various size and composition by superprecise calorimetry has revised several long-held beliefs concerning the forces responsible for the formation of the double helix. It was established that: 1) the enthalpy and the entropy of duplex unfolding are temperature dependent, increasing with temperature rise and having the same heat capacity increment for CG and AT pairs; 2) the enthalpy of AT melting is greater than that of the CG pair, so the stabilizing effect of the CG pair in comparison with AT results not from its larger enthalpic contribution (as expected from its extra hydrogen bond), but from the larger entropic contribution of the AT pair that results from its ability to fix ordered water in the minor groove and release it upon duplex unfolding; 3) the translation entropy, resulting from the appearance of a new kinetic unit on duplex dissociation, determines the dependence of duplex stability on its length and its concentration (it is an order-of-magnitude smaller than predicted from the statistical mechanics of gases and is fully expressed by the stoichiometric correction term); 4) changes in duplex stability on reshuffling the sequence (the “nearest-neighbor effect”) result from the immobilized water molecules fixed by AT pairs in the minor groove; and 5) the evaluated thermodynamic components permit a quantitative expression of DNA duplex stability.     

Differential hydration of dA.dT base pairing and dA and dT bulges in deoxyoligonucleotides

The role of water in the formation of stable duplexes of nucleic acids is being studied by determining the concurrent volume change, heats, and counterion uptake that accompany the duplexation process. The variability of the volume contraction that we have observed in the formation of a variety of homoduplexes suggests that sequence and conformation acutely affect the degree of hydration. We have used a combination of densimetric and calorimetric techniques to measure the change in volume and enthalpy resulting from the mixing of two complementary strands to form (a) fully paired duplexes with 10 or 11 base pairs and (b) bulged decameric duplexes with an extra dA or dT unmatched residue. We also monitored absorbance vs temperature profiles as a function of strand and salt concentration for all four duplexes. Relative to the decamer duplex, insertion of an extra dA.dT base pair to form an undecamer duplex results in a favorable enthalpy of -5.6 kcal/mol that is nearly compensated by an unfavorable entropy term of -5.1 kcal/mol. This enthalpy difference correlates with a differential uptake of water molecules, corresponding to an additional hydration of 16 mol of water molecules/mol of base pair. Relative to the fully paired duplexes, both bulged duplexes are 12-16 degrees C less stable and exhibit marginally larger counterion uptake on forming the duplex. The enthalpy change is slightly lower for the T-bulge duplex and less still for the A-bulge duplex. The volume change results indicate that an unmatched residue increases the amount of coulombic and/or structural hydration. The combined results strongly suggest that the destabilizing forces in bulged duplexes are partially compensated by an increase in hydration levels.     

Differential hydration of homopurine sequences relative to alternating purine/pyrimidine sequences.

The minor groove ligand distamycin A has been used to probe the relative hydration of the minor groove of eight synthetic polynucleotides of known sequence and composition. A combination of densimetric, calorimetric, and temperature-dependent spectroscopic techniques have been used to obtain complete thermodynamic profiles (delta Gzero, delta Hzero, delta Szero, and delta Vzero) for the association of distamycin A to all polymer duplexes. In 10 mM phosphate buffer, pH 7, binding of the drug to each of the polymeric duplexes resulted in characteristic negative changes in both the volume and enthalpy. Although the binding constants were found to be identical for pairs of isomer polynucleotides having identical compositions but different sequences, the values of delta Hzero, delta Szero, and delta Vzero of each such pair were remarkably different. The entropy changes were found to roughly parallel the volume changes; no such trend was seen between delta Hzero and delta Vzero. The data support the hypothesis that the volume changes observed for these systems reflect the coulombic-hydration contribution to the entropy. The heteropolymer duplexes generated much larger exothermic contributions, less favorable entropies and larger volume contractions than did the corresponding homopolymer duplexes of identical composition, and strongly suggest that polynucleotides with homopurine sequences are more hydrated than polynucleotides with alternating purine/pyrimidine sequences. In addition, it was found that duplexes containing guanine sharply reduced the affinity for the drug, also lowering the exothermicity but raising the entropy. This may be explained by the presence of an amino group in the minor groove that prevents hydrogen bonding. Substitution of the guanine with inosine reversed this trend in the thermodynamic properties. Furthermore, substitution of poly(dA) for poly(rA) in a duplex produced a similar reduction in the affinity, while raising the exothermic contribution and greatly reducing the favorable entropy effect in agreement with an apparent increase in the hydration state.     

The energetics of HMG box interactions with DNA: thermodynamic description of the target DNA duplexes

The thermal properties and energetics of formation of 10, 12 and 16 bp DNA duplexes, specifically interacting with the HMG box of Sox-5, have been studied by isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). DSC studies show that the partial heat capacity of these short duplexes increases considerably prior to the cooperative process of strand separation. Direct extrapolation of the pre and post-transition heat capacity functions into the cooperative transition zone suggests that unfolding/dissociation of strands results in no apparent heat capacity increment. In contrast, ITC measurements show that the negative enthalpy of complementary strand association increases in magnitude with temperature rise, implying that strand association proceeds with significant decrease of heat capacity. Furthermore, the ITC-measured enthalpy of strand association is significantly smaller in magnitude than the enthalpy of cooperative unfolding measured by DSC. To resolve this paradox, the heat effects upon heating and cooling of the separate DNA strands have been measured by DSC. This showed that cooling of the strands from 100 degrees C to -10 degrees C proceeds with significant heat release associated with the formation of intra and inter-molecular interactions. When the enthalpy of residual structure in the strands and the temperature dependence of the heat capacity of the duplexes and of their unfolded strands have been taken into account, the ITC and DSC results are brought into agreement. The analysis shows that the considerable increase in heat capacity of the duplexes with temperature rise is due to increasing fluctuations of their structure (e.g. end fraying and twisting) and this effect obscures the heat capacity increment resulting from the cooperative separation of strands, which in fact amounts to 200(+/-40) JK(-1) (mol bp)(-1). Using this heat capacity increment, the averaged standard enthalpy, entropy and Gibbs energy of formation of fully folded duplexes from fully unfolded strands have been determined at 25 degrees C as -33(+/-2) kJ (mol bp)(-1), -93(+/-4) J K(-1) (mol bp)(-1) and -5.0(+/-0.5) kJ (mol bp)(-1), respectively.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

The energetics of HMG box interactions with DNA: thermodynamics of the DNA binding of the HMG box from mouse sox-5

The energetics of the Sox-5 HMG box interaction with DNA duplexes, containing the recognition sequence AACAAT, were studied by fluorescence spectroscopy, isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). Fluorescence titration showed that the association constant of this HMG box with the duplexes is of the order 4x10(7) M(-1), increasing somewhat with temperature rise, i.e. the Gibbs energy is -40 kJ mol(-1) at 5 degrees C, decreasing to -48 kJ mol(-1) at 32 degrees C. ITC measurements of the enthalpy of association over this temperature range showed an endothermic effect below 17 degrees C and an exothermic effect above, suggesting a heat capacity change on binding of about -4 kJ K(-1) mol(-1), a value twice larger than expected from structural considerations. A straightforward interpretation of ITC data in heat capacity terms assumes, however, that the heat capacities of all participants in the association reaction do not change over the considered temperature range. Our previous studies showed that over the temperature range of the ITC experiments the HMG box of Sox-5 starts to unfold, absorbing heat and the heat capacities of the DNA duplexes also increase significantly. These heat capacity effects differ from that of the DNA/Sox-5 complex. Correcting the ITC measured binding enthalpies for the heat capacity changes of the components and complex yielded the net enthalpies which exhibit a temperature dependence of about -2 kJ K(-1) mol(-1), in good agreement with that predicted on the basis of dehydration of the protein-DNA interface. Using the derived heat capacity change and the enthalpy and Gibbs energy of association measured at 5 degrees C, the net enthalpy and entropy of association of the fully folded HMG box with the target DNA duplexes was determined over a broad temperature range. These functions were compared with those for other known cases of sequence specific DNA/protein association. It appears that the enthalpy and entropy of association of minor groove binding proteins are more positive than for proteins binding in the major groove. The observed thermodynamic characteristics of protein binding to the A+T-rich minor groove of DNA might result from dehydration of both polar and non-polar groups at the interface and release of counterions. The expected entropy of dehydration was calculated and found to be too large to be compensated by the negative entropy of reduction of translational/rotational freedom. This implies that DNA/HMG box association proceeds with significant decrease of conformational entropy, i.e. reduction in conformational mobility.     

Temperature dependence of thermodynamic properties for DNA/DNA and RNA/DNA duplex formation.

A clear difference in the enthalpy changes derived from spectroscopic and calorimetric measurements has recently been shown. The exact interpretation of this deviation varied from study to study, but it was generally attributed to the non-two-state transition and heat capacity change. Although the temperature-dependent thermodynamics of the duplex formation was often implied, systemic and extensive studies have been lacking in universally assigning the appropriate thermodynamic parameter sets. In the present study, the 24 DNA/DNA and 41 RNA/DNA oligonucleotide duplexes, designed to avoid the formation of hairpin or slipped duplex structures and to limit the base pair length less than 12 bp, were selected to evaluate the heat capacity changes and temperature-dependent thermodynamic properties of duplex formation. Direct comparison reveals that the temperature-independent thermodynamic parameters could provide a reasonable approximation only when the temperature of interest has a small deviation from the mean melting temperature over the experimental range. The heat capacity changes depend on the base composition and sequences and are generally limited in the range of -160 to approximately -40 cal.mol-1.K-1 per base pair. In contrast to the enthalpy and entropy changes, the free energy change and melting temperature are relatively insensitive to the heat capacity change. Finally, the 16 NN-model free energy parameters and one helix initiation at physiological temperature were extracted from the temperature-dependent thermodynamic data of the 41 RNA/DNA hybrids.     

Prediction of hybridization and melting for double-stranded nucleic acids

This article presents a general statistical mechanical approach to describe self-folding together with the hybridization between a pair of finite length DNA or RNA molecules. The model takes into account the entire ensemble of single- and double-stranded species in solution and their mole fractions at different temperatures. The folding and hybridization models deal with matched pairs, mismatches, symmetric and asymmetric interior loops, bulges, and single-base stacking that might exist at duplex ends or at the ends of helices. All possible conformations of the single- and double-stranded species are explored. Only intermolecular basepairs are considered in duplexes at this stage.In particular we focus on the role of stacking between neighboring nucleotide residues of single unfolded strands as an important source of enthalpy change on helix formation which has not been modeled computationally thus far. Changes in the states of the single strands with temperature are shown to lead to a larger heat effect at higher temperature. An important consequence of this is that predictions of enthalpies, which are based on databases of nearest-neighbor energy parameters determined for molecules or duplexes with lower melting temperatures compared with the melting temperatures of the oligos for which they are used as a predictive tool, will be underestimated.     

Oligonucleotide conjugates with minor groove ligands as probes for hybridization microarray chips

A possibility of using oligonucleotide conjugates with minor groove ligands as probes for hybridization microarray chips was studied. The oligonucleotide conjugates contain a hairpin ligand (MGB) composed of two tripyrrolcarboxamide residues with an aminocaproic acid residue as a linker and bound to the oligonucleotide duplex AT tract in a site-specific manner. We used as (5′-3′)-probes: GACAAGAp, GACAAAAp, GACAAGA-MGB, and GACAAAA-MGB. The oligonucleotides labeled with the Cy3 cyanine dye, Cy3-ACTAATTTTGTC and Cy3-ACTAATCTTGTC, were used as targets. The maximal MGB effect on the fluorescence level of microarray chip spots, which caused its fourfold increase as compared with the initial unmodified duplex, was observed for the duplex containing only AT pairs in the ligand binding site. The presence of AC and GT mutations in the binding site (imperfect duplexes) or a CG pair (perfect duplex) affect the change in fluorescence level to a considerably lesser degree.     

Solution structure and stability of the DNA undecamer duplexes containing oxanine mismatch

Solution structures of DNA duplexes containing oxanine (Oxa, O) opposite a cytosine (O:C duplex) and opposite a thymine (O:T duplex) have been solved by the combined use of (1)H NMR and restrained molecular dynamics calculation. One mismatch pair was introduced into the center of the 11-mer duplex of [d(GTGACO(6)CACTG)/d(CAGTGX(17)GTCAC), X = C or T]. (1)H NMR chemical shifts and nuclear Overhauser enhancement (NOE) intensities indicate that both the duplexes adopt an overall right-handed B-type conformation. Exchangeable resonances of C(17) 4-amino proton of the O:C duplex and of T(17) imino proton of O:T duplex showed unusual chemical shifts, and disappeared with temperature increasing up to 30 °C, although the melting temperatures were >50 °C. The O:C mismatch takes a wobble geometry with positive shear parameter where the Oxa ring shifted toward the major groove and the paired C(17) toward the minor groove, while, in the O:T mismatch pair with the negative shear, the Oxa ring slightly shifted toward the minor groove and the paired T(17) toward the major groove. The Oxa mismatch pairs can be wobbled largely because of no hydrogen bond to the O1 position of the Oxa base, and may occupy positions in the strands that optimize the stacking with adjacent bases.     

C-C-A-G-G-C-m5C-T-G-G. Helical fine structure, hydration, and comparison with C-C-A-G-G-C-C-T-G-G.

The self-complementary DNA duplex C-C-A-G-G-C-m5C-T-G-G has been refined against 1.75-A x-ray diffraction data to an R value of 17.4%. In the crystal of space group P6, 10-base pair DNA fragments with characteristic sequence-related fine structure stack end to end to form long antiparallel B-type double helices. As shown by a structure analysis at lower resolution (Heinemann, U., and Alings, C. (1991) EMBO J. 10, 35-43), the overall geometry of C-C-A-G-G-C-m5C-T-G-G is similar to that of the unmethylated analog C-C-A-G-G-C-C-T-G-G despite a different crystal environment. The present high resolution structure analysis permits a detailed comparison of the two duplexes and their hydration spheres. Helical parameters are significantly correlated between both molecules, with the exception of the base pair propeller. Sugar pucker and backbone torsion angles alpha, gamma, delta, and chi show similar mean values, but their individual values deviate significantly between duplexes. In contrast, torsion angles beta, epsilon, and zeta change along the strands of both duplexes in much the same way. The effect of single-site methylation on DNA conformation appears to be small and limited to the base pairs directly involved. Methylation tends to push base pairs toward the minor groove of the helix. A regular minor groove hydration pattern involves dual hydrogen bonding of water molecules to O-4' and base atoms of C-C-A-G-G-C-m5C-T-G-G.     

Thermodynamics of DNA branching.

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.     

Importance of minor groove functional groups for the stability of DNA duplexes

Eight oligonucleotide duplexes have been prepared with four pairs of selected complementary pairs of native/analogue heterocyclic bases incorporated at a selected test site. The base pairs vary in the nature of their functionality in the minor groove. Each pair has a minor groove purine amino group present or absent, and correspondingly has a minor grove pyrimidine carbonyl present or absent. Loss of duplex stability is most notable when the minor groove pyrimidine carbonyl is absent although in other respects normal Watson-Crick hydrogen bonding is maintained in these sequences. These differences in stability are discussed in terms of possible variations in minor groove hydration.     

Crystal structure of the highly distorted chimeric decamer r(C)d(CGGCGCCG)r(G).spermine complex--spermine binding to phosphate only and minor groove tertiary base-pairing.

The crystal structure of the self-complementary chimeric decamer duplex r(C)d(CGGCGCCG)r(G), with RNA base pairs at both termini, has been solved at 1.9 A resolution by the molecular replacement method and refined to an R value of 0.145 for 2,314 reflections. The C3'-endo sugar puckers of the terminal riboses apparently drive the entire chimeric duplex into an A-DNA conformation, in contrast to the B-DNA conformation adopted by the all-deoxy decamer of the same sequence. Five symmetry related duplexes encapsulate a spermine molecule which interacts with ten phosphate groups, both directly and through water molecules to form multiple ionic and hydrogen bonding interactions. The spermine interaction severely bends the duplexes by 31 degrees into the major groove at the fourth base pair G(4).C(17), jolts it and slides the 'base plate' into the minor groove. This base pair, together with the adjacent base pair in the top half and the corresponding pseudo two-fold related base pairs in the bottom half, form four minor groove base-paired multiples with the terminal base pairs of two neighboring duplexes.     

Enthalpy and heat capacity changes for formation of an oligomeric DNA duplex: interpretation in terms of coupled processes of formation and association of single-stranded helices.

The thermodynamics of self-assembly of a 14 base pair DNA double helix from complementary strands have been investigated by titration (ITC) and differential scanning (DSC) calorimetry, in conjunction with van't Hoff analysis of UV thermal scans of individual strands. These studies demonstrate that thermodynamic characterization of the temperature-dependent contributions of coupled conformational equilibria in the individual "denatured" strands and in the duplex is essential to understand the origins of duplex stability and to derive stability prediction schemes of general applicability. ITC studies of strand association at 293 K and 120 mM Na+ yield an enthalpy change of -73 +/- 2 kcal (mol of duplex)-1. ITC studies between 282 and 312 K at 20, 50, and 120 mM Na+ show that the enthalpy of duplex formation is only weakly salt concentration-dependent but is very strongly temperature-dependent, decreasing approximately linearly with increasing temperature with a heat capacity change (282-312 K) of -1.3 +/- 0.1 kcal K-1 (mol of duplex)-1. From DSC denaturation studies in 120 mM Na+, we obtain an enthalpy of duplex formation of -120 +/- 5 kcal (mol of duplex)-1 and an estimate of the corresponding heat capacity change of -0.8 +/- 0.4 kcal K-1 (mol of duplex)-1 at the Tm of 339 K. van't Hoff analysis of UV thermal scans on the individual strands indicates that single helix formation is noncooperative with a temperature-independent enthalpy change of -5.5 +/- 0.5 kcal at 120 mM Na+. From these observed enthalpy and heat capacity changes, we obtain the corresponding thermodynamic quantities for two fundamental processes: (i) formation of single helices from disordered strands, involving only intrastrand (vertical) interactions between neighboring bases; and (ii) formation of double helices by association (docking) of single helical strands, involving interstrand (horizontal and vertical) interactions. At 293 K and 120 mM Na+, we calculate that the enthalpy change for association of single helical strands is approximately -64 kcal (mol of duplex)-1 as compared to -210 kcal (mol of duplex)-1 calculated for duplex formation from completely unstructured single strands and to the experimental ITC value of -73 kcal (mol of duplex)-1. The intrinsic heat capacity change for association of single helical strands to form the duplex is found to be small and positive [ approximately 0.1 kcal K-1 (mol of duplex)-1], in agreement with the result of a surface area analysis, which also predicts an undetectably small heat capacity change for single helix formation.     

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

The introduction of cationic 5-(ω-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°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°C yielded a marginal differential ΔG° term, positive ΔΔH_ITC–Δ(TΔS) compensation, negative ΔΔV 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.     

Interaction between the Z-type DNA duplex and 1,3-propanediamine: crystal structure of d(CACGTG)2 at 1.2 A resolution

The crystal structure of a hexamer duplex d(CACGTG)(2) has been determined and refined to an R-factor of 18.3% using X-ray data up to 1.2 A resolution. The sequence crystallizes as a left-handed Z-form double helix with Watson-Crick base pairing. There is one hexamer duplex, a spermine molecule, 71 water molecules, and an unexpected diamine (Z-5, 1,3-propanediamine, C(3)H(10)N(2)) in the asymmetric unit. This is the high-resolution non-disordered structure of a Z-DNA hexamer containing two AT base pairs in the interior of a duplex with no modifications such as bromination or methylation on cytosine bases. This structure does not possess multivalent cations such as cobalt hexaammine that are known to stabilize Z-DNA. The overall duplex structure and its crystal interactions are similar to those of the pure-spermine form of the d(CGCGCG)(2) structure. The spine of hydration in the minor groove is intact except in the vicinity of the T5A8 base pair. The binding of the Z-5 molecule in the minor grove of the d(CACGTG)(2) duplex appears to have a profound effect in conferring stability to a Z-DNA conformation via electrostatic complementarity and hydrogen bonding interactions. The successive base stacking geometry in d(CACGTG)(2) is similar to the corresponding steps in d(CG)(3). These results suggest that specific polyamines such as Z-5 could serve as powerful inducers of Z-type conformation in unmodified DNA sequences with AT base pairs. This structure provides a molecular basis for stabilizing AT base pairs incorporated into an alternating d(CG) sequence.     

Probing the hydration of the minor groove of A.T synthetic DNA polymers by volume and heat changes

The minor-groove ligand netropsin provides a sensitive probe of the hydration difference between poly(dA).poly(dT) and poly[d(AT)].poly[d(AT)]. We have measured the volume change delta V accompanying binding of netropsin to these polymers, using an improved magnetic suspension densimeter. For poly(dA).poly(dT) we find delta V = +97 mL/mol of bound netropsin at pH 7.0 and 10 mM sodium phosphate buffer. For poly[d(AT)].poly[d(AT)] we find delta V = -16 mL/mol of bound netropsin. This striking differential effect suggests that the poly(dA).poly(dT) duplex compresses more water (or is more extensively hydrated). From our enthalpy and entropy results we estimate the approximately 10 water molecules, immobilized in the minor groove of this system, are displaced by each netropsin bound. The volume increase, however, is substantially larger than can be explained by a simple melting of these immobilized water molecules in the minor groove. A decompression of at least 40 water molecules must attend the complexation to the poly(dA).poly(dT) duplex. This suggests that the conformation change attending the binding of the drug to this polymer duplex causes a further dehydration, whereas no such change in dehydration and configuration for the heteropolymer system is indicated.     

Hydration pattern of A4T4 and T4A4 DNA: a molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.     

Hydration of B-DNA: comparison between the water network around poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) on the basis of Monte Carlo computations

A computational method is elaborated for studying the water environment around regular polynucleotide duplexes; it allows rigorous structural information on the hydration shell of DNA to be obtained. The crucial aspect of this Monte Carlo simulation is the use of periodical boundary conditions. The output data consists of local maxima of water density in the space near the DNA molecule and the properties of one- and two-membered water bridges as function of pairs of polar groups of DNA. In the present paper the results for poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) are presented. The differences in their hydration shells are of a purely structural nature and are caused by the symmetry of the polar groups of the polymers under study, the symmetry being reflected by the hydration shell. The homopolymer duplex hydration shell mirrors the mononucleotide repeat. The water molecules contacting the polynucleotide in the minor groove are located nearly in the plane midway between the planes of successive base pairs. One water molecule per base pair forms a water bridge facing two polar groups of bases from adjacent base pairs and on different strands making a "spine"-like structure. In contrast, the major groove hydration is stabilized exclusively by two-membered water bridges; the water molecules deepest in the groove are concentrated near the plane of the corresponding base pair. The alternating polymer is characterized by a marked dyad symmetry of the hydration shell corresponding to the axis between two successive base pairs. The minor groove hydration of the dCpdG step resembles the characteristic features of the homopolymer, but the bridge between the O2 oxygens of the other base-stacking type is formed by two water molecules. The major groove hydration is characterized by high probability of one-membered water bridges and by localization of a water molecule on the dyad axis of the dGpdC step. The found structural elements are discussed as reasonable invariants of a dynamic hydration shell.