Ab initio QM/MM study of class A beta-lactamase acylation: dual participation of Glu166 and Lys73 in a concerted base promotion of Ser70

Beta-lactamase acquisition is the most prevalent basis for Gram-negative bacteria resistance to the beta-lactam antibiotics. The mechanism used by the most common class A Gram-negative beta-lactamases is serine acylation followed by hydrolytic deacylation, destroying the beta-lactam. The ab initio quantum mechanical/molecular mechanical (QM/MM) calculations, augmented by extensive molecular dynamics simulations reported herein, describe the serine acylation mechanism for the class A TEM-1 beta-lactamase with penicillanic acid as substrate. Potential energy surfaces (based on approximately 350 MP2/6-31+G calculations) reveal the proton movements that govern Ser70 tetrahedral formation and then collapse to the acyl-enzyme. A remarkable duality of mechanism for tetrahedral formation is implicated. Following substrate binding, the pathway initiates by a low energy barrier (5 kcal mol(-1)) and an energetically favorable transfer of a proton from Lys73 to Glu166, through the catalytic water molecule and Ser70. This gives unprotonated Lys73 and protonated Glu166. Tetrahedral formation ensues in a concerted general base process, with Lys73 promoting Ser70 addition to the beta-lactam carbonyl. Moreover, the three-dimensional potential energy surface also shows that the previously proposed pathway, involving Glu166 as the general base promoting Ser70 through a conserved water molecule, exists in competition with the Lys73 process. The existence of two routes to the tetrahedral species is fully consistent with experimental data for mutant variants of the TEM beta-lactamase.     

The fragmentation-recombination mechanism of the enzyme glutamate mutase studied by QM/MM simulations

The radical mechanism of the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is studied with quantum mechanical/molecular mechanical (QM/MM) simulations based on density functional theory (DFT/MM). The hydrogen transfer between the substrate and the cofactor is found to be rate limiting with a barrier of 101.1 kJ mol(-1). A careful comparison to the uncatalyzed reaction in water is performed. The protein influences the reaction predominantly electrostatically and to a lesser degree sterically. Our calculations shed light on the atomistic details of the reaction mechanism. The well-known arginine claw and Glu 171 ( Clostridium cochlearium notation) are found to have the strongest influence on the reaction. However, a catalytic role of Glu 214, Lys 322, Gln 147, Glu 330, Lys 326, and Met 294 is found as well. The arginine claw keeps the intermediates in place and is probably responsible for the enantioselectivity. Glu 171 temporarily accepts a proton from the glutamyl radical intermediate and donates it back at the end of the reaction. We relate our results to experimental data when available. Our simulations lead to further understanding of how glutamate mutase catalyzes the carbon skeleton rearrangement of glutamate.     

A theoretical study on the catalytic mechanism of Mus musculus adenosine deaminase

The catalytic mechanism of Mus musculus adenosine deaminase (ADA) has been studied by quantum mechanics and two‐layered ONIOM calculations. Our calculations show that the previously proposed mechanism, involving His238 as the general base to activate the Zn‐bound water, has a high activation barrier of about 28 kcal/mol at the proposed rate‐determining nucleophilic addition step, and the corresponding calculated kinetic isotope effects are significantly different from the recent experimental observations. We propose a revised mechanism based on calculations, in which Glu217 serves as the general base to abstract the proton of the Zn‐bound water, and the protonated Glu217 then activates the substrate for the subsequent nucleophilic addition. The rate‐determining step is the proton transfer from Zn‐OH to 6‐NH₂ of the tetrahedral intermediate, in which His238 serves as a proton shuttle for the proton transfer. The calculated kinetic isotope effects agree well with the experimental data, and calculated activation energy is also consistent with the experimental reaction rate. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical study toward understanding the catalytic mechanism of pyruvate decarboxylase

Density functional calculations are employed to explore the mechanisms of all elementary reaction steps involved in the catalytic cycle of pyruvate decarboxylase (PDC). Different models are constructed for mimicking the involvement of some key residues in a certain step. The effect of the protein framework on the potential energy profiles of active site models is approximately modeled by fixing some freedoms, based on the crystal structure of the PDC enzyme from Saccharomyces cerevisiae (ScPDC). Our calculations confirm that Glu51 is the most important residue in the formation of the ylide and the release of acetaldehyde via the proton relay between Glu51, N1', and the 4'-amino group of thiamine diphosphate. The presence of Glu477 and Asp28 residues makes the decarboxylation of lactylthiamin diphosphate (LThDP) an endothermic process with a significant free energy barrier. The protonation of the alpha-carbanion to form 2-(1-hydroxyethyl)-thiamin diphosphate is found to go through a concerted double proton transfer transition state involving both Asp28 and His115 residues. The final step, acetaldehyde release, is likely to proceed through a concerted transition state involving carbon-carbon bond-breaking and the deprotonation of the alpha-hydroxyl group. The decarboxylation of LThDP and the protonation of the alpha-carbanion are two rate-limiting steps, relative to the facile occurrence of the ylide formation and acetaldehyde release. The catalytic roles of residues Glu51, Glu477, Asp28, and Gly417 in the active site of ScPDC in individual steps elucidated from the present study are in good agreement with those derived from site-directed mutagenesis.     

Catalytic mechanism of S-ribosylhomocysteinase (LuxS): direct observation of ketone intermediates by 13C NMR spectroscopy

S-Ribosylhomocysteinase (LuxS) catalyzes the cleavage of the thioether linkage of S-ribosylhomocysteine (SRH) to produce l-homocysteine and 4,5-dihydroxy-2,3-pentanedione (DHPD). This is a key step in the biosynthetic pathway of the type II autoinducer (AI-2) in both Gram-positive and Gram-negative bacteria. Previous studies demonstrated that LuxS contains a catalytically essential Fe2+ ion. The catalytic mechanism of LuxS was investigated using 2- and 3-13C-labeled SRH as substrate and 13C NMR spectroscopy. These studies revealed the presence of 2- and 3-keto intermediates in the catalytic pathway. The 2-keto intermediate was chemically synthesized, and its chemical and kinetic competence was demonstrated. The results support a catalytic mechanism in which the metal ion catalyzes an internal redox reaction, shifting the carbonyl group from the C-1 position to the C-3 position. Subsequent beta-elimination at the C-4 and C-5 positions releases homocysteine as a free thiol. The results also suggest that Cys-84 and Glu-57 are the possible general acids/bases for proton transfer during catalysis and that the keto intermediates are released from the enzyme active site before rebinding and completion of the reaction.     

Ab initio model study on acetylcholinesterase catalysis: potential energy surfaces of the proton transfer reactions

Ab initio molecular orbital (MO) and hybrid density functional theory (DFT) calculations have been applied to the initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE), which is the nucleophiric addition of Ser200 in catalytic triads to a neurotransmitter acetylcholine (ACh). We focus our attention mainly on the effects of oxyanion hole and Glu327 on the potential energy surfaces (PESs) for the proton transfer reactions in the catalytic triad Ser200-His440-Glu327. The activation barrier for the addition reaction of Ser200 to ACh was calculated to be 23.4 kcal/mol at the B3LYP/6-31G(d)//HF/3-21G(d) level of theory. The barrier height under the existence of oxyanion hole, namely, Ser200-His440-Glu327-ACh-(oxyanion hole) system, decreased significantly to 14.2 kcal/mol, which is in reasonable agreement with recent experimental value (12.0 kcal/mol). Removal of Glu327 from the catalytic triad caused destabilization of both energy of transition state for the reaction and tetrahedral intermediate (product). PESs calculated for the proton transfer reactions showed that the first proton transfer process is the most important in the stabilization of tetrahedral intermediate complex. The mechanism of addition reaction of ACh was discussed on the basis of theoretical results.     

A Theoretical Investigation on N7(H)→N9(H)Tautomerization of Isolated and Monohydrated 2,6-Dithiopurine Combined with IPCM Solvent Model

The tautomerization reaction mechanism has been reported between N7(H) and N9(H) of isolated and monohydrated 2,6‐dithiopurine using B3LYP/6‐311+G(d,p). The isodensity polarized continuum model (IPCM) in the self‐consistent reaction field (SCRF) method is employed to account for the solvent effect of water on the tautomerization reaction activation energies. The results show that the two pathways P(1) (via the carbene intermediate I1) and P(2) (via the sp³‐hybrid intermediate I2) are found in intramolecular proton transfer, and each pathway is composed by two primary steps. The calculated activation energy barriers of the rate‐determining steps in isolated 2,6‐dithiopurine N7(H)→N9(H) tautomerism are 308.2 and 220.0 kJ·mol^?1 in the two pathways, respectively. Interestingly, in one‐water molecule catalyst, it dramatically lowers the N7(H)→N9(H) energy barriers by the concerted double proton transfer mechanism in P(1), favoring the formation of 2,6‐dithiopurine N9(H). However, the single proton transfer mechanism assisted with out‐of‐plane water in the first step of P(2) increases the activation energy barrier from 220.0 to 232.3 kJ·mol^?1, while the second step is the out‐of‐plane concerted double proton transfer mechanism, indicating that they will be less preferable for proton transfer. Additionally, the results also show that all the pathways are put into the aqueous solution, and the activation energy barriers have no significant changes. Therefore, the long‐range electrostatic effect of bulk solvent has no significant impact on proton transfer reactions and the interaction with explicit water molecules will significantly influence proton transfer reactions.     

Computational study of phosphatase activity in soluble epoxide hydrolase: high efficiency through a water bridge mediated proton shuttle

Recently, a new branch of fatty acid metabolism has been opened by the novel phosphatase activity found in the N-terminal domain of the, hence bifunctional, soluble epoxide hydrolase (sEH). Importantly, this finding has also provided a new site for drug targeting in sEH's activity regulation. Classical MD and hybrid Car-Parrinello QM/MM calculations have been performed to investigate the reaction mechanism of the phosphoenzyme intermediate formation in the first step of the catalysis. The results support a concerted multi-event reaction mechanism: (1) a dissociative in-line nucleophilic substitution for the phosphoryl transfer reaction; (2) a double proton transfer involved in the formation of a good leaving group in the transition state. The presence of a water bridge in the substrate/enzyme complex allowed an efficient proton shuttle, showing its key role in speeding up the catalysis. The calculated free energy of the favored catalytic pathway is approximately 19 kcal/mol, in excellent agreement with experimental data.     

Evaluation of the catalytic mechanism of AICAR transformylase by pH-dependent kinetics, mutagenesis, and quantum chemical calculations.

The catalytic mechanism of 5-aminoimidazole-4-carboxamide ribonucleotide transformylase (AICAR Tfase) is evaluated with pH dependent kinetics, site-directed mutagenesis, and quantum chemical calculations. The chemistry step, represented by the burst rates, was not pH-dependent, which is consistent with our proposed mechanism that the 4-carboxamide of AICAR assists proton shuttling. Quantum chemical calculations on a model system of 5-amino-4-carboxamide imidazole (AICA) and formamide using the B3LYP/6-31G level of theory confirmed that the 4-carboxamide participated in the proton-shuttling mechanism. The result also indicated that the amide-assisted mechanism is concerted such that the proton transfers from the 5-amino group to the formamide are simultaneous with nucleophilic attack by the 5-amino group. Because the process does not lead to a kinetically stable intermediate, the intramolecular proton transfer from the 5-amino group through the 4-carboxamide to the formamide proceeds in the same transition state. Interestingly, the calculations predicted that protonation of the N3 of the imidazole of AICA would reduce the energy barrier significantly. However, the pK(a) of the imidazole of AICAR was determined to be 3.23 +/- 0.01 by NMR titration, and AICAR is likely to bind to the enzyme with its imidazole in the free base form. An alternative pathway was suggested by modeling Lys266 to have a hydrogen-bonding interaction with the N3 of the imidazole of AICAR. Lys266 has been implicated in catalysis based on mutagenesis studies and the recent X-ray structure of AICAR Tfase. The quantum chemical calculations on a model system that contains AICA complexed with CH3NH3+ as a mimic of the Lys residue confirmed that such an interaction lowered the activation energy of the reaction and likewise implicated the 4-carboxamide. To experimentally verify this hypothesis, we prepared the K266R mutant and found that its kcat is reduced by 150-fold from that of the wild type without changes in substrate and cofactor Km values. The kcat-pH profile indicated virtually no pH-dependence in the pH range 6-10.5. The results suggest that the ammonium moiety of Lys or Arg is important in catalysis, most likely acting as a general acid catalyst with a pK(a) value greater than 10.5. The H267A mutant was also prepared since His267 has been found in the active site and implicated in catalysis. The mutant enzyme showed no detectable activity while retaining its binding affinity for substrate, indicating that it plays a critical role in catalysis. We propose that His267 interacts with Lys266 to aid in the precise positioning of the general acid catalyst to the N3 of the imidazole of AICAR.     

Concerted or stepwise hydrogen transfer in the transfer hydrogenation of acetophenone catalyzed by ruthenium-acetamido complex: a theoretical mechanistic investigation

In this paper, the mechanism of transfer hydrogenation of acetophenone catalyzed by ruthenium-acetamido complex was studied using density function theory (DFT) method. The catalytic cycle of transfer hydrogenation consists of hydrogen transfer (HT) step and dehydrogenation (DH) step of isopropanol (IPA). Inner sphere mechanism (paths 1 and 7) and outer sphere mechanism (paths 2-6) in HT step are fully investigated. Calculated results indicate that DH step of IPA (from (i)1 to (i)2) is the rate-determining step in the whole catalytic cycle, which has a potential energy barrier of 16.2 kcal/mol. On the other hand, the maximum potential energy barriers of paths 1-7 in the HT step are 5.9, 12.7, 24.4, 16.8, 23.7, 7.2, and 6.1 kcal/mol, respectively. The inner sphere pathways (paths 1 and 7) are favorable hydrogen transfer modes compared with outer sphere pathways, and the proton transferred to the oxygen atom of acetophenone comes from the hydroxyl group but not from amino group of acetamido ligand. Those theoretical results are in agreement with experimental report. However, in view of this DFT study in the inner sphere mechanism of HT step, hydride transfer and proton transfer are concerted and asynchronous hydrogen transfer but not a stepwise one, and hydride transfer precedes proton transfer in this case.     

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Extensive combined quantum mechanical (B3LYP/6‐31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton‐transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2‐pyrimidinone. The rate‐determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.     

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.     

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.     

A QM/MM simulation study of transamination reaction at the active site of aspartate aminotransferase: Free energy landscape and proton transfer pathways

Transaminase is a key enzyme for amino acid metabolism, which reversibly catalyzes the transamination reaction with the help of PLP (pyridoxal 5^' -phosphate) as its cofactor. Here we have investigated the mechanism and free energy landscape of the transamination reaction involving the aspartate transaminase (AspTase) enzyme and aspartate-PLP (Asp-PLP) complex using QM/MM simulation and metadynamics methods. The reaction is found to follow a stepwise mechanism where the active site residue Lys258 acts as a base to shuttle a proton from α -carbon (CA) to imine carbon (C4A) of the PLP-Asp Schiff base. In the first step, the Lys258 abstracts the CA proton of the substrate leading to the formation of a carbanionic intermediate which is followed by the reprotonation of the Asp-PLP Schiff base at C4A atom by Lys258. It is found that the free energy barrier for the proton abstraction by Lys258 and that for the reprotonation are 17.85 and 3.57 kcal/mol, respectively. The carbanionic intermediate is 7.14 kcal/mol higher in energy than the reactant. Hence, the first step acts as the rate limiting step. The present calculations also show that the Lys258 residue undergoes a conformational change after the first step of transamination reaction and becomes proximal to C4A atom of the Asp-PLP Schiff base to favor the second step. The active site residues Tyr70* and Gly38 anchor the Lys258 in proper position and orientation during the first step of the reaction and stabilize the positive charge over Lys258 generated at the intermediate step.     

Class I ribonucleotide reductase revisited: the effect of removing a proton on Glu441

The substrate mechanism of class I ribonucleotide reductase has been revisited using the hybrid density functional B3LYP method. The molecular model used is based on the X-ray structure and includes all the residues of the R1 subunit commonly considered in the RNR substrate conversion scheme: Cys439 initiating the reaction as a thiyl radical, the redox-active cysteines Cys225 and Cys462, and the catalytically important Glu441 and Asn437. In contrast to previous theoretical studies of the overall mechanism, Glu441 is added as an anion. All relevant transition states have been optimized, including one where an electron is transferred 8 A from the disulfide to the substrate simultaneously with a proton transfer from Glu441. The calculated barrier for this step is 19.1 kcal/mol, which can be compared to the rate-limiting barrier indicated by experiments of about 17 kcal/mol. Even though the calculated barrier is somewhat higher than the experimental limit, the discrepancy is within the normal error bounds of B3LYP. The suggestion from the present modeling study is thus that a protonated Glu441 does not need to be present at the active site from the beginning of the catalytic cycle. However, the previously suggested mechanism with an initial protonation of Glu441 cannot be ruled out, because even with the cost added for protonation of Glu441 with a typical pK(a) of 4, the barrier for that mechanism is lower than the one obtained for the present mechanism. The results are compared to experimental results and suggestions.     

Computational Characterization of Reactive Intermediates of Carbon Monoxide Dehydrogenase

The catalytic oxidation of CO to CO2 by carbon monoxide dehydrogenases has been explored theoretically, and a large C-cluster model including the metal core [Ni-4Fe-4S] and surrounding residues and crystal water molecules was used in density functional calculations. The key species involved in the oxidation of CO at the C-cluster, Cred1, Cred2 and Cint, have been elucidated. On the basis of computational results, the plausible enzymatic mechanism for the CO oxidation was proposed. In the catalytic reaction, the first proton abstraction from the Fe（1）-bound water leads to a precursor to accommodate CO binding and the subsequently consecutive proton transfers from the metal-bound carboxylate to the amino acid residues facilitate the release of CO2. The hydrogen-bond network around the C-cluster formed by conserved residues His93, His96, Glu299, Lys563, and four water molecules in the active domain plays an important role in proton transfer and intermediate stabilization. Predicted geometries of key species show good agreement with the reported crystal structures.     

A QM/MM study of the racemization of vinylglycolate catalyzed by mandelate racemase enzyme.

The experimentally postulated mechanism for the interconversion between (S)-vinylglycolate and (R)-vinylglycolate catalyzed by mandelate racemase enzyme consists of a two-step quite symmetric process through a dianionic enolic intermediate that is formed after the abstraction of the alpha-proton of vinylglycolate by a basic enzymatic residue and is then reprotonated by another residue. The challenging problem behind this reaction is how the enzyme manages to stabilize such an intermediate, that is, how it lowers enough the high pK(a) of the alpha-proton for the reaction to take place. The QM/MM simulations performed in this paper indicate that catalysis is based on the stabilization of the negative charge developed on the substrate along the reaction. We have identified three different reaction mechanisms starting from different quasi-degenerate structures of the substrate-enzyme complex. In two of them the stabilizing role is done by means of a catalytic proton transfer that avoids the formation of a dianionic intermediate, and they involve six steps instead of the two experimentally proposed. On the contrary, the third mechanism passes through a dianionic species stabilized by the concerted approach of a protonated enzymatic residue during the proton abstraction. The potential energy barriers theoretically found along these mechanisms are qualitatively in good agreement with the experimental free energy barriers determined for racemization of vinylglycolate and mandelate. The theoretical study of the effect of the mutation of Glu317 by Gln317 in the kinetics of the reaction reveals the important role in the catalysis of the hydrogen bond formed by Glu317 in the native enzyme, as only one of the mechanisms, the slower one, is able to produce the racemization in the active site of the mutant. However, we have found that this hydrogen bond is not an LBHB within our model.     

Mechanism of RNase T1: concerted triester-like phosphoryl transfer via a catalytic three-centered hydrogen bond

BACKGROUND: The microscopic events of ribonuclease (RNase) catalyzed phosphoryl transfer reactions are still a matter of debate in which the contenders adhere to either the classical concerted acid-base mechanism or a more sequential triester-like mechanism. In the case of RNase A, small thio-effects of the nonbridging oxygens have been invoked in favor of the classical mechanism. However, the RNase T1 catalyzed transphosphorylation of phosphorothioate RNA is highly stereoselective. R(P) thio-substituted RNA is depolymerized 60000 times faster than S(P) thio-substituted RNA by this enzyme, whereas the uncatalyzed cleavage of both substrates occurs at comparable rates. We combined site-directed mutagenesis in the RNase active site and stereospecific thio-substitution of an RNA substrate to probe the intermolecular interactions of the enzyme with the nonbridging pro-S(P) oxygen that bring about this stereoselectivity of RNase T1. RESULTS: Thio-substitution of the nonbridging pro-S(P) oxygen in the substrate afflicts chemical turnover but not ground state binding whereas thio-substitution of the nonbridging pro-R(P) oxygen does not affect the kinetics of RNase T1. Site-directed mutagenesis of the catalytic base Glu58 impairs the enzyme's ability to discriminate both phosphorothioate diastereomers. Glu58Ala RNase T1 cleaves R(P) and S(P) phosphorothioate RNA with similar rates. The dependence of the pro-S(P) thio-effect on the presence of the Glu58 carboxylate evidences a strong rate-limiting interaction between the nonbridging pro-S(P) oxygen and the catalytic base Glu58 in the wild type enzyme. CONCLUSIONS: Based on these results, we put forward a new triester-like mechanism for the RNase T1 catalyzed reaction that involves a three-centered hydrogen bond between the 2'-OH group, the nonbridging pro-S(P) oxygen and one of the carboxylate oxygens of Glu58. This interaction allows nucleophilic attack on an activated phosphate to occur simultaneously with general base catalysis, ensuring concerted phosphoryl transfer via a triester-like mechanism.     

Theoretical studies on pyridoxal 5'-phosphate-dependent transamination of alpha-amino acids

Density functional methods have been applied to investigate the irreversible transamination between glyoxylic acid and pyridoxamine analog and the catalytic mechanism for the critical [1,3] proton transfer step in aspartate aminotransferase (AATase). The results indicate that the catalytic effect of pyridoxal 5'-phosphate (PLP) may be attributed to its ability to stabilize related transition states through structural resonance. Additionally, the PLP hydroxyl group and the carboxylic group of the amino acid can shuttle proton, thereby lowering the barrier. The rate-limiting step is the tautomeric conversion of the aldimine to ketimine by [1,3] proton transfer, with a barrier of 36.3 kcal/mol in water solvent. A quantum chemical model consisting 142 atoms was constructed based on the crystal structure of the native AATase complex with the product L-glutamate. The electron-withdrawing stabilization by various residues, involving Arg386, Tyr225, Asp222, Asn194, and peptide backbone, enhances the carbon acidity of 4'-C of PLP and Calpha of amino acid. The calculations support the proposed proton transfer mechanism in which Lys258 acts as a base to shuttle a proton from the 4'-C of PLP to Calpha of amino acid. The first step (proton transfer from 4'-C to lysine) is shown to be the rate-limiting step. Furthermore, we provided an explanation for the reversibility and specificity of the transamination in AATase.     

Theoretical investigation on the dissociation of (R)‐benzoin catalyzed by benzaldehyde lyase

Benzaldehyde lyase (BAL) is a versatile thiamin diphosphate (THDP)‐dependent enzyme with widespread synthetic applications in industry. Besides lyase activity, BAL also performs the functions as carboligase and decarboxylase. Unlike many other THDP‐dependent enzymes, the active center of BAL is devoid of any acid‐base amino acid residues except Glu50 and His29, and therefore, the catalytic mechanism of BAL is unusual. In this article, the dissociation mechanism of (R)‐benzoin to benzaldehyde catalyzed by BAL has been studied by using density functional theory method. The calculation results indicate that the whole reaction consists of four elementary steps, and at least two steps contribute to rate‐limiting. A big difference with other THDP‐dependent enzymes is that, in the first stage of the reaction, the ligation of substrate and THDP ylide is not companied by proton transfer, and in the subsequent transition states and intermediates, the carbonyl oxygen always exists in the form of anion. Gln113, His29, and 4′‐amino group of THDP are found to have the function to stabilize the transition states and intermediates. His29 acts as the proton acceptor in step 2 and proton donor in step 3 using one water molecule as mediator. © 2013 Wiley Periodicals, Inc.