Thymidylate synthase: structure, inhibition, and strained conformations during catalysis

A computational investigation of thymidylate synthase inhibitors through a combined approach of 3D-QSAR and pharmacophore modelling

Thymidylate synthase (TS) is a crucial target of cancer drug discovery and is mainly involved in the De novo synthesis of the DNA precursor thymine. In the present study, to generate reliable models and identify a few promising molecules, we combined QSAR modelling with the pharmacophore hypothesis-generating technique. Input molecules were clustered on their similarity, and a cluster of 74 molecules with a pyrimidine moiety was chosen as the set for 3D-QSAR and pharmacophore modelling. Atom-based and field-based 3D-QSAR models were generated and statistically validated with R2 > 0.90 and Q2 > 0.75. The common pharmacophore hypothesis(CPH) generation identified the best six-point model ADHRRR. Using these best models, a library of FDA-approved drugs was screened for activity and filtered via molecular docking, ADME profiling, and molecular dynamics simulations. The top ten promising TS-inhibiting candidates were identified, and their chemical features profitable for TS inhibitors were explored.Communicated by Ramaswamy H. Sarma.

Vitamins as Possible Cancer Biomarkers: Significance and Limitations

The Western-style diet, which is common in developed countries and spreading into developing countries, is unbalanced in many respects. For instance, micronutrients (vitamins A, B complex, C, D, E, and K plus iron, zinc, selenium, and iodine) are generally depleted in Western food (causing what is known as ‘hidden hunger’), whereas some others (such as phosphorus) are added beyond the daily allowance. This imbalance in micronutrients can induce cellular damage that can increase the risk of cancer. Interestingly, there is a large body of evidence suggesting a strong correlation between vitamin intake as well as vitamin blood concentrations with the occurrence of certain types of cancer. The direction of association between the concentration of a given vitamin and cancer risk is tumor specific. The present review summarized the literature regarding vitamins and cancer risk to assess whether these could be used as diagnostic or prognostic markers, thus confirming their potential as biomarkers. Despite many studies that highlight the importance of monitoring vitamin blood or tissue concentrations in cancer patients and demonstrate the link between vitamin intake and cancer risk, there is still an urgent need for more data to assess the effectiveness of vitamins as biomarkers in the context of cancer. Therefore, this review aims to provide a solid basis to support further studies on this promising topic.

[Expression of thymidylate synthase in salivary adenoid myoepithelial cells and its clinical significance]

OBJECTIVE

To evaluate the expression of thymidylate synthase (TS) in myoepithelial cells (MECs) of salivary adenoid tissues and explore its clinical significance.

METHODS

Immunohistochemical staining EnVision method was used to detect the expression of TS, P63, Calponin, CK5/6 and S-100 in 32 salivary gland specimens, including 10 non-neoplastic and salivary inflammation specimens, 11 mixed tumor specimens, 5 basal cell carcinoma specimens and 6 adenoid cyst carcinoma specimens. The specificity and sensitivity of TS as a specific molecular marker of salivary muscle epithelial cells were evaluated in comparison with P63, Calponin, CK5/6 and S-100.

RESULTS

The expression pattern of TS in all the salivary gland tissue specimens was identical with that of p63. TS and P63 both showed strong immunohistochemical expressions in MECs of salivary adenoid tissue specimens. Calponin, CK5/6, and S-100 showed cytoplasmic/membranous expressions in the MECs. In addition, TS exhibited weak or moderate cytoplasmic expression in a few salivary gland epithelial cells, cancer cells and scattered stromal cells, with negative expression in the cell nuclei. The expression of TS in the MECs of all the salivary adenoid specimens was highly consistent with those of P63, Calponin, CK5/6 and S-100 (P>0.05) Except for CK5/6 expression in Salivary inflammation and Salivary gland specimens. Kappa>0.75. The specificity and sensitivity of TS as a molecular marker of MECs were both 100%.

CONCLUSIONS

TS is a new specific marker of MECs for differential diagnosis of salivary gland tumors.

Crystal structure of the active form of native human thymidylate synthase in the absence of bound substrates

Human thymidylate synthase (hTS) provides the sole de novo intracellular source of thymidine 5'-monophosphate (dTMP). hTS is required for DNA replication prior to cell division, making it an attractive target for anticancer chemotherapy and drug discovery. hTS binds 2'-deoxyuridine 5'-monophosphate (dUMP) and the folate co-substrate N5,N10-methylenetetrahydrofolate (meTHF) in a pocket near the catalytic residue Cys195. The catalytic loop, which is composed of amino-acid residues 181-197, can adopt two distinct conformations related by a 180° rotation. In the active conformation Cys195 is close to the active site, while in the inactive conformation it is rotated and Cys195 is too distant from the active site for catalysis. Several hTS structures, either native or engineered, have been solved in the active conformation in complex with ligands or inhibitors and at different salt concentrations. However, apo hTS structures have been solved in an inactive conformation in high-salt and low-salt conditions (PDB entries 1ypv, 4h1i, 4gyh, 3egy and 3ehi). Here, the structure of apo hTS crystallized in the active form with sulfate ions coordinated by the arginine residue that binds dUMP is reported.

Inflammation-associated nitrotyrosination affects TCR recognition through reduced stability and alteration of the molecular surface of the MHC complex

Nitrotyrosination of proteins, a hallmark of inflammation, may result in the production of MHC-restricted neoantigens that can be recognized by T cells and bypass the constraints of immunological self-tolerance. Here we biochemically and structurally assessed how nitrotyrosination of the lymphocytic choriomeningitis virus (LCMV)-associated immunodominant MHC class I-restricted epitopes gp33 and gp34 alters T cell recognition in the context of both H-2D^b and H-2K^b. Comparative analysis of the crystal structures of H-2K^b/gp34 and H-2K^b/NY-gp34 demonstrated that nitrotyrosination of p3Y in gp34 abrogates a hydrogen bond interaction formed with the H-2K^b residue E152. As a consequence the conformation of the TCR-interacting E152 was profoundly altered in H-2K^b/NY-gp34 when compared to H-2K^b/gp34, thereby modifying the surface of the nitrotyrosinated MHC complex. Furthermore, nitrotyrosination of gp34 resulted in structural over-packing, straining the overall conformation and considerably reducing the stability of the H-2K^b/NY-gp34 MHC complex when compared to H-2K^b/gp34. Our structural analysis also indicates that nitrotyrosination of the main TCR-interacting residue p4Y in gp33 abrogates recognition of H-2D^b/gp33-NY complexes by H-2D^b/gp33-specific T cells through sterical hindrance. In conclusion, this study provides the first structural and biochemical evidence for how MHC class I-restricted nitrotyrosinated neoantigens may enable viral escape and break immune tolerance.

Affinity characterization-mass spectrometry methodology for quantitative analyses of small molecule protein binding in solution

Affinity characterization by mass spectrometry (AC-MS) is a novel LC-MS methodology for quantitative determination of small molecule ligand binding to macromolecules. Its most distinguishing feature is the direct determination of all three concentration terms of the equilibrium binding equation, i.e., (M), (L), and (ML), which denote the macromolecule, ligand, and the corresponding complex, respectively. Although it is possible to obtain the dissociation constant from a single mixing experiment, saturation analyses are still valuable for assessing the overall binding phenomenon based on an established formalism. In addition to providing the prerequisite dissociation constant and binding stoichiometry, the technique also provides valuable information about the actual solubility of both macromolecule and ligand upon dilution and mixing in binding buffers. The dissociation constants and binding mode for interactions of DNA primase and thymidylate synthetase (TS) with high and low affinity small molecule ligands were obtained using the AC-MS method. The data were consistent with the expected affinity of TS for these ligands based on dissociation constants determined by alternative thermal-denaturation techniques: TdF or TdCD, and also consistent enzyme inhibition constants reported in the literature. The validity of AC-MS was likewise extended to a larger set of soluble protein-ligand systems. It was established as a valuable resource for counter screen and structure-activity relationship studies in drug discovery, especially when other classical techniques could only provide ambiguous results.

Design, synthesis, biological evaluation and X-ray crystal structure of novel classical 6,5,6-tricyclic benzo[4,5]thieno[2,3-d]pyrimidines as dual thymidylate synthase and dihydrofolate reductase inhibitors

Classical antifolates (4-7) with a tricyclic benzo[4,5]thieno[2,3-d]pyrimidine scaffold and a flexible and rigid benzoylglutamate were synthesized as dual thymidylate synthase (TS) and dihydrofolate reductase (DHFR) inhibitors. Oxidative aromatization of ethyl 2-amino-4-methyl-4,5,6,7-tetrahydro-1-benzothiophene-3-carboxylate (±)-9 to ethyl 2-amino-4-methyl-1-benzothiophene-3-carboxylate 10 with 10% Pd/C was a key synthetic step. Compounds with 2-CH₃ substituents inhibited human (h) TS (IC₅₀ =0.26-0.8 μM), but not hDHFR. Substitution of the 2-CH₃ with a 2-NH₂ increases hTS inhibition by more than 10-fold and also affords excellent hDHFR inhibition (IC₅₀ = 0.09-0.1 μM). This study shows that the tricyclic benzo[4,5]thieno[2,3-d]pyrimidine scaffold is highly conducive to single hTS or dual hTS-hDHFR inhibition depending on the 2-position substituents. The X-ray crystal structures of 6 and 7 with hDHFR reveal, for the first time, that tricyclics 6 and 7 bind with the benzo[4,5]thieno[2,3-d]pyrimidine ring in the folate binding mode with the thieno S mimicking the 4-amino of methotrexate.

Structures of human thymidylate synthase R163K with dUMP, FdUMP and glutathione show asymmetric ligand binding

Thymidylate synthase (TS) is a well validated target in cancer chemotherapy. Here, a new crystal form of the R163K variant of human TS (hTS) with five subunits per asymmetric part of the unit cell, all with loop 181-197 in the active conformation, is reported. This form allows binding studies by soaking crystals in artificial mother liquors containing ligands that bind in the active site. Using this approach, crystal structures of hTS complexes with FdUMP and dUMP were obtained, indicating that this form should facilitate high-throughput analysis of hTS complexes with drug candidates. Crystal soaking experiments using oxidized glutathione revealed that hTS binds this ligand. Interestingly, the two types of binding observed are both asymmetric. In one subunit of the physiological dimer covalent modification of the catalytic nucleophile Cys195 takes place, while in another dimer a noncovalent adduct with reduced glutathione is formed in one of the active sites.

rTSβ as a Novel 5-fluorouracil Resistance Marker of Colorectal Cancer: A Preliminary Study

Introduction: Colorectal cancer is the most common form of malignancy in Taiwan and the third leading cause of death from cancer, preceded only by lung and hepatic cancers. Colorectal cancer is typically treated by surgical intervention and/or chemotherapy and radiotherapy, if necessary. To date, 5-fluorouracil (5-FU) is the most commonly used anti-cancer chemotherapy drug. However, patients commonly experience resistance to the drug therefore limiting its efficiency. In this study, we measured the expression of rTSβ in human colon cancer as a novel 5-FU resistance marker. Materials and Methods: We collected 172 colon cancer samples from 4 different hospitals (including 21 pairs of colon cancer biopsies and 151 pathologic slides of colon cancer). In vitro, we measured the cytotoxicity of 5-FU and 5-FU plus leucovorin in H630 and H630-1 colon cancer cell lines. Results: The results revealed that rTSβ was expressed in 115 (66.9 %) pathology samples and that tumour expression was higher than in corresponding normal tissue. Survival rates of up to 5 years following treatment was significantly higher for patients without rTSβ expression than for those with rTSβ expression (P = 0.0023). In vitro, H630-1 (with rTSβ overexpression) had significantly higher IC50 of 5-FU than did H630. IC50of 5-FU decreased when leucovorin was added. Conclusions: Results indicate a close relationship between rTSβ expression and resistance to the drug 5-FU in human colorectal cancer. These results provide further evidence for rTSβ expression as a novel 5-FU resistance marker of colorectal cancer. Key words: Chemotherapy, Immunohistochemistry, Prognosis

Variants of human thymidylate synthase with loop 181-197 stabilized in the inactive conformation

Loop 181-197 of human thymidylate synthase (hTS) populates two major conformations, essentially corresponding to the loop flipped by 180 degrees . In one of the conformations, the catalytic Cys195 residue lies distant from the active site making the enzyme inactive. Ligands stabilizing this inactive conformation may function as allosteric inhibitors. To facilitate the search for such inhibitors, we have expressed and characterized several mutants designed to shift the equilibrium toward the inactive conformer. In most cases, the catalytic efficiency of the mutants was only somewhat impaired with values of k(cat)/K(m) reduced by factors in a 2-12 range. One of the mutants, M190K, is however unique in having the value of k(cat)/K(m) smaller by a factor of approximately 7500 than the wild type. The crystal structure of this mutant is similar to that of the wt hTS with loop 181-197 in the inactive conformation. However, the direct vicinity of the mutation, residues 188-194 of this loop, assumes a different conformation with the positions of C(alpha) shifted up to 7.2 A. This affects region 116-128, which became ordered in M190K while it is disordered in wt. The conformation of 116-128 is however different than that observed in hTS in the active conformation. The side chain of Lys190 does not form contacts and is in solvent region. The very low activity of M190K as compared to another mutant with a charged residue in this position, M190E, suggests that the protein is trapped in an inactive state that does not equilibrate easily with the active conformer.

Direct role of nucleotide metabolism in C-MYC-dependent proliferation of melanoma cells

To identify C-MYC targets rate-limiting for proliferation of malignant melanoma, we stably inhibited C-MYC in several human metastatic melanoma lines via lentivirus-based shRNAs approximately to the levels detected in normal melanocytes. C-MYC depletion did not significantly affect levels of E2F1 protein reported to regulate expression of many S-phase specific genes, but resulted in the repression of several genes encoding enzymes rate-limiting for dNTP metabolism. These included thymidylate synthase (TS), inosine monophosphate dehydrogenase 2 (IMPDH2) and phosphoribosyl pyrophosphate synthetase 2 (PRPS2). C-MYC depletion also resulted in reduction in the amounts of deoxyribonucleoside triphosphates (dNTPs) and inhibition of proliferation. shRNA-mediated suppression of TS, IMPDH2 or PRPS2 resulted in the decrease of dNTP pools and retardation of the cell cycle progression of melanoma cells in a manner similar to that of C-MYC-depletion in those cells. Reciprocally, concurrent overexpression of cDNAs for TS, IMPDH2 and PRPS2 delayed proliferative arrest caused by inhibition of C-MYC in melanoma cells. Overexpression of C-MYC in normal melanocytes enhanced expression of the above enzymes and increased individual dNTP pools. Analysis of in vivo C-MYC interactions with TS, IMPDH2 and PRPS2 genes confirmed that they are direct C-MYC targets. Moreover, all three proteins express at higher levels in cells from several metastatic melanoma lines compared to normal melanocytes. Our data establish a novel functional link between C-MYC and dNTP metabolism and identify its role in proliferation of tumor cells.

The R163K mutant of human thymidylate synthase is stabilized in an active conformation: structural asymmetry and reactivity of cysteine 195

Loop 181-197 of human thymidylate synthase (hTS) populates two conformational states. In the first state, Cys195, a residue crucial for catalytic activity, is in the active site (active conformer); in the other conformation, it is about 10 A away, outside the active site (inactive conformer). We have designed and expressed an hTS variant, R163K, in which the inactive conformation is destabilized. The activity of this mutant is 33% higher than that of wt hTS, suggesting that at least one-third of hTS populates the inactive conformer. Crystal structures of R163K in two different crystal forms, with six and two subunits per asymmetric part of the unit cells, have been determined. All subunits of this mutant are in the active conformation while wt hTS crystallizes as the inactive conformer in similar mother liquors. The structures show differences in the environment of catalytic Cys195, which correlate with Cys195 thiol reactivity, as judged by its oxidation state. Calculations show that the molecular electrostatic potential at Cys195 differs between the subunits of the dimer. One of the dimers is asymmetric with a phosphate ion bound in only one of the subunits. In the absence of the phosphate ion, that is in the inhibitor-free enzyme, the tip of loop 47-53 is about 11 A away from the active site.

Sequence analysis and functional study of thymidylate synthase from zebrafish, Danio rerio

The thymidylate synthase (TS), an important target for many anticancer drugs, has been cloned from different species. But the cDNA property and function of TS in zebrafish are not well documented. In order to use zebrafish as an animal model for screening novel anticancer agents, we isolated TS cDNA from zebrafish and compared its sequence with those from other species. The open reading frame (ORF) of zebrafish TS cDNA sequence was 954 nucleotides, encoding a 318-amino acid protein with a calculated molecular mass of 36.15 kDa. The deduced amino acid sequence of zebrafish TS was similar to those from other organisms, including rat, mouse and humans. The zebrafish TS protein was expressed in Escherichia coli and purified to homogeneity. The purified zebrafish TS showed maximal activity at 28 degrees C with similar K(m) value to human TS. Western immunoblot assay confirmed that TS was expressed in all the developmental stages of zebrafish with a high level of expression at the 1-4 cell stages. To study the function of TS in zebrafish embryo development, a short hairpin RNA (shRNA) expression vector, pSilencer 4.1-CMV/TS, was constructed which targeted the protein-coding region of zebrafish TS mRNA. Significant change in the development of tail and epiboly was found in zebrafish embryos microinjected pSilencer4.1-CMV/TS siRNA expression vector.

The role of protein dynamics in thymidylate synthase catalysis: variants of conserved 2'-deoxyuridine 5'-monophosphate (dUMP)-binding Tyr-261

The enzyme thymidylate synthase (TS) catalyzes the reductive methylation of 2'-deoxyuridine 5'-monophosphate (dUMP) to 2'-deoxythymidine 5'-monophosphate. Using kinetic and X-ray crystallography experiments, we have examined the role of the highly conserved Tyr-261 in the catalytic mechanism of TS. While Tyr-261 is distant from the site of methyl transfer, mutants at this position show a marked decrease in enzymatic activity. Given that Tyr-261 forms a hydrogen bond with the dUMP 3'-O, we hypothesized that this interaction would be important for substrate binding, orientation, and specificity. Our results, surprisingly, show that Tyr-261 contributes little to these features of the mechanism of TS. However, the residue is part of the structural core of closed ternary complexes of TS, and conservation of the size and shape of the Tyr side chain is essential for maintaining wild-type values of kcat/Km. Moderate increases in Km values for both the substrate and cofactor upon mutation of Tyr-261 arise mainly from destabilization of the active conformation of a loop containing a dUMP-binding arginine. Besides binding dUMP, this loop has a key role in stabilizing the closed conformation of the enzyme and in shielding the active site from the bulk solvent during catalysis. Changes to atomic vibrations in crystals of a ternary complex of Escherichia coli Tyr261Trp are associated with a greater than 2000-fold drop in kcat/Km. These results underline the important contribution of dynamics to catalysis in TS.

Crystal structures of thymidylate synthase mutant R166Q: structural basis for the nearly complete loss of catalytic activity

Thymidylate synthase (TS) catalyzes the folate-dependent methylation of deoxyuridine monophosphate (dUMP) to form thymidine monophosphate (dTMP). We have investigated the role of invariant arginine 166, one of four arginines that contact the dUMP phosphate, using site-directed mutagenesis, X-ray crystallography, and TS from Escherichia coli. The R166Q mutant was crystallized in the presence of dUMP and a structure determined to 2.9 A resolution, but neither the ligand nor the sulfate from the crystallization buffer was found in the active site. A second structure determined with crystals prepared in the presence of dUMP and the antifolate 10-propargyl-5,8-dideazafolate revealed that the inhibitor was bound in an extended, nonproductive conformation, partially occupying the nucleotide-binding site. A sulfate ion, rather than dUMP, was found in the nucleotide phosphate-binding site. Previous studies have shown that the substitution at three of the four arginines of the dUMP phosphate-binding site is permissive; however; for Arg166, all the mutations lead to a near-inactive mutant. The present structures of TS R166Q reveal that the phosphate-binding site is largely intact, but with a substantially reduced affinity for phosphate, despite the presence of the three remaining arginines. The position of Cys146, which initiates catalysis, is shifted in the mutant and resides in a position that interferes with the binding of the dUMP pyrimidine moiety.

Analysis of target genes induced by IL-13 cytotoxin in human glioblastoma cells

IL-13 cytotoxin comprised of IL-13 and a mutated form of Pseudomonas exotoxin (fusion protein termed IL-13-PE38QQR) has been shown to inhibit protein synthesis leading to necrotic and apoptotic cell death in glioblastoma cells that express high levels of interleukin-13 receptors (IL-13R). To identify target genes of cell death and other cellular genes with IL-13 receptors in glioblastoma cells, we utilized the cDNA microarrays to analyze global gene expression profiles after IL-13 cytotoxin and IL-13 treatment. IL-13 cytotoxin mediated cytotoxicity to U251 cells in a dose-dependent manner. Hierarchical cluster analysis of differentially expressed genes in U251 glioma cells at different time points after IL-13 cytotoxin treatment showed three major groups, each representing a specific expression pattern. Randomly selected differentially expressed genes from each group were confirmed by RT-PCR analysis. Most down-regulated genes belong to cell adhesion, motility, angiogenesis, DNA repair, and metabolic pathways. While up-regulated genes belong to cell cycle arrest, apoptosis, signaling and various metabolic pathways. Unexpectedly, at early time points, both IL-13 and IL-13 cytotoxin induced several genes belonging to different pathways most notably IL-8, DIO2, END1, and ALDH1A3 indicating that these genes are early response genes and their products may be associated with IL-13R. In addition, IL-13 cytotoxin induced IL-13Ralpha2 mRNA expression during the treatment in glioma cells. Our results indicate that novel cellular genes are involved with IL-13 receptors and that IL-13 cytotoxin induced cell death involves various target genes in human glioblastoma cells. On going studies will determine the role of associated genes and their products in the IL-13R functions in glioma cells.

5,10-Methylene-5,6,7,8-tetrahydrofolate conformational transitions upon binding to thymidylate synthase: molecular mechanics and continuum solvent studies

We applied the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach to evaluate relative stability of the extended (flat) and C-shaped (bent) solution conformational forms of the 5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) molecule in aqueous solution. Calculations indicated that both forms have similar free energies in aqueous solution but detailed energy components are different. The bent solution form has lower intramolecular electrostatic and van der Waals interaction energies. The flat form has more favorable solvation free energy and lower contribution from the bond, angle and torsion angle molecular mechanical internal energies. We exploit these results and combine them with known crystallographic data to provide a model for the progressive binding of the mTHF molecule, a natural cofactor of thymidylate synthase (TS), to the complex forming in the TS-catalyzed reaction. We propose that at the time of initial weak binding in the open enzyme the cofactor molecule remains in a close balance between the flat and bent solution conformations, with neither form clearly favored. Later, thymidylate synthase undergoes conformational change leading to the closure of the active site and the mTHF molecule is withdrawn from the solvent. That effect shifts the thermodynamic equilibrium of the mTHF molecule toward the bent solution form. At the same time, burying the cofactor molecule in the closed active site produces numerous contacts between mTHF and protein that render change in the shape of the mTHF molecule. As a result, the bent solution conformer is converted to more strained L-shaped bent enzyme conformer of the mTHF molecule. The strain in the bent enzyme conformation allows for the tight binding of the cofactor molecule to the productive ternary complex that forms in the closed active site, and facilitates the protonation of the imidazolidine N10 atom, which promotes further reaction.

Catalytic mechanism of Chlamydia trachomatis flavin-dependent thymidylate synthase

Here we report on a Chlamydia trachomatis gene that complements the growth defect of a thymidylate synthase-deficient strain of Escherichia coli. The complementing gene encodes a 60.9-kDa protein that shows low level primary sequence homology to a new class of thymidylate-synthesizing enzymes, termed flavin-dependent thymidylate synthases (FDTS). Purified recombinant chlamydial FDTS (CTThyX) contains bound flavin. Results with site-directed mutants indicate that highly conserved arginine residues are required for flavin binding. Kinetic characterization indicates that CTThyX is active as a tetramer with NADPH, methylenetetrahydrofolate, and dUMP required as substrates, serving as source of reducing equivalents, methyl donor, and methyl acceptor, respectively. dTMP and H(4)folate are products of the reaction. Production of H(4)folate rather than H(2)folate, as in the classical thymidylate synthase reaction, eliminates the need for dihydrofolate reductase, explaining the trimethoprim-resistant phenotype displayed by thyA(-) E. coli-expressing CTThyX. In contrast to the extensively characterized thyA-encoded thymidylate synthases, which form a ternary complex with substrates dUMP and CH(2)H(4)folate and follow an ordered sequential mechanism, CTThyX follows a ping-pong kinetic mechanism involving a methyl enzyme intermediate. Mass spectrometry was used to localize the methyl group to a highly conserved arginine, and site-directed mutagenesis showed this arginine to be critical for thymidylate synthesizing activity. These differentiating characteristics clearly distinguish FDTS from ThyA, making this class of enzymes attractive targets for rational drug design.

Testing a Flexible-receptor Docking Algorithm in a Model Binding Site

Sampling receptor flexibility is challenging for database docking. We consider a method that treats multiple flexible regions of the binding site independently, recombining them to generate different discrete conformations. This algorithm scales linearly rather than exponentially with the receptor's degrees of freedom. The method was first evaluated for its ability to identify known ligands of a hydrophobic cavity mutant of T4 lysozyme (L99A). Some 200000 molecules of the Available Chemical Directory (ACD) were docked against an ensemble of cavity conformations. Surprisingly, the enrichment of known ligands from among a much larger number of decoys in the ACD was worse than simply docking to the apo conformation alone. Large decoys, accommodated in the larger cavity conformations sampled in the ensemble, were ranked better than known small ligands. The calculation was redone with an energy correction term that considered the cost of forming the larger cavity conformations. Enrichment improved, as did the balance between high-ranking large and small ligands. In a second retrospective test, the ACD was docked against a conformational ensemble of thymidylate synthase. Compared to docking against individual enzyme conformations, the flexible receptor docking approach improved enrichment of known ligands. Including a receptor conformational energy weighting term improved enrichment further. To test the method prospectively, the ACD database was docked against another cavity mutant of lysozyme (L99A/M102Q). A total of 18 new compounds predicted to bind this polar cavity and to change its conformation were tested experimentally; 14 were found to bind. The bound structures for seven ligands were determined by X-ray crystallography. The predicted geometries of these ligands all corresponded to the observed geometries to within 0.7A RMSD or better. Significant conformational changes of the cavity were observed in all seven complexes. In five structures, part of the observed accommodations were correctly predicted; in two structures, the receptor conformational changes were unanticipated and thus never sampled. These results suggest that although sampling receptor flexibility can lead to novel ligands that would have been missed when docking a rigid structure, it is also important to consider receptor conformational energy.

Purification and characterization of the aromatic desulfinase, 2-(2'-hydroxyphenyl)benzenesulfinate desulfinase

2-(2(')-Hydroxyphenyl)benzenesulfinate desulfinase (HPBS desulfinase) catalyzes the cleavage of the carbon-sulfur bond of 2-(2(')-hydroxyphenyl)benzenesulfinate (HPBS) to form hydroxybiphenyl and sulfite. This is the final step in the desulfurization of dibenzothiophene, the organosulfur compound used to study biodesulfurization of petroleum middle distillate. HPBS desulfinase was purified 1600-fold from Rhodococcus IGTS8. The purification was monitored using a spectrofluorimetric assay and SDS-PAGE. The pI of HPBS desulfinase is 5.6, the temperature optimum is 35 degrees C, and the pH optimum is 7.0. HPBS desulfinase has a K(m) of 0.90+/-0.15 microM and a k(cat) of 1.3+/-0.07 min(-1). Several analogs were tested for their ability to act as substrates or inhibitors of HPBS desulfinase. No alternative substrates and very few inhibitors were identified. HPBS desulfinase activity decreases in the presence of Cu(2+) and Zn(2+), while no metals significantly enhance enzyme activity. HPBS desulfinase is susceptible to tyrosine, tryptophan, and cysteine specific modification agents.

Approaches to solving the rigid receptor problem by identifying a minimal set of flexible residues during ligand docking

BACKGROUND

Using fixed receptor sites derived from high-resolution crystal structures in structure-based drug design does not properly account for ligand-induced enzyme conformational change and imparts a bias into the discovery and design of novel ligands. We sought to facilitate the design of improved drug leads by defining residues most likely to change conformation, and then defining a minimal manifold of possible conformations of a target site for drug design based on a small number of identified flexible residues.

RESULTS

The crystal structure of thymidylate synthase from an important pathogenic target Pneumocystis carinii (PcTS) bound to its substrate and the inhibitor, BW1843U89, is reported here and reveals a new conformation with respect to the structure of PcTS bound to substrate and the more conventional antifolate inhibitor, CB3717. We developed an algorithm for determining which residues provide 'soft spots' in the protein, regions where conformational adaptation suggests possible modifications for a drug lead that may yield higher affinity. Remodeling the active site of thymidylate synthase with new conformations for only three residues that were identified with this algorithm yields scores for ligands that are compatible with experimental kinetic data.

CONCLUSIONS

Based on the examination of many protein/ligand complexes, we develop an algorithm (SOFTSPOTS) for identifying regions of a protein target that are more likely to accommodate plastically to regions of a drug molecule. Using these indicators we develop a second algorithm (PLASTIC) that provides a minimal manifold of possible conformations of a protein target for drug design, reducing the bias in structure-based drug design imparted by structures of enzymes co-crystallized with inhibitors.

Enzyme-catalyzed therapeutic agent (ECTA) design: activation of the antitumor ECTA compound NB1011 by thymidylate synthase

The in vivo administration of enzyme-inhibiting drugs for cancer and infectious disease often results in overexpression of the targeted enzyme. We have developed an enzyme-catalyzed therapeutic agent (ECTA) approach in which an enzyme overexpressed within the resistant cells is recruited as an intracellular catalyst for converting a relatively non-toxic substrate to a toxic product. We have investigated the potential of the ECTA approach to circumvent the thymidylate synthase (TS) overexpression-based resistance of tumor cells to conventional fluoropyrimidine [i.e. 5-fluorouracil (5-FU)] cancer chemotherapy. (E)-5-(2-Bromovinyl)-2'-deoxy-5'-uridyl phenyl L-methoxyalaninylphosphoramidate (NB1011) is a pronucleotide analogue of (E)-5-(2-bromovinyl)-2'-deoxyuridine (BVdU), an antiviral agent known to be a substrate for TS when in the 5'-monophosphorylated form. NB1011 was synthesized and found to be at least 10-fold more cytotoxic to 5-FU-resistant, TS-overexpressing colorectal tumor cells than to normal cells. This finding demonstrates that the ECTA approach to the design of novel chemotherapeutics results in compounds that are selectively cytotoxic to tumor cell lines that overexpress the target enzyme, TS, and therefore may be useful in the treatment of fluoropyrimidine-resistant cancer.

The crystal structure of thymidylate synthase from Pneumocystis carinii reveals a fungal insert important for drug design

Thymidylate synthase from Pneumocystis carinii (PcTS) is an especially important drug target, since P. carinii is a fungus that causes opportunistic pneumonia infections in immune-compromised patients and is among the leading causes of death of AIDS patients. Thymidylate synthase (TS) is the sole enzyme responsible for the de novo production of deoxythymidine monophosphate and hence is crucial for DNA replication in every organism. Inhibitors selective for P. carinii TS over human TS would be greatly beneficial in combating this disease. The crystal structure of TS from P. carinii bound to its substrate, dUMP, and a cofactor mimic, CB3717, was determined to 2.6 A resolution. A comparison with other species of TS shows that the volume of the closed PcTS active-site is 20 % larger than that of five other TS closed active-sites. A two-residue proline insert that is strictly conserved among all fungal species of TS, and a novel C-terminal closing interaction involving a P. carinii-specific tyrosine residue are primarily responsible for this increase in volume. The structure suggests several options for designing an inhibitor specific to PcTS and avoiding interactions with human TS. Taking advantage of the residue substitutions of P. carinii TS over human TS enables the design of a selective inhibitor. Additionally, the larger volume of the active-site of PcTS is an important advantage for designing de novo inhibitors that will exclude the human TS active-site through steric hindrance.

Crystal structures of rat thymidylate synthase inhibited by Tomudex, a potent anticancer drug

Two crystal structures of rat thymidylate synthase (TS) complexed with dUMP and the anticancer drug Tomudex (ZD1694) have been determined to resolutions of 3.3 and 2.6 A. Tomudex is one of several new antifolates targeted to TS and the first to be approved for clinical use. The structures represent the first views of any mammalian TS bound to ligands and suggest that the rat protein undergoes a ligand-induced conformational change similar to that of the Escherichia coli protein. Surprisingly, Tomudex does not induce the "closed" conformation in rat TS that is seen on binding to E. coli TS, resulting in inhibitor atoms that differ in position by more than 1.5 A. Several species-specific differences in sequence may be the reason for this. Phe 74 shifts to a new position in the rat complex and is in van der Waals contact with the inhibitor, while in the E. coli protein the equivalent amino acid (His 51) hydrogen bonds to the glutamate portion of the inhibitor. Amino acids Arg 101, Asn 106, and Met 305 make no contacts with the inhibitor in the open conformation, unlike the equivalent residues in the E. coli protein (Thr 78, Trp 83, and Val 262). dUMP binding is similar in both proteins, except that there is no covalent adduct to the active site cysteine (Cys 189) in the rat structures. Two insertions in the rat protein are clearly seen, but the N-termini (residues 1-20) and C-termini (residues 301-307) are disordered in both crystal forms.

Resistance to 5-fluorouracil

1. Primary and secondary resistance to the widely used antimetabolite 5-fluorouracil (5-FU) are common phenomena in cancer chemotherapy. Because 5-FU still remains the agent of choice in the treatment of, for example, colorectal cancer, circumvention of resistance is of vital importance. 2. Resistance to fluoropyrimidines is a multifactorial event, which includes transport mechanisms, metabolism, molecular mechanisms, protection from apoptosis, and resistance via cell cycle kinetics. To date, the prediction of primary resistance to 5-FU in the clinic is limited to few studies focusing mainly on the key enzyme thymidylate synthase. To gain a deeper insight into the key events responsible for 5-FU resistance in vivo, the evaluation of additional parameters such as other (fluoro)pyrimidine converting enzymes, the mutational status of regulators of apoptosis, and tumour angiogenesis is currently under investigation. 3. Most studies on the circumvention of fluoropyrimidine resistance refer to preclinical investigations and were rarely confirmed in clinical trials. Although our understanding of resistance to 5-FU leaves many open questions, the fundamental insights accomplished during the last years provide a rational understanding to exceed the bounds of the actual therapeutic schedules.