In vitro efficacy of 2,N-bisarylated 2-ethoxyacetamides against Plasmodium falciparum.

Investigation of a series of 2,N-bisarylated 2-ethoxyacetamides resulted in the identification of four inhibitors 5, 20, 24, 29 with single-digit micromolar in vitro efficacy against two drug-resistant Plasmodium falciparum strains. These compounds are analogs of structurally-related 1,3-bisaryl-2-propen-1-ones (chalcones), the latter showing efficacy in vitro but not in a malaria-infected mouse. The 2,N-bisarylated 2-ethoxyacetamides (e.g., 2, 5, 20) were shown to possess significantly greater stability in the presence of metabolizing enzymes than the corresponding 1,3-bisaryl-2-propen-1-ones (e.g., 1, 3, 18).

Discovery of subtype selective muscarinic receptor antagonists as alternatives to atropine using in silico pharmacophore modeling and virtual screening methods.

Muscarinic acetylcholine receptors (mAChRs) have five known subtypes which are widely distributed in both the peripheral and central nervous system for regulation of a variety of cholinergic functions. Atropine is a well known muscarinic subtype non-specific antagonist that competitively inhibits acetylcholine (ACh) at postganglionic muscarinic sites. Atropine is used to treat organophosphate (OP) poisoning and resulting seizures in the warfighter because it competitively inhibits acetylcholine (ACh) at the muscarinic cholinergic receptors. ACh accumulates due to OP inhibition of acetylcholinesterase (AChE), the enzyme that hydrolyzes ACh. However, atropine produces several unwanted side-effects including dilated pupils, blurred vision, light sensitivity, and dry mouth. To overcome these side-effects, our goal was to find an alternative to atropine that emphasizes M1 (seizure prevention) antagonism but has minimum M2 (cardiac) and M3 (e.g., eye) antagonism so that an effective less toxic medical countermeasure may be developed to protect the warfighter against OP and other chemical warfare agents (CWAs). We adopted an in silico pharmacophore modeling strategy to develop features that are characteristics of known M1 subtype-selective compounds and used the model to identify several antagonists by screening an in-house (WRAIR-CIS) compound database. The generated model for the M1 selectivity was found to contain two hydrogen bond acceptors, one aliphatic hydrophobic, and one ring aromatic feature distributed in a 3D space. From an initial identification of about five hundred compounds, 173 compounds were selected through principal component and cluster analyses and in silico ADME/Toxicity evaluations. Next, these selected compounds were evaluated in a subtype-selective in vitro radioligand binding assay. Twenty eight of the compounds showed antimuscarinic activity. Nine compounds showed specificity for M1 receptors and low specificity for M3 receptors. The pK(i) values of the compounds range from 4.5 to 8.5 nM in comparison to a value of 8.7 nM for atropine. 2-(diethylamino)ethyl 2,2-diphenylpropanoate (ZW62841) was found have the best desired selectivity. None of the newly found compounds were previously reported to exhibit antimuscarinic specificity. Both theoretical and experimental results are presented.

In vitro efficacy of 7-benzylamino-1-isoquinolinamines against Plasmodium falciparum related to the efficacy of chalcones.

A series of 1,7-diaminoisoquinolinamines, that are expected to mediate antimalarial activity by the same mechanism employed by the chalcones, were produced. Six 7-benzylamino-1-isoquinolinamines were found to be submicromolar inhibitors in vitro of drug-resistant Plasmodium falciparum, with the best possessing activity comparable to chloroquine. Despite being developed from a lead that is a DHFR inhibitor, these compounds do not mediate their antimalarial effects by inhibition of DHFR.

In silico pharmacophore model for tabun-inhibited acetylcholinesterase reactivators: a study of their stereoelectronic properties.

Organophosphorus (OP) nerve agents that inhibit acetylcholinesterase (AChE; EC 3.1.1.7) function in the nervous system, causing acute intoxication. If untreated, death can result. Inhibited AChE can be reactivated by oximes, antidotes for OP exposure. However, OP intoxication caused by the nerve agent tabun (GA) is particularly resistant to oximes, which poorly reactivate GA-inhibited AChE. In an attempt to develop a rational strategy for the discovery and design of novel reactivators with lower toxicity and increased efficacy in reactivating GA-inhibited AChE, we developed the first in silico pharmacophore model for binding affinity of GA-inhibited AChE from a set of 11 oximes. Oximes were analyzed for stereoelectronic profiles and three-dimensional quantitative structure-activity relationship pharmacophores using ab initio quantum chemical and pharmacophore generation methods. Quantum chemical methods were sequentially used from semiempirical AM1 to hierarchical ab initio calculations to determine the stereoelectronic properties of nine oximes exhibiting affinity for binding to GA-inhibited AChE in vivo. The calculated stereoelectronic properties led us to develop the in silico pharmacophore model using CATALYST methodology. Specific stereoelectronic profiles including the distance between bisquarternary nitrogen atoms of the pyridinium ring in the oximes, hydrophilicity, surface area, nucleophilicity of the oxime oxygen, and location of the molecular orbitals on the isosurfaces have important roles for potencies for reactivating GA-inhibited AChE. The in silico pharmacophore model of oxime affinity for binding to GA-inhibited AChE was found to require a hydrogen bond acceptor, a hydrogen bond donor at the two terminal regions, and an aromatic ring in the central region of the oximes. The model was found to be well-correlated (R = 0.9) with experimental oxime affinity for binding to GA-inhibited AChE. Additional stereoelectronic features relating activity with the location of molecular orbitals and weak electrostatic potential field over the aromatic rings were found to be consistent with the pharmacophore model. These results provided the first predictive pharmacophore model of oxime affinity for binding toward GA-inhibited AChE. The model may be useful for virtual screening of compound libraries to discover and/or custom synthesize more efficacious and less toxic reactivators that may be useful for GA intoxication.

Computational Approaches for the Design of Mosquito Repellent Chemicals.

BACKGROUND: In view of many current mosquito-borne diseases there is a need for the design of novel repellents. OBJECTIVE: The objective of this article is to review the results of the researches carried out by the authors in the computer-assisted design of novel mosquito repellents. METHODS: Two methods in the computational design of repellents have been discussed: a) Quantitative Structure Activity Relationship (QSAR) studies from a set of repellents structurally related to DEET using computed mathematical descriptors, and b) Pharmacophore based modeling for design and discovery of novel repellent compounds including virtual screening of compound databases and synthesis of novel analogues. RESULTS: Effective QSARs could be developed using mathematical structural descriptors. The pharmacophore based method is an effective tool for the discovery of new repellent molecules. CONCLUSION: Results reviewed in this article show that both QSAR and pharmacophore based methods can be used to design novel repellent molecules.

3D-QSAR studies of 2,2-diphenylpropionates to aid discovery of novel potent muscarinic antagonists.

Muscarinic acetylcholine receptors (mAChRs) consisting of five known subtypes, are widely distributed in both central and peripheral nervous systems for regulation of a variety of critical functions. The present theoretical study describes correlations between experimental and calculated molecular properties of 15 alpha-substituted 2,2-diphenylpropionate antimuscarinics using quantum chemical and pharmacophore generation methods to characterize the drug mAChR properties and design new therapeutics. The calculated stereoelectronic properties, such as total energies, bond distances, valence angles, torsion angles, HOMO-LUMO energies, reactivity indices, vibrational frequencies of ether and carbonyl moieties, and nitrogen atom proton affinity were found to be well correlated when compared with experimentally determined inhibition constants from the literature using three muscarinic receptor assays: [(3)H]NMS receptor binding, alpha-amylase release from rat pancreas, and guinea pig ileum contraction. In silico predicted toxicity on rat oral LD(50) values correlated well with the [(3)H]NMS binding in N4TG1 cells and alpha-amylase release assays, but not the ileum contraction assay. Next, to explore the functional requirements for potent activity of the compounds, we developed a preliminary 3D pharmacophore model using the in silico techniques. The resulting model contained a hydrogen bond acceptor site on the carbonyl oxygen atom and a ring aromatic feature on one of the two aromatic rings in these compounds. This model was used as a template to search an in-house database for novel analogs. We found compounds equal in inhibition potency to atropine and, importantly, six not reported before as antimuscarinics. These results demonstrate that this QSAR approach not only provides a basis for understanding the molecular mechanism of action but a pharmacophore to aid in the discovery and design of novel potent muscarinic antagonists.

Computer-Assisted and Data Driven Approaches for Surveillance, Drug Discovery, and Vaccine Design for the Zika Virus.

Human life has been at the edge of catastrophe for millennia due diseases which emerge and reemerge at random. The recent outbreak of the Zika virus (ZIKV) is one such menace that shook the global public health community abruptly. Modern technologies, including computational tools as well as experimental approaches, need to be harnessed fast and effectively in a coordinated manner in order to properly address such challenges. In this paper, based on our earlier research, we have proposed a four-pronged approach to tackle the emerging pathogens like ZIKV: (a) Epidemiological modelling of spread mechanisms of ZIKV; (b) assessment of the public health risk of newly emerging strains of the pathogens by comparing them with existing strains/pathogens using fast computational sequence comparison methods; (c) implementation of vaccine design methods in order to produce a set of probable peptide vaccine candidates for quick synthesis/production and testing in the laboratory; and (d) designing of novel therapeutic molecules and their laboratory testing as well as validation of new drugs or repurposing of drugs for use against ZIKV. For each of these stages, we provide an extensive review of the technical challenges and current state-of-the-art. Further, we outline the future areas of research and discuss how they can work together to proactively combat ZIKV or future emerging pathogens.

Antiviral activities of selected antimalarials against dengue virus type 2 and Zika virus.

In a previous study, twelve antimalarial compounds, amodiaquine (AQ) and derivatives, were shown to have potent anti-dengue viral (DENV) activity by using the stable DENV2 Renilla luciferase reporter replicon expressing BHK-21 cells, infectivity (plaque), and the qRT-PCR assays. In this study, we performed molecular modeling on these compounds to determine their stereo-electronic properties required for optimal antiviral activity. Based on the similarity of calculated stereo-electronic profiles, specifically the electrostatic potential profiles of the compounds, and in silico screening of related compounds from literature, we identified three additional compounds, Quinacrine (QC), Mefloquine (MQ), and GSK369796. Analysis of their antiviral activities indicated that all three compounds have high anti-DENV activity in the DENV2 replicon expressing cells with EC50 values of 5.30 ± 1.31 µM (QC), 3.22 ± 0.37 µM (MQ), and 5.06 ± 0.86 µM (GSK369796). The infectivity assays revealed the EC50 values of 7.09 ± 1.67 µM (QC), 4.36 ± 0.31 µM (MQ) and 3.03 ± 0.35 µM (GSK369796). The mode of action of these compounds is through inhibition of autophagy, thereby affecting DENV2 replication. Moreover, these compounds also showed antiviral activity against the rapidly emerging Zika virus (ZIKV) with EC50 values of 2.27 ± 0.14 µM (QC), 3.95 ± 0.21 µM (MQ), and 2.57 ± 0.09 µM (GSK369796).

Identification of smooth and rough forms in cultures of Brucella melitensis strains by flow cytometry.

Brucella melitensis strains may occur as either smooth or rough variants depending on the expression of O-polysaccharides (OPS) as a component of the bacterial outer membrane lipopolysaccharide (LPS). The wboA gene, which codes for the enzyme glycosyl transferase, is essential for the assembly of O-chain in Brucella. Deletion of wboA in smooth virulent B. melitensis 16M results in a rough mutant designated WRR51. We developed a flow cytometric method to determine the proportion of B. melitensis cells displaying surface O-polysaccharide (OPS) in liquid culture. OPS was detected using polyclonal antibodies from rabbits immunized with smooth (S) or rough (R) Brucella LPS. First, we evaluated the binding of these antibodies to 16M (S), WRR51 (R) and complemented WRR51 expressing the wboA gene (S) as well as to their corresponding GFP-expressing derivative strains 16M/GFP, WRR51/GFP and WRR51/GFP+wboA. The rough mutants did not react with anti-S-LPS nor did the smooth strains react with anti-R-LPS. Second, using different ratios of 16M/GFP and WRR51/GFP, we were able to detect the presence of 1% rough bacteria spiked into a sample of smooth organisms. Third, we evaluated the purity of cultures of B. melitensis strains grown in a fermenter. These flow cytometric methods may be useful for quality control of process development for large-scale vaccine production.

Activity of and initial mechanistic studies on a novel antileishmanial agent identified through in silico pharmacophore development and database searching.

A 3D pharmacophore was generated to describe the antileishmanial activity of dinitroaniline sulfonamides by CATALYST 3D-QSAR methodology, and this pharmacophore was used to search the Maybridge database. Two compounds identified in this search, BTB 06237 and BTB 06256, were highly active with IC(50) values against L. donovani amastigotes of 0.5 +/- 0.2 and 2.3 +/- 0.8 microM, respectively. BTB 06237 also reduced parasite burdens in L. mexicana-infected J774 macrophages at low micromolar concentrations. Unlike the dinitroaniline sulfonamides, the active compounds did not display antimitotic effects against Leishmania. Transmission electron microscopy showed that the single parasite mitochondrion becomes dilated following incubation with BTB 06237, and fluorescence microscopy demonstrated that this organelle fragments into intensely staining spheres when treated with a mitochondrion-specific dye. The mitochondrial membrane potential was also dissipated in BTB 06237-treated parasites. These results indicate that BTB 06237 is an intriguing antileishmanial lead compound that likely interferes with mitochondrial function.

The medicinal chemistry of botulinum, ricin and anthrax toxins.

The potential use of weapons of mass destruction (nuclear, biological or chemical) by terrorist organizations represents a major threat to world peace and safety. Only a limited number of vaccines are available to protect the general population from the medical consequences of these weapons. In addition there are major health concerns associated with a pre-exposure mass vaccination of the general population. To reduce or eliminate the impact of these terrible threats, new drugs must be developed to safely treat individuals exposed to these agents. A review of all therapeutic agents under development for the treatment of the illnesses and injuries that result from exposure to nuclear, biological or chemical warfare agents is beyond the scope of any single article. The intent here is to provide a focused review for medicinal and organic chemists of three widely discussed and easily deployed biological warfare agents, botulinum neurotoxin and ricin toxins and the bacteria Bacillus anthracis. Anthrax will be addressed because of its similarity in both structure and mechanism of catalytic activity with botulinum toxin. The common feature of these three agents is that they exhibit their biological activity via toxin enzymatic hydrolysis of a specific bond in their respective substrate molecules. A brief introduction to the history of each of the biological warfare agents is presented followed by a discussion on the mechanisms of action of each at the molecular level, and a review of current potential inhibitors under investigation.

Analysis of molecular stereoelectronic similarity between N,N-diethyl-m-toluamide (DEET) analogs and insect juvenile hormone to develop a model pharmacophore for insect repellent activity.

Similarity analysis on molecular stereoelectronic properties of N,N-diethyl-m-toluamide (DEET), natural insect juvenile hormone (JH), a synthetic insect juvenile hormone mimic (JH-mimic, undecen-2-yl carbamate), and DEET compounds reveals remarkable similarities that lead to a reliable pharmacophore for the design of efficacious insect repellents and provide insights for understanding the mechanism of repellent action. The study involves an AM1 quantum chemical computational procedure enabling a conformational search for the lowest and most abundant energy conformers of JH, JH-mimic, and 15 DEET compounds and complete geometry optimization of the conformers. Similarity analyses of stereoelectronic properties such as structural parameters, atomic charges, dipole moments, molecular electrostatic potentials, and highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies were performed on JH, JH-mimic, and the DEET compounds. Similarity of stereoelectronic attributes of the amide/ester moiety, negative electrostatic potential regions beyond the molecular surface, and a large distribution of hydrophobic regions in the compounds appears to be the 3 important factors leading to a similar interaction with the JH receptor. The similarity of electrostatic profiles beyond the molecular surface is likely to play a crucial role toward molecular recognition interaction with the JH receptor from a distance which suggests a possible electrostatic bioisosterism of the amide group of the DEET compounds and JH-mimic and, thus, a model for molecular recognition at the JH receptor.

Discovery of non-oxime reactivators using an in silico pharmacophore model of reactivators for DFP-inhibited acetylcholinesterase.

Utilizing our previously reported in silico pharmacophore model for reactivation efficacy of oximes, we present here a discovery of twelve new non-oxime reactivators of diisopropylfluorophosphate (DFP)-inhibited acetylcholinesterase (AChE) obtained through virtual screening of an in-house compound database. Rate constant (kr) efficacy values of the non-oximes were found to be within ten-fold of pralidoxime (2-PAM) in an in vitro DFP inhibited eel AChE assay and one of them showed in vivo efficacy comparable to 2-PAM against brain symptoms for DFP induced neuropathology in guinea pigs. Short listing of the identified compounds were performed on the basis of in silico evaluations for favorable blood brain barrier penetrability, octanol-water partition (Clog P), toxicity (rat oral LD50) and binding affinity to the active site of the crystal structure of a OP- inhibited AChE.

The anthrax protective antigen (PA63) bound conformation of a peptide inhibitor of the binding of lethal factor to PA63: as determined by trNOESY NMR and molecular modeling.

Anthrax protective antigen (PA) is one of the three proteins produced by the gram positive bacteria Bacillus anthracis collectively known as the "anthrax toxin" (Ascenzi, P.; Visca, P.; Ippolito, G.; Spallarossa, A.; Bolognesi, M.; et al. Anthrax toxin: a tripartite lethal combination. FEBS Lett. 2002, 531, 384-388). The role played by PA in anthrax intoxication is to transport the two enzymes lethal factor (LF) and edema factor (EF) into the cell. Collier and co-workers (Mourez, M.; Kane, R. S.; Mogridge, J.; Metallo, S.; Deschatelets, P.; et al. Designing a polyvalent inhibitor of anthrax toxin. Nat. Biotechnol. 2001, 958). reported the isolation of two peptides via phage display that bind to the PA63 heptamer and inhibit its interaction with LF and EF, and thereby prevent the transport of LF and EF into the cell. One of these peptides, His-Thr-Ser-Thr-Try-Trp-Trp-Leu-Asp-Gly-Ala-Pro (P1), was selected for structural investigation on the basis of its ability to prevent the binding of LF to the PA63 heptamer bundle. Two-dimensional trNOESY experiments coupled with NOE restrained simulated annealing calculations were used to determine the PA63-bound conformation of P1. On binding to PA63, P1 adopts a helical conformation involving residues 3-9 while the C- and N-terminal residues exhibit dynamic fraying.

Carbon isosteres of the 4-aminopyridine substructure of chloroquine: effects on pK(a), hematin binding, inhibition of hemozoin formation, and parasite growth.

Unlike diprotic chloroquine (CQ), its two 4-aminoquinoline carbon isosteres (1, 2) are monoprotic at physiological pH. Compared to CQ, hematin binding affinity of 1 decreased 6.4-fold, and there was no measurable binding for 2. Although 1 was a weak inhibitor of hemozoin formation, neither isostere inhibited P. falciparum in vitro. Evidently, the CQ-hematin interaction is largely a function of its pyridine substructure, but inhibition of hemozoin formation and parasite growth depends on its 4-aminopyridine substructure.

A three-dimensional in silico pharmacophore model for inhibition of Plasmodium falciparum cyclin-dependent kinases and discovery of different classes of novel Pfmrk specific inhibitors.

The cell division cycle is regulated by a family of cyclin-dependent protein kinases (CDKs) that are functionally conserved among many eukaryotic species. The characterization of plasmodial CDKs has identified them as a leading antimalarial drug target in our laboratory. We have developed a three-dimensional QSAR pharmacophore model for inhibition of a Plasmodium falciparum CDK, known as Pfmrk, from a set of fifteen structurally diverse kinase inhibitors with a wide range of activity. The model was found to contain two hydrogen bond acceptor functions and two hydrophobic sites including one aromatic-ring hydrophobic site. Although the model was not developed from X-ray structural analysis of the known CDK2 structure, it is consistent with the structure-functional requirements for binding of the CDK inhibitors in the ATP binding pocket. Using the model as a template, a search of the in-house three-dimensional multiconformer database resulted in the discovery of sixteen potent Pfmrk inhibitors. The predicted inhibitory activities of some of these Pfmrk inhibitors from the molecular model agree exceptionally well with the experimental inhibitory values from the in vitro CDK assay.

Oxindole-based compounds are selective inhibitors of Plasmodium falciparum cyclin dependent protein kinases.

Cyclin dependent protein kinases (CDKs) have become attractive drug targets in an effort to identify effective inhibitors of the parasite Plasmodium falciparum, the causative agent of the most severe form of human malaria. We tested known CDK inhibitors for their ability to inhibit two malarial CDKs: Pfmrk and PfPK5. Many broad spectrum CDK inhibitors failed to inhibit Pfmrk suggesting that the active site differs from other CDKs in important ways. By screening compounds in the Walter Reed chemical database, we identified oxindole-based compounds as effective inhibitors of Pfmrk (IC(50) = 1.5 microM). These compounds have low cross-reactivity against PfPK5 and human CDK1 demonstrating selectivity for Pfmrk. Amino acid comparison of the active sites of Pfmrk and PfPK5 identified unique amino acid differences that may explain this selectivity and be exploited for further drug development efforts.

Inhibition of cholinesterases with cationic phosphonyl oximes highlights distinctive properties of the charged pyridine groups of quaternary oxime reactivators.

Oxime-induced reactivation of phosphonylated cholinesterases (ChEs) produces charged phosphonyl pyridine oxime intermediates (POXs) that are most potent organophosphate (OP) inhibitors of ChEs. To understand the role of cationic pyridine oxime leaving groups in the enhanced anti-ChE activity of POXs, the bimolecular rate constants for the inhibition (k(i)) of acetylcholinesterases (AChE) and butyrylcholinesterases (BChE), and the rate of decomposition (k(d)) of authentic O-alkyl methylphosphonyl pyridine oximes (AlkMeP-POXs) and N,N-dimethylamidophosphoryl pyridine oximes (EDMP-POXs), were studied. Stability ranking order in aqueous solutions correlated well with the electronic features and optimized geometries that were obtained by ab initio calculations at 6-31G(**) basis set level. AlkMeP-POXs of the 2-pyridine oxime series were found to be 4- to 8-fold more stable (t(1/2)=0.7 to 1.5 min) than the homologous O,O-diethylphosphoryl (DEP) oxime. Results suggest that re-inhibition of enzyme activity by POX is less likely during the reactivation of DEP-ChEs (obtained by use of DEP-containing pesticides) by certain oximes, compared to nerve agent-inhibited ChEs. The greatest inhibition was observed for the O-cyclohexyl methylphosphonyl-2PAM derivative (4.0 x 10(9)M(-1)min(-1); mouse AChE) and is 10-fold higher than the k(i) of cyclosarin. Increasing the size of the O-alkyl substituent of AlkMeP-POXs had only a small to moderate effect on the k(i) of ChEs, signifying a major role for the cationic pyridine oxime leaving group in the inhibition reaction. The shape of plots of logk(i) vs. pK(a) of the leaving groups for AlkMeP-PAMs and DEP-PAMs, could be used as a diagnostic tool to highlight and rationalize the unique properties of the cationic moiety of pyridine oxime reactivators.

Nuclear magnetic resonance and molecular modeling analysis of the interaction of the antimalarial drugs artelinic acid and artesunic acid with beta-cyclodextrin.

The artemisinin derivatives artelinic acid and artesunic acid are members of a class of compounds that have shown promise for the treatment of multidrug resistant strains of Plasmodium falciparum. Unfortunately, these compounds exhibit poor solubility and stability in aqueous solution. The research presented herein was conducted to determine whether complexation of artelinic acid or artesunic acid with beta-cyclodextrin would result in complexes with increased aqueous solubility while retaining the potent antimalarial activity of these compounds. Preliminary complexation studies with natural beta-cyclodextrins were conducted as a proof of concept, with a primary focus on understanding the electrostatic interactions that stabilize the resulting complexes. Complex formation was monitored using UV spectroscopy. The structures of the resulting complexes were determined using multidimensional nuclear magnetic resonance spectroscopy (NMR) and molecular modeling. NMR results are most consistent for artelinic acid and beta-cyclodextrin forming complexes in a ratio of 2:1; however, the presence of 1:1, 2:2, and 3:1 complexes in solution cannot be excluded based on the experimental data collected. The NMR data also indicate selective insertion of artelinic acid into the hydrophobic cavity of the beta-cyclodextrin via the primary face. NMR results indicate artesunic acid forms a similar complex with beta-cyclodextrin in a ratio of 1:1; again however, the presence of 1:1, 2:2, and 3:1 complexes in solution cannot be ruled out.

Structure-activity relationship study of antimalarial indolo [2,1-b]quinazoline-6,12-diones (tryptanthrins). Three dimensional pharmacophore modeling and identification of new antimalarial candidates.

A widely applicable three-dimensional QSAR pharmacophore model for antimalarial activity was developed from a set of 17 substituted antimalarial indolo[2,1-b]quinazoline-6,12-diones (tryptanthrins) that exhibited remarkable in vitro activity (below 100 ng/mL) against sensitive and multidrug-resistant Plasmodium falciparum malaria. The pharmacophore, which contains two hydrogen bond acceptors (lipid) and two hydrophobic (aromatic) features, was found to map well onto many well-known antimalarial drug classes including quinolines, chalcones, rhodamine dyes, Pfmrk cyclin dependent kinase inhibitors, malarial FabH inhibitors, and plasmepsin inhibitors. The phamacophore allowed searches for new antimalarial candidates from multiconformer 3D databases and enabled custom designed synthesis of new potent analogues.