Catalytic cycle of the phosphatidylcholine-preferring phospholipase C from Bacillus cereus. Solvent viscosity, deuterium isotope effects, and proton inventory studies

The phosphatidylcholine-preferring phospholipase C from Bacillus cereus (PLCBc) is a 28.5 kDa enzyme with three zinc ions in its active site. Although much is known about the roles that various PLCBc active site amino acids play in binding and catalysis, there is little information about the rate-determining step of the PLCBc-catalyzed hydrolysis of phospholipids and the catalytic cycle of the enzyme. To gain insight into these aspects of the hydrolysis, solvent viscosity variation experiments were conducted to determine whether an external step (substrate binding or product release) or an internal step (hydrolysis) is rate-limiting. The data indicate that the PLCBc-catalyzed reaction is unaffected by changes in solvent viscosity. This observation is inconsistent with the notion of substrate binding or product release being rate-determining and supports the hypothesis that a chemical step is rate-limiting. Furthermore, a deuterium isotope effect of 1.9 and a linear proton inventory plot indicate one proton is transferred in the rate-determining step. These data may be used to formulate a comprehensive catalytic cycle that is for the first time based on experimental evidence. In this mechanism, Asp55 of PLCBc activates an active site water molecule for attack on the phosphodiester bond, the hydrolysis of which is rate-limiting. The phosphorylcholine product is the first to leave the active site, followed by diacylglycerol.     

Influence of stereochemistry on proton transfer in protonated tripeptide models

Vectorial proton transfer among carbonyl oxygen atoms was studied in two models of tripeptide via quantum chemical calculations using the hybrid B3LYP functional and the 6-31++G** basis set. Two principal proton transfer pathways were found: a first path involving isomerization of the proton around the double bond of the carbonyl group, and a second based on the large conformational flexibility of the tripeptide model where all carbonyl oxygen atoms cooperate. The latter pathway has a rate-determining step energy barrier that is only around half of that for the first pathway. As conformational flexibility plays a crucial role in second pathway, the effect of attaching methyl groups to the alpha carbon atoms was studied. The results obtained are presented for all four possible stereochemical configurations.     

Exploring hydrogen bond in the excited state leading toward intramolecular proton transfer: detailed analysis of the structure and charge density topology along the reaction path using QTAIM

Excited state intramolecular proton transfer (ESIPT) reaction along the O-H[Symbol: see text][Symbol: see text][Symbol: see text][Symbol: see text]O hydrogen bond of o-hydroxy benzaldehyde (OHBA), methyl salicylate (MS) and salicylic acid (SA) was investigated by ab-initio quantum chemical calculation and theory of atoms and molecules (QTAIM) for the first time. Variation in several geometric as well as QTAIM parameters along the reaction coordinate was monitored in the fully relaxed excited state potential energy curve (PEC) obtained from intrinsic reaction coordinate (IRC) analysis. Although, the excited state barrier height for the forward reaction (?E (0) (#) ) reduces substantially in all the systems, MS and SA do not show any obvious asymmetry for proton transfer. For MS and SA, the crossover of the bond index as well as the lengths of the participating bonds at the saddle point is assigned due to this symmetry in accordance with bond energy - bond order (BEBO) model, which does not hold true in OHBA both in the ground and excited states. Bond ellipticity, covalent and metallic character were examined for different structures along the reaction path within the QTAIM framework. The QTAIM analysis was found to be able to uniquely distinguish between the ground and excited states of the OHBA molecule as well as both determining the effects on the bonding character of adding different substituent groups and differentiating between the ESIPT reactions in the SA and MS molecules.     

QM/MM investigation on the catalytic mechanism of Bacteroides thetaiotaomicron α-glucosidase BtGH97a

Bacteroides thetaiotaomicron α-glucosidase BtGH97a is an inverting enzyme. In this paper, the hydrolysis mechanism of p-nitro-phenyl α-d-glucopyranoside (pNP-Glc) catalyzed by BtGH97a was firstly studied by using quantum mechanical/molecular mechanical (QM/MM) approach. Two possible reaction pathways were considered. In the first pathway, a water molecule deprotonated by a nucleophilic base (here E439 or E508) attacks firstly on the anomeric carbon of pNP-Glc, then a proton from an acid residue (E532) attacks on the glycosidic oxygen to finish the hydrolysis reaction (named as nucleophilic attack-first pathway). In the second pathway, the proton from E532 attacks firstly on the glycosidic oxygen, then the water deprotonated by the nucleophilic base attacks on the anomeric carbon of pNP-Glc (named as proton attack-first pathway). Our calculation results indicate that the nucleophilic attack-first pathway is favorable in energy, in which the nucleophilic attack process is the rate-determining step with an energy barrier of 15.4kcal/mol in the case of residue E508 as nucleophilic base. In this rate-determining step, the deprotonation of water and the attack on the anomeric carbon are concerted. In the proton attack-first pathway, the proton attack on the glycosidic oxygen is the rate-determining step, and the energy barrier is 24.1kcal/mol. We conclude that the hydrolysis mechanism would follow nucleophilic attack-first pathway.     

Mechanism of proton transfer. II. Transfers to and from imidazole

A CNDO/2 calculation of an imidazole plus a proton in the plane perpendicular to the imidaxole plane has shown that the proton is lost most easilu in the direction defined by the N---H bond regardless of whether the hydrogen bond is bent or is not bent. However, bent hydrogen bonds transfer protons much easier than unstrained hydrogen bonds in agreement with the chymotrypsin mechanism of Wang.     

Slow-binding inhibition of chymotrypsin and cathepsin G by the peptide aldehyde chymostatin

The microbial, peptide-derived aldehyde chymostatin is a potent, competitive inhibitor of chymotrypsin and cathepsin G: Ki = 4 X 10(-10) and 1.5 X 10(-7) M, respectively. Et is "slow-binding inhibitor" of both proteases and, as such, allows determination of rate constants for its association with and dissociation from these proteases. Inhibition kinetics indicate second-order rate constants for the association of chymostatin with chymotrypsin and cathepsin G of 360,000 and 2000 M-1 S-1, respectively and a first-order rate constant for the dissociation of both protease-chymostatin complexes of approximately 0.0002 s-1. Thus, the extreme difference in potency of chymostatin as an inhibitor of chymotrypsin and cathepsin G originates entirely in Kon. Solvent deuterium isotope effects (SIE) were determined to probe the reaction step that rate limits Kon. For the reaction of chymotrypsin with chymostatin, the SIE for Kon is 1.6 +/- 0.1, while for the reaction of chymotrypsin with the peptide substrates Ala-Ala-Phe-pNA and Suc-Ala-Ala-Pro-Phe-pNA, the SIE's for Kc/Km are 2.8 +/- 0.2 and 1.9 +/- 0.1, respectively. These results suggest that Kon for the association of chymotrypsin with chymostatin is at least partially rate limited by a reaction step involving proton transfer. Combined with results for the inhibition of chymotrypsin by Bz-Phe-H [Kennedy, W.P., & Schultz, R. M. (1979) Biochemistry 18, 349-356], these data suggest a mechanism for inhibition by chymostatin involving the general-base-catalyzed formation of an enzyme-bound hemiacetal, followed by a conformational change of this intermediate that produces the final, stable complex of enzyme and inhibitor.(ABSTRACT TRUNCATED AT 250 WORDS)     

QM/MM investigation on the catalytic mechanism of Bacteroides thetaiotaomicron α-glucosidase BtGH97a

Bacteroides thetaiotaomicron α-glucosidase BtGH97a is an inverting enzyme. In this paper, the hydrolysis mechanism of p-nitro-phenyl α-d-glucopyranoside (pNP-Glc) catalyzed by BtGH97a was firstly studied by using quantum mechanical/molecular mechanical (QM/MM) approach. Two possible reaction pathways were considered. In the first pathway, a water molecule deprotonated by a nucleophilic base (here E439 or E508) attacks firstly on the anomeric carbon of pNP-Glc, then a proton from an acid residue (E532) attacks on the glycosidic oxygen to finish the hydrolysis reaction (named as nucleophilic attack-first pathway). In the second pathway, the proton from E532 attacks firstly on the glycosidic oxygen, then the water deprotonated by the nucleophilic base attacks on the anomeric carbon of pNP-Glc (named as proton attack-first pathway). Our calculation results indicate that the nucleophilic attack-first pathway is favorable in energy, in which the nucleophilic attack process is the rate-determining step with an energy barrier of 15.4 kcal/mol in the case of residue E508 as nucleophilic base. In this rate-determining step, the deprotonation of water and the attack on the anomeric carbon are concerted. In the proton attack-first pathway, the proton attack on the glycosidic oxygen is the rate-determining step, and the energy barrier is 24.1 kcal/mol. We conclude that the hydrolysis mechanism would follow nucleophilic attack-first pathway.     

Short-strong hydrogen bonds and a low barrier transition state for the proton transfer reaction in RNase A catalysis: a quantum chemical study

There is growing evidence that some enzymes catalyze reactions through the formation of short-strong hydrogen bonds as first suggested by Gerlt and Gassman. Support comes from several experimental and quantum chemical studies that include correlation energies on model systems. In the present study, the process of proton transfer between hydroxyl and imidazole groups, a model of the crucial step in the hydrolysis of RNA by the enzymes of the RNase A family, is investigated at the quantum mechanical level of density functional theory and perturbation theory at the MP2 level. The model focuses on the nature of the formation of a complex between the important residues of the protein and the hydroxyl group of the substrate. We have also investigated different configurations of the ground state that are important in the proton transfer reaction. The nature of bonding between the catalytic unit of the enzyme and the substrate in the model is investigated by Bader's atoms in molecule theory. The contributions of solvation and vibrational energies corresponding to the reactant, the transition state and the product configurations are also evaluated. Furthermore, the effect of protein environment is investigated by considering the catalytic unit surrounded by complete proteins--RNase A and Angiogenin. The results, in general, indicate the formation of a short-strong hydrogen bond and the formation of a low barrier transition state for the proton transfer model of the enzyme.     

Elucidation of the mechanism of selenoprotein glutathione peroxidase (GPx)-catalyzed hydrogen peroxide reduction by two glutathione molecules: a density functional study

The mechanism of the hydrogen peroxide reduction by two molecules of glutathione catalyzed by the selenoprotein glutatione peroxidase (GPx) has been computationally studied. It has been shown that the first elementary reaction of this process, (E-SeH) + H(2)O(2) --> (E-SeOH) + H(2)O (1), proceeds via a stepwise pathway with the overall barrier of 17.1 kcal/mol, which is in good agreement with the experimental barrier of 14.9 kcal/mol. During reaction 1, the Gln83 residue has been found to play a key role as a proton acceptor, which is consistent with experiments. The second elementary reaction, (E-SeOH) + GSH --> (E-Se-SG) + HOH (2), proceeds with the barrier of 17.9 kcal/mol. The last elementary reaction, (E-Se-SG) + GSH --> (E-SeH) + GS-SG (3), is initiated with the coordination of the second glutathione molecule. The calculations clearly suggest that the amide backbone of the Gly50 residue directly participates in this reaction and the presence of two water molecules is absolutely vital for the reaction to occur. This reaction proceeds with the barrier of 21.5 kcal/mol and is suggested to be a rate-determining step of the entire GPx-catalyzed reaction H(2)O(2) + 2GSH --> GS-SG + 2H(2)O. The results discussed in the present study provide intricate details of every step of the catalytic mechanism of the GPx enzyme and are in good general agreement with experimental findings and suggestions.     

Density functional theory studies on hydroxylamine mechanism of cyclohexanone ammoximation on titanium silicalite-1 catalyst

The hydroxylamine mechanism of cyclohexanone ammoximation on defective titanium active site of titanium silicalite-1 (TS-1) was simulated using two-layer ONIOM (M062X/6-31G**:PM6) method. A new energy favorable reaction route was found, which contained two parts: (1) the catalytic oxidation of adsorbed NH₃ to form hydroxylamine using the Ti-OOH as an active oxidant formed by reacting H₂O₂ with the defective Ti active site; (2) the subsequent noncatalytic oximation of desorbed hydroxylamine and cyclohexanone out of TS-1 pores to form cyclohexanone oxime. In the catalytic formation of hydroxylamine on the Ti active site of TS-1, the proposed mechanism of two-step single-proton transfer aided by a lattice oxygen atom bonded to Ti atom need a lower reaction energy than the mechanism proposed before. In the subsequent noncatalytic oximation of hydroxylamine and cyclohexanone, which contained two elementary reaction steps in total, the mechanisms of one-step double-proton transfer in the first elementary reaction step and the subsequent one-step three-proton transfer for the second elementary reaction step were proposed, in which the solvent water molecules played a very important role in assisting and stabilizing the proton transfer processes.     

The mechanism of the proton transfer: an outline

A brief summary of the principal notions of the quantum-mechanical theory of the charge transfer reactions has been presented. In the framework of this theory, the mechanism of the proton transfer consists in the classical medium reorganization that equalizes the proton energy levels in the initial and final states, and a consequent proton transfer via a quantum-mechanical underbarrier transition. On the basis of this mechanism, factors influencing the proton transfer probability, and hence kinetic isotope effect, have been discussed; among them are the optimum tunneling distance, the involvement of the excited vibrational states, etc. Semi-classical and quantum-mechanical treatments of the Swain-Schaad relations have been compared. Some applications to enzymatic proton-transfer reactions have been described.     

Novel mechanism of hydrolysis of therapeutic beta-lactams by Stenotrophomonas maltophilia L1 metallo-beta-lactamase.

Stopped-flow tryptophan fluorescence under single turnover and pseudo-first-order conditions has been used to investigate the kinetic mechanism of beta-lactam hydrolysis by the Stenotrophomonas maltophilia L1 metallo-beta-lactamase. For the cephalosporin substrates nitrocefin and cefaclor and the carbapenem meropenem, a substantial quench of fluorescence is observed on association of substrate with enzyme. We have assigned this to a rearrangement event subsequent to formation of an initial collision complex. For the colorimetric compound nitrocefin, decay of this dark inter- mediate represents the overall rate-determining step for the reaction and is equivalent to decay of a previously observed state in which the beta-lactam amide bond has already been cleaved. For both cefaclor and meropenem, the rate-determining step for hydrolysis is loss of a second, less quenched state, in which, however, the beta-lactam amide bond remains intact. We suggest, therefore, that the mechanism of hydrolysis of nitrocefin by binuclear metallo-beta-lactamases may be atypical and that cleavage of the beta-lactam amide bond is the rate-determining step for breakdown of the majority of beta-lactam substrates by the L1 enzyme.     

Quantum chemical studies for oxidation of morpholine by Cytochrome P450

Since morpholine oxidation has recently been shown to involve Cytochrome P450, the study on its mechanism at molecular level using quantum chemical calculations for the model of cytochrome active site is reported here. The reaction pathway is investigated for two electronic states, the doublet and the quartet, by means of density functional theory. The results show that morpholine hydroxylation occurs through hydrogen atom abstraction and rebound mechanism. However, in the low spin state, the reaction is concerted and hydrogen atom abstraction yields directly ferric-hydroxy morpholine complex without a distinct rebound step while in quartet state the reaction is stepwise. The presence of nitrogen in a morpholine heterocycle is postulated to greatly facilitate hydrogen abstraction. The hydroxylated product undergoes intramolecular hydrogen atom transfer from hydroxy group to nitrogen, leading to the cleavage of the C-N bond and the formation of 2-(2-aminoethoxy) acetaldehyde. The cleavage of the C-N bond is indicated as the rate-determining step for the studied reaction. The assistance of explicit water molecule is shown to lower the energy barrier for the C-N bond cleavage in enzymatic environment whereas solvent effects mimicked by COSMO solvent model have minor influence on relative energies along the pathway.     

Differences in the catalytic mechanisms of mesophilic and thermophilic indole-3-glycerol phosphate synthase enzymes at their adaptive temperatures

Thermophilic enzymes tend to be less catalytically-active at lower temperatures relative to their mesophilic counterparts, despite having very similar crystal structures. An often cited hypothesis for this general observation is that thermostable enzymes have evolved a more rigid tertiary structure in order to cope with their more extreme, natural environment, but they are also less flexible at lower temperatures, leading to their lower catalytic activity under mesophilic conditions. An alternative hypothesis, however, is that complementary thermophilic-mesophilic enzyme pairs simply operate through different evolutionary-optimized catalytic mechanisms. In this communication, we present evidence that while the steps of the catalytic mechanisms for mesophilic and thermophilic indole-3-glycerol phosphate synthase (IGPS) enzymes are fundamentally similar, the identity of the rate-determining step changes as a function of temperature. Our findings indicate that while product release is rate-determining at 25°C for thermophilic IGPS, near its adaptive temperature (75°C), a proton transfer event, involving a general acid, becomes rate-determining. The rate-determining steps for thermophilic and mesophilic IGPS enzymes are also different at their respective, adaptive temperatures with the mesophilic IGPS-catalyzed reaction being rate-limited before irreversible CO2 release, and the thermophilic IGPS-catalyzed reaction being rate limited afterwards.     

A theoretical model of the catalytic mechanism of the Delta5-3-ketosteroid isomerase reaction

The present paper describes a theoretical approach to the catalytic reaction mechanism involved in the conversion of 5-androstene-3,17-dione to 4-androstene-3,17-dione. The model incorporates the side chains of the residues tyrosine (Tyr(14)), aspartate (Asp(38)) and aspartic acid (Asp(99)) of the enzyme Delta(5)-3-ketosteroid isomerase (KSI; EC 5.3.3.1). The reaction involves two steps: first, Asp(38) acts as a base, abstracting the 4beta-H atom (proton) from C-4 of the steroid to form a dienolate as the intermediate; next, the intermediate is reketonized by proton transfer to the 6beta-position. Each step goes through its own transition state. Functional groups of the Tyr(14) and Asp(99) side chains act as hydrogen bond donors to the O1 atom of the steroid, providing stability along the reaction coordinate. Calculations were assessed at high level Hartree-Fock theory, using the 6-31G(*) basis set and the most important physicochemical properties involved in each step of the reaction, such as total energy, hardness, and dipole moment. Likewise, to explain the mechanism of reaction, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), atomic orbital contributions to frontier orbitals formation, encoded electrostatic potentials, and atomic charges were used. Energy minima and transition state geometries were confirmed by vibrational frequency analysis. The mechanism described herein accounts for all of the properties, as well as the flow of atomic charges, explaining both catalytic mechanism and proficiency of KSI.     

Photosystem II: The machinery of photosynthetic water splitting

This review summarizes our current state of knowledge on the structural organization and functional pattern of photosynthetic water splitting in the multimeric Photosystem II (PS II) complex, which acts as a light-driven water: plastoquinone-oxidoreductase. The overall process comprises three types of reaction sequences: (1) photon absorption and excited singlet state trapping by charge separation leading to the ion radical pair [Formula: see text] formation, (2) oxidative water splitting into four protons and molecular dioxygen at the water oxidizing complex (WOC) with P680+* as driving force and tyrosine Y(Z) as intermediary redox carrier, and (3) reduction of plastoquinone to plastoquinol at the special Q(B) binding site with Q(A)-* acting as reductant. Based on recent progress in structure analysis and using new theoretical approaches the mechanism of reaction sequence (1) is discussed with special emphasis on the excited energy transfer pathways and the sequence of charge transfer steps: [Formula: see text] where (1)(RC-PC)* denotes the excited singlet state (1)P680* of the reaction centre pigment complex. The structure of the catalytic Mn(4)O(X)Ca cluster of the WOC and the four step reaction sequence leading to oxidative water splitting are described and problems arising for the electronic configuration, in particular for the nature of redox state S(3), are discussed. The unravelling of the mode of O-O bond formation is of key relevance for understanding the mechanism of the process. This problem is not yet solved. A multistate model is proposed for S(3) and the functional role of proton shifts and hydrogen bond network(s) is emphasized. Analogously, the structure of the Q(B) site for PQ reduction to PQH(2) and the energetic and kinetics of the two step redox reaction sequence are described. Furthermore, the relevance of the protein dynamics and the role of water molecules for its flexibility are briefly outlined. We end this review by presenting future perspectives on the water oxidation process.     

Violet bioluminescence and fast kinetics from W92F obelin: structure-based proposals for the bioluminescence triggering and the identification of the emitting species

Obelin from the hydroid Obelia longissima and aequorin are members of a subfamily of Ca(2+)-regulated photoproteins that is a part of the larger EF-hand calcium binding protein family. On the addition of Ca(2+), obelin generates a blue bioluminescence emission (lambda(max) = 485 nm) as the result of the oxidative decarboxylation of the bound substrate, coelenterazine. The W92F obelin mutant is noteworthy because of the unusually high speed with which it responds to sudden changes of [Ca(2+)] and because it emits violet light rather than blue due to a prominent band with lambda(max) = 405 nm. Increase of pH in the range from 5.5 to 8.5 and using D(2)O both diminish the contribution of the 405 nm band, indicating that excited state proton transfer is involved. Fluorescence model studies have suggested the origin of the 485 nm emission as the excited state of an anion of coelenteramide, the bioluminescence reaction product, and 405 nm from the excited neutral state. Assuming that the dimensions of the substrate binding cavity do not change during the excited state formation, a His22 residue within hydrogen bonding distance to the 6-(p-hydroxy)-phenyl group of the excited coelenteramide is a likely candidate for accepting the phenol proton to produce an ion-pair excited state, in support of recent suggestions for the bioluminescence emitting state. The proton transfer could be impeded by removal of the Trp92 H-bond, resulting in strong enhancement of a 405 nm band giving the violet color of bioluminescence. Comparative analysis of 3D structures of the wild-type (WT) and W92F obelins reveals that there are structural displacements of certain key Ca(2+)-ligating residues in the loops of the two C-terminal EF hands as well as clear differences in hydrogen bond networks in W92F. For instance, the hydrogen bond between the side-chain oxygen atom of Asp169 and the main-chain nitrogen of Arg112 binds together the incoming alpha-helix of loop III and the exiting alpha-helix of loop IV in WT, providing probably concerted changes in these EF hands on calcium binding. But this linkage is not found in W92F obelin. These differences apparently do not change the overall affinity to calcium of W92F obelin but may account for the kinetic differences between the WT and mutant obelins. From analysis of the hydrogen bond network in the coelenterazine binding cavity, it is proposed that the trigger for bioluminescence reaction in these Ca(2+)-regulated photoproteins may be a shift of the hydrogen bond donor-acceptor separations around the coelenterazine-2-hydroperoxy substrate, initiated by small spatial adjustment of the exiting alpha-helix of loop IV.     

On the inactivity of thiol-subtilisin. The role of the intramolecular electric field

Based on computed proton affinities for several model systems, the energetics of proton transfer and the acidity of the catalytic triads Cys-His-Asn (papain). Cys-His-Asp (thiol-subtilisin) and Ser-His-Asp (subtilisin) are discussed. It is shown that in papain the ion-pair Cys--HisH+ exists owing to the intramolecular electric field, and that a similar situation is found in thiol-subtilisin. but not in subtilisin. Assuming similar reaction mechanisms for papain and thiol-subtilisin - i.e. proton transfer from HisH+ to the NH group of the scissile peptide bond - the inactivity of thil-subtilisin towards proteins is explained by the much greater basicity of His in the complex His-Asp- than in His-Asn. In order for this explanation to be consistent, it is tentatively concluded that the catalytic mechanism of the serine proteases is different from that of the cystein proteases, and involves direct transfer of the serine proton to the leaving group in the acylation step.     

On mechanisms of fluorescence quenching by water

Mechanisms of fluorescence quenching of aromatic chromophores by water are reviewed. The mechanisms include polarity of chromophore environment, proton or electron transfer between the excited chromophore and water. A hypothesis is proposed that the quenching can be a result of chromophore-solvent hydrogen bond breaking in the excited state.