Statistical analysis of plasma thermograms measured by differential scanning calorimetry.

Melting curves of human plasma measured by differential scanning calorimetry (DSC), known as thermograms, have the potential to markedly impact diagnosis of human diseases. A general statistical methodology is developed to analyze and classify DSC thermograms to analyze and classify thermograms. Analysis of an acquired thermogram involves comparison with a database of empirical reference thermograms from clinically characterized diseases. Two parameters, a distance metric, P, and correlation coefficient, r, are combined to produce a 'similarity metric,' ρ, which can be used to classify unknown thermograms into pre-characterized categories. Simulated thermograms known to lie within or fall outside of the 90% quantile range around a median reference are also analyzed. Results verify the utility of the methods and establish the apparent dynamic range of the metric ρ. Methods are then applied to data obtained from a collection of plasma samples from patients clinically diagnosed with SLE (lupus). High correspondence is found between curve shapes and values of the metric ρ. In a final application, an elementary classification rule is implemented to successfully analyze and classify unlabeled thermograms. These methods constitute a set of powerful yet easy to implement tools for quantitative classification, analysis and interpretation of DSC plasma melting curves.

Calf-thymus topoisomerase I equilibrates metastable secondary structure subsequent to relaxation of superhelical stress.

After relaxation of superhelical stress by various methods not involving topoisomerases, a long-lived metastable secondary structure with an anomalously low torsion elastic constant commonly prevails. The aim here is to ascertain whether such metastable secondary structure also results from the action of calf-thymus topoisomerase I (CT Topo I) on a native supercoiled DNA and, if so, whether the enzyme catalyzes its subsequent equilibration. The action of CT Topo I on supercoiled p30delta DNA was examined over a range of times from 10 min to 6 h. We verify that the enzyme operates in an almost completely processive manner, and at each time point determine the twist energy parameter, E(T), that governs the supercoiling free energy. E(T) is initially low, 533 +/- 60, and remains essentially constant up to at least 360 min, when no further CT Topo I is added. The activity of the rather dilute enzyme dies within approximately 60 min. During the 60 min after a second addition of fresh enzyme at either 60 or 120 min, E(T) rises up to a plateau at approximately 1100, which lies within the consensus equilibrium range, 1000 +/- 100. Over that same time period, the average peak spacing between the gel bands (corresponding to individual topoisomers) decreases somewhat with increasing time of exposure to active CT Topo I. After a third addition of fresh CT Topo I at 240 min, there is no further change in either E(T) or the average gel spacing. These and other observations indicate that active CT Topo I catalyzes the equilibration of a metastable secondary structure with abnormally low torsion and bending elastic constants that prevails after the initial release of superhelical stress. An observed temporal lag of this structural equilibration behind the relaxation of native superhelical DNAs suggests that it may require cleavage and religation events at multiple sites on the DNA. A novel analysis of the unwinding kinetics using literature data accounts for the almost complete processivity of the enzyme. The action of CT Topo I was also examined in the presence of 20 and 40 w/v% ethylene glycol (EG), which shift a secondary structure equilibrium toward an alternative state with altered torsion and bending elastic constants. The present results suggest that the usual metastable state coexists with the EG-induced state, and is equilibrated more rapidly than in the absence of EG.

Origins of the "nucleation" free energy in the hybridization thermodynamics of short duplex DNA.

Thermodynamic parameters deltaH(cal), deltaS(cal), and deltaG(cal) of the melting transitions for 19 short DNA/DNA duplexes ranging in length from 6 to 35 base pairs were systematically evaluated by differential scanning calorimetry melting experiments carried out at four salt concentrations from 85 mM to 1.0 M [Na+]. As expected, thermodynamic stabilities of the DNA duplexes increased with length and increasing [Na+]. From plots of deltaG25 versus duplex length, extrapolation to N = 0 provided estimates on values of deltaG(cal)25 (N = 0) as a function of [Na+], corresponding to the free-energy of the "hypothetical duplex" having zero base pairs, but occupying precisely the same molar volume as the fully base paired duplex. The values obtained for deltaG(cal)25 (N = 0) were 3.68, 5.59, 7.86, and 8.68 kcal/mol in 1.00, 0.60, 0.30, and 0.085 M Na+, respectively. These values are in reasonable agreement with published values of the nucleation or initiation free-energy, attributed to formation of the first base pair in a short duplex compared to formation of the remaining base pairs. A statistical thermodynamic formulation of the association of two strands accounting for displaced solvent was utilized to relate [Na+]-dependent deltaG(cal)25 (N = 0) values to configuration integrals for both single and duplex strands. Relative differences between two single strands in their standard states and the duplex (in its standard state), and solvent displaced during the annealing process was taken into account. This analysis provides a new vantage point to view what has historically been referred to as the helix initiation or nucleation parameter and provides an alternate interpretation and mechanism for the nucleation complex in duplex formation.

Conformational equilibria of bulged sites in duplex DNA studied by EPR spectroscopy.

Conformational flexibility in nucleic acids provides a basis for complex structures, binding, and signaling. One-base bulges directly neighboring single-base mismatches in nucleic acids can be present in a minimum of two distinct conformations, complicating the examination of the thermodynamics by calorimetry or UV-monitored melting techniques. To provide additional information about such structures, we demonstrate how electron paramagnetic resonance (EPR) active spin-labeled base analogues, base-specifically incorporated into the DNA, are monitors of the superposition of different bulge-mismatch conformations. EPR spectra provide information about the dynamic environments of the probe. This information is cast in terms of "dynamic signatures" that have an underlying basis in structural variations. By examining the changes in the equilibrium of the different states across a range of temperatures, the enthalpy and entropy of the interconversion among possible conformations can be determined. The DNA constructs with a single bulge neighboring a single-base mismatch ("bulge-mismatches") may be approximately modeled as an equilibrium between two possible conformations. This structural information provides insight into the local composition of the bulge-mismatch sequences. Experiments on the bulge-mismatches show that basepairing across the helix can be understood in terms of purine and pyrimidine interactions, rather than specific bases. Measurements of the enthalpy and entropy of formation for the bulge-mismatches by differential scanning calorimetry and UV-monitored melting confirm that the formation of bulge-mismatches is in fact more complicated than a simple two-state process, consistent with the base-specific spectral data that bulge-mismatches exist in multiple conformations in the premelting temperature region. We find that the calculations with the nearest-neighbor (NN) model for the two likely conformations do not correlate well with the populations of structures and thermodynamic parameters inferred from the base-specific EPR dynamics probe. We report that the base-specific spin probes are able to identify a bistable, temperature dependent, switching between conformations for a particular complex bulged construct.

Influence of buffer species on the thermodynamics of short DNA duplex melting: sodium phosphate versus sodium cacodylate.

Thermodynamic parameters of the melting transitions of 53 short duplex DNAs were experimentally evaluated by differential scanning calorimetry melting curve analysis. Solvents for the DNA solutions contained approximately 1 M Na+ and either 10 mM cacodylate or phosphate buffer. Thermodynamic parameters obtained in the two solvent environments were compared and quantitatively assessed. Thermodynamic stabilities (deltaG(o) (25 degrees C)) of the duplexes studied ranged from quite stable perfect match duplexes (approximately -30 kcal/mol) to relatively unstable mismatch duplexes (approximately -9 kcal/mol) and ranged in length from 18 to 22 basepairs. A significant difference in stability (average free energy difference of approximately 3 kcal/mol) was found for all duplexes melted in phosphate (greater stability) versus cacodylate buffers. Measured effects of buffer species appear to be relatively unaffected by duplex length or sequence content. The popular sets of published nearest-neighbor (n-n) stability parameters for Watson-Crick (w/c) and single-base mismatches were evaluated from melting studies performed in cacodylate buffer (SantaLucia and Hicks, Annu. Rev. Biophys. Biomol. Struct. 2004, 33, 415). Thus, when using these parameters to make predictions of sequence dependent stability of DNA oligomers in buffers other than cacodylate (e.g., phosphate) one should be mindful that in addition to sodium ion concentration, the type of buffer species also provides a minor but significant contribution to duplex stability. Such considerations could potentially influence results of sequence dependent analysis using published n-n parameters and impact results of thermodynamic calculations. Such calculations and analyses are typically employed in the design and interpretation of DNA multiplex hybridization experiments.

A structural transition in duplex DNA induced by ethylene glycol.

The twist energy parameter ( E T) that governs the supercoiling free energy, and the linking difference (Delta l) are measured for p30delta DNA in solutions containing 0-40 w/v % ethylene glycol (EG). A plot of E T vs -ln a w, where a w is the water activity, displays the full (reverse) sigmoidal profile of a discrete structural transition. A general theory for the effect of added osmolyte on a cooperative structural transition between two duplex states, 1 right arrow over left arrow 2, is formulated in terms of parameters applicable to individual base-pair subunits. The resulting fraction of base pairs in the 2-state ( f 2 (0)) is incorporated into expressions for the effective torsion and bending elastic constants, the effective twist energy parameter ( E T (eff)), and the change in intrinsic twist (delta l 0). Fitting the expression for E T (eff) to the measured E T values yields reasonably unambiguous estimates of E T 1 and E T 2 , the midpoint value (ln a w) 1/2, and the midpoint slope ( partial differential E T/ partial differential ln a w) 1/2, but does not yield unambiguous estimates of the equilibrium constant ( K 0), the difference in DNA-water preferential interaction coefficient (DeltaGamma), or the inverse cooperativity parameter ( J). Fitting a noncooperative model (assumed J = 1.0) to the data yields K 0 = 0.067 and DeltaGamma = -30.0 per base pair (bp). Essentially equivalent fits are provided by models with a wide range of correlated J, DeltaGamma, and K 0 values. Other results favor DeltaGamma in the range -1.0 to 0, which then requires K 0 > or = 0.914, and a cooperativity parameter, 1/ J > or = 30.0 bp. The measured delta l 0 and circular dichroism (CD) at 272 nm are found to be compatible with curves predicted using the same f 2 (0) values that best-fit the E T data. At least 7-10% of the base pairs are inferred to exist in the 2-state in 0.1 M NaCl in the complete absence of added osmolyte. Compared with the 1-state, the 2-state has a approximately 2.0- to 2.1-fold greater torsion elastic constant, a approximately 0.70-fold smaller bending elastic constant, a approximately 0.91-fold smaller E T value, a approximately 0.2% lower intrinsic twist, a somewhat lower CD near both 272 and 245 nm, and less water and/or more EG in its neighborhood. However, the relative change in preferential interaction coefficient associated with the transition is likely rather slight.

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 degrees C to 80 degrees 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.

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 DeltaH degrees , DeltaS degrees and DeltaG degrees 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.

Multiplex SNP discrimination.

Multiplex hybridization reactions of perfectly matched duplexes and duplexes containing a single basepair mismatch (SNPs) were investigated on DNA microarrays. Effects of duplex length, G-C percentage, and relative position of the SNP on duplex hybridization and SNP resolution were determined. Our theoretical model of multiplex hybridization accurately predicts observed results and implicates target concentration as a critical variable in multiplex SNP detection.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30delta DNA.

Upon increasing the concentration of ethylene glycol (EG) at 37 degrees C, the twist energy parameter, E(T), which governs the supercoiling free energy, was recently found to undergo a decreasing (or reverse) sigmoidal transition with a midpoint near 20 w/v % EG. In this study, the effects of adding 20 w/v % EG on the torsion elastic constant (alpha) of linear p30delta DNA and on the hydrodynamic radius (R(H)) of a synthetic 24 bp duplex DNA were examined at both 40 and 20 degrees C. The time-resolved fluorescence intensity and fluorescence polarization anisotropy (FPA) of intercalated ethidium were measured in order to assess the effects of 20 w/v % EG on: (1) alpha; (2) R(H); (3) the lifetimes of intercalated and non-intercalated dye; (4) the amplitude of dye wobble in its binding site; and (5) the binding constant for intercalation. The effects of 20 w/v % EG on the circular dichroism (CD) spectrum of the DNA and on the emission spectrum of the free dye were also measured. At 40 degrees C, addition of 20 w/v % EG caused a substantial (1.27- to 1.35-fold) increase in alpha, a significant change in the CD spectrum, and a very small, marginally significant increase in R(H), but little or no change in the amplitude of dye wobble in its binding site or the lifetime of intercalated dye. Together with previously reported measurements of E(T), these results imply that the bending elastic constant of DNA is significantly decreased by 20 w/v % EG at 40 degrees C. At 20 degrees C, addition of 20 w/v % EG caused a marginally significant decrease in alpha and very little change in any other measured properties. Also at 20 degrees C, addition of 30 w/v % betaine caused a marginally significant increase in alpha and significant but modest change in the CD spectrum, but very little change in any other properties.