Specificity in structure-based drug design: Identification of a novel,selective inhibitor of Pneumocystis carinii dihydrofolate reductase

Specificity is an important aspect of structure-based drug design. Distinguishing between related targets in different organisms is often the key to therapeutic success. Pneumocystis carinii is a fungal opportunist which causes a crippling pneumonia in immunocompromised individuals. We report the identification of novel inhibitors of P. cariniidihydrofolate reductase (DHFR) that are selective versus inhibition of human DHFR using computational molecular docking techniques. The Fine Chemicals Directory, a database of commercially available compounds, was screened with the DOCK program suite to produce a list of potential P. carinii DHFR inhibitors. We then used a postdocking refinement directed at discerning subtle structural and chemical features that might reflect species specificity. Of 40 compounds predicted to exhibit anti-PneumocystisDHFR activity, each of novel chemical framework, 13 (33%) show IC₅₀ values better than 150 μM in an enzyme assay. These inhibitors were further assayed against human DHFR: 10 of the 13 (77%) bind preferentially to the fungal enzyme. The most potent compound identified is a 7 μM inhibitor of P. carinii DHFR with 25-fold selectivity. The ability of molecular docking methods to locate selective inhibitors reinforces our view of structure-based drug discovery as a valuable strategy, not only for identifying lead compounds, but also for addressing receptor specificity. Proteins 29:59–67, 1997. © 1997 Wiley-Liss, Inc.     

Virtual screening,docking, and dynamics of potential new inhibitors of dihydrofolate reductase from Yersinia pestis

In the present work, we propose to design drugs that target the enzyme dihydrofolate redutase (DHFR) as a means of a novel drug therapy against plague. Potential inhibitors of DHFR from Yersinia pestis (YpDHFR) were selected by virtual screening and subjected to docking, molecular dynamics (MD) simulations, and Poisson–Boltzmann surface area method, in order to evaluate their interactions in the active sites of YpDHFR and human DHFR (HssDHFR). The results suggested selectivity for three compounds that were further used to propose the structures of six new potential selective inhibitors for YpDHFR.     

The influence of glycerol on the binding of methotrexate to Mycobacterium tuberculosis dihydrofolate reductase: a molecular modelling study

The enzyme, dihydrofolate reductase (DHFR), from Mycobacterium tuberculosis (mt-DHFR) is believed to be a potential drug target for the treatment of tuberculosis. The co-crystal structure of mt-DHFR bound with glycerol (GOL), NAPDH and methotrexate (MTX) reveals a GOL binding site on the enzyme. This GOL binding site could be very important for the design of novel, selective mt-DHFR inhibitors, because this binding site is absent on human DHFR (h-DHFR). We have performed molecular dynamic simulations and free energy calculations to evaluate the binding affinity of GOL and its free energy contribution to the binding of MTX to mt-DHFR. The results showed that GOL does not bind tightly to mt-DHFR. Although GOL itself contributed free energy on MTX binding to mt-DHFR, GOL also increased the flexibilities of MTX, so that MTX cannot maintain strong electronic interactions with ARG32 and ARG60, which caused the total binding free energy to decrease. These data suggest that GOL binding is weak and it could be expelled from the binding site, to allow inhibitors containing appropriate side chains to bind. This observation can be used to inform future drug design studies, especially those aimed at improving drug selectivity against h-DHFR.     

Multicomplex-based pharmacophore modeling in conjunction with multi-target docking and molecular dynamics simulations for the identification of PfDHFR inhibitors

Abstract

Plasmodium falciparum dihydrofolate reductase enzyme (PfDHFR) is counted as one of the attractive and validated antimalarial drug targets. However, the point mutations in the active site of wild-type PfDHFR have developed resistance against the well-known antifolates. Therefore, there is a dire need for the development of inhibitors that can inhibit both wild-type and mutant-type DHFR enzyme. In the present contribution, we have constructed the common feature pharmacophore models from the available PfDHFR. A representative hypothesis was prioritized and then employed for the screening of natural product library to search for the molecules with complementary features responsible for the inhibition. The screened candidates were processed via drug-likeness filters and molecular docking studies. The docking was carried out on the wild-type PfDHFR (3QGT); double-mutant PfDHFR (3UM5 and 1J3J) and quadruple-mutant PfDHFR (1J3K) enzymes. A total of eight common hits were obtained from the docking calculations that could be the potential inhibitors for both wild and mutant type DHFR enzymes. Eventually, the stability of these candidates with the selected proteins was evaluated via molecular dynamics simulations. Except for SPECS14, all the prioritized candidates were found to be stable throughout the simulation run. Overall, the strategy employed in the present work resulted in the retrieval of seven candidates that may show inhibitory activity against PfDHFR and could be further exploited as a scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Structural genomics for drug design against the pathogen Coxiella burnetii

Coxiella burnetii is a highly infectious bacterium and potential agent of bioterrorism. However, it has not been studied as extensively as other biological agents, and very few of its proteins have been structurally characterized. To address this situation, we undertook a study of critical metabolic enzymes in C. burnetii that have great potential as drug targets. We used high‐throughput techniques to produce novel crystal structures of 48 of these proteins. We selected one protein, C. burnetii dihydrofolate reductase (CbDHFR), for additional work to demonstrate the value of these structures for structure‐based drug design. This enzyme's structure reveals a feature in the substrate binding groove that is different between CbDHFR and human dihydrofolate reductase (hDHFR). We then identified a compound by in silico screening that exploits this binding groove difference, and demonstrated that this compound inhibits CbDHFR with at least 25‐fold greater potency than hDHFR. Since this binding groove feature is shared by many other prokaryotes, the compound identified could form the basis of a novel antibacterial agent effective against a broad spectrum of pathogenic bacteria. Proteins 2015; 83:2124–2136. © 2015 Wiley Periodicals, Inc.     

Preventing the return of smallpox: molecular modeling studies on thymidylate kinase from Variola virus

Smallpox was one of the most devastating diseases in the human history and still represents a serious menace today due to its potential use by bioterrorists. Considering this threat and the non-existence of effective chemotherapy, we propose the enzyme thymidylate kinase from Variola virus (VarTMPK) as a potential target to the drug design against smallpox. We first built a homology model for VarTMPK and performed molecular docking studies on it in order to investigate the interactions with inhibitors of Vaccinia virus TMPK (VacTMPK). Subsequently, molecular dynamics (MD) simulations of these compounds inside VarTMPK and human TMPK (HssTMPK) were carried out in order to select the most promising and selective compounds as leads for the design of potential VarTMPK inhibitors. Results of the docking and MD simulations corroborated to each other, suggesting selectivity towards VarTMPK and, also, a good correlation with the experimental data.     

Cloning,recombinant expression and inhibitor profiles of dihydrofolate reductase from the Australian sheep blow fly,Lucilia cuprina

While dihydrofolate reductase (DHFR) is an important drug target in mammals, bacteria and protozoa, no inhibitors of this enzyme have been developed as commercial insecticides. We therefore examined the potential of this enzyme as a drug target in an important ectoparasite of livestock, the Australian sheep blow fly, Lucilia cuprina (Diptera: Calliphoridae) (Wiedemann). The non‐specific DHFR inhibitors aminopterin and methotrexate significantly inhibited the growth of L. cuprina larvae, with IC₅₀ values at µg levels. Trimethoprim and pyrimethamine were 5–30‐fold less active. Relative IC₅₀ values for the inhibition of recombinant L. cuprina DHFR by various inhibitors were in accordance with their relative effects on larval growth. The active‐site amino acid residues of L. cuprina DHFR differed by between 34% and 50% when compared with two mammalian species, as well as two bacteria and two protozoa. There were significant charge and size differences in specific residues between the blow fly and human DHFR enzymes, notably the L. cuprina Asn21, Lys31 and Lys63 residues. This study provides bioassay evidence to highlight the potential of blow fly DHFR as an insecticide target, and describes differences in active site residues between blow flies and other organisms which could be exploited in the design of blow fly control chemicals.     

A molecular dynamics study of components of the ginger (Zingiber officinale) extract inside human acetylcholinesterase: implications for Alzheimer disease

Components of ginger (Zingiber officinale) extracts have been described as potential new drug candidates against Alzheimer disease (AD), able to interact with several molecular targets related to the AD treatment. However, there are very few theoretical studies in the literature on the possible mechanisms of action by which these compounds can work as potential anti-AD drugs. For this reason, we performed here docking, molecular dynamic simulations and mmpbsa calculations on four components of ginger extracts former reported as active inhibitors of human acetylcholinesterase (HssAChE), and compared our results to the known HssAChE inhibitor and commercial drug in use against AD, donepezil (DNP). Our findings points to two among the compounds studied: (E)-1,7-bis(4-hydroxy-3-methoxyphenyl)hept-4-en-3-on and 1-(3,4-dihydroxy-5-methoxyphenyl)-7-(4-hydroxy-3- ethoxyphenyl) heptane-3,5-diyl diacetate, as promising new HssAChE inhibitors that could be as effective as DNP. We also mapped the binding of the studied compounds in the different binding pockets inside HssAChE and established the preferred interactions to be favored in the design of new and more efficient inhibitors.     

Analysis of Bacillus anthracis nucleoside hydrolase via in silico docking with inhibitors and molecular dynamics simulation

As the enzyme nucleoside hydrolase (NH) is widely found in nature but has not yet been detected in mammals, it is considered an ideal target in the development of chemotherapy against parasitic diseases and bacterial infections like anthrax. Considering the risk that this biological warfare agent represents nowadays, the search for new drugs and new molecular targets in the development of chemotherapy against anthrax is imperative. On this basis, we performed docking studies of six known NH inhibitors at the active site of NH from Bacillus anthracis (BaNH). Subsequently, molecular dynamics (MD) simulations of these compounds inside BaNH were carried out in order to complement the docking studies and select the most promising compounds as leads for the design of potential BaNH inhibitors. Most of the docking and MD results obtained agreed well with each other and showed good correlation with experimental data.     

Mitogen activated protein kinase-1 and cell division control protein-42 are putative targets for the binding of novel natural lead molecules: a therapeutic intervention against Candida albicans

Abstract

Candida albicans, fungal yeast causes several lethal infections in immune-suppressed patients and recently emerged as drug-resistant pathogens worldwide. The present study aimed to screen putative drug targets of Candia albicans and to study the binding potential of novel natural lead compounds towards these targets by computational virtual screening and molecular dynamic (MD) simulation. Through extensive analysis of mitogen-activated protein kinase (MAPK) signalling pathways, mitogen-activated protein kinase-1 (HOG1) and cell division control protein-42 (CDC42) genes were prioritized as putative targets based on their virulent functions. The three-dimensional structures of these genes, not available in their native forms, were computationally modeled and validated. 76 lead molecules from various natural sources were screened and their drug likeliness and pharmacokinetic features were predicted. Among these ligands, two lead molecules that demonstrated ideal drug-likeliness and pharmacokinetic features were docked against HOG1 and CDC42 and their binding potential was compared with the binding of conventional drug Fluconazole with their usual target. The prediction was computationally validated by MD simulation. The current study revealed that Cudraxanthone-S present in Cudrania cochinchinensis and Scutifoliamide-B present in Piper scutifolium exhibited ideal drug likeliness, pharmacokinetics and binding potential to the prioritized targets in comparison with the binding of Fluconazole and their usual target. MD simulation showed that CDC42-Cudraxanthone-S and HOG1-Scutifoliamide-B complexes were exhibited stability throughout MD simulation. Thus, the study provides significant insight into employing HOG1 and CDC42 of MAPK as putative drug targets of C. albicans and Cudraxanthone-S and Scutifoliamide-B as potential inhibitors for drug discovery.

Communicated by Ramaswamy H. Sarma     

Comparative assessment of the therapeutic drug targets of C. botulinum ATCC 3502 and C. difficile str. 630 using in silico subtractive proteomics approach

Growing antimicrobial resistance of the pathogens against multiple drugs posed a serious threat to the human health worldwide. This fueled the need of identifying the novel therapeutic targets that can be used for developing new class of the drugs. Recently, there is a substantial rise in the rate of Clostridium infections as well as in the emergence of virulent and antibiotic resistant strains. Hence, there is an urgent need for the identification of potential therapeutic targets and the development of new drugs for the treatment and prevention of Clostridium infections. In the present study, a combinatorial approach involving systems biology and comparative genomics strategy was tested against Clostridium botulinum ATCC 3502 and Clostridium difficile str. 630 pathogens, to render potential therapeutic target at qualitative and quantitative level. This resulted in the identification of five common (present in both the pathogens, 34 in C. botulinum ATCC 3502 and 42 in C. difficile str. 630) drug targets followed by virtual screening–based identification of potential inhibitors employing molecular docking and molecular dynamics simulations. The identified targets will provide a solid platform for the designing of novel wide-spectrum lead compounds capable of inhibiting their catalytic activities against multidrug-resistant Clostridium in the near future.     

Lipases secreted by a gut bacterium inhibit arbovirus transmission in mosquitoes

Arboviruses are etiological agents of various severe human diseases that place a tremendous burden on global public health and the economy; compounding this issue is the fact that effective prophylactics and therapeutics are lacking for most arboviruses. Herein, we identified 2 bacterial lipases secreted by a Chromobacterium bacterium isolated from Aedes aegypti midgut, Chromobacterium antiviral effector-1 (CbAE-1) and CbAE-2, with broad-spectrum virucidal activity against mosquito-borne viruses, such as dengue virus (DENV), Zika virus (ZIKV), Japanese encephalitis virus (JEV), yellow fever virus (YFV) and Sindbis virus (SINV). The CbAEs potently blocked viral infection in the extracellular milieu through their lipase activity. Mechanistic studies showed that this lipase activity directly disrupted the viral envelope structure, thus inactivating infectivity. A mutation in the lipase motif of CbAE-1 fully abrogated the virucidal ability. Furthermore, CbAEs also exert lipase-dependent entomopathogenic activity in mosquitoes. The anti-arboviral and entomopathogenic properties of CbAEs render them potential candidates for the development of novel transmission control strategies against vector-borne diseases.     

Human dihydrofolate reductase and thymidylate synthase form a complex in vitro and co-localize in normal and cancer cells

Enzymes involved in thymidylate biosynthesis, thymidylate synthase (TS), and dihydrofolate reductase (DHFR) are well-known targets in cancer chemotherapy. In this study, we demonstrated for the first time, that human TS and DHFR form a strong complex in vitro and co-localize in human normal and colon cancer cell cytoplasm and nucleus. Treatment of cancer cells with methotrexate or 5-fluorouracil did not affect the distribution of either enzyme within the cells. However, 5-FU, but not MTX, lowered the presence of DHFR-TS complex in the nucleus by 2.5-fold. The results may suggest the sequestering of TS by FdUMP in the cytoplasm and thereby affecting the translocation of DHFR-TS complex to the nucleus. Providing a strong likelihood of DHFR-TS complex formation in vivo, the latter complex is a potential new drug target in cancer therapy. In this paper, known 3D structures of human TS and human DHFR, and some protozoan bifunctional DHFR-TS structures as templates, are used to build an in silico model of human DHFR–TS complex structure, consisting of one TS dimer and two DHFR monomers. This complex structure may serve as an initial 3D drug target model for prospective inhibitors targeting interfaces between the DHFR and TS enzymes.     

The systematic modeling studies and free energy calculations of the phenazine compounds as anti-tuberculosis agents

Abstract

Phenazine compounds have good activity against Mycobacterium tuberculosis (MTB). Based on the reported activities that were obtained in MTB H37Rv, a three-dimensional quantitative structure–activity relationship (3D-QSAR) model was built to design novel compounds against MTB. A fivefold cross-validation method and external validation were used to analyze the accuracy of forecasting. The model has a cross-validation coefficient q²=0.7 and a non-cross-validation coefficient r² = 0.903, indicating that the model has good predictive possibility. The design of anti-pneumococcus MTB compounds was guided by the obtained 3D-QSAR model, and several compounds with better activity were obtained. To test the activity of these compounds, molecular docking, molecular dynamics simulation, and post-simulation analysis of the already reported drug targets in MTB were carried out. Among the total 15 drug targets, only three targets (Rv2361c, Rv2965c, and Rv3048c) were selected based on the docking results. Initial results reported that these compounds possessed good inhibition activity for Rv2361c. The top nine complexes of Rv2361 ligands were only subjected to MD simulation which resulted in a stable dynamics of the structures and showed a residual fluctuation in inhibitors binding pocket. Free energy reported that overall, the derivatives hold strong energy against the protein target. Energetic contribution results showed that residues, Asp76, Arg80, Asn124, Arg127, Arg244, and Arg250, play a major role in total energy. Systems biology approach validates shortlisted drug effect on the entire system which might be useful to predict potential drug in wet lab as well.

Communicated by Ramaswamy H. Sarma     

Design,synthesis, biological evaluation and computational investigation of novel inhibitors of dihydrofolate reductase of opportunistic pathogens

The present work deals with design, synthesis and biological evaluation of novel, diverse compounds as potential inhibitors of dihydrofolate reductase (DHFR) from opportunistic microorganisms; Pneumocystis carinii (pc), Toxoplasma gondii (tg) and Mycobacterium avium (ma). A set of 14 structurally diverse compounds were designed with varying key pharmacophoric features of DHFR inhibitors, bulky distal substitutions and different bridges joining the distal part and 2,4-diaminopyrimidine nucleus. The designed compounds were synthesized and evaluated in enzyme assay against pc, tg and ma DHFR. The rat liver (rl) DHFR was used as mammalian standard. As the next logical step of the project, flexible molecular docking studies were carried out to predict the binding modes of these compounds in pcDHFR active site and the obtained docked poses were post processed using MM-GBSA protocol for prediction of relative binding affinity. The predicted binding modes were able to rationalize the experimental results in most cases. Of particular interest, both the docking scores and MM-GBSA predicted ΔG_bind were able to distinguish between the active and low active compounds. Furthermore, good correlation coefficient of 0.797 was obtained between the IC₅₀ values and MM-GBSA predicted ΔG_bind. Taken together, the current work provides not only a novel scaffold for further optimization of DHFR inhibitors but also an understanding of the specific interactions of inhibitors with DHFR and structural modifications that improve selectivity.     

Nanoparticle Fullerene (C60) demonstrated stable binding with antibacterial potential towards probable targets of drug resistant Salmonella typhi – a computational perspective and in vitro investigation

Salmonella typhi, a Gram negative bacterium, has become multidrug resistant (MDR) to wide classes of antibacterials which necessitate an alarming precaution. This study focuses on the binding potential and therapeutic insight of Nano-Fullerene C60 towards virulent targets of Salmonella typhi by computational prediction and preliminary in vitro assays. The clinical isolates of Salmonella typhi were collected and antibiotic susceptibility profiles were assessed. The drug targets of pathogen were selected by rigorous literature survey and gene network analysis by various metabolic network resources. Based on this study, 20 targets were screened and the 3D structures of few drug targets were retrieved from PDB and others were computationally predicted. The structures of nanoleads such as Fullerene C60, ZnO and CuO were retrieved from drug databases. The binding potential of these nanoleads towards all selected targets were predicted by molecular docking. The best docked conformations were screened and concept was investigated by preliminary bioassays. This study revealed that most of the isolates of Salmonella typhi were found to be MDR (p < .05). The theoretical models of selected drug targets showed high stereochemical validity. The molecular docking studies suggested that Fullerene C60 showed better binding affinity towards the drug targets when compared to ZnO and CuO. The preliminary in vitro assays suggested that 100 μg/L Fullerene C60 posses significant inhibitory activities and absence of drug resistance to this nanoparticle. This study suggests that Fullerene C60 can be scaled up as probable lead molecules against the major drug targets of MDR Salmonella typhi.     

Exploring insights for virulent gene inhibition of multidrug resistant Salmonella typhi,Vibrio cholerae,and Staphylococcus areus by potential phytoligands via in silico screening

In our recent studies on prevalence of multidrug resistant pathogens in Byramangala reservoir, Karnataka, India, we identified Salmonella typhi, Staphylococcus aureus, and Vibrio cholerae which had acquired multiple drug resistance (MDR) and emerged as superbugs. Hence, there is a pressing demand to identify alternative therapeutic remedies. Our study focused on the screening of herbal leads by structure-based virtual screening. The virulent gene products of these pathogens towards Kanamycin(aph), Trimethoprim(dfrA1), Methicillin (mecI), and Vancomycin (vanH) were identified as the probable drug targets and their 3D structures were predicted by homology modeling. The predicted models showed good stereochemical validity. By extensive literature survey, we selected 58 phytoligands and their drug likeliness and pharmacokinetic properties were computationally predicted. The inhibitory properties of these ligands against drug targets were studied by molecular docking. Our studies revealed that Baicalein from S. baicalensis (baikal skullcap) and Luteolin from Taraxacum officinale (dandelion) were identified as potential inhibitors against aph of S. typhi. Resveratrol from Vitis vinifera (grape vine) and Wogonin from S. baicalensis were identified as potential inhibitors against dfrA1 of S. typhi. Herniarin from Herniaria glabra (rupture worts) and Pyrocide from Daucus carota (Carrot) were identified as the best leads against dfrA1 of V. cholerae. Taraxacin of T. officinale (weber) and Luteolin were identified as potential inhibitors against Mec1. Apigenin from Coffee arabica (coffee) and Luteolin were identified as the best leads against vanH of S. aureus. Our findings pave crucial insights for exploring alternative therapeutics against MDR pathogens.     

Molecular docking studies on DMDP derivatives as human DHFR inhibitors

Molecular docking is routinely used for understanding drug‐receptor interaction in modern drug design. Here, we describe the docking of 2, 4-diamino-5-methyl-5-deazapteridine (DMDP) derivatives as inhibitors to human dihydrofolate reductase (DHFR). We docked 78 DMDP derivates collected from literature to DHFR and studied their specific interactions with DHFR. A new shape-based method, LigandFit, was used for docking DMDP derivatives into DHFR active sites. The result indicates that the molecular docking approach is reliable and produces a good correlation coefficient (r² = 0.499) for the 73 compounds between docking score and IC₅₀ values (Inhibitory Activity). The chloro substituted naphthyl ring of compound 63 makes significant hydrophobic contact with Leu 22, Phe 31 and Pro 61 of the DHFR active site leading to enhanced inhibition of the enzyme. The docked complexes provide better insights to design more potent DHFR inhibitors prior to their synthesis.     

Genome‐wide profiling of humoral immune response to Coxiella burnetii infection by protein microarray

Comprehensive evaluation of the humoral immune response to Coxiella burnetii may identify highly needed diagnostic antigens and potential subunit vaccine candidates. Here we report the construction of a protein microarray containing 1901 C. burnetii ORFs (84% of the entire proteome). This array was probed with Q‐fever patient sera and naïve controls in order to discover C. burnetii‐specific seroreactive antigens. Among the 21 seroreactive antigens identified, 13 were significantly more reactive in Q‐fever cases than naïve controls. The remaining eight antigens were cross‐reactive in both C. burnetii infected and naïve patient sera. An additional 64 antigens displayed variable seroreactivity in Q‐fever patients, and underscore the diversity of the humoral immune response to C. burnetii. Nine of the differentially reactive antigens were validated on an alternative immunostrip platform, demonstrating proof‐of‐concept development of a consistent, safe, and inexpensive diagnostic assay alternative. Furthermore, we report here the identification of several new diagnostic antigens and potential subunit vaccine candidates for the highly infectious category B alphaproteobacteria, C. burnetii.     

Designing of inhibitors against CTX-M-15 type β-lactamase: potential drug candidate against β-lactamases-producing multi-drug-resistant bacteria

CTX-M-15 are the most prevalent types of β-lactamases that hydrolyze almost all antibiotics of β-lactam group lead to multiple-antibiotic resistance in bacteria. Three β-lactam inhibitors are available for use in combination with different antibiotics of cephalosporine group against the CTX-M-15-producing strains. Therefore, strategies to identify novel anti β-lactamase agents with specific mechanisms of action are the need of an hour. In this study, we screened three novel non-β-lactam inhibitors against CTX-M-15 by multi-step virtual screening approach. The potential for virtually screened drugs was estimated through in vitro cell assays. Hence, we proposed a study to understand the binding mode of CTX-M-15 with inhibitors using bioinformatics and experimental approach. We calculated the dissociation constants (K_d), association constant (K_a), stoichiometry (n) and binding energies (ΔG) of compounds with the respective targets. Molecular dynamic simulation carried out for 25 ns, revealed that these complexes were found stable throughout the simulation with relative RMSD in acceptable range. Moreover, microbiological and kinetic studies further confirmed high efficacies of these inhibitors by reducing the minimum inhibitory concentration (MIC) and catalysis of antibiotics by β-lactamases in the presence of inhibitors. Therefore, we conclude that these potential inhibitors may be used as a lead molecule for future drug candidates against β-lactamases-producing bacteria.