Structure-based design of inhibitors specific for bacterial thymidylate synthase

Synthesis, Characterization, and Biological Evaluation of Novel 7-Oxo-7H-thiazolo[3,2-b]-1,2,4-triazine-2-carboxylic Acid Derivatives

A series of novel 7-oxo-7H-thiazolo[3,2-b]-1,2,4-triazine-2-carboxylic acid derivatives was synthesized in good yields by a multi-step procedure that included the generation of the S-alkylated derivatives from 6-substituted arylmethyl-3-mercapto-1,2,4-triazin-5-ones with ethyl 2-chloroacetoacetate, intramolecular cyclization with microwave irradiation, hydrolysis and amidation. All of the target compounds were fully characterized through ¹H-NMR, ¹³C-NMR and HRMS spectra. The intramolecular cyclization occurred regioselectively at the N2-position of 1,2,4-triazine ring, which was confirmed by compound 3e using single-crystal X-ray diffraction analysis. The antibacterial and antitubercular activities of the target compounds were evaluated. Compared with Ciprofloxacin and Rifampicin, compounds 5d, 5f and 5g containing the terminal amide fragment exhibited broad spectrum antibacterial activity, and carboxylic acid derivatives or its corresponding ethyl esters had less effect on antibacterial properties. The most potent compound 5f also displayed excellent in vitro antitubercular activity against Mycobacterium smegmatis (minimum inhibitory concentration (MIC) = 50 μg/mL) and better growth inhibition activity of leucyl-tRNA synthetase (78.24 ± 4.05% at 15 μg/mL).

FLIP: An assisting software in structure based drug design using fingerprint of protein-ligand interaction profiles

With the growing number of labor-intensive data in the pharmaceutical industries and public domain for protein-ligand complexes, a significant challenge is still remaining in managing and leveraging this vast information. Here, a standalone application is presented for analysis, organization, and illustration of structural data and molecular interactions for exploiting 3D-structures into simple 1D fingerprints. The utility of the approach was shown in unraveling a feasible solution for post-processing of docking results in parallel with providing fruitful analysis for users in order to investigate molecular interactions. Remarkably, all interaction possibilities including (hydrogen bond, water-bridged, electrostatic, and hydrophobic as well as π- π and cation-π interactions) are supported both in the form of fingerprints and compelling reports. These investigations are mainly considered based on right orientation, location, and geometry of the interacting pairs rather than the acquisition of the energy terms. The reasonable efficiency of our application in different models was comparable to recent methods It is clearly presented that FLIP provides a faster way to generate usable fingerprints for ligand and protein binding modes. FLIP is free for academic use and is available at: http://zistrayan.com/development/download/flip/package.zip.

Crystal structures of nematode (parasitic T. spiralis and free living C. elegans), compared to mammalian, thymidylate synthases (TS). Molecular docking and molecular dynamics simulations in search for nematode-specific inhibitors of TS

Three crystal structures are presented of nematode thymidylate synthases (TS), including Caenorhabditis elegans (Ce) enzyme without ligands and its ternary complex with dUMP and Raltitrexed, and binary complex of Trichinella spiralis (Ts) enzyme with dUMP. In search of differences potentially relevant for the development of species-specific inhibitors of the nematode enzyme, a comparison was made of the present Ce and Ts enzyme structures, as well as binary complex of Ce enzyme with dUMP, with the corresponding mammalian (human, mouse and rat) enzyme crystal structures. To complement the comparison, tCONCOORD computations were performed to evaluate dynamic behaviors of mammalian and nematode TS structures. Finally, comparative molecular docking combined with molecular dynamics and free energy of binding calculations were carried out to search for ligands showing selective affinity to T. spiralis TS. Despite an overall strong similarity in structure and dynamics of nematode vs mammalian TSs, a pool of ligands demonstrating predictively a strong and selective binding to TsTS has been delimited. These compounds, the E63 family, locate in the dimerization interface of TsTS where they exert species-specific interactions with certain non-conserved residues, including hydrogen bonds with Thr174 and hydrophobic contacts with Phe192, Cys191 and Tyr152. The E63 family of ligands opens the possibility of future development of selective inhibitors of TsTS and effective agents against trichinellosis.

X-ray crystal structures of Enterococcus faecalis thymidylate synthase with folate binding site inhibitors

Infections caused by Enterococcus faecalis (Ef) represent nowadays a relevant health problem. We selected Thymidylate synthase (TS) from this organism as a potential specific target for antibacterial therapy. We have previously demonstrated that species-specific inhibition of the protein can be achieved despite the relatively high structural similarity among bacterial TSs and human TS. We had previously obtained the EfTS crystal structure of the protein in complex with the metabolite 5-formyl-tetrahydrofolate (5-FTHF) suggesting the protein role as metabolite reservoir; however, protein-inhibitors complexes were still missing. In the present work we identified some inhibitors bearing the phthalimidic core from our in-house library and we performed crystallographic screening towards EfTS. We obtained two X-ray crystallographic structures: the first with a weak phthalimidic inhibitor bound in one subunit and 5-hydroxymethylene-6-hydrofolic acid (5-HMHF) in the other subunit; a second X-ray structure complex with methotrexate. The structural information achieved confirm the role of EfTS as an enzyme involved in the folate pool system and provide a structural basis for structure-based drug design.

Perturbation of genome integrity to fight pathogenic microorganisms

BACKGROUND

Resistance against antibiotics is unfortunately still a major biomedical challenge for a wide range of pathogens responsible for potentially fatal diseases.

SCOPE OF REVIEW

In this study, we aim at providing a critical assessment of the recent advances in design and use of drugs targeting genome integrity by perturbation of thymidylate biosynthesis.

MAJOR CONCLUSION

We find that research efforts from several independent laboratories resulted in chemically highly distinct classes of inhibitors of key enzymes within the routes of thymidylate biosynthesis. The present article covers numerous studies describing perturbation of this metabolic pathway in some of the most challenging pathogens like Mycobacterium tuberculosis, Plasmodium falciparum, and Staphylococcus aureus.

GENERAL SIGNIFICANCE

Our comparative analysis allows a thorough summary of the current approaches to target thymidylate biosynthesis enzymes and also include an outlook suggesting novel ways of inhibitory strategies. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

2'-Deoxyuridine 5'-monophosphate substrate displacement in thymidylate synthase through 6-hydroxy-2H-naphtho[1,8-bc]furan-2-one derivatives

Thymidylate synthase (TS) is a target for antifolate-based chemotherapies of microbial and human diseases. Here, ligand-based, synthetic, and X-ray crystallography studies led to the discovery of 6-(3-cyanobenzoyloxy)-2-oxo-2H-naphto[1,8-bc]furan, a novel inhibitor with a Ki of 310 nM against Pneumocystis carinii TS. The X-ray ternary complex with Escherichia coli TS revealed, for the first time, displacement of the substrate toward the dimeric protein interface, thus providing new opportunities for further design of specific inhibitors of microbial pathogens.

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.

Perspective: Alchemical free energy calculations for drug discovery

Computational techniques see widespread use in pharmaceutical drug discovery, but typically prove unreliable in predicting trends in protein-ligand binding. Alchemical free energy calculations seek to change that by providing rigorous binding free energies from molecular simulations. Given adequate sampling and an accurate enough force field, these techniques yield accurate free energy estimates. Recent innovations in alchemical techniques have sparked a resurgence of interest in these calculations. Still, many obstacles stand in the way of their routine application in a drug discovery context, including the one we focus on here, sampling. Sampling of binding modes poses a particular challenge as binding modes are often separated by large energy barriers, leading to slow transitions. Binding modes are difficult to predict, and in some cases multiple binding modes may contribute to binding. In view of these hurdles, we present a framework for dealing carefully with uncertainty in binding mode or conformation in the context of free energy calculations. With careful sampling, free energy techniques show considerable promise for aiding drug discovery.

The structure of Enterococcus faecalis thymidylate synthase provides clues about folate bacterial metabolism

Drug resistance to therapeutic antibiotics poses a challenge to the identification of novel targets and drugs for the treatment of infectious diseases. Infections caused by Enterococcus faecalis are a major health problem. Thymidylate synthase (TS) from E. faecalis is a potential target for antibacterial therapy. The X-ray crystallographic structure of E. faecalis thymidylate synthase (EfTS), which was obtained as a native binary complex composed of EfTS and 5-formyltetrahydrofolate (5-FTHF), has been determined. The structure provides evidence that EfTS is a half-of-the-sites reactive enzyme, as 5-FTHF is bound to two of the four independent subunits present in the crystal asymmetric unit. 5-FTHF is a metabolite of the one-carbon transfer reaction catalysed by 5-formyltetrahydrofolate cyclo-ligase. Kinetic studies show that 5-FTHF is a weak inhibitor of EfTS, suggesting that the EfTS-5-FTHF complex may function as a source of folates and/or may regulate one-carbon metabolism. The structure represents the first example of endogenous 5-FTHF bound to a protein involved in folate metabolism.

Folate metabolic pathways in Leishmania

Trypanosomatid parasitic protozoans of the genus Leishmania are autotrophic for both folate and unconjugated pteridines. Leishmania salvage these metabolites from their mammalian hosts and insect vectors through multiple transporters. Within the parasite, folates are reduced by a bifunctional DHFR (dihydrofolate reductase)-TS (thymidylate synthase) and by a novel PTR1 (pteridine reductase 1), which reduces both folates and unconjugated pteridines. PTR1 can act as a metabolic bypass of DHFR inhibition, reducing the effectiveness of existing antifolate drugs. Leishmania possess a reduced set of folate-dependent metabolic reactions and can salvage many of the key products of folate metabolism from their hosts. For example, they lack purine synthesis, which normally requires 10-formyltetrahydrofolate, and instead rely on a network of purine salvage enzymes. Leishmania elaborate at least three pathways for the synthesis of the key metabolite 5,10-methylene-tetrahydrofolate, required for the synthesis of thymidylate, and for 10-formyltetrahydrofolate, whose presumptive function is for methionyl-tRNAMet formylation required for mitochondrial protein synthesis. Genetic studies have shown that the synthesis of methionine using 5-methyltetrahydrofolate is dispensable, as is the activity of the glycine cleavage complex, probably due to redundancy with serine hydroxymethyltransferase. Although not always essential, the loss of several folate metabolic enzymes results in attenuation or loss of virulence in animal models, and a null DHFR-TS mutant has been used to induce protective immunity. The folate metabolic pathway provides numerous opportunities for targeted chemotherapy, with strong potential for 'repurposing' of compounds developed originally for treatment of human cancers or other infectious agents.

Identification of the binding modes of N-phenylphthalimides inhibiting bacterial thymidylate synthase through X-ray crystallography screening

To identify specific bacterial thymidylate synthase (TS) inhibitors, we exploited phenolphthalein (PTH), which inhibits both bacterial and human enzymes. The X-ray crystal structure of Lactobacillus casei TS (LcTS) that binds PTH showed multiple binding modes of the inhibitor, which prevented a classical structure-based drug design approach. To overcome this issue, we synthesized two phthalimidic libraries that were tested against TS enzymes and then we performed X-ray crystallographic screening of the active compounds. Compounds 6A, 8A, and 12A showed 40-fold higher affinity for bacterial TS than human TS. The X-ray crystallographic screening characterized the binding mode of six inhibitors in complexes with LcTS. Of these, 20A, 23A, and 24A showed a common unique binding mode, whereas 8A showed a different, unique binding mode. A comparative analysis of the LcTS X-ray complexes that were obtained with the pathogenic TS enabled the selection of compounds 8A and 23A as specific compounds and starting points to be exploited for the specific inhibition of pathogen enzymes.

Antifolate agents: a patent review (2006 - 2010)

INTRODUCTION

For > 50 years, drugs targeting the folate pathway have significantly impacted disease treatment as anticancer, antimicrobial and immunomodulatory agents. The discovery of novel antifolate agents with improved properties and superior activities remains an attractive strategy, both in academia and the pharmaceutical industry.

AREAS COVERED

This review surveys the patent literature from 2006 to 2010 for small molecule inhibitors of enzymatic targets in the folate biosynthetic pathway.

EXPERT OPINION

The pursuit of antifolates as anticancer and antimicrobial agents continues to be an active area of research. New patent disclosures reveal novel antifolate scaffolds, antifolates with improved drug-like properties and new strategies to effectively target cancer cells. The continued use of high resolution structural information has guided the discovery of several compounds. Owing to the need for high levels of potency and selectivity, especially in targeting pathogenic species, the use of high resolution crystal structures remains an important tool to guide the design of novel antifolates. Interestingly, the patents disclosing novel compounds were ones where X-ray crystallography was an integral component of the design process. Finally, a variety of new structures have been reported that may play an important role in the future development of therapeutic antifolates.

Mapping of ligand-binding cavities in proteins

The complex interactions between proteins and small organic molecules (ligands) are intensively studied because they play key roles in biological processes and drug activities. Here, we present a novel approach to characterize and map the ligand-binding cavities of proteins without direct geometric comparison of structures, based on Principal Component Analysis of cavity properties (related mainly to size, polarity, and charge). This approach can provide valuable information on the similarities and dissimilarities, of binding cavities due to mutations, between-species differences and flexibility upon ligand-binding. The presented results show that information on ligand-binding cavity variations can complement information on protein similarity obtained from sequence comparisons. The predictive aspect of the method is exemplified by successful predictions of serine proteases that were not included in the model construction. The presented strategy to compare ligand-binding cavities of related and unrelated proteins has many potential applications within protein and medicinal chemistry, for example in the characterization and mapping of "orphan structures", selection of protein structures for docking studies in structure-based design, and identification of proteins for selectivity screens in drug design programs.

Binding of small-molecule ligands to proteins: "what you see" is not always "what you get"

We review insights from computational studies of affinities of ligands binding to proteins. The power of structural biology is in translating knowledge of protein structures into insights about their forces, binding, and mechanisms. However, the complementary power of computer modeling is in showing "the rest of the story" (i.e., how motions and ensembles and alternative conformers and the entropies and forces that cannot be seen in single molecular structures also contribute to binding affinities). Upon binding to a protein, a ligand can bind in multiple orientations; the protein or ligand can be deformed by the binding event; waters, ions, or cofactors can have unexpected involvement; and conformational or solvation entropies can sometimes play large and otherwise unpredictable roles. Computer modeling is helping to elucidate these factors.

Kinetics and ligand-binding preferences of Mycobacterium tuberculosis thymidylate synthases, ThyA and ThyX

Background

Mycobacterium tuberculosis kills approximately 2 million people each year and presents an urgent need to identify new targets and new antitubercular drugs. Thymidylate synthase (TS) enzymes from other species offer good targets for drug development and the M. tuberculosis genome contains two putative TS enzymes, a conventional ThyA and a flavin-based ThyX. In M. tuberculosis, both TS enzymes have been implicated as essential for growth, either based on drug-resistance studies or genome-wide mutagenesis screens. To facilitate future small molecule inhibitors against these proteins, a detailed enzymatic characterization was necessary.

Methodology/Principal Findings

After cloning, overexpression, and purification, the thymidylate-synthesizing ability of ThyA and ThyX gene products were directly confirmed by HPLC analysis of reaction products and substrate saturation kinetics were established. 5-Fluoro-2′-deoxyuridine 5′-monophosphate (FdUMP) was a potent inhibitor of both ThyA and ThyX, offering important clues to double-targeting strategies. In contrast, the folate-based 1843U89 was a potent inhibitor of ThyA but not ThyX suggesting that it should be possible to find ThyX-specific antifolates. A turnover-dependent kinetic assay, combined with the active-site titration approach of Ackermann and Potter, revealed that both M. tuberculosis enzymes had very low k_cat values. One possible explanation for the low catalytic activity of M. tuberculosis ThyX is that its true biological substrates remain to be identified. Alternatively, this slow-growing pathogen, with low demands for TMP, may have evolved to down-regulate TS activities by altering the turnover rate of individual enzyme molecules, perhaps to preserve total protein quantities for other purposes. In many organisms, TS is often used as a part of larger complexes of macromolecules that control replication and DNA repair.

Conclusions/Significance

Thus, the present enzymatic characterization of ThyA and ThyX from M. tuberculosis provides a framework for future development of cell-active inhibitors and the biological roles of these TS enzymes in M. tuberculosis.

Role of an invariant lysine residue in folate binding on Escherichia coli thymidylate synthase: calorimetric and crystallographic analysis of the K48Q mutant

Thymidylate synthase (TS) catalyzes the reductive methylation of deoxyuridine monophosphate (dUMP) using methylene tetrahydrofolate (CH(2)THF) as cofactor, the glutamate tail of which forms a water-mediated hydrogen bond with an invariant lysine residue of this enzyme. To understand the role of this interaction, we studied the K48Q mutant of Escherichia coli TS using structural and biophysical methods. The k(cat) of the K48Q mutant was 430-fold lower than wild-type TS in activity, while the K(m) for the (R)-stereoisomer of CH(2)THF was 300 microM, about 30-fold larger than K(m) from the wild-type TS. Affinity constants were determined using isothermal titration calorimetry, which showed that binding was reduced by one order of magnitude for folate-like TS inhibitors, such as propargyl-dideazafolate (PDDF) or compounds that distort the TS active site like BW1843U89 (U89). The crystal structure of the K48Q-dUMP complex revealed that dUMP binding is not impaired in the mutant, and that U89 in a ternary complex of K48Q-nucleotide-U89 was bound in the active site with subtle differences relative to comparable wild-type complexes. PDDF failed to form ternary complexes with K48Q and dUMP. Thermodynamic data correlated with the structural determinations, since PDDF binding was dominated by enthalpic effects while U89 had an important entropic component. In conclusion, K48 is critical for catalysis since it leads to a productive CH(2)THF binding, while mutation at this residue does not affect much the binding of inhibitors that do not make contact with this group.

Thymidyl biosynthesis enzymes as antibiotic targets

The two long-known "classical" enzymes of uridyl-5-methylation, thymidylate synthase and ribothymidyl synthase, have been joined by two alternative methylation enzymes, flavin-dependent thymidylate synthase and folate-dependent ribothymidyl synthase. These two newly discovered enzymes have much in common: both contain flavin cofactors, utilize methylenetetrahydrofolate as a source of methyl group, and perform thymidylate synthesis via chemical pathways distinct from those of their classic counterparts. Several severe human pathogens (e.g., typhus, anthrax, tuberculosis, and more) depend on these "alternative" enzymes for reproduction. These and other distinctive properties make the alternative enzymes and their corresponding genes appealing targets for new antibiotics.

Antibacterial Agent Discovery Using Thymidylate Synthase Biolibrary Screening

Thymidylate synthase (TS, ThyA) catalyzes the reductive methylation of 2'-deoxyuridine 5'-monophosphate to 2'-deoxythymidine 5'-monophosphate, an essential precursor for DNA synthesis. A specific inhibition of this enzyme induces bacterial cell death. As a second round lead optimization design, new 1,2-naphthalein derivatives have been synthesized and tested against a TS-based biolibrary, including human thymidylate synthase (hTS). Docking studies have been performed to rationalize the experimentally observed affinity profiles of 1,2-naphthalein compounds toward Lactobacillus casei TS and hTS. The best TS inhibitors have been tested against a number of clinical isolates of Gram-positive-resistant bacterial strains. Compound 3,3-bis(3,5-dibromo-4-hydroxyphenyl)-1H,3H-naphtho[1,2-c]furan-1-one (5) showed significant antibacterial activity, no in vitro toxicity, and dose-response effects against Staphylococcus epidermidis (MIC=0.5-2.5 microg/mL) clinical isolate strains, which are resistant to at least 17 of the best known antibacterial agents, including vancomycin. So far this compound can be regarded as a leading antibacterial agent.

Structures of Leishmania major pteridine reductase complexes reveal the active site features important for ligand binding and to guide inhibitor design

Pteridine reductase (PTR1) is an NADPH-dependent short-chain reductase found in parasitic trypanosomatid protozoans. The enzyme participates in the salvage of pterins and represents a target for the development of improved therapies for infections caused by these parasites. A series of crystallographic analyses of Leishmania major PTR1 are reported. Structures of the enzyme in a binary complex with the cofactor NADPH, and ternary complexes with cofactor and biopterin, 5,6-dihydrobiopterin, and 5,6,7,8-tetrahydrobiopterin reveal that PTR1 does not undergo any major conformational changes to accomplish binding and processing of substrates, and confirm that these molecules bind in a single orientation at the catalytic center suitable for two distinct reductions. Ternary complexes with cofactor and CB3717 and trimethoprim (TOP), potent inhibitors of thymidylate synthase and dihydrofolate reductase, respectively, have been characterized. The structure with CB3717 reveals that the quinazoline moiety binds in similar fashion to the pterin substrates/products and dominates interactions with the enzyme. In the complex with TOP, steric restrictions enforced on the trimethoxyphenyl substituent prevent the 2,4-diaminopyrimidine moiety from adopting the pterin mode of binding observed in dihydrofolate reductase, and explain the inhibition properties of a range of pyrimidine derivates. The molecular detail provided by these complex structures identifies the important interactions necessary to assist the structure-based development of novel enzyme inhibitors of potential therapeutic value.

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.

Identification and Characterization of Small Molecule Inhibitors of the Calcium-Dependent S100B−p53 Tumor Suppressor Interaction

The binding of S100B to p53 down-regulates wild-type p53 tumor suppressor activity in cancer cells such as malignant melanoma, so a search for small molecules that bind S100B and prevent S100B-p53 complex formation was undertaken. Chemical databases were computationally searched for potential inhibitors of S100B, and 60 compounds were selected for testing on the basis of energy scoring, commercial availability, and chemical similarity clustering. Seven of these compounds bound to S100B as determined by steady state fluorescence spectroscopy (1.0 microM < or = K(D) < or = 120 microM) and five inhibited the growth of primary malignant melanoma cells (C8146A) at comparable concentrations (1.0 microM < or = IC(50) < or = 50 microM). Additionally, saturation transfer difference (STD) NMR experiments confirmed binding and qualitatively identified protons from the small molecule at the small molecule-S100B interface. Heteronuclear single quantum coherence (HSQC) NMR titrations indicate that these compounds interact with the p53 binding site on S100B. An NMR-docked model of one such inhibitor, pentamidine, bound to Ca(2+)-loaded S100B was calculated using intermolecular NOE data between S100B and the drug, and indicates that pentamidine binds into the p53 binding site on S100B defined by helices 3 and 4 and loop 2 (termed the hinge region).

Inhibitor Specificity via Protein Dynamics

Structure-based drug design of species-specific inhibitors generally exploits structural differences in proteins from different organisms. Here, we demonstrate how achieving specificity can be aided by targeting differences in the dynamics of proteins. Thymidylate synthase (TS) is a good target for anticancer agents and a potential target for antibacterial agents. Most inhibitors are folate-analogs that bind at the folate binding site and are not species specific. In contrast, alpha156 is not a folate-analog and is specific for bacterial TS; it has been shown crystallographically to bind in a nonconserved binding site. Docking calculations and crystal structure-based estimation of the essential dynamics of TSs from five different species show that differences in the dynamics of TSs make the active site more accessible to alpha156 in the prokaryotic than in the eukaryotic TSs and thereby enhance the specificity of alpha156.

Implications of protein flexibility for drug discovery

Proteins are in constant motion between different conformational states with similar energies. This has often been ignored in drug design. However, protein flexibility is fundamental to understanding the ways in which drugs exert biological effects, their binding-site location, binding orientation, binding kinetics, metabolism and transport. Protein flexibility allows increased affinity to be achieved between a drug and its target. This is crucial, because the lipophilicity and number of polar interactions allowed for an oral drug is limited by absorption, distribution, metabolism and toxicology considerations.

ortho-Halogen naphthaleins as specific inhibitors of Lactobacillus casei thymidylate synthase. Conformational properties and biological activity

Thymidylate synthase (TS) (EC 2.1.1.45), an enzyme involved in the DNA synthesis of both prokaryotic and eukaryotic cells, is a potential target for the development of anticancer and antinfective agents. Recently, we described a series of phthalein and naphthalein derivatives as TS inhibitors. These compounds have structures unrelated to the folate (Non-Analogue Antifolate Inhibitors, NAAIs) and were selective for the bacterial versus the human TS (hTS). In particular, halogen-substituted molecules were the most interesting. In the present paper the halogen derivatives of variously substituted 3,3-bis(4-hydroxyphenyl)-1H,3H-naphtho[2,3-c]furan-1-one (1-5) and 3,3-bis(4-hydroxyphenyl)-1H,3H-naphtho[1,8-c,d]pyran-1-one (6-14) were synthesized to investigate the biological effect of halogen substitution on the inhibition and selectivity for the TS enzymes. Conformational properties of the naphthalein series were explored in order to highlight possible differences between molecules that show species-specific biological profile with respect to non species-specific ones. With this aim, the conformational properties of the synthesized compounds were investigated by NMR, in various solvents and at different temperatures, and by computational analysis. The apparent inhibition constants (K(i)) for Lactobacillus casei TS (LcTS) were found to range from 0.7 to 7.0 microM, with the exception of the weakly active iodo-derivatives (4, 10, 13); all] the compounds were poorly active against hTS. The di-halogenated compounds 7, 8, 14 showed the highest specificity towards LcTS, their specificity index (SI) ranging between 40 and >558. The di-halogenated 1,8-naphthalein derivatives (7-10) exhibited different conformational properties with respect to the tetra-haloderivatives. Though a clear explanation for the observed specificity by means of conformational analysis is difficult to find, some interesting conformational effects are discussed in the context of selective recognition of the compounds investigated by the LcTS enzyme.

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.

Object-oriented approach to drug design enabled by NMR SOLVE: first real-time structural tool for characterizing protein-ligand interactions

As a result of genomics efforts, the number of protein drug targets is expected to increase by an order of magnitude. Functional genomics efforts are identifying these targets, while structural genomics efforts are determining structures for many of them. However, there is a significant gap in going from structural information for a protein target to a high affinity (K(d) < 100 nM) inhibitor, and the problem is multiplied by the sheer number of new targets now available. nature frequently designs proteins in classes that are related by the reuse, through gene duplication events, of cofactor binding domains. This reuse of functional domains is an efficient way to build related proteins in that it is object-oriented. There is a growing realization that the most efficient drug design strategies for attacking the mass of targets coming from genomics efforts will be systems-based approaches that attack groups of related proteins in parallel. We propose that the most effective drug design strategy will be one that parallels the object-oriented manner by which nature designed the gene families themselves. IOPE (Integrated Object-Oriented PharmacoEngineering) is such an approach. It is a three-step technology to build focused combinatorial libraries of potential inhibitors for major families and sub-families of enzymes, using cogent NMR data derived from representatives of these protein families. The NMR SOLVE (Structurally Oriented Library Valency Engineering) data used to design these libraries are gathered in days, and data can be obtained for large proteins (> 170 kDa). Furthermore, the process is fully object-oriented in that once a given bi-ligand is identified for a target, potency is retained if different cofactor mimics are swapped. This gives the drug design process maximum flexibility, allowing for the more facile transition from in vitro potency to in vivo efficacy.

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.

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.

Crystal structure of a deletion mutant of human thymidylate synthase Delta (7-29) and its ternary complex with Tomudex and dUMP

The crystal structures of a deletion mutant of human thymidylate synthase (TS) and its ternary complex with dUMP and Tomudex have been determined at 2.0 A and 2.5 A resolution, respectively. The mutant TS, which lacks 23 residues near the amino terminus, is as active as the wild-type enzyme. The ternary complex is observed in the open conformation, similar to that of the free enzyme and to that of the ternary complex of rat TS with the same ligands. This is in contrast to Escherichia coli TS, where the ternary complex with Tomudex and dUMP is observed in the closed conformation. While the ligands interact with each other in identical fashion regardless of the enzyme conformation, they are displaced by about 1.0 A away from the catalytic cysteine in the open conformation. As a result, the covalent bond between the catalytic cysteine sulfhydryl and the base of dUMP, which is the first step in the reaction mechanism of TS and is observed in all ternary complexes of the E. coli enzyme, is not formed. This displacement results from differences in the interactions between Tomudex and the protein that are caused by differences in the environment of the glutamyl tail of the Tomudex molecule. Despite the absence of the closed conformation, Tomudex inhibits human TS ten-fold more strongly than E. coli TS. These results suggest that formation of a covalent bond between the catalytic cysteine and the substrate dUMP is not required for effective inhibition of human TS by cofactor analogs and could have implications for drug design by eliminating this as a condition for lead compounds.

Induction of chromosome aberrations in vitro by phenolphthalein: mechanistic studies

Phenolphthalein induces tumors in rodents but because it is negative in assays for mutation in Salmonella and in mammalian cells, for DNA adducts and for DNA strand breaks, its primary mechanism does not seem to be DNA damage. Chromosome aberration (Ab) induction by phenolphthalein in vitro is associated with marked cytotoxicity. At very high doses, phenolphthalein induces weak increases in micronuclei (MN) in mouse bone marrow; a larger response is seen with chronic treatment. All this suggests genotoxicity is a secondary effect that may not occur at lower doses. In heterozygous TSG-p53((R)) mice, phenolphthalein induces lymphomas and also MN, many with kinetochores (K), implying chromosome loss. Induction of aneuploidy would be compatible with the loss of the normal p53 gene seen in the lymphomas. Here we address some of the postulated mechanisms of genotoxicity in vitro, including metabolic activation, inhibition of thymidylate synthetase, cytotoxicity, oxidative stress, DNA damage and aneuploidy. We show clearly that phenolphthalein does not require metabolic activation by S9 to induce Abs. Inhibition of thymidylate synthetase is an unlikely mechanism, since thymidine did not prevent Ab induction by phenolphthalein. Phenolphthalein dramatically inhibited DNA synthesis, in common with many non-DNA reactive chemicals that induce Abs at cytotoxic doses. Phenolphthalein strongly enhances levels of intracellular oxygen radicals (ROS). The radical scavenger DMSO suppresses phenolphthalein-induced toxicity and Abs whereas H(2)O(2) potentiates them, suggesting a role for peroxidative activation. Phenolphthalein did not produce DNA strand breaks in rat hepatocytes or DNA adducts in Chinese hamster ovary (CHO) cells. All the evidence points to an indirect mechanism for Abs that is unlikely to operate at low doses of phenolphthalein. We also found that phenolphthalein induces mitotic abnormalities and MN with kinetochores in vitro. These are also enhanced by H(2)O(2) and suppressed by DMSO. Our findings suggest that induction of Abs in vitro is a high-dose effect in oxidatively stressed cells and may thus have a threshold. There may be more than one mechanism operating in vitro and in vivo, possibly indirect genotoxicity at high doses and also chromosome loss, both of which would likely have a threshold.

Site-directed ligand discovery

We report a strategy (called "tethering") to discover low molecular weight ligands ( approximately 250 Da) that bind weakly to targeted sites on proteins through an intermediary disulfide tether. A native or engineered cysteine in a protein is allowed to react reversibly with a small library of disulfide-containing molecules ( approximately 1,200 compounds) at concentrations typically used in drug screening (10 to 200 microM). The cysteine-captured ligands, which are readily identified by MS, are among the most stable complexes, even though in the absence of the covalent tether the ligands may bind very weakly. This method was applied to generate a potent inhibitor for thymidylate synthase, an essential enzyme in pyrimidine metabolism with therapeutic applications in cancer and infectious diseases. The affinity of the untethered ligand (K(i) approximately 1 mM) was improved 3,000-fold by synthesis of a small set of analogs with the aid of crystallographic structures of the tethered complex. Such site-directed ligand discovery allows one to nucleate drug design from a spatially targeted lead fragment.

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.

Phthalein derivatives as a new tool for selectivity in thymidylate synthase inhibition

A new set of phthalein derivatives stemming from the lead compound, phenolphthalein, were designed to specifically complement structural features of a bacterial form of thymidylate synthase (Lactobacillus casei, LcTS) versus the human TS (hTS) enzyme. The new compounds were screened for their activity and their specificity against TS enzymes from different species, namely, L. casei (LcTS), Pneumocystis carinii (PcTS), Cryptococcus neoformans (CnTS), and human thymidylate synthase (hTS). Apparent inhibition constants (Ki) for all the compounds against LcTS were determined, and inhibition factors (IF, ratio between the initial rates of the enzymatic reaction in the presence and absence of each inhibitor) against each of the four TS species were measured. A strong correlation was found between the two activity parameters, IF and Ki, and therefore the simpler IF was used as a screening factor in order to accelerate biological evaluation. Compounds 5b, 5c, 5ba, and 6bc showed substantial inhibition of LcTS while remaining largely inactive against hTS, illustrating for the first time remarkable species specificity among TSs. Due to sequence homology between the enzymes, several compounds also showed high activity and specificity for CnTS. In particular, 3-hydroxy-3-(3-chloro-4-hydroxyphenyl)-6-nitro-1H, 3H-naphtho[1,8-c,d]pyran-1-one (6bc) showed an IF < 0.04 for CnTS (Ki = 0.45 microM) while remaining inactive in the hTS assay at the maximum solubility concentration of the compound (200 microM). In cell culture assays most of the compounds were found to be noncytotoxic to human cell lines but were cytotoxic against several species of Gram-positive bacteria. These results are consistent with the enzymatic assays. Intriguingly, several compounds also had selective activity against Cr. neoformans in cell culture assay. In general, the most active and selective compounds against the Gram-positive bacteria were those designed and found in the enzyme assay to be specific for LcTS versus hTS. The original lead compound was least selective against most of the cell lines tested. To our knowledge these compounds are the first TS inhibitors selective for bacterial TS with respect to hTS.

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.