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

A quantitative model of a cooperative two-state equilibrium in DNA: experimental tests, insights, and predictions

Quantitative parameters for a two-state cooperative transition in duplex DNAs were finally obtained during the last 5 years. After a brief discussion of observations pertaining to the existence of the two-state equilibrium per se, the lengths, torsion, and bending elastic constants of the two states involved and the cooperativity parameter of the model are simply stated. Experimental tests of model predictions for the responses of DNA to small applied stretching, twisting, and bending stresses, and changes in temperature, ionic conditions, and sequence are described. The mechanism and significance of the large cooperativity, which enables significant DNA responses to such small perturbations, are also noted. The capacity of the model to resolve a number of long-standing and sometimes interconnected puzzles in the extant literature, including the origin of the broad pre-melting transition studied by numerous workers in the 1960s and 1970s, is demonstrated. Under certain conditions, the model predicts significant long-range attractive or repulsive interactions between hypothetical proteins with strong preferences for one or the other state that are bound to well-separated sites on the same DNA. A scenario is proposed for the activation of the ilvPG promoter on a supercoiled DNA by integration host factor.

Polyelectrolyte effects in G-quadruplexes

The role of counterion condensation as a dominant force governing the stability of DNA duplexes and triplexes is well established. In contrast, the effect of counterion condensation on the stability of G-quadrupex conformations is poorly understood. Unlike other ordered nucleic acid structures, G-quadruplexes exhibit a specific binding of counterions (typically, Na(+) or K(+)) which are buried inside the central cavity and coordinated to the O6 carbonyls of the guanines forming the G-quartets. While it has been known that the G-quadruplex-to-coil transition temperature, TM, increases with an increase in the concentration of the stabilizing ion, the contributions of the specific (coordination in the central cavity) and nonspecific (condensation) ion binding have not been resolved. In this work, we separate the two contributions by studying the change in TM of preformed G-quadruplexes following the addition of nonstabilizing ions Li(+), Cs(+), and TMA(+) (tetramethylammonium). In our studies, we used two G-quadruplexes formed by the human telomeric sequences which are distinct with respect to the folding topology and the identity and the number of sequestered stabilizing ions. Our data suggest that the predominant ionic contribution to G-quadruplex stability comes from the specifically bound Na(+) or K(+) ions and not from counterion condensation. We offer molecular rationalizations to the observed insensitivity of G-quadruplex stability to counterion condensation and emphasize the need to expand such studies to assess the generality of our findings.

Electrokinetic characterization of superparamagnetic nanoparticle-aptamer conjugates: design of new highly specific probes for miniaturized molecular diagnostics

With the view of designing new nanoparticle (NP)-aptamer conjugates and proving their suitability as biorecognition tools for miniaturized molecular diagnostics, new maghemite-silica core-shell NP-aptamer conjugates were characterized for the first time in terms of grafting rate and colloidal stability under electrophoretic conditions using capillary electrophoresis. After the grafting rate (on the order of six to 50) of the lysozyme-binding aptamer had been estimated, the electrophoretic stability and peak dispersion of the resulting oligonucleotide-NP conjugates were estimated so as to determine the optimal separation conditions in terms of buffer pH, ionic strength and nature, as well as temperature and electric field strength. The effective surface charge density of the NPs was close to zero for pH lower than 5, which led to some aggregation. The NPs were stable in the pH range from 5 to 9, and an increase in electrophoretic mobility was evidenced with increasing pH. Colloidal stability was preserved at physiological pH for both non-grafted NPs and grafted NPs in the 10-100 mM ionic strength range and in the 15-60 °C temperature range. A strong influence of the nature of the buffer counterion on NP electrophoretic mobility and peak dispersion was evidenced, thus indicating some interactions between buffer components and NP-aptamer conjugates. Whereas an electric field effect (50-900 V cm(-1)) on NP electrophoretic mobility was evidenced, probably linked to counterion dissociation, temperature seems to have an appreciable effect on the zeta potential and aptamer configuration as well. This information is crucial for estimating the potentialities of such biorecognition tools in electrophoretic systems.

Synapsable quadruplex-mediated fibers

We have fabricated a DNA-based nanofiber created by self-assembly of guanine quadruplex (Hoogsteen base pairing) and double-stranded DNA (Watson-Crick base pairing). When duplexes containing a long stretch of contiguous guanines and single-stranded overhangs are incubated in potassium-containing buffer, the preformed duplexes create high molecular weight species that contain quadruplexes. In addition to observation of these larger species by gel electrophoresis, solutions were analyzed by atomic force microscopy to reveal nanofibers. Analysis of the atomic force microscopy images indicates that fibers form with lengths ranging from 250 to 2,000 nm and heights from 0.45 to 4.0 nm. This work is a first step toward the creation of new structurally heterogeneous (quadruplex/duplex), yet controllable, DNA-based materials exhibiting novel properties suitable for a diverse array of nanotechnology applications.

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.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30δ 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.

Torsional rigidities of weakly strained DNAs

Measurements on unstrained linear and weakly strained large (> or =340 bp) circular DNAs yield torsional rigidities in the range C = 170-230 fJ fm. However, larger values, in the range C = 270-420 fJ fm, are typically obtained from measurements on sufficiently small (< or =247 bp) circular DNAs, and values in the range C = 300-450 fJ fm are obtained from experiments on linear DNAs under tension. A new method is proposed to estimate torsional rigidities of weakly supercoiled circular DNAs. Monte Carlo simulations of the supercoiling free energies of solution DNAs, and also of the structures of surface-confined supercoiled plasmids, were performed using different trial values of C. The results are compared with experimental measurements of the twist energy parameter, E(T), that governs the supercoiling free energy, and also with atomic force microscopy images of surface-confined plasmids. The results clearly demonstrate that C-values in the range 170-230 fJ fm are compatible with experimental observations, whereas values in the range C > or = 269 fJ fm, are incompatible with those same measurements. These results strongly suggest that the secondary structure of DNA is altered by either sufficient coherent bending strain or sufficient tension so as to enhance its torsional rigidity.

Can reliable torsion elastic constants be determined from FPA data on 24 and 27 base-pair DNAs?

Torsion elastic constants obtained from fluorescence polarization anisotropy (FPA) measurements on fifty-three 24 and 27 base-pair (bp) DNAs were recently reported [F. Pedone, F. Mazzei, D. Santoni, Sequence-dependent DNA torsional rigidity: a tetranucleotide code, Biophys. Chem. 112 (2004) 77-88; F. Pedone, F. Mazzei, M. Matzeu, F. Barone, Torsional constant of 27-mer DNA oligomers of different sequences, Biophys. Chem. 94 (2001) 175-184]. The problem of extracting reliable torsion elastic constants (alpha) from FPA measurements on such short DNAs is examined in detail. The difficulty is illustrated by two (fictitious) 24 bp DNAs with approximately 5-fold different torsion elastic constants and 10% different initial anisotropies (r(0)), which exhibit practically indistinguishable anisotropy decays for all t>1 ns. FPA data were simulated for 24 bp DNAs with different input values of alpha and r(0) in the presence and absence of Poisson noise, and were fitted using different choices of the adjustable and fixed parameters. Experimental data for a 24 bp DNA were fitted in a similar manner. For either the simulated or experimental FPA data, it was not possible to determine both the initial anisotropy, r(0), and the torsion elastic constant, alpha, in a reliable (i.e. statistically significant) manner in the presence of Poisson noise. When r(0) is assumed to be fixed at any particular value in the fitting protocol, a unique best-fit value of alpha is obtained, but that best-fit alpha is extremely sensitive to small deviations of the assumed fixed value of r(0) away from the input r(0)-value of the simulated data. Pedone et al. fitted their FPA data by assuming that r(0)=0.360, and adjusting alpha, the hydrodynamic radius (R(H)), and effective length (L). In fact, the reported best-fit values of R(H) and L lay significantly outside their expected ranges. When this same fitting protocol is applied to simulated data for 27 bp DNAs, better overall agreement with the reported experimental values (alpha, R(H), and L) is obtained for a model, wherein all DNAs have the same typical input alpha=5.9 x 10(-12) dyn cm, R(H)=10.0 A, and L=27 (3.4)+2.7=94.5 A, but a 1.00- to 1.13-fold range of r(0)-values, than for the model of Pedone et al., wherein all DNAs have the same input r(0)=0.360, R(H)=10.0 A, and L=94.5 A, but a approximately 3-fold range of alpha-values. It is concluded that, in the absence of reliable independent estimates of r(0) for every DNA, the alpha-values reported for 24 and 27 bp DNAs cannot be regarded as experimentally justified. The reliability of the torsion elastic constants reported for the 136 distinct tetranucleotide steps, which are inferred from the values reported for the fifty-three 24 and 27 bp DNAs, is also briefly discussed.

Sequence-dependent DNA torsional rigidity: a tetranucleotide code

Using fluorescence polarization anisotropy (FPA), we measured the torsional constant of various DNA oligomers in different sequences and calculated the value for each of the 136 unique tetranucleotides. From these values, we obtained a "rigidity profile" for every double-stranded DNA sequence. We tested the code in the analysis of DNA sequences able to form nucleosomes. More than 50% of the sequences studied showed a common 20 and/or 30 bp modulation of the torsional constant. Many other profiles of rigidity were observed in the remaining sequences and this variety in torsional constant modulation may be related to functional differences between nucleosomes.

Thermodynamic treatment of oligonucleotide duplex-simplex equilibria

Thermodynamic formulations have been devised to obtain Delta G degrees values directly from spectroscopic data at a fixed common temperature in nucleic acid duplex-simplex melting curves. In addition, the dependence of melting on salt concentration has been expressed in terms of a stepwise stoichiometric representation, which leads to a specific equation for the partition of the added sodium ions between the different oligonucleotide forms.

A temperature-dependent interaction of neutral red with calf thymus DNA

Neutral red (NR) is used as a probe to study the temperature and concentration dependent interaction of a cationic dye with nucleic acid. A temperature-dependent interaction of NR with calf thymus DNA (CT DNA) has been studied by differential pulse voltammetry (DPV), UV-Visible absorption, circular dichroism (CD) and fluorescence spectroscopy. The experimental results of increasing peak current, changes in the UV-Visible absorption and fluorescence spectra of NR and decreasing the induced circular dichroism (ICD) intensity show that (i) the binding mode of NR molecules is changed from intercalating into DNA base pairs to aggregating along the DNA double helix and (ii) the orientation of NR chromophore in DNA double helix is also changed with the temperature.

Temperature dependence of the Raman spectrum of DNA. II. Raman signatures of premelting and melting transitions of poly(dA).poly(dT) and comparison with poly(dA-dT).poly(dA-dT)

The temperature dependence of the Raman spectrum of poly(dA).poly(dT) (dA: deoxyadenosine; dT: thymidine), a model for DNA containing consecutive adenine.thymine (A.T) pairs, has been analyzed using a spectrometer of high spectral precision and sensitivity. Three temperature intervals are distinguished: (a) premelting (10 < t < 70 degrees C), in which the native double helix is structurally altered but not dissociated into single strands; (b) melting (70 < t < 80 degrees C), in which the duplex is dissociated into single strands; and (c) postmelting (80 < t degrees C), in which no significant structural change can be detected. The distinctive Raman difference signatures observed between 10 and 70 degrees C and between 70 and 80 degrees C are interpreted in terms of the structural changes specific to premelting and melting transitions, respectively. Premelting alters the low-temperature conformation of the deoxyribose-phosphate backbone and eliminates base hydrogen bonding that is distinct from canonical Watson-Crick hydrogen bonding; these premelting perturbations occur without disruption of base stacking. Conversely, melting eliminates canonical Watson-Crick pairing and base stacking. The results are compared with those reported previously on poly(dA-dT).poly(dA-dT), the DNA structure consisting of alternating A.T and T.A pairs (L. Movileanu, J. M. Benevides, and G. J. Thomas, Jr. Journal of Raman Spectroscopy, 1999, Vol. 30, pp. 637-649). Poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) exhibit strikingly dissimilar temperature-dependent Raman profiles prior to the onset of melting. However, the two duplexes exhibit very similar melting transitions, including the same Raman indicators of ruptured Watson-Crick pairing, base unstacking and collapse of backbone order. A detailed analysis of the data provides a comprehensive Raman assignment scheme for adenosine and thymidine residues of B-DNA, delineates Raman markers diagnostic of consecutive A.T and alternating A.T/T.A tracts of DNA, and identifies the distinct Raman difference signatures for premelting and melting transitions in the two types of sequences.

Monte Carlo simulations of supercoiled DNAs confined to a plane

Recent advances in atomic force microscopy (AFM) have enabled researchers to obtain images of supercoiled DNAs deposited on mica surfaces in buffered aqueous milieux. Confining a supercoiled DNA to a plane greatly restricts its configurational freedom, and could conceivably alter certain structural properties, such as its twist and writhe. A program that was originally written to perform Monte Carlo simulations of supercoiled DNAs in solution was modified to include a surface potential. This potential flattens the DNAs to simulate the effect of deposition on a surface. We have simulated transfers of a 3760-basepair supercoiled DNA from solution to a surface in both 161 and 10 mM ionic strength. In both cases, the geometric and thermodynamic properties of the supercoiled DNAs on the surface differ significantly from the corresponding quantities in solution. At 161 mM ionic strength, the writhe/twist ratio is 1.20-1.33 times larger for DNAs on the surface than for DNAs in solution and significant differences in the radii of gyration are also observed. Simulated surface structures in 161 mM ionic strength closely resemble those observed by AFM. Simulated surface structures in 10 mM ionic strength are similar to a minority of the structures observed by AFM, but differ from the majority of such structures for unknown reasons. In 161 mM ionic strength, the internal energy (excluding the surface potential) decreases substantially as the DNA is confined to the surface. Evidently, supercoiled DNAs in solution are typically deformed farther from the minimum energy configuration than are the corresponding surface-confined DNAs. Nevertheless, the work (Delta A(int)) done on the internal coordinates, which include uniform rotations at constant configuration, during the transfer is positive and 2.6-fold larger than the decrease in internal energy. The corresponding entropy change is negative, and its contribution to Delta A(int) is positive and exceeds the decrease in internal energy by 3.6 fold. The work done on the internal coordinates during the solution-to-surface transfer is directed primarily toward reducing their entropy. Evidently, the number of configurations available to the more deformed solution DNA is vastly greater than for the less deformed surface-confined DNA.

Torsional constant of 27-mer DNA oligomers of different sequences

We have studied the torsional elastic constant (alpha) of short DNA (27mer) oligomers of various sequence by fluorescence polarization anysotropy (FPA) measurements. The lowest alpha values were found in samples with sequence rich in AA dinucleotides or containing the alternating d(A-T) x d(A-T) motif. The torsional rigidity of our DNA samples was compared to that calculated according to the current values of twist angle fluctuations derived for ten dinucleotide steps by recent analyses of DNA crystal structure database. The values of torsional rigidity derived from crystals are higher than our experimental ones, obtained by FPA analysis, suggesting that packing force in crystals may notably hinder the dinucleotide twist angle fluctuations that occur in solution. This behaviour is more evident for samples containing AA, TA and AT steps. In all the samples there is about a twofold change of the alpha value in the 10-40 degrees C range. An activation enthalpy (Delta H (#)) of about 17.4 kJ mol(-1), on average, was obtained for the temperature dependence of eight of the ten samples studied. A correlation with the stacking energy is discussed.

Effect of polyethylene glycol on the supercoiling free energy of DNA

The supercoiling free energy of pUC19 DNA [2686 base pairs (bp)] was measured in various concentrations of PEG 8000 (polyethylene glycol; molecular weight 8000) by the topoisomer distribution method. The effective twist energy parameter (E(T)) that governs the supercoiling free energy declined linearly by 1.9-fold with increasing w/v % PEG from 0 to 7.5%, which lies below the threshold for intermolecular condensation. In principle, PEG could affect E(T) either via an osmotic exclusion mechanism or by altering the torsion elastic constant, bending rigidity, or self-repulsions of the DNA. Possible alterations of the DNA secondary structure and torsion elastic constant were assessed by CD spectroscopy and time-resolved fluorescence polarization anisotropy of intercalated ethidium. Up to 7.5% PEG, the secondary structure of the DNA remained largely unaltered, as evidenced by (1) the absence of any significant change in the CD spectrum, (2) an extremely small relative decrease (-0.0013) in intrinsic twist, and (3) a negligibly small change in the torsion elastic constant. The observed reduction in E(T) cannot be ascribed primarily to a decrease in torsion elastic constant, and most likely does not stem from a decrease in bending rigidity either. The decrease in medium dielectric constant due to PEG should increase the self-repulsions, and thereby increase E(T), which is opposite to the observed trend. Instead, the observed decline in E(T) is attributed to an osmotic exclusion mechanism. The change in molar volume excluded to the PEG (Delta V(ex)), when the linking difference converts from Delta l = 0 to Delta l = +/-1, was determined from the observed E(T) value and PEG osmotic pressure at each concentration. The experimental Delta V(ex) values agree well with theoretical estimates reckoned for a simple osmotic exclusion model, in which PEG is excluded by hard-core interactions from a concentric cylindrical volume around every duplex segment. The difference in volume excluded to PEG between the Delta l = 0 and the Delta l = +/-1 topoisomers is attributed entirely to the approximately 0.7 additional writhe "crossing" of two duplex strands at roughly 90 degrees, which is known to occur in the latter species. When the separation between the duplex centers at the "crossing" was adjusted so that the theoretical estimate of Delta V(ex) matched the experimental value at each PEG concentration, a value near 5.7 nm was obtained in each case. The invariance and plausible magnitude of this mean separation at the crossing provide strong support for this simple osmotic exclusion model. An alternative model, in which the PEG is excluded from the entire coil envelope of the DNA out to its radius of gyration, perhaps because it decreases the local dielectric constant, was also considered. The estimated difference in excluded volume in that case exceeds the experimental value by a factor of nearly 10(4), and could be ruled out on that basis.

Conformational deformability of RNA: a harmonic mode analysis

The harmonic mode analysis method was used to characterize the conformational deformability of regular Watson-Crick paired, mismatch- and bulge-containing RNA. Good agreement between atomic Debye-Waller factors derived from x-ray crystallography of a regular RNA oligonucleotide and calculated atomic fluctuations was obtained. Calculated helical coordinate fluctuations showed a small sequence dependence of up to approximately 30-50%. A negative correlation between motions at a given base pair step and neighboring steps was found for most helical coordinates. Only very few calculated modes contribute significantly to global motions such as bending, twisting, and stretching of the RNA molecules. With respect to a local helical description of the RNA helix our calculations suggest that RNA bending is mostly due to a periodic change in the base pair step descriptors slide and roll. The presence of single guanine:uridine or guanine:adenine mismatches had little influence on the calculated RNA flexibility. In contrast, for tandem guanine:adenine base pairs the harmonic mode approach predicts a significantly reduced conformational flexibility in the case of a sheared arrangement and slightly enhanced flexibility for a face-to-face (imino proton) pairing relative to regular RNA. The presence of a single extra adenine bulge nucleotide stacked between flanking sequences resulted in an increased local atomic mobility around the bulge site (approximately 40%) and a slightly enhanced global bending flexibility. For an adenine bulge nucleotide in a looped-out conformation a strongly enhanced bulge nucleotide mobility but no increased bending flexibility compared to regular RNA was found.

Dynamic twisting correlations in a model DNA with uniform torsion elastic constant

The extent to which the twisting motions of two separate subunits (base pairs) in an elastic filament are correlated is discussed in terms of the two-point correlation function of their azimuthal angular displacements C(Delta,t) identical with <(phi(m)(t) - phi(m)(0))(phi(n)(t) - phi(n)(0))>, where m,n are the sequential subunit indices and Delta = |m - n| is their absolute difference. An approximate expression is derived for C(Delta,t) for an infinitely long model DNA from the analytical theory developed previously to treat the decay of the fluorescence polarization anisotropy of intercalated ethidium. C(Delta,t) is numerically evaluated as a function of Delta for a range of times (20, 40, 60, 80, 100, and 120 ns) for a model DNA with a typical torsion elastic constant. By t = 120 ns, significant dynamic correlations are observed to extend over a domain (Delta) several hundred base pairs. Copyright 1999 John Wiley & Sons, Inc.