Structure and dynamics studies of sterol 24-C-methyltransferase with mechanism based inactivators for the disruption of ergosterol biosynthesis

Diosgenin biosynthesis pathway and its regulation in Dioscorea cirrhosa L

Dioscorea cirrhosa L. (D. cirrhosa) tuber is a traditional medicinal plant that is abundant in various pharmacological substances. Although diosgenin is commonly found in many Dioscoreaceae plants, its presence in D. cirrhosa remained uncertain. To address this, HPLC-MS/MS analysis was conducted and 13 diosgenin metabolites were identified in D. cirrhosa tuber. Furthermore, we utilized transcriptome data to identify 21 key enzymes and 43 unigenes that are involved in diosgenin biosynthesis, leading to a proposed pathway for diosgenin biosynthesis in D. cirrhosa. A total of 3,365 unigenes belonging to 82 transcription factor (TF) families were annotated, including MYB, AP2/ERF, bZIP, bHLH, WRKY, NAC, C2H2, C3H, SNF2 and Aux/IAA. Correlation analysis revealed that 22 TFs are strongly associated with diosgenin biosynthesis genes (-r2- > 0.9, P < 0.05). Moreover, our analysis of the CYP450 gene family identified 206 CYP450 genes (CYP450s), with 40 being potential CYP450s. Gene phylogenetic analysis revealed that these CYP450s were associated with sterol C-22 hydroxylase, sterol-14-demethylase and amyrin oxidase in diosgenin biosynthesis. Our findings lay a foundation for future genetic engineering studies aimed at improving the biosynthesis of diosgenin compounds in plants.

Novel chemical scaffold as potential drug against Leishmania donovani: Integrated computational and experimental approaches

In this study, we have screened a large number of Food and Drug Administration-approved compounds for novel anti-leishmanial molecules targeting the citrate synthase enzyme of the parasite. Based on their docking and molecular dynamic simulation statistics, five compounds were selected. These compounds followed Lipinski's rule of five. Additionally, in vitro, antileishmanial and cytotoxicity studies were performed. The three compounds, Abemaciclib, Bazedoxifene, and Vorapaxar, had shown effective anti-leishmanial activities with IC50 values of 0.92 ± 0.02, 0.65 ± 0.09, and 6.1 ± 0.91 against Leishmania donovani promastigote and with EC50 values of 1.52 ± 0.37, 2.11 ± 0.38, 10.4 ± 1.27 against intramacrophagic amastigote without significantly harming macrophage cells. Among them, from in silico and antileishmanial activities studies, Abemaciclib had been selected based on their less binding energy, good antileishmanial activities, and also a significant difference in their binding energy with human citrate synthase for cell death mechanistic studies using flow cytometry and a DNA fragmentation assay. The action of this compound resulted in an increased reactive oxygen species production, depolarization of mitochondrial membrane potential, DNA damage, and an increase in the sub-G1 cell population. These properties are the hallmarks of apoptosis which were further confirmed by apoptotic assay. Based on the above result, this anticancer compound Abemaciclib could be employed as a potential treatment option for leishmaniasis after further confirmation.

Inhibiting Leishmania donovani Sterol Methyltransferase to Identify Lead Compounds Using Molecular Modelling

The recent outlook of leishmaniasis as a global public health concern coupled with the reportage of resistance and lack of efficacy of most antileishmanial drugs calls for a concerted effort to find new leads. The study combined In silico and in vitro approaches to identify novel potential synthetic small-molecule inhibitors targeting the Leishmania donovani sterol methyltransferase (LdSMT). The LdSMT enzyme in the ergosterol biosynthetic pathway is required for the parasite’s membrane fluidity, distribution of membrane proteins, and control of the cell cycle. The lack of LdSMT homologue in the human host and its conserved nature among all Leishmania parasites makes it a viable target for future antileishmanial drugs. Initially, six known inhibitors of LdSMT with IC50 < 10 μM were used to generate a pharmacophore model with a score of 0.9144 using LigandScout. The validated model was used to screen a synthetic library of 95,630 compounds obtained from InterBioScreen limited. Twenty compounds with pharmacophore fit scores above 50 were docked against the modelled three-dimensional structure of LdSMT using AutoDock Vina. Consequently, nine compounds with binding energies ranging from −7.5 to −8.7 kcal/mol were identified as potential hit molecules. Three compounds comprising STOCK6S-06707, STOCK6S-84928, and STOCK6S-65920 with respective binding energies of −8.7, −8.2, and −8.0 kcal/mol, lower than 22,26-azasterol (−7.6 kcal/mol), a known LdSMT inhibitor, were selected as plausible lead molecules. Molecular dynamics simulation studies and molecular mechanics Poisson–Boltzmann surface area calculations showed that the residues Asp25 and Trp208 were critical for ligand binding. The compounds were also predicted to have antileishmanial activity with reasonable pharmacological and toxicity profiles. When the antileishmanial activity of the three hits was evaluated in vitro against the promastigotes of L. donovani, mean half-maximal inhibitory concentrations (IC50) of 21.9 ± 1.5 μM (STOCK6S-06707), 23.5 ± 1.1 μM (STOCK6S-84928), and 118.3 ± 5.8 μM (STOCK6S-65920) were obtained. Furthermore, STOCK6S-84928 and STOCK6S-65920 inhibited the growth of Trypanosoma brucei, with IC50 of 14.3 ± 2.0 μM and 18.1 ± 1.4 μM, respectively. The identified compounds could be optimised to develop potent antileishmanial therapeutic agents.

Evaluation of the inhibitory potential of Valproic acid against histone deacetylase of Leishmania donovani and computational studies of Valproic acid derivatives

Valproic acid (VA) is a proven inhibitor of human histone deacetylases (HDACs). The homogenous HDAC has been associated with all major human parasitic pathogens and hence, it has been considered an attractive drug target for anti-leishmanial therapy. To assist in drug design endeavors for HDACs, an in-vitro study has been presented to investigate the VA inhibition on Leishmania donovani HDAC (LdHDAC). The regression analysis of VA by 24 hrs viability assay confirmed its activity against LdHDAC. Moreover, the toxicity of VA is also well documented. Thus, the in-silico experiments were also conducted to screen the non-toxic VA derivatives as anti-leishmanial drug candidates having potential as inhibitors of LdHDAC. For in-silico study, the 3D structure of target LdHDAC was developed by homology modeling. Based on their in-silico activity, we shortlisted 13 VA derivatives having maximum affinity for LdHDAC and identified four potential derivatives that can specifically bind to this protein. After that, these ligands were subjected to molecular dynamics simulation. These derivatives may be effective against L. donovani promastigotes since they followed Lipinski's RO5 and were non-toxic. Thus, screened derivatives can be considered as lead ligands for targeting LdHDAC and may be used as possible drug candidates to treat leishmaniasis and overcome the limitation of anti-leishmanial drugs. This is the first report of antileishmanial potential of VA and its derivatives targeting LdHDAC. Hence, the current investigation presents a search for novel target specific drugs to aid the anti-leishmanial drug development. Communicated by Ramaswamy H. Sarma.

Homology Modeling, de Novo Design of Ligands, and Molecular Docking Identify Potential Inhibitors of Leishmania donovani 24-Sterol Methyltransferase

The therapeutic challenges pertaining to leishmaniasis due to reported chemoresistance and toxicity necessitate the need to explore novel pathways to identify plausible inhibitory molecules. Leishmania donovani 24-sterol methyltransferase (LdSMT) is vital for the synthesis of ergosterols, the main constituents of Leishmania cellular membranes. So far, mammals have not been shown to possess SMT or ergosterols, making the pathway a prime candidate for drug discovery. The structural model of LdSMT was elucidated using homology modeling to identify potential novel 24-SMT inhibitors via virtual screening, scaffold hopping, and de-novo fragment-based design. Altogether, six potential novel inhibitors were identified with binding energies ranging from −7.0 to −8.4 kcal/mol with e-LEA3D using 22,26-azasterol and S1–S4 obtained from scaffold hopping via the ChEMBL, DrugBank, PubChem, ChemSpider, and ZINC15 databases. These ligands showed comparable binding energy to 22,26-azasterol (−7.6 kcal/mol), the main inhibitor of LdSMT. Moreover, all the compounds had plausible ligand efficiency-dependent lipophilicity (LELP) scores above 3. The binding mechanism identified Tyr92 to be critical for binding, and this was corroborated via molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The ligand A1 was predicted to possess antileishmanial properties with a probability of activity (Pa) of 0.362 and a probability of inactivity (Pi) of 0.066, while A5 and A6 possessed dermatological properties with Pa values of 0.205 and 0.249 and Pi values of 0.162 and 0.120, respectively. Structural similarity search via DrugBank identified vabicaserin, daledalin, zanapezil, imipramine, and cefradine with antileishmanial properties suggesting that the de-novo compounds could be explored as potential antileishmanial agents.

Computational study to discover potent phytochemical inhibitors against drug target, squalene synthase from Leishmania donovani

AIMS

The parasite, Leishmania donovani is responsible for lethal visceral leishmaniasis (VL) in humans. There is a need to investigate novel medicines as antileishmanial drugs, as medication currently introduced for leishmaniasis may cause resistance, serious side-effects, chemical instability and high cost. Therefore, this computational study was designed to explore potential phytochemical inhibitors against Leishmania donovani squalene synthase (LdSQS) enzyme, a drug target.

MAIN METHODS

Multiple sequence alignment was carried to detect conserved regions across squalene synthases from different Leishmania spp. Their evolutionary relationships were studied by generating phylogenetic tree. Homology modeling method was used to build a three dimensional model of the protein. The validated model was explored by docking simulation with the phytochemicals of interest to identify the most potent inhibitors. Two reported inhibitors were used as references in the virtual screening. The top hit compounds (binding energy less than -9 kcal/mol) were further subjected to intermolecular interaction analysis, pharmacophore modeling, pharmacokinetic and toxicity prediction.

KEY FINDINGS

Seven phytochemicals displayed binding energies less than -9 kcal/mol hence demonstrating ability to be strongly bound to the active site of LdSQS to inhibit the enzymatic activity. Ancistrotanzanine B demonstrated the lowest binding affinity of -9.83 kcal/mol superior to reported inhibitors in literature. Conserved two aspartate rich regions and two signatory motifs were found in the L. donovani squalene synthase by multiple sequence alignment. In addition, study of pharmacophore modeling confirmed that top hit phytochemicals and the reported inhibitor (E5700) share common chemical features for their biochemical interaction with LdSQS. Among seven phytochemicals, 3-O-methyldiplacol showed admissible physicochemical, pharmacokinetic and toxicity predictions compared to the reported inhibitors. All seven phytochemicals satisfied in silico prediction criteria for oral bioavailability.

SIGNIFICANCE

Based on the current study, these hits can be further structurally optimized and validated under laboratory conditions to develop antileishmanial drugs.

Repurposing of FDA-approved drugs as inhibitors of sterol C-24 methyltransferase of Leishmania donovani to fight against leishmaniasis

Leishmaniasis is a vector-borne disease caused by around 20 species of Leishmania. The main clinical forms of leishmaniasis are cutaneous leishmaniasis (CL) and visceral leishmaniasis (VL). VL is caused by Leishmania infantum in Central and South America, Mediterranean Basin, Middle East, and by L. donovani in Asia and Africa. Sterol C-24 methyltransferase (LdSMT) of L. donovani is a transferase enzyme of the sterol biosynthesis pathway. This pathway is one of the major targets for drug developments in Leishmania. Due to insufficient evidence about the exact function of SMT inside the cell and the uniqueness of the SMT enzyme in the Leishmania parasites made it a significant target for an effective drug development approach. We performed virtual screening of the Food and Drug Administration (FDA)-approved drug library against LdSMT and found simeprevir, an antiviral drug on top in the binding score. It showed a significant binding affinity with LdSMT. The binding was supported by hydrogen bonds and several other interactions. Simeprevir inhibited L. donovani growth of promastigotes with 50% inhibitory concentration (IC₅₀ ) of 51.49 ± 5.87 μM. Further studies showed that simeprevir induced ROS generation in 44.7% of parasites at 125-μM concentration. Here, we for the first time reported simeprevir as an antileishmanial lead molecule using a drug repurposing approach.

Screening Pipeline for Flavivirus Based Inhibitors for Zika Virus NS1

In-silico pipeline is applied for identifying and designing novel inhibitors against ZIKV NS1 protein. Comparative molecular docking studies are performed to explore the binding of structurally diverse compounds to ZIKV NS1 by AutoDock/Vina and GOLD. The Zika virus (ZIKV) is a flavivirus, responsible for life-threatening infections and transmitted by Aedes mosquitoes in other organisms. It is associated with Guillain Barre Syndrome (GBS) and microcephaly. This epidemic increase in GBS and microcephaly convoyed the World Health Organization to affirm ZIKV a public health crisis. To combat the ZIKV infections, non-structural protein 1 (NS1), a major host-interaction molecule contributing towards replication, pathogenesis and immune evasion is targeted in the current study. For this purpose, a comprehensive study is required to develop potential novel antiviral inhibitors. Three compounds were identified through docking programs exhibiting properties which are non-toxic to human host and could inhibit the elusive ZIKV. Significant interaction with active site residues and H-bond interactions with the key residues were analyzed for these compounds using molecular dynamics simulation. Free energy calculation predicted higher affinity of Deoxycalyxin-A for ZIKV NS1. This study contributes towards fighting ZIKV infections and can help researchers in designing drug for the treatment of ZIKV.

Deciphering the Inhibition Mechanism of under Trial Hsp90 Inhibitors and Their Analogues: A Comparative Molecular Dynamics Simulation

Heat shock protein 90 (Hsp90) performs functions in cellular activities together with other signaling pathways. Hsp90 is evolutionarily conserved and universally articulated as a human cancer-causing agent involved in lung cancer and breast cancer followed by colon and rectum cancers. It has emerged as an effective drug candidate, and inhibition may affect several signaling pathways associated with cancer spread. Therefore, in-silico approaches, molecular docking, molecular dynamics simulation, and binding free energy calculations were applied to create insights into the inhibition mechanism against Hsp90 to identify new cancer therapeutic drugs. Top-docked Hsp90-inhibitor complexes with their analogues were selected as the best complexes based on the GOLD fitness score and orientation. The significant interaction of Hsp90 inhibitors and their analogues were observed to be bound with active site residues as well as residing within the same cavity region. System stability factors RMSD, RMSF, beta-factor, and radius of gyration were analyzed for top-docked complexes and ensure strong binding interaction between inhibitors and the Hsp90 cavity. Cavity bound inhibitors were found to retain consistent hydrogen bonding during the simulation. The radial distribution function (RDF) illustrated that interacting active site residues drive the binding and stability of the inhibitors. Similarly, the axial frequency distribution, which is an indigenously developed analytical tool, produced noteworthy knowledge of the hydrogen-bonding pattern. Results yielded new insights into the design of cancer therapeutic drugs against Hsp90. This finding suggests that under trial Hsp90 inhibitors MPC-3100 could be a potential starting point into the development of potential anticancer agents with the possibility of future directions for the improvement of early existing Hsp90 inhibitors CNF-2024 and SNX-5422 as an anticancer agent.

Bioactive compounds induced in Physalis angulata L. by methyl-jasmonate: an investigation of compound accumulation patterns and biosynthesis-related candidate genes

We employed both metabolomic and transcriptomic approaches to explore the accumulation patterns of physalins, flavonoids and chlorogenic acid in Physalis angulata and revealed the genes associated with the biosynthesis of bioactive compounds under methyl-jasmonate (MeJA) treatment. Physalis angulata L. is an annual Solanaceae plant with a number of medicinally active compounds. Despite the potential pharmacological benefits of P. angulata, the scarce genomic information regarding this plant has limited the studies on the mechanisms of bioactive compound biosynthesis. To facilitate the basic understanding of the main chemical constituent biosynthesis pathways, we performed both metabolomic and transcriptomic approaches to reveal the genes associated with the biosynthesis of bioactive compounds under methyl-jasmonate (MeJA) treatment. Untargeted metabolome analysis showed that most physalins, flavonoids and chlorogenic acid were significantly upregulated. Targeted HPLC-MS/MS analysis confirmed variations in the contents of two important representative steroid derivatives (physalins B and G), total flavonoids, neochlorogenic acid, and chlorogenic acid between MeJA-treated plants and controls. Transcript levels of a few steroid biosynthesis-, flavonoid biosynthesis-, and chlorogenic acid biosynthesis-related genes were upregulated, providing a potential explanation for MeJA-induced active ingredient synthesis in P. angulata. Systematic correlation analysis identified a number of novel candidate genes associated with bioactive compound biosynthesis. These results may help to elucidate the regulatory mechanism underlying MeJA-induced active compound accumulation and provide several valuable candidate genes for further functional study.

Dual RNA-Seq Analysis of Trichophyton rubrum and HaCat Keratinocyte Co-Culture Highlights Important Genes for Fungal-Host Interaction

The dermatophyte Trichophyton rubrum is the major fungal pathogen of skin, hair, and nails that uses keratinized substrates as the primary nutrients during infection. Few strategies are available that permit a better understanding of the molecular mechanisms involved in the interaction of T. rubrum with the host because of the limitations of models mimicking this interaction. Dual RNA-seq is a powerful tool to unravel this complex interaction since it enables simultaneous evaluation of the transcriptome of two organisms. Using this technology in an in vitro model of co-culture, this study evaluated the transcriptional profile of genes involved in fungus-host interactions in 24 h. Our data demonstrated the induction of glyoxylate cycle genes, ERG6 and TERG_00916, which encodes a carboxylic acid transporter that may improve the assimilation of nutrients and fungal survival in the host. Furthermore, genes encoding keratinolytic proteases were also induced. In human keratinocytes (HaCat) cells, the SLC11A1, RNASE7, and CSF2 genes were induced and the products of these genes are known to have antimicrobial activity. In addition, the FLG and KRT1 genes involved in the epithelial barrier integrity were inhibited. This analysis showed the modulation of important genes involved in T. rubrum–host interaction, which could represent potential antifungal targets for the treatment of dermatophytoses.

Subtractive Genomics, Molecular Docking and Molecular Dynamics Simulation Revealed LpxC as a Potential Drug Target Against Multi-Drug Resistant Klebsiella pneumoniae

The emergence and dissemination of pan drug resistant clones of Klebsiella pneumoniae are great threat to public health. In this regard new therapeutic targets must be highlighted to pave the path for novel drug discovery and development. Subtractive proteomic pipeline brought forth UDP-3-O-[3-hydroxymyristoyl] N-acetylglucosamine deacetylase (LpxC), a Zn⁺² dependent cytoplasmic metalloprotein and catalyze the rate limiting deacetylation step of lipid A biosynthesis pathway. Primary sequence analysis followed by 3-dimensional (3-D) structure elucidation of the protein led to the detection of K. pneumoniae LpxC (KpLpxC) topology distinct from its orthologous counterparts in other bacterial species. Molecular docking study of the protein recognized receptor antagonist compound 106, a uridine-based LpxC inhibitory compound, as a ligand best able to fit the binding pocket with a Gold Score of 67.53. Molecular dynamics simulation of docked KpLpxC revealed an alternate binding pattern of ligand in the active site. The ligand tail exhibited preferred binding to the domain I residues as opposed to the substrate binding hydrophobic channel of subdomain II, usually targeted by inhibitory compounds. Comparison with the undocked KpLpxC system demonstrated ligand induced high conformational changes in the hydrophobic channel of subdomain II in KpLpxC. Hence, ligand exerted its inhibitory potential by rendering the channel unstable for substrate binding.

In silico Tools for Target Identification and Drug Molecular Docking in Leishmania

Neglected tropical diseases represent a significant health burden in large parts of the world. Drug discovery is currently a key bottleneck in the pipeline of these diseases. In this chapter, the in silico approaches used for the processes involved in drug discovery, identification and validation of druggable Leishmania targets, and design and optimisation of new anti-leishmanial drugs are discussed. We also provide a general view of the different computational tools that can be employed in pursuit of this aim, along with the most interesting cases found in the literature.