In Silico and Experimental ADAM17 Kinetic Modeling as Basis for Future Screening System for Modulators.

Understanding the mechanisms of modulators' action on enzymes is crucial for optimizing and designing pharmaceutical substances. The acute inflammatory response, in particular, is regulated mainly by a disintegrin and metalloproteinase (ADAM) 17. ADAM17 processes several disease mediators such as TNFα and APP, releasing their soluble ectodomains (shedding). A malfunction of this process leads to a disturbed inflammatory response. Chemical protease inhibitors such as TAPI-1 were used in the past to inhibit ADAM17 proteolytic activity. However, due to ADAM17's broad expression and activity profile, the development of active-site-directed ADAM17 inhibitor was discontinued. New 'exosite' (secondary substrate binding site) inhibitors with substrate selectivity raised the hope of a substrate-selective modulation as a promising approach for inflammatory disease therapy. This work aimed to develop a high-throughput screen for potential ADAM17 modulators as therapeutic drugs. By combining experimental and in silico methods (structural modeling and docking), we modeled the kinetics of ADAM17 inhibitor. The results explain ADAM17 inhibition mechanisms and give a methodology for studying selective inhibition towards the design of pharmaceutical substances with higher selectivity.

Design of SARS-CoV-2 PLpro Inhibitors for COVID-19 Antiviral Therapy Leveraging Binding Cooperativity.

Antiviral agents that complement vaccination are urgently needed to end the COVID-19 pandemic. The SARS-CoV-2 papain-like protease (PLpro), one of only two essential cysteine proteases that regulate viral replication, also dysregulates host immune sensing by binding and deubiquitination of host protein substrates. PLpro is a promising therapeutic target, albeit challenging owing to featureless P1 and P2 sites recognizing glycine. To overcome this challenge, we leveraged the cooperativity of multiple shallow binding sites on the PLpro surface, yielding novel 2-phenylthiophenes with nanomolar inhibitory potency. New cocrystal structures confirmed that ligand binding induces new interactions with PLpro: by closing of the BL2 loop of PLpro forming a novel "BL2 groove" and by mimicking the binding interaction of ubiquitin with Glu167 of PLpro. Together, this binding cooperativity translates to the most potent PLpro inhibitors reported to date, with slow off-rates, improved binding affinities, and low micromolar antiviral potency in SARS-CoV-2-infected human cells.

Discovery of novel antagonists targeting the DNA binding domain of androgen receptor by integrated docking-based virtual screening and bioassays.

Androgen receptor (AR), a ligand-activated transcription factor, is a master regulator in the development and progress of prostate cancer (PCa). A major challenge for the clinically used AR antagonists is the rapid emergence of resistance induced by the mutations at AR ligand binding domain (LBD), and therefore the discovery of novel anti-AR therapeutics that can combat mutation-induced resistance is quite demanding. Therein, blocking the interaction between AR and DNA represents an innovative strategy. However, the hits confirmed targeting on it so far are all structurally based on a sole chemical scaffold. In this study, an integrated docking-based virtual screening (VS) strategy based on the crystal structure of the DNA binding domain (DBD) of AR was conducted to search for novel AR antagonists with new scaffolds and 2-(2-butyl-1,3-dioxoisoindoline-5-carboxamido)-4,5-dimethoxybenzoicacid (Cpd39) was identified as a potential hit, which was competent to block the binding of AR DBD to DNA and showed decent potency against AR transcriptional activity. Furthermore, Cpd39 was safe and capable of effectively inhibiting the proliferation of PCa cell lines (i.e., LNCaP, PC3, DU145, and 22RV1) and reducing the expression of the genes regulated by not only the full-length AR but also the splice variant AR-V7. The novel AR DBD-ARE blocker Cpd39 could serve as a starting point for the development of new therapeutics for castration-resistant PCa.

Advances in Magnetic Microbead Affinity Selection Screening: Discovery of Natural Ligands to the SARS-CoV-2 Spike Protein.

Affinity selection-mass spectrometry, which includes magnetic microbead affinity selection-screening (MagMASS), is ideal for the discovery of ligands in complex mixtures that bind to pharmacological targets. Therapeutic agents are needed to prevent or treat COVID-19, which is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Infection of human cells by SARS-CoV-2 involves binding of the virus spike protein subunit 1 (S1) to the human cell receptor angiotensin converting enzyme-2 (ACE2). Like antibodies, small molecules have the potential to block the interaction of the viral S1 protein with human ACE2 and prevent SARS-CoV-2 infection. Therefore, a MagMASS assay was developed for the discovery of ligands to the S1 protein. Unlike previous MagMASS approaches, this new assay used robotics for 5-fold enhancement of throughput and sensitivity. The assay was validated using the SBP-1 peptide, which is identical to the ACE2 amino acid sequence recognized by the S1 protein, and then applied to the discovery of natural ligands from botanical extracts. Small molecule ligands to the S1 protein were discovered in extracts of the licorice species, Glycyrrhiza inflata. In particular, the licorice ligand licochalcone A was identified through dereplication and comparison with standards using HPLC with high-resolution tandem mass spectrometry.

Targeting novel coronavirus SARS-CoV-2 spike protein with phytoconstituents of Momordica charantia.

BACKGROUND: Infections by the SARS-CoV-2 virus causing COVID-19 are presently a global emergency. The current vaccination effort may reduce the infection rate, but strain variants are emerging under selection pressure. Thus, there is an urgent need to find drugs that treat COVID-19 and save human lives. Hence, in this study, we identified phytoconstituents of an edible vegetable, Bitter melon (Momordica charantia), that affect the SARS-CoV-2 spike protein. METHODS: Components of Momordica charantia were tested to identify the compounds that bind to the SARS-CoV-2 spike protein. An MTiOpenScreen web-server was used to perform docking studies. The Lipinski rule was utilized to evaluate potential interactions between the drug and other target molecules. PyMol and Schrodinger software were used to identify the hydrophilic and hydrophobic interactions. Surface plasmon resonance (SPR) was employed to assess the interaction between an extract component (erythrodiol) and the spike protein. RESULTS: Our in-silico evaluations showed that phytoconstituents of Momordica charantia have a low binding energy range, -5.82 to -5.97 kcal/mol. A docking study revealed two sets of phytoconstituents that bind at the S1 and S2 domains of SARS-CoV-2. SPR showed that erythrodiol has a strong binding affinity (KD = 1.15 µM) with the S2 spike protein of SARS-CoV-2. Overall, docking, ADME properties, and SPR displayed strong interactions between phytoconstituents and the active site of the SARS-CoV-2 spike protein. CONCLUSION: This study reveals that phytoconstituents from bitter melon are potential agents to treat SARS-CoV-2 viral infections due to their binding to spike proteins S1 and S2.

Flavonoids in Ampelopsis grossedentata as covalent inhibitors of SARS-CoV-2 3CLpro: Inhibition potentials, covalent binding sites and inhibitory mechanisms.

Coronavirus 3C-like protease (3CLpro) is a crucial target for treating coronavirus diseases including COVID-19. Our preliminary screening showed that Ampelopsis grossedentata extract (AGE) displayed potent SARS-CoV-2-3CLpro inhibitory activity, but the key constituents with SARS-CoV-2-3CLpro inhibitory effect and their mechanisms were unrevealed. Herein, a practical strategy via integrating bioactivity-guided fractionation and purification, mass spectrometry-based peptide profiling and time-dependent biochemical assay, was applied to identify the crucial constituents in AGE and to uncover their inhibitory mechanisms. The results demonstrated that the flavonoid-rich fractions (10-17.5 min) displayed strong SARS-CoV-2-3CLpro inhibitory activities, while the constituents in these fractions were isolated and their SARS-CoV-2-3CLpro inhibitory activities were investigated. Among all isolated flavonoids, dihydromyricetin, isodihydromyricetin and myricetin strongly inhibited SARS-CoV-2 3CLpro in a time-dependent manner. Further investigations demonstrated that myricetin could covalently bind on SARS-CoV-2 3CLpro at Cys300 and Cys44, while dihydromyricetin and isodihydromyricetin covalently bound at Cys300. Covalent docking coupling with molecular dynamics simulations showed the detailed interactions between the orthoquinone form of myricetin and two covalent binding sites (surrounding Cys300 and Cys44) of SARS-CoV-2 3CLpro. Collectively, the flavonoids in AGE strongly and time-dependently inhibit SARS-CoV-2 3CLpro, while the newly identified SARS-CoV-2 3CLpro inhibitors in AGE offer promising lead compounds for developing novel antiviral agents.

Novel STAT3 small-molecule inhibitors identified by structure-based virtual ligand screening incorporating SH2 domain flexibility.

Efforts to develop STAT3 inhibitors have focused on its SH2 domain starting with short phosphotyrosylated peptides based on STAT3 binding motifs, e.g. pY905LPQTV within gp130. Despite binding to STAT3 with high affinity, issues regarding stability, bioavailability, and membrane permeability of these peptides, as well as peptidomimetics such as CJ-887, have limited their further clinical development and led to interest in small-molecule inhibitors. Some small molecule STAT3 inhibitors, identified using structure-based virtual ligand screening (SB-VLS); while having favorable drug-like properties, suffer from weak binding affinities, possibly due to the high flexibility of the target domain. We conducted molecular dynamic (MD) simulations of the SH2 domain in complex with CJ-887, and used an averaged structure from this MD trajectory as an "induced-active site" receptor model for SB-VLS of 110,000 compounds within the SPEC database. Screening was followed by re-docking and re-scoring of the top 30% of hits, selection for hit compounds that directly interact with pY + 0 binding pocket residues R609 and S613, and testing for STAT3 targeting in vitro, which identified two lead hits with good activity and favorable drug-like properties. Unlike most small-molecule STAT3 inhibitors previously identified, which contain negatively-charged moieties that mediate binding to the pY + 0 binding pocket, these compounds are uncharged and likely will serve as better candidates for anti-STAT3 drug development. IMPLICATIONS: SB-VLS, using an averaged structure from molecular dynamics (MD) simulations of STAT3 SH2 domain in a complex with CJ-887, a known peptidomimetic binder, identify two highly potent, neutral, low-molecular weight STAT3-inhibitors with favorable drug-like properties.

Identification of novel anti-inflammatory Nur77 modulators by virtual screening.

Orphan nuclear receptor Nur77 is a unique member of the NR4A nuclear receptor subfamily, which is critical for cellular processes especially the inflammatory responses. Many efforts have been made to discover novel scaffold small molecules targeting Nur77. Herein, we evaluated the previously reported binding sites in crystal structures of Nur77 with small molecules, and then discovered compound 13 as a hit of Nur77 via virtual screening targeting the best-scored binding site. Based on the results of fluorescence titration assay, structure-activity relationship (SAR) analysis was summarized for compound 13 and its analogs. Among these analogs, compound 13e displayed the most potent binding affinity (0.54 ± 0.02 µM). The binding mode of compound 13e was predicted via molecule docking. Moreover, 13e exhibited significant anti-inflammation activity in TNF-α induced HepG2 cell model. Taken together, these results provided a new insight into the understanding the functions of specific binding sites on Nur77 for small molecular compounds, and the development of new scaffold Nur77 modulators.

Protein druggability assessment for natural products using in silico simulation: A case study with estrogen receptor and the flavonoid genistein.

Traditional herbal medicine (THM) comprises a vast number of natural compounds. Most of them are metabolized into different structures after administration, which makes the clarification of THM's mode of action more complicated. To evaluate the biological activities of those components and metabolites, in silico simulation technology is helpful. We focused on mixed-solvent molecular dynamics (MD) simulation for druggability assessment of natural products. Mixed-solvent MD is an in silico simulation method for the exploration of ligand-binding sites on target proteins, which uses water and an organic molecule mixture. The selection of organic small molecules is an important factor for predicting the characteristics of natural products. In this study, we used the known crystal structure of estrogen receptors with genistein as a test case and explored fragments reflecting the characteristics of natural products. We found that structures with a 4-pyrone structure are more often included in the natural products database compared with the DrugBank database, and we selectively detected the known-binding sites of estrogen receptor α and ß. The results indicate that the 4-pyrone structure might be promising for predicting the protein druggability of flavonoids. Additionally, mixed-solvent MD simulation discriminates the selectivity of genistein between estrogen receptor ß and α, indicating that the simulation can be evaluated using indices that differ from those of traditional ligand docking. Although this approach is still in its early stages, it has the potential to provide valuable information for understanding the diverse biological activities of natural products.

Bioassay-based Corchorus capsularis L. leaf-derived ß-sitosterol exerts antileishmanial effects against Leishmania donovani by targeting trypanothione reductase.

Leishmaniasis, a major neglected tropical disease, affects millions of individuals worldwide. Among the various clinical forms, visceral leishmaniasis (VL) is the deadliest. Current antileishmanial drugs exhibit toxicity- and resistance-related issues. Therefore, advanced chemotherapeutic alternatives are in demand, and currently, plant sources are considered preferable choices. Our previous report has shown that the chloroform extract of Corchorus capsularis L. leaves exhibits a significant effect against Leishmania donovani promastigotes. In the current study, bioassay-guided fractionation results for Corchorus capsularis L. leaf-derived ß-sitosterol (ß-sitosterolCCL) were observed by spectroscopic analysis (FTIR, 1H NMR, 13C NMR and GC-MS). The inhibitory efficacy of this ß-sitosterolCCL against L. donovani promastigotes was measured (IC50 = 17.7 ± 0.43 µg/ml). ß-SitosterolCCL significantly disrupts the redox balance via intracellular ROS production, which triggers various apoptotic events, such as structural alteration, increased storage of lipid bodies, mitochondrial membrane depolarization, externalization of phosphatidylserine and non-protein thiol depletion, in promastigotes. Additionally, the antileishmanial activity of ß-sitosterolCCL was validated by enzyme inhibition and an in silico study in which ß-sitosterolCCL was found to inhibit Leishmania donovani trypanothione reductase (LdTryR). Overall, ß-sitosterolCCL appears to be a novel inhibitor of LdTryR and might represent a successful approach for treatment of VL in the future.

Discovery of Potential Chemical Probe as Inhibitors of CXCL12 Using Ligand-Based Virtual Screening and Molecular Dynamic Simulation.

CXCL12 are small pro-inflammatory chemo-attractant cytokines that bind to a specific receptor CXCR4 with a role in angiogenesis, tumor progression, metastasis, and cell survival. Globally, cancer metastasis is a major cause of morbidity and mortality. In this study, we targeted CXCL12 rather than the chemokine receptor (CXCR4) because most of the drugs failed in clinical trials due to unmanageable toxicities. Until now, no FDA approved medication has been available against CXCL12. Therefore, we aimed to find new inhibitors for CXCL12 through virtual screening followed by molecular dynamics simulation. For virtual screening, active compounds against CXCL12 were taken as potent inhibitors and utilized in the generation of a pharmacophore model, followed by validation against different datasets. Ligand based virtual screening was performed on the ChEMBL and in-house databases, which resulted in successive elimination through the steps of pharmacophore-based and score-based screenings, and finally, sixteen compounds of various interactions with significant crucial amino acid residues were selected as virtual hits. Furthermore, the binding mode of these compounds were refined through molecular dynamic simulations. Moreover, the stability of protein complexes, Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), and radius of gyration were analyzed, which led to the identification of three potent inhibitors of CXCL12 that may be pursued in the drug discovery process against cancer metastasis.

Novel and potent inhibitors targeting DHODH are broad-spectrum antivirals against RNA viruses including newly-emerged coronavirus SARS-CoV-2.

Emerging and re-emerging RNA viruses occasionally cause epidemics and pandemics worldwide, such as the on-going outbreak of the novel coronavirus SARS-CoV-2. Herein, we identified two potent inhibitors of human DHODH, S312 and S416, with favorable drug-likeness and pharmacokinetic profiles, which all showed broad-spectrum antiviral effects against various RNA viruses, including influenza A virus, Zika virus, Ebola virus, and particularly against SARS-CoV-2. Notably, S416 is reported to be the most potent inhibitor so far with an EC50 of 17 nmol/L and an SI value of 10,505.88 in infected cells. Our results are the first to validate that DHODH is an attractive host target through high antiviral efficacy in vivo and low virus replication in DHODH knock-out cells. This work demonstrates that both S312/S416 and old drugs (Leflunomide/Teriflunomide) with dual actions of antiviral and immuno-regulation may have clinical potentials to cure SARS-CoV-2 or other RNA viruses circulating worldwide, no matter such viruses are mutated or not.

Virtual screening based on molecular docking of possible inhibitors of Covid-19 main protease.

Coronavirus (COVID-19) is an enveloped RNA virus that is diversely found in humans and that has now been declared a global pandemic by the World Health Organization. Thus, there is an urgent need to develop effective therapies and vaccines against this disease. In this context, this study aimed to evaluate in silico the molecular interactions of drugs with therapeutic indications for treatment of COVID-19 (Azithromycin, Baricitinib and Hydroxychloroquine) and drugs with similar structures (Chloroquine, Quinacrine and Ruxolitinib) in docking models from the SARS-CoV-2 main protease (M-pro) protein. The results showed that all inhibitors bound to the same enzyme site, more specifically in domain III of the SARS-CoV-2 main protease. Therefore, this study allows proposing the use of baricitinib and quinacrine, in combination with azithromycin; however, these computer simulations are just an initial step for conceiving new projects for the development of antiviral molecules.

The selenium-containing drug ebselen potently disrupts LEDGF/p75-HIV-1 integrase interaction by targeting LEDGF/p75.

Lens-epithelium-derived growth-factor (LEDGF/p75)-binding site on HIV-1 integrase (IN), is an attractive target for antiviral chemotherapy. Small-molecule compounds binding to this site are referred as LEDGF-IN inhibitors (LEDGINs). In this study, compound libraries were screened to identify new inhibitors of LEDGF/p75-IN interaction. Ebselen (2-phenyl-1,2-benzisoselenazol-3-one), a reported anti-HIV-1 agent, was identified as a moderate micromolar inhibitor of LEDGF/p75-IN interaction. Ebselen inhibited the interaction by binding to LEDGF/p75 and the ability of ebselen to inhibit the interaction could be reversed by dithiothreitol (DTT). BLI experiment showed that ebselen probably formed selenium-sulphur bonds with reduced thiols in LEDGF/p75. To the best of our knowledge, we showed for the first time that small-molecule compound, ebselen inhibited LEDGF/p75-IN interaction by directly binding to LEDGF/p75. The compound discovered here could be used as probe compounds to design and develop new disrupter of LEDGF/p75-IN interaction.

In Silico Inhibition of BACE-1 by Selective Phytochemicals as Novel Potential Inhibitors: Molecular Docking and DFT Studies.

BACKGROUND: Alzheimer's Disease (AD) has become the most common age-dependent disease of dementia. The trademark pathologies of AD are the presence of amyloid aggregates in neurofibrils. Recently phytochemicals being considered as potential inhibitors against various neurodegenerative, antifungal, antibacterial and antiviral diseases in human beings. OBJECTIVE: This study targets the inhibition of BACE-1 by phytochemicals using in silico drug discovery analysis. METHODS: A total of 3150 phytochemicals were collected from almost 25 different plants through literature assessment. The ADMET studies, molecular docking and density functional theory (DFT) based analysis were performed to analyze the potential inhibitory properties of these phytochemicals. RESULTS: The ADMET and docking results exposed seven compounds that have high potential as an inhibitory agent against BACE-1 and show binding affinity >8.0 kcal/mol against BACE-1. They show binding affinity greater than those of various previously reported inhibitors of BACE-1. Furthermore, DFT based analysis has shown high reactivity for these seven phytochemicals in the binding pocket of BACE- 1, based on ELUMO, EHOMO and Kohn-Sham energy gap. All seven phytochemicals were testified (as compared to experimental ones) as novel inhibitors against BACE-1. CONCLUSION: Out of seven phytochemicals, four were obtained from plant Glycyrrhiza glabra i.e. Shinflavanone, Glabrolide, Glabrol and PrenyllicoflavoneA, one from Huperzia serrate i.e. Macleanine, one from Uncaria rhynchophylla i.e. 3a-dihydro-cadambine and another one was from VolvalerelactoneB from plant Valeriana-officinalis. It is concluded that these phytochemicals are suitable candidates for drug/inhibitor against BACE-1, and can be administered to humans after experimental validation through in vitro and in vivo trials.

Investigation of Binding Characteristics of Phosphoinositide-dependent Kinase-1 (PDK1) Co-crystallized Ligands Through Virtual Pharmacophore Modeling Leading to Novel Anti-PDK1 Hits.

BACKGROUND: 3-Phosphoinositide Dependent Protein Kinase-1 (PDK1) is being lately considered as an attractive and forthcoming anticancer target. A Protein Data Bank (PDB) cocrystallized crystal provides not only rigid theoretical data but also a realistic molecular recognition data that can be explored and used to discover new hits. OBJECTIVE: This incited us to investigate the co-crystallized ligands' contacts inside the PDK1 binding pocket via a structure-based receptor-ligand pharmacophore generation technique in Discovery Studio 4.5 (DS 4.5). METHODS: Accordingly, 35 crystals for PDK1 were collected and studied. Every single receptorligand interaction was validated and the significant ones were converted into their corresponding pharmacophoric features. The generated pharmacophores were scored by the Receiver Operating Characteristic (ROC) curve analysis. RESULTS: Consequently, 169 pharmacophores were generated and sorted, 11 pharmacophores acquired good ROC-AUC results of 0.8 and a selectivity value above 8. Pharmacophore 1UU3_2_01 was used in particular as a searching filter to screen NCI database because of its acceptable validity criteria and its distinctive positive ionizable feature. Several low micromolar PDK1 inhibitors were revealed. The most potent hit illustrated anti-PDK1 IC50 values of 200 nM with 70% inhibition against SW480 cell lines. CONCLUSION: Eventually, the active hits were docked inside the PDK1 binding pocket and the recognition points between the active hits and the receptor were analyzed that led to the discovery of new scaffolds as potential PDK1 inhibitors.

Cilantro leaf harbors a potent potassium channel-activating anticonvulsant.

Herbs have a long history of use as folk medicine anticonvulsants, yet the underlying mechanisms often remain unknown. Neuronal voltage-gated potassium channel subfamily Q (KCNQ) dysfunction can cause severe epileptic encephalopathies that are resistant to modern anticonvulsants. Here we report that cilantro (Coriandrum sativum), a widely used culinary herb that also exhibits antiepileptic and other therapeutic activities, is a highly potent KCNQ channel activator. Screening of cilantro leaf metabolites revealed that one, the long-chain fatty aldehyde (E)-2-dodecenal, activates multiple KCNQs, including the predominant neuronal isoform, KCNQ2/KCNQ3 [half maximal effective concentration (EC50), 60 ± 20 nM], and the predominant cardiac isoform, KCNQ1 in complexes with the type I transmembrane ancillary subunit (KCNE1) (EC50, 260 ± 100 nM). (E)-2-dodecenal also recapitulated the anticonvulsant action of cilantro, delaying pentylene tetrazole-induced seizures. In silico docking and mutagenesis studies identified the (E)-2-dodecenal binding site, juxtaposed between residues on the KCNQ S5 transmembrane segment and S4-5 linker. The results provide a molecular basis for the therapeutic actions of cilantro and indicate that this ubiquitous culinary herb is surprisingly influential upon clinically important KCNQ channels.-Manville, R. W., Abbott, G. W. Cilantro leaf harbors a potent potassium channel-activating anticonvulsant.

Identification of New Potent Human Uncoupling Protein 1 (UCP1) Agonists Using Virtual Screening and in vitro Approaches.

Recent studies suggested that activation of Uncoupling Protein 1 (UCP1) has become an appealing therapeutic strategy against obesity and diabetes. In our research, the 3D structure of UCP1 was constructed through homology modelling, refined through molecular dynamics simulation, and evaluated by Ramachandran plot, the molecular docking of UCP1 activators brought about the proposal of an interaction mode inside the UCP1 active site. Remarkably, Reside Lys126 formed hydrogen bond; residues Pro121, Val125, Tyr146, Tyr149 and Arg150 formed hydrophobic interaction, which are key amino acids within UCP1 site. Then a pharmacophore model was generated consisting of three hydrophobic groups, a negative center and an additional hydrophobic group. Pharmacophore-based virutal screening of Specs database yield 5 hits. In vitro assay indicated ZINC 04660290 significantly increased the protein expression of UCP1 and decreased the fat droplet in a dose-dependent manner. Besides, pharmacokinetic properties were predicted for those five compounds through ADME/T prediction. All of these will guide us to design new UCP1 activators for the treatment of obesity and diabetes.

Prediction of Ligands Binding Acetylcholinesterase with Potential Antidotal Activity: A Virtual Screening Approach.

Intoxications caused by organophosphorus compounds (OPs) are associated with the reversible, and sometimes irreversible interaction with acetylcholinesterase (AChE). OPs are commonly used as pesticides mainly in developing countries, where the associated poisoning is a major health problem related to suicidal attempts, careless manipulation, and chemical warfare. The current antidotes are oxime-based drugs that can regenerate the AChE catalytic activity. Nevertheless, challenges associated with lack of efficiency and difficulties for crossing blood-brain barrier have motivated the design of novel alternatives. We used a validated molecular docking approach for the virtual screening of 579,890 synthetic ligands and 478 drugs against a human AChE in its apo conformation, and a murine AChE conjugated with the OP tabun. After filtering, 7 hits were selected as potential competitors due to the formation of key interactions within the active site gorge of the AChE structure, and potential reactivators based on interactions with amino acids of the catalytic triad in the presence of organophosphorus compounds. The selected candidates will be further evaluated through in vitro and in vivo assays.