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.

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).

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 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.

In silico stereo-electronic analysis of PMD (p-Menthane-3-8-Diol) and its derivatives for pharmacophore development may aid discovery of novel insect repellents.

PMD (p-menthane-3-8-diol) is an insect repellent that can be synthesized chemically or derived from a steam distillate residue of the leaves of lemon eucalyptus, Corymbia citriodora. It is one of the few natural product endorsed by the Center for Disease Control (USA) for topical application to protect against mosquitoes though it is not as effective as the common repellent DEET (N,N -diethyl-1,3-toluamide). However, DEET has several undesirable side effects and toxicity too. Thus, although PMDs are comparatively safer than DEET, no quantitative structure activity relationship (QSAR) and pharmacophore modeling studies have been reported in literature to improve efficacy and aid further development of more effective PMD analogues. In this study, we report results of quantum chemical analysis of stereoelectronic properties and pharmacophore modeling of PMD and eight of its synthetic derivatives to aid discovery and design of more effective PMD analogues. Stereo-electronic analysis indicates that lower aqueous stabilization (favorable lipophilicity) and larger separation of electrostatic potential energy together with a large localized negative electrostatic potential region by the oxygen atom play important roles for repellent activity. Consistent to these properties, the generated pharmacophore model of the PMDs showed two aliphatic hydrophobic and a hydrogen-bond donor features for potent activity. These results aided us to design more effective PMD repellents which are currently under further investigations.

In silico pharmacophore modeling on known pyridinium oxime reactivators of cyclosarin (GF) inhibited AChE to Aid discovery of potential, more efficacious novel non-oxime reactivators.

Cyclohexyl methylphosphonofluoridate (cyclosarin, cyclosin, GF) is a highly toxic organophosphorus (OP) nerve agent considered as potential warfare threats and known to be resistant to conventional oxime antidotal therapy. To aid discovery of novel antidotes for GF toxicity, a three-dimensional in silico pharmacophore model for reactivation efficacy against GF intoxication is presented. The model was generated from published experimental percentage reactivation data on oximes as changes of AChE/BuChE activities in the whole blood after cyclosarin intoxication and administration. The generated pharmacophore model was found to contain a hydrogen bond donor site and two ring aromatic sites as necessary optimal features for reactivation of GF intoxication. Stereo-electronic features of oximes reported by us earlier provided guidance to develop the model and were found to be consistent with the reported structure activity data. Furthermore, from virtual screening of two commercial databases, Maybridge and ChemNavigator using map-fitting of the model led us to identify two new non-oxime compounds showing reactivation efficacy within 10-fold range of 2-PAM for DFP-inhibited AChE. Since GF is a G simulator like DFP (diisopropylfluorophosphate), the model should have the potential for discovery of novel reactivators against GF intoxication.

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.

An in silico stereo-electronic comparison of conventional pyridinium oximes and K-oximes for organophosphate (OP) poisoning.

A comparative analysis of stereo-electronic properties of five cholinesterase reactivators (pralidoxime (2- PAM), trimedoxime, obidoxime, HI-6, and HLo-7) and six "K-oximes" was performed to assess their roles in reactivating OP-inhibited phosphorylated serine residue of mouse AChE. Quantum mechanical (QM) calculations starting from semiempirical to ab initio levels were sequentially performed with hierarchical basis sets to obtain the individual optimized geometry and stereo-electronic properties of the eleven oximes. Next, solvation effects were computed on the optimized structures using two different (PCM and COSMO) QM models. Results indicate that properties, such as the distance between the bisquarternary nitrogen atoms, surface area, molecular volume, and hydrophilicity have important roles in the reactivation of OP-inhibited AChE. Electronic attributes, such as the molecular electrostatic potentials and orbital energies were also found to be important parameters for reactivation. Nucleophilicity of the oxygen atoms at the terminal regions, electrophilicity in the central regions of the oximes, and location of the molecular orbitals on aromatic rings have significant roles for the experimentally observed reactivations in several OP agents inhibited mouse AChE. Analysis of solvation free energy indicates high solute polarization and dispersion energies of the oximes to be particularly critical for the tabun- inhibited mouse AChE, whereas lower values of these properties favor reactivation against other OP agents, such as soman, sarin and cyclosarin. Feature mappings of our recently reported pharmacophore model were also observed to be consistent with the above observed electronic properties. In silico toxicity evaluation on these oximes predicts the Koximes to have somewhat higher oral toxicity compared to the other bispyridinium oximes.

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

We earlier reported an in silico pharmacophore model for reactivation of oximes to tabun-inhibited AChE. Since DFP (diisopropylfluorophosphate) like tabun is a G-agent simulator, we utilized the model as a rational strategy to discover non-oxime reactivators of DFP-inhibited AChE in this study. The phramacophore was used for virtual screening of two commercial databases, Maybridge and ChemNavigator, to identify reactivators which lack the oxime functions. The procedure led us to identify several potent non-oxime compounds that reactivate DFP-inhibited AChE. These non-oxime reactivators contain a nucleophile group in lieu of the oxime moiety in the compound. Five of these novel non-oximes showed Kr values within ten-fold of 2-PAM in an in vitro assay. The pharmacophore model contained a hydrogen bond acceptor, a hydrogen bond donor, and an aromatic ring features distributed in a 3D space. Calculated stereoelectronic properties reported earlier with respect to the location of molecular orbitals and electrostatic potentials were consistent with the model and the newly identified compounds. Down selection of compounds after virtual screening was performed on the basis of fit score to the model, conformational energy, and in silico evaluations for favorable blood-brain barrier (BBB) penetrability, octanol-water partition (log P), and toxicity (rat oral LD(50)) assessments. In vitro reactivation efficacy of the compounds was evaluated in a DFP-inhibited eel acetylcholinesterase assay.

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.

Enhanced antibody responses to a detoxified lipopolysaccharide-group B meningococcal outer membrane protein vaccine are due to synergistic engagement of Toll-like receptors.

When given passively or elicited actively, antibodies induced by a detoxified Escherichia coli Rc chemotype (J5) mutant lipopolysaccharide (J5dLPS)-group B meningococcal outer membrane protein (OMP) complex vaccine protected animals from lethal sepsis. The protection from sepsis is believed to be dependent on high levels of antibodies against the core glycolipid (CGL), a region of LPS that is rather conserved among Enterobacteriaceae. The addition of unmethylated deoxycytidyl-deoxyguanosine dinucleotide (CpG)-containing oligodeoxynucleotides (ODN) was used as an immuno-adjuvant to improve antibody responses. In preparation for a Phase I human trial, we elucidated potential contributions by which the sepsis vaccine (J5dLPS-OMP) and CpG ODN might enhance the antibody response and provide evidence that the generation of immune responses is Toll-like receptor (TLR) dependent. Toll-like receptor 2, TLR4, and TLR9 were each essential for generating robust cytokine and antibody responses. The signature cytokine of dendritic cells, interleukin-12, was one of the cytokines that demonstrated synergy with the optimal TLR ligand/ engagement combination. We conclude that the involvement of multiple TLRs upon immunization was critical for the generation of optimal antibody responses. These observations provide further evidence for the inclusion of innate immune-based adjuvants during the development of next-generation vaccines.

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.

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.

Intranasal administration of a detoxified endotoxin vaccine protects mice against heterologous Gram-negative bacterial pneumonia.

When given passively or elicited actively, antibodies induced by a detoxified Escherichia coli J5 mutant lipopolysaccharide (J5dLPS)-group B meningococcal outer membrane protein (-OMP) vaccine previously protected animals from lethal sepsis. To assess the use of this vaccine for the treatment of Gram-negative bacillary pneumonia, we vaccinated mice, with or without the adjuvant CpG, by intranasal (i.n.) or intraperitoneal (i.p.) routes of administration. Local and systemic IgG levels were 2-3 logs higher following i.p. immunization compared to i.n. However, i.n. immunization elicited both local and systemic IgA, unlike i.p. administration. The addition of CpG to the vaccine, by either route of administration, elicited greater levels of antibody. Intranasal immunization protected mice against lethal heterologous Gram-negative bacillary pneumonia and post-immunization serum and broncho-alveolar lavage fluid mediated enhanced bacterial killing with peritoneal and alveolar macrophages in vitro. We conclude that further studies on the use of J5dLPS-OMP for the prevention of nosocomial pneumonia are warranted.

Novel semi-automated methodology for developing highly predictive QSAR models: application for development of QSAR models for insect repellent amides.

Conventional 3D-QSAR models are built using global minimum conformations or quantum-mechanics based geometry-optimized conformations as bioactive conformers. QSAR models developed using the global minima as bioactive conformers, employing the GFA, PLS and G/PLS methodologies, gave good non-validated r(2) (0.898, 0.868 and 0.922) and performed well on an internal validation test with leave-one-out correlation q(2) (LOO) (0.902, 0.726 and 0.924), leave-10%-out correlation q(2) (L10O) (0.874, 0.728 and 0.883) and leave-20%-out q(2) (L20O) (0.811, 0.716 and 0.907). However, they showed poor predictive ability on an external data set with best predictive r(2) (Pred-r(2)) of 0.349, 0.139 and 0.204 respectively. A novel methodology to mine bioactive conformers, from clusters of conformations with good 3D-spatial representation around pharmacophoric moiety, furnishes highly predictive 3D-QSAR models. The best QSAR model (model A) showed r(2) of 0.989, q(2) (LOO) of 0.989, q(2) (L10O) of 0.980, q(2) (L20O) of 0.963 and Pred-r(2) on eight test compounds of 0.845. The methodology is based on mimicking the multi-way Partial Least Squares (PLS) technique by performing several automated sequential PLS analyses. The poses/shapes of the mined bioactive conformers provide valuable insight into the mechanism of action of the insect repellents. All of the repetitive tasks were automated using Tcl-based Cerius2 scripts.