Product release is rate-limiting in the activation of the prodrug 5-fluorocytosine by yeast cytosine deaminase

Yeast cytosine deaminase (yCD), a zinc metalloenzyme, catalyzes the hydrolytic deamination of cytosine to uracil. The enzyme is of great biomedical interest because it also catalyzes the deamination of the prodrug 5-fluorocytosine (5FC) to form the anticancer drug 5-fluorouracil (5FU). yCD/5FC is one of the most widely used enzyme/prodrug combinations for gene-directed enzyme prodrug therapy for the treatment of cancers. A pH indicator assay has been developed for the measurement of the steady-state kinetic parameters for the deamination reaction. Transient kinetic studies have shown that the product release is a rate-limiting step in the activation of the prodrug 5FC by yCD. The rate constant of the chemical step for the forward reaction (250 s(-)(1)) is approximately 8 times that of the product release (31 s(-)(1)) and approximately 15 times k(cat) (17 s(-)(1)). The transient kinetic results are consistent with those of the steady-state kinetic analysis in the sense that the k(cat) and K(m) values calculated from the rate constants determined by the transient kinetic analysis are in close agreement with those measured by the steady-state kinetic analysis. NMR experiments have demonstrated that free 5FU is in slow exchange with its complex with yCD but has a low affinity for yCD. The transient kinetic and NMR results together suggest that the release of 5FU is rate-limiting in the activation of the prodrug 5FC by yCD and may involve multiple steps.     

Role of glutamate 64 in the activation of the prodrug 5-fluorocytosine by yeast cytosine deaminase

Yeast cytosine deaminase (yCD) catalyzes the hydrolytic deamination of cytosine to uracil as well as the deamination of the prodrug 5-fluorocytosine (5FC) to the anticancer drug 5-fluorouracil. In this study, the role of Glu64 in the activation of the prodrug 5FC was investigated by site-directed mutagenesis, biochemical, nuclear magnetic resonance (NMR), and computational studies. Steady-state kinetics studies showed that the mutation of Glu64 causes a dramatic decrease in k(cat) and a dramatic increase in K(m), indicating Glu64 is important for both binding and catalysis in the activation of 5FC. (19)F NMR experiments showed that binding of the inhibitor 5-fluoro-1H-pyrimidin-2-one (5FPy) to the wild-type yCD causes an upfield shift, indicating that the bound inhibitor is in the hydrated form, mimicking the transition state or the tetrahedral intermediate in the activation of 5FC. However, binding of 5FPy to the E64A mutant enzyme causes a downfield shift, indicating that the bound 5FPy remains in an unhydrated form in the complex with the mutant enzyme. (1)H and (15)N NMR analysis revealed trans-hydrogen bond D/H isotope effects on the hydrogen of the amide of Glu64, indicating that the carboxylate of Glu64 forms two hydrogen bonds with the hydrated 5FPy. ONIOM calculations showed that the wild-type yCD complex with the hydrated form of the inhibitor 1H-pyrimidin-2-one is more stable than the initial binding complex, and in contrast, with the E64A mutant enzyme, the hydrated inhibitor is no longer favored and the conversion has a higher activation energy, as well. The hydrated inhibitor is stabilized in the wild-type yCD by two hydrogen bonds between it and the carboxylate of Glu64 as revealed by (1)H and (15)N NMR analysis. To explore the functional role of Glu64 in catalysis, we investigated the deamination of cytosine catalyzed by the E64A mutant by ONIOM calculations. The results showed that without the assistance of Glu64, both proton transfers before and after the formation of the tetrahedral reaction intermediate become partially rate-limiting steps. The results of the experimental and computational studies together indicate that Glu64 plays a critical role in both the binding and the chemical transformation in the conversion of the prodrug 5FC to the anticancer drug 5-fluorouracil.     

Crystal structure of yeast cytosine deaminase. Insights into enzyme mechanism and evolution

Yeast cytosine deaminase is an attractive candidate for anticancer gene therapy because it catalyzes the deamination of the prodrug 5-fluorocytosine to form 5-fluorouracil. We report here the crystal structure of the enzyme in complex with the inhibitor 2-hydroxypyrimidine at 1.6-A resolution. The protein forms a tightly packed dimer with an extensive interface of 1450 A2 per monomer. The inhibitor was converted into a hydrated adduct as a transition-state analog. The essential zinc ion is ligated by the 4-hydroxyl group of the inhibitor together with His62, Cys91, and Cys94 from the protein. The enzyme shares similar active-site architecture to cytidine deaminases and an unusually high structural homology to 5-aminoimidazole-4-carboxamide-ribonucleotide transformylase and thereby may define a new superfamily. The unique C-terminal tail is involved in substrate specificity and also functions as a gate controlling access to the active site. The complex structure reveals a closed conformation, suggesting that substrate binding seals the active-site entrance so that the catalytic groups are sequestered from solvent. A comparison of the crystal structures of the bacterial and fungal cytosine deaminases provides an elegant example of convergent evolution, where starting from unrelated ancestral proteins, the same metal-assisted deamination is achieved through opposite chiral intermediates within distinctly different active sites.     

Possible ligand release pathway of dipeptidyl peptidase IV investigated by molecular dynamics simulations

Dipeptidyl peptidase IV (DPP4) is an important target for the treatment of Type II diabetes mellitus. The crystal structure of DPP4 demonstrates that there are two possible pathways to the active site, a side opening and a β propeller opening. However, it still lacks quantitative evidence to illustrate which pathway is more favorable for inhibitor to enter into or release from the active site. In this study, conventional and steered molecular dynamics simulations were performed to explore the details of inhibitor Q448 release from the active site of DPP4 via the two potential pathways. The comparisons of force and work together with potentials of mean force results suggested that the side opening might be more favorable for the inhibitor to pass through. Moreover, Glu205-Glu206 and Phe357 were recognized as two "key residues" in the active site for inhibitor binding. Accordingly, suggestions for further inhibitor design were provided.     

Exploring the role of the phospholipid ligand in endothelial protein C receptor: A molecular dynamics study

Endothelial protein C receptor (EPCR) is a CD1‐like transmembrane glycoprotein with important regulatory roles in protein C (PC) pathway, enhancing PC's anticoagulant, anti‐inflammatory, and antiapoptotic activities. Similarly to homologous CD1d, EPCR binds a phospholipid [phosphatidylethanolamine (PTY)] in a groove corresponding to the antigen‐presenting site, although it is not clear if lipid exchange can occur in EPCR as in CD1d. The presence of PTY seems essential for PC γ‐carboxyglutamic acid (Gla) domain binding. However, the lipid‐free form of the EPCR has not been characterized. We have investigated the structural role of PTY on EPCR, by multiple molecular dynamics (MD) simulations of ligand bound and unbound forms of the protein. Structural changes, subsequent to ligand removal, led to identification of two stable and folded ligand‐free conformations. Compared with the bound form, unbound structures showed a narrowing of the A′ pocket and a high flexibility of the helices around it, in agreement with CD1d simulation. Thus, a lipid exchange with a mechanism similar to CD1d is proposed. In addition, unbound conformations presented a reduced interaction surface for Gla domain, confirming the role of PTY in establishing the proper EPCR conformation for the interaction with its partner protein. Single MD simulations were also obtained for 29 mutant models with predicted structural stability and impaired binding ability. Ligand affinity calculations, based on linear interaction energy method, showed that substitution‐induced conformational changes affecting helices around the A′ pocket were associated to a reduced binding affinity. Mutants responsible for this effect may represent useful reagents for experimental tests. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Increased proliferation and chemosensitivity of human mesenchymal stromal cells expressing fusion yeast cytosine deaminase

Mesenchymal stromal cells (MSCs) are considered to be suitable vehicles for cellular therapy in various conditions. The expression of reporter and/or effector protein(s) enabled both the identification of MSCs within the organism and the exploitation in targeted tumor therapies. The aim of this study was to evaluate cellular changes induced by retrovirus-mediated transgene expression in MSCs in vitro. Human Adipose Tissue-derived MSCs (AT-MSCs) were transduced to express (i) the enhanced green fluorescent protein (EGFP) reporter transgene, (ii) the fusion yeast cytosine deaminase::uracil phosphoribosyltransferase (CDy::UPRT) enzyme along with the expression of dominant positive selection gene NeoR or (iii) the selection marker NeoR alone (MOCK). CDy::UPRT expression resulted in increased proliferation of CDy::UPRT-MSCs versus naïve AT-MSCs, MOCK-MSCs or EGFP-MSCs. Furthermore, CDy::UPRT-MSCs were significantly more sensitive to 5-fluorouracil (5FU), cisplatin, cyclophosphamide and cytosine arabinoside as determined by increased Caspase 3/7 activation and/or decreased relative proliferation. CDy::UPRT-MSCs in direct cocultures with breast cancer cells MDA-MB-231 increased tumor cell killing induced by low concentrations of 5FU. Our data demonstrated the changes in proliferation and chemoresistance in engineered MSCs expressing transgene with enzymatic function and suggested the possibilities for further augmentation of targeted MSC-mediated antitumor therapy.     

Inhibition of yeast cytosine deaminase by 5-bromo-2-pyrimidinone and its covalent hydrate

Yeast cytosine deaminase (EC 3.5.4.1) is inhibited by 5-bromo-2-pyrimidinone. In aqueous solution at neutral pH three forms of this compound (the anion, the parent, and the covalent hydrate) are in equilibrium. Experiments were undertaken in order to determine the relative contributions of these three forms to the observed inhibition. The anion makes little or no contribution. Both the parent and the covalent hydrate inhibit the enzyme, with the Ki for the hydrate being 0.2-0.02 times that of the parent. In the presence of stoichiometric concentrations of the enzyme, the equilibrium between parent and hydrate is displaced towards the hydrate; however, the hydration is not catalyzed by cytosine deaminase.     

A kinetic study of the inhibition of yeast AMP deaminase by polyphosphate

Inorganic pyrophosphate and polyphosphates have acted as potent inhibitors of purified AMP deaminase (EC 3.5.4.6) from yeast: the activity fell to a definite limit with the increase in the concentration of the inhibitor. The effect of polyphosphate was largely on the maximal velocity of the enzyme with some decrease in affinity. The cooperative effect of AMP, analyzed in terms of a Hill coefficient, remained at 2 in the absence and presence of polyphosphate. Binding of polyphosphate to the enzyme showed no cooperativity. The inhibition of AMP deaminase by polyphosphate can be qualitatively and quantitatively accounted for by the partial mixed-type inhibition mechanism. Both the Ki value for the inhibitor and the breakdown rate of the enzyme-substrate-inhibitor complex are dependent on the chain length of polyphosphate, suggesting that the breakdown rate of the enzyme-substrate-inhibitor complex is regulated by binding of polyphosphate to a specific inhibitory site.     

Loop motion and base release in purine-specific nucleoside hydrolase: A molecular dynamics study

Although various Trypanosoma vivax purine-specific inosine–adenosine–guanosine nucleoside hydrolase (IAG-NH) crystal structures have been determined and the chemical reaction mechanism of substrate hydrolysis has been studied recently, the mechanistic details for the release of base and ribose are still unclear. Herein molecular dynamics (MD) simulations combined with umbrella sampling technique were utilized to explore the regulation mechanisms of key residues and loops 1 and 2 for the base release. Our results have indicated that the base release process is not the rate-limiting step in the entire hydrolysis process, and the very low barrier of ~ 5.6 kcal/mol can be washed out easily by the notable exothermicity from the substrate hydrolysis step. Moreover, the MD simulations have revealed that Glu82/Trp83 in loop 1 and His247/Arg252 in loop 2 are important to modulate the base release. The partial helix-to-coil change of loop 2 along with the base release process has been observed, showing good agreement with the IAG-NH crystal structures. The local binding site around the ribose after the base release is also discussed.     

A negative selection scheme based on the expression of cytosine deaminase in plastids

The enzyme cytosine deaminase (CD) encoded by codA catalyzes deamination of cytosine to uracil. CD is present in prokaryotes and in many eukaryotic micro-organisms, but is absent in higher plants. 5-fluorocytosine (5FC) is metabolized in CD-expressing cells, causing cellular death. A chimeric codA has been introduced into the tobacco plastid genome and 5FC was used to select against tissue culture cells and seedlings expressing CD. This negative selection scheme will be useful in identifying nuclear genes which control plastid gene expression in higher plants.     

Location of the gene specifying cytosine deaminase in Escherichia coli

Molecular Genetics and Genomics - The gene specifying cytosine deaminase (cod) is shown to be located at approximately 86 minutes on the linkage map of E. coli. The corresponding gene in S....     

Recombinogenic phenotype of human activation-induced cytosine deaminase

Class switch recombination, gene conversion, and somatic hypermutation that diversify rearranged Ig genes to produce various classes of high affinity Abs are dependent on the enzyme activation-induced cytosine deaminase (AID). Evidence suggests that somatic hypermutation is due to error-prone DNA repair that is initiated by AID-mediated deamination of cytosine in DNA, whereas the mechanism by which AID controls recombination remains to be elucidated. In this study, using a yeast model system, we have observed AID-dependent recombination. Expression of human AID in wild-type yeast is mutagenic for G-C to A-T transitions, and as expected, this mutagenesis is increased upon inactivation of uracil-DNA glycosylase. AID expression also strongly induces intragenic mitotic recombination, but only in a strain possessing uracil-DNA glycosylase. Thus, the initial step of base excision repair is required for AID-dependent recombination and is a branch point for either hypermutagenesis or recombination.     

The 1.14 A crystal structure of yeast cytosine deaminase: evolution of nucleotide salvage enzymes and implications for genetic chemotherapy

Cytosine deaminase (CD) catalyzes the deamination of cytosine and is only present in prokaryotes and fungi, where it is a member of the pyrimidine salvage pathway. The enzyme is of interest both for antimicrobial drug design and gene therapy applications against tumors. The structure of Saccharomyces cerevisiae CD has been determined in the presence and absence of a mechanism-based inhibitor, at 1.14 and 1.43 A resolution, respectively. The enzyme forms an alpha/beta fold similar to bacterial cytidine deaminase, but with no similarity to the alpha/beta barrel fold used by bacterial cytosine deaminase or mammalian adenosine deaminase. The structures observed for bacterial, fungal, and mammalian nucleic acid deaminases represent an example of the parallel evolution of two unique protein folds to carry out the same reaction on a diverse array of substrates.     

AMP deaminase from baker's yeast. Kinetic and molecular properties.

The kinetic and molecular properties of AMP deaminase [AMP aminohydrolase, EC 3.5.4.6] purified from baker's yeast (saccharomyces cerevisiae) were investigated. The enzyme was activated by ATP and dATP, but inhibited by Pi and GTP in an allosteric manner. Alkali metal ions and alkaline earth metal ions activated the enzyme to various extent. Kinetic negative cooperativity was observed in the binding of nucleoside triphosphates. Kinetic analysis showed that the number of interaction sites for AMP (substrate) and Pi (inhibitor) is two each per enzyme molecule. The molecular weight of the native enzyme was estimated to be 360,000 by sedimentation equilibrium studies. On polyacrylamide gel electrophoresis in the presence of sodium dodecyl sulfate, the enzyme gave a single polypeptide band with a molecular weight of 83,000, suggesting that the native enzyme has a tetrameric structure. Baker's yeast AMP deaminase was concluded to consist of two "promoter" units which each consist of two polypeptide chains with identical molecular weight.     

Proton transfer dynamics in the propionic acid dimer from path integral molecular dynamics calculations

The double proton transfer process in the cyclic dimer of propionic acid in the gas phase was studied using a path integral molecular dynamics method. Structures, energies and proton trajectories were determined. Very large amplitude motions of the skeleton of a propionic acid molecule were observed during the simulations, and almost free rotation of the C₂H₅ group around the C_α-C bond. A double-well symmetric potential with a very small energy barrier was determined from the free energy profile for the proton motions. Infrared spectra for different isotopomers were calculated, and comparative vibrational analysis was performed. The vibrational results from CPMD appear to be in qualitative agreement with the experimental ones.     

Location of gene specifying cytosine deaminase in Escherichia coli

Summary The gene specifying cytosine deaminase (cod) is shown to be located at approximately 86 minutes on the linkage map of E. coli. The corresponding gene in S. typhimurium has been reported to have a different location.     

Investigation of ligand selectivity in CYP3A7 by molecular dynamics simulations

Cytochrome P450 (CYP) 3A7 plays a crucial role in the biotransformation of the metabolized endogenous and exogenous steroids. To compare the metabolic capabilities of CYP3A7–ligands complexes, three endogenous ligands were selected, namely dehydroepiandrosterone (DHEA), estrone, and estradiol. In this study, a three-dimensional model of CYP3A7 was constructed by homology modeling using the crystal structure of CYP3A4 as the template and refined by molecular dynamics simulation (MD). The docking method was adopted, combined with MD simulation and the molecular mechanics generalized born surface area method, to probe the ligand selectivity of CYP3A7. These results demonstrate that DHEA has the highest binding affinity, and the results of the binding free energy were in accordance with the experimental conclusion that estrone is better than estradiol. Moreover, several key residues responsible for substrate specificity were identified on the enzyme. Arg372 may be the most important residue due to the low interaction energies and the existence of hydrogen bond with DHEA throughout simulation. In addition, a cluster of Phe residues provides a hydrophobic environment to stabilize ligands. This study provides insights into the structural features of CYP3A7, which could contribute to further understanding of related protein structures and dynamics.