In silico docking and molecular dynamics simulation of 3-dehydroquinate synthase (DHQS) from <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis>

The shikimate pathway is as an attractive target because it is present in bacteria, algae, fungi, and plants but does not occur in mammals. In Mycobacterium tuberculosis (MTB), the shikimate pathway is integral to the biosynthesis of naphthoquinones, menaquinones, and mycobactin. In these study, novel inhibitors of 3-dehydroquinate synthase (DHQS), an enzyme that catalyzes the second step of the shikimate pathway in MTB, were determined. 12,165 compounds were selected from two public databases through virtual screening and molecular docking analysis using PyRx 8.0 and Autodock 4.2, respectively. A total of 18 compounds with the best binding energies (?13.23 to ?8.22 kcal/mol) were then selected and screened for absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis, and nine of those compounds were found to satisfy all of the ADME and toxicity criteria. Among those nine, the three compounds—ZINC633887?(binding energy =??10.29 kcal/mol), ZINC08983432?(?9.34 kcal/mol), and PubChem73393?(?8.61 kcal/mol)—with the best binding energies were further selected for molecular dynamics (MD) simulation analysis. The results of the 50-ns MD simulations showed that the two compounds ZINC633887 and PubChem73393 formed stable complexes with DHQS and that the structures of those two ligands remained largely unchanged at the ligand-binding site during the simulations. These two compounds identified through docking and MD simulation are potential candidates for the treatment of TB, and should undergo validation in vivo and in vitro.     

Virtual screening for novel COX-2 inhibitors using the ZINC database

Cyclooxygenase-2 (COX-2) enzyme binds to arachidonic acid and releases metabolites that are used to induce pain and inflammation. COX-2 selective inhibitors such as celecoxib, rofecoxib and valdecoxib are currently used to reduce inflammatory response. However, they lack anti-thrombotic activity and hence lead to cardiovascular and renal liabilities apart from gastrointestinal irritation. Therefore, there is still a need to develop more potent COX-2 inhibitors. In this paper, we report the screening of various compounds from the ZINC database (contains 4.6 million small molecule compounds) using the eHiTS (electronic High Throughput Screening) software tool against the COX-2 protein. The strategy employed can be conveniently split into two categories, viz. screening and docking, respectively. Screening was performed using molecular constraints tool to filter compounds with physico-chemical properties similar to the 6COX bound ligand SC-558. The analysis resulted in 1042 Lipinski compliant hits which are docked and scored to identify structurally novel ligands that make similar interactions to those of known ligands or may have different interactions with other parts of the binding site. Our screening approach identified two molecules ZINC00663976 (eHITS score of -7.135 kcal/mol) and ZINC02062094 (eHITS score of -7.242 kcal/mol) from the ZINC database. Their energy scores are better than the 6COX bound co-crystallized ligand SC-558 with an eHiTS score of -6.559 kcal/mol. Both the ligands were docked within the binding pocket forming interactions with Leu352, Phe518, Met522, Val523, Ala527 and Ser353. Visual inspection suggested similar orientation and binding mode for ZINC02062094 with SC-558 ligand. The NH group of the ligand formed hydrogen bond interactions with the backbone NH of Ala527.     

Computational finding of potential inhibitor for Cytochrome P450 Mono-oxygenases Enzyme of Mycobacterium tuberculosis

Cytochrome P450 mono-oxygenases (2UUQ) enzyme from Mycobacterium tuberculosis catalyzes oxidation of organic compounds such as lipids and steroidal hormones therefore remain as potential drug target. Currently available first line anti-tuberculosis drugs have been caused several side effects in the body as well as resistance development by mycobacterium against these drugs, necessitates the considerable need for finding new drugs. Therefore, we propose a structure based computational method to find a new potential inhibitor for cytochrome P450 mono-oxygenases enzyme. Compounds from several ligand databases were docked against the functional sites of 2UUQ (A) through the standard GEMDOCK v2.0 and AUTODOCK4.0 molecular docking tools. Commercially available chemical compound ZINC00004165 (5-[3-(2-nitroimidazol-1-yl) propyl] phenanthridine) has produced top rank with lowest interaction energy of -113.2 (via GEMDOCK) and lowest docking energy of -9.80 kcal/mol (via AUTODOCK) as compared to first line anti TB compounds. Z score and normal distribution analysis verified that the ZINC00004165 compound has more affinity towards 2UUQ in comparison to large number of random population of compounds. ZINC00004165 is also in agreement with the drug likeness properties of Lipinski rule of five without any violation. Therefore, our finding concludes that the commercial compound ZINC00004165 can act as a potential inhibitor against cytochrome P450 mono-oxygenases enzyme of Mycobacterium tuberculosis.     

Molecular docking analysis of HER-2 inhibitor from the ZINC database as anticancer agents

The human epidermal growth factor (HER2) is a transmembrane receptor that is highly expressed in breast cancer and in different other cancers. Therefore, it is of interest to identify the new HER2 inhibitors from a selected 300 compounds in the ZINC database. The top two hit compounds (ZINC000014780728 (-11.0 kcal/mol) and ZINC000014762512 (-10.8 kcal/mol)) showed a high affinity with HER2 relative to the reference compound (lapatinib (-10.2 kcal/mol)) for further consideration.     

In-silico analysis of Sirt2 from Schistosoma mansoni: structures,conformations, and interactions with inhibitors

Sirtuins are NAD+-dependent lysine deacetylases member of the class III HDAC family. These are demonstrated to be therapeutic targets in parasitic diseases like schistosomiasis. Observations suggested that sirtuin enzyme is necessary for the functionality of fe/male reproductive system, due to which SmSirt2 is treated as a potential therapeutic target. There are no structural and molecular features of SmSirt2 have been reported yet. In this study, homology modeling has been used to determine the three-dimensional features of the SmSITRT2. Further, structure validation has been performed by energy minimization and Ramachandran plot. Validated structures are further subjected to molecular docking and virtual screening to find the best lead molecules for downstream analysis. Ten lead molecules were selected while comparing virtual screening of hSirt2 and SmSirt2 both. These leads are further compared with AKG2 which is known inhibitor of hSirt2 (?8.8 kcal/mol). Out of selected 10 leads, four of them (ZINC23995485 (?9.5 kcal/mol), ZINC53298162 (?9.4 kcal/mol), ZINC70927268 (?10.0 kcal/mol), ZINC89878705 (?11.2 kcal/mol)) have shown better interaction with SmSirt2, in which ZINC89878705 (?11.2 kcal/mol) shows a more compact packing as compared to AKG2 and rest of ligands. These molecules could be further subject to in vitro study and model of SmSirt2 has been proposed for further structure-based drug design projects concerning sirtuins from Schistosoma mansoni.     

Virtual screening using the ligand ZINC database for novel lipoxygenase-3 inhibitors

The leukotrienes constitute a group of arachidonic acid-derived compounds with biologic activities suggesting important roles in inflammation and immediate hypersensitivity. Epidermis-type lipoxygenase-3 (ALOXE3), a distinct subclass within the multigene family of mammalian lipoxygenases, is a novel isoenzyme involved in the metabolism of leukotrienes and plays a very important role in skin barrier functions. Lipoxygenase selective inhibitors such as azelastine and zileuton are currently used to reduce inflammatory response. Nausea, pharyngolaryngeal pain, headache, nasal burning and somnolence are the most frequently reported adverse effects of these drugs. Therefore, there is still a need to develop more potent lipoxygenase inhibitors. In this paper, we report the screening of various compounds from the ZINC database (contains over 21 million compounds) using the Molegro Virtual Docker software against the ALOXE3 protein. Screening was performed using molecular constraints tool to filter compounds with physico-chemical properties similar to the 1N8Q bound ligand protocatechuic acid. The analysis resulted in 4319 Lipinski compliant hits which are docked and scored to identify structurally novel ligands that make similar interactions to those of known ligands or may have different interactions with other parts of the binding site. Our screening approach identified four molecules ZINC84299674; ZINC76643455; ZINC84299122 & ZINC75626957 with MolDock score of -128.901, -120.22, -116.873 & - 102.116 kcal/mol, respectively. Their energy scores were better than the 1N8Q bound co-crystallized ligand protocatechuic acid (with MolDock score of -77.225 kcal/mol). All the ligands were docked within the binding pocket forming interactions with amino acid residues.     

In silico analyses of major active constituents of fingerroot (Boesenbergia rotunda) unveils inhibitory activities against SARS-CoV-2 main protease enzyme

Boesenbergia rotunda (L.) Mansf., commonly known as fingerroot is a perennial herb in the Zingiberaceae family with anticancer, anti-leptospiral, anti-inflammatory, antioxidant, antiulcer, and anti-herpes viral activities. While the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) inhibitory activity of B. rotunda extract has been recently found, the active compounds contributing to this activity are yet unknown. The main protease (M^pro) enzyme is one of the most well established therapeutic targets among coronaviruses which plays a vital role in the maturation and cleavage of polyproteins during viral replication. The current work aims to identify active phytochemical substances from B. rotunda extract that can inhibit the replication of SARS-CoV-2 by using a combined molecular docking and dynamic simulation approaches. The virtual screening experiment revealed that fifteen molecules out of twenty-three major active compounds in the plant extract have acceptable drug-like characteristics. Alpinetin, Pinocembrin, and Pinostrobin have binding energies of ?7.51 kcal/mol, ?7.21 kcal/mol, and ?7.18 kcal/mol, respectively, and can suppress M^pro activity. The stability of the simulated complexes of the lead compounds with the drug-receptor is demonstrated by 100-ns MD simulations. The binding free energies study utilizing molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA) show that the compounds and M^pro enzyme have favourable thermodynamic interactions, which are majorly driven by van der Waals forces. Thus, the selected bioactive phytochemicals from B. rotunda might be used as anti-SARS-CoV-2 candidates that target the M^pro enzyme.     

Molecular docking analysis of known flavonoids as duel COX-2 inhibitors in the context of cancer

Cyclooxygenase-2 (COX-2) catalyzed synthesis of prostaglandin E2 and it associates with tumor growth, infiltration, and metastasis in preclinical experiments. Known inhibitors against COX-2 exhibit toxicity. Therefore, it is of interest to screen natural compounds like flavanoids against COX-2. Molecular docking using 12 known flavanoids against COX-2 by FlexX and of ArgusLab were performed. All compounds showed a favourable binding energy of >-10 KJ/mol in FlexX and > -8 kcal/mol in ArgusLab. However, this data requires in vitro and in vivo verification for further consideration.     

Identifications small molecules inhibitor of p53-mortalin complex for cancer drug using virtual screening

Mortalin was over expressed in tumor cells and bind to p53 protein. This interaction was suggested to promote sequestration of p53 in the cytoplasm, thereby inhibiting its nuclear activity. The p53 is a tumor suppressor that is essential for the prevention of cancer development and loss of p53 function is one of the early events in immortalization of human cells. Therefore, abrogation p53-mortalin interaction using small molecule is guaranteed stop cancer cell grow. However study interaction of p53-mortalin, and its inhibition using small molecule is still challenging because specific site of mortalin that bind to p53, vice versa, is still debatable. This study has aims to analyze the p53-binding site of mortalin using molecular docking and to screen drug-like compounds that have potential as inhibitors of p53-mortalin interaction using virtual screening. The result showed that the lowest energy binding of p53-mortalin complex is -31.89 kcal/mol, and p53 protein bind to substrate binding domain of mortalin (THR433; VAL435; LEU436; LEU437; PRO442; ILE558; LYS555). Furthermore, the p53-binding domain of mortalin was used as receptor to screen 9000 drug-like compounds from ZINC database using molecular docking program Auto Dock Vina in PyRx 0.8 (Virtual Screening Tools). Here, we have identified three drug-like compounds that are ZINC01019934, ZINC00624418 and ZINC00664532 adequate to interrupt stability of p53-mortalin complex that warrant for anticancer agent.     

Anti-methanogenic effect of rhubarb (Rheum spp.) – An in silico docking studies on methyl-coenzyme M reductase (MCR)

The present study explored anti-methanogenic properties of rhubarb compounds using in silico analysis on methyl-coenzyme M reductase (MCR) for identifying its anti-methanogen mechanism. To identify pharmacokinetics of 35 compounds from rhubarb, molecular docking and ADME analysis were performed against MCR using AutoDockVina, FAFDrugs3 and PROTOX programs. Docking results successfully indicated three possible candidate compounds 9,10-anthracenedione, 1,8-dihydroxy-3-methyl (?6.92 kcal/mol); phthalic acid isobutyl octadecyl ester (?5.26 kcal/mol); and diisooctyl phthalate (?5.61 kcal/mol) showed minimum binding energy (kcal/mol) with the target protein MCR which catalyze the biosynthesis of rumen methane. In conclusion, the identified compounds showed the most docking fitness score against the target methyl-coenzyme M reductase and the decrease in ruminal methane emission by rhubarb might be a result of these compounds by inhibition of methanogenesis.     

Virtual screening to identify <Emphasis Type="Italic">Leishmania braziliensis N</Emphasis>-myristoyltransferase inhibitors: pharmacophore models,docking, and molecular dynamics

Leishmaniasis is caused by several protozoa species belonging to genus Leishmania that are hosted by humans and other mammals. Millions of new cases are recorded every year and the drugs available on the market do not show satisfactory efficacy and safety. A hierarchical virtual screening approach based on the pharmacophore model, molecular docking, and molecular dynamics was conducted to identify possible Leishmania braziliensis N-misristoyltransferase (LbNMT) inhibitors. The adopted pharmacophore model had three main features: four hydrophobic centers, four hydrogen-bond acceptor atoms, and one positive nitrogen center. The molecules (n=15,000) were submitted to alignment with the pharmacophore model and only 27 molecules aligned to model. Six molecules were submitted to molecular docking, using receptor PDB ID 5A27. After docking, the ZINC35426134 was a top-ranked molecule (? 64.61 kcal/mol). The molecule ZINC35426134 shows hydrophobic interactions with Phe82, Tyr209, Val370, and Leu391 and hydrogen bonds with Asn159, Tyr318, and Val370. Molecular dynamics simulations were performed with the protein in its APO and HOLO forms for 37 ns in order to assess the stability of the protein–ligand complex. Results showed that the HOLO form was more stable than the APO one, and it suggests that the ZINC35426134 binding stabilizes the enzyme. Therefore, the selected molecule has the potential to meet the herein proposed target.     

Molecular modeling of Ruellia tuberosa L compounds as a-amylase inhibitor: an in silico comparation between human and rat enzyme model

Inhibition of α-amylase is an important strategy to control post-prandial hyperglycemia. The present study on Ruellia tuberosa, known as traditional anti-diabetic agent, is being provided in silico study to identify compounds inhibiting α-amylase in rat and human. Compounds were explored from PubChem database. Molecular docking was studied using the autodock4. The interactions were further visualized and analyzed using the Accelrys Discovery Studio version 3.5. Binding energy of compounds to α-amylase was varying between -1.92 to -6.66 kcal/mol in rat pancreatic alpha amylase and -3.06 to -8.42kcal/mol in human pancreatic alpha amylase, and inhibition konstanta (ki) was varying between 13.12- 39460µM in rat and 0.67-5600µM in human. The docking results verify that betulin is the most potential inhibitor of all towards rat model alpha amylase and human alpha amylase. Further analysis reveals that betulin could be a potential inhibitor with non-competitive pattern like betulinic acid. In comparison, betulin has smaller Ki (0.67µM) than acarbose (2.6 µM), which suggesting that betulin is more potential as inhibitor than acarbose, but this assumption must be verified in vitro.     

Comparative docking of dual conformations in human fatty acid synthase thioesterase domain reveals potential binding cavity for virtual screening of ligands

Human fatty acid synthase (hFASN), a homo dimeric lipogenic enzyme with seven catalytic domains, is an important clinical target in cancer, metabolic syndrome and infections. Here, molecular modelling and docking methods were implemented to examine the inter-molecular interactions of thioesterase (TE) domain in hFASN with its physiological substrate, and to identify potential chemical inhibitors. TE catalyses the hydrolysis of thioester bond between palmitate and the 4’ phosphopantetheine of acyl carrier protein, releasing 16-carbon palmitate. The crystal structure of hFASN TE in two inhibitory conformations (A and B) were geometry-optimized and used for molecular docking with palmitate, orlistat (a known FASN inhibitor) and virtual screening against compounds from National Cancer Institute (NCI) database. Relatively, low binding affinity was observed during the complex formation of palmitate with A (?.164 kcal/mol) and B (?.332 kcal/mol) forms of TE, when compared with orlistat-docked TE (A form: ?5.872 kcal/mol and B form: ?5.484 kcal/mol), clearly indicating that the native inhibited conformation (crystal structure) was unfavourable for substrate binding. We used these orlistat dual binding modes as positive controls for prioritizing the ligands during virtual screening. From 2, 31,617 molecules in the NCI database, 916 high-scoring compounds (hit ligands) were obtained for A-form and 4582 for B-form of the TE-domain, which were then ranked according to glide docking score, XP H bond score, absorption, distribution, metabolism and excretion and binding free energy (Prime/MM-GBSA). Consequently, two top scoring ligands (NSC: 319661 and NSC: 153166) emerged as promising drug candidates that may be tested in FASN-over-expressing diseases.     

A molecular modeling based screening for potential inhibitors to alpha hemolysin from Staphylococcus aureus

Staphylococcus aureus, a Gram-positive bacterium is pathogenic in nature. It is known that secreted toxins remain active after antibiotic treatment. The alpha hemolysin or alpha toxin damages cell membrane and induces apoptosis and degradation of DNA. The titer of alphahemolysin increases and causes hemostasis disturbances, thrombocytopenia, and pulmonary lesions during staphylococcal infection. Therefore, it is of interest to inhibit alpha hemolysin using novel compounds. We used the structure of alpha hemolysin(PDB: 7AHL) to screen structures for 100,000 compounds from the ZINC database using molecular docking with AutoDock VINA. Nine (9) successive hits were then subjected for pharmacokinetic and toxicity properties by PROTOX (a webserver for the prediction of oral toxicities of small molecules) and FAFDrugs (a tool for prediction of ADME and Toxicity). This exercise further identified hit #1 ({[3a-(Dihydroxymethyl)-6-hydroxy-2,2-dimethyl-1,3,4-trioxatetrahydro-2H-pentalen-5- yl]methyl}amino(9H-fluoren-9-yl)acetate with binding affinity: -10.3 kcal/mol) and hit #2 (6-(Dihydroxymethyl)-2-{2-[3- (methylamino)propyl]-2-azatricyclo[9.4.0.03,8]pentadeca-1(11),3,5,7,12,14-hexaen-6-yloxy}tetrahydro-2H-pyran-3,4,5-triol with binding affinity: -9.6 kcal/mol) with acceptable toxicity and ADME properties for potential predicted hemolysin inhibition. These compounds should then be evaluated in vitro using inhibitory studies.     

Identification of Potential Carbonic Anhydrase Inhibitors for Glaucoma Treatment Through an In-Silico Approach

Glaucoma is a neurodegenerative disease and second leading cause of blindness in western world. The disease is characterized by an elevated intraocular pressure. Carbonic anhydrase plays a major role by forming aqueous humor and its inhibition can reduce intraocular pressure by partially suppressing the secretion of aqueous humor. Thus in this study, we proposed to identify the potential novel compounds targeting the carbonic anhydrase. The diversity set-II molecules library consisting of 1880 compounds from National Cancer Institute were virtually screened (molecular docking) against human carbonic anhydrase protein. For the obtained best compounds, the nature of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), which determine nucleophilic and electrophilic activity, were calculated by using density functional theory (DFT). The in silico screening suggested 5 best compounds that are effective in comparison to the dorzolamide, a widely used carbonic anhydrase inhibitor for glaucoma treatment. Of the five compounds, 4-nitro-7-[(1-oxidopyridin-1-ium-2-yl) thio] benzofurazan (ZINC01757986) exhibited the better binding affinity (??9.2 cal/mol) in comparison to dorzolamide (??7.2 kcal/mol). The DFT studies on novel identified compound, ZINC01757986 exhibited less HOMO–LUMO energy gap, low hardness and more softness (0.2305 eV, 0.1152 eV and 8.6805 eV) when compared to dorzalamide (0.9536 eV, 0.4768 eV and 2.0973 eV). These studies emphasize that ZINC01757986 can be used as potential carbonic anhydrase inhibitor and lead compounds for the development of an effective anti-glaucoma drug. The results emphasize that these compounds could be potential lead molecules for further structure-based discovery of antiglaucoma drugs.

     

Molecular docking analysis of Clostridium perfringens beta toxin model with potential inhibitors from the ZINC database

Beta toxin from Clostridium perfringens after being secreted in gut is capable of causing necrotic enteritis in humans and several other animal species and does not respond to routinely used antibiotics. Therefore, there is a need to design an effective inhibitor for the Clostridium perfringens beta toxin (CPB) using cutting edge drug discovery technologies. Hence, potential CPB inhibitors were identified using computer aided screening of compounds from the ZINC database. Further, we document the molecular docking analysis of Clostridium perfringens beta toxin model (that revealed 4 binding pockets, A-D) with the identified potential inhibitors. We show that ZINC291192 [N-[(1-methylindol-3-yl) methyl eneamino]-7,10-dioxabicyclo[4.4.0]deca-2,4,11-triene-8- carboxamide] has optimal binding features with calculated binding energy of -10.38 kcal/mol and inhibition constant of 24.76 nM for further consideration.     

Molecular modelling studies unveil potential binding sites on human serum albumin for selected experimental and in silico COVID-19 drug candidate molecules

Human serum albumin (HSA) is the most prevalent protein in the blood plasma which binds an array of exogenous compounds. Drug binding to HSA is an important consideration when developing new therapeutic molecules, and it also aids in understanding the underlying mechanisms that govern their pharmacological effects. This study aims to investigate the molecular binding of coronavirus disease 2019 (COVID-19) therapeutic candidate molecules to HSA and to identify their putative binding sites. Binding energies and interacting residues were used to evaluate the molecular interaction. Four drug candidate molecules (β-D-N4-hydroxycytidine, Chloroquine, Disulfiram, and Carmofur) demonstrate weak binding to HSA, with binding energies ranging from −5 to −6.7 kcal/mol. Ivermectin, Hydroxychloroquine, Remdesivir, Arbidol, and other twenty drug molecules with binding energies ranging from −6.9 to −9.5 kcal/mol demonstrated moderate binding to HSA. The strong HSA binding drug candidates consist of fourteen molecules (Saquinavir, Ritonavir, Dihydroergotamine, Daclatasvir, Paritaprevir etc.) with binding energies ranging from −9.7 to −12.1 kcal/mol. All these molecules bind to different HSA subdomains (IA, IB, IIA, IIB, IIIA, and IIIB) through molecular forces such as hydrogen bonds and hydrophobic interactions. Various pharmacokinetic properties (gastrointestinal absorption, blood-brain barrier permeation, P-glycoprotein substrate, and cytochrome P450 inhibitor) of each molecule were determined using SwissADME program. Further, the stability of the HSA-ligand complexes was analyzed through 100 ns molecular dynamics simulations considering various geometric properties. The binding free energy between free HSA and compounds were calculated using Molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) approach. The findings of this study might be useful in understanding the mechanism of COVID-19 drug candidates binding to serum albumin protein, as well as their pharmacodynamics and pharmacokinetics.Keyword: Human serum albumin, HAS, Serum protein, COVID-19, Molecular docking, Molecular dynamics simulation, Pharmacokinetics, Pharmacodynamics     

Biochemical evaluation and molecular docking assessment of Cymbopogon citratus as a natural source of acetylcholine esterase (AChE)- targeting insecticides

Acetylcholinesterase (AChE) has been an effective target for insecticide development which is a very important aspect of the global fight against insect-borne diseases. The drastic reduction in the sensitivity of insects to AChE-targeting insecticides like organophosphates and carbamates have increased the need for insecticides of natural origin. In this study, we used Drosophila melanogaster as a model to investigate the insecticidal and AChE inhibitory potentials of Cymbopogon citratus and its bioactive compounds. Flies were exposed to 100 and 200 mg/mL C. citratus leaf extract for a 3-h survival assay followed by 45 min exposure for negative geotaxis and biochemical assays. Molecular docking analysis of 45 bioactive compounds of the plant was conducted against Drosophila melanogaster AChE (DmAChE). Exposure to C. citratus significantly reduced the survival rate of flies throughout the exposure period and this was accompanied by a significant decrease in percentage negative geotaxis, AChE activity, catalase activity, total thiol level and a significant increase in glutathione-S-transferase (GST) activity. The bioactive compounds of C. citratus showed varying levels of binding affinities for the enzyme. (+)-Cymbodiacetal scored highest (?9.407 kcal/mol) followed by proximadiol (?8.253 kcal/mol), geranylacetone (?8.177 kcal/mol), and rutin (?8.148 kcal/mol). The four compounds occupied the same binding pocket and interacted with important active site amino acid residues as the co-crystallized ligand (1qon). These compounds could be responsible for the insecticidal and AChE inhibitory potentials of C. citratus and they could be further explored in the development of AChE-targeting insecticides.