Computer-assisted combinatorial design of bicyclic thymidine analogs as inhibitors of Mycobacterium tuberculosis thymidine monophosphate kinase

Thymidine monophosphate kinase (TMPKmt) is an essential enzyme for nucleotide metabolism in Mycobacterium tuberculosis, and thus an attractive target for novel antituberculosis agents. In this work, we have explored the chemical space around the 2′,3′-bicyclic thymidine nucleus by designing and in silico screening of a virtual focused library selected via structure based methods to identify more potent analogs endowed with favorable ADME-related properties. In all the library members we have exchanged the ribose ring of the template with a cyclopentane moiety that is less prone to enzymatic degradation. In addition, we have replaced the six-membered 2′,3′-ring by a number of five-membered and six-membered heterocyclic rings containing alternative proton donor and acceptor groups, to exploit the interaction with the carboxylate groups of Asp9 and Asp163 as well as with several cationic residues present in the vicinity of the TMPKmt binding site. The three-dimensional structure of the TMPKmt complexed with 5-hydroxymethyl-dUMP, an analog of dTMP, was employed to develop a QSAR model, to parameterize a scoring function specific for the TMPKmt target and to select analogues which display the highest predicted binding to the target. As a result, we identified a small highly focused combinatorial subset of bicyclic thymidine analogues as virtual hits that are predicted to inhibit the mycobacterial TMPK in the submicromolar concentration range and to display favorable ADME-related properties.Graphical abstract

[1]  Ray A. Jarvis,et al.  Clustering Using a Similarity Measure Based on Shared Near Neighbors , 1973, IEEE Transactions on Computers.

[2]  S. Miertus,et al.  Design and in silico screening of combinatorial library of antimalarial analogs of triclosan inhibiting Plasmodium falciparum enoyl-acyl carrier protein reductase. , 2009, European journal of medicinal chemistry.

[3]  Thierry Langer,et al.  Impact of Scoring Functions on Enrichment in Docking-Based Virtual Screening: An Application Study on Renin Inhibitors , 2004, J. Chem. Inf. Model..

[4]  R. Sclafani,et al.  Yeast gene CDC8 encodes thymidylate kinase and is complemented by herpes thymidine kinase gene TK. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

[6]  I. Muegge Effect of ligand volume correction on PMF scoring , 2001, J. Comput. Chem..

[7]  Piet Herdewijn,et al.  3'-C-branched-chain-substituted nucleosides and nucleotides as potent inhibitors of Mycobacterium tuberculosis thymidine monophosphate kinase. , 2003, Journal of medicinal chemistry.

[8]  H. Munier-Lehmann,et al.  Discovery of bicyclic thymidine analogues as selective and high-affinity inhibitors of Mycobacterium tuberculosis thymidine monophosphate kinase. , 2004, Journal of medicinal chemistry.

[9]  S. Miertus,et al.  Design of peptidomimetic inhibitors of aspartic protease of HIV-1 containing -Phe Psi Pro- core and displaying favourable ADME-related properties. , 2008, Journal of molecular graphics & modelling.

[10]  Ming-Jing Hwang,et al.  Derivation of class II force fields. I. Methodology and quantum force field for the alkyl functional group and alkane molecules , 1994, J. Comput. Chem..

[11]  S. Miertus,et al.  Combinatorial design of nonsymmetrical cyclic urea inhibitors of aspartic protease of HIV-1. , 2005, Bioorganic & medicinal chemistry.

[12]  G. Labesse,et al.  Thymidylate kinase of Mycobacterium tuberculosis: A chimera sharing properties common to eukaryotic and bacterial enzymes , 2001, Protein science : a publication of the Protein Society.

[13]  H. Munier-Lehmann,et al.  Thymidine and thymidine-5'-O-monophosphate analogues as inhibitors of Mycobacterium tuberculosis thymidylate kinase. , 2003, Bioorganic & medicinal chemistry letters.

[14]  Maria Kontoyianni,et al.  Evaluation of library ranking efficacy in virtual screening , 2005, J. Comput. Chem..

[15]  Y. Martin,et al.  A general and fast scoring function for protein-ligand interactions: a simplified potential approach. , 1999, Journal of medicinal chemistry.

[16]  W. L. Jorgensen,et al.  Prediction of drug solubility from structure. , 2002, Advanced drug delivery reviews.

[17]  E. De Clercq,et al.  Novel bicyclic sugar modified nucleosides: synthesis, conformational analysis and antiviral evaluation. , 2004, Bioorganic & medicinal chemistry.

[18]  M R Lee,et al.  Use of MM‐PB/SA in estimating the free energies of proteins: Application to native, intermediates, and unfolded villin headpiece , 2000, Proteins.

[19]  J. Balzarini,et al.  Thymidine phosphorylase inhibitors: recent developments and potential therapeutic applications. , 2005, Mini reviews in medicinal chemistry.

[20]  Jeremy C. Smith,et al.  Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal‐gas entropy correction? , 2002, J. Comput. Chem..

[21]  H. Munier-Lehmann,et al.  Synthesis and evaluation of thymidine-5'-O-monophosphate analogues as inhibitors of Mycobacterium tuberculosis thymidylate kinase. , 2002, Bioorganic & medicinal chemistry letters.

[22]  T. Rungrotmongkol,et al.  Design of oseltamivir analogs inhibiting neuraminidase of avian influenza virus H5N1. , 2009, Antiviral research.

[23]  H. Munier-Lehmann,et al.  Exploring Acyclic Nucleoside Analogues as Inhibitors of Mycobacterium tuberculosis Thymidylate Kinase , 2008, ChemMedChem.

[24]  G. Labesse,et al.  Nucleoside Analogues as Inhibitors of Thymidylate Kinases: Possible Therapeutic Applications , 2002, Chembiochem : a European journal of chemical biology.

[25]  L Wang,et al.  Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[26]  Renxiao Wang,et al.  Comparative evaluation of 11 scoring functions for molecular docking. , 2003, Journal of medicinal chemistry.

[27]  Jan Balzarini,et al.  Rational design of 5'-thiourea-substituted alpha-thymidine analogues as thymidine monophosphate kinase inhibitors capable of inhibiting mycobacterial growth. , 2007, Journal of medicinal chemistry.

[28]  P A Kollman,et al.  Continuum solvent studies of the stability of RNA hairpin loops and helices. , 1998, Journal of biomolecular structure & dynamics.

[29]  Michael K. Gilson,et al.  The inclusion of electrostatic hydration energies in molecular mechanics calculations , 1991, J. Comput. Aided Mol. Des..

[30]  Chandra Verma,et al.  Dissecting the Vibrational Entropy Change on Protein/Ligand Binding: Burial of a Water Molecule in Bovine Pancreatic Trypsin Inhibitor , 2001 .

[31]  Marc Delarue,et al.  Enzymatic and Structural Analysis of Inhibitors Designed against Mycobacterium tuberculosis Thymidylate Kinase , 2003, The Journal of Biological Chemistry.

[32]  Markus H J Seifert,et al.  Targeted scoring functions for virtual screening. , 2009, Drug discovery today.

[33]  Vladimir Frecer,et al.  Design, structure-based focusing and in silico screening of combinatorial library of peptidomimetic inhibitors of Dengue virus NS2B-NS3 protease , 2010, J. Comput. Aided Mol. Des..

[34]  P. Kollman,et al.  Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA. , 2001, Journal of the American Chemical Society.

[35]  C. Frömmel,et al.  The automatic search for ligand binding sites in proteins of known three-dimensional structure using only geometric criteria. , 1996, Journal of molecular biology.

[36]  Hans-Joachim Böhm,et al.  On the use of LUDI to search the Fine Chemicals Directory for ligands of proteins of known three-dimensional structure , 1994, J. Comput. Aided Mol. Des..

[37]  H. Munier-Lehmann,et al.  Substituted benzyl-pyrimidines targeting thymidine monophosphate kinase of Mycobacterium tuberculosis: synthesis and in vitro anti-mycobacterial activity. , 2008, Bioorganic & medicinal chemistry.

[38]  P. Hajduk,et al.  Evaluation of PMF scoring in docking weak ligands to the FK506 binding protein. , 1999, Journal of medicinal chemistry.

[39]  K. Duncan,et al.  Current strategies for identifying and validating targets for new treatment-shortening drugs for TB. , 2007, Current molecular medicine.

[40]  M. Delarue,et al.  X-ray structure of TMP kinase from Mycobacterium tuberculosis complexed with TMP at 1.95 A resolution. , 2001, Journal of molecular biology.

[41]  Gennady M Verkhivker,et al.  Molecular recognition of the inhibitor AG-1343 by HIV-1 protease: conformationally flexible docking by evolutionary programming. , 1995, Chemistry & biology.

[42]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings , 1997 .

[43]  A. Gilles,et al.  Improved spectrophotometric assay of nucleoside monophosphate kinase activity using the pyruvate kinase/lactate dehydrogenase coupling system. , 1994, Analytical biochemistry.

[44]  B. Barrell,et al.  Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence , 1998, Nature.

[45]  G. Labesse,et al.  Comparative Study of Purine and Pyrimidine Nucleoside Analogues Acting on the Thymidylate Kinases of Mycobacterium tuberculosis and of Humans , 2003, Chembiochem : a European journal of chemical biology.

[46]  W. L. Jorgensen,et al.  Prediction of Properties from Simulations: Free Energies of Solvation in Hexadecane, Octanol, and Water , 2000 .

[47]  W L Jorgensen,et al.  Prediction of drug solubility from Monte Carlo simulations. , 2000, Bioorganic & medicinal chemistry letters.

[48]  C. Peloquin,et al.  Tuberculosis pharmacotherapy: strategies to optimize patient care , 2009, Expert opinion on pharmacotherapy.

[49]  Gennady Verkhivker,et al.  Deciphering common failures in molecular docking of ligand-protein complexes , 2000, J. Comput. Aided Mol. Des..

[50]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[51]  W. Parker,et al.  Structure-activity relationship for nucleoside analogs as inhibitors or substrates of adenosine kinase from Mycobacterium tuberculosis. I. Modifications to the adenine moiety. , 2006, Biochemical pharmacology.

[52]  M. Spigelman,et al.  Challenges in tuberculosis drug research and development , 2007, Nature Medicine.

[53]  Hans-Joachim Böhm,et al.  The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure , 1994, J. Comput. Aided Mol. Des..

[54]  Alison R. Gregro,et al.  Evaluating scoring functions for docking and designing β-secretase inhibitors , 2007 .

[55]  H. Munier-Lehmann,et al.  Synthesis and Biological Evaluation of Bicyclic Nucleosides as Inhibitors of M. tuberculosis Thymidylate Kinase , 2006, ChemMedChem.

[56]  M. Delarue,et al.  Crystallization and preliminary X-ray analysis of the thymidylate kinase from Mycobacterium tuberculosis. , 2000, Acta crystallographica. Section D, Biological crystallography.