The role of methylation in the intrinsic dynamics of B- and Z-DNA

Methylation of cytosine at the 5-carbon position (5 mC) is observed in both prokaryotes and eukaryotes. In humans, DNA methylation at CpG sites plays an important role in gene regulation and has been implicated in development, gene silencing, and cancer. In addition, the CpG dinucleotide is a known hot spot for pathologic mutations genome-wide. CpG tracts may adopt left-handed Z-DNA conformations, which have also been implicated in gene regulation and genomic instability. Methylation facilitates this B-Z transition but the underlying mechanism remains unclear. Herein, four structural models of the dinucleotide d(GC)(5) repeat sequence in B-, methylated B-, Z-, and methylated Z-DNA forms were constructed and an aggregate 100 nanoseconds of molecular dynamics simulations in explicit solvent under physiological conditions was performed for each model. Both unmethylated and methylated B-DNA were found to be more flexible than Z-DNA. However, methylation significantly destabilized the BII, relative to the BI, state through the Gp5mC steps. In addition, methylation decreased the free energy difference between B- and Z-DNA. Comparisons of α/γ backbone torsional angles showed that torsional states changed marginally upon methylation for B-DNA, and Z-DNA. Methylation-induced conformational changes and lower energy differences may contribute to the transition to Z-DNA by methylated, over unmethylated, B-DNA and may be a contributing factor to biological function.     

Effects of base substituents on the hydration of B- and Z-DNA: correlations to the B- to Z-DNA transition.

We present a study of how substituent groups of naturally occurring and modified nucleotide bases affect the degree of hydration of right-handed B-DNA and left-handed Z-DNA. A comparison of poly(dG-dC) and poly(dG-dm5C) titrations with the lipotropic salts of the Hofmeister series infers that the methyl stabilization of cytosines as Z-DNA is primarily a hydrophobic effect. The hydration free energies of various alternating pyrimidine-purine sequences in the two DNA conformations were calculated as solvent free energies from solvent accessible surfaces. Our analysis focused on the N2 amino group of purine bases that sits in the minor groove of the double helix. Removing this amino group from guanine to form inosine (I) destabilizes Z-DNA, while adding this group to adenines to form 2-aminoadenine (A') stabilizes Z-DNA. These predictions were tested by comparing the salt concentrations required to crystallize hexanucleotide sequences that incorporate d(CG), d(CI), d(TA) and d(TA') base pairs as Z-DNA. Combining the current results with our previous analysis of major groove substituents, we derived a thermodynamic cycle that relates the systematic addition, deletion, or substitution of each base substituent to the B- to Z-DNA transition free energy.     

Structural requirements for the action of steroids as quenchers of albumin fluorescence.

1. Androgens, corticoids, gestagens, estrogens and related steroids are effective quenchers of the intrinsic fluorescence of bovine serum albumin. The quenching effect involves the formation of a steroid albumin complex which formation constant (Kf) and free energy of formation (delta G 0) can be determined by fluorescence titration. The fluorimetrically determined delta G 0 values range from -6.5 to -7.5 kcal/mol. 2. 5 alpha-Androstane and 5 alpha-pregnane are effective quenchers of albumin fluorescence, in accord with the essentially hydrophobic nature of the steroid-albumin interaction. Introduction of hydroxy or oxo groups in 5 alpha-androstane decreases the fluorescence quenching action, but the effect of each group declines when other polar groups are present in the steroid molecule. Similar effects occur with 5 alpha-pregnane except that 20-hydroxy (or oxo) duo-polar derivatives are more effective than the parent hydrocarbon. 3. Comparison of delta G 0 values for steroids differing in a single grouping shows that the steroid-albumin interaction is increased by (a) the benzenoid A-ring; (b) sulfate or carboxylate ions in the vicinity of C-3; (c) the 3-oxo group in place of the 3 alpha-hydroxyl (with 5 beta-pregnane derivatives; not with 5 alpha-androstane derivatives); (d) 17 beta-acetyl or 17 beta-hydroxyethyl residues; (e) acetylated or propionated 17 beta-hydroxy groups; (f) acetylated or methylated hydroxy groups at the C-3 of estrogens; (g) delta 5 and delta 6 double bonds; and (h) the 19 beta-methyl group. The maximal variation of delta G 0 determined by affinity-enhancing groups is -0.8 kcal/mol. Conversely, the steroid-albumin interaction is decreased by introduction of (i) oxygen atoms at C-3, C-6, C-11, C-16, and C-17; (j) 17 alpha-ethynyl and 17 alpha-acetoxyl residues; (k) benzoylated or hexahydro-benzoylated beta-hydroxy groups at C-17; (l) acetylated and benzoylated hydroxy groups at C-3; and delta 1 (conjugated) double bond. Oxo groups at C-3, C-6, C-16 and the 16 alpha, 17 alpha-epoxy group are more effective than the corresponding alpha-hydroxyl in decreasing affinity, while at C-11 and C-17, the alpha-hydroxyl is more effective than the beta-hydroxyl and the oxo group. The effect of substituents is influenced by the whole molecular structure, particularly, by the stereostructure at the A/B juncture, and the presence of an oxo group at C-17. 4. The stereospecific effect of substituents at different positions in the steroid molecule suggests that with non-aromatic, A/B trans (planar) steroids, binding to albumin primarily involves the (alpha) rear surface of the B-, C- and D-ring, and possibly, the 17 beta-side chain. With estrogens and A/B cis (dihedral) steroids, the benzenoid A-ring and electron attracting groups at C-3, respectively, may participate in binding.     

Energetics of Z-DNA formation in poly d(A-T), poly d(G-C), and poly d(A-C) poly d(G-T).

The conformational change for the alternating purine-pyrimidine polydeoxyribonucleotides i.e. poly d(A-T), poly d(G-C), and poly d(A-C) poly d(G-T) from a right-handed conformation at room temperature to the left-handed Z-DNA like double helix at elevated temperatures has been studied by UV spectroscopy, Raman spectroscopy, and by adiabatic differential scanning microcalorimetry (DSC) in the presence of Na+ and Mg2+ or Ni2+ respectively as counterions. The differential UV spectra reveal through a hyperchromic shift at around 280nm and a hypochromic shift at 260nm that a conformational change to the left-handed conformation occurs. The Raman spectra clearly show characteristic changes, a drastic decrease of the band at 680cm-1 and the appearance of a new band at 628cm-1, due to the change of the purine bases to the syn conformation upon inversion of the helix-handedness. The course of the transition as function of temperature can be followed quantitatively by plotting the change in the excess heat capacity vs. temperature. The transition enthalpy delta H for the B- to Z-DNA transition per mole base pairs (mbp) amounts to 2.0 +/- 0.2kcal for poly d(G-C), to 4.0 +/- 0.4kcal for poly d(A-T), and to 3.1 +/- 0.3kcal for poly d(A-C) poly d(G-T). The enthalpy change due to the Z-DNA to coil transitions (per mole base pairs) amounts to 11kcal for poly d(G-C), 10.5kcal for poly d(A-T) and 11.3kcal for poly d(A-C) poly d(G-T).     

Analysis of trimethylpsoralen photoreactivity to Z-DNA provides a general in vivo assay for Z-DNA: analysis of the hypersensitivity of (GT)n B-Z junctions

We have described an exonuclease III/photoreversal procedure to map, with base pair resolution, the bases which have photoreacted with 4,5',8-trimethylpsoralen (Me3-psoralen) forming either monoadducts or interstrand cross-links in DNA (20). This assay allows quantitation of relative rates of Me3-psoralen photobinding to bases in DNA at levels as low as one cross-link per 8,000 base pairs. This assay should be useful for a wide variety of applications of Me3-psoralen photobinding to DNA. Here, we demonstrate the applicability of the Me3-psoralen exo III assay for analysis of the conformation of the Z forming sequences (GT)12ATGT and GAATTC(TG)6TA(TG)6. We have shown previously that Me3-psoralen forms crosslinks in the 5'TA within the (CG)6TA(CG)6 sequence when it exists in the B conformation but not when it exists in the Z conformation (34). More recently we have confirmed this result with the exo III assay and have shown at least a hundred fold increase in Me3-psoralen photoreactivity at the 5'AT sequence within the EcoR I sites (GAATTC) which presumably represent B-Z junctions flanking (CG)6TA(CG)6 (20). Here we demonstrate both the characteristic decrease in psoralen photobinding to 5'TAs within (GT)12ATGT and (TG)6TA(TG)6 and the hyperreactivity of B-Z junctions. These characteristic properties of Me3-psoralen photobinding provide an assay for Z-DNA that is applicable in vivo. The general applicability of this approach for assaying Z-DNA in vivo is discussed.     

Effect of N4-methylcytosine and 5-methylcytosine on the stability of the DNA helix

The thermodynamic parameters (delta H, delta S) of the helix-coil transition of self-complementary oligonucleotides d(CGCGCGCG), d(CG5mCGCGCG), d(CG4mCGCGCG), d(GGACCCGGGTCC), d(GGA5mCCCGGGTCC), and d(GGA4mCCCGGGTCC) were determined. The substitution of 4mC for C was found to decrease the melting temperature of the oligonucleotides. The destabilization effect of the two substitutions is equivalent to the change of A.T for G.C pair. The free energy decrease of the helix-coil transition due to the introduction of two 4mC into an octanucleotide was estimated to be 1,24 kcal/mol.     

The tendency to recreate ancestral CG dinucleotides in the human genome

Background  

The CG dinucleotides are known to be deficient in the human genome, due to a high mutation rate from 5-methylated CG to TG and its complementary pair CA. Meanwhile, many cellular functions rely on these CG dinucleotides, such as gene expression controlled by cytosine methylation status. Thus, CG dinucleotides that provide essential functional substrates should be retained in genomes. How these two conflicting processes regarding the fate of CG dinucleotides - i.e., high mutation rate destroying CG dinucleotides, vs. functional processes that require their preservation remains an unsolved question.     

Stabilization of Z-DNA by demethylation of thymine bases: 1.3-A single-crystal structure of d(m5CGUAm5CG)

Methylation of cytosine bases at the C5 position has been known to stabilize Z-DNA. We had previously predicted from calculations of solvent-accessible surfaces that the methyl group at the same position of thymine has a destabilizing effect on Z-DNA. In the current studies, the sequence d(m5CGUAm5CG) has been crystallized and its structure solved as Z-DNA to 1.3-A resolution. A well-defined octahedral hexaaquomagnesium complex was observed to bridge the O4 oxygens of the adjacent uridine bases at the major groove surface, and four well-structured water molecules were found in the minor groove crevice at the d(UA) dinucleotide. These solvent interactions were not observed in the previously published Z-DNA structure of the analogous d(m5CGTAm5CG) sequence. A comparison of the thymine and uridine structures supports our prediction that demethylation of thymine bases helps to stabilize Z-DNA. A comparison of this d(UA)-containing Z-DNA structure with the analogous d(TA) structure shows that access of the O4 position is hindered by the C5 methyl of thymine due to steric and hydrophobic inhibition. In the absence of the methyl group, a magnesium-water complex binds to and slightly affects the structure of the Z-DNA major groove surface. This perturbation of the solvent structure at the major groove surface is translated into a much larger 1.41-A widening of the minor groove crevice, thereby allowing the specific binding of two water molecules at well-defined sites of each internal d(UA) base pair. Possible mechanisms by which modifications at the major groove surface of Z-DNA can affect the solvent properties of the minor groove crevice are discussed.     

The B-Z conformational transition in folded oligodeoxynucleotides: loop size and stability of Z-hairpins

The capacity to assume a left-handed conformation and the thermodynamics of loop formation in concentrated aqueous NaClO4 have been investigated for the following palindromic sequences: d-(CGCGCGAAAAACGCGCG) (A5), d(CGCGCGTTTTTCGCGCG) (T5), d(CGCGCGTACGCGCG) (TA), and d(CGCGCGATCGCGCG) (AT). The results show that (a) each oligomer assumes a Z conformation upon exposure to increasing NaClO4 concentrations; the salt concentration at the transition midpoint is 1.8 M for both A5 and T5 and 3 and 3.5 M for TA and AT, respectively; (b) in high salt the four oligomers exist, over a wide range of nucleotide concentrations (up to 10(-3) M) and of temperature (greater than 0 degrees C), as unimolecular hairpin structures; (c) hairpins TA and AT exhibit, in buffer A, a lower thermal stability with respect to A5 and T5 (delta T about 16 degrees C), contrary to what is observed at low ionic strength; (d) on hairpin formation, the enthalpic term is about -52 kcal/mol for the two 17-mers and -38 kcal/mol for the two 14-mers, while the change in entropy is found to be around -150 eu for A5 and T5 and -115 eu for TA and AT. This thermodynamic picture suggests that a two-residue loop for TA and AT, found at low ionic strength [see preceding paper (Xodo, L.E., Manzini, G., Quadrifoglio, F., van der Marel, G.A., & van Boom, J.H. (1988) Biochemistry (preceding paper in this issue)], is substituted by a longer one including two additional residues from a missing dC.dG base pairing at the top of the stem.     

In vivo assessment of the Z-DNA-forming potential of d(CA.GT)n and d(CG.GC)n sequences cloned into SV40 minichromosomes

Alternating repeated d(CA.GT)n and d(CG.GC)n sequences constitute a significant proportion of the simple repeating elements found in eukaryotic genomic DNA. These sequences are known to form left-handed Z-DNA in vitro. In this paper, we have addressed the question of the in vivo determination of the Z-DNA-forming potential of such sequences in eukaryotic chromatin. For this purpose, we have investigated the ability of a d(CA.GT)30 sequence and a d(CG.GC)5 sequence to form left-handed Z-DNA when cloned into simian virus 40 (SV40) minichromosomes at two different positions: the TaqI site, which occurs in the intron of the T-antigen gene, and the HpaII site, which is located in the late promoter region within the SV40 control region. Formation of Z-DNA at the inserted repeated sequences was analyzed through the change in DNA linkage associated with the B to Z transition. Our results indicate that regardless of: (1) the site of insertion (either TaqI or HpaII), (2) the precise moment of the viral lytic cycle (from 12 h to 48 h postinfection) and (3) the condition of incorporation of the SV40 recombinants to the host cells (either as minichromosomes or as naked DNA, relaxed or negatively supercoiled), neither the d(CA.GT)30 nor the d(CG.GC)5 sequence are stable in the left-handed Z-DNA conformation in the SV40 minichromosome. The biological relevance of these results is discussed.     

DNA binding properties of key sandramycin analogues: systematic examination of the intercalation chromophore

The examination of a key series of chromophore analogues of sandramycin (1) is detailed employing surface plasmon resonance to establish binding constants within a single high affinity bis-intercalation binding site 5'-d(GCATGC)2, and to establish the preference for sandramycin binding to 5'-d(GCXXGC)2 where XX=AT, TA, GC, and CG. From the latter studies, sandramycin was found to exhibit a preference that follows the order: 5'-d(GCATGC)2 > 5'-d(GCGCGC)2, delta deltaG(o)= 0.4 kcal/mol > 5'-d(GCTAGC)2, delta deltaG(o) = 0.9 kcal/mol> or =5'-d(GCCGGC)2, delta delta G(o) = 1.0 kcal/mol although it binds with high affinity to all four deoxyoligonucleotides. The two highest affinity sequences constitute repeating 5'-PuPy motifs with each intercalation event occurring at a 5'-PyPu step. The most effective sequence constitutes the least stable duplex, contains the sterically most accessible minor groove central to the bis-intercalation site, and the ability to accept two gly-NH/T C2 carbonyl H-bonds identified in prior NMR studies. Similarly, the contribution of the individual structural features of the chromophore were assessed with the high affinity duplex sequence 5'-d(GCATGC)2. In addition to the modest affinity differences, one of the most distinguishing features of the high affinity versus lower affinity bis-intercalation or mono-intercalation directly observable by surface plasmon resonance was the temporal stability of the complexes characterized by the exceptionally slow off-rates.     

Synthesis, structure and thermodynamic properties of 8-methylguanine-containing oligonucleotides: Z-DNA under physiological salt conditions.

Various oligonucleotides containing 8-methylguanine (m8G) have been synthesized and their structures and thermodynamic properties investigated. Introduction Of M8G into DNA sequences markedly stabilizes the Z conformation under low salt conditions. The hexamer d(CGC[M8G]CG)2 exhibits a CD spectrum characteristic of the Z conformation under physiological salt conditions. The NOE-restrained refinement unequivocally demonstrated that d(CGC[m8G]CG)2 adopts a Z structure with all guanines in the syn conformation. The refined NMR structure is very similar to the Z form crystal structure of d(CGCGCG)2, with a root mean square deviation of 0.6 between the two structures. The contribution of m8G to the stabilization of Z-DNA has been estimated from the mid-point NaCl concentrations for the B-Z transition of various m8G-containing oligomers. The presence of m8G in d(CGC[m8G]CG)2 stabilizes the Z conformation by at least deltaG = -0.8 kcal/mol relative to the unmodified hexamer. The Z conformation was further stabilized by increasing the number of m8Gs incorporated and destabilized by incorporating syn-A or syn-T, found respectively in the (A,T)-containing alternating and non-alternating pyrimidine-purine sequences. The results suggest that the chemically less reactive m8G base is a useful agent for studying molecular interactions of Z-DNA or other DNA structures that incorporate syn-G conformation.     

Synthesis and high stability of complementary complexes of N-(2-hydroxyethyl)phenazinium derivatives of oligonucleotides.

Two simple methods for the synthesis of oligonucleotides bearing a N-(2-hydroxyethyl)phenazinium (Phn) residue at the 5'- and/or 3'-terminal phosphate groups are proposed. By forming complexes between a dodecanucleotide d(pApApCpCpTpGpTpTpTpGpGpC), a heptanucleotide d(pCpCpApApApCpA), and Phn derivatives of the latter, it is shown that the introduction of a dye at the end of an oligonucleotide chain strongly stabilizes its complementary complexes. The Tmax and the thermodynamic parameters (delta H, delta S, delta G) of complex formation were determined. According to these data, coupling of a dye with the 5'-terminal phosphate group is the most advantageous: delta G(37 degrees C) is increased by 3.59 +/- 0.04 kcal/mol compared to 2.06 +/- 0.04 kcal/mol for 3'-Phn derivatives. The elongation of the linker, which connects the dye to the oligonucleotide, from a dimethylene up to a heptamethylene usually leads to destabilization of the oligonucleotide complex. The complementary complex formed by the 3',5'-di-Phn derivative of the heptanucleotide was found to be the most stable among all duplexes investigated. Relative to the unmodified complex the increase in free energy was 4.96 +/- 0.04 kcal/mol. The association constant of this modified complex at 37 degrees C is 9.5 x 10(6) M-1, whereas the analogous value for the unmodified complex is only 3 x 10(3) M-1.     

Stochastic distribution of a short region of Z-DNA within a long repeated sequence in negatively supercoiled plasmids

We have analyzed, at nucleotide resolution, the progress of the B-to-Z transition as a function of superhelical density in a 2.2-kilobase plasmid containing the sequence d(C-A)31.d(T-G)31. The transition was monitored by means of reactivity to two chemical probes: diethyl pyrocarbonate, which is sensitive to the presence of Z-DNA, and hydroxylamine, which detects B-Z junctions. At a threshold negative superhelical density between about 0.048 and 0.056, hyper-reactivity to diethyl pyrocarbonate appears throughout the CA/TG repeat and remains as the superhelical density is further increased. However, there is no reactivity characteristic of B-Z junctions until the superhelical density reaches 0.084, when single cytosines at each end of the repeat become hyper-reactive to hydroxylamine. A two-dimensional gel analysis of this system by others (Haniford, D. B., and Pulleyblank, D. E. (1983) Nature 302, 632-634) indicates that only about half of the 62 base pairs of the CA/TG repeat undergo the initial transition at omega = 0.056. Our results indicate that this region of Z-DNA is free to exist anywhere along the CA/TG repeat and is probably constantly in motion. Well defined B-Z junctions are seen only when there is sufficient supercoiling to convert the entire CA/TG sequence to Z-DNA. The implications for possible B-Z transitions in chromosomal domains of different sizes are discussed.     

CG sites and their nucleotide environment in the human gamma-delta-beta-globin gene cluster: distribution, frequency and possible frame for differential methylation

Available restriction endonucleases including CG dinucleotides in their target sequences (most of them being unable to cut the DNA when the cytosine of the CG sequence is methylated) have been used to map cloned DNA covering the human gamma-delta-beta globin gene cluster. Since the human DNA fragments were cloned in Escherichia coli, only the internal cytosine in the sequence CCAT GG could be methylated. Thus, any recognized "CG enzyme" site can be detected since they are unmethylated. Results show that frequencies of "CG enzyme" sites regularly decrease from the gamma-globin region to the beta-globin region, the latter being very poor in "CG enzyme"' sites. The array of enzymes used here detects 4 times more CG sites than the classical MspI/HpaII system. Examination of previously sequenced parts of the gamma-delta-beta globin gene cluster shows that CG dinucleotides correspond to an average frequency of 1 out of 104 nucleotides in the gamma-globin region and 1 out of 138 nucleotides in the beta-globin region. In the gamma-globin region, 1 CG out of 4 or 5 may be detected by the enzymes used; the detected frequency is less than 1 out of 10 CG in the beta-region. Analysis of nucleotide environment around CG dinucleotides shows occurrence of local differences, the main sequences being CGG in the 5' side flanking the gamma genes and ACG in the corresponding area of the beta gene. The results presented introduce some new considerations about analysis of cytosine methylation which has been previously proposed as playing a role in the control of the activity of gamma, delta and beta genes respectively.     

The Z-Z junction: the boundary between two out-of-phase Z-DNA regions

The boundary between two segments of Z-DNA that differ in the phase of their syn-anti alternation about the glycosidic bond is termed a Z-Z junction. Using chemical probes and two-dimensional gel electrophoresis, we examined a Z-Z junction consisting of the sequence d[(CG)8C(CG)8] inserted into a plasmid and used energy minimization techniques to devise a three-dimensional model that is consistent with the available data. We show that both alternating CG segments undergo the B-Z transition together to form a Z-Z junction. The junction is very compact, displaying a distinctive reactivity signature at the two base pairs at the junction. In particular, the 5' cytosine of the CC dinucleotide at the junction is hyperreactive toward hydroxylamine, and the two guanines of the GG dinucleotide on the complementary strand are less reactive toward diethyl pyrocarbonate than are the surrounding Z-DNA guanines. Statistical mechanical treatment of the 2-D gel data yields a delta G for forming the Z-Z junction equal to 3.5 kcal, significantly less than the cost of a B-Z junction and approximately equal to the cost of a base out of alternation (i.e., a Z-DNA pyrimidine in the syn conformation). The computer-generated model shows little distortion of the Z helix outside of the central two base pairs, and the energy of the structure and the steric accessibility of the reactive groups are consistent with the data.     

Stereochemical Effects of Methylphosphonate in B- and Z-DNA Helices: Variation in Hydrophobicity and Effective Widths of Grooves

Abstract

Stereochemical effects of methylphosphonate (MP) in B-DNA and Z-DNA duplexes are studied through molecular mechanics approach. Duplexes of different lengths, tetramers, hexamers, dodecamers are examined to assess the interstrand and intrastrand electrostatic effects due to MPs vis-a-vis phosphates. A variety of models which include duplexes with alternating S-MP and R-MP, alternating phosphate and MP and, duplexes posessing MPs in only one of the strands, are examined by considering both the S- and R-stereoisomers. Majority of the calculations are performed with CG sequences to delineate factors responsible for the stability of B- and Z-DNA as well as B × Z-DNA transition under nonionic conditions. The results show that both B- and Z-DNA duplexes are energetically favoured in the presence of MP due to overwhelming reduction in intrastrand as well as interstrand electrostatic repulsive interactions. The effect is distinct in oligomers longer than tetramers. Comparison of energetics of MP B- and Z-DNA duplexes suggests that an oligodeoxynucleotide such as d(CG)₆ with all phosphates replaced by MPs may favour equally both B- and Z-DNA conformations. The analysis further provides an estimate of electrostatic interactions, operating at the grooves under a variety of conditions. Several specific and localised effects due to S-MP and R-MP are seen at CG and GC steps in various B-DNA and Z-DNA models. S-MP in B- DNA reduces the effective major groove width by nearly 3 Å hence denying access to the functional groups of endonucleases thereby enhancing the resistance of MP-DNA to enzymatic digestion. Further, methyl groups of MP render the surface of the DNA helix to be significantly hydrophobic which may explain higher permeability of MP-DNA in membranes as well as its less soluble nature in aqueous media.     

The thermodynamics of bovine and porcine insulin and proinsulin association determined by concentration difference spectroscopy

Difference spectroscopy was used to determine the equilibrium constants and thermodynamic parameters for the monomer-dimer association of bovine and porcine insulin and bovine proinsulin at pH 2.0 and 7.0. At pH 2 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine insulin were found to be -6.6 kcal/mol, -18 cal/mol-deg, and -12 kcal/mol, respectively. Porcine insulin behaved similarly to bovine insulin in its dimerization properties in that delta G degree 25, delta S degree, and delta H degree were found to be -6.8 kcal/mol, -14 cal/mol-deg, and -11 kcal/mol, respectively. At pH 7 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine insulin were found to be -7.2 kcal/mol, -16 cal/mol/deg, and -12 kcal/mol, respectively. At pH 7.0 delta G degree 25, delta S degree, and delta H degree for dimerization of porcine insulin were -6.7 kcal/mol, -11.6 cal/mol-deg, and -10 kcal/mol, respectively. The similarity in the thermodynamic parameters of both insulin species at the different pH's suggests that there are minimal structural changes at the monomer-monomer contact site over this pH range. The dimerization of both insulin species is under enthalpic control. This may suggest that the formation of the insulin dimer is not driven by hydrophobic bonding but, rather, is driven by the formation between subunits of four hydrogen bonds in an apolar environment. At pH 2 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine proinsulin were found to be -5.3 kcal/mol, -26 cal/mol-deg, and -13 kcal/mol, respectively. At pH 7 delta G degree 25, delta S degree, and delta H degree for dimerization of proinsulin were -5.9 kcal/mol, -4.2 cal/mol-deg, and -7.2 kcal/mol, respectively. Although the presence of the C-peptide on proinsulin does not drastically affect the overall free energy change of dimer formation (as compared to insulin), the other thermodynamic parameters are rather drastically altered. This may be because of electrostatic interactions of groups on the C-peptide with groups on the B-chain which are near the subunit contact site in the insulin dimer.