A transition state in pieces: major contributions of entropic effects to ligand binding by adenosine deaminase.

Nebularine undergoes hydration at the active site of adenosine deaminase, in a reaction analogous to a partial reaction in the displacement of ammonia from adenosine by water, to generate an inhibitory complex that captures much of the binding affinity expected of an ideal transition-state analogue. Enzyme affinities of several compounds related to nebularine 1,6-hydrate, and to its stable analog 2'-deoxycoformycin, were compared in an effort to identify the structural origins of strong binding. Binding of the stable transition-state analog inhibitor 2'-deoxycoformycin was rendered 9.8 kcal/mol less favorable by removal of substituent ribose, 9.7 kcal/mol less favorable by inversion of the 8-hydroxyl substituent of the diazepine ring, and 10.0 kcal/mol less favorable by removal of atoms 4-6 of the diazepine ring. Binding of the unstable transition-state analog nebularine hydrate was rendered at least 9.9 kcal/mol less favorable by removal of the 6-hydroxyl group and 10.2 kcal/mol less favorable by removal of atoms 1-3 of the pyrimidine ring. In each case, the enzyme exhibited only modest affinity (Kd greater than or equal to 10(-2) M) for the "missing piece", indicating that incorporation of 2 binding determinants within a single molecule permits an additional 7-12 kcal/mol of intrinsic binding energy to be manifested as observed binding energy. These results are consistent with earlier indications that adenosine deaminase may use 10.5 kcal/mol of the intrinsic free energy of binding of the two substrates to place them in positions appropriate for reaction at the active site, overcoming the unfavorable entropy change of -35 eu for the equilibrium of 1,6-hydration of purine ribonucleoside and reducing the equilibrium constant for attainment of the transition state in deamination of adenosine. Thus, adenosine deaminase may achieve up to 8 orders of magnitude of its catalytic power by converting the nonenzymatic, bimolecular, hydration reaction to a monomolecular reaction at its active site. Several new 6-substituted 1,6-dihydropurine ribonucleosides, prepared by photoaddition of formate and by low-temperature addition of organolithium reagents to a derivative of purine ribonucleoside, exhibited Ki values of 9-1400 microM against adenosine deaminase, in accord with the active site's considerable tolerance of bulky leaving groups in substrates. Inhibition by one diastereomer of 6-carboxy-1,6-dihydropurine ribonucleoside was found to be time-dependent, progressing from a weakly bound to a more strongly bound complex.     

Energetic contribution of non-essential 5' sequence to catalysis in a hepatitis delta virus ribozyme

Hepatitis delta virus (HDV) ribozymes employ multiple catalytic strategies to achieve overall rate enhancement of RNA cleavage. These strategies include general acid-base catalysis by a cytosine side chain and involvement of divalent metal ions. Here we used a trans-acting form of the antigenomic ribozyme to examine the contribution of the 5' sequence in the substrate to HDV ribozyme catalysis. The cleavage rate constants increased for substrates with 5' sequence alterations that reduced ground-state binding to the ribozyme. Quantitatively, a plot of activation free energy of chemical conversion versus Gibb's free energy of substrate binding revealed a linear relationship with a slope of -1. This relationship is consistent with a model in which components of the substrate immediately 5' to the cleavage site in the HDV ribozyme-substrate complex destabilize ground-state binding. The intrinsic binding energy derived from the ground-state destabilization could contribute up to 2 kcal/mol toward the total 8.5 kcal/mol reduction in activation free energy for RNA cleavage catalyzed by the HDV ribozyme.     

Binding energetics of phosphorus-containing inhibitors of thermolysin

The importance of a specific hydrogen bond between thermolysin and a phosphonamidate inhibitor, Z-NHCH2-PO(O-)-Leu-Leu (1) [Bartlett, P. A., & Marlowe, C. K. (1987) Science (Washington D.C.) 235, 569-571], has been reevaluated. We have determined the inhibition constants (binding free energies) for thermolysin of phosphonamidate n-hexyl-P(O)(O-)-Leu-Trp-NHMe (4), phosphonate n-hexyl-P-(O)(O-)OCH(iBu)CO-Trp-NHMe (5), and phosphinates n-hexyl-P(O)(O-)CH2CH(iBu)CO-Trp-NHMe (6) and Z-NHCH2PO(O-)CH2CH(iBu)CO-Leu (3). Replacement of the P-NH group by P-CH2 (1----3 and 4----6) weakens the overall binding free energy by about 1.5 kcal/mol. A negligible difference in solvation energy has been measured for these pairs, and the basicity of the P-O- ligand for zinc in each pair remains nearly unchanged as determined by pH titration of their 31P NMR resonances. Therefore, this value of 1.5 kcal/mol can be assigned to the specific hydrogen bond known to exist between the P-NH of 1 and thermolysin [Tronrud, D. E., Holden, H. M., & Matthews, B. W. (1987) Science (Washington, D.C.) 235, 871-574] and inferred to exist between 4 and the enzyme. Substitution of P-O for P-NH (1----2 [Bartlett, P. A., & Marlowe, C. K. (1987) Science (Washington, D.C.) 235, 569-571] and 4----5) weakens the overall binding free energy by 4.1 kcal/mol for each pair as the basicity of the P-O- ligand decreases by about 1.6 pH units. The measured solvation energy difference between 4 and 5 (and by inference between 1 and 2) is negligible.(ABSTRACT TRUNCATED AT 250 WORDS)     

Study of the insulin dimerization: binding free energy calculations and per-residue free energy decomposition

A calculation of the binding free energy for the dimerization of insulin has been performed using the molecular mechanics-generalized Born surface area approach. The calculated absolute binding free energy is -11.9 kcal/mol, in approximate agreement with the experimental value of -7.2 kcal/mol. The results show that the dimerization is mainly due to nonpolar interactions. The role of the hydrogen bonds between the 2 monomers appears to give the direction of the interactions. A per-atom decomposition of the binding free energy has been performed to identify the residues contributing most to the self association free energy. Residues B24-B26 are found to make the largest favorable contributions to the dimerization. Other residues situated at the interface between the 2 monomers were found to make favorable but smaller contributions to the dimerization: Tyr B16, Val B12, and Pro B28, and to an even lesser extent, Gly B23. The energy decomposition on a per-residue basis is in agreement with experimental alanine scanning data. The results obtained from a single trajectory (i.e., the dimer trajectory is also used for the monomer analysis) and 2 trajectories (i.e., separate trajectories are used for the monomer and dimer) are similar.     

Thermodynamic analysis of inducer binding to the lactose repressor protein: contributions of galactosyl hydroxyl groups and beta substituents

Kinetic and equilibrium studies of the binding of modified beta-D-galactoside sugars to the lac repressor were carried out to generate thermodynamic data for protein-inducer interactions. The energetic contributions of the galactosyl hydroxyl groups to binding were assessed by using a series of methyl deoxyfluoro-beta-D-galactosides. The C-3 and C-6 hydroxyls contributed less than or equal to -2.3 and -1.7 +/- 0.3 kcal/mol to the binding free energy change, respectively, whereas the C-4 hydroxyl provided only a nominal contribution (-0.1 +/- 0.2 kcal/mol). Favorable contributions to the total binding free energy change were observed for replacement of O-methyl by S-methyl at the beta-anomeric position and for S-methyl by S-isopropyl. Negative delta H degrees values characteristic of protein-sugar complexes [Quiocho, F. A. (1986) Annu. Rev. Biochem. 55, 287-315] were observed for a series of beta-D-galactosides differing at the beta-glycosidic position. A decrease in delta H degrees of approximately 6 kcal/mol upon replacement of the O-methyl substituent by S-methyl indicates a substantial increase in van der Waals' interactions and/or hydrogen bonding in this region of the ligand binding site. The more favorable free energy change for the binding of the S-isopropyl vs S-methyl compound is due mainly to more positive entropic contributions, consistent with an increase in apolar interactions.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

8 Using pseudorotation as a reaction coordinate in free energy simulations of nucleic acids

Backbone sugar groups are central components of nucleic acids. The conformations of the ribose/deoxyribose can be elegantly described using the concept of pseudorotation (Altona and Sundaralingam, 1972), and are dominated by the C2′- and C3′-endo conformers. The free energy barrier of the transition between these two major puckering modes can be probed by NMR relaxation experiments (Johnson and Hoogstraten, 2008), but an atomic picture of the transition path per se is only available for several truncated nucleoside analogues (Brameld & Goddard III, 1999). Here, we implemented a new free energy simulation method for Molecular Dynamics simulations using pseudorotation as the reaction coordinate (Cremer and Pople, 1975). This allowed us to compute the free energy landscape of a complete pseudorotation cycle. The free energy landscape revealed not only the relative stability of C2′- and C3′-endo conformers, but also the main transition path and its free energy barrier. As a validation of our new approach, we calculated free energy surface of the pseudorotation of guanosine monophosphate. The free energy surface revealed that the C2′-endo conformation is ?1?kcal/mol that is more stable and the free energy barrier for the transition is 4.5–5?kcal/mol. These are in excellent agreement with previous NMR measurements (Zhang et al., 2012; Röder et al., 1975). We have further applied this method to other systems that are important in pre-biotic chemistry, including an RNA duplex with unique 2′, 5′-phosphodiester linkages.     

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.     

Cooperativity of papain-substrate interaction energies in the S2 to S2' subsites

Enzyme-substrate contacts in the hydrolysis of ester substrates by the cysteine protease papain were investigated by systematically altering backbone hydrogen-bonding and side-chain hydrophobic contacts in the substrate and determining each substrate's kinetic constants. The observed specificity energies [defined as delta delta G obs = -RT ln [(kcat/KM)first/(kcat/KM)second)]] of the substrate backbone hydrogen bonds were -2.7 kcal/mol for the P2 NH and -2.6 kcal/mol for the P1 NH when compared against substrates containing esters at those sites. The observed binding energies were -4.0 kcal/mol for the P2 Phe side chain, -1.0 kcal/mol for the P1' C=O, and -2.3 kcal/mol for the P2' NH. The latter three values probably all significantly underestimate the incremental binding energies. The P2 NH, P2 Phe side-chain, and P1 NH contacts display a strong interdependence, or cooperativity, of interaction energies that is characteristic of enzyme-substrate interactions. This interdependence arises largely from the entropic cost of forming the enzyme-substrate transition state. As favorable contacts are added successively to a substrate, the entropic penalty associated with each decreases and the free energy expressed approaches the incremental interaction energy. This is the first report of a graded cooperative effect. Elucidation of favorable enzyme-substrate contacts remote from the catalytic site will assist in the design of highly specific cysteine protease inhibitors.     

Site-bound water and the shortcomings of a less than perfect transition state analogue

Kinetic measurements have shown that substantial enthalpy changes accompany substrate binding by cytidine deaminase, increasing markedly as the reaction proceeds from the ground state (1/K(m), DeltaH = -13 kcal/mol) to the transition state (1/K(tx), DeltaH = -20 kcal/mol) [Snider, M. J., et al. (2000) Biochemistry 39, 9746-9753]. In the present work, we determined the thermodynamic changes associated with the equilibrium binding of inhibitors by cytidine deaminase by isothermal titration calorimetry and van't Hoff analysis of the temperature dependence of their inhibition constants. The results indicate that the binding of the transition state analogue 3,4-dihydrouridine DeltaH = -21 kcal/mol), like that of the transition state itself (DeltaH = -20 kcal/mol), is associated with a large favorable change in enthalpy. The significantly smaller enthalpy change that accompanies the binding of 3,4-dihydrozebularine (DeltaH = -10 kcal/mol), an analogue of 3,4-dihydrouridine in which a hydrogen atom replaces this inhibitor's 4-OH group, is consistent with the view that polar interactions with the substrate at the site of its chemical transformation play a critical role in reducing the enthalpy of activation for substrate hydrolysis. The entropic shortcomings of 3,4-dihydrouridine, in capturing all of the free energy involved in binding the actual transition state, may arise from its inability to displace a water molecule that occupies the binding site normally occupied by product ammonia.     

Thermodynamics of ion-induced RNA folding in the hammerhead ribozyme: an isothermal titration calorimetric study

The hammerhead ribozyme undergoes a well-defined two-stage conformational folding process, induced by the binding of magnesium ions. In this study, we have used isothermal titration calorimetry to analyze the thermodynamics of magnesium binding and magnesium ion-induced folding of the ribozyme. Binding to the natural sequence ribozyme is strongly exothermic and can be analyzed in terms of sequential interaction at two sites with association constants K(A) = 480 and 2840 M(-1). Sequence variants of the hammerhead RNA give very different isothermal titration curves. An A14G variant that cannot undergo ion-induced folding exhibits endothermic binding. By contrast, a deoxyribose G5 variant that can undergo only the first of the two folding transitions gives a complex titration curve. However, despite these differences the ITC data for all three species can be analyzed in terms of the sequential binding of magnesium ions at two sites. While the binding affinities are all in the region of 10(3) M(-1), corresponding to free energies of Delta G degrees = -3.5 to -4 kcal mol(-1), the enthalpic and entropic contributions show much greater variation. The ITC experiments are in good agreement with earlier conformational studies of the folding of the ion-induced folding of the hammerhead ribozyme.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.     

Novel mechanism-based substrates of dihydrofolate reductase and the thermodynamics of ligand binding: A comparison of theory and experiment for 8-methylpterin and 6,8-dimethylpterin

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of the designed mechanism-based pterins, 8-methylpterin and 6,8-dimethylpterin, to dihydrofolate reductase (DHFR), with co-factor nicotinamide adenine dinucleotide phosphate (NADPH). The calculated free energy differences suggest that DHFR.NADPH.6,8-dimethylpterin is thermodynamically more stable than DHFR.NADPH.8-methylpterin by 2.4 kcal/mol when the substrates are protonated and by 1.3 kcal/mol when neutral. The greater binding strength of 6,8-dimethylpterin may be attributed largely to hydration effects. In terms of an appropriate model for the pH-dependent kinetic mechanism, these differences can be interpreted consistently with experimental data obtained from previous kinetic studies, i.e., 6,8-dimethylpterin is a more efficient substrate of vertebrate DHFRs than 8-methylpterin. The kinetic data suggest a value of 6.6 ± 0.2 for the pK_a of the active site Glu-30 in DHFR.NADPH. We have also used experimental data to estimate absolute values for thermodynamic dissociation constants of the active (i.e., protonated) forms of the substrates: these are of the same order as for the binding of folate (0.1–10 μM). The relative binding free energy calculated from the empirically derived dissociation constants for the protonated forms of 8-methylpterin and 6,8-dimethylpterin is 1.4 kcal/mol, a value which compares reasonably well with the theoretical value of 2.4 kcal/mol. © 1993 Wiley-Liss, Inc.     

N-acetylcolchinol O-methyl ether and thiocolchicine, potent analogs of colchicine modified in the C ring. Evaluation of the mechanistic basis for their enhanced biological properties

Two colchicine analogs with modifications only in the C ring are better inhibitors than colchicine of cell growth and tubulin polymerization. Radiolabeled thiocolchicine (with a thiomethyl instead of a methoxy group at position C-10) and N-acetylcolchinol O-methyl ether (NCME) (with a methoxy-substituted benzenoid instead of the methoxy-substituted tropone C ring) were prepared for comparison with colchicine. Scatchard analysis indicated a single binding site with KD values of 1.0-2.3 microM. Thiocolchicine was bound 2-4 times as rapidly as colchicine, but the activation energies of the reactions were nearly identical (18 kcal/mol for colchicine, 20 kcal/mol for thiocolchicine). NCME bound to tubulin in a biphasic reaction. The faster phase was 60 times as fast as colchicine binding at 37 degrees C, and a substantial reaction occurred at 0 degrees C. The rate of the faster phase of NCME binding changed relatively little as a function of temperature, so the activation energy was only 7.0 kcal/mol. Dissociation reactions were also evaluated, and at 37 degrees C the half-lives of the tubulin-drug complexes were 11 min for NCME, 24 h for thiocolchicine, and 27 h for colchicine. Relative dissociation rates as a function of temperature varied little among the drug complexes. Activation energies for the dissociation reactions were 30 kcal/mol for thiocolchicine, 27 kcal/mol for NCME, and 24 kcal/mol for colchicine. Comparison of the activation energies of association and dissociation yielded free energies for the binding reactions of -20 kcal/mol for NCME, -10 kcal/mol for thiocolchicine, and -6 kcal/mol for colchicine. The greater effectiveness of NCME and thiocolchicine as compared with colchicine in biological assays probably derives from their more rapid binding to tubulin and the lower free energies of their binding reactions.     

Predicting the binding properties of cibacron blue F3GA in affinity separation systems

The binding properties of cibacron blue F3GA (CB-F3GA) bound to a model NAD(P)H/FAD(H2)-dependent protein system, namely cytosolic quinone reductase (QR), was characterized by AMBER in an attempt to address the binding properties of immobilized CB-F3GA used in the separation of serum albumin. A favorable binding free energy of -4.52kcal/mol (KD=5.09 x 10(-4)kcal/mol) was determined for CB-F3GA binding by MM-PBSA method, which was found to be a ballpark estimate of empirical values reported in literature (DeltaG approximately -6kcal/mol). We propose that CB-F3GA primarily follows a class III binding motif in presence of FAD in the binding site of QR in solution, while a class II binding motif is observed in the crystal form. It was found that favorable van der Waals/hydrophobic interactions take place in the binding site making a major contribution to a favorably dominating enthalpy of binding (DeltaHtot=-25.87kcal/mol) as compared to a disfavorable binding entropy term (TDeltaStot=-21.35kcal/mol). Additional MM-PBSA experiments in the absence of FAD gave rise to a disfavorable binding free energy for CB in complex with QR, suggesting that FAD is an essential determinant of CB-F3GA binding. This is in contrast to an earlier observation of Denizli et al. on separation of human serum albumin (HSA) by immobilized CB-F3GA in the absence of FAD. Therefore, a class I binding model for CB-F3GA is proposed here to account for the efficient separation of HSA in affinity chromatography systems.     

The standard Gibbs free energy change of hydrolysis of alpha-D-ribose 1-phosphate

The standard Gibbs free energy change of hydrolysis of α-d-ribose 1-phosphate has been measured at pH 7.0, ionic strength 0.1 m, and 25 °C by combining the corresponding values of the two following reactions: adenosine + H₂O ág adenine + ribose (ΔG^0′ = ?2.3 ± 0.1 kcal/mol), catalyzed by adenosine nucleosidase, and ribose 1-phosphate + adenine ág adenosine + P_i (ΔG^0′ = ?3.1 ± 0.1 kcal/mol), catalyzed by adenosine phosphorylase. The standard Gibbs free energy changes were calculated for both reactions from the equilibrium constant. A value of -5.4 ± 0.15 kcal/mol, comparable to that of other hemiacetal phosphoric esters, was obtained for the hydrolysis of ribose 1-phosphate.     

Thermodynamics of glycophorin A transmembrane helix dimerization in C14 betaine micelles

We have used sedimentation equilibrium analytical ultracentrifugation to measure the free energy change for the glycophorin A transmembrane helix-helix dimerization in C14 betaine micelles. By varying the amount of micellar C14 betaine, we show that the protein association reaction in the micellar C14 phase behaves as an ideal-dilute solution. In this hydrophobic environment, the mole-fraction standard state free energy change for self-association of the SNGpA99 glycophorin A construct is -5.7 (+/-0.3, N=5) kcal mol(-1) at 25 degrees C. Compared with previous results carried out in C(8)E(5) micellar solutions, the free energy of dimerization is 1.3 kcal mol(-1) less favorable in C14 betaine micelles. In contrast, when considered on a per-interface basis, the formation of the glycophorin A transmembrane dimer in C14 betaine micelles may be more favorable than the association of several designed transmembrane peptides.     

Temperature dependence of anion transport inhibitor binding to human red cell membranes

The binding characteristics of the inhibitor of anion transport in human red cells, 4,4'-dibenzamido-2,2'-disulfonic stilbene (DBDS), to the anion transport protein of red cell ghost membranes in buffer containing 150 mM NaCl have been measured over the temperature range 0-30 degrees C by equilibrium and stopped-flow fluorescence methods. The equilibrium dissociation constant Keq, increased with temperature. No evidence of a 'break' in the ln(Keq) vs. 1/T plot was found. The standard dissociation enthalpy and entropy changes calculated from the temperature dependence are 9.1 +/- 0.9 kcal/mol and 3.2 +/- 0.3 e.u., respectively. Stopped-flow kinetic studies resolve the overall binding into two steps: a bimolecular association of DBDS with the anion transport protein, followed by a unimolecular rearrangement of the DBDS-protein complex. The rate constants for the individual steps in the binding mechanism can be determined from an analysis of the concentration dependence of the binding time course. Arrhenius plots of the rate constants showed no evidence of a break. Activation energies for the individual steps in the binding mechanism are 11.6 +/- 0.9 kcal/mol (bimolecular, forward step), 17 +/- 2 kcal/mol (bimolecular, reverse step), 6.4 +/- 2.3 kcal/mol (unimolecular, forward step), and 10.6 +/- 1.9 kcal/mol (unimolecular, reverse step). Our results indicate that there is an appreciable enthalpic energy barrier for the bimolecular association of DBDS with the transport protein, and appreciable enthalpic and entropic barriers for the unimolecular rearrangement of the DBDS-protein complex.     

Binding to tubulin of the colchicine analog 2-methoxy-5-(2', 3', 4'-trimethoxyphenyl)tropone. Thermodynamic and kinetic aspects

The thermodynamics and kinetics of the binding to tubulin of the colchicine analog 2-methoxy-5-(2', 3', 4'-trimethoxyphenyl) tropone (termed AC because it lacks the B-ring of colchicine) have been characterized by fluorescence techniques. The fluorescence of AC is weak in aqueous solution and is enhanced 250-fold upon binding to tubulin. The following thermodynamic values were obtained for the interaction at 37 degrees C: K = 3.5 X 10(5) M-1; delta G0 = -7.9 kcal/mol; delta H0 = -6.8 kcal/mol; delta S0 = 3.6 entropy units. The AC-tubulin complex is 1-2 kcal/mol less stable than the colchicine-tubulin complex. The change in fluorescence of AC was employed to measure the kinetics of the association process, and quenching of protein fluorescence was used to measure both association and dissociation. The association process, like that of colchicine, could be resolved into a major fast phase and a minor slow phase. The apparent second order rate constant for the fast phase was found to be 5.2 X 10(4) M-1 S-1 at 37 degrees C, and the activation energy was 13 kcal/mol. This activation energy is 7-11 kcal/mol less than that for the binding of colchicine to tubulin. The difference in activation energies can most easily be rationalized by a mechanism involving a tubulin-induced conformational change in the ligand ( Detrich , H. W., III, Williams, R. C., Jr., Macdonald, T. L., Wilson, L., and Puett , D. (1981) Biochemistry 20, 5999-6005). Such a change would be expected to have a small activation energy in AC because it possesses a freely rotating single bond in place of the B-ring of colchicine.