The free energy of DNA supercoiling is enthalpy-determined

Development and evaluation of a simple and effective real time PCR assay for mitochondrial quantification in racing camels

Camel racing is a popular sport in the Middle East region, where the demand is high for racing camels with higher stamina and endurance. Devising a technique to measure oxidative capacity and endurance in camels should be useful. Mitochondria are highly specialized organelles involved in metabolism in all higher organisms for sustaining life and providing energy for physical functions. The ratio of mitochondrial DNA (mtDNA) to nuclear DNA (nDNA) is often used as an estimate for the metabolic status of the tissue. A greater quantity of mitochondria per unit of tissue translates into greater oxidative capacity and endurance. In this report, we describe a simple, sensitive and efficient real-time PCR assay for the quantification of blood mitochondria in racing camels. The primer sequences selected for the SYBR green-based PCR assay included mitochondrial D-loop region, mitochondrial ATP6ase gene and the nuclear β-actin gene. The assay was validated using two groups of camels comprising racing and dairy camels. The racing camels demonstrated a higher mtDNA/nDNA ratio compared with dairy camels based on the ΔΔCt values, with a higher variability among racing camels. The mean ΔΔCt values of adult and young racing camels did not vary considerably. The findings show that the present assay can be used as an evaluative tool for racing camels.

Determining DNA supercoiling enthalpy by isothermal titration calorimetry

DNA supercoiling plays a critical role in certain essential DNA transactions, such as DNA replication, recombination, and transcription. For this reason, exploring energetics of DNA supercoiling is fundamentally important for understanding its biological functions. In this paper, using a unique property of DNA intercalators, such as ethidium bromide and daunorubicin, which bind to supercoiled, nicked, and relaxed DNA templates with different DNA-binding enthalpies, we determined DNA supercoiling enthalpy of plasmid pXXZ6, a 4.5 kb plasmid to be about 11.5 kcal/mol per linking number change. This determination allowed us to partition the DNA supercoiling free energy into enthalpic and entropic contributions where the unfavorable DNA supercoiling free energy exclusively originated from the large positive supercoiling enthalpy and was compensated by a large, favorable entropy term (TΔS).

Tension management in the kinetochore

The kinetochore is the protein machine built at the centromere that integrates mechanical force and chemical energy from dynamic microtubules into directed chromosome motion. The kinetochore also provides a powerful signaling function that is able to alter the properties of the spindle checkpoint and initiate a signal transduction cascade that leads to inhibition of the anaphase promoting complex and cell cycle arrest. Together, the kinetochore accomplishes the feat of chromosome segregation with unparalleled accuracy. Errors in segregation lead to Down's syndrome, the most frequent inherited birth defect, pregnancy loss, and cancer. Over a century after the discovery of the kinetochore, an architectural map comprising greater than 100 proteins is emerging. Understanding the architecture and physical biology of the key components provides new insights into how this fascinating machine moves genomes.

Temperature dependence of DNA persistence length

We have determined the temperature dependence of DNA persistence length, a, using two different methods. The first approach was based on measuring the j-factors of short DNA fragments at various temperatures. Fitting the measured j-factors by the theoretical equation allowed us to obtain the values of a for temperatures between 5°C and 42°C. The second approach was based on measuring the equilibrium distribution of the linking number between the strands of circular DNA at different temperatures. The major contribution into the distribution variance comes from the fluctuations of DNA writhe in the nicked circular molecules which are specified by the value of a. The computation-based analysis of the measured variances was used to obtain the values of a for temperatures up to 60°C. We found a good agreement between the results obtained by these two methods. Our data show that DNA persistence length strongly depends on temperature and accounting for this dependence is important in quantitative comparison between experimental results obtained at different temperatures.

The effect of changes in the bulk dielectric constant on the DNA torsional properties

The effect of changes in the bulk dielectric constant on the DNA torsional properties was evaluated from plasmid circularization reactions. In these reactions, pUC18 previously linearized by EcoRI digestion was recircularized with T4 DNA ligase. The bulk dielectric constant of the reaction medium was decreased by the addition of different concentrations of neutral solutes: ethylene glycol, glycerol, sorbitol, and sucrose, or increased by the addition of glycine. The topoisomers generated by the ligase reaction were resolved by agarose-gel electrophoresis. The DNA twist energy parameter (kappa), which is an apparent torsional constant, was determined by linearization of the Gaussian topoisomers' distribution. It was observed that the twist energy parameter for the given solutes is almost linearly dependent on the bulk dielectric constant. In the reaction buffer, the twist energy parameter was determined to be 1100 +/- 100. By decreasing the dielectric constant to 74 with the addition of sorbitol, the value of the parameter reaches kappa = 900 +/- 100, whereas the addition of ethylene glycol leads to kappa = 400 +/- 50. Upon addition of glycine, which resulted in a dielectric constant equal to 91, the value of the twist energy parameter increased to kappa = 1750 +/- 100.

DNA Mechanics

We review the history of DNA mechanics and its analysis. We evaluate several methods to analyze the structures of superhelical DNA molecules, each predicated on the assumption that DNA can be modeled with reasonable accuracy as an extended, linearly elastic polymer. Three main approaches are considered: mechanical equilibrium methods, which seek to compute minimum energy conformations of topologically constrained molecules; statistical mechanical methods, which seek to compute the Boltzmann distribution of equilibrium conformations that arise in a finite temperature environment; and dynamic methods, which seek to compute deterministic trajectories of the helix axis by solving equations of motion. When these methods include forces of self-contact, which prevent strand passage and preserve the topological constraint, each predicts plectonemically interwound structures. On the other hand, the extent to which these mechanical methods reliably predict energetic and thermodynamic properties of superhelical molecules is limited, in part because of their inability to account explicitly for interactions involving solvent. Monte Carlo methods predict the entropy associated with supercoiling to be negative, in conflict with a body of experimental evidence that finds it is large and positive, as would be the case if superhelical deformations significantly disrupt the ordering of ambient solvent molecules. This suggests that the large-scale conformational properties predicted by elastomechanical models are not the only ones determining the energetics and thermodynamics of supercoiling. Moreover, because all such models that preserve the topological constraint correctly predict plectonemic interwinding, despite these and other limitations, this constraint evidently dominates energetic and thermodynamic factors in determining supercoil geometry. Therefore, agreement between predicted structures and structures obtained experimentally, for example, by electron microscopy, does not in itself provide evidence for the correctness or completeness of any given model of DNA mechanics.

The question of long-range allosteric transitions in DNA

The question of long-range allosteric transitions of DNA secondary structure and their possible involvement in transcriptional activation is discussed in the light of new results. A variety of recent evidence strongly supports a fluctuating long-range description of DNA secondary structure. Balanced equilibria between two or more different secondary structures, and the occurrence of very large domain sizes, have been documented in several instances. Long-range allosteric effects stemming from changes in sequence or secondary structure over a small region of the DNA have been observed to extend over distances up to hundreds of base pairs in some cases. The discovery that coherent bending strain beyond a threshold level in small (N < or = 250 base pairs (bp)] circular DNAs significantly alters the DNA secondary structure has important implications, especially for transcriptional activators that either bend the DNA directly or are involved in the formation of DNA loops of sufficiently small size (N < or = 250 bp). Whether the RNA polymerase is activated primarily via protein: protein contacts, as is widely believed, or instead via a bend-induced allosteric transition of the DNA in such a small loop, is now an open question. Binding of the transcriptional activator Sp1 to linear DNA induces a remarkably long-range change in its secondary structure, and catabolite activator protein binding to a supercoiled DNA behaves similarly, though possibly for different reasons. Compelling evidence for a bend-induced long-range structural transmission effect of the transcriptional activator integration host factor on RNA polymerase activity was recently reported. These results may augur a new paradigm in which allosteric transitions of duplex DNA, as well as of the proteins, are involved in the regulation of transcription.

Electrostatic-undulatory theory of plectonemically supercoiled DNA

We present an analytical calculation of the electrostatic interaction in a plectonemic supercoil within the Poisson-Boltzmann approximation. Undulations of the supercoil strands arising from thermal motion couple nonlinearly with the electrostatic interaction, giving rise to a strong enhancement of the bare interaction. In the limit of fairly tight winding, the free energy of a plectonemic supercoil may be split into an elastic contribution containing the bending and torsional energies and an electrostatic-undulatory free energy. The total free energy of the supercoil is minimized according to an iterative scheme, which utilizes the special symmetry inherent in the usual elastic free energy of the plectoneme. The superhelical radius, opening angle, and undulation amplitudes in the radius and pitch are obtained as a function of the specific linking difference and the concentration of monovalent salt. Our results compare favorably with the experimental values for these parameters of Boles et al. (1990. J. Mol. Biol. 213:931-951). In particular, we confirm the experimental observation that the writhe is a virtually constant fraction of the excess linking number over a wide range of superhelical densities. Another important prediction is the ionic strength dependence of the plectonemic parameters, which is in reasonable agreement with the results from computer simulations.

Effect of temperature on DNA secondary structure in the absence and presence of 0.5 M tetramethylammonium chloride

Changes in the average secondary structures of three different linear DNAs over the premelting region from 5 to 60 degrees C were investigated by measuring their CD spectra and also their torsion elastic constants (<alpha>) by time-resolved fluorescence polarization anisotropy. For one of these DNAs, the Haell fragment of pBR322, the apparent diffusion coefficients [Dapp(k)] at small and large scattering vectors (k) were also measured by dynamic light scattering. With increasing temperature, all three DNAs exhibited typical premelting changes in their CD spectra, and these were accompanied by 1.4- to 1.7-fold decreases in <alpha>. Also for the 1876 base pair fragment, Dapp(k) at large scattering vectors, which is sensitive to the dynamic bending rigidity, decreased by 17%, even though there was no change at small scattering vectors, where Dapp(k) = D0 is the translational diffusion coefficient of the center-of-mass. These observations demonstrate conclusively that the premelting CD changes of these DNAs are associated with a significant change in average secondary structure and mechanical properties, though not in persistence length. In the presence of 0.5 M tetramethylammonium chloride (TMA-Cl) the premelting change in CD is largely suppressed, and the corresponding changes in <alpha> and Dapp(k) at large scattering vectors are substantially diminished. These observations suggest that TMA-Cl, which binds preferentially to A.T-rich regions and stabilizes those regions (relative to G.C-rich regions) against melting, effectively stabilizes the prevailing low-temperature secondary structure sufficiently that the DNA is effectively trapped in that state over the temperature range observed.

On the origin of the temperature dependence of the supercoiling free energy

Monte Carlo simulations using temperature-invariant torsional and bending rigidities fail to predict the rather steep decline of the experimental supercoiling free energy with increasing temperature, and consequently fail to predict the correct sign and magnitude of the supercoiling entropy. To illustrate this problem, values of the twist energy parameter (E(T)), which governs the supercoiling free energy, were simulated using temperature-invariant torsion and bending potentials and compared to experimental data on pBR322 over a range of temperatures. The slope, -dE(T)/dT, of the simulated values is also compared to the slope derived from previous calorimetric data. The possibility that the discrepancies arise from some hitherto undetected temperature dependence of the torsional rigidity was investigated. The torsion elastic constant of an 1876-bp restriction fragment of pBR322 was measured by time-resolved fluorescence polarization anisotropy of intercalated ethidium over the range 278-323 K, and found to decline substantially over that interval. Simulations of a 4349-bp model DNA were performed using these measured temperature-dependent torsional rigidities. The slope, -dE(T)/dT, of the simulated data agrees satisfactorily with the slope derived from previous calorimetric measurements, but still lies substantially below that of Duguet's data. Models that involve an equilibrium between different secondary structure states with different intrinsic twists and torsion constants provide the most likely explanation for the variation of the torsion constant with T and other pertinent observations.

Monte Carlo simulations of supercoiling free energies for unknotted and trefoil knotted DNAs

A new Monte Carlo (MC) algorithm is proposed for simulating inextensible circular chains with finite twisting and bending rigidity. This new algorithm samples the relevant Riemann volume elements in a uniform manner, when the constraining potential vanishes. Simulations are performed for filaments comprising 170 subunits, each containing approximately 28 bp, which corresponds to a DNA length of 4770 bp. The bending rigidity is chosen to yield a persistence length, P = 500 A, and the intersubunit potential is taken to be a hard-cylinder potential with diameter d = 50 A. This value of d yields the same second virial coefficient as the electrostatic potential obtained by numerical solution of the Poisson-Boltzmann equation for 150 mM salt. Simulations are performed for unknotted circles and also for trefoil knotted circles using two different values of the torsional rigidity, C = (2.0 and 3.0) x 10(-19) dyne cm2. In the case of unknotted circles, the simulated supercoiling free energy varies practically quadratically with linking difference delta l. The simulated twist energy parameter ET, its slope dET/dT, and the mean reduced writhe <w>/delta l for C = 3 x 10(-19) dyne cm2 all agree well with recent simulations for unknotted circles using the polygon-folding algorithm with identical P, d, and C. The simulated ET vs. delta l data for C = 2.0 x 10(-19) dyne cm2 agree rather well with recent experimental data for p30 delta DNA (4752 bp), for which the torsional rigidity, C = 2.07 x 10(-19) dyne cm2, was independently measured. The experimental data for p30 delta are enormously more likely to have arisen from C = 2.0 x 10(-19) than from C = 3.0 x 10(-19) dyne cm2. Serious problems with the reported experimental assessments of ET for pBR322 and their comparison with simulated data are noted. In the case of a trefoil knotted DNA, the simulated value, (ET)tre, exceeds that of the unknotted DNA, (ET)unk, by approximately equal to 1.40-fold at magnitude of delta l = 1.0, but declines to a plateau about 1.09-fold larger than (ET)unk when magnitude of delta l > or = 15. Although the predicted ratio, (ET)tre/(ET)unk approximately equal to 1.40, agrees fairly well with recent experimental measurements on a 5600-bp DNA, the individual measured ET values, like some of those reported for pBR322, are so large that they cannot be simulated using P = 500 A, d = 50 A, and any previous experimental estimate of C.

Thermodynamics of condensation of nuclear chromatin. A differential scanning calorimetry study of the salt-dependent structural transitions

We present a detailed thermodynamic investigation of the conformational transitions of chromatin in calf thymus nuclei. Differential scanning calorimetry was used as the leading method, in combination with infrared spectroscopy, electron microscopy, and techniques for the molecular characterization of chromatin components. The conformational transitions were induced by changes in the counterion concentration. In this way, it was possible to discriminate between the interactions responsible for the folding of the higher order structure and for the coiling of nucleosomal DNA. Our experiments confirm that the denaturation of nuclear chromatin at physiological ionic strength occurs at the level of discrete structural domains, the linker and the core particle, and we were able to rule out that the actual denaturation pattern might be determined by dissociation of the nucleohistone complex and successive migration of free histones toward native regions, as recently suggested. The sequence of the denaturation events is (1) the conformational change of the histone complement at 66 degrees C, (2) the unstacking of the linker DNA at 74 degrees C, and (3) the unstacking of the core particle DNA, that can be observed either at 90 or at 107 degrees C, depending on the degree of condensation of chromatin. Nuclear chromatin unfolds in low-salt buffers, and can be refolded by increasing the ionic strength, in accordance with the well-known behavior of short fragments. The process is athermal, therefore showing that the stability of the higher order structure depends on electrostatic interactions. The transition between the folded conformation and the unfolded one proceeds through an intermediate condensation state, revealed by an endotherm at 101 degrees C. The analysis of the thermodynamic parameters of denaturation of the polynucleosomal chain demonstrates that the wrapping of the DNA around the histone octamer involves a large energy change. The most striking observation concerns the linker segment, which melts a few degrees below the peak temperature of naked DNA. This finding is in line with previous thermal denaturation investigations on isolated chromatin at low ionic strength, and suggests that a progressive destabilization of the linker occurs in the course of the salt-induced coiling of DNA in the nucleosome.

Energy-structure correlations of plasmid DNA in different topological forms

Differential scanning microcalorimetry (DSC), UV absorption and circular dichroism (CD) have been used to study structure and stability of linear (lin), open circular (oc), supercoiled (cd) and relaxed circular duplex (rd) DNA and calf thymus (CT) DNA. Investigations were made in low salt buffer and in the presence of 7.2 M NaClO4. The chaotropic action of perchlorate promotes a reduction of the overall stability of DNA, which permits a direct determination of the transition enthalpies of all four DNA configurations. The stabilities against thermal denaturation have been found to increase in the series lin approximately oc less than cd less than rd. These relative stabilities can be rationalized on the basis of the linkage between supercoiling and secondary structural changes in topologically constrained duplex DNA. On the basis of these studies, a model of the melting process could be suggested that is consistent with the energetic and spectroscopic data.

Melting transition of covalently closed DNA with supercoil-induced cruciforms

The melting curve for covalently closed supercoiled DNA has been studied by assuming the existence of cruciforms as significant structural perturbations in the pre-melting region. The statistical mechanical treatment used incorporates these cruciform structures through an appropriate sequence generating function. The variation of the effective hydrogen bond energy with temperature is taken into account by an empirical procedure. The results obtained are in close agreement with the corresponding experimental data in TEA solution where the effect of heterogeneity of the base pairs is minimized.

The flexibility of alternating dA-dT sequences

The flexibility of alternating poly (dA-dT) has been investigated by the technique of transient electric dichroism. Rotational relaxation times, which are very sensitive to changes in the end-to-end length of flexible polymers, are determined from the field free dichroism decay curves of four, well defined fragments of poly (dA-dT) ranging in size from 136 to 270 base pairs. Persistence lengths, calculated from the results of Hagerman and Zimm (Biopolymers (1981) 29, 1481-1502), are in the range 200-250 A. This makes alternating dA-dT sequences about twice as flexible as naturally occurring, "random" sequence DNA. Considering a bend around a nucleosome, for example, this difference in persistence length translates to an energy difference between poly (dA-dT) and random sequence DNA of 0.17 kT/base pair or 1 kcal per 10 base pair stretch. This energy difference is sufficiently large to suggest that dA-dT sequences could serve as markers in DNA packaging, for example, at sites where DNA must tightly bend to accommodate structures.