Insight into the process of product expulsion in cellobiohydrolase Cel6A from Trichoderma reesei by computational modeling

Glycoside hydrolase cellulase family 6 from Trichoderma reesei (TrCel6A) is an important cellobiohydrolase to hydrolyze cellooligosaccharide into cellobiose. The knowledge of enzymatic mechanisms is critical for improving the conversion efficiency of cellulose into ethanol or other chemicals. However, the process of product expulsion, a key component of enzymatic depolymerization, from TrCel6A has not yet been described in detail. Here, conventional molecular dynamics and steered molecular dynamics (SMD) were applied to study product expulsion from TrCel6A. Tyr103 may be a crucial residue in product expulsion given that it exhibits two different posthydrolytic conformations. In one conformation, Tyr103 rotates to open the ?3 subsite. However, Tyr103 does not rotate in the other conformation. Three different routes for product expulsion were proposed on the basis of the two different conformations. The total energy barriers of the three routes were calculated through SMD simulations. The total energy barrier of product expulsion through Route 1, in which Tyr103 does not rotate, was 22.2 kcal·mol^?1. The total energy barriers of product expulsion through Routes 2 and 3, in which Tyr103 rotates to open the ?3 subsite, were 10.3 and 14.4 kcal·mol^?1, respectively. Therefore, Routes 2 and 3 have lower energy barriers than Route 1, and Route 2 is the thermodynamically optimal route for product expulsion. Consequently, the rotation of Tyr103 may be crucial for product release from TrCel6A. Results of this work have potential applications in cellulase engineering.     

Product binding varies dramatically between processive and nonprocessive cellulase enzymes

Cellulases hydrolyze β-1,4 glycosidic linkages in cellulose, which are among the most prevalent and stable bonds in Nature. Cellulases comprise many glycoside hydrolase families and exist as processive or nonprocessive enzymes. Product inhibition negatively impacts cellulase action, but experimental measurements of product-binding constants vary significantly, and there is little consensus on the importance of this phenomenon. To provide molecular level insights into cellulase product inhibition, we examine the impact of product binding on processive and nonprocessive cellulases by calculating the binding free energy of cellobiose to the product sites of catalytic domains of processive and nonprocessive enzymes from glycoside hydrolase families 6 and 7. The results suggest that cellobiose binds to processive cellulases much more strongly than nonprocessive cellulases. We also predict that the presence of a cellodextrin bound in the reactant site of the catalytic domain, which is present during enzymatic catalysis, has no effect on product binding in nonprocessive cellulases, whereas it significantly increases product binding to processive cellulases. This difference in product binding correlates with hydrogen bonding between the substrate-side ligand and the cellobiose product in processive cellulase tunnels and the additional stabilization from the longer tunnel-forming loops. The hydrogen bonds between the substrate- and product-side ligands are disrupted by water in nonprocessive cellulase clefts, and the lack of long tunnel-forming loops results in lower affinity of the product ligand. These findings provide new insights into the large discrepancies reported for binding constants for cellulases and suggest that product inhibition will vary significantly based on the amount of productive binding for processive cellulases on cellulose.     

The role of binding subsite A in reactions catalyzed by hen egg-white lysozyme

The role of binding subsite A, located at the terminal of the six binding subsites of hen egg-white lysozyme, in substrate binding and catalytic reactions was investigated by kinetic studies using a chemical modification method. Computer simulation showed that, although subsite A participates in the binding of the substrate, a decrease in the affinity of subsite A to the sugar residue does not cause a lowering of the rate of substrate consumption but changes the mode of the reaction by changing the distribution of the products formed. The binding free energies of subsites for Asp-101-modified lysozymes were estimated by data-fitting from the experimental time-courses. The contribution of Asp-101 in hen egg-white lysozyme to the substrate binding at subsite A was estimated to correspond to a binding free energy of about -3 kJ/mol, 30% of the total binding free energy of subsite A. Modification of Asp-101 affected not only the binding free energy of subsite A but also that of subsite C.     

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone,androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

Accurate PDZ/Peptide Binding Specificity with Additive and Polarizable Free Energy Simulations

PDZ domains contain 80–100 amino acids and bind short C-terminal sequences of target proteins. Their specificity is essential for cellular signaling pathways. We studied the binding of the Tiam1 PDZ domain to peptides derived from the C-termini of its Syndecan-1 and Caspr4 targets. We used free energy perturbation (FEP) to characterize the binding energetics of one wild-type and 17 mutant complexes by simulating 21 alchemical transformations between pairs of complexes. Thirteen complexes had known experimental affinities. FEP is a powerful tool to understand protein/ligand binding. It depends, however, on the accuracy of molecular dynamics force fields and conformational sampling. Both aspects require continued testing, especially for ionic mutations. For six mutations that did not modify the net charge, we obtained excellent agreement with experiment using the additive, AMBER ff99SB force field, with a root mean square deviation (RMSD) of 0.37 kcal/mol. For six ionic mutations that modified the net charge, agreement was also good, with one large error (3 kcal/mol) and an RMSD of 0.9 kcal/mol for the other five. The large error arose from the overstabilization of a protein/peptide salt bridge by the additive force field. Four of the ionic mutations were also simulated with the polarizable Drude force field, which represents the first test of this force field for protein/ligand binding free energy changes. The large error was eliminated and the RMS error for the four mutations was reduced from 1.8 to 1.2 kcal/mol. The overall accuracy of FEP indicates it can be used to understand PDZ/peptide binding. Importantly, our results show that for ionic mutations in buried regions, electronic polarization plays a significant role.     

Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.     

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.     

Enzymatic hydrolysis of cellulose to glucose lung immobilized β-glucosidase

The production of sugars by enzymatic hydrolysis of cellulose is a multistep process which includes conversion of the intermediate cellobiose to glucose by β-glucosidase. Aside from its role as an intermediate, cellobiose inhibits the endoglucanase components of typical cellulase enzyme systems. Because these enzyme systems often contain insufficient concentrations of β-glucosidase to prevent accumulation of inhibitory cellobiose, this research investigated the use of supplemental immobilized β-glucosidase to increase yield of glucose. Immobilized β-glucosidase from Aspergillus phoenicis was produced by sorption at controlled-pore alumina with about 90% activity retention. The product lost only about 10% of the original activity during an on-stream reaction period of 500 hr with cellobiose as substrate; maximum activity occurred near pH 3.5 and the apparent activation energy was about 11 kcal/mol. The immobilized β-glucosidase was used together with Trichoderma reesei cellulase to hydrolyze cellulosic materials, such as Solka Floc, corn stove and exploded wood. Increased yields of glucose and greater conversions of cellobiose of glucose were observed when the reaction systems contained supplemental immobilized β-glucosidase.     

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.     

Multiple functions of aromatic-carbohydrate interactions in a processive cellulase examined with molecular simulation

Proteins employ aromatic residues for carbohydrate binding in a wide range of biological functions. Glycoside hydrolases, which are ubiquitous in nature, typically exhibit tunnels, clefts, or pockets lined with aromatic residues for processing carbohydrates. Mutation of these aromatic residues often results in significant activity differences on insoluble and soluble substrates. However, the thermodynamic basis and molecular level role of these aromatic residues remain unknown. Here, we calculate the relative ligand binding free energy by mutating tryptophans in the Trichoderma reesei family 6 cellulase (Cel6A) to alanine. Removal of aromatic residues near the catalytic site has little impact on the ligand binding free energy, suggesting that aromatic residues immediately upstream of the active site are not directly involved in binding, but play a role in the glucopyranose ring distortion necessary for catalysis. Removal of aromatic residues at the entrance and exit of the Cel6A tunnel, however, dramatically impacts the binding affinity, suggesting that these residues play a role in chain acquisition and product stabilization, respectively. The roles suggested from differences in binding affinity are confirmed by molecular dynamics and normal mode analysis. Surprisingly, our results illustrate that aromatic-carbohydrate interactions vary dramatically depending on the position in the enzyme tunnel. As aromatic-carbohydrate interactions are present in all carbohydrate-active enzymes, these results have implications for understanding protein structure-function relationships in carbohydrate metabolism and recognition, carbon turnover in nature, and protein engineering strategies for biomass utilization. Generally, these results suggest that nature employs aromatic-carbohydrate interactions with a wide range of binding affinities for diverse functions.     

Structures of Phanerochaete chrysosporium Cel7D in complex with product and inhibitors

The cellobiohydrolase Pc_Cel7D is the major cellulase produced by the white-rot fungus Phanerochaete chrysosporium, constituting approximately 10% of the total secreted protein in liquid culture on cellulose. The enzyme is classified into family 7 of the glycoside hydrolases and, like other family members, catalyses cellulose hydrolysis with net retention of the anomeric carbon configuration. Previous work described the apo structure of the enzyme. Here we investigate the binding of the product, cellobiose, and several inhibitors, i.e. lactose, cellobioimidazole, Tris/HCl, calcium and a thio-linked substrate analogue, methyl 4-S-beta-cellobiosyl-4-thio-beta-cellobioside (GG-S-GG). The three disaccharides bind in the glucosyl-binding subsites +1 and +2, close to the exit of the cellulose-binding tunnel/cleft. Pc_Cel7D binds to lactose more strongly than cellobiose, while the opposite is true for the homologous Trichoderma reesei cellobiohydrolase Tr_Cel7A. Although both sugars bind Pc_Cel7D in a similar fashion, the different preferences can be explained by varying interactions with nearby loops. Cellobioimidazole is bound at a slightly different position, displaced approximately 2 A toward the catalytic centre. Thus the Pc_Cel7D complexes provide evidence for two binding modes of the reducing-end cellobiosyl moiety; this conclusion is confirmed by comparison with other available structures. The combined results suggest that hydrolysis of the glycosyl-enzyme intermediate may not require the prior release of the cellobiose product from the enzyme. Further, the structure obtained in the presence of both GG-S-GG and cellobiose revealed electron density for Tris at the catalytic centre. Inhibition experiments confirm that both Tris and calcium are effective inhibitors at the conditions used for crystallization.     

Comparison of binding of mixed ribose-deoxyribose analogues of CUCU to a ribozyme and to GGAGAA by equilibrium dialysis: evidence for ribozyme specific interactions with 2' OH groups

Dissociation constants at 15 degrees C were measured by equilibrium dialysis for the binding of rCrUrCrU, dCrUrCrU, rCdUrCrU, rCrUdCrU, and rCrUrCdU to the L-21 ScaI form of the self-splicing group I LSU intron from Tetrahymena thermophila. Substitution of deoxyribose for ribose in each of the middle two positions makes the free energy change for binding 1-2 kcal/mol less favorable, compared to about 0.3 kcal/mol less favorable for each of the terminal positions. Dissociation constants for binding of the same oligomers to rGGAGAA were measured by optical melting methods. Substitution of a single deoxyribose for ribose makes the free energy change for binding less favorable by 0.4-0.9 kcal/mol for this simple duplex formation. Comparison of the effects for binding to ribozyme and to rGGAGAA indicate that ribozyme-specific tertiary interactions dependent on the middle two 2' OH groups of rCrUrCrU add about 2 kcal/mol of favorable free energy for binding to L-21 ScaI. Comparisons are made with results from gel retardation studies [Pyle, A. M., McSwiggen, J. A., & Cech, T. R. (1990) Proc. Natl. Acad. Sci. U.S.A. 87, 8187-8191; Pyle, A. M., & Cech, T. R. (1991) Nature (London) 350, 628-631].     

次氯酸与不饱和脂肪酸氧化机制和转化产物的理论研究

摘要 目的：揭示次氯酸与不饱和脂肪酸的氧化反应机制及转化产物。方法：运用Gaussian 16软件包,采用密度泛函方法M06-2X(D3),结合6-31+G(d)基组,在SMD液相水模型水平下进行计算。结果：次氯酸与单不饱和脂肪酸油酸的氧化反应是先形成氯鎓离子中间体,氯鎓离子再与水分子反应生成氯醇,第一步氯鎓离子的形成是控速步骤,其反应活化自由能~8 kcal/mol。环氧化合物和短链的醛是两种转化产物,前者由氯醇脱氯化氢而来,而后者由环氧化合物和氯醇通过系列与次氯酸根的反应而得到,生成它们的控速步骤的反应活化自由能分别为23 和24 kcal/mol。选取两个乙基为取代基的乙烯为油酸模型,其与次氯酸反应的活化自由能仅比油酸高1 kcal/mol。计算得到次氯酸与亚油酸、顺-9,反-11 亚油酸、梓树酸和花生四烯酸模型氧化反应生成氯醇的活化自由能分别是~10、13、16和14 kcal/mol。结论：氯鎓离子中间体机制是次氯酸与不饱和脂肪酸氧化反应的主要机制,反应的活化自由能通常低于15 kcal/mol,意味着此氧化反应动力学上容易发生。氧化产物氯醇能转化为环氧化合物和短链的醛,但活化自由能较高,约23和24 kcal/mol。选取距离双键3个碳以内的结构为不饱和脂肪酸模型,它能够很好地反映不饱和脂肪酸的反应活性。     

Thermodynamics and stoicheiometry of the binding of substrate analogues to uricase.

The subunit composition, metal content, substrate-analogue binding and thermal stability of Aspergillus flavus uricase were determined. A. flavus uricase is a tetramer and contains no copper, iron or any other common prosthetic group. Analytical-gel-filtration and equilibrium-dialysis experiments showed one binding site per subunit for urate analogues. The free energy of xanthine binding was -30.5 kJ (-7.3 kcal)/mol of subunit by equilibrium dialysis and -30.1 kJ (-7.2 kcal)/mol of subunit by microcalorimetry. The enthalpy change for xanthine binding was -15.9 kJ (-3.8 kcal)/mol of subunit when determined from the temperature-dependence of the equilibrium constant and -18.0 kJ (-4.3 kcal)/mol of subunit when measured microcalorimetrically. The thermal inactivation rate of A. flavus uricase increases as protein concentration is decreased. This concentration-dependent instability is not due to subunit dissociation.     

Turkey egg white lysozyme. Free energy, enthalpy, and steady state kinetics of reaction with N-acetylglucosamine oligosaccharides.

Equilibrium and calorimetric studies of substrate binding to turkey egg white (TEW) lysozyme were carried out at 30degrees as a function of pH (2 to 9) and ligand size (monosaccharide to hexasaccharide of N-acetylglucosamine). Steady state kinetic measurements using the N-acetylglucosamine hexasaccharide were carried out as a function of pH (2 to 9) and temperature (20-60degrees). These experiments allow comparison of the properties of TEW lysozyme with those of the hen egg white (HEW) enzyme reported previously (Banerjee, S. K., Holler, E., Hess, G. P., and Rupley, J. A. (1975) J. Biol. Chem. 250, 4355-4367, and references therein). The free energies and enthalpies of oligosaccharide binding are the same for TEW and HEW lysozymes at pH 2 but are less negative for TEW lysozyme at pH 5. The pH dependence of the binding of (GlcNAc)3 and higher oligomers to TEW lysozyme is like that for the binding of beta-methyl-N-acetylglucosaminide to TEW lysozyme. These data indicate that oligosaccharide ligands bind identically with HEW and TEW lysozymes, except for the interactions of residue 101, which is aspartic acid in the HEW protein and glycine in the TEW protein (Larue, J. N., and Speck, J. C., Jr. (1970) J. Biol. Chem. 245, 1985-1991). The pH dependence of kcat is described by apparent pK values of 3.9 and 6.8 and a maximum value of kcat of 0.135 s-1. A value of 21.0 kcal/mol was calculated for deltaH from the temperature dependence of kcat. These values and the dependence of the transglycosylation reaction on acceptor concentration are within experimental error the same as those for HEW lysozyme. The more acid pK seen in the pH rate profile reflects the ionization of Asp-52 in the lysozyme-(GlcNAc)6 complex. The pK of Asp-52 in the free protein is 0.3 pK unit lower. The essential identity of the active sites of the HEW and TEW enzymes, except for the Asp-101 interactions, allows estimation of the thermodynamic properties associated with formation of the two hydrogen bonds between Asp-101 and substrate as deltaG0 = -1.2 kcal/mol, DeltaH0 = -3.6 kcal/mol, and deltaS0 = -7.9 e.u.     

High free energy of lipid/protein interaction in biological membranes

The free energy of lipid/protein interaction in biological membranes is still unknown although extensive partitioning and modelling studies have revealed many partial energetic increments. Multiple site binding kinetics are now applied to four well-studied functional membrane proteins, and mean free energy values (+/-S.D.) of -4.23+/-0.49 kcal/mol for single lipid binding sites and of -89.7+/-35.4 kcal/mol for complete lipid substitution are obtained. These high free energy values point to an important bioenergetic role of lipid/protein interaction in membrane functions.     

Structural basis for enantiomer binding and separation of a common beta-blocker: crystal structure of cellobiohydrolase Cel7A with bound (S)-propranolol at 1.9 A resolution

Cellobiohydrolase Cel7A (previously called CBH 1), the major cellulase produced by the mould fungus Trichoderma reesei, has been successfully exploited as a chiral selector for separation of stereo-isomers of some important pharmaceutical compounds, e.g. adrenergic beta-blockers. Previous investigations, including experiments with catalytically deficient mutants of Cel7A, point unanimously to the active site as being responsible for discrimination of enantiomers.In this work the structural basis for enantioselectivity of basic drugs by Cel7A has been studied by X-ray crystallography. The catalytic domain of Cel7A was co-crystallised with the (S)-enantiomer of a common beta-blocker, propranolol, at pH 7, and the structure of the complex was determined and refined at 1. 9 A resolution. Indeed, (S)-propranolol binds at the active site, in glucosyl-binding subsites -1/+1. The catalytic residues Glu212 and Glu217 make tight salt links with the secondary amino group of (S)-propranolol. The oxygen atom attached to the chiral centre of (S)-propranolol forms hydrogen bonds to the nucleophile Glu212 O(epsilon1) and to Gln175 N(epsilon2), whereas the aromatic naphthyl moiety stacks with the indole ring of Trp376 in site +1. The bidentate charge interaction with the catalytic glutamate residues is apparently crucial, since no enantioselectivity has been obtained with the catalytically deficient mutants E212Q and E217Q.Activity inhibition experiments with wild-type Cel7A were performed in conditions close to those used for crystallisation. Competitive inhibition constants for (R)- and (S)-propranolol were determined at 220 microM and 44 microM, respectively, corresponding to binding free energies of 20 kJ/mol and 24 kJ/mol, respectively. The K(i) value for (R)-propranolol was 57-fold lower than the highest concentration, 12.5 mM, used in co-crystallisation experiments. Still several attempts to obtain a complex with the (R)-enantiomer have failed.By using cellobiose as a selective competing ligand, the retention of the enantiomers of propranolol on the chiral stationary phase (CSP) based on Cel7A mutant D214N were resolved into enantioselective and non- selective binding. The enantioselective binding was weaker for both enantiomers on D214N-CSP than on wild-type-CSP.     

The thermodynamics of hapten and antigen binding by rabbit anti-dinitrophenyl antibody.

Rabbit anti-dinitrophenyl antibody from a serum pool was obtained as five fractions of purified specific antibody by a limiting antigen precipitation method. Each fraction had a different binding affinity for epsilon-N-2,4-dinitrophenyl-L-lysine. The free energy changes for hapten binding to the five antibody fractions varied from -8.35 to -10.0 kcal/mol. An average deltaH of -13.9 kcal/mol was measured for the fractions with a batch calorimeter. The results indicate no significant correlation between enthalpy changes and free energy changes. However, a statistically significant correlation between the free energy changes and the entropy changes was found. The enthalpy of the anti-dinitrophenyl antibody interaction with multivalent dinitrophenyl human serum albumin was determined. These are the first enthalpy measurements of an antibody antigen reaction in which the intrinsic binding enthalpy between the antibody and the determinant group is known. The deltaH for the antigen binding reaction was -10.1 kcal/mol which is 3.8 kcal/mol less exothermic than the deltaH for the hapten binding reaction. The interactions that could lead to such a difference in enthalpy are discussed.