In search of the representative pharmacophore hypotheses of the enzymatic proteome of Plasmodium falciparum: a multicomplex-based approach

Analysis of the interaction of antimalarial agents with Plasmodium falciparum glutathione reductase through molecular mechanical calculations

CONTEXT

Malaria remains a significant global health challenge with emerging resistance to current treatments. Plasmodium falciparum glutathione reductase (PfGR) plays a critical role in the defense mechanisms of malaria parasites against oxidative stress. In this study, we investigate the potential of targeting PfGR with conventional antimalarials and dual drugs combining aminoquinoline derivatives with GR inhibitors, which reveal promising interactions between PfGR and studied drugs. The naphthoquinone Atovaquone demonstrated particularly high affinity and potential dual-mode binding with the enzyme active site and cavity. Furthermore, dual drugs exhibit enhanced binding affinity, suggesting their efficacy in inhibiting PfGR, where the aliphatic ester bond (linker) is essential for effective binding with the enzyme's active site. Overall, this research provides important insights into the interactions between antimalarial agents and PfGR and encourages further exploration of its role in the mechanisms of action of antimalarials, including dual drugs, to enhance antiparasitic efficacy.

METHODS

The drugs were tested as PfGR potential inhibitors via molecular docking on AutoDock 4, which was performed based on the preoptimized structures in HF/3-21G-PCM level of theory on ORCA 5. Drug-receptor systems with the most promising binding affinities were then studied with a molecular dynamic's simulation on AMBER 16. The molecular dynamics simulations were performed with a 100 ns NPT ensemble employing GAFF2 forcefield in the temperature of 310 K, integration time step of 2 fs, and non-bond cutoff distance of 6.0 Å.

Insights into the selective mechanism of PDE2/9a inhibitors from silico aspects

The selective design of competitive enzyme inhibitors is an extremely difficult task but necessary work for certain types of systems, such as the phosphodiesterase (PDE) system addressed in this article. In the PDE family, PDE2A and PDE9 respectively target the central nervous system and heart failure, and share many conserved amino acids at their binding sites. Therefore, gaining a deep understanding of the selective mechanisms of PDE2A/9A is crucial for designing highly selective drugs. In this study, various computer-aided drug design (CADD) methods, including molecular docking, molecular dynamics simulations (MD), and binding free energy calculations, are employed to explore the selective mechanisms of PDE2A/9A. Overall, our research results indicate a selective design strategy for PDE2A, which involves incorporating hydrophobic or aromatic moieties into the molecular structure to better accommodate the hydrophobic pocket of PDE2A. Additionally, it is recommended to introduce functional groups capable of forming connections with selective residues, such as Phe830 and Gln812 for PDE2A, or Ala452 and Tyr424 for PDE9A, to enhance the selectivity of inhibitors targeting PDE2A/9A. This achievement is anticipated to pave the way for the development of innovative and selective small molecules targeting PDE2A/9A.Communicated by Ramaswamy H. Sarma.

Targeting allosteric binding site in methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) to identify natural product inhibitors via structure-based computational approach

Cancer has been viewed as one of the deadliest diseases worldwide. Among various types of cancer, breast cancer is the most common type of cancer in women. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is a promising druggable target and is overexpressed in cancerous cells, like, breast cancer. We conducted structure-based modeling on the allosteric site of the enzyme. Targeting the allosteric site avoids the problem of drug resistance. Pharmacophore modeling, molecular docking, HYDE assessment, drug-likeness, ADMET predictions, simulations, and free-energy calculations were performed. The RMSD, RMSF, RoG, SASA, and Hydrogen-bonding studies showed that seven candidates displayed stable behaviour. As per the literature, average superimposed simulated structures revealed a similar protein conformational change in the αE'-βf' loop, causing its displacement away from the allosteric site. The MM-PBSA showed tight binding of six compounds with the allosteric pocket. The effect of inhibitors interacting in the allosteric site causes a decrease in the binding energy of J49 (active-site inhibitor), suggesting the effect of allosteric binding. The PCA and FEL analysis revealed the significance of the docked compounds in the stable behaviour of the complexes. The outcome can contribute to the development of potential natural products with drug-like properties that can inhibit the MTHFD2 enzyme.

Computational studies to explore inhibitors against the cyclin-dependent kinase 12/13 enzyme: an insilco pharmacophore modeling, molecular docking and dynamics approach

Cancer is enlisted among the deadliest disease all over the world. The cyclin-dependent kinases 12 and 13 have been identified as cell cycle regulators. They conduct transcription and co-transcriptional processes by phosphorylating the C-terminal of RNA polymerase-II. Inhibition of CDK12 and 13 selectively presents a novel strategy to treat triple-negative breast cancer, but dual inhibitors are still lacking. Here, we report the screening of the natural product compound class against the dual CDK12/13 enzyme by employing various in silico methods. Complexes of CDK12 enzymes are used to form common feature pharmacophore models, whereas we perform receptor-based pharmacophore modelling on CDK13 enzyme owing to the availability of a single PDB. On conducting screening over the representative pharmacophores, the common drug-like screened natural products were shortlisted for conducting molecular docking studies. After molecular docking calculations, the candidates that showed crucial interaction with CDK12 and CDK13 enzymes were shortlisted for simulation studies. Five common docked candidates were selected for molecular dynamics simulations and free energy calculations. Based on the cut-off criteria of free energy calculations, one common hit was selected as the dual CDK12/13 inhibitor. The outcome concluded that the hit with ID CNP0386383 possesses drug-like properties, displays crucial interaction in the binding pocket, and shows stable dynamic behaviour and higher binding energy than the experimentally reported inhibitor of both CDK12 and CDK13 enzymes.Communicated by Ramaswamy H. Sarma.

Identification of natural products against enoyl-acyl-carrier-protein reductase in malaria via combined pharmacophore modeling, molecular docking and simulations studies

Plasmodium falciparum is counted as one of the deadly species causing malaria. In that respect, enoyl acyl carrier protein reductase is recognized as one of the attractive druggable targets for the identification of antimalarials. Thus, from the structural proteome of ENR, common feature pharmacophores were constructed. To identify the representative models, all the hypotheses were subjected to validation methods, like, test set, enrichment factor, and Güner-Henry method, and the selected representative hypotheses were used to screen out the drug-like natural products. Further, the screened candidates were advanced to molecular docking calculations. Based on the docking score criteria and presence of essential interaction with Tyr277, seven candidates were shortlisted to conduct the HYDE and QSAR assessment. Further, the stability of these complexes was evaluated by employing molecular dynamics simulations, molecular mechanics-generalized born surface area approach-based free binding energy calculations with the residue-wise contribution of PfENR to the total binding free energy of the complex. On comparing the root mean square deviation, and fluctuation plots of the docked candidates with the reference, all the candidates displayed stable behavior, and the same outcome was depicted from the secondary structure element. However, from the free energy calculations, and residue-wise contribution conducted after dynamics, it was observed that out of seven, only five candidates sustain the binding with Tyr277 and cofactor of PfENR. Therefore, in the current work, the hybrid study of screening and stability lead to the identification of five structurally diverse candidates that can be employed for the design of novel antimalarials.Communicated by Ramaswamy H. Sarma.

Structure- and ligand-based drug design methods for the modeling of antimalarial agents: a review of updates from 2012 onwards

Malaria still persists as one of the deadliest infectious disease having a huge morbidity and mortality affecting the higher population of the world. Structure and ligand-based drug design methods like molecular docking and MD simulations, pharmacophore modeling, QSAR and virtual screening are widely used to perceive the accordant correlation between the antimalarial activity and property of the compounds to design novel dominant and discriminant molecules. These modeling methods will speed-up antimalarial drug discovery, selection of better drug candidates for synthesis and to achieve potent and safer drugs. In this work, we have extensively reviewed the literature pertaining to the use and applications of various ligand and structure-based computational methods for the design of antimalarial agents. Different classes of molecules are discussed along with their target interactions pattern, which is responsible for antimalarial activity. Communicated by Ramaswamy H. Sarma.

Computational studies to identify the common type-I and type-II inhibitors against the CDK8 enzyme

In this study, multicomplex-based pharmacophore modeling was conducted on the structural proteome of the two states of CDK8 protein, that is, DMG-in and out. Three pharmacophores having six, five, and four features were selected as the representative models to conduct the virtual screening process using the prepared drug-like natural product database. The screened candidates were subjected to molecular docking studies on DMG-in (5XS2) and out (4F6U) conformation of the CDK8 protein. Subsequently, the common four docked candidates of 5XS2 and 4F6U were selected to perform the molecular dynamics simulation studies. Apart from one of the complexes of DMG-in (5XS2-UNPD163102), all other complexes displayed stable dynamic behavior. The interaction and stability studies of the docked complexes were compared with the references selected from the two conformations (DMG-in and out) of the protein. The current work leads to the identification of three common DMG-in and out hits with diverse scaffolds which can be employed as the initial leads for the design of the novel CDK8 inhibitors.

Application of ensemble pharmacophore-based virtual screening to the discovery of novel antimitotic tubulin inhibitors

Tubulin is a well-validated target for herbicides, fungicides, anti-parasitic, and anti-tumor drugs. Many of the non-cancer tubulin drugs bind to its colchicine site but no colchicine-site anticancer drug is available. The colchicine site is composed of three interconnected sub-pockets that fit their ligands and modify others' preference, making the design of molecular hybrids (that bind to more than one sub-pocket) a difficult task. Taking advantage of the more than eighty published X-ray structures of tubulin in complex with ligands bound to the colchicine site, we generated an ensemble of pharmacophore representations that flexibly sample the interactional space between the ligands and target. We searched the ZINC database for scaffolds able to fit several of the subpockets, such as tetrazoles, sulfonamides and diarylmethanes, selected roughly ~8000 compounds with favorable predicted properties. A Flexi-pharma virtual screening, based on ensemble pharmacophore, was performed by two different methodologies. Combining the scaffolds that best fit the ensemble pharmacophore-representation, we designed a new family of ligands, resulting in a novel tubulin modulator. We synthesized tetrazole 5 and tested it as a tubulin inhibitor in vitro. In good agreement with the design principles, it demonstrated micromolar activity against in vitro tubulin polymerization and nanomolar anti-proliferative effect against human epithelioid carcinoma HeLa cells through microtubule disruption, as shown by immunofluorescence confocal microscopy. The integrative methodology succedes in the design of new scaffolds for flexible proteins with structural coupling between pockets, thus expanding the way in which computational methods can be used as significant tools in the drug design process.

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 EC₅₀ = 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 (ΔG_binding=-57.337 kJ/mol).Communicated by Ramaswamy H. Sarma.

Identification of the natural compound inhibitors against Plasmodium falciparum plasmepsin-II via common feature based screening and molecular dynamics simulations

Malaria is counted amongst the deadly disease caused by Plasmodium falciparum. Recently, plasmepsin-II enzyme has gained much importance as an attractive drug target for the exploration of antimalarials. Therefore, the common feature pharmacophore models were generated from the crystallized complexes of the plasmepsin-II proteome. These models were subjected to a series of validation procedures, i.e. test set and Güner Henry studies to enlist the representative models. The selected representative hypotheses incorporating the most essential chemical features (common ZHHA) were screened against the natural product database to retrieve the potential candidates. To ensure the selection of the drug-like candidates, prior to screening, filtering steps (Drug-likeness and ADMET filters) were employed on the selected database. To study the interaction pattern of the candidates within the protein, these molecules were advanced to the molecular docking studies. Subsequently, based on the selected cut-off criteria obtained via redocking of the reference (4Z22), 15 compounds showed higher docking score (> -16.05 kcal/mol), and displayed the presence of hydrogen bonding with the crucial amino acids, i.e. Asp34 and Asp214. Further, the stability of the docked molecules was scrutinized via molecular dynamics simulations, and the results were compared with the reference compound 4Z22. All the docked compounds showed stable dynamics behaviour. Thus, in the present contribution, the combination of screening and stability procedures resulted in the identification of 15 hits that can serve as a new chemical space in the designing of the novel antimalarials.Communicated by Ramaswamy H. Sarma.