A quantum chemical study of the catalysis for cytidine deaminase: contribution of the extra water molecule

Cytidine deaminase is known as an important enzyme responsible for the hydrolytic deamination of cytidine, which is applied as a key step to the conversion of the precursor of the cancer drug to an active form in the living body. Cytidine with water is efficiently converted to uridine with ammonia in the cleft of cytidine deaminase. In this work, the catalysis of cytidine deaminase for the hydrolytic deamination was examined using cytosine as a model of cytidine and the model molecules for the active site of cytidine deaminase by means of the quantum chemical method. We especially investigated the contribution of the water molecule from the solvent to the catalysis, because the X-ray diffraction analysis of a crystal structure has revealed the existence of the water molecule in the vicinity of the substrate bound to the active site inside the cleft. Our computations showed that the extra water molecule from the solvent has a possibility to support the catalysis of cytidine deaminase.     

Computational Insights on the Hydride and Proton Transfer Mechanisms of D-Arginine Dehydrogenase

D-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) is an amine oxidase which catalyzes the conversion of D-arginine into iminoarginine. It contains a non-covalent FAD cofactor that is involved in the oxidation mechanism. Based on substrate, solvent, and multiple kinetic isotope effects studies, a stepwise hydride transfer mechanism is proposed. It was shown that D-arginine binds to the active site of enzyme as α-amino group protonated, and it is deprotonated before a hydride ion is transferred from its α-C to FAD. Based on a mutagenesis study, it was concluded that a water molecule is the most likely catalytic base responsible from the deprotonation of α-amino group. In this study, we formulated computational models based on ONIOM method to elucidate the oxidation mechanism of D-arginine into iminoarginine using the crystal structure of enzyme complexed with iminoarginine. The calculations showed that Arg222, Arg305, Tyr249, Glu87, His 48, and two active site water molecules play key roles in binding and catalysis. Model systems showed that the deprotonation step occurs prior to hydride transfer step, and active site water molecule(s) may have participated in the deprotonation process.     

Solvent structure and hammerhead ribozyme catalysis

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.     

Theoretical investigation of the origins of catalysis of a retro-Diels-Alder reaction by antibody 10F11

The antibody 10F11 catalyzes a retro-Diels-Alder reaction that forms HNO. Deductions about the mechanism of catalysis were made by Reymond, Baumann et al. from X-ray crystal structures and from kinetic measurements for mutated antibodies. We report a study of these reactions with quantum mechanical methods and a study of the substrate and transition state binding to the active site of the antibody 10F11 using density functional theory and empirical docking methods. We have quantitated the likely contributions to catalysis of three residues identified as possible causes of catalysis: Trp H104, Phe H101, and Ser H100. Trp H104 can make a significant contribution to catalysis through dispersive interactions (pi-stacking aromatic-aromatic stabilization). On its own, Phe H101 makes only a small contribution to catalysis. When both aromatic residues are present, they act cooperatively and can make greater contributions to catalysis than expected for each residue alone. Ser H100 and the backbone NH of Phe H101 are expected to act through hydrogen bonding to speed up the reaction, but our calculations suggest that they make only a small contribution to catalysis. Reymond's studies suggest that the hydrogen-bonding network may be mediated through a water molecule in the binding site.     

How do azoles inhibit cytochrome P450 enzymes? A density functional study

To examine how azole inhibitors interact with the heme active site of the cytochrome P450 enzymes, we have performed a series of density functional theory studies on azole binding. These are the first density functional studies on azole interactions with a heme center and give fundamental insight into how azoles inhibit the catalytic function of P450 enzymes. Since azoles come in many varieties, we tested three typical azole motifs representing a broad range of azole and azole-type inhibitors: methylimidazolate, methyltriazolate, and pyridine. These structural motifs represent typical azoles, such as econazole, fluconazole, and metyrapone. The calculations show that azole binding is a stepwise mechanism whereby first the water molecule from the resting state of P450 is released from the sixth binding site of the heme to create a pentacoordinated active site followed by coordination of the azole nitrogen to the heme iron. This process leads to the breaking of a hydrogen bond between the resting state water molecule and the approaching inhibitor molecule. Although, formally, the water molecule is released in the first step of the reaction mechanism and a pentacoordinated heme is created, this does not lead to an observed spin state crossing. Thus, we show that release of a water molecule from the resting state of P450 enzymes to create a pentacoordinated heme will lead to a doublet to quartet spin state crossing at an Fe-OH(2) distance of approximately 3.0 A, while the azole substitution process takes place at shorter distances. Azoles bind heme with significantly stronger binding energies than a water molecule, so that these inhibitors block the catalytic cycle of the enzyme and prevent oxygen binding and the catalysis of substrate oxidation. Perturbations within the active site (e.g., a polarized environment) have little effect on the relative energies of azole binding. Studies with an extra hydrogen-bonded ethanol molecule in the model, mimicking the active site of the CYP121 P450, show that the resting state and azole binding structures are close in energy, which may lead to chemical equilibrium between the two structures, as indeed observed with recent protein structural studies that have demonstrated two distinct azole binding mechanisms to P450 heme.     

Mechanism of OMP decarboxylation in orotidine 5'-monophosphate decarboxylase

Despite extensive experimental and theoretical studies, the detailed catalytic mechanism of orotidine 5'-monophosphate decarboxylase (ODCase) remains controversial. In particular simulation studies using high level quantum mechanics have failed to reproduce experimental activation free energy. One common feature of many previous simulations is that there is a water molecule in the vicinity of the leaving CO2 group whose presence was only observed in the inhibitor bound complex of ODCase/BMP. Various roles have even been proposed for this water molecule from the perspective of stabilizing the transition state and/or intermediate state. We hypothesize that this water molecule is not present in the active ODCase/OMP complex. Based on QM/MM minimum free energy path simulations with accurate density functional methods, we show here that in the absence of this water molecule the enzyme functions through a simple direct decarboxylation mechanism. Analysis of the interactions in the active site indicates multiple factors contributing to the catalysis, including the fine-tuned electrostatic environment of the active site and multiple hydrogen-bonding interactions. To understand better the interactions between the enzyme and the inhibitor BMP molecule, simulations were also carried out to determine the binding free energy of this special water molecule in the ODCase/BMP complex. The results indicate that the water molecule in the active site plays a significant role in the binding of BMP by contributing approximately -3 kcal/mol to the binding free energy of the complex. Therefore, the complex of BMP plus a water molecule, instead of the BMP molecule alone, better represents the tight binding transition state analogue of ODCase. Our simulation results support the direct decarboxylation mechanism and highlight the importance of proper recognition of protein bound water molecules in the protein-ligand binding and the enzyme catalysis.     

Water-assisted reaction mechanism of monozinc beta-lactamases

Hybrid Car-Parrinello QM/MM calculations are used to investigate the reaction mechanism of hydrolysis of a common beta-lactam substrate (cefotaxime) by the monozinc beta-lactamase from Bacillus cereus (BcII). The calculations suggest a fundamental role for an active site water in the catalytic mechanism. This water molecule binds the zinc ion in the first step of the reaction, expanding the zinc coordination number and providing a proton donor adequately oriented for the second step. The free energy barriers of the two reaction steps are similar and consistent with the available experimental data. The conserved hydrogen bond network in the active site, defined by Asp120, Cys221, and His263, not only contributes to orient the nucleophile (as already proposed), but it also guides the second catalytic water molecule to the zinc ion after the substrate is bound. The hydrolysis reaction in water has a relatively high free energy barrier, which is consistent with the stability of cefotaxime in water solution. The modeled Michaelis complexes for other substrates are also characterized by the presence of an ordered water molecule in the same position, suggesting that this mechanism might be general for the hydrolysis of different beta-lactam substrates.     

Efficient intramolecular general acid catalysis of nucleophilic attack on a phosphodiester

The hydrolysis of methyl 8-dimethylamino-1-naphthyl phosphate 4 and its reactions with a representative range of nucleophiles are catalyzed by the dimethylammonium group at acidic pH with rate accelerations of the order of 106. The reaction persists up to pH 7 because the strong intramolecular hydrogen bond, which is the key to efficient general acid catalysis, is present also in the reactant. The sensitivity to the basicity of the nucleophile (Br?nsted beta(nuc) = 0.29) lies between values measured previously for mono- and triesters. The comparisons suggest that general acid catalyzed reactions of phosphate mono- or diesters with strongly basic oxyanion nucleophiles (like those derived from a serine oxygen or a bound water molecule in an enzyme active site) will be fastest when their negative charges are neutralized by protonation. Reactions with NH2OH and its N-methylated derivatives show an apparent alpha-effect, but NH2OMe reacts no faster than a primary amine of similar basicity. It is suggested that the reaction involving NH2OH as an oxygen nucleophile proceeds through the pre-equilibrium formation of the tautomer H3N+-O- as the active nucleophile: ab initio calculations support this idea.     

First computational evidence for a catalytic bridging hydroxide ion in a phosphodiesterase active site.

Phosphodiesterases are clinical targets for a variety of biological disorders, because this superfamily of enzymes regulates the intracellular concentration of cyclic nucleotides that serve as the second messengers playing a critical role in a variety of physiological processes. Understanding the structure and mechanism of a phosphodiesterase will provide a solid basis for rational design of the more efficient therapeutics. Although a three-dimensional X-ray crystal structure of the catalytic domain of human phosphodiesterase 4B2B was recently reported, it is uncertain whether a critical bridging ligand in the active site is a water molecule or a hydroxide ion. The identity of this bridging ligand is theoretically determined by performing first-principles quantum chemical calculations on models of the active site. All the results obtained indicate that this critical bridging ligand in the active site of the reported X-ray crystal structure is a hydroxide ion, rather than a water molecule, expected to serve as the nucleophile to initialize the catalytic degradation of the intracellular second messengers.     

Elucidating the mechanism for the reduction of nitrite by copper nitrite reductase--a contribution from quantum chemical studies

Density functional methods have been applied to investigate the properties of the active site of copper-containing nitrite reductases and possible reaction mechanisms for the enzyme catalysis. The results for a model of the active site indicate that a hydroxyl intermediate is not formed during the catalytic cycle, but rather a state with a protonated nitrite bound to the reduced copper. Electron affinity calculations indicate that reduction of the T2 copper site does not occur immediately after nitrite binding. Proton affinity calculations are indicative of substantial pK(a) differences between different states of the T2 site. The calculations further suggest that the reaction does not proceed until uptake of a second proton from the bulk solution. They also indicate that Asp-92 may play both a key role as a proton donor to the substrate, and a structural role in promoting catalysis. In the D92N mutant another base, presumably a nearby histidine (His-249) may take the role as the proton donor. On the basis of these model calculations and available experimental evidence, an ordered reaction mechanism for the reduction of nitrite is suggested. An investigation of the binding modes of the nitric oxide product and the nitrite substrate to the model site has also been made, indicating that nitric oxide prefers to bind in an end-on fashion to the reduced T2 site.     

Theoretical examination of Mg(2+)-mediated hydrolysis of a phosphodiester linkage as proposed for the hammerhead ribozyme

The hammerhead ribozyme is an RNA molecule capable of self-cleavage at a unique site within its sequence. Hydrolysis of this phosphodiester linkage has been proposed to occur via an in-line attack geometry for nucleophilic displacement by the 2'-hydroxyl on the adjoining phosphorus to generate a 2',3'-cyclic phosphate ester with elimination of the 5'-hydroxyl group, requiring a divalent metal ion under physiological conditions. The proposed S(N)2(P) reaction mechanism was investigated using density functional theory calculations incorporating the hybrid functional B3LYP to study this metal ion-dependent reaction with a tetraaquo magnesium (II)-bound hydroxide ion. For the Mg(2+)-catalyzed reaction, the gas-phase geometry optimized calculations predict two transition states with a kinetically insignificant, yet clearly defined, pentacoordinate intermediate. The first transition state located for the reaction is characterized by internal nucleophilic attack coupled to proton transfer. The second transition state, the rate-determining step, involves breaking of the exocyclic P-O bond where a metal-ligated water molecule assists in the departure of the leaving group. These calculations demonstrate that the reaction mechanism incorporating a single metal ion, serving as a Lewis acid, functions as a general base and can afford the necessary stabilization to the leaving group by orienting a water molecule for catalysis.     

A new understanding on how heme metabolism occurs in heme oxygenase: water-assisted oxo mechanism

Heme metabolism by heme oxygenase (HO) is investigated with quantum mechanical/molecular mechanical (QM/MM) calculations. A mechanism assisted by water is proposed: (1) an iron-oxo species and a water molecule are generated by the heterolytic cleavage of the O-O bond of an iron-hydroperoxo species in a similar way to P450-mediated reactions, (2) a hydrogen atom abstraction by the iron-oxo species from the generated water molecule and the C-O bond formation between the water molecule and the α-meso carbon take place simultaneously. The water molecule is hydrogen-bonded to the oxo ligand and to the water cluster in the active site of HO. The water cluster can control the position of the generated water molecule to ensure the regioselective oxidation of heme at the α-meso position, at the same time, can facilitate the oxidation by stabilizing a positive charge on the water molecule in the transition state. A key difference between HO and P450 is observed in the structure of the active site; Thr252 in P450 blocks the access of the water molecule to the α-meso position, and can thus suppress the undesired heme oxidation for P450.     

Active site labeling of G8 in the hairpin ribozyme: implications for structure and mechanism

There is mounting evidence that suggests that general acid/base catalysis is operative in the hairpin ribozyme, with analogy to the protein enzyme RNaseA. Nevertheless, the extent of general base catalysis as well as the identity of the specific chemical groups responsible remains the subject of some controversy. An affinity label has previously been used to alkylate histidine 12 (His12), the active general base in RNaseA. To date, no such experiment has been applied to a ribozyme. We have synthesized the analogous affinity label for the hairpin ribozyme with an electrophilic 2'-bromoacetamide group in lieu of the 2'-hydroxyl (2'OH) at the substrate cleavage site and show that guanosine 8 (G8) of the hairpin ribozyme is specifically alkylated, most likely at the N1 position. This evidence strongly implicates N1 of G8 in active site chemistry. By direct analogy to RNase A, these findings could be consistent with the hypothesis that deprotonated G8 residue functions as a general base in the hairpin ribozyme. Other mechanistic possibilities for N1 of G8 such as indirect general base catalysis mediated by a water molecule or transition state stabilization could also be consistent with our findings.     

A molecular dynamics exploration of the catalytic mechanism of yeast cytosine deaminase

Yeast cytosine deaminase (yCD), a zinc metalloenzyme of significant biomedical interest, is investigated by a series of molecular dynamics simulations in its free form and complexed with its reactant (cytosine), product (uracil), several reaction intermediates, and an intermediate analogue. Quantum chemical calculations, used to construct a model for the catalytic Zn ion with its ligands (two cysteines, a histidine, and one water) show, by comparison with crystal structure data, that the cysteines are deprotonated and the histidine is monoprotonated. The simulations suggest that Glu64 plays a critical role in the catalysis by yCD. The rotation of the Glu64 side-chain carboxyl group that can be protonated or deprotonated permits it to act as a proton shuttle between the Zn-bound water and cytosine and subsequent reaction intermediates. Free energy methods are used to obtain the barriers for these rotations, and they are sufficiently small to permit rotation on a nanosecond time scale. In the course of the reaction, cytosine reorients to a geometry to favor nucleophilic attack by a Zn-bound hydroxide. A stable position for a reaction product, ammonia, was located in the active site, and the free energy of exchange with a water molecule was evaluated. The simulations also reveal small motions of the C-terminus and the loop that contains Phe114 that may be important for reactant binding and product release.     

Molecular orbital studies of enzyme mechanisms. II. Catalytic oxidation of alcohols by liver alcohol dehydrogenase

Semiempirical (AM1) molecular orbital theory has been used to investigate the oxidation of alcohols at the active site of liver alcohol dehydrogenase (LADH). The model active site consists of a zinc dication coordinated to two methyl-mercaptans (Cys-46, Cys-176), an imidazole (His-67), and a water. An imidazole (His-51) hydrogen bonded to a hydroxy-acetate (Ser-48) forms the remote base. AM1 calculations that address the two distinct steps in the catalytic mechanism of ethanol oxidation by LADH are reported. These two steps are: (1) the deprotonation of ethanol by imidazole (His-51) via hydrogen-bonded hydroxy-acetate (Ser-48), creating a proton relay system; and (2) the rate-limiting hydride transfer step from ethanol C1 to nicotinamide adenine dinucleotide (NAD⁺), leading to product formation. Detailed calculations have been used to resolve the unsolved problems of mechanisms that have been suggested on the basis of kinetic data and crystal structures of several LADH complexes. We investigated two possible mechanisms for the deprotonation of ethanol, by zinc-bound OH^? and by direct deprotonation of zinc-bound ethanol by imidazole via hydroxyacetate (Ser-48). Our calculations show that there is no need for LADH to activate a water molecule at the active site as in many other zinc enzymes. This result agrees with experimental evidence. Our calculations also indicate that substrates are bound in an inner-sphere-pentacoordinated complex to the active site zincion. In this case, spectroscopic investigations agree with our results but crystallographic data do not. The highest activation energy is found for the hydride transfer, in agreement with the experiment. Finally, we proposed an alternative mechanism for the mode of action of LADH based upon our results. © 1993 John Wiley & Sons, Inc.     

Study on the effects of intermolecular interactions on firefly multicolor bioluminescence

Firefly luciferase exhibits a color-tuning mechanism based on pH-induced changes in the structure of the active site. These changes increase the polarity of the active site, and thus modulate the intermolecular interactions between the light emitter and active site molecules. In this study, the effects exerted by adenosine monophosphate (AMP), water molecules, and amino acids of Luciola cruciata luciferase active site on the emission wavelength of oxyluciferin were assessed by TD-DFT calculations. The redshift results mainly from decreased interaction of oxyluciferin with AMP and increased interaction of the emitter with a water molecule and Phe249. Breaking of a hydrogen bond between the benzothiazole oxygen atom with formation of a similar bond to the thiazolone oxygen atom is also instrumental.     

Electronic aspects of enzymatic catalysis

Enzymes are large aperiodic structures and this hinders both ab initio molecular orbital and Bloch-type band theory of calculations. A frontier orbital perturbation theory of catalysis is applied to enzymes. Reasons are given for proposing that the induced-fit conformational changes, essential to enzymatic catalysis, leads to an increase in the electronic eigenvalue density at the active site, enhancing the necessary catalytic orbital perturbation.     

Spectroscopic methods in bioinorganic chemistry: blue to green to red copper sites

A wide variety of spectroscopic methods are now available that provide complimentary insights into the electronic structures of transition-metal complexes. Combined with calculations, these define key bonding interactions, enable the evaluation of reaction coordinates, and determine the origins of unique spectroscopic features/electronic structures that can activate metal centers for catalysis. This presentation will summarize the contributions of a range of spectroscopic methods combined with calculations in elucidating the electronic structure of an active site using the blue copper site as an example. The contribution of electronic structure to electron-transfer reactivity will be considered in terms of anisotropic covalency, electron-transfer pathways, reorganization energy, and protein contributions to the geometric and electronic structures of blue-copper-related active sites.     

Investigation of the mechanism of the cell wall DD-carboxypeptidase reaction of penicillin-binding protein 5 of Escherichia coli by quantum mechanics/molecular mechanics calculations

Penicillin-binding protein 5 (PBP 5) of Escherichia coli hydrolyzes the terminal D-Ala-D-Ala peptide bond of the stem peptides of the cell wall peptidoglycan. The mechanism of PBP 5 catalysis of amide bond hydrolysis is initial acylation of an active site serine by the peptide substrate, followed by hydrolytic deacylation of this acyl-enzyme intermediate to complete the turnover. The microscopic events of both the acylation and deacylation half-reactions have not been studied. This absence is addressed here by the use of explicit-solvent molecular dynamics simulations and ONIOM quantum mechanics/molecular mechanics (QM/MM) calculations. The potential-energy surface for the acylation reaction, based on MP2/6-31+G(d) calculations, reveals that Lys47 acts as the general base for proton abstraction from Ser44 in the serine acylation step. A discrete potential-energy minimum for the tetrahedral species is not found. The absence of such a minimum implies a conformational change in the transition state, concomitant with serine addition to the amide carbonyl, so as to enable the nitrogen atom of the scissile bond to accept the proton that is necessary for progression to the acyl-enzyme intermediate. Molecular dynamics simulations indicate that transiently protonated Lys47 is the proton donor in tetrahedral intermediate collapse to the acyl-enzyme species. Two pathways for this proton transfer are observed. One is the direct migration of a proton from Lys47. The second pathway is proton transfer via an intermediary water molecule. Although the energy barriers for the two pathways are similar, more conformers sample the latter pathway. The same water molecule that mediates the Lys47 proton transfer to the nitrogen of the departing D-Ala is well positioned, with respect to the Lys47 amine, to act as the hydrolytic water in the deacylation step. Deacylation occurs with the formation of a tetrahedral intermediate over a 24 kcal x mol(-1) barrier. This barrier is approximately 2 kcal x mol(-1) greater than the barrier (22 kcal x mol(-1)) for the formation of the tetrahedral species in acylation. The potential-energy surface for the collapse of the deacylation tetrahedral species gives a 24 kcal x mol(-1) higher energy species for the product, signifying that the complex would readily reorganize and pave the way for the expulsion of the product of the reaction from the active site and the regeneration of the catalyst. These computational data dovetail with the knowledge on the reaction from experimental approaches.