The enhancement of electrostriction caused by lowering the solvent dielectric constant leads to the decrease of activation energy in trypsin catalysis.

The contribution of electrostriction of the solvent to the stabilization of the negatively charged tetrahedral transition state of a trypsin-catalyzed reaction was probed by means of kinetic studies involving high-pressure and solvent dielectric constant. A good correlation was observed between the increased catalytic efficiency of trypsin and the decreased solvent dielectric constant. When the dielectric constant of the solvents was lowered by 4.68 units, the loss of activation energy and that of free energy of activation were 2.26 kJ/mol and 3.09 kJ/mol, respectively. The activation volume for k(cat) decreased significantly as the dielectric constant of the solvent decreased, indicating that the degree of electrostriction of the solvent around the charged tetrahedral transition state has been enhanced. These observations demonstrate that the increase in the catalytic efficiency of the trypsin reaction with decreasing dielectric constant resulted from the stabilization of electrostatic energy for the formation of an oxyanion hole, and this stabilization was caused by the increase of electrostricted water around the charged tetrahedral transition state. Therefore, we conclude that control of the solvent dielectric constant can stabilize the tetrahedral transition state, and this lowers the activation energy.     

Stabilization of Thermolysin Induced by a Decrease in the PH and Dielectric Constant of the Reaction Medium Resulting from a Temperature Increase in Z-Phe-Phe-OMe Synthesis

Thermolysin was stabilized by increasing the reaction temperature from 90°C to 110°C during peptide synthesis of N-(benzyloxycarbonyl)-_L-phenylalanyl-_L-phenylalanine methyl ester (Z-Phe-Phe-OMe). The stabilization energy was acquired from the drop in both pH and dielectric constant due to the temperature increase. The acquired stabilization energy was as high as ca. 42 kJ/mol (corresponding to 20°C). This acquired stabilization energy did not result from a single event such as a change in electrical charges. It was evaluated as the overall stabilization energy at and around the active site area of the enzyme using an electrostatic potential equation.     

Thermal stabilization of trypsin with glycol chitosan

Glycol chitosan was evaluated as thermoprotectant additive for trypsin in aqueous solutions. Maximal stabilization was achieved by using a polymer/protein ratio of 2 (w/w). The catalytic properties of trypsin were not affected by the presence of the polysaccharide. The enzyme thermostability was increased from 49 °C to 93 °C in the presence of the additive. Trypsin was also 37-fold more stable against incubation at 55 °C and its activation free energy of thermal inactivation was increased by 9.9 kJ/mol when adding glycol chitosan.     

Transition state stabilization of subtilisins in organic media

Electrostatic forces are among the stabilizing interactions that contribute to the high degree of enzyme-transition state complementarity. The active-site polarity, which can differ substaintially from that of water, is thus an important determinant of transition state stabilization. Here we pose the question of whether the rate of an enzymatic reaction proceeding through a charged transition state can be increased by increasing the active-site polarity in an organic solvent. The active-site polarity of subtilisin has been reduced by dehydration and suspension in a nonpolar solvent (tetrahydrofuran), and then increased by adding water to the solvent. Enhancing the local polarity substantially increasing the rate of catalysis, implicating polarity as an important factor in stabilizing the charged tetrahedral transition state. Studies with subtilisins whose active sites have been modified by site-directed mutagenesis support the role of polarity in transition state stabilization. (c) 1994 John Wiley & Sons, Inc.     

Powering DNA repair through substrate electrostatic interactions

The reaction catalyzed by the DNA repair enzyme uracil DNA glycosylase (UDG) proceeds through an unprecedented stepwise mechanism involving a positively charged oxacarbenium ion sugar and uracil anion leaving group. Here we use a novel approach to evaluate the catalytic contribution of electrostatic interactions between four essential phosphodiester groups of the DNA substrate and the cationic transition state. Our strategy was to substitute each of these phosphate groups with an uncharged (R)- or (S)-methylphosphonate linkage (MeP). We then compared the damaging effects of these methylphosphonate substitutions on catalysis with their damaging effects on binding of a cationic 1-azadeoxyribose (1-aza-dR(+)) oxacarbenium ion analogue to the UDG-uracil anion binary complex. A plot of log k(cat)/K(m) for the series of MeP-substituted substrates against log K(D) for binding of the 1-aza-dR(+) inhibitors gives a linear correlation of unit slope, confirming that the electronic features of the transition state resemble that of the 1-aza-dR(+), and that the anionic backbone of DNA is used in transition state stabilization. We estimate that all of the combined phosphodiester interactions with the substrate contribute 6-8 kcal/mol toward lowering the activation barrier, a stabilization that is significant compared to the 16 kcal/mol catalytic power of UDG. However, unlike groups of the enzyme that selectively stabilize the charged transition state by an estimated 7 kcal/mol, these phosphodiester groups also interact strongly in the ground state. To our knowledge, these results provide the first experimental evidence for electrostatic stabilization of a charged enzymatic transition state and intermediate using the anionic backbone of DNA.     

Catalytic mechanism of class B2 metallo-beta-lactamase

The initial nucleophilic substitution step of biapenem hydrolysis catalyzed by a subclass B2 metallo-beta-lactamase (CphA from Aeromonas hydrophila) is investigated using hybrid quantum mechanical/molecular mechanical methods and density functional theory. We focused on a recently proposed catalytic mechanism that involves a non-metal-binding water nucleophile in the active site of the monozinc CphA. Both theoretical models identified a single transition state featuring nearly concomitant nucleophilic addition and elimination steps, and the activation free energy from the potential of mean force calculations was estimated to be approximately 14 kcal/mol. The theoretical results also identified the general base for activating the water nucleophile to be the metal-binding Asp-120 rather than His-118, as suggested earlier. The protonation of Asp-120 leads to cleavage of the O(delta2)-Zn coordination bond, whereas the negatively charged nitrogen leaving group resulting from the ring opening replaces Asp-120 as the fourth ligand of the sole zinc ion. The electrophilic catalysis by the metal ion provides sufficient stabilization for the leaving group to avoid a tetrahedral intermediate. The theoretical studies provided detailed insights into the catalytic strategy of this unique metallo-beta-lactamase.     

Studies on the aminopeptidase activity of rat cathepsin H.

Three synthetic substrates H-Arg-NH-Mec, Bz-Arg-NH-Mec and H-Cit-NH-Mec (Bz, Benzoyl; NH-Mec, 4-methylcoumaryl-7-amide; Cit, citrulline) were used to characterize specificity requirements for the P1-S1 interaction of cathepsin H from rat liver. From rapid equilibrium kinetic studies it was shown that Km, kcat and the specificity constants kcat/Km are quite similar for substrates with a free alpha-amino group. In contrast, a 25-fold decrease of kcat/Km was observed for the N-terminal-blocked substrate Bz-Arg-NH-Mec. The activation energies for H-Arg-NH-Mec and Bz-Arg-NH-Mec were determined to be 37 kJ/mol and 55 kJ/mol, respectively, and the incremental binding energy delta delta Gb of the charged alpha-amino group was estimated to -8.1 kJ/mol at pH 6.8. The shown preference of cathepsin H for the unblocked substrates H-Arg-NH-Mec and H-Cit-NH-Mec was further investigated by inspection of the pH dependence of kcat/Km. The curves of the two substrates with a charged alpha-amino group showed identical bell-shaped profiles which both exhibit pKa1 and pKa2 values of 5.5 and 7.4, respectively, at 30 degrees C. The residue with a pKa1 of 5.5 in the acid limb of the activity profile of H-Arg-NH-Mec was identified by its ionization enthalpy delta Hion = 21 kJ/mol as a beta-carboxylate or gamma-carboxylate of the enzyme, whereas the residue with a pKa2 of 7.4 was assigned to the free alpha-amino group of the substrate with a delta Hion of 59 kJ/mol. Bz-Arg-NH-Mec showed a different pH-activity profile with a pKa1 of 5.4 and a pKa2 of 6.6 at 30 degrees C. Cathepsin H exhibits no preference for a basic P1 side chain as has been shown by the similar kinetics of H-Arg-NH-Mec and the uncharged, isosteric substrate H-Cit-NH-Mec. In summary, specific interactions of an anionic cathepsin H active site residue with the charged alpha-amino group of substrates caused transition state stabilization which proves the enzyme to act preferentially as an aminopeptidase.     

Homotropic cooperative binding of organic solvent vapors by solid trypsin

Homotropic cooperative binding was observed at vapor sorption of organic solvents (acetonitrile, propionitrile, ethanol, 1-propanol, 2-propanol, nitroethane) by dried solid trypsin from porcine pancreas (0.05 g H2O/g protein). The vapor sorption isotherms were obtained by the static method of gas chromatographic headspace analysis at 298 K for 'vapor solvent+solid trypsin' systems in the absence of the liquid phase. All isotherms have a sigmoidal shape with significant sorbate uptake only above the threshold of sorbate thermodynamic activity. On the sorption isotherms of non-hydroxylic sorbates the saturation of trypsin by organic solvent was observed above the sorbate threshold activity. The formation of inclusion compounds with phase transition between solvent-free and solvent-saturated trypsin is supposed. Approximation of obtained isotherms by the Hill equation gives the inclusion stoichiometry S, inclusion free energy, and the Hill constant N of clathrates. The inclusion stoichiometry S depends significantly on the size and shape of sorbate molecules and changes from S=31 mol of sorbate per mol of trypsin for ethanol to S=6 for nitroethane. The inclusion free energies determined for the standard states of pure liquid sorbate and infinitely dilute solution in toluene are in the range from -0.5 to -1.2 kJ/mol and from -3.1 to -8.1 kJ/mol, respectively, per 1 mol of sorbate. The Hill constants are relatively high: from N=5.6 for 1-propanol to N approximately equal to 10(3) for nitroethane. The implication of the obtained results for the interpretation of solvent effects on the enzyme activity and stability in low-water medium is discussed.     

Autooxidation and oxygenation of human hemoglobin]

The processes of reversible oxygen binding and nonreversible autoxidation of human hemoglobin were studied. The activation energy of the oxygen binding, as determined by the temperature dependence of the P50 parameter, was 26 +/- 4 kJ/mol, the activation energy of the autoxidation, as determined by the temperature dependence of the apparent rate constant of autoxidation, was 120 +/- 15 kJ/mol. Pyridoxal phosphate decreased the oxygen affinity of hemoglobin, slightly diminished the cooperativity of the oxygenation process and unaffected the activation energy of the oxygen binding. Pyridoxal phosphate slightly reduced the Bohr coefficient value from 0.70 to 0.65. Pyridoxal phosphate, but not pyridoxal, raised the apparent rate constant of autoxidation reaction. The rate of autoxidation significantly increased as the pH value of the medium decreased, reflecting, probably, protonation of the distal histidine of the hemoglobin. The activation energy of autoxidation was independent of pH. Aliphatic alcohols also increased the rate of the autoxidation process, probably, either by stabilization of the hemoglobin T-state, or by direct nucleophilic displacement of the oxygen molecule.     

Continuum dielectric modelling of the protein-solvent system, and calculation of the long-range electrostatic field of the enzyme phosphoglycerate mutase

The numerical continuum electrostatic method presented previously (Warwicker, J. & Watson, H. C. (1982) J. Mol. Biol., 157, 671-679), is developed with an improved analysis of the protein-solvent system. Inclusion in the model of saturable solvent dielectric, and counterions is discussed and presented. A number of long-range electrostatic field calculations are made on bovine pancreatic trypsin inhibitor to demonstrate the differences between various solvent and counterion models. The long-range potential field, due to polar side-chain and alpha-helix dipole charge, is calculated for the glycolytic enzyme phosphoglycerate mutase. The positive potential in and around the catalytic cleft region is sufficiently large to suggest that it may play a role in long-range attraction of the enzyme's negatively charged substrates. Analogous systems with charge-charge interactions in solvent water are considered. It is suggested that a long-range enzyme-substrate attractive force-field may, in part, offset the repulsive energy arising from overlap of hydration shells between enzyme and substrate.     

Thermodynamic and kinetic processes during the unfolding of BSA in the presence of the mycotoxin patulin

The effects of the mycotoxin patulin on the thermodynamics and kinetics of the transition of bovine serum albumin (BSA) in aqueous solution were studied by Differential Scanning Calorimetry and Photoluminescence methods. Results show that in the presence of patulin, the free enthalpy change during the transition of BSA was decreased by an average of ～ 46 kJ/mol, the free energy change was decreased by ～ 4 kJ/mol, and the activation energy fell from ～ 1546 to ～ 840 kJ/mol. These results indicate that the bioactivity of patulin is based on the kinetic rather than on the thermodynamic properties of the transition. This is the first evidence of the direct interaction of patulin with the free thiol-containing BSA, a process which could contribute to the adverse cyto- and genotoxic effects induced by patulin.     

蛋白质与极性配体复合物的绝对结合自由能计算

介绍了用分子动力学模拟与热力学积分法相结合,模拟蛋白质与配体的绝对结合自由能的方法.通过分子转换法,使蛋白质分子（包括水分子）与配体小分子之间的相互作用逐渐减弱 （或增强）至完全消失（或完全出现）. 运用体约束方法,计算了配体与受体结合后平动、转动自由度的丧失即熵效应所引起的自由能变化.以胰蛋白酶双突变体（D189G/G226D）与极性配体苯甲脒为例,研究了蛋白质活性部位与极性配体的相互作用对结合自由能的影响,该复合物绝对结合自由能的模拟结果(－15.5 kJ/mol)与实验值(－10.5 kJ/mol)相近.     

Irreversible thermal denaturation of elongation factor Ts from Thermus thermophilus effect of the residual structure and intermonomer disulfide bond

The homodimeric wild-type elongation factor Ts, EF-Ts(wt), and its C190A mutant, EF-Ts(C190A), from Thermus thermophilus goes through thermal denaturation in a way consistent with a two state irreversible model with a relatively high activation energy, approximately 530 kJ/mol (Supplemental materials provides a list of 98 activation energies from 54 proteins in various solvent conditions). Removing the intermonomeric disulfide bond by substituting alanine for cysteine 190 affects the rate constant of the irreversible thermal transition. At physiological temperatures, the half-life of the native conformations was estimated to be approximately 21 days for wt and 1.3 days for C190A. Thermally denatured EF-Ts refolds into a molten-globule-like state as indicated by its native-like circular dichroism spectrum in the far UV region and the enhanced fluorescence of the hydrophobic probe, 1-anilinonaphtalene-8-sulphonate. The residual secondary structure observed in the thermally denatured state of EF-Ts at high temperatures affects its apparent temperature of thermal transition, T(trs), independent of the presence or absence of the intermonomeric disulfide bond. The effect of the GdmHCl concentration on the activation energy, E(a), and the temperature, T*, i.e., the temperature at which the rate of the irreversible step is 1 min(-1), indicates that the intermonomeric disulfide bond contributes to the irreversibility of thermal transition of EF-Ts.     

Do enzymes change the nature of transition states? Mapping the transition state for general acid-base catalysis of a serine protease

The properties of the transition state for serine protease-catalyzed hydrolysis of an amide bond were determined for a series of subtilisin variants from Bacillus lentus. There is no significant change in the structure of the enzyme upon introduction of charged mutations S156E/S166D, suggesting that changes in catalytic activity reflect global properties of the enzyme. The effect of charged mutations on the pK(a) of the active site histidine-64 N(epsilon)(2)-H was correlated with changes in the second-order rate constant k(cat)/K(m) for hydrolysis of tetrapeptide anilides at low ionic strength with a Br?nsted slope alpha = 1.1. The solvent isotope effect (D)2(O)(k(cat)/K(m))(1) = 1.4 +/- 0.2. These results are consistent with a rate-limiting breakdown of the tetrahedral intermediate in the acylation step with hydrogen bond stabilization of the departing amine leaving group. There is an increase in the ratio of hydrolysis of succinyl-Ala-Ala-Pro-Phe-anilides for p-nitroaniline versus aniline leaving groups with variants with more basic active site histidines that can be described by the interaction coefficient p(xy) = delta beta(lg)/delta pK(a) (H64) = 0.15. This is attributed to increased hydrogen bonding of the active site imidazolium N-H to the more basic amine leaving group as well as electrostatic destabilization of the transition state. A qualitative characterization of the transition state is presented in terms of a reaction coordinate diagram that is defined by the structure-reactivity parameters.     

Effect of molecular structure and water content on the dielectric relaxation behaviour of amorphous low molecular weight carbohydrates above and below their glass transition

The dielectric relaxation behaviour of several amorphous low molecular weight carbohydrates and their 10% w/w water mixtures has been studied in the supercooled liquid and glassy regions in the frequency range 100 Hz to 100 kHz. The dry carbohydrates show a primary alpha-relaxation (activation energy 250-405 kJ mol(-1)) at temperatures above the calorimetric glass transition temperature, Tg, and, in most cases, a secondary beta-relaxation (activation energy 42-55 kJ mol(-1)) at sub-Tg temperatures. Whilst D-mannose showed a beta-relaxation similar in strength to D-glucose, its deoxy sugar, L-rhamnose showed a relatively weak beta-relaxation. This indicates that the hydroxymethyl group influences relaxation in carbohydrate glasses. Addition of water shifted the alpha-relaxations to lower temperatures and increased the strength of the beta-relaxations. In glucitol this resulted in a merging of the alpha- and beta-relaxations. The beta-relaxation increased in strength and decreased in temperature for the series of water mixtures: D-glucose, maltose, and maltotriose.     

Influence of key residues on the reaction mechanism of the cAMP-dependent protein kinase

The reaction mechanism of the catalytic phosphoryl transfer of cAMP-dependent protein kinase (cAPK) was investigated by semi-empirical AM1 molecular orbital computations of an active site model system derived from the crystal structure of the catalytic subunit of the enzyme. The activation barrier is calculated as 20.7 kcal mol(-1) and the reaction itself to be exothermic by 12.2 kcal mol(-1). The active site residue Asp166, which was often proposed to act as a catalytic base, does not accept a proton in any of the reaction steps. Instead, the hydroxyl hydrogen of serine is shifted to the simultaneously transferred phosphate group of ATP. Although the calculated transition state geometry indicates an associative phosphoryl transfer, no concentration of negative charge is found. To study the influence of protein mutations on the reaction mechanism, we compared two-dimensional energy hypersurfaces of the protein kinase wild-type model and a corresponding mutant in which Asp166 was replaced by alanine. Surprisingly, they show similar energy profiles despite the experimentally known decrease of catalytic activity for corresponding mutants. Furthermore, a model structure was examined, where the charged NH3 group of Lys168 was replaced by a neutral methyl group. The energetic hypersurface of this hypothetical mutant shows two possible pathways for phosphoryl transfer, which both require significantly higher activation energies than the other systems investigated, while the energetic stabilization of the reaction product is similar in all systems. As the position of the amino acid side chains and the substrate peptide is virtually unchanged in all model systems, our results suggest that the exchange of Asp166 by other amino acid is less important to the phosphoryl transfer itself, but crucial to maintain the configuration of the active site in vivo. The positively charged side chain of Lys168, however, is necessary to stabilize the intermediate reaction states, particularly the side chain of the substrate peptide.     

Effects of diffusion limitation on immobilized nitrifying microorganisms at low temperatures

Activation energies of suspended and immobilized nitrifying bacteria were determined and compared to determine if diffusion limitation results in decreased sensitivity for temperature. The activation energy for the respiration activity of suspended Nitrosomonas europaea and Nitrobacter agilis was found to be 86.4 and 58.4 kJ mol(-1), respectively. The activation energy for oxygen diffusion in the support material, kappa-carrageenan, determined from the effect of temperature on the effective diffusion coefficient (D), was 17.2 kJ mol(-1). Consequently, the apparent actvation energy of diffusion limited cells should be lower. It was indeed shown that due to the effect of diffusion limitation and to temperature effects on the Monod constant K(s), the immobilized-cell activity was less sensitive to temperature. The apparent activation energy for immobilized Ns. europaea was between 28.6 and 94.2 kJ mol(-1) and for immobilized Nb. agilis between 1.4 and 72.9 kJ mol(-1), depending on the oxygen concentration and temperature. (c) 1995 John Wiley & Sons, Inc.     

Arginine 127 stabilizes the transition state in carboxypeptidase

Crystallographic studies suggest that Arg-127 is a key amino acid in the hydrolysis of peptides and esters by carboxypeptidase A. The guanidinium group of Arg-127 is hypothesized to stabilize the oxyanion of the tetrahedral intermediate formed by the attack of water on the scissile carbonyl bond. We have replaced this amino acid in rat carboxypeptidase A1 with lysine (R127K), methionine (R127M), and alanine (R127A), in order to define the role of Arg-127 in carboxypeptidase catalyzed hydrolysis. The wild-type and mutant enzymes were expressed in yeast and purified. Kinetic studies show that Arg-127 substitution decreases kcat for both ester and amide substrates, whereas Km is relatively unchanged; for R127M and R127A this corresponds to a 6 kcal/mol decrease in transition state stabilization of the rate-limiting step. The binding affinity for the phosphonate transition state analog, Cbz-Phe-Ala(P)-OAla, was decreased by 5.4 kcal/mol, whereas binding affinity for the ground state inhibitor, DL-benzylsuccinic acid, was decreased by only 1.7 kcal/mol for R127M. Electrostatic calculations employing a finite difference solution to the Poisson-Boltzmann equation predict that the positive charge of Arg-127 should stabilize the transition state by 6-8 kcal/mol. Therefore, the experimental and theoretical data suggest that the primary role of Arg-127 is stabilization of the transition state through electrostatic interaction with the oxyanion.     

Proteolytic degradation of ribonuclease A in the pretransition region of thermally and urea-induced unfolding.

The method of limited proteolysis has proven to be appropriate for the determination of unfolding rate constants (k(U)) of ribonuclease A in the transition region of thermal denaturation [Arnold, U. & Ulbrich-Hofmann, R. (1997) Biochemistry 36, 2166-2172]. The aim of the present paper was to extend this procedure to the pretransition region of thermally and urea-induced denaturation where spectroscopic methods do not allow direct measurement of k(U). The results show that the approach can be applied successfully to denaturing (free energy of unfolding Delta G < 10 kJ.mol(-1)) and to marginally native conditions (Delta G = 10-25 kJ.mol(-1)). Under moderately (Delta G = 25-30 kJ.mol(-1)) and strongly native conditions (Delta G > 30 kJ.mol(-1)), however, the determination of kU was not possible in this way as the proteolytic degradation of ribonuclease A by thermolysin or trypsin was no longer determined by global unfolding. Here, proteolysis proceeds via the native RNase A. In the presence of low concentrations of urea, the rate constants of proteolysis were, surprisingly, smaller than in the absence of urea. As the protease activity has been taken into account, this result points to a local stabilization of the RNase A molecule.     

Direct proton magnetic resonance determination of the pKa of the active center histidine in thiolsubtilisin

The serine proteases constitute a group of endopeptidases whose members owe their catalytic activity to the presence of a catalytic triad of amino acids consisting of a serine, a histidine and an aspartate. The pK(a) values for this histidine have been determined for several cases in which there is a negative charge installed at the serine to mimic the oxyanionic intermediate and related transition state for the catalytic pathway. Instances from this laboratory include (1) replacement of the serine by a cysteine in subtilisin to create a thiolate; (2) formation of monoisopropylphosphoryl-Ser 195 monoanionic phosphodiesters (in trypsin and chymotrypsin, Ser 221 in subtilisins); and (3) tetrahedral boronates formed with peptide boronic acids. The nuclear magnetic resonance (NMR) signals pertinent to this histidine, or signals indirectly reflecting the state of ionization of this histidine, have been used effectively to monitor changes in the active center ionization state. In every case studied, there is elevation of the pK(a) at the histidine when the negative charge is installed at the serine position. Herein is reported the first NMR measurement of the active center His 63 pK(a) in thiolsubtilisin Carlsberg; it is elevated by 3 units compared with the parent enzyme. Using a numerical solution (finite difference) of the Poisson-Boltzmann equation, a protein dielectric constant of 4 provides a good estimate of the experimentally observed pK(a) elevations. Very significantly, a very low protein dielectric constant (epsilon(p) = 3-5) is required in all of the comparisons, and for all three enzymes used (chymotrypsin, trypsin, and subtilisin). Finally, we discuss why the electrostatic perturbation sensed at His of the active center is more amplified by a negative charge on the Ser side than the same charge on the Asp side. A plausible explanation is that the positive charge on the imidazolium ring of the His is localized, with the N(delta 1) carrying a smaller fraction, the N(epsilon 2) carrying the bulk of the positive charge.