Mechanism for ribonucleotide reductase inactivation by the anticancer drug gemcitabine

Gemcitabine (2',2'-difluoro-2'-deoxycytidine, dFdC) is a very promising anticancer drug, already approved for clinical use in three therapeutic indications. It is metabolized intracellularly to 5'-diphosphate (dFdCDP), which is known to be a potent inhibitor of ribonucleotide reductase (RNR). Although several nucleotide analogs show in vitro capacity of RNR inactivation, none has shown the in vivo efficacy of dFdCDP. Accordingly, the experimental data suggests that its mechanism of inhibition is different from the other known RNR suicide inhibitors. Enzyme inhibition in the absence of reductive species leads to complete loss of the essential radical in subunit R2, and formation of a new nucleotide-based radical. Interestingly, however, the presence of the reductants does not prevent inhibition--the radical is not lost but the targeted subunit of RNR becomes R1, which is inactivated possibly by alkylation. We have conducted a theoretical study, which led us to the first proposal of a possible mechanism for RNR inhibition by dFdCDP in the absence of reductants. This mechanism turned out to be very similar to the natural substrate reduction pathway and only deviates from the natural course after the formation of the well-known disulphide bridge. This deviation is caused precisely by the F atom in the beta-face, only present in this inhibitor. The essential radical in R2 is lost, and so is the enzyme catalytic activity. The nucleotide-based radical that constitutes the end product of our mechanism has been suggested in the literature as a possible candidate for the one detected experimentally. In fact, all experimental data available has been reproduced by the theoretical calculations performed here.     

Theoretical study on the inhibition of ribonucleotide reductase by 2'-mercapto-2'-deoxyribonucleoside-5'-diphosphates

Ribonucleotide reductase (RNR) is responsible for the reduction of ribonucleotides into the correspondent 2'-deoxyribonucleotides in the only physiological process that yields the monomers of DNA. The enzyme has thus become an attractive target for chemotherapies that fight proliferation-based diseases, specifically cancer and infections by some viruses and parasites. 2'-Mercapto-2'-deoxyribonucleoside-5'-diphosphates (SHdNDP) are mechanism-based inhibitors of RNR and therefore potential chemotherapeutic agents for those indications. Previous experimental studies established the in vitro and in vivo activity of SHdNDP. In the in vitro studies, it was observed that the activity was dependent on the oxidative status of the medium, with the inactivation of RNR only occurring when molecular oxygen was available. To better understand the mechanism involved in RNR inactivation by SHdNDP, we performed theoretical calculations on the possible reactions between the inhibitors and the RNR active site. As a result, we propose the possible mechanistic pathways for the chemical events that occur in the absence and in the presence of O2. They correspond to a refinement and a complement of those proposed in the literature.     

Theoretical studies on the mode of inhibition of ribonucleotide reductase by 2'-substituted substrate analogues

Several 2'-substituted-2'-deoxyribonucleotides are potent time-dependent inactivators of the enzyme ribonucleotide reductase (RNR), which function by destructing its essential tyrosil radical and/or by performing covalent addition to the enzyme. The former leads to inhibition of the R2 dimer of RNR and the latter to inhibition of the R1 dimer. Efforts to elucidate the mechanism of inhibition have been undertaken in the last decades, and a general mechanistic scheme has emerged. Accordingly, two alternative pathways lead either to the inhibition of R1 or R2, for which the 2'-chloro-2'-deoxynucleotides serve as the model for the inhibition of R1 and the 2'-azido-2'-deoxynucleotides the model for the inhibition of R2. However, the underlying reason for the different behavior of the inhibitors has remained unknown until now. Moreover, a fundamental mechanistic alternative has been proposed, based on results from biomimetic reactions, in which the 2'-substituents would be eliminated as radicals, and not as anions, as previously assumed. This would lead to further reactions not predicted by the existing mechanistic scheme. To gain a better understanding we have performed high-level theoretical calculations on the active site of RNR. Results from this work support the general Stubbe's paradigm, although some changes to that mechanism are necessary. In addition, a rational explanation of the factors that determine which of the dimers (R1 or R2) will be inactivated is provided for the first time. It has been demonstrated also that the 2'-substituents are indeed eliminated as anions, and not as radicals. Biomimetic experiments have led to different results because they lack a basic group capable of deprotonating the 3'-HO group of the substrate. It has been found here that the chemical character of the leaving group (radical or anionic) can be manipulated by controlling the protonation state of the 3'-HO group.     

Ribonucleotide activation by enzyme ribonucleotide reductase: understanding the role of the enzyme

This article focuses on the first step of the catalytic mechanism for the reduction of ribonucleotides catalyzed by the enzyme Ribonucleotide Reductase (RNR). This corresponds to the activation of the substrate. In this work a large model of the active site region involving 130 atoms was used instead of the minimal gas phase models used in previous works. The ONIOM method was employed to deal with such a large system. The results gave additional information, which previous small models could not provide, allowing a much clearer evaluation of the role of the enzyme in this step. Enzyme-substrate interaction energies, specific transition state stabilization, and substrate steric strain energies were obtained. It was concluded that the transition state is stabilized in 4.0 kcal/mol by specific enzyme-substrate interactions. However, this stabilization is cancelled by the cost in conformational energy for the enzyme to adopt the transition state geometry; the overall result is that the enzyme machinery does not lead to a rate enhancement in this step. It was also found that the substrate binds to the active site with almost no steric strain, emphasizing the complementarity and specificity of the RNR active site for nucleotide binding. The main role of the enzyme at the very beginning of the catalytic cycle was concluded to be to impose stereospecifity upon substrate activation and to protect the enzyme radical from the solvent, rather than to be an reaction rate enhancement.     

Understanding ribonucleotide reductase inactivation by gemcitabine

This paper focuses on the inhibition of ribonucleotide reductase (RNR) by gemcitabine, 2',2'-difluoro-2'-deoxycytidine (dFdC), a deoxycytidine analogue that is a very active drug against solid tumors and is currently commercialized as gemzar. RNR inactivation is reductant-dependent and occurs in a very different way from that of other known substrate analogues. In the presence of reductants monomer R1 of RNR is inhibited, whereas in the absence of reductants the radical is lost and monomer R2 is inhibited. As inside the cell reductants are available, it is likely that R1 inactivation is the most favorable mechanism responsible for drug cytotoxicity. This inhibition pathway has been unknown to date, but we have conducted a theoretical study that has led us to the first proposal of a mechanism for RNR inhibition by dFdC in the presence of reductants.     

Theoretical studies on the mechanism of inhibition of Ribonucleotide Reductase by (E)-2'-Fluoromethylene-2'-deoxycitidine-5'-diphosphate

(E)-2'-Fluoromethylene-2'-deoxycitidine-5'-diphosphate (FMCDP) is a potent time-dependent inactivator of the enzyme Ribonucleotide Reductase, which operates by destructing an essential tyrosil radical and performing a covalent addition to an active site residue. Considerable effort to elucidate the inhibition mechanism has been undertaken in recent years, and some insights have been obtained. Although a mechanistic proposal has been put forward, based on a general paradigm of inhibition of RNR by 2' substituted substrate analogues, details about the mechanism have remained elusive. Namely, the exact residue that adds to the inhibitor is still not identified, although mutagenesis experiments suggest that it should correspond to the E441 residue. In this work, we performed an extensive theoretical exploration of the potential energy surface of a model system representing the active site of RNR with FMCDP, using Density Functional Theory. This study establishes the detailed mechanism of inhibition, which is considerably different from the one proposed earlier. The proposed mechanism is fully consistent with available experimental data. Energetic results reveal unambiguously that the residue adding to the inhibitor is a cysteine, most probably C439, and exclude the possibility of the addition of E441. However, the E441 residue is shown to be essential for inhibition, catalyzing both the major and minor inhibition pathways, in agreement with mutagenesis results. It is shown also that the major mode of inactivation mimics the early stages of the natural substrate pathway.     

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.     

Rational polynomial representation of ribonucleotide reductase activity

Background  

Recent data suggest that ribonucleotide reductase (RNR) exists not only as a heterodimer R1₂R2₂ of R1₂ and R2₂ homodimers, but also as tetramers R1₄R2₄ and hexamers R1₆R2₆. Recent data also suggest that ATP binds the R1 subunit at a previously undescribed hexamerization site, in addition to its binding to previously described dimerization and tetramerization sites. Thus, the current view is that R1 has four NDP substrate binding possibilities, four dimerization site binding possibilities (dATP, ATP, dGTP, or dTTP), two tetramerization site binding possibilities (dATP or ATP), and one hexamerization site binding possibility (ATP), in addition to possibilities of unbound site states. This large number of internal R1 states implies an even larger number of quaternary states. A mathematical model of RNR activity which explicitly represents the states of R1 currently exists, but it is complicated in several ways: (1) it includes up to six-fold nested sums; (2) it uses different mathematical structures under different substrate-modulator conditions; and (3) it requires root solutions of high order polynomials to determine R1 proportions in mono-, di-, tetra- and hexamer states and thus RNR activity as a function of modulator and total R1 concentrations.     

Structure of the nitrogen-centered radical formed during inactivation of E. coli ribonucleotide reductase by 2'-azido-2'-deoxyuridine-5'-diphosphate: trapping of the 3'-ketonucleotide

Ribonucleotide reductases (RNRs) catalyze the conversion of nucleotides to deoxynucleotides providing the monomeric precursors required for DNA replication and repair. The class I RNRs are composed of two homodimeric subunits: R1 and R2. R1 has the active site where nucleotide reduction occurs, and R2 contains the diiron tyrosyl radical (Y*) cofactor essential for radical initiation on R1. Mechanism-based inhibitors, such as 2'-azido-2'-deoxyuridine-5'-diphosphate (N(3)UDP), have provided much insight into the reduction mechanism. N(3)UDP is a stoichiometric inactivator that, upon interaction with RNR, results in loss of the Y* in R2 and formation of a nitrogen-centered radical (N*) covalently attached to C225 (R-S-N*-X) in the active site of R1. N(2) is lost prior to N* formation, and after its formation, stoichiometric amounts of 2-methylene-3-furanone, pyrophosphate, and uracil are also generated. On the basis of the hyperfine interactions associated with N*, it was proposed that N* is also covalently attached to the nucleotide through either the oxygen of the 3'-OH (R-S-N*-O-R') or the 3'-C (R-S-N*-C-OH). To distinguish between the proposed structures, the inactivation was carried out with 3'-[(17)O]-N(3)UDP and N* was examined by 9 and 140 GHz EPR spectroscopy. Broadening of the N* signal was detected and the spectrum simulated to obtain the [(17)O] hyperfine tensor. DFT calculations were employed to determine which structures are in best agreement with the simulated hyperfine tensor and our previous ESEEM data. The results are most consistent with the R-S-N*-C-OH structure and provide evidence for the trapping of a 3'-ketonucleotide in the reduction process.     

The reaction mechanism of xanthine oxidase: evidence for two-electron chemistry rather than sequential one-electron steps

Current research on xanthine oxidase has favored a mechanism involving base-catalyzed proton abstraction from a Mo-OH group, allowing nucleophilic attack on the substrate and hydride transfer from the substrate to Mo=S group in the active site. During the course of this reaction mechanism, the molybdenum redox cycles from MoVI to MoIV, with reoxidation of the MoIV speices to form the EPR active MoV intermediate. However, it has also been suggested that the reaction occurs in two subsequent one-electron steps. We have determined kinetic parameters kred and kred/Kd for a variety of plausible substrates as well as the one-electron reduction potentials for these substrates. Our data indicate no correlation between these kinetic parameters and their one-electron reduction potentials, as would be expected if the enzyme were using two subsequent one-electron reduction steps. Our results provide additional support to current evidence for the favored two-electron reduction mechanism.     

Conformational remodeling of femtomolar inhibitor-acetylcholinesterase complexes in the crystalline state

The active center of acetylcholinesterase (AChE), a target site for competitive inhibitors, resides centrosymmetric to the subunit at the base of a deep, narrow gorge lined by aromatic residues. At the gorge entry, a peripheral site encompasses overlapping binding loci for noncompetitive inhibitors, which alter substrate access to the gorge. The click-chemistry inhibitor TZ2PA6 links the active center ligand, tacrine, to the peripheral site ligand, propidium, through a biorthogonal reaction of an acetylene and an azide that forms either a syn1 or an anti1 triazole. Compared with wild-type mouse AChE, a Tyr337Ala mutant displays full catalytic activity, albeit with 2-3 orders of magnitude higher affinities for the TZ2PA6 syn1 and anti1 regioisomers, reflected in low femtomolar K(d) values, diffusion-limited association, and dissociation half-times greater than 1 month and 1 week, respectively. Three structures of each of the co-crystallized syn1 and anti1 complexes of the Tyr337Ala mutant were solved at three distinct times of crystal maturation, consistent with or exceeding the half-lives of the complexes in solution, while crystalline complexes obtained from soaked Tyr337Ala crystals led to picturing "freshly formed" complexes. The structures, at 2.55-2.75 ? resolution, reveal a range of unprecedented conformations of the bound regioisomers, not observed in the wild-type AChE complexes, associated with concerted positional rearrangements of side chains in the enzyme gorge. Moreover, time-dependent conformational remodeling of the crystalline complexes appears to correlate with the dissociation half-times of the solution complexes. Hence, for the tight-binding TZ2PA6 inhibitors, the initial complexes kinetically driven in solution slowly form more stable complexes governed by thermodynamic equilibrium and observable in mature crystals.     

Synthesis and characterisation of a ligand that forms a stable tetrahedral intermediate in the active site of the Aureobacterium species (-) gamma-lactamase

The crystal structure of a (-) gamma-lactamase from an Aureobacterium species showed a molecule bound covalently to the active site serine residue. This enzyme complex represented the first structure of a stably bound tetrahedral intermediate for an alpha/beta hydrolase fold enzyme. The structural elucidation of tetrahedral intermediates is important for the understanding of enzymatic mechanism, substrate recognition and enzyme inhibition. In this paper, we report the synthesis and subsequent characterisation of (3aR,7aS)-3a,4,7,7a-tetrahydrobenzo-[1,3]-dioxol-2-one (BD1), the molecule modelled into the Aureobacterium (-) gamma-lactamase active site. This molecule has been confirmed to be an inhibitor and to be displaced from the enzyme by the racemic gamma-lactam substrate.     

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.     

Computational study of ketosteroid isomerase: insights from molecular dynamics simulation of enzyme bound substrate and intermediate

Delta(5)-3-Ketosteroid Isomerase (KSI) catalyzes the isomerization of 5,6-unsaturated ketosteroids to their 4,5-unsaturated isomers at a rate approaching the diffusion limit. The isomerization reaction follows a two-step general acid-base mechanism starting with Asp38-CO(2)(-) mediated proton abstraction from a sp(3)-hybridized carbon atom, alpha to carbonyl group, providing a dienolate intermediate. In the second step, Asp38-CO(2)H protonates the C6 of the intermediate providing a 4,5-unsaturated ketosteroid. The details of the mechanism have been highly controversial despite several experimental and computational studies of this enzyme. The general acid-base catalysis has been proposed to involve either a catalytic diad or a cooperative hydrogen bond mechanism. In this paper, we report our results from the 1.5 nanosecond molecular dynamics (MD) simulation of enzyme bound natural substrate (E.S) and enzyme bound intermediate (E.In) solvated in a TIP3P water box. The final coordinates from our MD simulation strongly support the cooperative hydrogen bond mechanism. The MD simulation of E.S and E.In shows that both Tyr14 and Asp99 are hydrogen bonded to the O3 of the substrate or intermediate. The average hydrogen bonding distance between Tyr14-OH and O3 becomes shorter and exhibits less fluctuation on E.S --> E.In. We also observe dynamic motions of water moving in and out of the active site in the E.S structures. This free movement of water disappears in the E.In structures. The active site is shielded by hydrophobic residues, which come together and squeeze out the waters from the active site in the E.In complex.     

C-H bond activation by a ferric methoxide complex: modeling the rate-determining step in the mechanism of lipoxygenase.

Lipoxygenases are mononuclear non-heme iron enzymes that regio- and stereospecifcally convert 1,4-pentadiene subunit-containing fatty acids into alkyl peroxides. The rate-determining step is generally accepted to be hydrogen atom abstraction from the pentadiene subunit of the substrate by an active ferric hydroxide species to give a ferrous water species and an organic radical. Reported here are the synthesis and characterization of a ferric model complex, [Fe(III)(PY5)(OMe)](OTf)(2), that reacts with organic substrates in a manner similar to the proposed enzymatic mechanism. The ligand PY5 (2,6-bis(bis(2-pyridyl)methoxymethane)pyridine) was developed to simulate the histidine-dominated coordination sphere of mammalian lipoxygenases. The overall monoanionic coordination provided by the endogenous ligands of lipoxygenase confers a strong Lewis acidic character to the active ferric site with an accordingly positive reduction potential. Incorporation of ferrous iron into PY5 and subsequent oxidation yields a stable ferric methoxide species that structurally and chemically resembles the proposed enzymatic ferric hydroxide species. Reactivity with a number of hydrocarbons possessing weak C-H bonds, including a derivative of the enzymatic substrate linoleic acid, scales best with the substrates' bond dissociation energies, rather than pK(a)'s, suggesting a hydrogen atom abstraction mechanism. Thermodynamic analysis of [Fe(III)(PY5)(OMe)](OTf)(2) and the ferrous end-product [Fe(II)(PY5)(MeOH)](OTf)(2) estimates the strength of the O-H bond in the metal bound methanol in the latter to be 83.5 +/- 2.0 kcal mol(-1). The attenuation of this bond relative to free methanol is largely due to the high reduction potential of the ferric site, suggesting that the analogously high reduction potential of the ferric site in LO is what allows the enzyme to perform its unique oxidation chemistry. Comparison of [Fe(III)(PY5)(OMe)](OTf)(2) to other coordination complexes capable of hydrogen atom abstraction shows that, although a strong correlation exists between the thermodynamic driving force of reaction and the rate of reaction, other factors appear to further modulate the reactivity.     

Spectroscopic and electronic structure studies of intermediate X in ribonucleotide reductase R2 and two variants: a description of the FeIV-oxo bond in the FeIII-O-FeIV dimer

Spectroscopic and electronic structure studies of the class I Escherichia coli ribonucleotide reductase (RNR) intermediate X and three computationally derived model complexes are presented, compared, and evaluated to determine the electronic and geometric structure of the FeIII-FeIV active site of intermediate X. Rapid freeze-quench (RFQ) EPR, absorption, and MCD were used to trap intermediate X in R2 wild-type (WT) and two variants, W48A and Y122F/Y356F. RFQ-EPR spin quantitation was used to determine the relative contributions of intermediate X and radicals present, while RFQ-MCD was used to specifically probe the FeIII/FeIV active site, which displayed three FeIV d-d transitions between 16,700 and 22,600 cm(-1), two FeIV d-d spin-flip transitions between 23,500 and 24,300 cm(-1), and five oxo to FeIV and FeIII charge transfer (CT) transitions between 25,000 and 32,000 cm(-1). The FeIV d-d transitions were perturbed in the two variants, confirming that all three d-d transitions derive from the d-pi manifold. Furthermore, the FeIV d-pi splittings in the WT are too large to correlate with a bis-mu-oxo structure. The assignment of the FeIV d-d transitions in WT intermediate X best correlates with a bridged mu-oxo/mu-hydroxo [FeIII(mu-O)(mu-OH)FeIV] structure. The mu-oxo/mu-hydroxo core structure provides an important sigma/pi superexchange pathway, which is not present in the bis-mu-oxo structure, to promote facile electron transfer from Y122 to the remote FeIV through the bent oxo bridge, thereby generating the tyrosyl radical for catalysis.     

Single-molecule Förster resonance energy transfer reveals an innate fidelity checkpoint in DNA polymerase I

Enzymatic reactions typically involve complex dynamics during substrate binding, conformational rearrangement, chemistry, and product release. The noncovalent steps provide kinetic checkpoints that contribute to the overall specificity of enzymatic reactions. DNA polymerases perform DNA replication with outstanding fidelity by actively rejecting noncognate nucleotide substrates early in the reaction pathway. Substrates are delivered to the active site by a flexible fingers subdomain of the enzyme, as it converts from an open to a closed conformation. The conformational dynamics of the fingers subdomain might also play a role in nucleotide selection, although the precise role is currently unknown. Using single-molecule F?rster resonance energy transfer, we observed individual Escherichia coli DNA polymerase I (Klenow fragment) molecules performing substrate selection. We discovered that the fingers subdomain actually samples through three distinct conformations--open, closed, and a previously unrecognized intermediate conformation. We measured the overall dissociation rate of the polymerase-DNA complex and the distribution among the various conformational states in the absence and presence of nucleotide substrates, which were either correct or incorrect. Correct substrates promote rapid progression of the polymerase to the catalytically competent closed conformation, whereas incorrect nucleotides block the enzyme in the intermediate conformation and induce rapid dissociation from DNA. Remarkably, incorrect nucleotide substrates also promote partitioning of DNA to the spatially separated 3'-5' exonuclease domain, providing an additional mechanism to prevent misincorporation at the polymerase active site. These results reveal the existence of an early innate fidelity checkpoint, rejecting incorrect nucleotide substrates before the enzyme encloses the nascent base pair.     

Direct studies on 5-coordinate intermediates formed during substitution at tetrahedral Fe sites: role of bound nucleophile in labilisation of leaving group

The substitution reactions of the tetrahedral Fe sites in [FeCl(4)](-), [Fe(2)S(2)Cl(4)](2-), [Fe(4)S(4)Cl(4)](2-) and [{MoFe(3)S(4)Cl(3)}(2)(micro-SEt)(3)](3-) with 4-RC(6)H(4)S(-) (R = MeO, Me, H, Cl or NO(2)) all involve rapid binding of the thiolate to a Fe site and formation of a kinetically and spectroscopically detectable intermediate. Kinetic studies allow calculation of the rate of Fe-Cl dissociation from the 5-coordinate site of the intermediate (k(2)(R)). The rate of Fe-Cl dissociation from the intermediate exhibits a marked dependence on the nature of the bound thiolate with log(10)(k(2)(R)) increasing in a linear manner with the calculated NBO charge on the sulfur atom of the coordinated thiolate. This behaviour indicates that Fe-Cl bond dissociation at the 5-coordinate intermediate involves a process in which Fe-thiolate bond shortening occurs prior to movement of the Fe-Cl bond.     

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

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

Activation mechanisms of butyrylcholinesterase by 2,4,6-trinitrotoluene, 3,3-dimethylbutyl-N-n-butylcarbamate, and 2-trimethylsilyl-ethyl-N-n-butylcarbamate

The goal of this work was to propose a possible mechanism for the butyrylcholinesterase activation by 2,4,6-trinitrotoluene (TNT), 3,3-dimethylbutyl-N-n-butylcarbamate (1), and 2-trimethylsilyl-ethyl-N-n-butylcarbamate (2). Kinetically, TNT, and compounds 1 and 2 were characterized as the nonessential activators of butyrylcholinesterase. TNT, and compounds 1 and 2 were hydrophobic compounds and were proposed to bind to the hydrophobic activator binding site, which was located outside the active site gorge of the enzyme. The conformational change from a normal active site gorge to a more accessible active site gorge of the enzyme was proposed after binding of TNT, and compounds 1 and 2 to the activator binding site of the enzyme. Therefore, TNT, and compounds 1 and 2 may act as the excess of butyrylcholine in the substrate activator for the butyrylcholinesterase catalyzed reactions.