Albert S. Benight

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.

Studies of DNA dumbbells VII: Evaluation of the next-nearest-neighbor sequence-dependent interactions in duplex DNA

Melting experiments were conducted on 22 DNA dumbbells as a function of solvent ionic strength from 25–115 mM Na+. The dumbbell molecules have short duplex regions comprised of 16–20 base pairs linked on both ends by T4 single‐strand loops. Only the 4–8 central base pairs of the dumbbell stems differ for different molecules, and the six base pairs on both sides of the central sequence and adjoining loops on both ends are the same in every molecule. Results of melting analysis on the 22 new DNA dumbbells are combined with our previous results on 17 other DNA dumbbells, with stem lengths containing from 14–18 base pairs, reported in the first article of this series (Doktycz, Goldstein, Paner, Gallo, and Benight, Biopoly 32, 1992, 849–864). The combination of results comprises a database of optical melting parameters for 39 DNA dumbbells in ionic strengths from 25–115 mM Na+. This database is employed to evaluate the thermodynamics of singlet, doublet, and triplet sequence‐dependent interactions in duplex DNA. Analysis of the 25 mM Na+ data reveals the existence of significant sequence‐dependent triplet or next‐nearest‐neighbor interactions. The enthalpy of these interactions is evaluated for all possible triplets. Some of the triplet enthalpy values are less than the uncertainty in their evaluation, indicating no measurable interaction for that particular sequence. This finding suggests that the thermodynamic stability of duplex DNA depends on solvent ionic strength in a sequence‐dependent manner. As a part of the analysis, the nearest‐neighbor (base pair doublet) interactions in 55, 85, and 115 mM Na+ are also reevaluated from the larger database.

Electrical detection of the temperature induced melting transition of a DNA hairpin covalently attached to gold interdigitated microelectrodes

The temperature induced melting transition of a self-complementary DNA strand covalently attached at the 5′ end to the surface of a gold interdigitated microelectrode (GIME) was monitored in a novel, label-free, manner. The structural state of the hairpin was assessed by measuring four different electronic properties of the GIME (capacitance, impedance, dissipation factor and phase angle) as a function of temperature from 25°C to 80°C. Consistent changes in all four electronic properties of the GIME were observed over this temperature range, and attributed to the transition of the attached single-stranded DNA (ssDNA) from an intramolecular, folded hairpin structure to a melted ssDNA. The melting curve of the self-complementary single strand was also measured in solution using differential scanning calorimetry (DSC) and UV absorbance spectroscopy. Temperature dependent electronic measurements on the surface and absorbance versus temperature values measured in solution experiments were analyzed assuming a two-state process. The model analysis provided estimates of the thermodynamic transition parameters of the hairpin on the surface. Two-state analyses of optical melting data and DSC measurements provided evaluations of the thermodynamic transition parameters of the hairpin in solution. Comparison of surface and solution measurements provided quantitative evaluation of the effect of the surface on the thermodynamics of the melting transition of the DNA hairpin.

The helix-coil transition in heterogeneous double stranded DNA: Microcanonical method

A microscopic Potts-like one-dimensional model with many-particle interactions is developed to construct a statistical mechanical description of the melting of heterogeneous sequence duplex DNA. For this model, referred as the generalized model of polypeptide chains (GMPC), a closed-form expression for the free energy is derived. The characteristic equation of the model enables estimates on the melting temperature and transition interval, consistent with results obtained from more classical approaches. From the characteristic equation of the model, the temperature-dependent statistical weight parameter for helical states is evaluated. This parameter is shown to change throughout the transition from a harmonic form in early regions of the transition to an arithmetic form in later stages. The GMPC is extended to consider the influence of sequence heterogeneity in the melting of duplex DNA.

Competition for hydrogen-bond formation in the helix-coil transition and protein folding

The problem of the helix-coil transition of biopolymers in explicit solvents, such as water, with the ability for hydrogen bonding with a solvent is addressed analytically using a suitably modified version of the Generalized Model of Polypeptide Chains. Besides the regular helix-coil transition, an additional coil-helix or reentrant transition is also found at lower temperatures. The reentrant transition arises due to competition between polymer-polymer and polymer-water hydrogen bonds. The balance between the two types of hydrogen bonding can be shifted to either direction through changes not only in temperature, but also by pressure, mechanical force, osmotic stress, or other external influences. Both polypeptides and polynucleotides are considered within a unified formalism. Our approach provides an explanation of the experimental difficulty of observing the reentrant transition with pressure and underscores the advantage of pulling experiments for studies of DNA. Results are discussed and compared with those reported in a number of recent publications with which a significant level of agreement is obtained.

Stacking heterogeneity: A model for the sequence dependent melting cooperativity of duplex DNA

A microscopic Potts-like one-dimensional model with many particle interactions [referred as the generalized model of polypeptide chains (GMPCs)] is developed to investigate cooperativity of DNA sequence dependent melting. For modeling sequence, regular homogeneous sequences were arranged in heterogeneous blocks of various lengths. Within the framework of the GMPC the authors show that the inclusion of stacking interaction heterogeneity relative to homogeneous hydrogen bond interactions leads to an unexpected and quite remarkable increase in melting cooperativity for small blocks. In some cases this tendency persists for long blocks having sharp sequence heterogeneity.

Intersegment interactions and helix-coil transition within the generalized model of polypeptide chains approach

The generalized model of polypeptide chains is extended to describe the helix-coil transition in a system comprised of two chains interacting side-by-side. The Hamiltonian of the model takes into account four possible types of interactions between repeated units of the two chains, i.e., helix-helix, helix-coil, coil-helix, and coil-coil. Analysis reveals when the energy I(hh)+I(cc) of (h-h, c-c) interactions overwhelms the energy I(hc)+I(ch) of mixed (h-c, c-h) interactions, the correlation length rises substantially, resulting in narrowing of the transition interval. In the opposite case, when I(hh)+I(cc)<I(hc)+I(ch), nontrivial behavior of the system is predicted where an intermediate plateau appears on the denaturation curve. For the latter case, calculations of the number of junctions and the average length of helical segments indicate rearrangement of helical segments at the transition point. Conceptual links are established with experimentally oriented theories of Ghosh and Dill [J. Am. Chem. Soc. 131, 2306 (2009)] and Skolnick and Holtzer [Biochemistry 25, 6192 (1986)], providing a potential explanation for both favorable helix formation and disfavored intersegment interactions from the same theoretical perspective.

Two scale generalized model of polypeptide chains

The generalized model of polypeptide chains (GMPC) is expanded to simultaneously consider two types of interactions occurring over different scales. This new two scale GMPC is applied in several specific cases to examine: The combined influence of stacking or antistacking and hydrogen bonding, or spatial restrictions on the length of helical segments, on the cooperativity and temperature interval of the helix-coil transition of duplex DNA. For the cases of stacking or antistacking in combination with hydrogen bonding the model reduces to the basic uniscale model with a redefined scaling parameter Delta. Antistacking increases the cooperativity, while stacking decreases it. In each case, explanations are given in terms of different lengths of helical segments. Restrictions on the length of helical regions result in the appearance of antiferromagnetic-type correlations where there is no apparent link between cooperativity and transition interval.

The paradox of multiplex DNA melting on a surface

Under equilibrium conditions, there are two regimes of target capture on a surface--target limited and probe limited. In the probe limited regime, the melting curve from multiplex target dissociation from the surface exhibits a single transition due to a reverse displacement mechanism of the low affinity species. The melting curve cannot be used in analytical methods to resolve heteroduplexes; only with the simplex system can proper thermodynamics be obtained.