Towards development of new antimalarial compounds through in silico and in vitro assays.

Malaria is one of most widespread infectious disease in world. The antimalarial therapy presents a series of limitations, such as toxicity and the emergence of resistance, which makes the search for new drugs urgent. Thus, it becomes necessary to explore essential and exclusive therapeutic targets of the parasite to achieve selective inhibition. Enoyl-ACP reductase is an enzyme of the type II fatty acid biosynthetic pathway and is responsible for the rate-limiting step in the fatty acid elongation cycle. In this work, we use hierarchical virtual screening and drug repositioning strategies to prioritize compounds for phenotypic assays and molecular dynamics studies. The molecules were tested against chloroquine-resistant W2 strain of Plasmodium falciparum (EC50 between 330.05 and 13.92 µM). Nitrofurantoin was the best antimalarial activity at low micromolar range (EC50 = 13.92 µM). However, a hit compound against malaria must have a biological activity value below 1 µM. A large number of molecules present problems with permeability in biological membranes and reaching an effective concentration in their target's microenvironment. Nitrofurantoin derivatives with inclusions of groups which confer increased lipid solubility (methyl groups, halogens and substituted and unsubstituted aromatic rings) have been proposed. These derivatives were pulled through the lipid bilayer in molecular dynamics simulations. Molecules 14, 18 and 21 presented lower free energy values than nitrofurantoin when crossing the lipid bilayer.

Identification of Potential Multitarget Compounds against Alzheimer's Disease through Pharmacophore-Based Virtual Screening.

Alzheimer's disease (AD) is a neurodegenerative disease characterized by progressive loss of cognitive functions, and it is the most prevalent type of dementia worldwide, accounting for 60 to 70% of cases. The pathogenesis of AD seems to involve three main factors: deficiency in cholinergic transmission, formation of extracellular deposits of ß-amyloid peptide, and accumulation of deposits of a phosphorylated form of the TAU protein. The currently available drugs are prescribed for symptomatic treatment and present adverse effects such as hepatotoxicity, hypertension, and weight loss. There is urgency in finding new drugs capable of preventing the progress of the disease, controlling the symptoms, and increasing the survival of patients with AD. This study aims to present new multipurpose compounds capable of simultaneously inhibiting acetylcholinesterase (AChE), butyrylcholinesterase (BChE)-responsible for recycling acetylcholine in the synaptic cleft-and beta-secretase 1 (BACE-1)-responsible for the generation of amyloid-ß plaques. AChE, BChE, and BACE-1 are currently considered the best targets for the treatment of patients with AD. Virtual hierarchical screening based on a pharmacophoric model for BACE-1 inhibitors and a dual pharmacophoric model for AChE and BChE inhibitors were used to filter 214,446 molecules by QFITBACE > 0 and QFITDUAL > 56.34. The molecules selected in this first round were subjected to molecular docking studies with the three targets and further evaluated for their physicochemical and toxicological properties. Three structures: ZINC45068352, ZINC03873986, and ZINC71787288 were selected as good fits for the pharmacophore models, with ZINC03873986 being ultimately prioritized for validation through activity testing and synthesis of derivatives for SAR studies.

Molecular Multi-Target Approach for Human Acetylcholinesterase, Butyrylcholinesterase and ß-Secretase 1: Next Generation for Alzheimer's Disease Treatment.

Alzheimer's Disease (AD) is a neurodegenerative condition characterized by progressive memory loss and other affected cognitive functions. Pharmacological therapy of AD relies on inhibitors of the enzymes acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), offering only a palliative effect and being incapable of stopping or reversing the neurodegenerative process. However, recent studies have shown that inhibiting the enzyme ß-secretase 1 (BACE-1) may be able to stop neurodegeneration, making it a promising target. Considering these three enzymatic targets, it becomes feasible to apply computational techniques to guide the identification and planning of molecules capable of binding to all of them. After virtually screening 2119 molecules from a library, 13 hybrids were built and further screened by triple pharmacophoric model, molecular docking, and molecular dynamics (t = 200 ns). The selected hybrid G meets all stereo-electronic requirements to bind to AChE, BChE, and BACE-1 and offers a promising structure for future synthesis, enzymatic testing, and validation.

Antiplasmodial activity of coumarins isolated from Polygala boliviensis: in vitro and in silico studies.

Polygala boliviensis is found in the Brazilian semiarid region. This specie is little chemically and biologically studied. Polygala spp. have different metabolites, especially coumarins. Studies indicate that coumarins have antimalarial potential, denoting the importance of researching new active compounds from plants, since the resistance of Plasmodium strains to conventional therapy has increased. The present study aimed to evaluate the antiplasmodial activity of auraptene and poligalen against a chloroquine-resistant strain of Plasmodium falciparum. Coumarins were isolated from P. boliviensis by open column chromatography and identified by Nuclear Magnetic Resonance Spectroscopy. A cytotoxicity assay was carried out using MTT test, and the in vitro antiplasmodial activity was evaluated using the W2 strain. The antiplasmodial activity results found were IC50=0.171 ± 0.016 for auraptene and 0.164 ± 0.012 for poligalen; the selectivity indexes were 78.71 and 609.76, respectively. Inverse virtual screening in the BRAMMT database by OCTOPUS 1.2 was applied to coumarins to find potential P. falciparum targets and showed higher affinity energy of auraptene for purine nucleoside phosphorylase (PfPNP) and of poligalen for dihydroorotate dehydrogenase (PfDHODH). Molecular Dynamics studies (MD and MM-GBSA) approach were applied to calculate binding energies against selected P. falciparum targets and showed that all coumarins were stable at the binding site during simulations. Furthermore, energies were favorable for complexation. This is the first report of auraptene in P. boliviensis species and of in vitro antiplasmodial activity of auraptene and poligalen. In silico studies indicated that the mechanism of action of coumarins is the inhibition of PfPNP and PfDHODH.Communicated by Ramaswamy H. Sarma.

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.

Identification of novel Aedes aegypti odorant-binding protein 1 modulators by ligand and structure-based approaches and bioassays.

Arboviruses are a group of viruses (e.g. Dengue, Chikungunya and Yellow fever virus) that are transmitted by arthropod vectors, which Aedes aegipty is the vector of main viruses in Americas. This vector is responsible to 2.4 millions of arboviruses cases in Brazil with less than a thousand deaths annually. Despite of epidemiological data, arboviruses treatment is symptomatic and the vaccine control is not effective, which makes the vector control against A. aegipty a promising strategy to diseases control. One way to achieve this goal is to development of A. aegipty sensitive olfactory modulators. Odorant binding protein 1 from A. aegypti (AaegOBP1) is essential in sensory communication, and is the first filter in odorant selection, which makes this target promising to development of new repellents. For this reason, hierarchical virtual screening (ligand-based pharmacophore model and molecular docking) together volatility filter was applied at Sigma-Aldrich database (n = 126.851) to prioritize potential molecules to repellency assays. Three compounds showed adequate stereo-electronic requirements (QFIT> 81.53), score to AaegOBP1 binding site (Score > 36.0) and volatile properties and it was chosen for repellency assays. ZINC00170981 and ZINC00131924 showed a dose-response behavior, while ZINC01621824 did not showed activity in repellency assays. Finally, Molecular Dynamics (MD) was employed to hypothesize the stability of protein-ligand complexes. According to RMSD, RMSF and binding free energy data, ZINC00170981 and ZINC00131924 were able to stabilize AaegOBP1 binding-site during the trajectory by interactions with key residues such as His77, Leu89 and Trp114). Communicated by Ramaswamy H. Sarma.

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.

Functional annotation and comparative modeling of ligninolytic enzymes from Trametes villosa (SW.) Kreisel for biotechnological applications.

Functional annotation of Trametes villosa genome was performed to search Class II peroxidase proteins in this white-rot fungus, which can be valuable for several biotechnological processes. After sequence identification and manual curation, five proteins were selected to build 3 D models by comparative modeling. Analysis of sequential and structural sequences from selected targets revealed the presence of two putative Lignin Peroxidase and three putative Manganese Peroxidase on this fungal genome. All 3 D models had a similar folding pattern from selected 3 D structure templates. After minimization and validation steps, the best 3 D models were subjected to docking studies and molecular dynamics to identify structural requirements and the interactions required for molecular recognition. Two reliable 3 D models of Class II peroxidases, with typical catalytic site and architecture, and its protein sequences are indicated to recombinant production in biotechnological applications, such as bioenergy.Communicated by Ramaswamy H. Sarma.

A Review on the In Vitro Evaluation of the Anticholinesterase Activity Based on Ellman's Method.

Inhibition of cholinesterases is a common strategy for the treatment of several disorders, especially Alzheimer´s disease. In vitro assays represent a critical step towards identifying molecules with potential anticholinesterase effect. This study aimed at providing a comprehensive review of the methodologies used in vitro for the anticholinesterase activity based on the spectrophotometry of Ellman's method. This work used two databases (PubMed and ScienceDirect) to search for original articles and selected publications between 1961 and 2019, which reported in vitro spectrophotometry assays for anticholinesterase activity. After the search process and the selection of publications, the final sample consisted of 146 articles published in several journals submitted by researchers from different countries. Although the studies analyzed in this work are all within the same conception of in vitro tests based on Ellman's method, one can observe a wide divergence in the origin and concentration of enzyme, the choice and pH of the buffer, the concentration of the substrate, the sample diluent, incubation time, temperature, and time of the spectrophotometric reading interval. There is no consensus in the methodology of studies with in vitro tests for anticholinesterase assessment. The methodological variations related to kinetic parameters may interfere in the characterization of cholinesterase inhibitors.

In silico development of adenosine A2B receptor antagonists for sickle cell disease.

Sickle cell disease (SCD) is a disease resulting from mutation in the globin portion of hemoglobin caused by the replacement of adenine for thymine in the codon of the ß globin gene. In Brazil, SCD affects about 0.3% of the black and Caucasian population. Until now, there is no specific treatment and the available drugs have several serious adverse effects which makes the search for new drugs an emergently need. The use of computational techniques can accelerate the drug development process by prioritization of molecules with affinity against essential targets. Adenosine A2b receptor (rA2b) has been studied in SCD due to its relationship with red blood cells concentration of 2,3-diphosphoglycerate which reduces the hemoglobin affinity for oxygen (O2), facilitating its availability for the tissues. Then, development of rA2b antagonists could be helpful for the treatment of SCD. However, there is still no 3D structure of rA2b and to overcome this limitation, homology modeling should be applied. In this scenario, this study aims to build a suitable 3D model of rA2b by SWISS MODEL and to evaluate the structural aspects of rA2b with known antagonists that may be useful for the identification of new potential antagonists by molecular dynamics on a lipid bilayer environment using GROMACS 5.1.4. The complexes with antagonists ZINC223070016 and ZINC17974526 interacted with key residues by hydrophobic contacts and hydrogen bonds which stabilized them at the rA2b binding site. This intermolecular profile can contribute to the development of more potent rA2b antagonists. 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.

Pharmacophore-based virtual screening and molecular docking to identify promising dual inhibitors of human acetylcholinesterase and butyrylcholinesterase.

The dual inhibition of human acetylcholinesterase (hAChE) and butyrylcholinesterase (hBuChE) plays an important role in Alzheimer's disease treatment. Thus, this study aims identify promising dual inhibitors against hAChE and hBuChE by in silico approaches (pharmacophore-based virtual screening and molecular docking). Ten 3 D pharmacophore models for dual inhibitors using default genetic parameters were built by GALAHAD™ available on SYBYL-X 2.0. Validation steps were carried out according to Energy (<100.0 kcal/mol), Pareto = 0, Area under the ROC Curve (>0.70), Boltzmann-Enhanced Discrimination of ROC curve (BEDROC >0.50) and structure-activity relationship (SAR) for known inhibitors. The best dual pharmacophore model based on internal/external statistical parameters and SAR data (one hydrogen bond acceptor, two hydrogen bond donors and four hydrophobic centers) was employed in virtual screening at Sigma-Aldrich® subset (n = 214,446) of ZINC database by UNITY module of SYBYL-X 2.0. According to superposition values (QFIT), the best ranked compounds were prioritized for molecular docking and partition coefficient analysis (clog p < 5.0). 37 top-ranked compounds (QFIT > 64.22) from pharmacophore model showed affinity in hAChE (-10.2 < Affinity energy < -6.3 kcal/mol) and hBuChE (-10.9 < Affinity energy < -2.3 kcal/mol) binding sites. Next, liposolubity prediction and commercially available showed that ZINC43198636, ZINC43198637 and ZINC00390718 can be potential dual inhibitors against hAChE and hBuChE.Communicated by Ramaswamy H. Sarma.

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.

Identification of potential modulator of Anopheles gambiae odorant binding protein 1 by hierarchical virtual screening and molecular dynamics.

Malaria is a protozoan infection transmitted by the bite of the infected female mosquito belonging to the genus Anopheles spp., which causes more than 445 million annual deaths worldwide. Available drugs have serious adverse effects (e.g. blurred vision, hypotension and headache) and species-dependent efficacy. An alternative to overcome these problems involve the use of molecules with affinity to the Anopheles gambiae mosquito odor receptors, minimizing the reinfection process as well as reducing the problems related to pharmacological therapy. The vector control can interrupt the epidemiological cycle and, therefore, control the malaria incidence. In the olfactory pathway, odorant binding protein 1 acts on the first level of odor recognition on malarial vector and thus can be used to modulate mosquito behavior and development of new attracts or repellents. Thus, this study applied ligand-based (2D-chemical similarity) and structure-based (docking and molecular dynamics) computational approaches to prioritize potential olfactory modulators on natural products catalogs at ZINC15 database (n = 98,379). Hierarchical virtual screening prioritized a potential olfactory modulator (Z8217) against Anopheles gambiae odorant binding protein 1 (AgOBP1). Next, it was submitted to molecular dynamics routine to identify structural requirements and the interactions profile required for binding-site affinity. This promising natural compound can interact like experimental ligand and will be used in repellency assay to confirm its sensorial behavior.Communicated by Ramaswamy H. Sarma.

Identification of potential Leishmania chagasi superoxide dismutase allosteric modulators by structure-based computational approaches: homology modelling, molecular dynamics and pharmacophore-based virtual screening.

The visceral form of Leishmaniasis, also known as kala-azar, caused by Leishmania chagasi is the main etiological agent of this form in Brazil responsible for 30,000 annual deaths. Despite its epidemiological impact, treatment of the disease is limited by resistance, species-dependent efficacy and serious adverse effects. The application of computational tools to prioritize potential bioactive molecules based on 3D structural of biological target is a viable alternative. Among the L. chagasi validated targets, Fe + 2 superoxide dismutase B2 (LcFeSODB2) is the first parasite enzyme against oxidative stress and it is involved in essential metabolic processes for its survival. Due to substrate binding-site volume (superoxide ion) and consequent difficulty in its active site modulation for small molecules, the search for allosteric sites at LcFeSODB2 3D structure is a promising strategy. As there are no 3D structures of LcFeSODB2, comparative modeling was applied to build 3D models by SWISS-MODEL and MODELLER version 9.19. Next, the best 3D model was used in molecular dynamics (MD) routines with multiple probes on GROMACS version 5.1.2. In addition, potential allosteric sites predicted by FTMap and Metapocket web servers were used with probe occupancy maps from MD to select an allosteric binding site and propose a pharmacophore model. Next, it was used as a template in virtual screening by UNITY® module available on SYBYL-X version 2.1.1 at Sigma-Aldrich CPR™ subset of ZINC12 database. The pharmacophore-based virtual screening resulted in the selection of two potential allosteric LcFeSOD compounds with partial pharmacophoric requirements, drug-like properties and commercial availability for enzymatic assays. Communicated by Ramaswamy H. Sarma.

Antileishmanial activity evaluation of thiazolidine-2,4-dione against Leishmania infantum and Leishmania braziliensis.

Leishmaniasis is responsible for approximately 65,000 annual deaths. Despite the mortality data, drugs available for the treatment of patients are insufficient and have moderate therapeutic efficacy in addition to serious adverse effects, which makes the development of new drugs urgent. To achieve this goal, the integration of kinetic and DSF assays against parasitic validated targets, along with phenotypic assays, can help the identification and optimization of bioactive compounds. Pteridine reductase 1 (PTR1), a validated target in Leishmania sp., is responsible for the reduction of folate and biopterin to tetrahydrofolate and tetrahydrobiopterin, respectively, both of which are essential for cell growth. In addition to the in vitro evaluation of 16 thiazolidine-2,4-dione derivatives against Leishmania major PTR1 (LmPTR1), using the differential scanning fluorimetry (ThermoFluor®), phenotypic assays were employed to evaluate the compound effect over Leishmania braziliensis (MHOM/BR/75/M2903) and Leishmania infantum (MHOM/BR/74/PP75) promastigotes viability. The ThermoFluor® results show that thiazolidine-2,4-dione derivatives have micromolar affinity to the target and equivalent activity on Leishmania cells. 2b is the most potent compound against L. infantum (EC50 = 23.45 ± 4.54 µM), whereas 2a is the most potent against L. braziliensis (EC50 = 44.16 ± 5.77 µM). This result suggests that lipophilic substituents on either-meta and/or-para positions of the benzylidene ring increase the potency against L. infantum. On the other hand, compound 2c (CE50 = 49.22 ± 7.71 µM) presented the highest selectivity index.

Pharmacophore modeling, docking and molecular dynamics to identify Leishmania major farnesyl pyrophosphate synthase inhibitors.

Leishmaniasis is caused by protozoa of the genus Leishmania spp. and is considered the second most important protozoa in the world due to the number of cases and mortality. Despite its importance in terms of public health, the treatment of patients is limited and has mostly low levels of efficacy and safety. Farnesyl pyrophosphate synthase (FPPS) acts in the early stages of isoprenoid synthesis, and is important for maintaining the integrity of the lipid bilayer of the parasite that causes the disease. The aim of this work was to identify one potential inhibitor of the FPPS of Leishmania major through virtual screening by pharmacophore modeling and docking. A total of 85,000 compounds from a natural products database (ZINC15) was submitted for virtual hierarchical screening, and the top ranked molecule in both methods was analyzed by intermolecular interaction profile and 20 ns molecular dynamics simulations. These results showed a promising compound from natural products that mimic the major interactions present in the substrate/inhibitor.

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.

An integrated approach towards the discovery of novel non-nucleoside Leishmania major pteridine reductase 1 inhibitors.

Despite the fact that Leishmania ssp are pteridine auxotrophs, Dihydrofolate Reductase-Thymidylate Synthase (DHFR-TS) inhibitors are ineffective against Leishmania major. On the other hand Pteridine Reductase 1 (PTR1) inhibitors proved to be lethal to the parasite. Aiming at identifying hits that lie outside the chemical space of known PTR1 inhibitors, pharmacophore models that differentiate true-binders from decoys and explain the structure-activity relationships of known inhibitors were employed to virtually screen the lead-like subset of ZINC database. This approach leads to the identification of Z80393 (IC50 = 32.31 ± 1.18 µM), whose inhibition mechanism was investigated by Thermal Shift Assays. This experimental result supports a competitive mechanism and was crucial to establish the docking search space as well as select the best pose, which was then investigated by molecular dynamics studies that corroborate the hit putative binding profile towards LmPTR1. The information gathered from such studies shall be useful to design more potent non-nucleoside LmPTR1 inhibitors.

Caffeoylquinic acids from antiplasmodial active extract of Xanthium cavanillesii fruits and their molecular modelling studies.

The antiplasmodial active extract of Xanthium cavanillesii contains 3,4-dicaffeoyl quinic acid (3,4-DCQA), 3,5-dicaffeoyl quinic acid (3,5-DCQA) and 1,3,5-tricaffeoyl quinic acid (1,3,5-TCQA). These results inspired us to investigate the interaction of these molecules with a promising validated target of Plasmodium, PfATP6 orthologue of mammalian Ca+2-ATPase. Models of this receptor were obtained through comparative modelling. Afterwards, molecular docking studies were used to identify possible interaction modes of these caffeoyl quinic derivatives on the binding site. The 1,3,5-TCQA had the best energy, but all of these had better energy than thapsigargin, a non-competitive inhibitor of the sarco/endoplasmatic reticulum Ca+2-ATPase (SERCA).