Melting of Cross-Linked DNA V. Cross-Linking Effect Caused by Local Stabilization of the Double Helix

Abstract

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this “cross-linking effect” (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the “cross-linking effect” (18–20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20–40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

Melting of cross-linked DNA v. cross-linking effect caused by local stabilization of the double helix

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this "cross-linking effect" (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the "cross-linking effect" (18 - 20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20- 40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

The flexibility of alternating dA-dT sequences

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

Striking stabilization of Arc repressor by an engineered disulfide bond

A solvent-exposed Cys11-Cys11' disulfide bond was designed to link the antiparallel strands of the beta sheet both in the Arc repressor dimer and in a single-chain variant in which the Arc subunits are connected by a 15-residue peptide tether. In both proteins, the presence of the disulfide bond increased the T(m) by approximately 40 degrees C. In the single-chain background, the disulfide bond stabilized Arc by 8.5 kcal/mol relative to the reduced form, a significantly larger degree of stabilization than caused by other engineered disulfides and most natural disulfides. This exceptional stabilization arises from a modest effective concentration of the Cys11-Cys11' disulfide in the native state (71 M) and an anomalously low effective concentration in the denatured state (40 microM). Disulfide cross-linking of the two beta strands in the single-chain Arc background accelerated refolding by a factor of 170 into the sub-microsecond time scale. However, the major energetic effect of the disulfide occurs after the transition state for Arc refolding, slowing unfolding by 200 000-fold.     

Sequence-dependent mechanics of single DNA molecules

Atomic force microscope-based single-molecule force spectroscopy was employed to measure sequence-dependent mechanical properties of DNA by stretching individual DNA double strands attached between a gold surface and an AFM tip. We discovered that in lambda-phage DNA the previously reported B-S transition, where 'S' represents an overstretched conformation, at 65 pN is followed by a nonequilibrium melting transition at 150 pN. During this transition the DNA is split into single strands that fully recombine upon relaxation. The sequence dependence was investigated in comparative studies with poly(dG-dC) and poly(dA-dT) DNA. Both the B-S and the melting transition occur at significantly lower forces in poly(dA-dT) compared to poly(dG-dC). We made use of the melting transition to prepare single poly(dG-dC) and poly(dA-dT) DNA strands that upon relaxation reannealed into hairpins as a result of their self-complementary sequence. The unzipping of these hairpins directly revealed the base pair-unbinding forces for G-C to be 20 +/- 3 pN and for A-T to be 9 +/- 3 pN.     

A model for the LexA repressor DNA complex

A structural model for the interaction of the LexA repressor DNA binding domain (DBD) with operator DNA is derived by means of Monte Carlo docking. Protein–DNA complexes were generated by docking the LexA repressor DBD NMR solution structure onto both rigid and bent B-DNA structures while giving energy bonuses for contacts in agreement with experimental data. In the resulting complexes, helix III of the LexA repressor DBD is located in the major groove of the DNA and residues Asn-41, Glu-44, and Glu-45 form specific hydrogen bonds with bases of the CTGT DNA sequence. Ser-39, Ala-42, and Asn-41 are involved in a hydrophobic interaction with the methyl group of the first thymine base. Residues in the loop region connecting the two β-sheet strands are involved in nonspecific contacts near the dyad axis of the operator. The contacts observed in the docked complexes cover the entire consensus CTGT half-site DNA operator, thus explaining the specificity of the LexA repressor for such sequences. In addition, a large number of nonspecific interactions between protein and DNA is observed. The agreement between the derived model for the LexA repressor DBD/DNA complex and experimental biochemical results is discussed. © 1995 Wiley-Liss, Inc.     

Mechanical stability of single DNA molecules

Using a modified atomic force microscope (AFM), individual double-stranded (ds) DNA molecules attached to an AFM tip and a gold surface were overstretched, and the mechanical stability of the DNA double helix was investigated. In lambda-phage DNA the previously reported B-S transition at 65 piconewtons (pN) is followed by a second conformational transition, during which the DNA double helix melts into two single strands. Unlike the B-S transition, the melting transition exhibits a pronounced force-loading-rate dependence and a marked hysteresis, characteristic of a nonequilibrium conformational transition. The kinetics of force-induced melting of the double helix, its reannealing kinetics, as well as the influence of ionic strength, temperature, and DNA sequence on the mechanical stability of the double helix were investigated. As expected, the DNA double helix is considerably destabilized under low salt buffer conditions (相似文献   

Binding of polynucleotides to fibroblasts. Effects of complex formation with vinyl analogs of nucleic acids

The nitrocellulose filter assay was used to study the effect of the DNA denaturants glycerol and dimethylsulfoxide (Me2SO) on the lac repressor-operator interaction. Both glycerol and Me2SO decrease the rate of dissociation (kb) of the repressor-operator complex but do not significantly alter the rate of association of repressor and operator. In the presence of 10% Me2SO an almost 10-fold increase of affinity of repressor for operator is observed. A small increase in affinity of repressor for Escherichia coli DNA, chicken blood DNA, and poly(dA-dT) is also found. The results lead to the conclusion that lac repressor when interacting with the operator causes local destabilization of the DNA.     

Structure and Dynamics of Netropsin-Poly(dA-dT)· Poly(dA-dT) Complex: 500 MHz 1H NMR Studies

Abstract

Antibiotic netropsin is known to bind specifically to A and T regions in DNA; the mode of binding being non-intercalative. Obviously, H-bonding between the proton donors of netropsin and acceptors N3 of A and 02 of T comes as a strong possibility which might render this specificity. In netropsin there could be 8 proton donors: four terminal amino groups and four internal imino groups. However, methylation of the terminal amino groups does not alter the binding affinity of netropsin to DNA—but the modification of the internal imino groups significantly lowers the binding affinity. Hence, the logical conclusion is that netropsin may specifically interact with A and T through H-bonding and in order to do so, it should approach the helix from the minor groove. The present paper provides experimental data which verify the conclusion mentioned above.

Using poly(dA-dT)? poly(dA-dT) as a model system it was observed following a thorough theoretical stereochemical analysis that netropsin could bind to -(T-A-T) sequence of the polymer in the B-form through the minor groove by forming specific B-bonding. Models could be either right or left-handed B-DNA with a mono or dinucleotide repeat.

By monitoring the ³¹P signals of free poly(dA-dT) ? poly(dA-dT) and netropsin-poly(dA-dT)? poly(dA-dT) complex we show that the drug changes the DNA structure from essentially a mononucleotide repeat to that of very dominant dinucleotide repeat; however the base- pairing in the DNA-drug complex remain to be Watson-Crick. Whether H-bonding is the specific mode of interaction was judged by monitoring the imino protons of netropsin in the presence of poly(dA-dT) ? poly(dA-dT). This experiment was conducted in 90% H₂O + 10% D₂O Using the time-shared long pulse. It was found that exchangeable imino protons of netropsin appear in the drug-DNA complex and disappear upon increasing the D₂O content; thus confirming that H-bonding is indeed the specific mode of interaction. From these and several NOE measurements, we propose a structure for poly(dA-dT)? poly(dA-dT(-netropsin complex.

In summary, experimental data indicate that netropsin binds to poly(dA-dT)? poly(dA-dT) by forming specific hydrogen bonds and that the binding interaction causes the structure to adopt a Watson-Crick paired dinucleotide repeat motif. The proposed hydrogen bonds can form only if the drug approaches the DNA from the minor groove. Within the NMR time scale the interaction between the ligand and DNA is a fast one. From the NOE experimental data, it appears that poly(dA-dT)? poly(dA-dT) in presence of netropsin exists as an equilibrium mixture of right- and left-handed B-DNA duplexes with a dinucleotide repeat—with a predominance of the left-handed form. The last conclusion is a soft one because it was very difficult to make sure the absence of spin diffusion. In a 400 base pairs long DNA duplex- drug complex (as used in this study), equilibrium between right and left-handed helices can also mean the existence of both helical domains in the same molecule with fast interchange between these domains or/and unhindered motion/propagation of these domains along the helix axis.     

Effect of locked nucleic acid modifications on the thermal stability of noncanonical DNA structure

We studied the kinetic and thermodynamic effects of locked nucleic acid (LNA) modifications on parallel and antiparallel DNA duplexes. The LNA modifications were introduced at cytosine bases of the pyrimidine strand. Kinetic parameters evaluated from melting and annealing curves showed that the association and dissociation rate constants for the formation of the LNA-modified parallel duplex at 25.0 °C were 3 orders of magnitude larger and 6 orders of magnitude smaller, respectively, than that of the unmodified parallel duplex. The activation energy evaluated from the temperature-dependent rate constants was largely altered by the LNA modifications, suggesting that the LNA modifications affected a prenucleation event in the folding process. Moreover, thermodynamic parameters showed that the extent of stabilization by the LNA modification for parallel duplexes (3.6 kcal mol(-1) per one modification) was much more significant than that of antiparallel duplexes (1.6 kcal mol(-1)). This large stabilization was due to the decrease in ΔH° that was more favorable than the decrease in TΔS°. These quantitative parameters demonstrated that LNA modification specifically stabilized the noncanonical parallel duplex. On the basis of these observations, we succeeded to stabilize the parallel duplex by LNA modification at the physiological pH. These results can be useful in the rational design of functional molecules such as more effective antisense and antigene strands, more sensitive strands for detection of target DNA and RNA strands, and molecular switches responding to solution pH.     

Calculation of the stability of mini-duplexes of DNA in an electrolyte solution using a molecular mechanics method in approximating a statistical model of the environment

Approximation of the statistical model of environment (SME) to estimate the energy of the macromolecule in electrolyte solution has been developed and used for calculating the conformational energy of nucleic acids by means of the molecular mechanics method. Calculation of base pairs opening delta Hcalop enthalpies and enthalpies of DNA miniduplexes dissociation delta Hcaldis were performed for 10 types of diduplexes. The approximation SME enable to perform calculations of the absolute base-dependent values of delta Hcalop which coincide in the range of 1 kcal/mol with the experimental base-dependent values of the helix-coil transition enthalpies. Values of dissociation enthalpies delta Hcalop greater than delta Hcaldis for all miniduplexes, the difference of delta Hcalop--delta Hcaldis determine the base-dependent energy of the helix-coil boundary. The values of activation barriers for the strands dissociation delta H d not equal to congruent to 8 kcal/mol and association delta H not equal to as congruent to 4 kcal/mol were obtained for GG/CC and AA/TT duplexes. It is concluded that the approximation SME enables to increase substantially the accuracy of the calculations of the macromolecule conformational rearrangement enthalpies in the electrolytes solution.     

Flanking DNA-sequences contribute to the specific binding of cI-repressor and OR1.

The binding of cI-repressor to a series of mutant operators containing OR1 of the right operator of bacteriophage lambda was investigated. Sites OR2 and/or OR3 were inactivated by either point or deletion mutations. The free energy of binding repressor to OR1 in the wildtype operator, delta G1, is -13.7 +/- 0.3 kcal/mol. delta G1 determined for an OR2- operator created by a single point mutation in OR2 is -13.6 +/- 0.2 kcal/mol. In contrast, delta G1 for the binding of repressor to a cloned synthetic OR1 operator containing only 24 bp of lambda sequence is -12.2 +/- 0.1 kcal/mol. When sequence 5' to OR1 is present, the binding affinity increases to -13.0 +/- 0.1 kcal/mol. In addition, the proximity of OR1 to a fragment-end decreases delta G1 from -13.7 to -12.3 +/- 0.1 kcal/mol. These results suggest that the DNA sequence outside the 17 bp OR1 binding-site contributes to the specific binding of cI-repressor.     

Interactions of Arg2 in the Mnt N-terminal arm with the central and flanking regions of the Mnt operator

Arg2, in the N-terminal arm of the Mnt repressor, plays an important role in determining operator-binding specificity. In the complex of the Mnt tetramer with the 21 base-pair mnt operator, there are four potential sites for Arg2 interactions, two in the central region of the operator and two on the outer flanks of the operator. Single-chain variants of the dimeric N-terminal domain of Mnt containing one Arg2 residue and one Lys2 or Met2 residue were constructed and interactions with operator DNA were probed using Fe. EDTA affinity cleavage. The results of these orientation studies show that the majority of the energetically significant interactions mediated by Arg2 occur in the central region of the mnt operator. The RK2, RA2, and RM2 mutations reduce the free energy of operator binding by 1.7 kcal/mol, 3.3 kcal/mol, and 4.9 kcal/mol, respectively. Double-mutant thermodynamic cycle analyses using the RA2, RM2, and operator variants also reveal interaction free energies between Arg2 and operator base-pairs 9, 10, 11, 12 and 13, which in aggregate account for most of the Arg2 contribution to operator binding.     

Thermodynamics of DNA branching.

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.     

Characterization of 1,25-dihydroxyvitamin D3-receptor complex interactions with DNA by a competitive assay

To identify and assess the specificity of the 1,25-dihydroxyvitamin D₃ chick intestinal cytoplasmic receptor's nucleotide binding site, a competitive DNA-cellulose binding assay was utilized. Unlike other steroid hormone receptors, the 1,25-dihydroxyvitamin D₃-receptor complex binds homologous DNA at 4 °C and does not appear to undergo thermal- or salt-induced activation. Data are presented which suggest that receptor binding discriminates between double-stranded DNA and RNA but is not specific with respect to DNA base sequences. However, DNA base sequence selectivity by 1,25-dihydroxyvitamin D₃-receptor complexes is observed using synthetic polydeoxyribonucleotides, particularly, poly(dA-dT) · poly(dA-dT) and poly(dA) · poly(dT). Preference for double-stranded over single-stranded DNA was also observed. Consistent with this finding, both actinomycin D and ethidium bromide caused a dose-dependent inhibition of receptor binding to DNA-cellulose. It is concluded that the 1,25-dihydroxyvitamin D₃-receptor complex has specificity for AT-rich segments of double-stranded DNA and that this interaction is not merely electrostatic, but also involves hydrophobic interaction with the major and/or minor grooves of the DNA helix.     

Inhibition of DNA polymerase activity by thymine hydrates.

Ultraviolet irradiation of DNA results in various pyrimidine modifications. We have demonstrated formation of both cis-thymine hydrate and trans-thymine hydrate (6-hydroxy-5,6-dihydrothymine) in UV-irradiated poly(dA-dT):poly(dA-dT). Both are released from DNA as free bases by bacterial and human glycosylases. Thymine hydrates are stable in DNA and can be detected in control, unirradiated substrates. We examined the effects of thymine hydrates in UV-irradiated substrate poly(dA-dT):poly(dA-dT) on E. coli DNA polymerase I activity. Enzymic incorporation of labeled thymidine-5'-monophosphate significantly decreased with increasing UV dose. Reversal of DNA thymine hydrates to thymines by mild heating of the substrate prior to enzymic reaction resulted in partial recovery of nucleotide incorporation. Cyclobutane thymine dimers are formed between non-adjacent thymines in UV-irradiated poly(dA-dT):poly(dA-dT). These are responsible for the incomplete recovery of DNA polymerase activity following heating due to their heat stability. Analyses of the irradiated and hydrolyzed substrate also demonstrated formation of minor yields of photoproducts formed by covalent linkage of adjacent thymines and adenines by UV-irradiation. Therefore, the thymine hydrates formed in UV-irradiated DNA partially inhibit polymerase activity during DNA synthesis and thus could be potentially lethal if unrepaired.     

Sequence Specific Binding of Beta Carboline Alkaloid Harmalol with Deoxyribonucleotides: Binding Heterogeneity,Conformational, Thermodynamic and Cytotoxic Aspects

Background

Base dependent binding of the cytotoxic alkaloid harmalol to four synthetic polynucleotides, poly(dA).poly(dT), poly(dA-dT).poly(dA-dT), poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) was examined by various photophysical and calorimetric studies, and molecular docking.

Methodology/Principal Findings

Binding data obtained from absorbance according to neighbor exclusion model indicated that the binding constant decreased in the order poly(dG-dC).poly(dG-dC)>poly(dA-dT).poly(dA-dT)>poly(dA).poly(dT)>poly(dG).poly(dC). The same trend was shown by the competition dialysis, change in fluorescence steady state intensity, stabilization against thermal denaturation, increase in the specific viscosity and perturbations in circular dichroism spectra. Among the polynucleotides, poly(dA).poly(dT) and poly(dG).poly(dC) showed positive cooperativity where as poly(dG-dC).poly(dG-dC) and poly(dA-dT).poly(dA-dT) showed non cooperative binding. Isothermal calorimetric data on the other hand showed enthalpy driven exothermic binding with a hydrophobic contribution to the binding Gibbs energy with poly(dG-dC).poly(dG-dC), and poly(dA-dT).poly(dA-dT) where as harmalol with poly(dA).poly(dT) showed entropy driven endothermic binding and with poly(dG).poly(dC) it was reported to be entropy driven exothermic binding. The study also tested the in vitro chemotherapeutic potential of harmalol in HeLa, MDA-MB-231, A549, and HepG2 cell line by MTT assay.

Conclusions/Significance

Studies unequivocally established that harmalol binds strongly with hetero GC polymer by mechanism of intercalation where the alkaloid resists complete overlap to the DNA base pairs inside the intercalation cavity and showed maximum cytotoxicity on HepG2 with IC₅₀ value of 14 µM. The results contribute to the understanding of binding, specificity, energetic, cytotoxicity and docking of harmalol-DNA complexation that will guide synthetic efforts of medicinal chemists for developing better therapeutic agents.     

Reciprocal cooperative effects of multiple ligand binding to pyruvate kinase

The formation of multiple ligand complexes with muscle pyruvate kinase was measured in terms of dissociation constants and the standard free energies of formation were calculated. The binding of Mn2+ to the enzyme (KA = 55 +/- 5 X 10(-6) M; deltaF degrees = -5.75 +/- 0.05 kcal/mol) and to the enzyme saturated with phosphoenolpyruvate (conditional free energy) KA' = 0.8 +/- 0.4 X 10(-6) M; deltaF degrees = -8.22 +/- 0.34 kcal/mol) has been measured under identical conditions giving a free energy of coupling, delta(deltaF degrees) = -2.47 +/- 0.34 kcal/mol. Such a large negative free energy of coupling is diagnostic of a strong positively cooperative effect in ligand binding. The binding of the substrate phosphoenolpyruvate to free enzyme and the enzyme-Mn2+ complex was, by necessity, measured by different methods. The free energy of phosphoenolpyruvate binding to free enzyme (KS = 1.58 +/- 0.10 X 10(-4)M; deltaF degrees = -5.13 +/- 0.04 kcal/mol) and to the enzyme-Mn2+ complex (K3 = 0.75 +/- 0.10 X 10(-6)M; deltaF degrees = -8.26 +/- 0.07 kcal/mol) also gives a large negative free energy of coupling, delta(deltaF degrees) = -3.16 +/- 0.08 kcal/mol. Such a large negative value confirms reciprocal binding effects between the divalent cation and the substrate phosphoenolpyruvate. The binding of Mn2+ to the enzyme-ADP complex was also investigated and a free energy of coupling, delta(deltaF degrees) = -0.08 +/- 0.08 kcal/mol, was measured, indicative of little or no cooperativity in binding. The free energy of coupling with Mn2+ and pyruvate was measured as -1.52 +/- 0.14 kcal/mol, showing a significant amount of cooperativity in ligand binding but a substantially smaller effect than that observed for phosphoenolpyruvate binding. The magnitude of the coupling free energy may be related to the role of the divalent cation in the formation of the enzyme-substrate complexes. In the absence of the activating monovalent cation, the coupling free energies for phosphoenolpyruvate and pyruvate binding decrease by 40-60% and 25%, respectively, substantiating a role for the monovalent cation in the formation of enzyme-substrate complexes with phosphoenolpyruvate and with pyruvate.     

Energetics of the binding of calcium and troponin I to troponin C from rabbit skeletal muscle.

We determined the free energy of interaction between rabbit skeletal troponin I (TNI) and troponin C (TNC) at 10 degrees and 20 degrees C with fluorescently labeled proteins. The sulfhydryl probe 5-iodoacetamidoeosin (IAE) was attached to cysteine (Cys)-98 of TNC and to Cys-133 of TNI, and each of the labeled proteins was titrated with the other unlabeled protein. The association constant for formation of the complex between labeled TNC (TNC*) and TNI was 6.67 X 10(5) M-1 in 0.3 M KCl, and pH 7.5 at 20 degrees C. In the presence of bound Mg2+, the binding constant increased to 4.58 X 10(7) M-1 and in the presence of excess of Ca2+, the association constant was 5.58 X 10(9) M-1. Very similar association constants were obtained when labeled TNI was titrated with unlabeled TNC. The energetics of Ca2+ binding to TNC* and the complex TNI X TNC* were also determined at 20 degrees C. The two sets of results were used to separately determine the coupling free energy for binding TNI and Mg2+, or Ca2+ to TNC. The results yielded a total coupling free energy of -5.4 kcal. This free energy appeared evenly partitioned into the two species: TNI X TNC(Mg)2 or TNI X TNC(Ca)2, and TNI X TNC(Ca)4. The first two species were each stabilized by -2.6 kcal, with respect to the Ca2+ free TNI X TNC complex, and TNI X TNC(Ca)4 was stabilized by -2.8 kcal, respect to TNI X TNC(Ca)2 or TNI X TNC(Mg)2. The coupling free energy was shown to produce cooperatively complexes formed between TNI and TNC in which the high affinity sites were initially saturated as a function of free Ca2+ to yield TNI X TNC(Ca)4. This saturation occurred in the free Ca2+ concentration range 10(-7) to 10(-5) M. The cooperative strengthening of the linkage between TNI and TNC induced by Ca2+ binding to the Ca2+-specific sites of TNC may have a direct relationship to activation of actomyosin ATPase. The nature of the forces involved in the Ca2+-induced strengthening of the complex is discussed.