Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has received global attention due to the serious threat it poses to public health. Since the outbreak in December 2019, millions of people have been affected and its rapid global spread has led to an upsurge in the search for treatment. To discover hit compounds that can be used alone or in combination with repositioned drugs, we first analyzed the pharmacokinetic and toxicological properties of natural products from Brazil’s semiarid region. After, we analyzed the site prediction and druggability of the SARS-CoV-2 main protease (Mpro), followed by docking and molecular dynamics simulation. The best SARS-CoV-2 Mpro complexes revealed that other sites were accessed, confirming that our approach could be employed as a suitable starting protocol for ligand prioritization, reinforcing the importance of catalytic cysteine-histidine residues and providing new structural data that could increase the antiviral development mainly against SARSCoV-2. Here, we selected 10 molecules that could be in vitro assayed in response to COVID-19. Two compounds (b01 and b02) suggest a better potential for interaction with SARS-CoV-2 Mpro and could be further studied.Research Dean and Graduate Studies of the Federal University of Pará (PROPESP/UFPA)Brazilian National Council for Scientific and Technological Development (CNPq)Brazilian Coordination for Improvement of Personnel Higher Education (CAPES)Bahia Research Foundation (FAPESB, grant numbers APP071/2011, JCB-0039/2013, and RED-008/2013

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.PPBE and PPGCF/UEFS; Fundação de Amparo à Pesquisa do Estado de Minas Gerais—FAPEMIG, grants APQ-02741-17, APQ-00855-19, APQ-01733-21, and APQ-04559-22Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq-Brazil, grants 305117/2017-3, 426261/2018-6Fellowship of 2021 (grant 310108/2020-9

Identification of New Rofecoxib-Based Cyclooxygenase-2 Inhibitors: A Bioinformatics Approach

The cyclooxygenase-2 receptor is a therapeutic target for planning potential drugs with anti-inflammatory activity. The selective cyclooxygenase-2 (COX-2) inhibitor rofecoxib was selected as a pivot molecule to perform virtual ligand-based screening from six commercial databases. We performed the search for similarly shaped Rapid Overlay of Chemical Structures (ROCS) and electrostatic (EON) compounds. After, we used pharmacokinetic and toxicological parameters to determine the best potential compounds, obtained through the softwares QikProp and Derek, respectively. Then, the compounds proceeded to the molecular anchorage study, which showed promising results of binding affinity with the hCOX-2 receptor: LMQC72 (∆G = −11.0 kcal/mol), LMQC36 (∆G = −10.6 kcal/mol), and LMQC50 (∆G = −10.2 kcal/mol). LMQC72 and LMQC36 showed higher binding affinity compared to rofecoxib (∆G = −10.4 kcal/mol). Finally, molecular dynamics (MD) simulations were used to evaluate the interaction of the compounds with the target hCOX-2 during 150 ns. In all MD simulation trajectories, the ligands remained interacting with the protein until the end of the simulation. The compounds were also complexing with hCOX-2 favorably. The compounds obtained the following affinity energy values: rofecoxib: ΔGbind = −45.31 kcal/mol; LMQC72: ΔGbind = −38.58 kcal/mol; LMQC36: ΔGbind = −36.10 kcal/mol; and LMQC50: ΔGbind = −39.40 kcal/mol. The selected LMQC72, LMQC50, and LMQC36 structures showed satisfactory pharmacokinetic results related to absorption and distribution. The toxicological predictions of these compounds did not display alerts for possible toxic groups and lower risk of cardiotoxicity compared to rofecoxib. Therefore, future in vitro and in vivo studies are needed to confirm the anti-inflammatory potential of the compounds selected here with bioinformatics approaches based on rofecoxib ligand

Identification of a Novel Dual Inhibitor of Acetylcholinesterase and Butyrylcholinesterase: In Vitro and In Silico Studies

The enhancement of cholinergic functions via acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) inhibition is considered a valuable therapeutic strategy for the treatment of Alzheimer’s disease. This study aimed to evaluate the in vitro effect of ZINC390718, previously filtered using computational approaches, on both cholinesterases and to characterize, using a molecular dynamics (MD) simulation, the possible binding mode of this compound inside the cholinesterase enzymes. The in vitro cytotoxicity effect was also investigated using a primary astrocyte-enriched glial cell culture. ZINC390718 presented in vitro dual inhibitory activity against AChE at a high micromolar range (IC50 = 543.8 µM) and against BuChE (IC50 = 241.1 µM) in a concentration-dependent manner, with greater activity against BuChE. The MD simulation revealed that ZINC390718 performed important hydrophobic and H-bond interactions with the catalytic residue sites on both targets. The residues that promoted the hydrophobic interactions and H-bonding in the AChE target were Leu67, Trp86, Phe123, Tyr124, Ser293, Phe295, and Tyr341, and on the BuChE target, they were Asp70, Tyr332, Tyr128, Ile442, Trp82, and Glu197. The cytotoxic effect of Z390718, evaluated via cell viability, showed that the molecule has low in vitro toxicity. The in vitro and in silico results indicate that ZINC390718 can be used as chemotype for the optimization and identification of new dual cholinesterase inhibitors

Identification of Potential New Aedes aegypti Juvenile Hormone Inhibitors from N-Acyl Piperidine Derivatives: A Bioinformatics Approach

Aedes aegypti mosquitoes transmit several human pathogens that cause millions of deaths worldwide, mainly in Latin America. The indiscriminate use of insecticides has resulted in the development of species resistance to some such compounds. Piperidine, a natural alkaloid isolated from Piper nigrum, has been used as a hit compound due to its larvicidal activity against Aedes aegypti. In the present study, piperidine derivatives were studied through in silico methods: pharmacophoric evaluation (PharmaGist), pharmacophoric virtual screening (Pharmit), ADME/Tox prediction (Preadmet/Derek 10.0®), docking calculations (AutoDock 4.2) and molecular dynamics (MD) simulation on GROMACS-5.1.4. MP-416 and MP-073 molecules exhibiting ΔG binding (MMPBSA −265.95 ± 1.32 kJ/mol and −124.412 ± 1.08 kJ/mol, respectively) and comparable to holo (ΔG binding = −216.21 ± 0.97) and pyriproxyfen (a well-known larvicidal, ΔG binding= −435.95 ± 2.06 kJ/mol). Considering future in vivo assays, we elaborated the theoretical synthetic route and made predictions of the synthetic accessibility (SA) (SwissADME), lipophilicity and water solubility (SwissADME) of the promising compounds identified in the present study. Our in silico results show that MP-416 and MP-073 molecules could be potent insecticides against the Aedes aegypti mosquitoes

Alkylated Sesamol Derivatives as Potent Antioxidants

Sesamol is a phenolic derivative. Its antioxidant activity is low than that of Trolox and depends on benzodioxole moiety. Thus, a molecular modification strategy through alkylation, inspired by natural and synthetic antioxidants, was studied by molecular modeling at the DFT/B3LYP level of theory by comparing the 6-31+G(d,p) and 6-311++G(2d,2p) basis sets. All proposed derivatives were compared to classical related antioxidants such as Trolox, t-butylated hydroxytoluene (BHT) and t-butylated hydroxyanisole (BHA). According to our results, molecular orbitals, single electron or hydrogen-atom transfers, spin density distributions, and alkyl substitutions at the ortho positions related to phenol moiety were found to be more effective than any other positions. The trimethylated derivative was more potent than Trolox. t-Butylated derivatives were stronger than all other alkylated derivatives and may be new alternative forms of modified antioxidants from natural products with applications in the chemical, pharmaceutical, and food industries

Toward of Safer Phenylbutazone Derivatives by Exploration of Toxicity Mechanism

A drug design for safer phenylbutazone was been explored by reactivity and docking studies involving single electron transfer mechanism, as well as toxicological predictions. Several approaches about its structural properties were performed through quantum chemistry calculations at the B3LYP level of theory, together with the 6-31+G(d,p) basis sets. Molecular orbital and ionization potential were associated to electron donation capacity. The spin densities contribution showed a preferential hydroxylation at the para-positions of phenyl ring when compared to other positions. In addition, on electron abstractions the aromatic hydroxylation has more impact than alkyl hydroxylation. Docking studies indicate that six structures 1, 7, 8 and 13–15 have potential for inhibiting human as well as murine COX-2, due to regions showing similar intermolecular interactions to the observed for the control compounds (indomethacin and refecoxib). Toxicity can be related to aromatic hydroxylation. In accordance to our calculations, the derivatives here proposed are potentially more active as well safer than phenylbutazone and only structures 8 and 13–15 were the most promising. Such results can explain the biological properties of phenylbutazone and support the design of potentially safer candidates

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)