Competitive electrostatic binding of charged ligands to polyelectrolytes: practical approach using the non-linear Poisson-Boltzmann equation

We have developed a practical analytical treatment of the non-linear Poisson-Boltzmann (P-B) equation to characterize the strong but non-specific binding of charged ligands to DNA and other highly charged macromolecules. These reactions are notable for their strong salt dependence and anti-cooperativity, features which the theory fully explains. We summarize analytical results for concentration profiles and ion binding in various regimes of surface curvature and ionic strength, and show how counterion size and charge distribution may influence competitive binding. We present several practical applications of the formalism, showing how to estimate the ligand concentration needed to effectively compete with a given buffer salt, and how to calculate the amounts of counterion species bound at various distances from the DNA surface under given bulk solution conditions. We cast our results into the form of a Scatchard binding isotherm, showing how the apparent binding constant K(obs) and S = -dlog K (obs )d log[M (+)] can be predicted from the basic theory. Anti-cooperativity arises naturally without steric repulsion, and binding curves can be fitted with K(obs) and effective charge as the only free parameters. We extend the analytical P-B analysis to an arbitrary number of counterion species, and apply the results to fit and predict three-ion competition data.

Thermodynamics of the hydrophobic effect. II. Calorimetric measurement of enthalpy, entropy, and heat capacity of aggregation of alkylamines and long aliphatic chains

The thermodynamics of long aliphatic chain alkylamine aggregation in aqueous solution was studied by isothermal titration calorimetry (ITC). Protonated alkylammonium cations with linear aliphatic chains of 10-14 carbon atoms were fully soluble in aqueous solution at the beginning of titration, but practically insoluble after deprotonation by titrating with sodium hydroxide. The alkylamines aggregated and precipitated during the reaction, enabling direct measurement of the enthalpy of aggregation. The enthalpy of aggregation became increasingly exothermic upon increasing the chain length. Hydrophobic aggregation was enthalpy-driven and entropy-opposed for alkylamines with 12-14 carbon atoms at room temperature. Direct observation of hydrophobic aggregation by ITC at constant temperature and pressure provided more accurate thermodynamic parameters than obtainable from van't Hoff analysis. Aggregation into liquid or solid phases could be distinguished by ITC, but not by van't Hoff analysis of alkylamine solubility data.

Thermodynamics of the hydrophobic effect. I. Coupling of aggregation and pK(a) shifts in solutions of aliphatic amines

Long aliphatic hydrocarbon chains aggregate in aqueous solution due to the hydrophobic effect, forming structures such as micelles and membranes, while amino groups titrate at basic pH. These two biologically important behaviors are linked in alkylamines, in which the pK(a) of the amino group is shifted downward by aggregation. In this paper we study the thermodynamics of these coupled processes, following aggregation by observing alkylamine pH titration behavior. The magnitude of the shift depended on the aliphatic chain length and on the concentration of alkylamine: longer chains and higher concentrations lowered the pK(a) to a greater extent. Gibbs free energies of protonation and aggregation were calculated from the pK(a) shifts. Enthalpies, entropies, and heat capacities were estimated by van't Hoff analysis from the pK(a) shift dependencies on temperature. However, the results were less precise than the calorimetrically measured values, as described in the following article. A model to calculate titration curves, pK(a) shifts, and aggregation of uncharged alkylamines as a function of aliphatic chain length, concentration, and temperature is presented.

Hexamminecobalt(III)-induced condensation of calf thymus DNA: circular dichroism and hydration measurements

The interaction of hexamminecobalt(III), Co(NH(3))(6)(3+), with 160 and 3000-8000 bp length calf thymus DNA has been investigated by circular dichroism, acoustic and densimetric techniques. The acoustic titration curves of 160 bp DNA revealed three stages of interaction: (i) Co(NH(3))(6)(3+) binding up to the molar ratio [Co(NH(3))(6)(3+)]/[P] = 0.25, prior to DNA condensation; (ii) a condensation process between [Co(NH(3))(6)(3+)]/[P] = 0.25 and 0.30; and (iii) precipitation after [Co(NH(3))(6)(3+)]/[P] = 0.3. In the case of 3000-8000 bp DNA only two processes were observed: (i) binding up to [Co(NH(3))(6)(3+)]/[P] = 0.3; and (ii) precipitation after this point. In agreement with earlier observations, long DNA aggregates without changes in its B-form circular dichroism spectrum, while short DNA demonstrates a positive B-->Psi transition after [Co(NH(3))(6)(3+)]/[P] = 0.25. From ultrasonic and densimetric measurements the effects of Co(NH(3))(6)(3+) binding on volume and compressibility have been obtained. The binding of Co(NH(3))(6)(3+) to both short and long DNA is characterized by similar changes in volume and compressibility calculated per mole Co(NH(3))(6)(3+): DeltaV = 9 cm(3) mol(-1) and Deltakappa = 33 x 10(-4) cm(3) mol(-1) bar(-1). The positive sign of the parameters indicates dehydration, i.e. water release from Co(NH(3))(6)(3+) and the atomic groups of DNA. This extent of water displacement would be consistent with the formation of two direct, hydrogen bonded contacts between the cation and the phosphates of DNA.

Mechanism for nucleic acid chaperone activity of HIV-1 nucleocapsid protein revealed by single molecule stretching

The nucleocapsid protein (NC) of HIV type 1 is a nucleic acid chaperone that facilitates the rearrangement of nucleic acids into conformations containing the maximum number of complementary base pairs. We use an optical tweezers instrument to stretch single DNA molecules from the helix to coil state at room temperature in the presence of NC and a mutant form (SSHS NC) that lacks the two zinc finger structures present in NC. Although both NC and SSHS NC facilitate annealing of complementary strands through electrostatic attraction, only NC destabilizes the helical form of DNA and reduces the cooperativity of the helix-coil transition. In particular, we find that the helix-coil transition free energy at room temperature is significantly reduced in the presence of NC. Thus, upon NC binding, it is likely that thermodynamic fluctuations cause continuous melting and reannealing of base pairs so that DNA strands are able to rapidly sample configurations to find the lowest energy state. The reduced cooperativity allows these fluctuations to occur in the middle of complex double-stranded structures. The reduced stability and cooperativity, coupled with the electrostatic attraction generated by the high charge density of NC, is responsible for the nucleic acid chaperone activity of this protein.

Entropy and heat capacity of DNA melting from temperature dependence of single molecule stretching

When a single molecule of double-stranded DNA is stretched beyond its B-form contour length, the measured force shows a highly cooperative overstretching transition. We have measured the force at which this transition occurs as a function of temperature. To do this, single molecules of DNA were captured between two polystyrene beads in an optical tweezers apparatus. As the temperature of the solution surrounding a captured molecule was increased from 11 degrees C to 52 degrees C in 500 mM NaCl, the overstretching transition force decreased from 69 pN to 50 pN. This reduction is attributed to a decrease in the stability of the DNA double helix with increasing temperature. These results quantitatively agree with a model that asserts that DNA melting occurs during the overstretching transition. With this model, the data may be analyzed to obtain the change in the melting entropy DeltaS of DNA with temperature. The observed nonlinear temperature dependence of DeltaS is a result of the positive change in heat capacity of DNA upon melting, which we determine from our stretching measurements to be DeltaC(p) = 60 +/- 10 cal/mol K bp, in agreement with calorimetric measurements.

Force-induced melting of the DNA double helix. 2. Effect of solution conditions

In this paper, we consider the implications of the general theory developed in the accompanying paper, to interpret experiments on DNA overstretching that involve variables such as solution temperature, pH, and ionic strength. We find the DNA helix-coil phase boundary in the force-temperature space. At temperatures significantly below the regular (zero force) DNA melting temperature, the overstretching force, f(ov)(T), is predicted to decrease nearly linearly with temperature. We calculate the slope of this dependence as a function of entropy and heat-capacity changes upon DNA melting. Fitting of the experimental f(ov)(T) dependence allows determination of both of these quantities in very good agreement with their calorimetric values. At temperatures slightly above the regular DNA melting temperature, we predict stabilization of dsDNA by moderate forces, and destabilization by higher forces. Thus the DNA stretching curves, f(b), should exhibit two rather than one overstretching transitions: from single stranded (ss) to double stranded (ds) and then back at the higher force. We also predict that any change in DNA solution conditions that affects its melting temperature should have a similar effect on DNA overstretching force. This result is used to calculate the dependence of DNA overstretching force on solution pH, f(ov)(pH), from the known dependence of DNA melting temperature on pH. The calculated f(ov)(pH) is in excellent agreement with its experimental determination (M. C. Williams, J. R. Wenner, I. Rouzina, and V. A. Bloomfield, Biophys. J., accepted for publication). Finally, we quantitatively explain the measured dependence of DNA overstretching force on solution ionic strength for crosslinked and noncrosslinked DNA. The much stronger salt dependence of f(ov) in noncrosslinked DNA results from its lower linear charge density in the melted state, compared to crosslinked or double-stranded overstretched S-DNA.

Force-induced melting of the DNA double helix 1. Thermodynamic analysis

The highly cooperative elongation of a single B-DNA molecule to almost twice its contour length upon application of a stretching force is interpreted as force-induced DNA melting. This interpretation is based on the similarity between experimental and calculated stretching profiles, when the force-dependent free energy of melting is obtained directly from the experimental force versus extension curves of double- and single-stranded DNA. The high cooperativity of the overstretching transition is consistent with a melting interpretation. The ability of nicked DNA to withstand forces greater than that at the transition midpoint is explained as a result of the one-dimensional nature of the melting transition, which leads to alternating zones of melted and unmelted DNA even substantially above the melting midpoint. We discuss the relationship between force-induced melting and the B-to-S transition suggested by other authors. The recently measured effect on T7 DNA polymerase activity of the force applied to a ssDNA template is interpreted in terms of preferential stabilization of dsDNA by weak forces approximately equal to 7 pN.

Effect of pH on the overstretching transition of double-stranded DNA: evidence of force-induced DNA melting

When a single molecule of double-stranded DNA is stretched beyond its B-form contour length, the measured force shows a highly cooperative overstretching transition. We have investigated the source of this transition by altering helix stability with solution pH. As solution pH was increased from pH 6.0 to pH 10.6 in 250 mM NaCl, the overstretching transition force decreased from 67.0 +/- 0.8 pN to 56.2 +/- 0.8 pN, whereas the transition width remained nearly constant. As the pH was lowered from pH 6.0 to pH 3.1, the overstretching force decreased from 67.0 +/- 0.8 pN to 47.0 +/- 1.0 pN, but the transition width increased from 3.0 +/- 0.6 pN to 16.0 +/- 3 pN. These results quantitatively agree with a model that asserts that DNA strand dissociation, or melting, occurs during the overstretching transition.

The greater negative charge density of DNA in tris-borate buffers does not enhance DNA condensation by multivalent cations

As indicated by recent measurements of the electrophoretic free solution mobility, DNA appears to have a greater helical charge density in Tris-borate-EDTA (TBE) buffers than in Tris-acetate-EDTA (TAE) buffers. Since electrostatic forces play a major role in DNA packaging processes, we have investigated the condensation of closed circular plasmid DNA using total intensity and dynamic light scattering in Tris-borate, Tris-acetate, and Tris-cacodylate buffers with cobaltic hexa-amine (III) [Co(NH(3))(3+)(6)]. We find that neither the critical concentration of Co(NH(3))(3+)(6) nor the hydrodynamic radii of the resulting condensates vary significantly in the buffer systems studied here despite the prediction that DNA condensation should occur at significantly lower Co(NH(3))(3+)(6) concentrations in Tris-borate buffers. Assuming a persistence length behavior similar to B-DNA in the presence of multivalent cations, a decrease in the attractive counterion correlation pressure decay length in Tris-borate buffers does not account for our observations. It is possible that the binding of multivalent cations to DNA may hinder borate association with the DNA double helix.

Excluded volume in solvation: sensitivity of scaled-particle theory to solvent size and density

Changes in solvent environment greatly affect macromolecular structure and stability. To investigate the role of excluded volume in solvation, scaled-particle theory is often used to calculate delta G(tr)(ev), the excluded-volume portion of the solute transfer free energy, delta G(tr). The inputs to SPT are the solvent radii and molarities. Real molecules are not spheres. Hence, molecular radii are not uniquely defined and vary for any given species. Since delta G(tr)(ev) is extremely sensitive to solvent radii, uncertainty in these radii causes a large uncertainty in delta G(tr)(ev)-several kcal/mol for amino acid solutes transferring from water to aqueous mixtures. This uncertainty is larger than the experimental delta G(tr) values. Also, delta G(tr)(ev) can be either positive or negative. Adding neutral crowding molecules may not necessarily reduce solubility. Lastly, delta G(tr)(ev) is very sensitive to solvent density, rho. A few percent error in rho may even cause qualitative deviations in delta G(tr)(ev). For example, if rho is calculated by assuming the hard-sphere pressure to be constant, then delta G(tr)(ev) values and uncertainties are now only tenths of a kcal/mol and are positive. Because delta G(tr)(ev) values calculated by scaled-particle theory are strongly sensitive to solvent radii and densities, determining the excluded-volume contribution to transfer free energies using SPT may be problematic.

Static and dynamic light scattering from aggregating particles

We use standard hydrodynamic and light scattering theories to calculate the total intensity and dynamic light scattering properties of random aggregates of spherical particles containing up to ten spheres. When the aggregates have dimensions comparable to the wavelength of light, intraaggregate interference effects can dramatically reduce the apparent size of the aggregates. These results could be significant in interpreting DNA condensation, protein polymerization, and other biomolecular aggregation reactions.

Structural basis of polyamine-DNA recognition: spermidine and spermine interactions with genomic B-DNAs of different GC content probed by Raman spectroscopy

Four genomic DNAs of differing GC content (Micrococcus luteus, 72% GC; Escherichia coli, 50% GC; calf thymus, 42% GC; Clostridium perfringens, 27% GC) have been employed as targets of interaction by the cationic polyamines spermidine ([H(3)N(CH(2))(3)NH(2)(CH(2))(4)NH(3)](3+)) and spermine ([(CH(2))(4)(NH(2)(CH(2))(3)NH(3))(2)](4+)). In solutions containing 60 mM DNA phosphate (approximately 20 mg DNA/ml) and either 1, 5 or 60 mM polyamine, only Raman bands associated with the phosphates exhibit large spectral changes, demonstrating that B-DNA phosphates are the primary targets of interaction. Phosphate perturbations, which are independent of base composition, are consistent with a model of non-specific cation binding in which delocalized polyamines diffuse along DNA while confined by the strong electrostatic potential gradient perpendicular to the helix axis. This finding provides experimental support for models in which polyamine-induced DNA condensation is driven by non-specific electrostatic binding. The Raman spectra also demonstrate that major groove sites (guanine N7 and thymine C5H(3)) are less affected than phosphates by polyamine-DNA interactions. Modest dependence of polyamine binding on genome base composition suggests that sequence context plays only a secondary role in recognition. Importantly, the results demonstrate that polyamine binding has a negligible effect on the native B-form secondary structure. The capability of spermidine or spermine to bind and condense genomic B-DNA without disrupting the native structure must be taken into account when considering DNA organization within bacterial nucleoids or cell nuclei.

DNA condensation by multivalent cations

In the presence of multivalent cations, high molecular weight DNA undergoes a dramatic condensation to a compact, usually highly ordered toroidal structure. This review begins with an overview of DNA condensation: condensing agents, morphology, kinetics, and reversibility, and the minimum size required to form orderly condensates. It then summarizes the statistical mechanics of the collapse of stiff polymers, which shows why DNA condensation is abrupt and why toroids are favored structures. Various ways to estimate or measure intermolecular forces in DNA condensation are discussed, all of them agreeing that the free energy change per base pair is very small, on the order of 1% of thermal energy. Experimental evidence is surveyed showing that DNA condensation occurs when about 90% of its charge is neutralized by counterions. The various intermolecular forces whose interplay gives rise to DNA condensation are then reviewed. The entropy loss upon collapse of the expanded wormlike coil costs free energy, and stiffness sets limits on tight curvature. However, the dominant contributions seem to come from ions and water. Electrostatic repulsions must be overcome by high salt concentrations or by the correlated fluctuations of territorially bound multivalent cations. Hydration must be adjusted to allow a cooperative accommodation of the water structure surrounding surface groups on the DNA helices as they approach. Undulations of the DNA in its confined surroundings extend the range of the electrostatic forces. The condensing ions may also subtly modify the local structure of the double helix.

Condensation of DNA by multivalent cations: experimental studies of condensation kinetics

DNA in viruses and cells exists in highly condensed, tightly packaged states. We have undertaken an in vitro study of the kinetics of DNA condensation by the trivalent cation hexaammine cobalt (III) with the aim of formulating a quantitative, mechanistic model of the condensation process. Experimental approaches included total intensity and dynamic light scattering, electron microscopy, and differential sedimentation. We determined the average degree of condensation, the distribution of condensate sizes, and the fraction of uncondensed DNA as a function of reaction time for a range of [DNA] and [Co(NH(3))(3+)(6)]. We find the following: (1) DNA condensation occurs only above a critical [Co(NH(3))(3+)(6)] for a given DNA and salt concentration. At the onset of condensation, [Co(NH(3))(3+)(6)]/[DNA-phosphate] is close to the average value of 0.54, which reflects the 89-90% charge neutralization criterion for condensation. (2) The equilibrium weight average hydrodynamic radius <R(H) > of the condensates first decreases, then increases with increasing [Co(NH(3))(3+)(6)] as they undergo a transition from intramolecular (monomolecular) to intermolecular (multimolecular) condensation. However, <R(H) > is insensitive to [DNA]. (3) The uncondensed DNA fraction decays approximately exponentially with time. The equilibrium uncondensed DNA fraction and relaxation time decrease with increasing [Co(NH(3))(3+)(6)] but are insensitive to [DNA]. (4) The condensation rate in its early stages is insensitive to [DNA] but proportional to [Co(NH(3))(3+)(6)](xs) = [Co(NH(3))(3+)(6)] - [Co(NH(3))(3+)(6)](crit). (5) Data for low [DNA] and low [Co(NH(3))(3+)(6)] at early stages of condensation are most reliable for kinetic modeling since under these conditions there is minimal clumping and network formation among separate condensates. A mechanism with initial monomolecular nucleation and subsequent bimolecular association and unimolecular dissociation steps with rate constants that depend on the number of DNA molecules in the condensate, accounts reasonably well for these observations.

Stretching of single collapsed DNA molecules

The elastic response of single plasmid and lambda phage DNA molecules was probed using optical tweezers at concentrations of trivalent cations that provoked DNA condensation in bulk. For uncondensed plasmids, the persistence length, P, decreased with increasing spermidine concentration before reaching a limiting value 40 nm. When condensed plasmids were stretched, two types of behavior were observed: a stick-release pattern and a plateau at approximately 20 pN. These behaviors are attributed to unpacking from a condensed structure, such as coiled DNA. Similarly, condensing concentrations of hexaammine cobalt(III) (CoHex) and spermidine induced extensive changes in the low and high force elasticity of lambda DNA. The high force (5-15 pN) entropic elasticity showed worm-like chain (WLC) behavior, with P two- to fivefold lower than in low monovalent salt. At lower forces, a 14-pN plateau abruptly appeared. This corresponds to an intramolecular attraction of 0.083-0.33 kT/bp, consistent with osmotic stress measurements in bulk condensed DNA. The intramolecular attractive force with CoHex is larger than with spermidine, consistent with the greater efficiency with which CoHex condenses DNA in bulk. The transition from WLC behavior to condensation occurs at an extension about 85% of the contour length, permitting looping and nucleation of condensation. Approximately half as many base pairs are required to nucleate collapse in a stretched chain when CoHex is the condensing agent.

Linkage analysis of chromosome 2q in osteoarthritis

BACKGROUND

In independent linkage studies chromosome 2q11-q24 and chromosome 2q23-35 have previously been implicated as regions potentially harbouring susceptibility loci for osteoarthritis (OA).

OBJECTIVE

To test chromosome 2q for linkage to idiopathic osteoarthritis.

METHODS

Using a cohort of 481 OA families that each contained at least one affected sibling pair with severe end-stage disease (ascertained by hip or knee joint replacement surgery), we conducted a linkage analysis of chromosome 2q using 16 polymorphic microsatellite markers at an average spacing of one marker every 8.5 cM.

RESULTS

Our results provide suggestive evidence for a locus at 2q31 with a maximum multipoint logarithm of the odds score (MLS) of 1.22 which increased to 2.19 in those families concordant for hip-only disease (n = 311). This suggestive linkage was greater in male-hip families (MLS = 1.57, n = 71) than in female-hip families (MLS = 0.71, n = 132).

CONCLUSIONS

Chromosome 2q is likely to contain at least one susceptibility locus for OA.

Thermodynamics of DNA binding and condensation: isothermal titration calorimetry and electrostatic mechanism

The thermodynamics of binding of the trivalent cations cobalt hexammine and spermidine to plasmid DNA was studied by isothermal titration calorimetry. Two stages were observed in the course of titration, the first attributed to cation binding and the second to DNA condensation. A standard calorimetric data analysis was extended by applying an electrostatic binding model, which accounted for most of the observed data. Both the binding and condensation reactions were entropically driven (TDeltaS approximately +10 kcal/mol cation) and enthalpically opposed (DeltaH approximately +1 kcal/mol cation). As predicted from their relative sizes, the binding constants of the cations were indistinguishable, but cobalt hexammine had a much greater DNA condensing capacity because it is more compact than spermidine. The dependence of both the free energy of cobalt hexammine binding and the critical cobalt hexammine concentration for DNA condensation on temperature and monovalent cation concentration followed the electrostatic model quite precisely. The heat capacity changes of both stages were positive, perhaps reflecting both the temperature dependence of the dielectric constant of water and the burial of polar surfaces. DNA condensation occurred when about 67 % of the DNA phosphate charge was neutralized by cobalt hexammine and 87 % by spermidine. During condensation, the remaining DNA charge was neutralized.

Linkage analysis of candidate genes as susceptibility loci for osteoarthritis-suggestive linkage of COL9A1 to female hip osteoarthritis

OBJECTIVE

To examine 11 candidate genes as susceptibility loci for osteoarthritis (OA).

METHODS

A total of 481 families have been ascertained in which at least two siblings have had joint replacement surgery of the hip, or knee, or hip and knee for idiopathic OA. Each candidate gene was targeted using one or more intragenic or closely linked microsatellite marker. The linkage data were analysed unstratified and following stratification by sex and by joint replaced (hip or knee).

RESULTS

The analyses revealed suggestive linkage of the type IX collagen gene COL9A1 (6q12-q13) to a subset of 132 families that contained affected females who were concordant for hip OA (female-hip) with a P-value of 0.00053 and logarithm of the odds (LOD) score of 2.33 [corrected P-value of 0. 0016, corrected LOD score of 1.85].

CONCLUSIONS

COL9A1 may therefore be a susceptibility locus for female hip OA. In addition, there was weak evidence of linkage to HLA/COL11A2 (6p21.3) in female hip OA with a corrected P-value of 0.016.

Association analysis of the vitamin D receptor gene, the type I collagen gene COL1A1, and the estrogen receptor gene in idiopathic osteoarthritis

OBJECTIVE

Evidence has accumulated supporting a role for genes in the etiology of osteoarthritis (OA). Several candidates have been targeted as potential susceptibility loci including genes that are involved in the regulation of bone density. Genetic association analysis has suggested a role for the vitamin D receptor gene (VDR) and the estrogen receptor gene (ER) in susceptibility. Such findings must be tested in additional independent cohorts. We tested for association of these 2 genes, plus a third gene implicated in bone density, COL1A1, with idiopathic OA.

METHODS

A case-control cohort of 371 affected probands and 369 unaffected spouses was used. Association was tested using 4 intragenic single nucleotide polymorphisms (SNP), one each for the VDR and COL1A1 genes, and 2 for the ER gene. The VDR and ER SNP are the same SNP that have been associated with OA. All 4 SNP affect restriction enzyme sites and were genotyped using polymerase chain reaction and enzyme digestion. Allele and genotype distributions for each SNP were compared between cases and controls and analyzed using Fisher's exact test.

RESULTS

There was no evidence of association of the VDR or the ER gene SNP to OA. There was weak evidence of association of the COL1A1 SNP in female cases (p = 0.017), reflected by a difference in the distribution of genotypes at this SNP between female cases and controls (p = 0.027). However, when corrected for multiple testing, these results were not significant.

CONCLUSION

If the VDR, ER, or COL1A1 genes do encode predisposition to OA then the 4 SNP tested are not associated with major susceptibility alleles at these 3 loci.

Heat capacity effects on the melting of DNA. 1. General aspects

In this paper we analyze published data on DeltaH and DeltaS values for the DNA melting transition under various conditions. We show that there is a significant heat capacity increase DeltaC(p) associated with DNA melting, in the range of 40-100 cal/mol K per base pair. This is larger than the transition entropy per base pair, DeltaS(0) approximately 25 cal/mol K. The ratio of DeltaC(p)/DeltaS(0) determines the importance of heat capacity effects on melting. For DNA this ratio is 2-4, larger than for many proteins. We discuss how DeltaC(p) values can be extracted from experimental data on the dependence of DeltaH and DeltaS on the melting temperature T(m). We consider studies of DNA melting as a function of ionic strength and show that while polyelectrolyte theory provides a good description of the dependence of T(m) on salt, electrostatics alone cannot explain the accompanying strong variation of DeltaH and DeltaS. While T(m) is only weakly affected by DeltaC(p), its dependence on one parameter (e.g., salt) as a function of another (e.g., DNA composition) is determined by DeltaC(p). We show how this accounts for the stronger stabilization of AT relative to GC base pairs with increasing ionic strength. We analyze the source of discrepancies in DeltaH as determined by calorimetry and van't Hoff analysis and discuss ways of analyzing data that yield valid van't Hoff DeltaH. Finally, we define a standard state for DNA melting, the temperature at which thermal contributions to DeltaH and DeltaS vanish, by analyzing experimental data over a broad range of stabilities.

Osmotic pressure effects on EcoRV cleavage and binding

Investigations of DNA binding proteins frequently measure pH and salt dependence, but relatively few studies measure protein binding in high concentrations of small molecules often found in vivo. We have measured kinetics of the restriction enzyme EcoRV in concentrated solutions of three small cosolvents that produce osmotic pressures from 0.25 to 2.5 mol/kg (6 to 62 atm or water activity of 0.995 to 0.956). We have correlated DNA cleavage and binding parameters with four solution parameters (dielectric constant, viscosity, water concentration, and water activity). We found that the responses of maximum velocity (Vmax) and the dissociation constant for nonspecific binding (Kd,ns) best correlate with water activity. The Michaelis constant (Km) correlates with both water activity and solution viscosity, the latter due to the highly dilute reactant concentrations, which make enzyme-substrate combination diffusion limited. Dielectric constant does not influence any of the kinetic parameters, which is consistent with a view that protein and DNA are preferentially hydrated, and excluded solutes cannot affect the local dielectric constant.

Crowding effects on EcoRV kinetics and binding

The cytosol of the cell contains high concentrations of small and large macromolecules, but experimental data are often obtained in dilute solutions that do not reflect in vivo conditions. We have studied the crowding effect that large macromolecules have on EcoRV cleavage by adding high-molecular-weight Ficoll 70 to reaction solutions. Results indicate that Ficoll has surprisingly little effect on overall EcoRV reaction velocity because of offsetting increases in V(max) and K(m), and stronger nonspecific binding. The changes in measured parameters can largely be attributed to the excluded volume effects on reactant activities and the slowing of protein diffusion. Covolume reduction upon binding appears to reinforce nonspecific binding strength, and k(cat) appears to be slowed by stronger nonspecific binding, which slows product release. The data also suggest that effective Ficoll particle volume decreases as its concentration increases above a few weight percent, which may be due to Ficoll interpenetration or compression.

Heat capacity effects on the melting of DNA. 2. Analysis of nearest-neighbor base pair effects

The stability of a DNA double helix of any particular sequence is conventionally estimated as the average of the stabilities of the 10 different nearest-neighbor (NN) base pair doublets that it contains. Therefore, much effort has been devoted to the experimental characterization and tabulation of the enthalpy, entropy, and free energy of melting for each of the NN doublets. Although data from different research groups generally agree for the NN free energies and melting temperatures, there are major disagreements for the enthalpies and entropies. The largest differences are between the parameters obtained on oligomeric relative to polymeric DNA. This disagreement interferes with the practical application of NN thermodynamic parameters. It also raises doubts regarding several fundamental assumptions about DNA melting, such as the absence of longer range interactions, the length dependence of DNA melting parameters per base pair, the applicability of polyelectrolyte theory to the description of salt effects on oligomers, and the purely enthalpic difference between NN doublets. Here we show that if one takes into account the significant heat capacity increase associated with DNA melting, all of the above assumptions are self-consistently reconciled with experiment.

Dependence of the Raman signature of genomic B-DNA on nucleotide base sequence

The vibrational spectra of four genomic and two synthetic DNAs, encompassing a wide range in base composition [poly(dA-dT). poly(dA-dT), 0% G + C; Clostridium perfringens DNA, 27% G + C; calf thymus DNA, 42% G + C; Escherichia coli DNA, 50% G + C; Micrococcus luteus DNA, 72% G + C; poly(dG-dC).poly(dG-dC), 100% G + C] (dA: deoxyadenosine; dG: deoxyguanosine; dC: deoxycytidine; dT: thymidine), have been analyzed using Raman difference methods of high sensitivity. The results show that the Raman signature of B DNA depends in detail upon both genomic base composition and sequence. Raman bands assigned to vibrational modes of the deoxyribose-phosphate backbone are among the most sensitive to base sequence, indicating that within the B family of conformations major differences occur in the backbone geometry of AT- and GC-rich domains. Raman bands assigned to in-plane vibrations of the purine and pyrimidine bases-particularly of A and T-exhibit large deviations from the patterns expected for random base distributions, establishing that Raman hypochromic effects in genomic DNA are also highly sequence dependent. The present study provides a basis for future use of Raman spectroscopy to analyze sequence-specific DNA-ligand interactions. The demonstration of sequence dependency in the Raman spectrum of genomic B DNA also implies the capability to distinguish genomic DNAs by means of their characteristic Raman signatures.

Structural effects of cobalt-amine compounds on DNA condensation

Light scattering and electron microscopy have been used to investigate the structural effects of the trivalent complexes hexaammine cobalt (III) chloride (Cohex), tris(ethylenediamine) cobalt(III) chloride (Coen), and cobalt(III) sepulchrate chloride (Cosep) on DNA condensation. These cobalt-amine compounds have similar ligand coordination geometries but differ slightly in size. Their hydrophobicity is in the order Cosep > Coen > Cohex, according to the numbers of methylene groups in these ligands. All of these compounds effectively precipitate DNA at high concentrations; but despite a lower surface charge density, Cosep condenses DNA twice as effectively as Coen or Cohex. UV and CD measurements of the supernatants of cobalt-amine/DNA solutions reveal a preferential binding of Delta-Coen over Lambda-Coen to the precipitated DNA, but there is no chiral selectivity for Cosep. Competition experiments show that the binding strengths of these three cobalt-amine compounds to aggregated DNA are comparable. A charge neutralization of 88-90% is required for DNA condensation. Our data indicate that 1) electrostatic interaction is the main driving force for binding of multivalent cations to DNA; 2) DNA condensation is dependent on the structure of the condensing agent; and 3) the hydration pattern or polarization of water molecules on the surface of condensing agents plays an important role in DNA condensation and chiral recognition.

Osteoarthritis-susceptibility locus on chromosome 11q, detected by linkage

We present a two-stage genomewide scan for osteoarthritis-susceptibility loci, using 481 families that each contain at least one affected sibling pair. The first stage, with 272 microsatellite markers and 297 families, involved a sparse map covering 23 chromosomes at intervals of approximately 15 cM. Sixteen markers that showed evidence of linkage at nominal P</=.05 were then taken through to the second stage, with an additional 184 families. This second stage confirmed evidence of linkage for markers on chromosome 11q. Additional markers from this region were then typed to create a denser map. We obtained a maximum single-point LOD score, at D11S901, of 2.40 (P=.0004) and a maximum multipoint-LOD score of 3.15, between markers D11S1358 and D11S35. A subset of 196 of the 481 families, comprising affected female sibling pairs, generated a corrected LOD score of 2.54 (P=.0003) for marker D11S901, with evidence for linkage extending 12 cM proximal to this marker. When we stratified for affected male sibling pairs there was no evidence of linkage to chromosome 11. Our data suggest that a female-specific susceptibility gene for idiopathic osteoarthritis is located on chromosome 11q.

1-anilino-8-naphthalene sulfonate as a protein conformational tightening agent

1-Anilino-8-naphthalene sulfonate (ANS) anion is conventionally considered to bind to preexisting hydrophobic (nonpolar) surfaces of proteins, primarily through its nonpolar anilino-naphthalene group. Such binding is followed by an increase in ANS fluorescence intensity, similar to that occurring when ANS is dissolved in organic solvents. It is generally assumed that neither the negative sulfonate charge on the ANS, nor charges on the protein, participate significantly in ANS-protein interaction. However, titration calorimetry has demonstrated that most ANS binding to a number of proteins occurs through electrostatic forces, in which ion pairs are formed between ANS sulfonate groups and cationic groups on the proteins (D. Matulis and R. E. Lovrien, Biophys. J., 1998, Vol. 74, pp. 1-8). Here we show by viscometry and diffusion coefficient measurements that bovine serum albumin and gamma-globulin, starting from their acid-expanded, most hydrated conformations, undergo extensive molecular compaction upon ANS binding. As the cationic protein binds negatively charged ANS anion it also takes up positively charged protons from water to compensate the effect of the negative charge, and leaves the free hydroxide anions in solution thus shifting pH upward (the Scatchard-Black effect). These results indicate that ANS is not always a definitive reporter of protein molecular conformation that existed before ANS binding. Instead, ANS reports on a conformationally tightened state produced by the interplay of ionic and hydrophobic characters of both protein and ligand.

Buffer effects on EcoRV kinetics as measured by fluorescent staining and digital imaging of plasmid cleavage

We have developed a protocol to quantify polymer DNA cleavage which replaces the traditional radiolabeling and scintillation counting with fluorescent staining and digital imaging. This procedure offers high sensitivity, speed, and convenience, while avoiding waste and error associated with traditional 32P radiolabeling. This protocol was used to measure cleavage of pBR322 plasmid DNA by EcoRV, a type II restriction enzyme. EcoRV was found to exhibit an order of magnitude difference in binding in two apparently similar buffers used in previous investigations. To determine the origin of this effect, we measured reaction kinetics in buffers of different chemical nature and concentration: Tris, bis-Tris propane, Tes, Hepes, and cacodylate. We found that buffer concentration and identity had significant effects on EcoRV reaction velocity through large changes in specific binding and nonspecific binding (reflected in the Michaelis constant Km and the dissociation constant for nonspecific binding Kns). There were only small changes in Vmax. The source of the buffer effect is the protonated amines common to many pH buffers. These buffer cations likely act as counterions screening DNA phosphates, where both the protonated buffer structure and concentration affect enzyme binding strength. It appears that by choosing anionic buffers or zwitterionic buffers with a buried positive charge, buffer influence on the protein binding to DNA can be largely eliminated.

The homo-dimeric form of ADP-ribosyl cyclase in solution

ADP-ribosyl cyclase is a multi-functional enzyme that catalyzes the formation of two Ca2+ signaling molecules, cyclic ADP-ribose (cADPR) and nicotinic acid adenine dinucleotide phosphate (NAADP). X-ray crystallography of three different crystal forms shows that it is a non-covalent dimer. Chemical cross-linking and dynamic light scattering were used in this study to determine if the cyclase is also a non-covalent dimer in solution. Treatment of the cyclase in dilute solution (0.05 mg/ml) with dimethylsuberimidate resulted in complete conversion to a species with molecular weight about twice that of the monomeric cyclase. Prolonged cross-linking of the cyclase at four times higher concentration produced also only the covalently linked dimers and no multimer formation was observed. The cross-linked dimer retained full enzymatic activity and readily catalyzed the formation of cADPR from NAD, NAADP from NADP, cyclic ADP-ribose phosphate from NADP, and cyclic GDP-ribose from nicotinamide guanine dinucleotide. Analysis of the autocorrelation functions obtained from dynamic light scattering measurements indicated the cyclase solution (2 mg/ml) was composed of a single molecular species and its diffusion coefficient was measured to be 7. 4x10-7 cm2/s. Computer modeling using the crystallographic dimensions of the non-covalent cyclase dimer, a donut shaped molecule with a central cavity and overall dimensions of 7x6x3 nm, gave a value for the diffusion coefficient essentially the same as that measured. These results indicate the cyclase is a non-covalent dimer in solution.

DNA bending by small, mobile multivalent cations

We propose a purely electrostatic mechanism by which small, mobile, multivalent cations can induce DNA bending. A multivalent cation binds at the entrance to the B-DNA major groove, between the two phosphate strands, electrostatically repelling sodium counterions from the neighboring phosphates. The unscreened phosphates on both strands are strongly attracted to the groove-bound cation. This leads to groove closure, accompanied by DNA bending toward the cationic ligand. We explicitly treat the dynamic character of the cation-DNA interaction using an adiabatic approximation, noting that DNA bending is much slower than the diffusion of nonspecifically bound, mobile cations. We make semiquantitative estimates of the free energy components of bending-electrostatic (with a sigmoidal distance-dependent dielectric function), elastic, and entropic cation localization-and find that the equilibrium state is bent B-DNA stabilized with a self-localized cation. This is a bending polaron, formation of which should be critically dependent on the strength of electrostatic interaction and the concentration of highly mobile cations available for self-localization. We predict that the resultant bend will be large (approximately 20-40 degrees), smooth (because it is spread over 6 bp), and infrequent. The stability of such a bend can be variable, from transient to highly stable (static) bending, observable with standard curvature-measuring techniques. We further predict that this bending mechanism will have an unusual sequence dependence: sequences with less binding specificity will be more bent, unless the specific binding site is in the major groove.

Insertion of telomere repeat sequence decreases plasmid DNA condensation by cobalt (III) hexaammine

Telomere repeat sequence (TRS) DNA is found at the termini of most eukaryotic chromosomes. The sequences are highly repetitive and G-rich (e.g., [C(1-3)A/TG(1-3)]n for the yeast Saccharomyces cerevisiae) and are packaged into nonnucleosomal protein-DNA structures in vivo. We have used total intensity light scattering and electron microscopy to monitor the effects of yeast TRS inserts on in vitro DNA condensation by cobalt (III) hexaammine. Insertion of 72 bp of TRS into a 3.3-kb plasmid depresses condensation as seen by light scattering and results in a 22% decrease in condensate thickness as measured by electron microscopy. Analysis of toroidal condensate dimensions suggests that the growth stages of condensation are inhibited by the presence of a TRS insert. The depression in total light scattering intensity is greater when the plasmid is linearized with the TRS at an end (39-49%) than when linearized with the TRS in the interior (18-22%). Circular dichroism of a 95-bp fragment containing the TRS insert gives a spectrum that is intermediate between the A-form and B-form, and the anomalous condensation behavior of the TRS suggests a noncanonical DNA structure. We speculate that under conditions in which the plasmid DNA condenses, the telomeric insert assumes a helical geometry that is similar to the A-form and is incompatible with packing into the otherwise B-form lattice of the condensate interior.

Ionic effects on the elasticity of single DNA molecules

We used a force-measuring laser tweezers apparatus to determine the elastic properties of lambda-bacteriophage DNA as a function of ionic strength and in the presence of multivalent cations. The electrostatic contribution to the persistence length P varied as the inverse of the ionic strength in monovalent salt, as predicted by the standard worm-like polyelectrolyte model. However, ionic strength is not always the dominant variable in determining the elastic properties of DNA. Monovalent and multivalent ions have quite different effects even when present at the same ionic strength. Multivalent ions lead to P values as low as 250-300 A, well below the high-salt "fully neutralized" value of 450-500 A characteristic of DNA in monovalent salt. The ions Mg2+ and Co(NH3)63+, in which the charge is centrally concentrated, yield lower P values than the polyamines putrescine2+ and spermidine3+, in which the charge is linearly distributed. The elastic stretch modulus, S, and P display opposite trends with ionic strength, in contradiction to predictions of macroscopic elasticity theory. DNA is well described as a worm-like chain at concentrations of trivalent cations capable of inducing condensation, if condensation is prevented by keeping the molecule stretched. A retractile force appears in the presence of multivalent cations at molecular extensions that allow intramolecular contacts, suggesting condensation in stretched DNA occurs by a "thermal ratchet" mechanism.

DNA melting investigated by differential scanning calorimetry and Raman spectroscopy

Thermal denaturation of the B form of double-stranded DNA has been probed by differential scanning calorimetry (DSC) and Raman spectroscopy of 160 base pair (bp) fragments of calf thymus DNA. The DSC results indicate a median melting temperature Tm = 75.5 degrees C with calorimetric enthalpy change delta Hcal = 6.7 kcal/mol (bp), van't Hoff enthalpy change delta HVH = 50.4 kcal/mol (cooperative unit), and calorimetric entropy change delta Scal = 19.3 cal/deg.mol (bp), at the experimental conditions of 55 mg DNA/ml in 5 mM sodium cacodylate at pH 6.4. The average cooperative melting unit (nmelt) comprises 7.5 bp. The Raman signature of 160 bp DNA is highly sensitive to temperature. Analyses of several conformation-sensitive Raman bands indicate the following ranges for thermodynamic parameters of melting: 43 < delta HVH < 61 kcal/mol (cooperative unit), 75 < Tm < 80 degrees C and 6 < (nmelt) < 9 bp, consistent with the DSC results. The changes observed in specific Raman band frequencies and intensities as a function of temperature reveal that thermal denaturation is accompanied by disruption of Watson-Crick base pairs, unstacking of the bases and disordering of the B form backbone. These three types of structural change are highly correlated throughout the investigated temperature range of 20 to 93 degrees C. Raman bands diagnostic of purine and pyrimidine unstacking, conformational rearrangements in the deoxyribose-phosphate moieties, and changes in environment of phosphate groups have been identified. Among these, bands at 834 cm-1 (due to a localized vibration of the phosphodiester group), 1240 cm-1 (thymine ring) and 1668 cm-1 (carbonyl groups of dT, dG and dC), are shown by comparison with DSC results to be the most reliable quantitative indicators of DNA melting. Conversely, the intensities of Raman marker bands at 786 cm-1 (cytosine ring), 1014 cm-1 (deoxyribose ring) and 1092 cm-1 (phosphate group) are largely invariant to melting and are proposed as appropriate standards for intensity normalizations.

Electrostatic effects on the stability of condensed DNA in the presence of divalent cations

Cylindrical cell model Poisson-Boltzmann (P-B) calculations are used to evaluate the electrostatic contributions to the relative stability of various DNA conformations (A, B, C, Z, and single-stranded (ss) with charge spacings of 3.38 and 4.2 A) as a function of interhelix distance in a concentrated solution of divalent cations. The divalent ion concentration was set at 100 mM, to compare with our earlier reports of spectroscopic and calorimetric experiments, which demonstrate substantial disruption of B-DNA geometry. Monovalent cations neutralize the DNA phosphates in two ways, corresponding to different experimental situations: 1) There is no significant contribution to the ionic strength from the neutralizing cations, corresponding to DNA condensation from dilute solution and to osmotic stress experiments in which DNA segments are brought into close proximity to each other in the presence of a large excess of buffer. 2) The solution is uniformly concentrated in DNA, so that the neutralizing cations add significantly to those in the buffer at close DNA packing. In case 1), conformations with lower charge density (Z and ssDNA) have markedly lower electrostatic free energies than B-DNA as the DNA molecules approach closely, due largely to ionic entropy. If the divalent cations bind preferentially to single-stranded DNA or a distorted form of B-DNA, as is the case with transition metals, the base pairing and stacking free energies that stabilize the double helix against electrostatic denaturation may be overcome. Strong binding to the bases is favored by the high concentration of divalent cations at the DNA surface arising from the large negative surface potential; the surface concentration increases sharply as the interhelical distance decreases. In case 2), the concentration of neutralizing monovalent cations becomes very large and the electrostatic free energy difference between secondary structures becomes small as the interhelical spacing decreases. Such high ionic concentrations will be expected to modify the stability of DNA by changing water activity as well as by screening electrostatic interactions. This may be the root of the decreased thermal stability of DNA in the presence of high concentrations of magnesium ions.

DNA condensation

Recent progress in our understanding of DNA condensation includes the observation of the collapse of single DNA molecules, greater insights into the intermolecular forces driving condensation, the recognition of helix-structure perturbation in condensed DNA, and the increasing recognition of the likely biological consequences of condensation. DNA condensed with cationic liposomes is an efficient agent for the transfection of eukaryotic cells, with considerable potential interest for gene therapy.

Brownian dynamics simulation of probe diffusion in DNA: effects of probe size, charge and DNA concentration

We have used Brownian dynamics simulation to study probe diffusion in solutions of short chain DNA using our previously developed simulation algorithm. We have examined the effect of probe size, charge, and DNA concentration on the probe diffusion coefficient, with the aim of gaining insight into the diffusion of proteins in a concentrated DNA environment. In these simulations, DNA was modeled as a worm-like chain of hydrodynamically equivalent spherical frictional elements while probe particles were modeled as spheres of given charge and hydrodynamic radius. The simulations allowed for both short range Lennard-Jones interactions and long ranged electrostatic interactions between charged particles. For uncharged systems, we find that the effects of probe size and DNA concentration on the probe diffusion coefficient are consistent with excluded volume models and we interpret our results in terms of both empirical scaling laws and the predictions of scaled particle theory. For charged systems, we observe that the effects of probe size and charge are most pronounced for the smallest probes and interpret the results in terms of the probe charge density. For an ionic strength of 0.1 M we find that, below a critical probe surface charge density, the probe diffusion coefficient is largely independent of probe charge and only weakly dependent on the DNA charge. These effects are discussed in terms of the interactions between the probe and the DNA matrix and are interpreted in terms of both the underlying physics of transport in concentrated solutions and the assumptions of the simulation model.

Raman spectroscopy of DNA-metal complexes. II. The thermal denaturation of DNA in the presence of Sr2+, Ba2+, Mg2+, Ca2+, Mn2+, Co2+, Ni2+, and Cd2+

Differential scanning calorimetry, laser Raman spectroscopy, optical densitometry, and pH potentiometry have been used to investigate DNA melting profiles in the presence of the chloride salts of Ba2+, Sr2+, Mg2+, Ca2+, Mn2+, Co2+, Ni2+, and Cd2+. Metal-DNA interactions have been observed for the molar ratio [M2+]/[PO2-] = 0.6 in aqueous solutions containing 5% by weight of 160 bp mononucleosomal calf thymus DNA. All of the alkaline earth metals, plus Mn2+, elevate the melting temperature of DNA (Tm > 75.5 degrees C), whereas the transition metals Co2+, Ni2+, and Cd2+ lower Tm. Calorimetric (delta Hcal) and van't Hoff (delta HVH) enthalpies of melting range from 6.2-8.7 kcal/mol bp and 75.6-188.6 kcal/mol cooperative unit, respectively, and entropies from 17.5 to 24.7 cal/K mol bp. The average number of base pairs in a cooperative melting unit (<nmelt>) varied from 11.3 to 28.1. No dichotomy was observed between alkaline earth and transition DNA-metal complexes for any of the thermodynamic parameters other than their effects on Tm. These results complement Raman difference spectra, which reveal decreases in backbone order, base unstacking, distortion of glycosyl torsion angles, and rupture of hydrogen bonds, which occur after thermal denaturation. Raman difference spectroscopy shows that transition metals interact with the N7 atom of guanine in duplex DNA. A broader range of interaction sites with single-stranded DNA includes ionic phosphates, the N1 and N7 atoms of purines, and the N3 atom of pyrimidines. For alkaline earth metals, very little interaction was observed with duplex DNA, whereas spectra of single-stranded complexes are very similar to those of melted DNA without metal. However, difference spectra reveal some metal-specific perturbations at 1092 cm-1 (nPO2-), 1258 cm-1 (dC, dA), and 1668 cm-1 (nC==O, dNH2 dT, dG, dC). Increased spectral intensity could also be observed near 1335 cm-1 (dA, dG) for CaDNA. Optical densitometry, employed to detect DNA aggregation, reveals increased turbidity during the melting transition for all divalent DNA-metal complexes, except SrDNA and BaDNA. Turbidity was not observed for DNA in the absence of metal. A correlation was made between DNA melting, aggregation, and the ratio of Raman intensities I1335/I1374. At room temperature, DNA-metal interactions result in a pH drop of 1.2-2.2 units for alkaline earths and more than 2.5 units for transition metals. Sr2+, Ba2+, and Mg2+ cause protonated sites on the DNA to become thermally labile. These results lead to a model that describes DNA aggregation and denaturation during heating in the presence of divalent metal cations; 1) The cations initially interact with the DNA at phosphate and/or base sites, resulting in proton displacement. 2) A combination of metal-base interactions and heating disrupts the base pairing within the DNA duplex. This allows divalent metals and protons to bind to additional sites on the DNA bases during the aggregation/melting process. 3) Strands whose bases have swung open upon disruption are linked to neighboring strands by metal ion bridges. 4) Near the midpoint of the melting transition, thermal energy breaks up the aggregate. We have no evidence to indicate whether metal ion cross-bridges or direct base-base interactions rupture first. 5) Finally, all cross-links break, resulting in single-stranded DNA complexed with metal ions.

Aggregation of melted DNA by divalent metal ion-mediated cross-linking

In an accompanying paper we reported the use of differential scanning calorimetry and optical densitometry to characterize the melting and aggregation of 160 bp fragments of calf thymus DNA during heating in the presence of divalent metal cations. Aggregation is observed as thermal denaturation begins and becomes more extensive with increasing temperature until the melting temperature Tm is reached, after which the aggregates dissolve extensively. The order of effectiveness of the metals in inducing aggregation is generally consistent with their ability to induce melting: Cd > Ni > Co > Mn approximately Ca > Mg. Under our experimental conditions (50 mg/ml DNA, 100 mM MCl2, [metal]/[DNA phosphate] approximately 0.6), no measurable aggregates were observed for BaDNA or SrDNA. In this paper we show that the Shibata-Schurr theory of aggregation in the thermal denaturation region provides a good model for our observations. Free energies of cross-linking, induced by the divalent cations, are estimated to be between 34% and 38% of the free energies of base stacking. The ability of a divalent metal cation to induce DNA aggregation can be attributed to its ability to disrupt DNA base pairing and simultaneously to link two different DNA sites.

DNA condensation by cobalt hexaammine (III) in alcohol-water mixtures: dielectric constant and other solvent effects

DNA molecules condense into compact structures in the presence of a critical concentration of multivalent cations. To probe the contribution of electrostatic forces to condensation, we used mixtures of water with methanol (MeOH), ethanol (EtOH), and isopropanol (iPrOH) to vary the dielectric constant epsilon from 80 to 50. The condensation of pUC18 plasmids by hexaammine cobalt (III), Co(NH3)(3+)6, was monitored by total intensity and dynamic light scattering, electron microscopy, and CD. The total scattering intensity increased as epsilon went from 80 to 70, and the decreased as epsilon decreased further. Ultraviolet spectrophotometry confirmed that the loss of intensity at low epsilon was not due to the particles' settling out of solution. The rate as well as the extent of condensation increased as epsilon was lowered from 80 to 70, and also depended on the species of alcohol (MeOH < EtOH < iPrOH). The hydrodynamic radii RH of the particles, however, remained roughly the same at 300-350 A and was independent of the species of alcohol. RH increased below epsilon = 70. The critical concentration of Co(NH3)6(3+) required to induce DNA condensation decreased from 21 microM to about 16 microM as the dielectric constant decreased from 80 to 70, and decreased moderately with the nonpolarity of the alcohol. The fraction of DNA charge neutralized at the onset of DNA condensation was calculated by a modification of Manning's two-variable counterion condensation theory to be 0.90 +/- 0.01, independent of epsilon. By electron microscopy we observed that the condensed particles changed from about 93% toroids at epsilon = 80 to 89% rods at epsilon = 70 and 98% rods at epsilon = 65. At epsilon lower than 65, DNA collapsed into a network of multistranded fibers. The morphology of condensed DNA particles, whether toroids, rods, or fibers, was independent of the alcohol species. CD spectra in ethanol-water mixtures indicated that both closed circular and linearized plasmids were in the B conformation when condensed with Co(NH3)6(3+) at epsilon > or = 70, although the closed circular molecules exhibited a weak psi-DNA spectrum. A transition from the B to A form took place between epsilon = 70 and 60, well above the normal dielectric constant of epsilon = 40 for this transition, indicating that ethanol and Co(NH3)6(3+) synergistically promote the B-A transition. We interpret these results to mean that alcohols have both electrostatic and structural effects on DNA, leading to three regimes of condensation. At the lowest alcohol concentrations the B conformation is stable and condensation is relatively slow, allowing time for the packing adjustments necessary to form toroids.(ABSTRACT TRUNCATED AT 400 WORDS)

Condensation of plasmids enhanced by Z-DNA conformation of d(CG)n inserts

DNA molecules collapse into compact structures in the presence of multivalent cations. To probe the possible importance of supercoiling and conformational effects, pUC18 plasmids (2686 bp) were modified by inserting 12-bp and 20-bp alternating d(CG)n sequences, which are capable of converting to a left-handed Z-conformation under appropriate conditions, into the polycloning region. Condensation was induced by rapid addition of hexaammine cobalt(III) [Co(NH3)6(3+)] and monitored by laser light scattering and electron microscopy. Light scattering shows that plasmids with longer d(CG)n inserts condense more extensively at natural superhelical densities. Electron microscopy indicates that the morphological distribution of condensed d(CG)n-containing plasmids changes as a function of Co(NH3)6(3+) concentration. At lower Co(NH3)6(3+) concentration, the proportion of rods is higher, and at higher Co-(NH3)6(3+) concentration, most of the condensates have the form of toroids. In addition, the inner radii of the toroids are much smaller relative to condensed pUC18 under the same conditions. Enzymatic analysis and chemical probing show that the d(CG)n inserts in naturally supercoiled plasmids have extensively converted from B-form to Z-form in the presence of Co(NH3)6(3+) at the upper range of concentrations under which condensation occurs. To determine whether the enhanced condensation of d(CG)n-containing plasmids results from the change of superhelical density due to the B-Z transition, we treated wild-type pUC18 molecules with topoisomerase I and varying amounts of ethidium bromide to generate a range of supercoil densities. Light scattering indicates that supercoiling did not affect the condensation process.(ABSTRACT TRUNCATED AT 250 WORDS)

Gel electrophoresis measurement of counterion condensation on DNA

We used agarose gel electrophoresis to measure the effective charge neutralization of DNA by counterions of different structure and valence, including Na+, Mg2+, Co(NH3)3+6, and spermidine3+, which competed for binding with an excess of Tris acetate buffer. Linear DNA molecules ranged in size from 1 to 5 kilobases, and supercoiled plasmid pUC18 was also measured. In all cases, the results were in good agreement with theoretical predictions from counterion condensation theory for two-counterion mixtures.

Condensation of supercoiled DNA induced by MnCl2

Multivalent cations condense DNA in vitro, but it had been thought that a valence of at least + 3 was required in aqueous solution. We have found that Mn2+ can produce toroidal condensates of supercoiled plasmid DNA, but not of linearized plasmid. Mg2+ does not cause condensation, and neither MgCl2 nor NaCl can negate the effect of MnCl2, indicating that the condensation mechanism with Mn is not primarily electrostatic. Supercoiled MnDNA is more extensively digested than the linear form by S1 nuclease. Supercoiling appears to cooperate with Mn2+ in stabilizing helix distortions and also provides a "pressure" that enhances lateral association.

Brownian dynamics simulations of probe and self-diffusion in concentrated protein and DNA solutions

We have developed a Brownian dynamics algorithm for simulating probe and self-diffusion in concentrated solutions of DNA and protein. In these simulations, proteins are represented as spheres with radii given by their hydrodynamic radii, while DNA is modeled as a wormlike chain of hydrodynamically equivalent spherical frictional elements. The molecular interaction potentials employed by the program allow for intramolecular stretching and bending motions of the DNA chains, short-range Lennard-Jones interactions, and long-range electrostatic interactions. To test the program, we have carried out simulations of bovine serum albumin (BSA) probe diffusion and DNA self-diffusion in solutions of short-chain DNA as a function of both DNA concentration and solution ionic strength. In addition, we report on simulations of BSA self-diffusion as a function of BSA concentration and ionic strength. Based on a comparison to available experimental data, we find that our simulations accurately predict these transport properties under conditions of physiological salt concentration and predict the stronger concentration dependence observed at lower salt concentrations. These results are discussed in light of the nature of the intermolecular interactions in such systems and the approximations and limitations of the simulation algorithm.

Raman spectroscopy of DNA-metal complexes. I. Interactions and conformational effects of the divalent cations: Mg, Ca, Sr, Ba, Mn, Co, Ni, Cu, Pd, and Cd

Interactions of divalent metal cations (Mg2+, Ca2+, Ba2+, Sr2+, Mn2+, Co2+, Ni2+, Cu2+, Pd2+, and Cd2+) with DNA have been investigated by laser Raman spectroscopy. Both genomic calf-thymus DNA (> 23 kilobase pairs) and mononucleosomal fragments (160 base pairs) were employed as targets of metal interaction in solutions containing 5 weight-% DNA and metal:phosphate molar ratios of 0.6:1. Raman difference spectra reveal that transition metal cations (Mn2+, Co2+, Ni2+, Cu2+, Pd2+, and Cd2+) induce the greatest structural changes in B-DNA. The Raman (vibrational) band differences are extensive and indicate partial disordering of the B-form backbone, reduction in base stacking, reduction in base pairing, and specific metal interaction with acceptor sites on the purine (N7) and pyrimidine (N3) rings. Many of the observed spectral changes parallel those accompanying thermal denaturation of B-DNA and suggest that the metals link the bases of denatured DNA. While exocyclic carbonyls of dT, dG, and dC may stabilize metal ligation, correlation plots show that perturbations of the carbonyls are mainly a consequence of metal-induced denaturation of the double helix. Transition metal interactions with the DNA phosphates are weak in comparison to interactions with the bases, except in the case of Cu2+, which strongly perturbs both base and phosphate group vibrations. On the other hand, the Raman signature of B-DNA is largely unperturbed by Mg2+, Ca2+, Sr2+, and Ba2+, suggesting much weaker interactions of the alkaline earth metals with both base and phosphate sites. A notable exception is a moderate perturbation by alkaline earths of purine N7 sites in 160-base pair DNA, with Ca2+ causing the greatest effect. Correlation plots demonstrate a strong interrelationship between perturbations of Raman bands assigned to ring vibrations of the bases and those of bands assigned to exocyclic carbonyls and backbone phosphodiester groups. However, strong correlations do not occur between the Raman phosphodioxy band (centered near 1092 cm-1) and other Raman bands, suggesting that the former is not highly sensitive to the structural changes induced by divalent metal cations. The structural perturbations induced by divalent cations are much greater for > 23-kilobase pair DNA than for 160-base pair DNA, as evidenced by both the Raman difference spectra and the tendency toward the formation of insoluble aggregates. In the presence of transition metals, aggregation of high-molecular-weight DNA is evident at temperatures as low as 11 degrees C. A relationship between DNA melting and aggregation is proposed in which initial metal binding at major groove sites locally destabilizes the B-DNA double helix, causing displacement of the bases away from one another and exposing additional metal binding sites. Metal cation linkage of two displaced bases would allow separate DNA strands to crosslink. Aggregation is proposed to result from the formation of an extended network of these crosslinks.

Separation of bile vesicles and micelles by gel filtration chromatography: the importance of the intermicellar bile salt concentration

Micelles and vesicles coexist in native bile. Mixed micelles are composed of bile salt, phospholipid, and cholesterol. Micellar bile salt is in equilibrium with the aqueous phase bile salt (intermicellar bile salt), and mixed micelles can be converted to cholesterol-phospholipid vesicles by depletion of bile salt. To determine the amount of cholesterol carried in vesicles and micelles, these two populations must be separated without altering the relative proportion of each. Based on the size difference between micelles and vesicles, gel filtration chromatography has been used to accomplish this separation. We reasoned that to maintain the proportion of micelles and vesicles in bile, the column must be equilibrated and eluted with buffer containing the intermicellar bile salt concentration (IMBC) and species. To test this hypothesis we created a model bile composed exclusively of micelles, a solution containing micelles and vesicles, and a model bile containing all vesicles, as demonstrated by quasielastic light scattering. Gel filtration on Sepharose 4B demonstrated that model vesicles and micelles could be separated on a column eluted with buffer containing bile salt at the IMBC. However, a modest decrease in the buffer bile salt concentration (less than 1 mmol/L) resulted in complete conversion of micelles to vesicles. A comparable increase in the buffer bile salt concentration converted vesicles to micelles. Using only taurocholate in the eluting buffer at the IMBC caused a complete shift of micelles to vesicles, whereas using only taurochenodeoxycholate resulted in conversion of vesicles to micelles. An initial collection of rat bile separated on a column equilibrated with the measured IMBC demonstrated that 94% of the cholesterol was in the micellar fractions.(ABSTRACT TRUNCATED AT 250 WORDS)