DNA Duplex Length and Salt Concentration Dependence of Enthalpy−Entropy Compensation Parameters for DNA Melting

Replacing salt correction factors with optimized RNA nearest-neighbour enthalpy and entropy parameters

We calculate the nearest-neighbour enthalpies and entropies at 5 salt concentrations of 18 RNA sequences, each for at least 9 different species concentrations, totalling 757 melting temperatures, using a melting temperature optimization method. These new parameters do not need to be salt-corrected and are shown to provide overall improved melting temperature predictions. They show a marked quadratic dependence with salt concentrations which are compensated to form linear Gibbs free energies. Two different parameter schemes were tested, with fixed or variable initial parameters. We have found that using variable initial parameters provides better predictive results than using salt correction factors and that the prediction uncertainty is considerably reduced for a validation set of independent sequences. An interpolation scheme is introduced to generate model parameters for arbitrary salt concentrations which performs better against a validation set than predictions using salt corrections.

Molecular rigidity and enthalpy-entropy compensation in DNA melting

Enthalpy-entropy compensation (EEC) is observed in diverse molecular binding processes of importance to living systems and manufacturing applications, but this widely occurring phenomenon is not sufficiently understood from a molecular physics standpoint. To gain insight into this fundamental problem, we focus on the melting of double-stranded DNA (dsDNA) since measurements exhibiting EEC are extensive for nucleic acid complexes and existing coarse-grained models of DNA allow us to explore the influence of changes in molecular parameters on the energetic parameters by using molecular dynamics simulations. Previous experimental and computational studies have indicated a correlation between EEC and changes in molecular rigidity in certain binding-unbinding processes, and, correspondingly, we estimate measures of DNA molecular rigidity under a wide range of conditions, along with resultant changes in the enthalpy and entropy of binding. In particular, we consider variations in dsDNA rigidity that arise from changes of intrinsic molecular rigidity such as varying the associative interaction strength between the DNA bases, the length of the DNA chains, and the bending stiffness of the individual DNA chains. We also consider extrinsic changes of molecular rigidity arising from the addition of polymer additives and geometrical confinement of DNA between parallel plates. All our computations confirm EEC and indicate that this phenomenon is indeed highly correlated with changes in molecular rigidity. However, two distinct patterns relating to how DNA rigidity influences the entropy of association emerge from our analysis. Increasing the intrinsic DNA rigidity increases the entropy of binding, but increases in molecular rigidity from external constraints decreases the entropy of binding. EEC arises in numerous synthetic and biological binding processes and we suggest that changes in molecular rigidity might provide a common origin of this ubiquitous phenomenon in the mutual binding and unbinding of complex molecules.

Impact of thymine glycol damage on DNA duplex energetics: Correlations with lesion-induced biochemical and structural consequences

The magnitude and nature of lesion-induced energetic perturbations empirically correlate with mutagenicity/cytotoxicity profiles and can be predictive of lesion outcomes during polymerase-mediated replication in vitro. In this study, we assess the sequence and counterbase-dependent energetic impact of the Thymine glycol (Tg) lesion on a family of deoxyoligonucleotide duplexes. Tg damage arises from thymine and methyl-cytosine exposure to oxidizing agents or radiation-generated free-radicals. The Tg lesion blocks polymerase-mediated DNA replication in vitro and the unrepaired site elicits cytotoxic lethal consequences in vivo. Our combined calorimetric and spectroscopic characterization correlates Tg -induced energetic perturbations with biological and structural properties. Specifically, we incorporate a 5R-Tg isomer centered within the tridecanucleotide sequence 5'-GCGTACXCATGCG-3' (X = Tg or T) which is hybridized with the corresponding complementary sequence 5'-CGCATGNGTACGC-3' (N = A, G, T, C) to generate families of Tg -damaged (Tg ·N) and lesion-free (T·N) duplexes. We demonstrate that the magnitude and nature of the Tg destabilizing impact is dependent on counterbase identity (i.e., A ∼ G < T < C). The observation that a Tg lesion is less destabilizing when positioned opposite purines suggests that favorable counterbase stacking interactions may partially compensate lesion-induced perturbations. Moreover, the destabilizing energies of Tg ·N duplexes parallel their respective lesion-free T·N mismatch counterparts (i.e., G < T < C). Elucidation of Tg-induced destabilization relative to the corresponding undamaged mismatch energetics allows resolution of lesion-specific and sequence-dependent impacts. The Tg-induced energetic perturbations are consistent with its replication blocking properties and may serve as differential recognition elements for discrimination by the cellular repair machinery.

Hydrogen bonding or stacking interactions in differentiating duplex stability in oligonucleotides containing synthetic nucleoside probes for alkylated DNA

Oligonucleotides that hybridize to modified DNA are useful chemical tools to probe the noncovalent interactions that stabilize DNA duplexes. In an effort to better understand the interactions that influence the specificity of hybridization probes for O(6)-alkylguanine lesions, we examined a series of synthetic nucleoside analogues (BIM, Benzi, and Peri) with respect to their ability to stabilize duplex DNA comprised of native or damaged DNA oligonucleotides. The base-modified nucleoside analogues contained systematically varied hydrogen-bonding and π-stacking properties. The nucleoside probes were incorporated into DNA and paired opposite canonical bases (A, T, C, or G), O(6) -methylguanine (O(6)-MeG), O(6)-benzylguanine (O(6)-BnG), or a stable abasic site analogue (tetrahydrofuran, THF). On the basis of the free energy of duplex formation, the highest degree of stabilization was observed when Peri was paired opposite O(6)-MeG. The thermodynamic data suggest that the smaller probes stabilize DNA duplexes more through hydrogen bonding, whereas the larger probes, with a greater capacity to π stack, contribute to duplex stabilization more on the basis of base stacking. These results demonstrate that increased helix stability could be achieved when BIM, Benzi, or Peri were paired opposite damage-containing DNA rather than unmodified DNA (that is, O(6)-MeG rather than G). This knowledge is expected to be useful in the design and development of nucleoside analogues for uses in DNA-based technologies.

Formation and analysis of topographical domains between lipid membranes tethered by DNA hybrids of different lengths

We recently described a strategy to prepare DNA-tethered lipid membranes either to fixed DNA on a surface or to DNA displayed on a supported bilayer [Boxer et al., J. Struct. Biol., 2009, 168, 190; Boxer et al., Langmuir, 2011, 27, 5492]. With the latter system, the DNA hybrids are laterally mobile; when orthogonal sense-antisense pairs of different lengths are used, the DNA hybrids segregate by height and the tethered membrane deforms to accommodate the height difference. This architecture is particularly useful for modelling interactions between membranes mediated by molecular recognition and resembles cell-to-cell junctions. The length, affinity and population of the DNA hybrids between the membranes are completely controllable. Interesting patterns of height segregation are observed by fluorescence interference contrast microscopy. Diverse behavior is observed in the segregation and pattern forming process and possible mechanisms are discussed. This model system captures some of the essential physics of synapse formation and is a step towards understanding lipid membrane behaviour in cell-to-cell junctions.

Entropy-enthalpy Compensation of Biomolecular Systems in Aqueous Phase: a Dry Perspective

We survey thermodynamic measurements on processes involving biological macromolecules in aqueous solution, which illustrate well the ubiquitous phenomenon of entropy-enthalpy compensation. The processes include protein folding/unfolding and ligand binding/unbinding, with compensation temperatures varying by about 50 K around an average near 293 K. We show that incorporating both near-exact entropy-enthalpy compensation (due to solvent relaxation) and multi-excitation entropy (from vibrational quanta) leads to a compensation temperature in water of about 230 K. We illustrate a general procedure for subtracting solvent and environment-related terms to determine the bare Gibbs free energy changes of chemical processes.