Computational studies on pterins and speculations on the mechanism of action of dihydrofolate reductase

The significance of the enol form of the pterin ring in enzymatic reduction of dihydrofolate by DHFR is discussed on the basis of the results of ab initio calculations carried out on the keto/enol tautomers of 6-methyl-7, 8-dihydropterin as the model compound for the natural substrate, dihydrofolate.     

Computational studies of the domain movement and the catalytic mechanism of thymidine phosphorylase

Thymidine phosphorylase (TP) is a dual substrate enzyme with two domains. Each domain binds a substrate. In the crystal structure of Escherichia coli TP, the two domains are arranged so that the two substrate binding sites are too far away for the two substrates to directly react. Molecular dynamics simulations reveal a different structure of the enzyme in which the two domains have moved to place the two substrates in close contact. This structure has a root-mean-square deviation from the crystal structure of 4.1 A. Quantum mechanical calculations using this structure find that the reaction can proceed by a direct nucleophilic attack with a low barrier. This mechanism is not feasible in the crystal structure environment and is consistent with the mechanism observed for other N-glycosidic enzymes. Important catalytic roles are found for the three highly conserved residues His 85, Arg 171, and Lys 190.     

Computational and experimental studies of the catalytic mechanism of Thermobifida fusca cellulase Cel6A (E2)

Mutagenesis experiments suggest that Asp79 in cellulase Cel6A (E2) from Thermobifida fusca has a catalytic role, in spite of the fact that this residue is more than 13 A from the scissile bond in models of the enzyme-substrate complex built upon the crystal structure of the protein. This suggests that there is a substantial conformational shift in the protein upon substrate binding. Molecular mechanics simulations were used to investigate possible alternate conformations of the protein bound to a tetrasaccharide substrate, primarily involving shifts of the loop containing Asp79, and to model the role of water in the active site complex for both the native conformation and alternative low-energy conformations. Several alternative conformations of reasonable energy have been identified, including one in which the overall energy of the enzyme-substrate complex in solution is lower than that of the conformation in the crystal structure. This conformation was found to be stable in molecular dynamics simulations with a cellotetraose substrate and water. In simulations of the substrate complexed with the native protein conformation, the sugar ring in the -1 binding site was observed to make a spontaneous transition from the (4)C(1) conformation to a twist-boat conformer, consistent with generally accepted glycosidase mechanisms. Also, from these simulations Tyr73 and Arg78 were found to have important roles in the active site. Based on the results of these various MD simulations, a new catalytic mechanism is proposed. Using this mechanism, predictions about the effects of changes in Arg78 were made which were confirmed by site-directed mutagenesis.     

Theoretical studies on the activation of the pterin cofactor in the catalytic mechanism of dihydrofolate reductase

Two mechanisms for facilitating hydride ion transfer from NADPH involving preprotonation of the pteridine rings of the dihydrofolate reductase substrates folate and dihydrofolate have been investigated by ab initio quantum mechanical methods. Protonation energies and effective solution pKas have been calculated for four protonated forms, three of which are nonpreferred in aqueous solution and therefore not directly accessible to experimental study. The pattern and degree of redistribution of the positive charge over the component rings of the N-heterobicyclic pi-system in these protonated forms have been analyzed in terms of changes in the electron populations of the ring atoms and total ring charges. The effects of such changes in promoting hydride ion transfer to C7 in folate and C6 in dihydrofolate have been evaluated by considering the extent of development of partial carbonium ion character at these carbon atoms and also the degree of electron deficiency in the pyrazine ring as a whole. The results illustrate that perturbations due, for instance, to protonation may be propagated by pi-electron coupling effects over medium-range distances of 4-6 A across the pteridine ring. The two mechanisms have been assessed in terms of the calculated absolute and relative pKas of the protonated species taking into account experimental information regarding possible stabilization of these forms in the enzyme active site and also the effectiveness of the various protonations in assisting the hydride ion transfer step. Judged against these criteria, the theoretical results favor the generally proposed mechanism involving preprotonation of N8 in folate and N5 in dihydrofolate. However, some support was also found for the alternative novel mechanism involving O4-protonation of both folate and dihydrofolate.     

The catalytic mechanism of yeast thiosulfate reductase

Thiosulfate reductase catalyzes the desulfuration of thiosulfonates while oxidizing GSH to GSSG. Kinetic studies of the enzyme-catalyzed reaction between GSH and benzenethiosulfonate have been carried out, and direct evidence for the occurrence of glutathione persulfide as an immediate product of the reaction has been obtained. The formal mechanism of this enzymic reaction has been shown to be rapid equilibrium-ordered with GSH as the leading substrate.     

Further studies on the catalytic mechanism of human liver alpha-L-fucosidase

Radiolabeling of human liver alpha-L-fucosidase (alpha-L-fucoside fucohydrolase, EC 3.2.1.51) with [1-3H]conduritol C trans-epoxide revealed that there are four active sites per tetrameric enzyme complex. Solvent isotope effect experiments give evidence for a proton transfer at the rate-limiting step in catalysis. Transglycosylase activity was observed using methanol as an alternative glycone acceptor to produce methyl alpha-L-fucoside, suggesting that alpha-L-fucose is formed when water is the acceptor. Initial burst kinetics experiments suggest that a glycosyl-enzyme intermediate is formed, although the magnitude of the burst is not stoichiometric with the number of active sites. These data, along with previous results, suggest a general acid-general base catalytic mechanism involving double inversion of stereochemistry at C-1 of fucose, as well as the formation of either a covalent glycosyl-enzyme intermediate or a tight ion pair between a charged active-site residue and a hypothetical fucosyl oxocarbonium ion intermediate.     

Open questions in ferredoxin-NADP+ reductase catalytic mechanism.

Ferredoxin (flavodoxin)-NADP(H) reductases (FNR) are ubiquitous flavoenzymes that deliver NADPH or low potential one-electron donors (ferredoxin, flavodoxin) to redox-based metabolisms in plastids, mitochondria and bacteria. The plant-type reductase is also the basic prototype for one of the major families of flavin-containing electron transferases that display common functional and structural properties. Many aspects of FNR biochemistry have been extensively characterized in recent years using a combination of site-directed mutagenesis, steady-state and transient kinetic experiments, spectroscopy and X-ray crystallography. Despite these considerable advances, various key features in the enzymology of these important reductases remain yet to be explained in molecular terms. This article reviews the current status of these open questions. Measurements of electron transfer rates and binding equilibria indicate that NADP(H) and ferredoxin interactions with FNR result in a reciprocal decrease of affinity, and that this induced-fit step is a mandatory requisite for catalytic turnover. However, the expected conformational movements are not apparent in the reported atomic structures of these flavoenzymes in the free state or in complex with their substrates. The overall reaction catalysed by FNR is freely reversible, but the pathways leading to NADP+ or ferredoxin reduction proceed through entirely different kinetic mechanisms. Also, the reductases isolated from various sources undergo inactivating denaturation on exposure to NADPH and other electron donors that reduce the FAD prosthetic group, a phenomenon that might have profound consequences for FNR function in vivo. The mechanisms underlying this reductive inhibition are so far unknown. Finally, we provide here a rationale to interpret FNR evolution in terms of catalytic efficiency. Using the formalism of the Albery-Knowles theory, we identified which parameter(s) have to be modified to make these reductases even more proficient under a variety of conditions, natural or artificial. Flavoenzymes with FNR activity catalyse a number of reactions with potential importance for biotechnological processes, so that modification of their catalytic competence is relevant on both scientific and technical grounds.     

Studies on the specificity of the tetrapyrrole substrate for human biliverdin-IXalpha reductase and biliverdin-IXbeta reductase. Structure-activity relationships define models for both active sites

A comparison of the initial rate kinetics for human biliverdin-IXalpha reductase and biliverdin-IXbeta reductase with a series of synthetic biliverdins with propionate side chains "moving" from a bridging position across the central methene bridge (alpha isomers) to a "gamma-configuration" reveals characteristic behavior that allows us to propose distinct models for the two active sites. For human biliverdin-IXalpha reductase, as previously discussed for the rat and ox enzymes, it appears that at least one "bridging propionate" is necessary for optimal binding and catalytic activity, whereas two are preferred. All other configurations studied were substrates for human biliverdin-IXalpha reductase, albeit poor ones. In the case of mesobiliverdin-XIIIalpha, extending the propionate side chains to hexanoate resulted in a significant loss of activity, whereas the butyrate derivative retained high activity. For human biliverdin-IXalpha reductase, we suggest that a pair of positively charged side chains play a key role in optimally binding the IXalpha isomers. In the case of human biliverdin-IXbeta reductase, the enzyme cannot tolerate even one propionate in the bridging position, suggesting that two negatively charged residues on the enzyme surface may preclude productive binding in this case. The flavin reductase activity of biliverdin-IXbeta reductase is potently inhibited by mesobiliverdin-XIIIalpha and protohemin, which is consistent with the hypothesis that the tetrapyrrole and flavin substrate bind at a common site.     

Rapid-scan stopped-flow studies of the flavoenzyme mercuric reductase during catalytic turnover

Time-resolved absorption spectra of the FAD-containing enzyme mercuric reductase were recorded during the catalytic reaction at 25 degrees C, pH 7.3. With an excess of NADPH over Hg2+ there was a rapid (k = 43 s-1) initial formation of a spectral species similar to that previously assigned to an NADPH complex of two-electron-reduced enzyme, EH2-NADPH. This spectrum persisted during the quasisteady-state phase of the reaction suggesting that EH2-NADPH is a true catalytic intermediate and that the rate of catalysis is limited by the oxidation of EH2-NADPH by Hg2+. Also with an excess of Hg2+ over NADPH a spectrum similar to that of EH2-NADPH was rapidly formed. As the NADPH was exhausted, the spectrum of oxidized enzyme, E, did not reappear but rather a spectrum similar to that previously assigned to an NADP+ complex of two-electron-reduced enzyme, EH2-NADP+. These results suggest that EH2-HADP+ cannot rapidly reduce the Hg2+ substrate. However, eventually all reducing equivalents from NADPH added to oxidized, activated enzyme are utilized for the reduction of Hg2+. A mechanism model is proposed that does not involve the free, oxidized enzyme in the catalytic cycle.     

The catalytic mechanism of Pseudomonas aeruginosa cd1 nitrite reductase

The cd1 NiRs (nitrite reductases) are enzymes catalysing the reduction of nitrite to NO (nitric oxide) in the bacterial energy conversion denitrification process. These enzymes contain two distinct redox centres: one covalently bound c-haem, which is reduced by external electron donors, and another peculiar porphyrin, the d1-haem (3,8-dioxo-17-acrylate-porphyrindione), where nitrite is reduced to NO. In the present paper, we summarize the most recent results on the mechanism of nitrite reduction by the cd1 NiR from Pseudomonas aeruginosa. We discuss the essential catalytic features of this enzyme, with special attention to the allosteric regulation of the enzyme's activity and to the mechanism employed to avoid product inhibition, i.e. trapping of the active-site reduced haem by the product NO. These results shed light on the reactivity of cd1 NiRs and assign a central role to the unique d1-haem, present only in this class of enzymes.     

Investigations of the catalytic mechanism of thioredoxin glutathione reductase from Schistosoma mansoni

Thioredoxin glutathione reductase from Schistosoma mansoni (SmTGR) catalyzes the reduction of both thioredoxin and glutathione disulfides (GSSG), thus playing a crucial role in maintaining redox homeostasis in the parasite. In line with this role, previous studies have demonstrated that SmTGR is a promising drug target for schistosomiasis. To aid in the development of efficacious drugs that target SmTGR, it is essential to understand the catalytic mechanism of SmTGR. SmTGR is a dimeric flavoprotein in the glutathione reductase family and has a head-to-tail arrangement of its monomers; each subunit has the components of both a thioredoxin reductase (TrxR) domain and a glutaredoxin (Grx) domain. However, the active site of the TrxR domain is composed of residues from both subunits: FAD and a redox-active Cys-154/Cys-159 pair from one subunit and a redox-active Cys-596'/Sec-597' pair from the other; the active site of the Grx domain contains a redox-active Cys-28/Cys-31 pair. Via its Cys-28/Cys-31 dithiol and/or its Cys-596'/Sec-597' thiol-selenolate, SmTGR can catalyze the reduction of a variety of substrates by NADPH. It is presumed that SmTGR catalyzes deglutathionylation reactions via the Cys-28/Cys-31 dithiol. Our anaerobic titration data suggest that reducing equivalents from NADPH can indeed reach the Cys-28/Cys-31 disulfide in the Grx domain to facilitate reductions effected by this cysteine pair. To clarify the specific chemical roles of each redox-active residue with respect to its various reactivities, we generated variants of SmTGR. Cys-28 variants had no Grx deglutathionylation activity, whereas Cys-31 variants retained partial Grx deglutathionylation activity, indicating that the Cys-28 thiolate is the nucleophile initiating deglutathionylation. Lags in the steady-state kinetics, found when wild-type SmTGR was incubated at high concentrations of GSSG, were not present in Grx variants, indicating that this cysteine pair is in some way responsible for the lags. A Sec-597 variant was still able to reduce a variety of substrates, albeit slowly, showing that selenocysteine is important but is not the sole determinant for the broad substrate tolerance of the enzyme. Our data show that Cys-520 and Cys-574 are not likely to be involved in the catalytic mechanism.     

QM/MM studies on the catalytic mechanism of Phenylethanolamine N-methyltransferase

Epinephrine is a naturally occurring adrenomedullary hormone that transduces environmental stressors into cardiovascular actions. As the only route in the catecholamine biosynthetic pathway, Phenylethanolamine N-methyltransferase (PNMT) catalyzes the synthesis of epinephrine. To elucidate the detailed mechanism of enzymatic catalysis of PNMT, combined quantum-mechanical/molecular-mechanical (QM/MM) calculations were performed. The calculation results reveal that this catalysis contains three elementary steps: the deprotonation of protonated norepinphrine, the methyl transferring step and deprotonation of the methylated norepinphrine. The methyl transferring step was proved to be the rate-determining step undergoing a SN2 mechanism with an energy barrier of 16.4kcal/mol. During the whole catalysis, two glutamic acids Glu185 and Glu219 were proved to be loaded with different effects according to the calculations results of the mutants. These calculation results can be used to explain the experimental observations and make a good complementarity for the previous QM study.     

QM/MM studies on the catalytic mechanism of Phenylethanolamine N-methyltransferase

Epinephrine is a naturally occurring adrenomedullary hormone that transduces environmental stressors into cardiovascular actions. As the only route in the catecholamine biosynthetic pathway, Phenylethanolamine N-methyltransferase (PNMT) catalyzes the synthesis of epinephrine. To elucidate the detailed mechanism of enzymatic catalysis of PNMT, combined quantum-mechanical/molecular-mechanical (QM/MM) calculations were performed. The calculation results reveal that this catalysis contains three elementary steps: the deprotonation of protonated norepinphrine, the methyl transferring step and deprotonation of the methylated norepinphrine. The methyl transferring step was proved to be the rate-determining step undergoing a SN2 mechanism with an energy barrier of 16.4 kcal/mol. During the whole catalysis, two glutamic acids Glu185 and Glu219 were proved to be loaded with different effects according to the calculations results of the mutants. These calculation results can be used to explain the experimental observations and make a good complementarity for the previous QM study.     

An essential histidine residue in the catalytic mechanism of mammalian glutathione reductase

Glutathione reductase from human erythrocytes was inactivated by ethoxyformic anhydride, and > 95% activity was lost by modification of about 1–1.5 histidine residues per flavin (or subunit), as measured by the increased absorbance at 240 nm. Full reactivation was obtained with hydroxylamine. The rate of inactivation increased with pH and an apparent pK = 5.9 was obtained for the protolytic dissociation. The modified enzyme was inactive with NADPH and GSSG as substrates, but almost fully active in catalysis of a transhydrogenase reaction involving pyridine nucleotides. The visible absorption spectrum of oxidized or two-electron-reduced enzyme was not changed, but the flavin fluorescence of oxidized enzyme increased 2-fold after the modification. NADPH or NADP⁺ did not protect the enzyme against inactivation. It is concluded that the modification affects a histidine involved in the second half-reaction of the catalysis, i.e. reduction of GSSG by the dithiol of reduced enzyme. Glutathione reductase from three additional mammalian sources was similarly inactivated, but enzyme from yeast was much less inactivated by the corresponding treatment with ethoxyformic anhydride.     

P-glycoprotein catalytic mechanism: studies of the ADP-vanadate inhibited state

Kinetics of inhibition of ATPase activity of pure mouse Mdr3 P-glycoprotein upon incubation with MgADP and vanadate were studied along with the trapping of [14C]ADP in presence of vanadate. The presence of verapamil strongly magnified both effects. Inhibition of ATPase was also increased by several other drugs known to bind to drug-binding sites. Inhibition by ADP-vanadate was slow and depended cooperatively on nucleotide binding. Stoichiometry of [14C]ADP trapping by vanadate was 1 mol/mol P-glycoprotein at full inhibition. Catalytic site mutants prevented [14C]ADP trapping, whereas interdomain signal communication mutants reduced it in approximate correlation with their effects upon drug stimulation of ATPase. In explanation of the results, we propose that a "closed conformation" involving dimerization and interdigitation of the two nucleotide-binding domains is necessary to allow inhibition by ADP-vanadate. The results suggest that such a conformation occurs naturally during ATP hydrolysis. It is proposed that in order for the catalytic transition state to form, the two nucleotide-binding domains dimerize to form an integrated single entity containing two bound ATP with just one of the two ATP being hydrolyzed per dimerization event.     

Structural, spectroscopic and catalytic activity studies on glutathione reductase reconstituted with FAD analogues

FAD-modified human glutathione reductases were reconstituted from apoenzyme using the FAD analogues 6-SH-FAD, 6-SCN-FAD, 6-OH-FAD, 6-NH2-FAD and 8-OH-FAD. The catalytic activities of the modified enzymes were substantially lower than for the native enzyme. All five species could be crystallized, but only those containing 6-SH-FAD, 6-OH-FAD and 6-NH2-FAD yielded crystals that could be analyzed. X-ray analyses and structural refinements were performed at 0.27 nm and 0.30 nm resolution resulting in R factors around 13.5%. The crystal structures showed the additional non-hydrogen atoms and small conformational changes of the polypeptide that were obviously induced by the substituents of the FAD analogues. The observed changes together with spectroscopic and activity data permit some conclusions about the chemical nature of the substituents.     

Kinetic characterization of the catalytic mechanism of methionine sulfoxide reductase B from Neisseria meningitidis

Methionine sulfoxide reductases catalyze the thioredoxin-dependent reduction of methionine sulfoxide back to methionine. The methionine sulfoxide reductases family is composed of two structurally unrelated classes of enzymes named MsrA and MsrB, which display opposite stereoselectivities toward the sulfoxide function. Both enzymes are monomeric and share a similar three-step chemical mechanism. First, in the reductase step, a sulfenic acid intermediate is formed with a concomitant release of 1 mol of methionine per mol of enzyme. Then, an intradisulfide bond is formed. Finally, Msrs return back to reduced forms via reduction by thioredoxin. In the present study, it is shown for the Neisseria meningitidis MsrB that (1) the reductase step is rate-determining in the process leading to formation of the disulfide bond and (2) the thioredoxin-recycling process is rate-limiting. Moreover, the data suggest that within the thioredoxin-recycling process, the rate-limiting step takes place after the two-electron chemical exchange and thus is associated with the release of oxidized thioredoxin.     

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.     

Complementary perturbation of the kinetic mechanism and catalytic effectiveness of dihydrofolate reductase by side-chain interchange.

The variable residue Leu-28 of Escherichia coli dihydrofolate reductase (DHFR) and the corresponding residue Phe-31 in murine DHFR were interchanged, and the impact on catalysis was evaluated by steady-state and pre-steady-state analysis. The E. coli L28F mutant increased the pH-independent kcat from 11 to 50 s-1 but had little effect on Km(H2F). An increase in the rate constant for dissociation of H4F from E.H4F.NH (from 12 to 80 s-1) was found to be largely responsible for the increase in kcat. Unexpectedly, the rate constant for hydride transfer increased from 950 to 4000 s-1 with little perturbation of NADPH and NADP+ binding to E. Consequently, the flux efficiency of the E. coli L28F mutant rose from 15% to 48% and suggests a role in genetic selection for this variable side chain. The murine F31L mutant decreased the pH-independent kcat from 28 to 4.8 s-1 but had little effect on Km(H2F). A decrease in the rate constant for dissociation of H4F from E.H4F.NH (from 40 to 22 s-1) and E.H4F (from 15 to 0.4 s-1) was found to be mainly responsible for the decrease in kcat. The rate constant for hydride transfer decreased from 9000 to 5000 s-1 with minor perturbation of NADPH binding. Thus, the free energy differences along the kinetic pathway were generally similar in magnitude but opposite in direction to those incurred by the E. coli L28F mutant. This conclusion implies that DHFR hydrophobic active-site side chains impart their characteristics individually and not collectively.     

Roles of N-terminal active cysteines and C-terminal cysteine-selenocysteine in the catalytic mechanism of mammalian thioredoxin reductase

Mammalian thioredoxin reductase [EC 1.6.4.5], a homodimeric flavoprotein, has a marked similarity to glutathione reductase. The two cysteines in the N-terminal FAD domain (-Cys59-x-x-x-x-Cys64-) and histidine (His472) are conserved between them at corresponding positions, but the mammalian thioredoxin reductase contains a C-terminal extension of selenocysteine (Sec or U) at the penultimate position and a preceding cysteine (-Gly-Cys497-Sec498-Gly). Introduction of mutations into the cloned rat thioredoxin reductase gene revealed that residues Cys59, Cys64, His472, Cys497, and Sec498, as well as the sequence of Cys497 and Sec498 were essential for thioredoxin-reducing activity. To analyze the catalytic mechanism of the mammalian thioredoxin reductase, the wild-type, U498C, U498S, C59S, and C64S were overproduced in a baculovirus/insect cell system and purified. The wild-type thioredoxin reductase produced in this system, designated as WT, was found to lack the Sec residue and to terminate at Cys497. A Sec-containing thioredoxin reductase, which was purified from COS-1 cells transfected with the wild-type cDNA, was designated as SecWT and was used as an authentic enzyme. Among mutant enzymes, only U498C retained a slight thioredoxin-reducing activity at about three orders magnitude lower than SecWT. WT, U498C, and U498S showed some 5,5'-dithiobis(2-nitrobenzoic acid)-reducing activity and transhydrogenase activity, and C59S and C64S had substantially no such activities. These data and spectral analyses of these enzymes suggest that Cys59 and Cys64 at the N-terminus, in conjunction with His472, function as primary acceptors for electrons from NADPH via FAD, and that the electrons are then transferred to Cys497-Sec498 at the C-terminus for the reduction of oxidized thioredoxin in the mammalian thioredoxin reductase.