Design and synthesis of novel 6,7-imidazotetrahydroquinoline inhibitors of thymidylate synthase using iterative protein crystal structure analysis

Development and binding mode assessment of N-[4-[2-propyn-1-yl[(6S)-4,6,7,8-tetrahydro-2-(hydroxymethyl)-4-oxo-3H-cyclopenta[g]quinazolin-6-yl]amino]benzoyl]-l-γ-glutamyl-D-glutamic acid (BGC 945), a novel thymidylate synthase inhibitor that targets tumor cells

N-[4-[2-Propyn-1-yl[(6S)-4,6,7,8-tetrahydro-2-(hydroxymethyl)-4-oxo-3H-cyclopenta[g]quinazolin-6-yl]amino]benzoyl]-l-γ-glutamyl-d-glutamic acid 1 (BGC 945, now known as ONX 0801), is a small molecule thymidylate synthase (TS) inhibitor discovered at the Institute of Cancer Research in London. It is licensed by Onyx Pharmaceuticals and is in phase 1 clinical studies. It is a novel antifolate drug resembling TS inhibitors plevitrexed and raltitrexed that combines enzymatic inhibition of thymidylate synthase with α-folate receptor-mediated targeting of tumor cells. Thus, it has potential for efficacy with lower toxicity due to selective intracellular accumulation through α-folate receptor (α-FR) transport. The α-FR, a cell-surface receptor glycoprotein, which is overexpressed mainly in ovarian and lung cancer tumors, has an affinity for 1 similar to that for its natural ligand, folic acid. This study describes a novel synthesis of 1, an X-ray crystal structure of its complex with Escherichia coli TS and 2'-deoxyuridine-5'-monophosphate, and a model for a similar complex with human TS.

Synthesis and biological activity of splitomicin analogs targeted at human NAD(+)-dependent histone deacetylases (sirtuins)

Small molecules interfering with posttranslational modification of histones are of interest as tools to study epigenetic regulation of gene transcription. Specifically, drugs that interfere with histone deacetylation could be useful to induce differentiation, growth arrest as well as apoptotic cell death in tumor cells. One class of histone deacetylases is known as sirtuins some of which (Saccharomyces cerevisiae Sir2) are for example inhibited by the lactone splitomicin leading to telomeric silencing in yeast. However, splitomicin is only a micromolar inhibitor of yeast Sir2 and does not inhibit human subtypes and the lactone is prone to hydrolytic ring opening. In preliminary SAR-studies, splitomicin analogs lacking this hydrolytically labile ring were described as inactive while the naphthalene moiety could successfully be replaced by smaller aromatic rings in a fragment-like dihydrocoumarin. Here we report the synthesis and biological activity of a series of hydrolytically stable analogs with activity against human SIRT1 and 2. These comparatively small compounds characterized by high ligand efficiency are used as a starting point toward the development of specific inhibitors of histone deacetylases from the class of sirtuins.

A molecular docking strategy identifies Eosin B as a non-active site inhibitor of protozoal bifunctional thymidylate synthase-dihydrofolate reductase

Protozoal parasites are unusual in that their thymidylate synthase (TS) and dihydrofolate reductase (DHFR) enzymes exist on a single polypeptide. In an effort to probe the possibility of substrate channeling between the TS and DHFR active sites and to identify inhibitors specific for bifunctional TS-DHFR, we used molecular docking to screen for inhibitors targeting the shallow groove connecting the two active sites. Eosin B is a 100 microm non-active site inhibitor of Leishmania major TS-DHFR identified by molecular docking. Eosin B slows both the TS and DHFR reaction rates. When Arg-283, a key residue to which eosin B is predicted to bind, is mutated to glutamate, however, eosin B only minimally inhibits the TS-DHFR reaction. Additionally, eosin B was found to be a 180 microm inhibitor of Toxoplasma gondii in both biochemical and cell culture assays.

Structure-based studies on species-specific inhibition of thymidylate synthase

Thymidylate synthase (TS) is a well-recognized target for anticancer chemotherapy. Due to its key role in the sole de novo pathway for thymidylate synthesis and, hence, DNA synthesis, it is an essential enzyme in all life forms. As such, it has been recently recognized as a valuable new target against infectious diseases. There is also a pressing need for new antimicrobial agents that are able to target strains that are drug resistant toward currently used drugs. In this context, species specificity is of crucial importance to distinguish between the invading microorganism and the human host, yet thymidylate synthase is among the most highly conserved enzymes. We combine structure-based drug design with rapid synthetic techniques and mutagenesis, in an iterative fashion, to develop novel antifolates that are not derived from the substrate and cofactor, and to understand the molecular basis for the observed species specificity. The role of structural and computational studies in the discovery of nonanalog antifolate inhibitors of bacterial TS, naphthalein and dansyl derivatives, and in the understanding of their biological activity profile, are discussed.

Predicting and harnessing protein flexibility in the design of species-specific inhibitors of thymidylate synthase

BACKGROUND

Protein plasticity in response to ligand binding abrogates the notion of a rigid receptor site. Thus, computational docking alone misses important prospective drug design leads. Bacterial-specific inhibitors of an essential enzyme, thymidylate synthase (TS), were developed using a combination of computer-based screening followed by in-parallel synthetic elaboration and enzyme assay [Tondi et al. (1999) Chem. Biol. 6, 319-331]. Specificity was achieved through protein plasticity and despite the very high sequence conservation of the enzyme between species.

RESULTS

The most potent of the inhibitors synthesized, N,O-didansyl-L-tyrosine (DDT), binds to Lactobacillus casei TS (LcTS) with 35-fold higher affinity and to Escherichia coli TS (EcTS) with 24-fold higher affinity than to human TS (hTS). To reveal the molecular basis for this specificity, we have determined the crystal structure of EcTS complexed with DDT and 2'-deoxyuridine-5'-monophosphate (dUMP). The 2.0 A structure shows that DDT binds to EcTS in a conformation not predicted by molecular docking studies and substantially differently than other TS inhibitors. Binding of DDT is accompanied by large rearrangements of the protein both near and distal to the enzyme's active site with movement of C alpha carbons up to 6 A relative to other ternary complexes. This protein plasticity results in novel interactions with DDT including the formation of hydrogen bonds and van der Waals interactions to residues conserved in bacterial TS but not hTS and which are hypothesized to account for DDT's specificity. The conformation DDT adopts when bound to EcTS explains the activity of several other LcTS inhibitors synthesized in-parallel with DDT suggesting that DDT binds to the two enzymes in similar orientations.

CONCLUSIONS

Dramatic protein rearrangements involving both main and side chain atoms play an important role in the recognition of DDT by EcTS and highlight the importance of incorporating protein plasticity in drug design. The crystal structure of the EcTS/dUMP/DDT complex is a model system to develop more selective TS inhibitors aimed at pathogenic bacterial species. The crystal structure also suggests a general formula for identifying regions of TS and other enzymes that may be treated as flexible to aid in computational methods of drug discovery.

A novel approach to the design of inhibitors of human secreted phospholipase A2 based on native peptide inhibition

Human Type IIA secreted phospholipase A(2) (sPLA(2)-IIA) is an important modulator of cytokine-dependent inflammatory responses and a member of a growing superfamily of structurally related phospholipases. We have previously shown that sPLA(2)-IIA is inhibited by a pentapeptide sequence comprising residues 70-74 of the native sPLA(2)-IIA protein and that peptides derived from the equivalent region of different sPLA(2)-IIA species specifically inhibit the enzyme from which they are derived. We have now used an analogue screen of the human pentapeptide (70)FLSYK(74) in which side-chain residues were substituted, together with molecular docking approaches that modeled low-energy conformations of (70)FLSYK(74) bound to human sPLA(2)-IIA, to generate inhibitors with improved potency. Importantly, the modeling studies showed a close association between the NH(2) and COOH termini of the peptide, predicting significant enhancement of the potency of inhibition by cyclization. Cyclic compounds were synthesized and indeed showed 5-50-fold increased potency over the linear peptide in an Escherichia coli membrane assay. Furthermore, the potency of inhibition correlated with steady-state binding of the cyclic peptides to sPLA(2)-IIA as determined by surface plasmon resonance studies. Two potential peptide interaction sites were identified on sPLA(2)-IIA from the modeling studies, one in the NH(2)-terminal helix and the other in the beta-wing region, and in vitro association assays support the potential for interaction of the peptides with these sites. The inhibitors were effective at nanomolar concentrations in blocking sPLA(2)-IIA-mediated amplification of cytokine-induced prostaglandin synthesis in human rheumatoid synoviocytes in culture. These studies provide an example where native peptide sequences can be used for the development of potent and selective inhibitors of enzyme function.

Design, synthesis, and X-ray crystal structure of a potent dual inhibitor of thymidylate synthase and dihydrofolate reductase as an antitumor agent

A novel N-¿2-amino-4-methyl[(pyrrolo[2, 3-d]pyrimidin-5-yl)ethyl]benzoyl¿-L-glutamic acid (3a) was designed and synthesized as a potent dual inhibitor of thymidylate synthase (TS) and dihydrofolate reductase (DHFR) and as an antitumor agent. Compound 3b, the N7-benzylated analogue of 3a, was also synthesized as an antitumor agent. The synthesis of 3a was accomplished via a 12-step sequence which involved the synthesis of 2-amino-4-methylpyrrolo[2,3-d]pyrimidine (10) in 5 steps from 2-acetylbutyrolactone. Protection of the 2-amino group of 10 and regioselective iodination at the 5-position followed by palladium-catalyzed coupling afforded intermediate 14 which was converted to 3a by reduction and saponification. Similar synthetic methodology was used for 3b. X-ray crystal structure of the ternary complex of 3a, DHFR, and NADPH showed that the pyrrolo[2, 3-d]pyrimidine ring binds in a "2,4-diamino mode" in which the pyrrole nitrogen mimics the 4-amino moiety of 2,4-diaminopyrimidines. This is the first example of a classical pyrrolo[2,3-d]pyrimidine antifolate shown to have this alternate mode of binding to DHFR. Compounds 3a and 3b were more inhibitory than LY231514 against TS from Lactobacillus casei and Escherichia coli. Analogue 3a was also more inhibitory against DHFR from human, Toxoplasma gondii, and Pneumocystis carinii. Evaluation of 3a against methotrexate (MTX)-resistant cell lines with defined mechanisms indicated that cross-resistance of 3a was much lower than that of MTX. Metabolite protection studies and folylpoly-gamma-glutamate synthetase studies suggest that the antitumor activity of 3a against the growth of tumor cells in culture is a result of dual inhibition of TS and DHFR. Compound 3a inhibited the growth of CCRF-CEM and FaDu cells in culture at ED(50) values of 12.5 and 7.0 nM, respectively, and was more active against FaDu cells than MTX. In contrast, compound 3b was inactive against both cell lines. Compound 3a was evaluated in the National Cancer Institute in vitro preclinical antitumor screening program and afforded IG(50) values in the nanomolar range against a number of tumor cell lines.

Ab initio studies of some amino acid residue complexes with 4-mercaptopyridine as a model for thymitaq (AG337), an inhibitor of thymidylate synthase

The complex formed by isopentane, as a model for the isoleucine residue present in the wild-type thymidylate synthase, with 4-mercaptopyridine as a fragment of the thymidylate synthase inhibitor Thymitaq (AG337) is investigated with ab initio quantum chemical calculations at Hartree-Fock and MP2 levels, using the 3-21G* basis set. The binding energy is compared with the binding energies of 4-mercaptopyridine with amino acid residues found in mutant thymidylate synthase enzymes. As compared with isoleucine, alanine and glycine do not show binding, in agreement with enzyme-inhibition results.

Structure-based discovery and in-parallel optimization of novel competitive inhibitors of thymidylate synthase

BACKGROUND

The substrate sites of enzymes are attractive targets for structure-based inhibitor design. Two difficulties hinder efforts to discover and elaborate new (nonsubstrate-like) inhibitors for these sites. First, novel inhibitors often bind at nonsubstrate sites. Second, a novel scaffold introduces chemistry that is frequently unfamiliar, making synthetic elaboration challenging.

RESULTS

In an effort to discover and elaborate a novel scaffold for a substrate site, we combined structure-based screening with in-parallel synthetic elaboration. These techniques were used to find new inhibitors that bound to the folate site of Lactobacillus casei thymidylate synthase (LcTS), an enzyme that is a potential target for proliferative diseases, and is highly studied. The available chemicals directory was screened, using a molecular-docking computer program, for molecules that complemented the three-dimensional structure of this site. Five high-ranking compounds were selected for testing. Activity and docking studies led to a derivative of one of these, dansyltyrosine (Ki 65 microM). Using solid-phase in-parallel techniques 33 derivatives of this lead were synthesized and tested. These analogs are dissimilar to the substrate but bind competitively with it. The most active analog had a Ki of 1.3 microM. The tighter binding inhibitors were also the most specific for LcTS versus related enzymes.

CONCLUSIONS

TS can recognize inhibitors that are dissimilar to, but that bind competitively with, the folate substrate. Combining structure-based discovery with in-parallel synthetic techniques allowed the rapid elaboration of this series of compounds. More automated versions of this approach can be envisaged.

Structure-based design of inhibitors specific for bacterial thymidylate synthase

Thymidylate synthase is an attractive target for antiproliferative drug design because of its key role in the synthesis of DNA. As such, the enzyme has been widely targeted for anticancer applications. In principle, TS should also be a good target for drugs used to fight infectious disease. In practice, TS is highly conserved across species, and it has proven to be difficult to develop inhibitors that are selective for microbial TS enzymes over the human enzyme. Using the structure of TS from Lactobacillus casei in complex with the nonsubstrate analogue phenolphthalein, inhibitors were designed to take advantage of features of the bacterial enzyme that differ from those of the human enzyme. Upon synthesis and testing, these inhibitors were found to be up to 40-fold selective for the bacterial enzyme over the human enzyme. The crystal structures of two of these inhibitors in complex with TS suggested the design of further compounds. Subsequent synthesis and testing showed that these second-round compounds inhibit the bacterial enzyme at sub-micromolar concentrations, while the human enzyme was not inhibited at detectable levels (selectivities of 100-1000-fold or greater). Although these inhibitors share chemical similarities, X-ray crystal structures reveal that the analogues bind to the enzyme in substantially different orientations. Site-directed mutagenesis experiments suggest that the individual inhibitors may adopt multiple configurations in their complexes with TS.

Probing the folate-binding site of human thymidylate synthase by site-directed mutagenesis. Generation of mutants that confer resistance to raltitrexed, Thymitaq, and BW1843U89

Human thymidylate synthase (TS) contains three highly conserved residues Ile-108, Leu-221, and Phe-225 that have been suggested to be important for cofactor and antifolate binding. To elucidate the role of these residues and generate drug-resistant human TS mutants, 14 variants with multiple substitutions of these three hydrophobic residues were created by site-directed mutagenesis and transfected into mouse TS-negative cells for complementation assays and cytotoxicity studies, and the mutant proteins expressed and characterized. The I108A mutant confers resistance to raltitrexed and Thymitaq with respective IC50 values 54- and 80-fold greater than wild-type but less resistance to BW1843U89 (6-fold). The F225W mutant displays resistance to BW1843U89 (17-fold increase in IC50 values), but no resistance to raltitrexed and Thymitaq. It also confers 8-fold resistance to fluorodeoxyuridine. Both the kinetic characterization of the altered enzymes and formation of antifolate-resistant colonies in mouse bone marrow cells that express mutant TS are in accord with the IC50 values for cytotoxicity noted above. The human TS mutants (I108A and F225W), by virtue of their desirable properties, including good catalytic function and resistance to antifolate TS inhibitors, confirm the importance of amino acid residues Ile-108 and Phe-225 in the binding of folate and its analogues. These novel mutants may be useful for gene transfer experiments to protect hematopoietic progenitor cells from the toxic effects of these drugs.

Thymidylate synthase inhibition: a structure-based rationale for drug design

Thymidylate synthase (TS) is a very interesting target in antiproliferative diseases. Its inhibition causes thimineless death of the cells and compounds inhibiting TS are widely used in anticancer therapy. The classical antifolate TS inhibitors are structural analogs of the folate cofactor; they often share the same metabolic pathways and this causes the development of resistance inside the cells. A detailed analysis of the available x-ray crystal structures of the complexes of the enzyme with different substrates and inhibitors support the finding of a structural basis of their biological activity. TS inhibitors nonstructural analog of folate, non-analog antifolate inhibitors (NAAI), are welcome as a new interesting research topic. Among the most recent and interesting ones, compounds from Agouron related to the indole structure, are independent on the folate metabolism, highly active and specific for human TS. Other compounds, phthalein derivatives, can inhibit TS enzymes from various sources and show an interesting biological activity profile: they inhibit better bacterial and fungal TS than human TS. The x-ray crystal structures of some of these inhibitors with TS show that they bind in a different binding site from that of the classical folate TS inhibitors. This indicates a potential allosteric binding site useful for future drug discovery studies.

Intracellular location of thymidylate synthase and its state of phosphorylation

Thymidylate synthase (TS), an enzyme that is essential for DNA synthesis, was found to be associated mainly with the nucleolar region of H35 rat hepatoma cells, as determined both by immunogold electron microscopy and by autoradiography. In the latter case, the location of TS was established through the use of [6-3H]5-fluorodeoxyuridine, which forms a tight ternary complex of TS with 5-fluorodeoxyuridylate (FdUMP) and 5, 10-methylenetetrahydrofolylpolyglutamate within the cell. However, with H35 cells containing 50-100-fold greater amounts of TS than unmodified H35 cells, the enzyme, although still in the nucleus, was located primarily in the cytoplasm as shown by autoradiography and immunohistochemistry. In addition, TS was also present in mitochondrial extracts of both cell lines, as determined by enzyme activity measurements and by ternary complex formation with [32P]FdUMP and 5,10-methylenetetrahydrofolate. Another unique observation is that the enzyme appears to be a phosphoprotein, similar to that found for other proteins associated with cell division and signal transduction. The significance of these findings relative to the role of TS in cell division remains to be determined, but suggest that this enzyme's contribution to the cell cycle may be more complex than believed previously.

Cloning, expression and characterization of thymidylate synthase from Cryptococcus neoformans

The thymidylate synthase (TS)-encoding gene from Cryptococcus neoformans (Cn) has been isolated from cDNA and genomic libraries. The 1127-bp gene contains three introns and a 951-bp open reading frame encoding a 35,844-Da protein. The cDNA clones lack 324 bp of the 5' coding region of the gene. The complete coding sequence was assembled as an expression cassette in pUC19 using parts of the coding sequence from the cDNA and genomic DNA and completing the sequence using synthetic DNA. Production of active TS from Cn (CnTS) was first demonstrated by complementation of a thymine(Thy)-requiring Escherichia coli strain. The expression cassette was subsequently subcloned into the T7 polymerase vector pET15-b. In this construct, CnTS is produced as approximately 10% of the total soluble protein in E. coli. Homogeneous enzyme was obtained at a 36% yield after consecutive chromatography on DEAE-cellulose, Q-Sepharose, phenyl-Sepharose and Affi-Gel Blue. Steady-state kinetic analysis showed that the Km values for dUMP and CH2H4.folate were 2.7 +/- 0.5 microM and 38.2 +/- 2.5 microM, respectively, and the kcat was 5.1 s-1. The enzyme was stable upon storage at -80 degrees C in Tris.HCl pH 7.4 and thiol.

5-Fluoro-2'-deoxycytidine 5'-monophosphate is a mechanism-based inhibitor of thymidylate synthase

Thymidylate synthase (TS) is inhibited by 5-fluoro-2'-deoxycytidine 5'-monophosphate (FdCMP). From initial velocity measurements, the apparent Ki for the binary FdCMP-enzyme complex was about 20 microM. In the presence of 5,10-methylene-5,6,7,8-tetrahydrofolate (CH2H4folate), FdCMP causes a time-dependent inactivation of the enzyme and formation of a TS-FdCMP-CH2H4 folate complex. The ternary complex contains one mol of inhibitor per monomer of enzyme, and can be readily isolated on nitrocellulose filters. Dissociation of the ternary complex is quite slow (t1/2 approximately 16 h), and yields unchanged FdCMP. As with the corresponding complex formed with 5-fluoro-2'-deoxyuridine 5'-monophosphate (FdUMP), the TS-FdCMP-CH2H4 folate complex shows a differential absorbance maximum at 326 nm, and is stable to SDS-PAGE. Taken together, these results indicated that FdCMP is a slow, tight binding inhibitor of TS and has a mechanism of inhibition similar to that of FdUMP.

Flexibases: a way to enhance the use of molecular docking methods

Specially expanded databases containing three-dimensional structures are created to enhance the utility of docking methods to find new leads, i.e., active compounds of pharmacological interest. The expansion is based on the automatic generation of a set of maximally dissimilar conformations. The ligand receptor system of methotrexate and dihydrofolate reductase is used to demonstrate the feasibility of creating flexibases and their utility in docking studies.

On the use of LUDI to search the Fine Chemicals Directory for ligands of proteins of known three-dimensional structure

It is shown that the computer program LUDI can be used to search large database of three-dimensional structures for putative ligands of proteins with known 3D structure. As an example, a subset of approximately 30,000 small molecules (with less than 40 atoms and 0-2 rotatable bonds) from the Fine Chemicals Directory has been used in the search for possible novel ligands for four different proteins (trypsin, streptavidin, purine nucleoside phosphorylase and HIV protease). For trypsin and streptavidin, known ligands or substructures of known ligands are retrieved as top-scoring hits. In addition, a number of new interesting structures are found in all considered cases. Therefore, the method holds promise to retrieve automatically protein ligands from a 3D database if the 3D structure of the target protein is known.

Recent advances in drug design methods: where will they lead?

Drug design methods have made significant new advances over the last ten years, mainly in the areas of molecular modelling. In more recent times important developments in theory have led to a different type of modelling becoming possible, the so-called de novo or automated design algorithms. In this new method the programs perform much of the chemist's thinking, in finding appropriately sized chemical groups to fit into a target site. However this is a combinatoric problem which has no general analytical solution; it is ripe for optimization. Other advances, such as combinatorial chemical synthesis and screening, will dramatically influence the search for new lead structures for target sites, which at present are poorly understood. Already these methods are being applied to peptide libraries. Peptides do not make good drug compounds because of their poor bioavailability; further, their flexibility reduces their affinity. In some cases peptide backbones can be removed and replaced with rigid non-peptide scaffolds.

A novel computational tool for automated structure-based drug design

The computer program LUDI for automated structure-based drug design is described. The program constructs possible new ligands for a given protein of known three-dimensional structure. This novel approach is based upon rules about energetically favourable non-bonded contact geometries between functional groups of the protein and the ligand which are derived from a statistical analysis of crystal packings of organic molecules. In a first step small fragments are docked into the protein binding site in such a way that hydrogen bonds and ionic interactions can be formed with the protein and hydrophobic pockets are filled with lipophilic groups of the ligand. The program can then append further fragments onto a previously positioned fragment or onto an already existing ligand (e.g., a lead structure that one seeks to improve). It is also possible to link several fragments together by bridge fragments to form a complete molecule. All putative ligands retrieved or constructed by LUDI are scored. We use a simple scoring function that was fitted to experimentally determined binding constants of protein-ligand complexes. LUDI is a very fast program with typical execution times of 1-5 min on a work station and is therefore suitable for interactive usage.

Structure-based drug design: applications in immunopharmacology and immunosuppression

Structure-based drug design (SBDD) combines the power of many scientific disciplines, such as X-ray crystallography, nuclear magnetic resonance, medicinal chemistry, molecular modeling, biology, enzymology and biochemistry, in a functional paradigm of drug development. The current strength of SBDD lies in parlaying enzyme inhibitors into drugs, but a variety of technological advances over the past few years now makes it possible to address complex biological targets, such as those regulating immunosuppression and immunoactivation. Manual Navia and Debra Peattie discuss the SBDD paradigm and consider several of its achievements and challenges in immunopharmacology, particularly as these apply to the design of novel, potent immunosuppressants.

Structure-based drug design: applications in immunopharmacology and immunosuppression

Structure-based drug design (SBDD) combines the power of many scientific disciplines, such as X-ray crystallography, nuclear magnetic resonance, medicinal chemistry, molecular modeling, biology, enzymology and biochemistry, in a functional paradigm of drug development. The current strength of SBDD lies in parlaying enzyme inhibitors into drugs, but a variety of technological advances over the past few years now makes it possible to address complex biological targets, such as those regulating immunosuppression and immunoactivation. Manuel Navia and Debra Peattie discuss the SBDD paradigm and consider several of its achievements and challenges in immunopharmacology, particularly as these apply to the design of novel, potent immunosuppressants.