Design,synthesis and biological evaluation of 4-aminoquinoline-guanylthiourea derivatives as antimalarial agents

Guanylthiourea (GTU) has been identified as an important antifolate antimalarial pharmacophore unit, whereas, 4-amino quinolones are already known for antimalarial activity. In the present work molecules carrying 4-aminoquinoline and GTU moiety have been designed using molecular docking analysis with PfDHFR enzyme and heme unit. The docking results indicated that the necessary interactions (Asp54 and Ile14) and docking score (−9.63 to −7.36 kcal/mmol) were comparable to WR99210 (−9.89 kcal/mol). From these results nine molecules were selected for synthesis. In vitro analysis of these synthesized compounds reveal that out of the nine molecules, eight show antimalarial activity in the range of 0.61–7.55 μM for PfD6 strain and 0.43–8.04 μM for PfW2 strain. Further, molecular dynamics simulations were performed on the most active molecule to establish comparative binding interactions of these compounds and reference ligand with Plasmodium falciparum dihydrofolate reductase (PfDHFR).     

Targeting plant DIHYDROFOLATE REDUCTASE with antifolates and mechanisms for genetic resistance

The folate biosynthetic pathway and its key enzyme dihydrofolate reductase (DHFR) is a popular target for drug development due to its essential role in the synthesis of DNA precursors and some amino acids. Despite its importance, little is known about plant DHFRs, which, like the enzymes from the malarial parasite Plasmodium, are bifunctional, possessing DHFR and thymidylate synthase (TS) domains. Here using genetic knockout lines we confirmed that either DHFR‐TS1 or DHFR‐TS2 (but not DHFR‐TS3) was essential for seed development. Screening mutated Arabidopsis thaliana seeds for resistance to antimalarial DHFR‐inhibitor drugs pyrimethamine and cycloguanil identified causal lesions in DHFR‐TS1 and DHFR‐TS2, respectively, near the predicted substrate‐binding site. The different drug resistance profiles for the plants, enabled by the G137D mutation in DHFR‐TS1 and the A71V mutation in DHFR‐TS2, were consistent with biochemical studies using recombinant proteins and could be explained by structural models. These findings provide a great improvement in our understanding of plant DHFR‐TS and suggest how plant‐specific inhibitors might be developed, as DHFR is not currently targeted by commercial herbicides.     

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     

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.     

Compound prioritization from inverse docking experiment using receptor‐centric and ligand‐centric methods: a case study on Plasmodium falciparum Fab enzymes

Prioritization of compounds using inverse docking approach is limited owing to potential drawbacks in its scoring functions. Classically, molecules ranked by best or lowest binding energies and clustering methods have been considered as probable hits. Mining probable hits from an inverse docking approach is very complicated given the closely related protein targets and the chemically similar ligand data set. To overcome this problem, we present here a computational approach using receptor‐centric and ligand‐centric methods to infer the reliability of the inverse docking approach and to recognize probable hits. This knowledge‐driven approach takes advantage of experimentally identified inhibitors against a particular protein target of interest to delineate shape and molecular field properties and use a multilayer perceptron model to predict the biological activity of the test molecules. The approach was validated using flavone derivatives possessing inhibitory activities against principal antimalarial molecular targets of fatty acid biosynthetic pathway, FabG, FabI and FabZ, respectively. We propose that probable hits can be retrieved by comparing the rank list of docking, quantitative‐structure activity relationship and multilayer perceptron models. Copyright © 2014 John Wiley & Sons, Ltd.     

Ligand‐based pharmacophore modelling and screening of DNA minor groove binders targeting Staphylococcus aureus

The recognition of DNA by small molecules is of special importance in the design of new drugs. Many natural and synthetic compounds have the ability to interact with the minor groove of DNA. In the present study, identification of minor groove binding compounds was attained by the combined approach of pharmacophore modelling, virtual screening and molecular dynamics approach. Experimentally reported 32 minor groove binding compounds were used to develop the pharmacophore model. Based on the fitness score, best three pharmacophore hypotheses were selected and used as template for screening the compounds from drug bank database. This pharmacophore‐based screening provides many compounds with the same pharmacological properties. All these compounds were subjected to four phases of docking protocols with combined Glide‐quantum‐polarized ligand docking approach. Molecular dynamics results indicated that selected compounds are more active and showed good interaction in the binding site of DNA. Based on the scoring parameters and energy values, the best compounds were selected, and antibacterial activity of these compounds was identified using in vitro antimicrobial techniques. Copyright © 2014 John Wiley & Sons, Ltd.     

Exploring dual inhibitory role of febrifugine analogues against Plasmodium utilizing structure-based virtual screening and molecular dynamic simulation

Malaria is an endemic disease caused by the protozoan parasite Plasomodium falciparum. Febrifugine analogues are natural compound obtained from the traditional Chinese herbs have shown significant antimalarial and anticancerous efficacy in experimental model. Development of resistance against the existing antimalarial drug has alarmed the scientific innovators to find a potential antimalarial molecule which can be further used by endemic countries for the elimination of this disease. In this study, structure-based virtual screening and molecular dynamics (MD) base approaches were used to generate potential antimalarial compound against plasmepsin II and prolyl-tRNA synthetase of Plasmodium. Here, we have docked series of febrifugine analogues (n = 11,395) against plasmepsin II in three different docking modes and then it was compared with previously reported target prolyl-tRNA synthetase. Extra precision docking resulted into 235 ligands having better docking score were subject for QikProp analysis. Better ligands (n = 39) obtained from QikProp analysis were subject for ADMET prediction and docking protocol validation through the estimation of receiver operator characteristics. In the later stage, 24 ligands obtained from ADMET study were subject for the estimation of binding energy through MM-GBSA and same were also docked against prolyl-tRNA synthetase to get compounds with dual inhibitor role. Finally, MD simulation and 2D fingerprint MACCS study of two best ligands have shown significant interaction with plasmepsin II and homology against known active ligand with noteworthy MACCS index, respectively. This study concludes that FA12 could be potential drug candidate to fight against Plasmodium falciparum parasites.     

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.     

Novel pyrazole–pyrazoline hybrids endowed with thioamide as antimalarial agents: their synthesis and 3D-QSAR studies

Abstract

One of the most viable options to tackle the growing resistance to the antimalarial drugs is hybrid molecules. It involves combination of different scaffolds in one frame that may lead to compounds with diverse biological profiles. In this context, new hybrids of three different scaffolds viz pyrazole, pyrazoline and thiosemicarbazone moiety were incorporated into one single compound and evaluated for their in vitro schizontocidal activity against the CQ-sensitive 3D7 strain of Plasmodium falciparum. Compounds with significant in vitro antimalarial activity were further evaluated for cytotoxicity against VERO cell lines. The best active compound 48 exhibited an IC₅₀ of 1.13?µM. The in vitro results were further validated by quantitative structure–activity relationship (QSAR).     

Discovery of a series of hydroximic acid derivatives as potent histone deacetylase inhibitors

Abstract

To develop potent histone deacetylase inhibitors as antitumor agents, structural modification was performed. The synthesized molecules were tested by enzymatic inhibition assay and anti-proliferation assay. Several molecules show improved activities in the enzymatic inhibition assay. However, in the MTT assays, all these derived molecules have limited performance compared with SAHA. The IC50 values of molecule ((S)-N-(6-(hydroxyamino)-6-oxohexyl)-4-(3-(2-oxo-1-phenyl-2-((3-(trifluoromethyl)phenyl)amino)ethyl)ureido)benzamide, L8) which has the best enzymatic inhibition activity (with an IC50 value of 11.7?nm and 967?nm against Hela nucleus extract and HDAC8, respectively) were calculated compared with SAHA. Molecular docking was performed to predict the binding mode of molecule L8 in the active site of HDAC2 and HDAC8. Hydrophobic interaction, chelate binding, electrostatic attraction and H-bond interaction in combination make contribution to the ligand–receptor interactions.     

Structure based discovery of small molecule suppressors targeting bacterial lysozyme inhibitors

The production of lysozyme inhibitors, competitively binding to the lysozyme active site, is a bacterial strategy to prevent the lytic activity of host lysozymes. Therefore, suppression of the lysozyme–inhibitor interaction is an interesting new approach for drug development since restoration of the bacterial lysozyme sensitivity will support bacterial clearance from the infected sites. Using molecular modelling techniques the interaction of the Salmonella PliC inhibitor with c-type lysozyme was studied and a protein–protein interaction based pharmacophore model was created. This model was used as a query to identify molecules, with potential affinity for the target, and subsequently, these molecules were filtered using molecular docking. The retained molecules were validated as suppressors of lysozyme inhibitory proteins using in vitro experiments revealing four active molecules.     

New Fluoroquinolone-1,2,4-triazoles as Potent Antibacterial Agents: Design,Synthesis, Docking Studies and in Silico ADME Profiles

Our current work is aimed at synthesizing novel substituted 1,2,4-triazolyl-fluoroquinolone analogs and study of their biological activity to find active promising molecules. The structural elucidation of the products was demonstrated by a variety of spectroscopic methods such as IR, ¹H-NMR, ¹³C-NMR, mass and elemental analysis. The newly synthesized 1,2,4-triazole derivatives were tested in vitro for their ability to inhibit the growth of seven different microbes including S. epidermidis, S. pneumoniae, S. aureus, B. subtilis, K. pneumoniae, E. coli, and P. aeruginosa. Five FQ derivatives 5d , 5e , 5h , 5j , and 5b have demonstrated good antibacterial activity against S. pneumoniae with MICs ranging from 2.5–22.0 μg/mL, while 5c , 5g reported comparable activity against P. aeruginosa with respect to the standard drugs moxifloxacin and ciprofloxacin. The possible mechanism of antibacterial activity of fluoroquinolones was investigated via molecular docking by using DNA gyrase of S. pneumoniae (3RAE). The pefloxacin derivatives also tended a good antibacterial ability based on the results of the molecular docking, ligand 5h with good binding affinity (−9.92 Kcal/mol) and binding site interactions via ValA:86, SerA:79, TyrA:82, MetA:116, AspA:78, AlaA:63, ArgA:117, ProA:112, ProA:113, AlaA:115, AlaA:114. These scaffolds were further evaluated for their ADMET and physicochemical properties by using SwissADME, ADMETlab2.0 web server as a good oral bioavailability.     

New pyrazolylpyrazoline derivatives as dual acting antimalarial-antileishamanial agents: synthesis,biological evaluation and molecular modelling simulations

Promising inhibitory activities of the parasite multiplication were obtained upon evaluation of in vivo antimalarial activities of new pyrazolylpyrazoline derivatives against Plasmodium berghei infected mice. Further evaluation of 5b and 6a against chloroquine-resistant strain (RKL9) of P. falciparum showed higher potency than chloroquine. In vitro antileishmanial activity testing against Leishmania aethiopica promastigote and amastigote forms indicated that 5b, 6a and 7b possessed promising activity compared to miltefosine and amphotericin B deoxycholate. Moreover, antileishmanial activity reversal of the active compounds via folic and folinic acids showed comparable results to the positive control trimethoprim, indicating an antifolate mechanism via targeting leishmanial DHFR and PTR1. The compounds were non-toxic at 125, 250 and 500 mg/kg. In addition, docking of the most active compound against putative malarial target Pf-DHFR-TS and leishmanial PTR1 rationalised the observed activities. Molecular dynamics simulations confirmed a stable and high potential binding of 7a against leishmanial PTR1.     

Synthesis,Antimicrobial Activity and Molecular Docking Study of Novel α‐(Diphenylphosphoryl)‐ and α‐(Diphenylphosphorothioyl)cycloalkanone Oximes

A series of novel α‐(diphenylphosphoryl)‐ and α‐(diphenylphosphorothioyl)cycloalkanone oximes have been synthesized in search for novel bioactive molecules. Their structures were characterized by various spectroscopic methods including IR, NMR (¹H, ³¹P, ¹³C), mass spectrometry and single crystal X‐ray diffraction. The newly synthesized phosphorus‐containing oximes were screened for their in vitro antimicrobial activity against Gram‐positive bacteria (Staphylococcus aureus and Bacillus subtilis), Gram‐negative bacteria (Escherichia coli and Salmonella typhimurium) and fungal strains (Candida albicans and Candida glabrata). The biological assays showed that all the studied compounds exhibited high antibacterial and antifungal activities at only 0.1–2.1 μg/mL. In silico molecular docking studies in FabH enzyme active site were performed in order to predict the possible interaction modes and binding energies of the drug candidates at the molecular level.     

First report on exploring structural requirements of alpha and beta thymidine analogs for PfTMPK inhibitory activity using in silico studies

With the emergence of multi-drug resistance of the currently available antimalarial drugs including the “magic bullet” artemisinin derivatives in the market, there is an urgent need for discovery and development of new potent antimalarial molecules. The present work deals with quantitative structure–activity relationship (QSAR) modeling, pharmacophore mapping and docking studies of a series of 35 thymidine analogs as inhibitors of Plasmodium falciparum thymidylate kinase (PfTMPK), an enzyme that catalyzes phosphorylation of thymidine monophosphate (TMP) to thymidine diphosphate (TDP). The models were validated both internally and externally and significant statistical results were obtained, indicating the robustness and reliability of the developed models. The docking study was performed using the LigandFit option of receptor–ligand interactions protocol section available in Discovery Studio 2.1 where lower RMSD values (0.6931 Å) between the co-crystallized ligand and re-docked ligand assured that the ligand was bound in the same binding pocket. The QSAR, pharmacophore mapping and docking studies provide an understanding of important structural requirements or essential molecular properties, or features of molecules, and important binding interactions, and provide an important guidance for the chemist to synthesis of new molecules with improved PfTMPK inhibitory activity profile. This work revealed the importance of –NH-fragment, electrophilicity of the molecules and the number of oxygen atom towards the PfTMPK inhibitory activity of the molecules. To the best of our knowledge, this work presents the first QSAR and pharmacophore report for thymidine analogs which may serve as an efficient tool for the design and synthesis of potent molecules as PfTMPK inhibitors to address the increasing threat of multi-drug resistance against P. falciparum.     

Role of Molecular Electronic Properties of Some Novel Antimalarial Cyclic Peroxy Ketals in Relation to Potency and Potential Toxicity

Calculated molecular electronicproperties of 20 recently reported cyclic peroxy ketals have been rationalized with their in vitro antimalarial activity with the overall goal to guide the design for safer and effective peroxide-containing antimalarial agents. Stereoelectronicproperties were calculated on the optimized geometry of each compound using the ab initio 3-21G* quantum chemical basis set. Potency appears to be related to electronic properties rather than structural properties such as bond lengths and bond angles though an aliphatic cyclic ring seems to be a structural requirement for potent activity. Electronic properties such as differences in molecular polarity, in electrostatic potential profiles about the peroxide bond and the aromatic ring, in peroxide bond strength, and in the LUMO orbital energy of the molecules are all associated with potency. The three-dimensionalisopotential profile beyondthe van der Waals surface at –10 kcal/mol for the more potent analogs has a distinct large negative potential region by the aromatic ring extending to the methoxy moiety, suggesting a site for initial recognition interaction with the receptor away from the peroxide bond. The HOMO and LUMO isodensity surfaces for all molecules are located on the aromatic ring. The peroxide bond strength of the compounds is around 100 kcal/mol greater than the peroxide-containing clinically used antimalarial artemisinin compounds. In addition, density of the peroxy ketals also appears related to potent activity. The above features of the cyclic peroxy ketals are consistent with these compounds being less potent than the artemisinin compounds, but at the same time are less likely to be as neurotoxic as the artemisinins.     

Tamoxifen,an anticancer drug,is an activator of human aldehyde dehydrogenase 1A1

The modulation of aldehyde dehydrogenase (ALDH) activity has been suggested as a promising option for the prevention or treatment of many diseases. To date, only few activating compounds of ALDHs have been described. In this regard, N‐(1,3‐benzodioxol‐5‐ylmethyl)?2,6‐dichlorobenzamide has been used to protect the heart against ischemia/reperfusion damage. In the search for new modulating ALDH molecules, the binding capability of different compounds to the active site of human aldehyde dehydrogenase class 1A1 (ALDH1A1) was analyzed by molecular docking, and their ability to modulate the activity of the enzyme was tested. Surprisingly, tamoxifen, an estrogen receptor antagonist used for breast cancer treatment, increased the activity and decreased the K_m for NAD⁺ by about twofold in ALDH1A1. No drug effect on human ALDH2 or ALDH3A1 was attained, showing that tamoxifen was specific for ALDH1A1. Protection against thermal denaturation and competition with daidzin suggested that tamoxifen binds to the aldehyde site of ALDH1A1, resembling the interaction of N‐(1,3‐benzodioxol‐5‐ylmethyl)?2,6‐dichlorobenzamide with ALDH2. Further kinetic analysis indicated that tamoxifen activation may be related to an increase in the K_d for NADH, favoring a more rapid release of the coenzyme, which is the rate‐limiting step of the reaction for this isozyme. Therefore, tamoxifen might improve the antioxidant response, which is compromised in some diseases. Proteins 2015; 83:105–116. © 2014 Wiley Periodicals, Inc.     

Towards in silico lead optimization: scores from ensembles of protein/ligand conformations reliably correlate with biological activity

Accurately ranking protein/ligand interactions and distinguishing subtle differences between homologous compounds in a virtual focused library in silico is essential in a structure-based drug discovery program. In order to establish a predictive model to design novel inhibitors of dihydrofolate reductase (DHFR) from the parasitic protozoa, Cryptosporidium hominis, we docked a series of 30 DHFR inhibitors with measured inhibition constants against the crystal structure of the protein. By including protein flexibility and averaging the energies of the 25 lowest protein/ligand conformers we obtained more accurate total nonbonded energies from which we calculated a predicted biological activity. The calculated and measured biological activities showed reliable correlations of 72.9%. Additionally, visual analysis of the ensemble of protein/ligand conformations revealed alternative ligand binding pockets in the active site. Using the same principles we then created a homology model of DHFR from Toxoplasma gondii and docked 11 inhibitors. A correlation of 50.2% between docking score and activity validates both the method and the model. The correlations presented here are particularly compelling considering the high structural similarity of the ligands and the fact that we have used structures derived from crystallographic data and homology modeling. These docking principles may be useful in any lead optimization study where accurate ranking of similar compounds is desired.     

2-Alkyl(aryl)-quinazolin-4(3H)-thiones, 2-R-(quinazolin-4(3H)-ylthio)carboxylic acids and amides: synthesis,molecular docking,antimicrobial and anticancer properties

In this study, a series of novel 2-alkyl(aryl)-quinazolin-4(3H)-thiones, 2-R-(quinazolin-4(3H)-ylthio)carboxylic acids and amides were synthesized and evaluated for antimicrobial and anticancer activities. Their structure was confirmed by elemental analysis and spectral data (FT-IR, LC-MS, ¹H-NMR). Antimicrobial activity was tested in vitro against Staphylococcus aureus, Enterococcus faecalis, Enterobacter aerogenes, Pseudomonas aeruginosa, Escherichia coli, Klebsiella pneumonia, Candida albicans and NCI in vitro preliminary anticancer activity against nine different cancer types. The most active antibacterial and antifungal compounds were: 2.1, 2.2 and 2.4. The introduction of the carboxylic acid or amide residue into the fourth position of quinazolin-4(3H)-thione resulted in the absence of antimicrobial activity. Substance 3.8 inhibited renal cancer UO-31 line and 2.18 – leukemia CCRF-CEM. The results of in silico molecular docking for DHFR and CK2 kinase had no correlation with in vitro properties, proposing the presence of other biological activity pathways.     

Selective peptide inhibitors of bifunctional thymidylate synthase‐dihydrofolate reductase from Toxoplasma gondii provide insights into domain–domain communication and allosteric regulation

The bifunctional enzyme thymidylate synthase–dihydrofolate reductase (TS–DHFR) plays an essential role in DNA synthesis and is unique to several species of pathogenic protozoans, including the parasite Toxoplasma gondii. Infection by T. gondii causes the prevalent disease toxoplasmosis, for which TS–DHFR is a major therapeutic target. Here, we design peptides that target the dimer interface between the TS domains of bifunctional T. gondii TS–DHFR by mimicking β‐strands at the interface, revealing a previously unknown allosteric target. The current study shows that these β‐strand mimetic peptides bind to the apo‐enzyme in a species‐selective manner to inhibit both the TS and distal DHFR. Fluorescence spectroscopy was used to monitor conformational switching of the TS domain and demonstrate that these peptides induce a conformational change in the enzyme. Using structure‐guided mutagenesis, nonconserved residues in the linker between TS and DHFR were identified that play a key role in domain–domain communication and in peptide inhibition of the DHFR domain. These studies validate allosteric inhibition of apo‐TS, specifically at the TS–TS interface, as a potential target for novel, species‐specific therapeutics for treating T. gondii parasitic infections and overcoming drug resistance.