Docking, QM/MM, and molecular dynamics simulations of the hexose transporter from Plasmodium falciparum (PfHT)

In vitro and in silico assessment of new beta amino ketones with antiplasmodial activity

BACKGROUND

Based on the current need for new drugs against malaria, our study evaluated eight beta amino ketones in silico and in vitro for potential antimalarial activity.

METHODS

Using the Brazilian Malaria Molecular Targets (BraMMT) and OCTOPUS® software programs, the pattern of interactions of beta-amino ketones was described against different proteins of P. falciparum and screened to evaluate their physicochemical properties. The in vitro antiplasmodial activities of the compounds were evaluated using a SYBR Green-based assay. In parallel, in vitro cytotoxic data were obtained using the MTT assay.

RESULTS

Among the eight compounds, compound 1 was the most active and selective against P. falciparum (IC50 = 0.98 µM; SI > 60). Six targets were identified in BraMMT that interact with compounds exhibiting a stronger binding energy than the crystallographic ligand: P. falciparum triophosphate phosphoglycolate complex (1LYX), P. falciparum reductase (2OK8), PfPK7 (2PML), P. falciparum glutaredoxin (4N0Z), PfATP6, and PfHT.

CONCLUSIONS

The physicochemical properties of compound 1 were compatible with the set of criteria established by the Lipinski rule and demonstrated its potential as a drug prototype for antiplasmodial activity.

An overview of the Plasmodium falciparum hexose transporter and its therapeutic interventions

Despite intense elimination efforts, human malaria, caused by the infection of five Plasmodium species, remains the deadliest parasitic disease in the world. Even worse, with the emergence and spreading of the first-line drug-resistant Plasmodium parasites, therapeutic interventions based on novel plasmodial drug targets are more necessary than ever. Given that the blood-stage parasites primarily rely on glycolysis for their energy supply, blocking glucose uptake, the rate-limiting step of ATP generation, was considered a promising approach to kill these parasites. To achieve this goal, characterization of the plasmodial hexose transporter and development of selective inhibitors have been pursued for decades. Here, we review the identification and characterization of the Plasmodium falciparum hexose transporter (PfHT1) and summarize current advances in its inhibitor development.

Molecular docking, simulation and binding free energy analysis of small molecules as PfHT1 inhibitors

Antimalarial drug resistance has thrown a spanner in the works of malaria elimination. New drugs are required for ancillary support of existing malaria control efforts. Plasmodium falciparum requires host glucose for survival and proliferation. On this basis, P. falciparum hexose transporter 1 (PfHT1) protein involved in hexose permeation is considered a potential drug target. In this study, we tested the antimalarial activity of some compounds against PfHT1 using computational techniques. We performed high throughput virtual screening of 21,352 small-molecule compounds against PfHT1. The stability of the lead compound complexes was evaluated via molecular dynamics (MD) simulation for 100 nanoseconds. We also investigated the pharmacodynamic, pharmacokinetic and physiological characteristics of the compounds in accordance with Lipinksi rules for drug-likeness to bind and inhibit PfHT1. Molecular docking and free binding energy analyses were carried out using Molecular Mechanics with Generalized Born and Surface Area (MMGBSA) solvation to determine the selectivity of the hit compounds for PfHT1 over the human glucose transporter (hGLUT1) orthologue. Five important PfHT1 inhibitors were identified: Hyperoside (CID5281643); avicularin (CID5490064); sylibin (CID5213); harpagoside (CID5481542) and quercetagetin (CID5281680). The compounds formed intermolecular interaction with the binding pocket of the PfHT1 target via conserved amino acid residues (Val314, Gly183, Thr49, Asn52, Gly183, Ser315, Ser317, and Asn48). The MMGBSA analysis of the complexes yielded high free binding energies. Four (CID5281643, CID5490064, CID5213, and CID5481542) of the identified compounds were found to be stable within the PfHT1 binding pocket throughout the 100 nanoseconds simulation run time. The four compounds demonstrated higher affinity for PfHT1 than the human major glucose transporter (hGLUT1). This investigation demonstrates the inhibition potential of sylibin, hyperoside, harpagoside, and avicularin against PfHT1 receptor. Robust preclinical investigations are required to validate the chemotherapeutic properties of the identified compounds.

Evaluation of antiplasmodial activity in silico and in vitro of N-acylhydrazone derivatives

N-acylhydrazones are considered privileged structures in medicinal chemistry, being part of antimicrobial compounds (for example). In this study we show the activity of N-acylhydrazone compounds, namely AH1, AH2, AH4, AH5 in in vitro tests against the chloroquine-resistant strain of Plasmodium falciparum (W2) and against WI26 VA-4 human cell lines. All compounds showed low cytotoxicity (LC₅₀ > 100 µM). The AH5 compound was the most active against Plasmodium falciparum, with an IC₅₀ value of 0.07 μM. AH4 and AH5 were selected among the tested compounds for molecular docking calculations to elucidate possible targets involved in their mechanism of action and the SwissADME analysis to predict their pharmacokinetic profile. The AH5 compound showed affinity for 12 targets with low selectivity, while the AH4 compound had greater affinity for only one target (3PHC). These compounds met Lipinski's standards in the ADME in silico tests, indicating good bioavailability results. These results demonstrate that these N-acylhydrazone compounds are good candidates for future preclinical studies against malaria.

A simple quinoline salt derivative is active in vitro against plasmodium Faciparum asexual blood stages and inhibits the development of cerebral malaria in murine model

Chloroquine (CQ) was the most effective and widely used drug for the prophylaxis and treatment of severe and non-severe malaria. Although its prophylactic use has led to resistance to P. falciparum in all endemic countries, CQ still remains the drug of choice for the treatment of vivax malaria. Otherwise, the speed in which parasite resistance to available antimalarials rises and spreads in endemic regions points to the urgent need for the development of new antimalarials. Quinoline derivatives have been used as a tool in the search for new drugs and were investigated in the present study in an attempt to produce a HIT compound to avoid the cerebral malarial (CM). Seven compounds were synthesized, including three quinoline derivate salts. The cytotoxicity and antiplasmodial activity were assayed in vitro, highlighting compound 3 as a HIT, which also showed interaction with ferriprotoporphyrin IX similarly to CQ. Physicochemical and pharmacokinetic properties of absorption were found to be favorable when analyzed in silico. The in vivo assays, using the experimental cerebral malaria (ECM) model, showed important values of parasite growth inhibition on the 7th day-post infection (Q15 15 mg/kg: 76.9%, Q30 30 mg/kg: 90,1% and Q50 50 mg/kg: 92,9%). Compound 3 also showed significant protection against the development of CM, besides hepatic and renal parameters better than CQ. In conclusion, this quinoline derivative demonstrated promising activity for the treatment of malaria and was able to avoid the development of severe malaria in mice.

Dehydrobufotenin extracted from the Amazonian toad Rhinella marina (Anura: Bufonidae) as a prototype molecule for the development of antiplasmodial drugs

Background:

The resistance against antimalarial drugs represents a global challenge in the fight and control of malaria. The Brazilian biodiversity can be an important tool for research and development of new medicinal products. In this context, toxinology is a multidisciplinary approach on the development of new drugs, including the isolation, purification, and evaluation of the pharmacological activities of natural toxins. The present study aimed to evaluate the cytotoxicity, as well as the antimalarial activity in silico and in vitro of four compounds isolated from Rhinella marina venom as potential oral drug prototypes.

Methods:

Four compounds were challenged against 35 target proteins from P. falciparum and screened to evaluate their physicochemical properties using docking assay in Brazilian Malaria Molecular Targets (BraMMT) software and in silico assay in OCTOPUS® software. The in vitro antimalarial activity of the compounds against the 3D7 Plasmodium falciparum clones were assessed using the SYBR Green I based assay (IC₅₀). For the cytotoxic tests, the LD₅₀ was determined in human pulmonary fibroblast cell line using the [3(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide] (MTT) assay.

Results:

All compounds presented a ligand-receptor interaction with ten Plasmodium falciparum-related protein targets, as well as antimalarial activity against chloroquine resistant strain (IC₅₀ = 3.44 μM to 19.11 μM). Three of them (dehydrobufotenine, marinobufagin, and bufalin) showed adequate conditions for oral drug prototypes, with satisfactory prediction of absorption, permeability, and absence of toxicity. In the cell viability assay, only dehydrobufotenin was selective for the parasite.

Conclusions:

Dehydrobufotenin revealed to be a potential oral drug prototype presenting adequate antimalarial activity and absence of cytotoxicity, therefore should be subjected to further studies.

In silico approach: biological prediction of nordentatin derivatives as anticancer agent inhibitors in the cAMP pathway

A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene. The inhibition of the cAMP pathway is the main target of anti-cancer drugs, which is responsible for uncontrolled cell division in cancer. Modeling was done using a combination of semi-empirical methods and the density functional theory (PM3-DFT/6-31G*/B3LYP) to obtain the optimal structure of a small ligand that could be modeled. Studies on the interaction of the ligands and amino acid residues on protein targets were carried out using a combination of molecular docking and molecular dynamic simulation. Molecular docking based on functional grid scores showed a very good native ligand pose with an RMSD of 0.93 Å in determining the initial coordinates of the ligand-receptor interactions. Furthermore, the amino acid residues responsible for interaction through H-bonds were Tyr103, His104, His177, Met217, and Gln313. The binding free energy (kcal mol-1) results of the candidates were PS-1 (-36.84 ± 0.31), PS-2 (-35.34 ± 0.28), PS-3 (-26.65 ± 0.30), PS-5 (-42.66 ± 0.26), PS-7 (-35.33 ± 0.23), and PS-9 (-32.57 ± 0.20), which are smaller than that of the native ligand Z72 (-24.20 ± 0.19), and thus these have good potential as drugs that can inhibit the cAMP pathway. These results provide theoretical information for the efficient inhibition of the cAMP pathway in the future.

Brazilian malaria molecular targets (BraMMT): selected receptors for virtual high-throughput screening experiments

BACKGROUND

Owing to increased spending on pharmaceuticals since 2010, discussions about rising costs for the development of new medical technologies have been focused on the pharmaceutical industry. Computational techniques have been developed to reduce costs associated with new drug development. Among these techniques, virtual high-throughput screening (vHTS) can contribute to the drug discovery process by providing tools to search for new drugs with the ability to bind a specific molecular target.

OBJECTIVES

In this context, Brazilian malaria molecular targets (BraMMT) was generated to execute vHTS experiments on selected molecular targets of Plasmodium falciparum.

METHODS

In this study, 35 molecular targets of P. falciparum were built and evaluated against known antimalarial compounds.

FINDINGS

As a result, it could predict the correct molecular target of market drugs, such as artemisinin. In addition, our findings suggested a new pharmacological mechanism for quinine, which includes inhibition of falcipain-II and a potential new antimalarial candidate, clioquinol.

MAIN CONCLUSIONS

The BraMMT is available to perform vHTS experiments using OCTOPUS or Raccoon software to improve the search for new antimalarial compounds. It can be retrieved from www.drugdiscovery.com.br or download of Supplementary data.

Engineering of lysine cyclodeaminase conformational dynamics for relieving substrate and product inhibitions in the biosynthesis of l-pipecolic acid

Efficient biocatalytic process construction by relieving substrate and product inhibitions via identification and engineering of enzyme conformational dynamics.

Structural dynamics of Casein Kinase I (CKI) from malarial parasite Plasmodium falciparum (Isolate 3D7): Insights from theoretical modelling and molecular simulations

The protein kinases (PKs), belonging to serine/threonine kinase (STKs), are important drug targets for a wide spectrum of diseases in human. Among protein kinases, the Casein Kinases (CKs) are vastly expanded in various organisms, where, the malarial parasite Plasmodium falciparum possesses a single member i.e., PfCKI, which can phosphorylate various proteins in parasite extracts in vitro condition. But, the structure-function relationship of PfCKI and dynamics of ATP binding is yet to be understood. Henceforth, an attempt was made to study the dynamics, stability, and ATP binding mechanisms of PfCKI through computational modelling, docking, molecular dynamics (MD) simulations, and MM/PBSA binding free energy estimation. Bi-lobed catalytic domain of PfCKI shares a high degree of secondary structure topology with CKI domains of rice, human, and mouse indicating co-evolution of these kinases. Molecular docking study revealed that ATP binds to the active site where the glycine-rich ATP-binding motif (G16-X-G18-X-X-G21) along with few conserved residues plays a crucial role maintaining stability of the complex. Structural superposition of PfCKI with close structural homologs depicted that the location and length of important loops are different, indicating the dynamic properties of these loops among CKIs, which is consistent with principal component analysis (PCA). PCA displayed that the overall global motion of ATP-bound form is comparatively higher than that of apo form. The present study provides insights into the structural features of PfCKI, which could contribute towards further understanding of related protein structures, dynamics of catalysis and phosphorylation mechanism in these important STKs from malarial parasite in near future.