Active Site of alpha-Chymotrypsin Activation by Association-Desolvation

Analyzing the Reaction of Orotidine 5'-Phosphate Decarboxylase as a Way to Examine Some Key Catalytic Proposals

This study analyzes the origin of enzyme catalysis by focusing on the reaction of orotidine 5'-phosphate decarboxylase (ODCase). This reaction involves an enormous catalytic effect of 23 kcal/mol that has been attributed to reactant state destabilization associated with the use of binding energy through the so-called Circe effect. However, our early studies and subsequent key experiments have shown that the presumed effect of the binding energy (namely, the strain exerted by a bond to a phosphate group) does not contribute to the catalysis. In this study, we perform quantitative empirical valence bond calculations that reproduce the catalytic effect of ODCase and the effect of removing the phosphate side chain. The calculations demonstrate that the effect of the phosphate is due to a change in reorganization energy and should not be described as an induced fit effect. Similarly, we show that the overall catalytic effect is due to electrostatic transition state stabilization, which again reflects the smaller reorganization energy in the enzyme than in water. We also elaborate on the problems with the induced fit proposal, including the fact that it does not serve to tell us what the actual origin of the action of the catalytic effect is. In addition to the above points, we use this paper to discuss misconceptions about the meaning of the preorganization effect, as well as other misunderstandings of what is being done in consistent calculations of enzyme catalysis.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

A smallest 6 kda metalloprotease, mini-matrilysin, in living world: a revolutionary conserved zinc-dependent proteolytic domain- helix-loop-helix catalytic zinc binding domain (ZBD)

Background

The Aim of this study is to study the minimum zinc dependent metalloprotease catalytic folding motif, helix B Met loop-helix C, with proteolytic catalytic activities in metzincin super family. The metzincin super family share a catalytic domain consisting of a twisted five-stranded β sheet and three long α helices (A, B and C). The catalytic zinc is at the bottom of the cleft and is ligated by three His residues in the consensus sequence motif, HEXXHXXGXXH, which is located in helix B and part of the adjacent Met turn region. An interesting question is - what is the minimum portion of the enzyme that still possesses catalytic and inhibitor recognition?”

Methods

We have expressed a 60-residue truncated form of matrilysin which retains only the helix B-Met turn-helix C region and deletes helix A and the five-stranded β sheet which form the upper portion of the active cleft. This is only 1/4 of the full catalytic domain. The E. coli derived 6 kDa MMP-7 ZBD fragments were purified and refolded. The proteolytic activities were analyzed by Mca-Pro-Leu-Gly-Leu-Dpa-Ala-Arg-NH2 peptide assay and CM-transferrin zymography analysis. SC44463, BB94 and Phosphoramidon were computationally docked into the 3day structure of the human MMP7 ZBD and TAD and thermolysin using the docking program GOLD.

Results

This minimal 6 kDa matrilysin has been refolded and shown to have proteolytic activity in the Mca-Pro-Leu-Gly-Leu-Dpa-Ala-Arg-NH2 peptide assay. Triton X-100 and heparin are important factors in the refolding environment for this mini-enzyme matrilysin. This minienzyme has the proteolytic activity towards peptide substrate, but the hexamer and octamer of the mini MMP-7 complex demonstrates the CM-transferrin proteolytic activities in zymographic analysis. Peptide digestion is inhibited by SC44463, specific MMP7 inhibitors, but not phosphorimadon. Interestingly, the mini MMP-7 can be processed by autolysis and producing ~ 6 ~ 7 kDa fragments. Thus, many of the functions of the enzyme are retained indicating that the helix B-Met loop-helix C is the minimal functional “domain” found to date for the matrixin family.

Conclusions

The helix B-Met loop-helix C folding conserved in metalloprotease metzincin super family is able to facilitate proteolytic catalysis for specific substrate and inhibitor recognition. The autolysis processing and producing 6 kDa mini MMP-7 is the smallest metalloprotease in living world.

Comparative thermodynamic studies on substrate and product binding of O-acetylserine sulfhydrylase reveals two different ligand recognition modes

BACKGROUND

The importance of understanding the detailed mechanism of cysteine biosynthesis in bacteria is underscored by the fact that cysteine is the only sulfur donor for all cellular components containing reduced sulfur. O-acetylserine sulfhydrylase (OASS) catalyzes this crucial last step in the cysteine biosynthesis and has been recognized as an important gene for the survival and virulence of pathogenic bacteria. Structural and kinetic studies have contributed to the understanding of mechanistic aspects of OASS, but details of ligand recognition features of OASS are not available. In the absence of any detailed study on the energetics of ligand binding, we have studied the thermodynamics of OASS from Salmonella typhimurium (StOASS), Haemophilus influenzae (HiOASS), and Mycobacterium tuberculosis (MtOASS) binding to their substrate O-acetylserine (OAS), substrate analogue (methionine), and product (cysteine).

RESULTS

Ligand binding properties of three OASS enzymes are studied under defined solution conditions. Both substrate and product binding is an exothermic reaction, but their thermodynamic signatures are very different. Cysteine binding to OASS shows that both enthalpy and entropy contribute significantly to the binding free energy at all temperatures (10-30°C) examined. The analyses of interaction between OASS with OAS (substrate) or methionine (substrate analogue) revealed a completely different mode of binding. Binding of both OAS and methionine to OASS is dominated by a favorable entropy change, with minor contribution from enthalpy change (ΔH(St-Met) = -1.5 ± 0.1 kJ/mol; TΔS(St-Met) = 8.2 kJ/mol) at 20°C. Our salt dependent ligand binding studies indicate that methionine binding affinity is more sensitive to [NaCl] as compared to cysteine affinity.

CONCLUSIONS

We show that OASS from three different pathogenic bacteria bind substrate and product through two different mechanisms. Results indicate that predominantly entropy driven methionine binding is not mediated through classical hydrophobic binding, instead, may involve desolvation of the polar active site. We speculate that OASS in general, may exhibit two different binding mechanisms for recognizing substrates and products.

A midgut-specific chymotrypsin cDNA (Slctlp1) from Spodoptera litura: cloning, characterization, localization and expression analysis

Serine proteases play important roles in food digestion and immune response during insect development. A full-length cDNA (Slctlp1) encoding a chymotrypsin-like serine protease was cloned from Spodoptera litura and characterized for its cDNA structure, developmental and induced expression and localization. The deduced protein of the Slctlp1 cDNA contains a catalytic triad and a substrate specificity pocket found in most of the serine proteases. Both the transcripts and protein of Slctlp1 were predominately expressed in the midgut at the feeding stages during the larval development. Immunohistochemistry analysis indicated that the SLCTLP1 protein was predominately present in the midgut of the 6th instar feeding larvae. Starvation suppressed the expression of Slctlp1 gene and protein in 6th instar larvae and the protein expression was increased again by re-feeding the insect. The results suggest that the cloned Slctlp1 cDNA may be involved in food protein digestion at the feeding stages during larval development.

The hydrolysis of phosphate diesters in cyclohexane and acetone

The hydrolysis of phosphate diesters is one of the most difficult reactions known. Here we show that in acetone or cyclohexane, at 25 degrees C, phosphodiesters undergo hydrolysis 5 x 10(5) and 2 x 10(9)-fold more rapidly than in water, respectively, and that this rate enhancement is achieved by lowering the enthalpy of activation.

Phosphate monoester hydrolysis in cyclohexane

The hydrolysis of simple phosphate monoesters is among the most difficult reactions that are subject to catalysis by enzymes, and it has been suggested that extraction of the substrates from solvent water may contribute to the catalytic effects of phosphohydrolases. Here, we show that the tetrabutylammonium salt of neopentyl phosphate enters wet cyclohexane at concentrations sufficient to allow determination of its rate of hydrolysis. The second-order rate constant for hydrolysis of the phosphomonoester dianion is enhanced approximately 2 x 10(12)-fold by transfer from water to cyclohexane. That rate enhancement arises from an increase in the entropy of activation.

Catalytic proficiency: the unusual case of OMP decarboxylase

Enzymes are called upon to differ greatly in the difficulty of the tasks that they perform. The catalytic proficiency of an enzyme can be evaluated by comparing the second-order rate constant (kcat/Km) with the rate of the spontaneous reaction in neutral solution in the absence of a catalyst. The proficiencies of enzymes, measured in this way, are matched by their affinity constants for the altered substrate in the transition state. These values vary from approximately approximately 10(9) M(-1) for carbonic anhydrase to approximately 10(23) M(-1) for yeast orotidine 5'-phosphate decarboxylase (ODCase). ODCase turns its substrate over with a half-time of 18 ms, in a reaction that proceeds in its absence with a half-time of 78 million years in neutral solution. ODCase differs from other decarboxylases in that its catalytic activity does not depend on the presence of metals or other cofactors, or on the formation of a covalent bond to the substrate. Several mechanisms of transition state stabilization are considered in terms of ODCase crystal structures observed in the presence and absence of bound analogs of the substrate, transition state, and product. Very large connectivity effects are indicated by the results of experiments testing how transition state stabilization is affected by the truncation of binding determinants of the substrate and the active site.

How much do enzymes really gain by restraining their reacting fragments?

The steric effect, exerted by enzymes on their reacting substrates, has been considered as a major factor in enzyme catalysis. In particular, it has been proposed that enzymes catalyze their reactions by pushing their reacting fragments to a catalytic configuration which is sometimes called near attack configuration (NAC). This work uses computer simulation approaches to determine the relative importance of the steric contribution to enzyme catalysis. The steric proposal is expressed in terms of well defined thermodynamic cycles that compare the reaction in the enzyme to the corresponding reaction in water. The S(N)2 reaction of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10, which was used in previous studies to support the strain concept is chosen as a test case for this proposal. The empirical valence bond (EVB) method provides the reaction potential surfaces in our studies. The reliability and efficiency of this method make it possible to obtain stable results for the steric free energy. Two independent strategies are used to evaluate the actual magnitude of the steric effect. The first applies restraints on the substrate coordinates in water in a way that mimics the steric effect of the protein active site. These restraints are then released and the free energy associated with the release process provides the desired estimate of the steric effect. The second approach eliminates the electrostatic interactions between the substrate and the surrounding in the enzyme and in water, and compares the corresponding reaction profiles. The difference between the resulting profiles provides a direct estimate of the nonelectrostatic contribution to catalysis and the corresponding steric effect. It is found that the nonelectrostatic contribution is about -0.7 kcal/mol while the full "apparent steric contribution" is about -2.2 kcal/mol. The apparent steric effect includes about -1.5 kcal/mol electrostatic contribution. The total electrostatic contribution is found to account for almost all the observed catalytic effect ( approximately -6.1 kcal/mol of the -6.8 calculated total catalytic effect). Thus, it is concluded that the steric effect is not the major source of the catalytic power of haloalkane dehalogenase. Furthermore, it is found that the largest component of the apparent steric effect is associated with the solvent reorganization energy. This solvent-induced effect is quite different from the traditional picture of balance between the repulsive interaction of the reactive fragments and the steric force of the protein.

Computer simulations of enzyme catalysis: finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Enzymes work by solvation substitution rather than by desolvation

Considerable attention has recently been drawn to the hypothesis that enzymes catalyze their reactions by displacing solvent and creating an environment similar to the gas phase for the reacting substrates. This "desolvation hypothesis" is reexamined in this paper by defining a common reference energy for reactions in various environments. It is argued that consistent attempts to describe the actual energetics of enzymatic reactions, taking either gas phase or solution as a reference, would contradict the above hypothesis. That is, the enzyme does remove water molecules from its substrate, but substitutes these molecules for another polar environment (namely, its active site). By taking amide hydrolysis as an example, we use experimentally estimated solvation energies and analyze the reaction profile in the gas phase, in solution, and in enzyme active sites. We show that the gas-phase reaction is characterized by an enormous activation barrier (associated with forming the charged nucleophile from neutral fragments), although the nucleophilic attack is essentially barrierless. On the other hand, the enzyme and solution reactions are found to have similar reaction profiles, with a lower activation barrier for the enzymatic reaction. Presumably, the fact that previous analyses of this problem did not involve the construction of the relevant thermodynamic cycles (and quantitative estimates of the corresponding solvation energies) led to the desolvation hypothesis. Our conclusion is that enzyme active sites provide specific polar environments that do not resemble the gas phase but that are designed for electrostatic stabilization of ionic transition states and that "solvate" these states more than water does.

Glycolamide esters as biolabile prodrugs of carboxylic acid agents: synthesis, stability, bioconversion, and physicochemical properties

Benzoic acid esters of various substituted 2-hydroxyacetamides (glycolamides) were found to be hydrolyzed extremely rapidly in human plasma solutions, the half-lives of hydrolysis being less than 5 s in 50% plasma solutions for some N,N-disubstituted glycolamide esters. The rapid rate of hydrolysis could be largely attributed to cholinesterase (also called pseudocholinesterase) present in plasma. From a study of a variety of substituted glycolamide esters and structurally related esters, the most prominent structural requirement needed for a rapid rate of hydrolysis was found to be the glycolamide ester structure combined with the presence of two substituents on the amide nitrogen atom. A structural similarity of such esters with benzoylcholine, a good substrate for cholinesterase, was put forward. Esters of N,N-disubstituted glycolamides are suggested to be a useful biolabile prodrug type for several carboxylic acid agents. The esters combine a high susceptibility to undergo enzymatic hydrolysis in plasma with a high stability in aqueous solution. Furthermore, as demonstrated with the benzoic acid model esters, it is feasible to obtain ester derivatives with almost any desired water solubility or lipophilicity with retainment of marked lability to enzymatic hydrolysis.

Waterlogged molecules

Measurements of vapor pressures over their aqueous solutions indicate that organic compounds show profound differences in hydrophilic character. These differences are of such magnitude as to suggest an important role for changing solvation in determining free energy changes associated with metabolic transformations in water, and in governing structural equilibria of proteins and other large molecules in water. When two or more functional groups are present within the same solute molecule, their combined effects on its free energy of solvation are commonly additive. Striking departures from additivity, observed in certain cases, indicate the existence of special interactions between different parts of a solute molecule and the water that surrounds it. Similar considerations presumably apply to activated intermediates in the interconversion of biological materials.

Viscosity-dependent structural fluctuations in enzyme catalysis

The effect of viscosity on the rate of catalysis of carboxypeptidase A has been tested. By use of the tripeptide carbobenzoxy-l-alanyl-l-alanyl-l-alanine [Z(L-Ala)3] as substrate, it was shown that most of the effect on the hydrolysis rate caused by the presence of 30 or 40% methanol or glycerol in aqueous solution can be ascribed to a contribution of viscosity to the catalytic rate constant, kcat. Arrhenius plots of kcat in 30 and 40% glycerol or methanol are linear and almost parallel. When the rate constants are "corrected" for the viscosity of various media, the difference between the various Arrhenius plots is considerably reduced; it vanishes, within experimental error, when the effect of the dielectric constant of the solutions is taken into account as well. It is proposed that the viscosity of the medium can influence the rate-limiting step of the enzymic reaction, which is the rate of transitions over the energy barrier preceding product formation. According to the suggested mechanism, the enzyme--substrate complex can overcome this energy barrier by viscosity-dependent structural fluctuations. The quantitative agreement between the theory and the experimental results suggests that (a) due to the temperature dependence of the viscosity of the solution, the potential energy barrier of the reaction is about 5 kcal/mol lower than the observed activation energy and (b) information about the structural flexibility of the complex can be obtained by kinetic measurements.

Substrate anchoring and the catalytic power of enzymes

Evidence from published nuclear magnetic resonance studies of enzyme-substrate and enzyme-inhibitor systems systems shows that substrates are confined at the active site of the enzyme, have a relatively long residence time, and tumble in solution as an enzyme-substrate complex. The consequences of this "substrate anchoring" with regard to the catalytic power of enzymes are considered. It appears that as a result of the relatively slow motions, the probability of existence of the "activated complex" in the enzymatic reaction may be increased and thus the reaction rate may be accelerated by factors of between 10(6) and 10(9). It is estimated that the combined effect of substrate anchoring, desolvation, and charge compensation may lead to reaction-rate accelerations by factors of between 10(9) and 10(8).

Evidence for an extended active center in elastase

Large increases in the esterase, amidase, and peptidase activities of elastase are observed on increasing the length of peptide substrates. The results suggest that at least five contiguous peptide units of the substrate (four N-terminal and one C-terminal to the scissile bond) can interact concurrently with the enzyme. When analyzed in terms of an acyl-enzyme mechanism, the increased effectiveness of the enzyme with larger substrates appears mainly in the acylation reaction and, to a lesser extent, in the binding and deacylation steps.