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.

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.