Development of Potential Multi-Target Inhibitors for Human Cholinesterases and Beta-Secretase 1: A Computational Approach.

Alzheimer's disease causes chronic neurodegeneration and is the leading cause of dementia in the world. The causes of this disease are not fully understood but seem to involve two essential cerebral pathways: cholinergic and amyloid. The simultaneous inhibition of AChE, BuChE, and BACE-1, essential enzymes involved in those pathways, is a promising therapeutic approach to treat the symptoms and, hopefully, also halt the disease progression. This study sought to identify triple enzymatic inhibitors based on stereo-electronic requirements deduced from molecular modeling of AChE, BuChE, and BACE-1 active sites. A pharmacophore model was built, displaying four hydrophobic centers, three hydrogen bond acceptors, and one positively charged nitrogen, and used to prioritize molecules found in virtual libraries. Compounds showing adequate overlapping rates with the pharmacophore were subjected to molecular docking against the three enzymes and those with an adequate docking score (n = 12) were evaluated for physicochemical and toxicological parameters and commercial availability. The structure exhibiting the greatest inhibitory potential against all three enzymes was subjected to molecular dynamics simulations (100 ns) to assess the stability of the inhibitor-enzyme systems. The results of this in silico approach indicate ZINC1733 can be a potential multi-target inhibitor of AChE, BuChE, and BACE-1, and future enzymatic assays are planned to validate those results.

Identification of potential human beta-secretase 1 inhibitors by hierarchical virtual screening and molecular dynamics.

Alzheimer's disease (AD) is a neurodegenerative pathology responsible for 70% of dementia cases worldwide. Despite its relevance, the few drugs available for the treatment of this disease offer only symptomatic relief, with limited efficacy and serious adverse effects. The most accepted hypothesis about the pathogenesis involves the aggregation and deposition of ß-amyloid peptides, mainly in the cerebral cortex and hippocampus, through the catalytic action of beta-secretase 1 (BACE-1), making this enzyme a promising target for the development of new drugs. In order to prioritize candidates for BACE-1 inhibitors, a hierarchical virtual screening by pharmacophore model and molecular docking was performed against the 216,833 molecules contained in several databases. Our previously built pharmacophore model was used for the first filtering step, which resulted in the selection of 399 molecules. The remaining molecules were filtered through molecular docking with GOLD 5.4.0. In this step, molecules with scoring values ​​greater than the mean plus standard deviation were evaluated for commercial availability and absence of asymmetric centers. Four molecules were selected and evaluated for mutagenic potential by the AMES test with the help of the pkCSM server. Finally, they were tested against the descriptors on Lipinski and Veber rules, and ZINC01589617 (QFIT = 56.52/Score = 44.95) satisfied all requirements, being subjected to molecular dynamics simulations (t = 100 ns) in order to obtain robust data on the mode of bonding and profile of intermolecular interactions. Those in silico strategies demonstrated that ZINC01589617 is a potential candidate for biological tests.Communicated by Ramaswamy H. Sarma.

Hierarchical Virtual Screening and Binding Free Energy Prediction of Potential Modulators of Aedes Aegypti Odorant-Binding Protein 1.

The Aedes aegypti mosquito is the main hematophagous vector responsible for arbovirus transmission in Brazil. The disruption of A. aegypti hematophagy remains one of the most efficient and least toxic methods against these diseases and, therefore, efforts in the research of new chemical entities with repellent activity have advanced due to the elucidation of the functionality of the olfactory receptors and the behavior of mosquitoes. With the growing interest of the pharmaceutical and cosmetic industries in the development of chemical entities with repellent activity, computational studies (e.g., virtual screening and molecular modeling) are a way to prioritize potential modulators with stereoelectronic characteristics (e.g., pharmacophore models) and binding affinity to the AaegOBP1 binding site (e.g., molecular docking) at a lower computational cost. Thus, pharmacophore- and docking-based virtual screening was employed to prioritize compounds from Sigma-Aldrich® (n = 126,851) and biogenic databases (n = 8766). In addition, molecular dynamics (MD) was performed to prioritize the most potential potent compounds compared to DEET according to free binding energy calculations. Two compounds showed adequate stereoelectronic requirements (QFIT > 81.53), AaegOBP1 binding site score (Score > 42.0), volatility and non-toxic properties and better binding free energy value (∆G < −24.13 kcal/mol) compared to DEET ((N,N-diethyl-meta-toluamide)) (∆G = −24.13 kcal/mol).

Mechanism of action involved in the anxiolytic-like effects of Hibalactone isolated from Hydrocotyle umbellata L.

Background and aim: Hibalactone (HB) is a lignan related to the anxiolytic-like effects of Hydrocotyle umbellata L. However, there is a need to understand better the mechanism of action of this lignan to support the ethnopharmacological uses of the species. This work aimed to evaluate by in vivo and in silico analysis the mechanism of action of HB involved in its anxiolytic-like effects. Experimental procedure: The effects of HB in mice were evaluated on light-dark box (LDB) and elevated plus maze (EPM) tests. The participation of 5-HT1A receptor and the benzodiazepine site of GABAA receptor was evaluated to investigate the possible mechanism of action. In silico tools were used to better elucidate the anxiolytic-like effects of HB. Results: Oral treatment with HB at a dose of 33 mg/kg showed an anxiolytic-like effect in the LDB and EPM tests. Besides that, the treatment altered the ethological parameters, frequency of head dips, and stretched-attend postures (SAP), important to better describe the anxiolytic profile of HB. Pretreatment with flumazenil (2 mg/kg) reverted the anxiolytic-like effect of HB on LDB and EPM tests. On the other hand, pretreatment with NAN-190 (0.5 mg/kg) not reverted the activity observed. In silico predictions revealed the potential of HB to increase GABAergic neurotransmission. Pharmacophore modelling and docking simulations showed that HB might interact with the α1ß2γ2 GABAA receptor. Conclusion: Together, the results presented herein suggest that activation of the benzodiazepine site of the GABAA receptor contributes to the anxiolytic-like effect of HB.

A novel antiplasmodial compound: integration of in silico and in vitro assays.

Malaria is a disease caused by Plasmodium genus. which P. falciparum is responsible for the most severe form of the disease, cerebral malaria. In 2018, 405,000 people died of malaria. Antimalarial drugs have serious adverse effects and limited efficacy due to multidrug-resistant strains. One way to overcome these limitations is the use of computational approaches for prioritizing candidates to phenotypic assays and/or in vitro assays against validated targets. Plasmodium falciparum Enoyl-ACP reductase (PfENR) is noteworthy because it catalyzes the rate-limiting step of the biosynthetic pathway of fatty acid. Thus, the study aimed to identify potential PfENR inhibitors by ligand (2D molecular similarity and pharmacophore models) and structure-based virtual screening (molecular docking). 2D similarity-based virtual screening using Tanimoto Index (> 0.45) selected 29,236 molecules from natural products subset available in ZINC database (n = 181,603). Next, 10 pharmacophore models for PfENR inhibitors were generated and evaluated based on the internal statistical parameters from GALAHAD™ and ROC/AUC curve. These parameters selected a suitable pharmacophore model with one hydrophobic center and two hydrogen bond acceptors. The alignment of the filtered molecules on best pharmacophore model resulted in the selection of 10,977 molecules. These molecules were directed to the docking-based virtual screening by AutoDock Vina 1.1.2 program. These strategies selected one compound to phenotypic assays against parasite. ZINC630259 showed EC50 = 0.12 ± 0.018 µM in antiplasmodial assays and selective index similar to other antimalarial drugs. Finally, MM/PBSA method showed stability of molecule within PfENR binding site (ΔGbinding=-57.337 kJ/mol).Communicated by Ramaswamy H. Sarma.

Identification of Bioactive Compounds against Aedes aegypti (Diptera: Culicidae) by Bioassays and in Silico Assays.

Most of the hematophagous insects act as disease vectors, including Aedes aegypti, responsible for transmitting some of the most critical arboviruses globally, such as Dengue. The use of repellents based on natural products is a promising alternative for personal protection compared to industrial chemical repellents. In this study, the repellent effect of essential oils extracted from Lippia thymoides, Lippia alba, Cymbopogon winterianus, and Eucalyptus globulus leaves was evaluated. Essential oils used showed repellent activity against Ae. aegypti in laboratory bioassays, obtaining protection rates above 70 % from 3.75 mg/mL and higher concentration for all analyzed oils. GC/MS identified 57 constituents, which were used in the ligand-based pharmacophore model to expose compounds with requirements for repellents that modulate mosquitoes behavior through odorant-binding protein 1 Ae. aegypti. Ligand-based pharmacophore model approach results suggested that repellent activity from C. winterianus, L. alba, and L. thymoides essential oils' metabolites is related to Citronelal (QFIT=26.77), Citronelol (QFIT=11.29), Citronelol acetate (QFIT=52.22) and Geranil acetate (QFIT=10.28) with synergistic or individual activity. E. globulus essential oil's repellent activity is associated with Ledol (0.94 %; QFIT=41.95). Molecular docking was applied to understand the binding mode and affinity of the essential oils' data set at the protein binding site. According to molecular docking, Citronelol (ChemPLP=60.98) and geranyl acetate (ChemPLP=60.55) were the best-classified compounds compared to the others and they can be explored to develop new repellents.

Repositioning of Etravirine as a Potential CK1ε Inhibitor by Virtual Screening.

CK1ε is a key regulator of WNT/ß-catenin and other pathways that are linked to tumor progression; thus, CK1ε is considered a target for the development of antineoplastic therapies. In this study, we performed a virtual screening to search for potential CK1ε inhibitors. First, we characterized the dynamic noncovalent interactions profiles for a set of reported CK1ε inhibitors to generate a pharmacophore model, which was used to identify new potential inhibitors among FDA-approved drugs. We found that etravirine and abacavir, two drugs that are approved for HIV infections, can be repurposed as CK1ε inhibitors. The interaction of these drugs with CK1ε was further examined by molecular docking and molecular dynamics. Etravirine and abacavir formed stable complexes with the target, emulating the binding behavior of known inhibitors. However, only etravirine showed high theoretical binding affinity to CK1ε. Our findings provide a new pharmacophore for targeting CK1ε and implicate etravirine as a CK1ε inhibitor and antineoplastic agent.

Identification of novel antiplasmodial compound by hierarquical virtual screening and in vitro assays.

Malaria is an infectious disease caused by protozoa of the genus Plasmodium spp. with approximately 219 million cases in 2017. P. falciparum is main responsible for the most severe form of the disease, cerebral malaria. Despite of public health impacts, chemotherapy against malaria is still limited due to the emergence of drug resistance cases used in monotherapy and combination therapies. Thus, the development of new antimalarial drugs becomes emergency. One way of achieve this goal is to explore essential and/or unique therapeutic targets of the parasite, or at least sufficiently different to ensure selective inhibition. Enoil-ACP reductase (ENR) is a NADH-dependent enzyme responsible for the limiting step of the type II fatty acid biosynthetic pathway (FAS II). Thus, pharmacophore and docking based virtual screening were applied to prioritize molecules for in vitro assays against P. falciparum W2 strain. The application of successive filters at OOCC database (n = 618) resulted in the identification of one molecule (13) (EC50 = 0.098 ± 0.021 µM) with similar biological activity to artemether. The molecule 13 is a typical drug repurposing case due to previous other approved therapeutic uses on Chinese medicine as a non-specific cholinergic antagonist, thus it could be accelerated the drug development process. Additionally, molecular dynamics studies were used to confirm stability of the molecular interactions identified by molecular docking. Thus, representative structures of P. falciparum ENR can be used in a study to propose new derivatives for evaluation of biological activity in vitro and in vivo. Communicated by Ramaswamy H. Sarma.

Synthesis of Pyrrolo[3,4-b]pyridin-5-ones via Multicomponent Reactions and In Vitro-In Silico Studies Against SiHa, HeLa, and CaSki Human Cervical Carcinoma Cell Lines.

A series of 12 polysubstituted pyrrolo[3,4-b]pyridin-5-ones were synthesized via a one-pot cascade process (Ugi-3CR/aza Diels-Alder/N-acylation/decarboxylation/dehydration) and studied in vitro using human epithelial cervical carcinoma SiHa, HeLa, and CaSki cell line cultures. Three compounds of the series exhibited significative cytotoxicity against the three cell lines, with HeLa being the most sensitive one. Then, based on these results, in silico studies by docking techniques were performed using Paclitaxel as a reference and αß-tubulin as the selected biological target. Worth highlighting is that strong hydrophobic interactions were observed between the three active molecules and the reference drug Paclitaxel, to the αß-tubulin. In consequence, it was determined that hydrophobic-aromatic moieties of bioactive compounds and Paclitaxel play a key role in making stronger interactions to the ligand-target complex. A quantitative structure activity relationship (QSAR) study revealed that the six membered rings are the most significant molecular frameworks, being present in all proposed models for the in vitro-studied cell lines. Finally, also from the docking interpretation, a ligand-based pharmacophore model is proposed in order to find further potential polyheterocyclic candidates to bind stronger to the αß-tubulin.

In Silico Study to Identify New Antituberculosis Molecules from Natural Sources by Hierarchical Virtual Screening and Molecular Dynamics Simulations.

Tuberculosis (TB) is an infection caused by Mycobacterium tuberculosis, responsible for 1.5 million documented deaths in 2016. The increase in reported cases of M. tuberculosis resistance to the main drugs show the need for the development of new and efficient drugs for better TB control. Based on these facts, this work aimed to use combined in silico techniques for the discovery of potential inhibitors to ß-ketoacyl-ACP synthase (MtKasA). Initially compounds from natural sources present in the ZINC database were selected, then filters were sequentially applied by virtual screening, initially with pharmacophoric modeling, and later the selected compounds (based on QFIT scores) were submitted to the DOCK 6.5 program. After recategorization of the variables (QFIT score and GRID score), compounds ZINC35465970 and ZINC31170017 were selected. These compounds showed great hydrophobic contributions and for each established system 100 ns of molecular dynamics simulations were performed and the binding free energy was calculated. ZINC35465970 demonstrated a greater capacity for the KasA enzyme inhibition, with a ΔGbind = -30.90 kcal/mol and ZINC31170017 presented a ΔGbind = -27.49 kcal/mol. These data can be used in other studies that aim at the inhibition of the same biological targets through drugs with a dual action.

Virtual screening to identify Leishmania braziliensis N-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.

Identification of new promising Plasmodium falciparum superoxide dismutase allosteric inhibitors through hierarchical pharmacophore-based virtual screening and molecular dynamics.

Malaria is the world's most widespread protozoan infection, being responsible for more than 445,000 annual deaths. Among the malaria parasites, Plasmodium falciparum is the most prevalent and lethal. In this context, the search for new antimalarial drugs is urgently needed. P. falciparum superoxide dismutase (PfSOD) is an important enzyme involved in the defense mechanism against oxidative stress. The goal of this study was to identify through hierarchical screening on pharmacophore models and molecular dynamics (MD), promising allosteric PfSOD inhibitors that do not show structural requirements for human inhibition. MD simulations of 1000 ps were performed on PfSOD using GROMACS 5.1.2. For this, the AMBER99SB-ILDN force field was adapted to describe the metal-containing system. The simulations indicated stability in the developed system. Therefore, a covariance matrix was generated, in which it was possible to identify residues with correlated and anticorrelated movements with the active site. These results were associated with the results found in the predictor of allosteric sites, AlloSitePro, which affirmed the ability of these residues to delimit an allosteric site. Then, after successive filtering of the Sigma-Aldrich® compounds database for HsSOD1 and PfSOD pharmacophores, 152 compounds were selected, also obeying Lipinski's rule of 5. Further filtering of those compounds based on molecular docking results, toxicity essays, availability, and price filtering led to the selection of a best compound, which was then submitted to MD simulations of 20,000 ps on the allosteric site. The study concludes that the ZINC00626080 compound could be assayed against SODs. Graphical Abstract Plasmodium falciparum superoxide dismutase.

Nanomolar anti-sickling compounds identified by ligand-based pharmacophore approach.

Adenosine receptors are considered as potential targets for drug development against several diseases. The discovery of subtype 2B adenosine receptors role in erythrocyte sickling process proved its importance to neglected diseases such as sickle cell anemia, which affects approximately 29.000 people around the world, but whose treatment is restricted to just one FDA approved drug (hydroxyurea). In order to widen the therapeutic arsenal available to treat sickle cell anemia patients, it is imperative to identify new lead compounds that modify the sickling course and not just its symptoms. In order to accomplish this goal, ligand-based pharmacophore models that differentiate true ligands from decoys and enlighten the structure-activity relationship of known RA2B antagonists were employed screen the lead-like subset of the ZINC database. Following a chemical diversity analysis, 18 compounds were selected for biological evaluation. Among them, one molecule Z1139491704 (pEC50 = 7.77 ± 0.17) has shown better anti-sickling activity than MRS1754 (pEC50 = 7.63 ± 0.12), a commercial RA2B antagonist. Moreover, these compounds exhibited no cytotoxic effect at low micromolar range on mammalian cells. In conclusion, the sound development of validated ligand-based pharmacophore models proved essential to identify novel chemical scaffolds that might be useful to develop anti-sickling drugs.