Tanaji T. Talele

Successful Applications of Computer Aided Drug Discovery: Moving Drugs from Concept to the Clinic

Drug discovery and development is an interdisciplinary, expensive and time-consuming process. Scientific advancements during the past two decades have changed the way pharmaceutical research generate novel bioactive molecules. Advances in computational techniques and in parallel hardware support have enabled in silico methods, and in particular structure-based drug design method, to speed up new target selection through the identification of hits to the optimization of lead compounds in the drug discovery process. This review is focused on the clinical status of experimental drugs that were discovered and/or optimized using computer-aided drug design. We have provided a historical account detailing the development of 12 small molecules (Captopril, Dorzolamide, Saquinavir, Zanamivir, Oseltamivir, Aliskiren, Boceprevir, Nolatrexed, TMI-005, LY-517717, Rupintrivir and NVP-AUY922) that are in clinical trial or have become approved for therapeutic use.

Identification and characterization of coumestans as novel HCV NS5B polymerase inhibitors

The hepatitis C virus (HCV) NS5B is essential for viral RNA replication and is therefore a prime target for development of HCV replication inhibitors. Here, we report the identification of a new class of HCV NS5B inhibitors belonging to the coumestan family of phytoestrogens. Based on the in vitro NS5B RNA-dependent RNA polymerase (RdRp) inhibition in the low micromolar range by wedelolactone, a naturally occurring coumestan, we evaluated the anti-NS5B activity of four synthetic coumestan analogues bearing different patterns of substitutions in their A and D rings, and observed a good structure-activity correlation. Kinetic characterization of coumestans revealed a noncompetitive mode of inhibition with respect to nucleoside triphosphate (rNTP) substrate and a mixed mode of inhibition towards the nucleic acid template, with a major competitive component. The modified order of addition experiments with coumestans and nucleic acid substrates affected the potencies of the coumestan inhibitors. Coumestan interference at the step of NS5B–RNA binary complex formation was confirmed by cross-linking experiments. Molecular docking of coumestans within the allosteric site of NS5B yielded significant correlation between their calculated binding energies and IC50 values. Coumestans thus add to the diversifying pool of anti-NS5B agents and provide a novel scaffold for structural refinement and development of potent NS5B inhibitors.

3D QSAR and Molecular Docking Studies of Benzimidazole Derivatives as Hepatitis C Virus NS5B Polymerase Inhibitors

The urgent need for novel HCV antiviral agents has provided an impetus for understanding the structural requisites of NS5B polymerase inhibitors at the molecular level. Toward this objective, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of 67 HCV NS5B polymerase inhibitors were performed using two methods. First, ligand-based 3D QSAR studies were performed based on the lowest energy conformations employing the atom fit alignment method. Second, receptor-based 3D QSAR models were derived from the predicted binding conformations obtained by docking all NS5B inhibitors at the allosteric binding site of NS5B (PDB ID: 2dxs). Results generated from the ligand-based model were found superior (r2cv values of 0.630 for CoMFA and 0.668 for CoMSIA) to those obtained by the receptor-based model (r2cv values of 0.536 and 0.561 for CoMFA and CoMSIA, respectively). The predictive ability of the models was validated using a structurally diversified test set of 22 compounds that had not been included in a preliminary training set of 45 compounds. The predictive r2 values for the ligand-based CoMFA and CoMSIA models were 0.734 and 0.800, respectively, while the corresponding predictive r2 values for the receptor-based CoMFA and CoMSIA models were 0.538 and 0.639, respectively. The greater potency of the tryptophan derivatives over that of the tyrosine derivatives was interpreted based on CoMFA steric and electrostatic contour maps. The CoMSIA results revealed that for a NS5B inhibitor to have appreciable inhibitory activity it requires hydrogen bond donor and acceptor groups at the 5-position of the indole ring and an R substituent at the chiral carbon, respectively. Interpretation of the CoMFA and CoMSIA contour maps in context of the topology of the allosteric binding site of NS5B provided insight into NS5B-inhibitor interactions. Taken together, the present 3D QSAR models were found to accurately predict the HCV NS5B polymerase inhibitory activity of structurally diverse test set compounds and to yield reliable clues for further optimization of the benzimidazole derivatives in the data set.

Nilotinib potentiates anticancer drug sensitivity in murine ABCB1-, ABCG2-, and ABCC10-multidrug resistance xenograft models.

A panel of clinically used tyrosine kinase inhibitors were compared and nilotinib was found to most potently sensitize specific anticancer agents by blocking the functions of ABCB1/P-glycoprotein, ABCG2/BCRP and ABCC10/MRP7 transporters involved in multi-drug resistance. Nilotinib appreciably enhanced the antitumor response of (1) paclitaxel in the ABCB1- and novel ABCC10-xenograft models, and (2) doxorubicin in a novel ABCG2-xenograft model. With no apparent toxicity observed in the above models, nilotinib attenuated tumor growth synergistically and increased paclitaxel concentrations in ABCB1-overexpressing tumors. The beneficial actions of nilotinib warrant consideration as viable combinations in the clinic with agents that suffer from MDR-mediated insensitivity.

Neratinib reverses ATP-binding cassette B1-mediated chemotherapeutic drug resistance in vitro, in vivo, and ex vivo

Neratinib, an irreversible inhibitor of epidermal growth factor receptor and human epidermal receptor 2, is in phase III clinical trials for patients with human epidermal receptor 2-positive, locally advanced or metastatic breast cancer. The objective of this study was to explore the ability of neratinib to reverse tumor multidrug resistance attributable to overexpression of ATP-binding cassette (ABC) transporters. Our results showed that neratinib remarkably enhanced the sensitivity of ABCB1-overexpressing cells to ABCB1 substrates. It is noteworthy that neratinib augmented the effect of chemotherapeutic agents in inhibiting the growth of ABCB1-overexpressing primary leukemia blasts and KBv200 cell xenografts in nude mice. Furthermore, neratinib increased doxorubicin accumulation in ABCB1-overexpressing cell lines and Rhodamine 123 accumulation in ABCB1-overexpressing cell lines and primary leukemia blasts. Neratinib stimulated the ATPase activity of ABCB1 at low concentrations but inhibited it at high concentrations. Likewise, neratinib inhibited the photolabeling of ABCB1 with [125I]iodoarylazidoprazosin in a concentration-dependent manner (IC50 = 0.24 μM). Neither the expression of ABCB1 at the mRNA and protein levels nor the phosphorylation of Akt was affected by neratinib at reversal concentrations. Docking simulation results were consistent with the binding conformation of neratinib within the large cavity of the transmembrane region of ABCB1, which provides computational support for the cross-reactivity of tyrosine kinase inhibitors with human ABCB1. In conclusion, neratinib can reverse ABCB1-mediated multidrug resistance in vitro, ex vivo, and in vivo by inhibiting its transport function.

The phosphodiesterase-5 inhibitor vardenafil is a potent inhibitor of ABCB1/P-glycoprotein transporter.

One of the major causes of chemotherapy failure in cancer treatment is multidrug resistance (MDR) which is mediated by the ABCB1/P-glycoprotein. Previously, through the use of an extensive screening process, we found that vardenafil, a phosphodiesterase 5 (PDE-5) inhibitor significantly reverses MDR in ABCB1 overexpressing cancer cells, and its efficacy was greater than that of tadalafil, another PDE-5 inhibitor. The present study was designed to determine the reversal mechanisms of vardenafil and tadalafil on ABC transporters-mediated MDR. Vardenafil or tadalafil alone, at concentrations up to 20 µM, had no significant toxic effects on any of the cell lines used in this study, regardless of their membrane transporter status. However, vardenafil when used in combination with anticancer substrates of ABCB1, significantly potentiated their cytotoxicity in ABCB1 overexpressing cells in a concentration-dependent manner, and this effect was greater than that of tadalafil. The sensitivity of the parenteral cell lines to cytotoxic anticancer drugs was not significantly altered by vardenafil. The differential effects of vardenafil and tadalafil appear to be specific for the ABCB1 transporter as both vardenafil and tadalafil had no significant effect on the reversal of drug resistance conferred by ABCC1 (MRP1) and ABCG2 (BCRP) transporters. Vardenafil significantly increased the intracellular accumulation of [3H]-paclitaxel in the ABCB1 overexpressing KB-C2 cells. In addition, vardenafil significantly stimulated the ATPase activity of ABCB1 and inhibited the photolabeling of ABCB1 with [125I]-IAAP. Furthermore, Western blot analysis indicated the incubation of cells with either vardenafil or tadalafil for 72 h did not alter ABCB1 protein expression. Overall, our results suggest that vardenafil reverses ABCB1-mediated MDR by directly blocking the drug efflux function of ABCB1.

Identification of an Allosteric Pocket on Human Hsp70 Reveals a Mode of Inhibition of This Therapeutically Important Protein

Hsp70s are important cancer chaperones that act upstream of Hsp90 and exhibit independent anti-apoptotic activities. To develop chemical tools for the study of human Hsp70, we developed a homology model that unveils a previously unknown allosteric site located in the nucleotide binding domain of Hsp70. Combining structure-based design and phenotypic testing, we discovered a previously unknown inhibitor of this site, YK5. In cancer cells, this compound is a potent and selective binder of the cytosolic but not the organellar human Hsp70s and has biological activity partly by interfering with the formation of active oncogenic Hsp70/Hsp90/client protein complexes. YK5 is a small molecule inhibitor rationally designed to interact with an allosteric pocket of Hsp70 and represents a previously unknown chemical tool to investigate cellular mechanisms associated with Hsp70.

Multiple Transport-Active Binding Sites Are Available for a Single Substrate on Human P-Glycoprotein (ABCB1)

P-glycoprotein (Pgp, ABCB1) is an ATP-Binding Cassette (ABC) transporter that is associated with the development of multidrug resistance in cancer cells. Pgp transports a variety of chemically dissimilar amphipathic compounds using the energy from ATP hydrolysis. In the present study, to elucidate the binding sites on Pgp for substrates and modulators, we employed site-directed mutagenesis, cell- and membrane-based assays, molecular modeling and docking. We generated single, double and triple mutants with substitutions of the Y307, F343, Q725, F728, F978 and V982 residues at the proposed drug-binding site with cys in a cysless Pgp, and expressed them in insect and mammalian cells using a baculovirus expression system. All the mutant proteins were expressed at the cell surface to the same extent as the cysless wild-type Pgp. With substitution of three residues of the pocket (Y307, Q725 and V982) with cysteine in a cysless Pgp, QZ59S-SSS, cyclosporine A, tariquidar, valinomycin and FSBA lose the ability to inhibit the labeling of Pgp with a transport substrate, [125I]-Iodoarylazidoprazosin, indicating these drugs cannot bind at their primary binding sites. However, the drugs can modulate the ATP hydrolysis of the mutant Pgps, demonstrating that they bind at secondary sites. In addition, the transport of six fluorescent substrates in HeLa cells expressing triple mutant (Y307C/Q725C/V982C) Pgp is also not significantly altered, showing that substrates bound at secondary sites are still transported. The homology modeling of human Pgp and substrate and modulator docking studies support the biochemical and transport data. In aggregate, our results demonstrate that a large flexible pocket in the Pgp transmembrane domains is able to bind chemically diverse compounds. When residues of the primary drug-binding site are mutated, substrates and modulators bind to secondary sites on the transporter and more than one transport-active binding site is available for each substrate.

Structure-based virtual screening, synthesis and SAR of novel inhibitors of hepatitis C virus NS5B polymerase.

Hepatitis C virus (HCV) NS5B polymerase is a key target for the development of therapeutic agents aimed at the treatment of HCV infections. Here we report on the identification of novel allosteric inhibitors of HCV NS5B through a combination of structure-based virtual screening, synthesis and structure-activity relationship (SAR) optimization approach. Virtual screening of 260,000 compounds from the ChemBridge database against the tetracyclic indole inhibitor binding pocket of NS5B (allosteric pocket-1, AP-1), sequentially down-sized the library by 4 orders of magnitude to yield 23 candidates. In vitro evaluation of the NS5B inhibitory activity of the in-silico selected compounds resulted in 17% hit rate, identifying two novel chemotypes. Of these, compound 3, bearing the rhodanine scaffold, proved amenable for productive SAR exploration and synthetic modification. As a result, 25 derivatives that exhibited IC₅₀ values ranging from 7.7 to 68.0 μM were developed. Docking analysis of lead compound 28 within the tetracyclic indole- and benzylidene-binding allosteric pockets (AP-1 and AP-3, respectively) of NS5B revealed topological similarities between these two pockets. Compound 28, a novel rhodanine analog with NS5B inhibitory potency in the low micromolar level range may be a promising lead for future development of more potent NS5B inhibitors.

GW583340 and GW2974, human EGFR and HER-2 inhibitors, reverse ABCG2- and ABCB1-mediated drug resistance

The overexpression of ATP binding cassette (ABC) transporters often leads to the development of multidrug resistance (MDR) and results in a suboptimal response to chemotherapy. Previously, we reported that lapatinib (GW572016), a human epidermal growth factor receptor (EGFR) and HER-2 tyrosine kinase inhibitor (TKI), significantly reverses MDR in cancer cells by blocking the efflux function of ABC subfamily B member 1 (ABCB1) and ABC subfamily G member 2 (ABCG2). In the present study, we conducted in vitro experiments to evaluate if GW583340 and GW2974, structural analogues of lapatinib, could reverse ABCB1- and ABCG2-mediated MDR. Our results showed that GW583340 and GW2974 significantly sensitized ABCB1 and ABCG2 overexpressing MDR cells to their anticancer substrates. GW583340 and GW2974 significantly increased the intracellular accumulation of [(3)H]-paclitaxel in ABCB1 overexpressing cells and [(3)H]-mitoxantrone in ABCG2 overexpressing cells respectively. In addition, GW583340 and GW2974 significantly inhibited ABCG2-mediated transport of methotrexate in ABCG2 overexpressing membrane vesicles. There was no significant change in the expression levels of ABCB1 and ABCG2 in the cell lines exposed to 5μM of either GW583340 or GW2974 for 3 days. In addition, a docking model predicted the binding conformation of GW583340 and GW2974 to be within the transmembrane region of homology modeled human ABCB1 and ABCG2. We conclude that GW583340 and GW2974, at clinically achievable plasma concentrations, reverse ABCB1- and ABCG2-mediated MDR by blocking the drug efflux function of these transporters. These findings may be useful in developing combination therapy for cancer treatment with EGFR TKIs.