Computational evidence for the catalytic mechanism of glutaminyl cyclase. A DFT investigation

The results of a DFT theoretical investigation on the catalytic mechanism of the QC enzyme are presented. A rather large model-system is used. It includes the most important residues that are believed to play a key-role in the catalysis. The computational results show that the rate-determining step of the catalytic process is not the nucleophilic attack leading to the cycle formation (a very easy and fast process with a negligible barrier of 0.8 kcal mol(-1)), but a proton transfer, which is assisted by the Glu201 residue acting as a proton shuttle (general base and general acid). A complex network of hydrogen bonds (involving Asp248 and other residues) contribute to lower the activation barrier for the proton shift which affords the formation of an ammonia molecule bonded to the substrate. The ammonia molecule is a good leaving group which is easily expelled from the substrate in the last step of the catalytic cycle, but remains anchored to the enzyme as a ligand of the zinc cation. The metal plays a key-role in assisting the nucleophilic attack (electrostatic catalysis) since it polarizes the substrate gamma-amide carbonyl group (its electrophilic character increases). Also, the strength of the nucleophilic nitrogen (substrate alpha-amino group) is enhanced by hydrogen bonds involving the Glu201 residue. The computations outline the important role of Trp329 in helping the substrate binding process and stabilizing the cyclization transition state.     

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.     

Catalytic role of the metal ion of carboxypeptidase A in ester hydrolysis.

The mechanism of action of bovine pancreatic carboxypeptidase. Aalpha (peptidyl-L-amino acid hydrolase; EC 3.4.12.2) has been investigated by application of cryoenzymologic methods. Kinetic studies of the hydrolysis of the specific ester substrate O-(trans-p-chlorocinnamoyl)-L-beta-phenyllactate have been carried out with both the native and the Co2+-substituted enzyme in the 25 to --45 degrees C temperature range. In the --25 to --45 degrees C temperature range with enzyme in excess, a biphasic reaction is observed for substrate hydrolysis characterized by rate constants for the fast (kf) and the slow (ks) processes. In Arrhenius plots, ks extrapolates to kcat at 25 degrees C for both enzymes in aqueous solution, indicating that the same catalytic rate-limiting step is observed. The slow process is analyzed for both metal enzymes, as previously reported (Makinen, M. W., Yamamura, K., and Kaiser, E. T. (1976) Proc Natl. Acad. Sci. U. S. A. 73, 3882-3886), to involve the deacylation of a mixed anhydride acyl-enzyme intermediate. Near --60 degrees C the acyl-enzyme intermediate of both metal enzymes can be stabilized for spectral characterization. The pH and temperature dependence of ks reveals a catalytic ionizing group with a metal ion-dependent shift in pKa and an enthalpy of ionization of 7.2 kcal/mol for the native enzyme and 6.2 kcal/mol for the Co2+ enzyme. These parameters identify the ionizing catalytic group as the metal-bound water molecule. Extrapolation of the pKa data to 25 degrees C indicates that this ionization coincides with that observed in the acidic limb of the pH profile of log(kcat/Km(app)) for substrate hydrolysis under steady state conditions. The results indicate that in the esterolytic reaction of carboxypeptidase. A deacylation of the mixed anhydride intermediate is catalyzed by a metal-bound hydroxide group.     

Oxygen exchange between acetate and the catalytic glutamate residue in glutaconate CoA-transferase from Acidaminococcus fermentans. Implications for the mechanism of CoA-ester hydrolysis.

The exchange of oxygen atoms between acetate, glutaryl-CoA, and the catalytic glutamate residue in glutaconate CoA-transferase from Acidaminococcus fermentans was analyzed using [(18)O(2)]acetate together with matrix-assisted laser desorption/ionization time of flight mass spectrometry of an appropriate undecapeptide. The exchange reaction was shown to be site-specific, reversible, and required both glutaryl-CoA and [(18)O(2)]acetate. The observed exchange is in agreement with the formation of a mixed anhydride intermediate between the enzyme and acetate. In contrast, with a mutant enzyme, which was converted to a thiol ester hydrolyase by replacement of the catalytic glutamate residue by aspartate, no (18)O uptake from H(2)(18)O into the carboxylate was detectable. This result is in accord with a mechanism in which the carboxylate of aspartate acts as a general base in activating a water molecule for hydrolysis of the thiol ester intermediate. This mechanism is further supported by the finding of a significant hydrolyase activity of the wild-type enzyme using acetyl-CoA as substrate, whereas glutaryl-CoA is not hydrolyzed. The small acetate molecule in the substrate binding pocket may activate a water molecule for hydrolysis of the nearby enzyme-CoA thiol ester.     

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).     

Studies on the mechanism of action of xanthine oxidase

Recent studies of the reaction mechanism of the molybdenum-containing enzyme xanthine oxidase are presented. The pH-dependence of both the steady-state and rapid reaction kinetics of the enzyme exhibits is bell-shaped, with pK(a)s for the acid and alkaline limbs of 6.6 and 7.4, respectively. These are assigned to ionizations of an active site base and substrate, respectively, with the implication that enzyme acts on the neutral rather than monoanionic form of the purine substrate. A computational study provides evidence that in the course of the reaction tautomerization of substrate occurs, with a proton moving from N-3 to N-9 in the course of the reaction - enzyme facilitation of this tautomerization may contribute as much as 24 kcal/mol in transition state stabilization for the reaction. Electron spin echo (ESEEM) and electron-nuclear double resonance (ENDOR) studies of the so-called "very rapid" Mo(V) intermediate of the reaction, the latter work using a newly synthesized form of the substrate 2-hydroxy-6-methylpurine that has been selectively isotopically labeled at C-8, indicates that product is bound to the molybdenum of the active site in a simple, end-on fashion, consistent with a reaction mechanism involving nucleophilic attack of a (deprotonated) Mo-OH on the C-8 position of substrate. A kinetic study using a series of purines has failed to identify a correlation between the one-electron reduction potential for substrate and catalytic effectiveness, indicating that a reaction mechanism initiated by one-electron, outer-sphere electron transfer is unlikely. Finally, a consideration of the active site structure in the context of the above work suggests specific amino acid residues to target for site-directed mutagenesis studies. Preliminary experiments with two such mutants are entirely consistent with the proposed catalytic roles of two active site glutamate residues.     

A comparison of dogfish and bovine chymotrypsins

Bovine and dogfish chymotrypsins were compared to determine if chymotrypsin from a poikilothermic organism (spiny dogfish (Squalus acanthias] adapted to low temperatures possessed catalytic properties different from those of the same enzyme from a warm-blooded animal. An improved procedure was developed for isolating dogfish pancreatic chymotrypsin. The least hydrophobic and smallest substrate used, p-nitrophenyl acetate, had similar enthalpies of association (delta Ha) with both enzymes, whereas larger, more hydrophobic substrates had delta Ha values that were of opposite sign for the two enzymes. As the temperature increased, the association constants (1/Ks) for p-nitrophenyl valerate and p-nitrophenyltrimethyl acetate increased for dogfish chymotrypsin and decreased for bovine chymotrypsin, while the free energies of association (delta Ga) remained relatively constant. Acylation of chymotrypsin was 1.5-2.5 times slower in the dogfish enzyme than in the bovine enzyme except below 15 degrees C with p-nitrophenyltrimethyl acetate. delta H++ for acylation by p-nitrophenyltrimethyl acetate were 2.0 kcal/mol for the dogfish enzyme and 10.2 kcal/mol for the bovine, whereas delta H++ values were only slightly lower in the dogfish enzyme for the other two substrates. For all substrates, the deacylation rate constant (kcat) was greater with dogfish chymotrypsin than bovine. However, the free energies of activation (delta G++) for deacylation were nearly equal between the two enzymes for each of the substrates.     

Potential of mean force calculation for the proton and hydride transfer reactions catalyzed by medium-chain acyl-CoA dehydrogenase: effect of mutations on enzyme catalysis

Potential of mean force calculations have been performed on the wild-type medium-chain acyl-CoA dehydrogenase (MCAD) and two of its mutant forms. Initial simulation and analysis of the active site of the enzyme reveal that an arginine residue (Arg256), conserved in the substrate-binding domain of this group of enzymes, exists in two alternate conformations, only one of which makes the enzyme active. This active conformation was used in subsequent computations of the enzymatic reactions. It is known that the catalytic alpha,beta-dehydrogenation of fatty acyl-CoAs consists of two C-H bond dissociation processes: a proton abstraction and a hydride transfer. Energy profiles of the two reaction steps in the wild-type MCAD demonstrate that the reaction proceeds by a stepwise mechanism with a transient species. The activation barriers of the two steps differ by only approximately 2 kcal/mol, indicating that both may contribute to the rate-limiting process. Thus this may be a stepwise dissociation mechanism whose relative barriers can be tuned by suitable alterations of the substrate and/or enzyme. Analysis of the structures along the reaction path reveals that Arg256 plays a key role in maintaining the reaction center hydrogen-bonding network involving the thioester carbonyl group, which stabilizes transition states as well as the intervening transient species. Mutation of this arginine residue to glutamine increases the activation barrier of the hydride transfer reaction by approximately 5 kcal/mol, and the present simulations predict a substantial loss of catalytic activity for this mutant. Structural analysis of this mutant reveals that the orientation of the thioester moiety of the substrate has been changed significantly as compared to that in the wild-type enzyme. In contrast, simulation of the active site of the Thr168Ala mutant shows no significant change in the relative orientation of the substrate and the cofactor in the active site; as a result, this mutation has very little effect on the overall reaction barrier, and this is consistent with the experimental data. This study demonstrates that significant insights into the catalytic mechanism can be obtained from simulation studies, and the results can be used to design novel mechanistic probes for the enzyme.     

A quantum mechanics/molecular mechanics study of the catalytic mechanism and product specificity of viral histone lysine methyltransferase

There are three reaction steps in the S-adenosylmethionine (AdoMet) methylation of lysine-NH2 catalyzed by a methyltransferase. They are (i) combination of enzyme.Lys-NH3+ with AdoMet, (ii) substrate ionization to provide enzyme.AdoMet.Lys-NH2, and (iii) methyl transfer providing enzyme.AdoHcy.Lys-N(Me)H2+ and the dissociation of AdoHcy. In this study of the viral histone methyltransferase (vSET), we find that substrate ionization of vSET.Lys27-NH3+, vSET.Lys27-N(Me)H2+, and vSET.Lys27-N(Me)2H+ takes place upon combination with AdoMet. The presence of a water channel allows dissociation of a proton to the solvent. There is no water channel in the absence of AdoMet. That the formation of a water channel is combined with AdoMet binding was first discovered in our investigation of Rubisco large subunit methyltransferase. Via a quantum mechanics/molecular mechanics (QM/MM) approach, the calculated free energy barrier (DeltaG++) of the first methyl transfer reaction catalyzed by vSET [Lys27-NH2 + AdoMet --> Lys27-N(Me)H2+ + AdoHcy] equals 22.5 +/- 4.3 kcal/mol, which is in excellent agreement with the free energy barrier (21.7 kcal/mol) calculated from the experimental rate constant (0.047 min-1). The calculated DeltaG++ of the second methyl transfer reaction [AdoMet + Lys27-N(Me)H --> AdoHcy + Lys27-N(Me)2H+] at the QM/MM level is 22.6 +/- 3.6 kcal/mol, which is in agreement with the value of 22.4 kcal/mol determined from the experimental rate constant (0.015 min-1). The third methylation [Lys27-N(Me)2 + AdoMet --> Lys27-N(Me)3+ + AdoHcy] is associated with a DeltaG++ of 23.1 +/- 4.0 kcal/mol, which is in agreement with the value of 23.0 kcal/mol determined from the experimental rate constant (0.005 min-1). Our computations establish that the first, second, and third methyl transfer steps catalyzed by vSET are linear SN2 reactions with the bond making being approximately 50% associative.     

A computational study of the deacylation mechanism of human butyrylcholinesterase

To investigate the mechanism of the deacylation reaction in the active site of human butyrylcholinesterase (BuChE), we carried out quantum mechanical (QM) calculations on cluster models of the active site built from a crystallographic structure. The models consisted of the substrate butyrate moiety, the catalytic triad of residues (Ser198, Glu325, and His438), the "oxy-anion hole" (Gly116, Gly117, and Ala199), the side chain of Glu197, four water molecules, the side chain of Ser225, and the peptide linkage between Val321 and Asn322. Analyses of the equilibrium geometries, electronic properties, and energies of the QM models gave insights into the catalytic mechanism. In addition, the QM calculations provided the data required to build a molecular mechanics representation of the reactive BuChE region that was employed in molecular dynamics simulations followed by molecular-mechanics-Poisson-Boltzmann (MM-PB) calculations. Subsequently, we combined the QM energies with average MM-PB energies to estimate the free energy of the reactive structures in the enzyme. The rate-determining step corresponds to the formation of a tetrahedral intermediate with a free-energy barrier of approximately 14.0 kcal/mol. The modulation of the BuChE activity, exerted by either neutral molecules (glycerol, GOL) or a second butyrylcholine (CHO) molecule bound to the cation-pi site, does not involve any significant allosteric effect. Interestingly, the presence of GOL or CHO stabilizes a product complex formed between a butyric acid molecule and BuChE. These results are in consonance with the crystallographic structure of BuChE, in which the catalytic Ser198 interacts with a butyric fragment, while the cation-pi site is occupied by one GOL molecule.     

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

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

QM/MM investigation of the catalytic mechanism of angiotensin-converting enzyme

Angiotensin-converting enzyme (ACE) converts angiotensin I to angiotensin II and degrades bradykinin and other vasoactive peptides. ACE inhibitors are used to treat diseases such as hypertension and heart failure. It is thus highly desirable to understand the catalytic mechanism of ACE, as this should facilitate the design of more powerful and selective ACE inhibitors. ACE exhibits two different active domains, the C-domain and the N-domain. In this work, we systematically investigated the inhibitor- and substrate-binding patterns in the N-domain of human ACE using a combined quantum mechanical and molecular mechanical approach. The hydrolysis of hippuryl–histidyl–leucine (HHL) as catalyzed by the N-domain of human somatic ACE was explored, and the effects of chloride ion on the overall reaction were also investigated. Two models, one with and one without a chloride ion at the first binding position, were then designed to examine the chloride dependence of inhibitor–substrate binding and the catalytic mechanism. Our calculations indicate that the hydrolysis reaction follows a stepwise general base/general acid catalysis path. The estimated mean free energy barrier height in the two models is about 15.6 kcal/mol, which agrees very well with the experimentally estimated value of 15.8 kcal/mol. Our simulations thus suggest that the N-domain is in a mixed form during ACE-catalyzed hydrolysis, with the single-chloride-ion and the double-chloride-ion forms existing simultaneously.

Open image in new window

Graphical Abstract Superposition of ACE C- and N- domains

     

Quantum mechanics/molecular mechanics modeling of substrate-assisted catalysis in family 18 chitinases: conformational changes and the role of Asp142 in catalysis in ChiB

Family 18 chitinases catalyze the hydrolysis of β-1,4-glycosidic bonds in chitin. The mechanism has been proposed to involve the formation of an oxazolinium ion intermediate via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile (instead of an enzyme residue). Here, we have modeled the first step of the chitin hydrolysis catalyzed by Serratia marcescens chitinase B for the first time using a combined quantum mechanics/molecular mechanics approach. The calculated reaction barriers based on multiple snapshots are 15.8-19.8 kcal mol(-1) [B3LYP/6-31+G(d)//AM1-CHARMM22], in good agreement with the activation free energy of 16.1 kcal mol(-1) derived from experiment. The enzyme significantly stabilizes the oxazolinium intermediate. Two stable conformations ((4)C(1)-chair and B(3,O)-boat) of the oxazolinium ion intermediate in subsite -1 were unexpectedly observed. The transition state structure has significant oxacarbenium ion-like character. The glycosyl residue in subsite -1 was found to follow a complex conformational pathway during the reaction ((1,4)B → [(4)H(5)/(4)E](++) → (4)C(1) ? B(3,O)), indicating complex conformational behavior in glycoside hydrolases that utilize a substrate-assisted catalytic mechanism. The D142N mutant is found to follow the same wild-type-like mechanism: the calculated barriers for reaction in this mutant (16.0-21.1 kcal mol(-1)) are higher than in the wild type, in agreement with the experiment. Asp142 is found to be important in transition state and intermediate stabilization.     

RNA catalytic properties of the minimum (-)sTRSV sequence

We have identified an RNA catalytic domain within the sequence of the 359 base long negative-strand satellite RNA of tobacco ringspot virus. The catalytic domain contains two minimal sequences of satellite RNA, a 50-base catalytic RNA sequence, and a 14-base substrate RNA sequence. The catalytic complex of catalytic RNA/substrate RNA represents a structure not previously found in any RNA catalytic reaction described to date. The reaction is truly catalytic since the catalytic RNA has multiple substrate cleavage events and is not consumed during the course of the reaction. A linear relationship is seen between reaction rate and catalytic RNA concentration. The reaction has a Km of 0.03 microM, a kcat of 2.1/min, a temperature optimum of near 37 degrees C, and an energy of activation of 19 kcal/mol.     

Comparative studies of the catalytic mechanisms of two chorismatases: CH‐fkbo and CH‐Hyg5

Chorismatase is an important enzyme involved in Shikimate pathway, which catalyzes the conversion of chorismate into pyruvate and (dihydro)‐benzoic acid derivatives. According to the outcomes of catalytic reactions, chorismatases can be divided into three subfamilies: CH‐Fkbo, CH‐Hyg5 and CH‐XanB2. Recently, the crystal structures of CH‐Fkbo and CH‐Hyg5 from Streptomyces hygroscopicus have been successfully obtained, allowing us to perform QM/MM calculations to explore the reaction details. Our calculation results support the proposal that CH‐Fkbo and CH‐Hyg5 employ different catalytic mechanisms and gave the mechanistic details. Fkbo follows a typical hydrolytic mechanism, which contains three consecutive steps, including the protonation step of the methylene group of substrate, the nucleophilic attack of the resulted carbocation by activated water and cleavage of C2′‐O8 bond of tetrahedral intermediate (hemiketal). The protonation of methylene group and the C2′‐O8 cleavage correspond to similar energy barriers (26.5 and 24.8 kcal/mol), suggesting both steps to be rate‐limiting. Whereas Hyg5 employs an intramolecular mechanism, in which the oxygen from C4 migrates to C3 via an arene oxide intermediate. The first step of Hyg5, which corresponds to the concerted protonation of methylene group and the cleavage of C3‐O8, is calculated to be rate‐limiting with an energy barrier of 26.3 kcal/mol. The nonconserved active site residue G240^Hyg5 (or A244^Fkb°) is suggested to be responsible for leading to different reaction mechanism in CH‐Fkbo and CH‐Hyg5. During the catalytic reaction, residue C327 plays an important role in directing the product selectivity in Hyg5 enzyme. Proteins 2017; 85:1146–1158. © 2017 Wiley Periodicals, Inc.     

Quantum chemical studies of the catalytic mechanism of N-terminal nucleophile hydrolase

Modeling of the catalytic mechanism of penicillin acylase, a member of the N-terminal nucleophile hydrolase superfamily, is for the first time conducted at ab initio quantum chemistry level. The uniqueness of this family of enzymes is that their active site lacks His and Asp (Glu) residues, comprising together with a serine residue the classical catalytic triad. The current investigation confirms that the amino group of the N-terminal serine residue in N-terminal hydrolases is capable of activating its own hydroxyl group. Using the MP2/RHF method with the 6−31+G** basis set, stationary points on the potential energy surface of the considered molecular system were located, corresponding to local minima (complexes of reagents, products, intermediate) and to saddle points (transition states). It turned out that the stage of acyl-serine formation proceeds via two transition states; the first one, which separates reagents from the so-called tetrahedral intermediate, has the highest relative energy (30 kcal/mol). In contrast to recently proposed empiric suggestions, we have found that participation of a bridging water molecule in proton shuttling is not necessary for the catalysis. The quantum chemical calculations showed a crucial role of a specific solvation in decreasing the activation barrier of the reaction by approximately 10 kcal/mol. Published in Russion in Biokhimiya, 2007, Vol. 72, No. 5, pp. 615–621.     

Crystal structure of a type II dihydrofolate reductase catalytic ternary complex

Type II dihydrofolate reductase (DHFR) is a plasmid-encoded enzyme that confers resistance to bacterial DHFR-targeted antifolate drugs. It forms a symmetric homotetramer with a central pore which functions as the active site. Its unusual structure, which results in a promiscuous binding surface that accommodates either the dihydrofolate (DHF) substrate or the NADPH cofactor, has constituted a significant limitation to efforts to understand its substrate specificity and reaction mechanism. We describe here the first structure of a ternary R67 DHFR.DHF.NADP+ catalytic complex, resolved to 1.26 A. This structure provides the first clear picture of how this enzyme, which lacks the active site carboxyl residue that is ubiquitous in Type I DHFRs, is able to function. In the catalytic complex, the polar backbone atoms of two symmetry-related I68 residues provide recognition motifs that interact with the carboxamide on the nicotinamide ring, and the N3-O4 amide function on the pteridine ring. This set of interactions orients the aromatic rings of substrate and cofactor in a relative endo geometry in which the reactive centers are held in close proximity. Additionally, a central, hydrogen-bonded network consisting of two pairs of Y69-Q67-Q67'-Y69' residues provides an unusually tight interface, which appears to serve as a "molecular clamp" holding the substrates in place in an orientation conducive to hydride transfer. In addition to providing the first clear insight regarding how this extremely unusual enzyme is able to function, the structure of the ternary complex provides general insights into how a mutationally challenged enzyme, i.e., an enzyme whose evolution is restricted to four-residues-at-a-time active site mutations, overcomes this fundamental limitation.     

Structure and possible catalytic residues of Taka-amylase A

A complete molecular model of Taka-amylase A consisting of 478 amino acid residues was built with the aid of amino acid sequence data. Some typical structural features of the molecule are described. A model fitting of an amylose chain in the catalytic site of the enzyme showed a possible productive binding mode between substrate and enzyme. On the basis of the difference Fourier analysis and the model fitting study, glutamic acid (Glu230) and aspartic acid (Asp297), which are located at the bottom of the cleft, were concluded to be the catalytic residues, serving as the general acid and base, respectively.     

Temperature dependence of binding and catalysis for the Cdc25B phosphatase

Using a combination of steady-state and single-turnover kinetics, we probe the temperature dependence of substrate association and chemistry for the reaction of Cdc25B phosphatase with its Cdk2-pTpY/CycA protein substrate. The transition state for substrate association is dominated by an enthalpic barrier (DeltaH(++) of 13 kcal/mol) and has a favorable entropic contribution of 4 kcal/mol at 298 K. Phosphate transfer from Cdk2-pTpY/CycA to enzyme (DeltaH(++) of 12 kcal/mol) is enthalpically more favorable than for the small molecule substrate p-nitrophenyl phosphate (DeltaH(++) of 18 kcal/mol), yet entropically less favorable (TDeltaS(++) of 2 vs. -6 kcal/mol at 298 K, respectively). By measuring the temperature dependence of binding and catalysis for several hotspot mutants involved in binding of protein substrate, we determine the enthalpy-entropy compensations for changes in rates of association and phosphate transfer compared to the wild type system. We conclude that the transition state for enzyme-substrate association involves tight and specific contacts at the remote docking site and that phospho-transfer from Cdk2-pTpY/CycA to the pre-organized active site of the enzyme is accompanied by unfavorable entropic rearrangements that promote rapid product dissociation.     

Theoretical insights into the protonation states of active site cysteine and citrullination mechanism of Porphyromonas gingivalis peptidylarginine deiminase

Porphyromonas gingivalis peptidylarginine deiminase (PPAD) catalyzes the citrullination of peptidylarginine, which plays a critical role in the rheumatoid arthritis (RA) and gene regulation. For a better understanding of citrullination mechanism of PPAD, it is required to establish the protonation states of active site cysteine, which is still a controversial issue for the members of guanidino‐group‐modifying enzyme superfamily. In this work, we first explored the transformation between the two states: State N (both C351 and H236 are neutral) and State I (both residues exist as a thiolate–imidazolium ion pair), and then investigated the citrullination reaction of peptidylarginine, using a combined QM/MM approach. State N is calculated to be more stable than State I by 8.46 kcal/mol, and State N can transform to State I via two steps of substrate‐assisted proton transfer. Citrullination of the peptidylarginine contains deamination and hydrolysis. Starting from State N, the deamination reaction corresponds to an energy barrier of 18.82 kcal/mol. The deprotonated C351 initiates the nucleophilic attack to the substrate, which is the key step for deamination reaction. The hydrolysis reaction contains two chemical steps. Both the deprotonated D238 and H236 can act as the bases to activate the hydrolytic water, which correspond to similar energy barriers (～17 kcal/mol). On the basis of our calculations, C351, D238, and H236 constitute a catalytic triad, and their protonation states are critical for both the deamination and hydrolysis processes. In view of the sequence similarity, these findings may be shared with human PAD1–PAD4 and other guanidino‐group‐modifying enzymes. Proteins 2017; 85:1518–1528. © 2017 Wiley Periodicals, Inc.