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.

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.

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 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 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, 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.

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.

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.

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.