L . donovani XPRT: Molecular characterization and evaluation of inhibitors

Identification of putative volatile biomarkers of canine leishmaniasis in dog's breath and hair employing a novel algorithm for automated chromatographic peak detection and matching

The analysis of the volatile compounds released by biological samples represents a promising approach for the non-invasive diagnosis of a disease. The present study, focused on a population of dogs infected with canine leishmaniasis, aimed to decipher the volatolomic profile associated with this disease in dogs, which represent the main animal reservoir for Leishmania pathogen transmission to humans. The volatiles emitted by the breath and hair of dogs were analysed employing the gas chromatography-mass spectrometry (GC-MS) technique. The acquired chromatograms were investigated using a novel algorithm developed in this study for automated chromatographic peak detection and matching in untargeted GC-MS analysis, which includes various steps that comprise noise reduction, m/z filtering, background subtraction, peak detection, peak matching, and generation of a peak table for compounds identification. The results revealed one tentative breath volatile biomarker and five tentative hair volatile biomarkers for the cutaneous form of the disease, which is characterised by skin ulcerations. Additionally, nine tentative breath volatile biomarkers and four tentative hair volatile biomarkers were found for the visceral form of the disease, which affects internal organs such as spleen, liver and bone marrow. All tentative biomarkers identified in this study were upregulated in cutaneous leishmaniasis, while in visceral leishmaniasis, all tentative biomarkers were upregulated in the breath and only one out of four in the hair. Only one compound (glyceryl monooleate) was identified as tentative volatile biomarker for both forms of the disease, in the hair of dogs.

Leishmania donovani adenylosuccinate synthetase requires IMP for dimerization and organization of the active site

Adenylosuccinate synthetase (AdSS), which catalyses the GTP-dependent conversion of inosine monophosphate (IMP) and aspartic acid to succinyl-AMP, plays a major role in purine biosynthesis. In some bacterial AdSS, it is implicated that IMP binding is important to organize the active site, but in certain plant AdSS, GTP performs this role. Here, we report that in Leishmania donovani AdSS, IMP binding favoured dimerization, induced greater conformational change and improved the protein stability more than GTP binding. IMP binding, which resulted in a network of hydrogen bonds, stabilized the conformation of active site loops and brought the switch loop to a closed conformation, which then facilitated GTP binding. Our results provide a basis for designing better inhibitors of leishmanial AdSS.

Insights into the ATP / GTP selectivity of a GTPase, adenylosuccinate synthetase from Leishmania donovani

Many GTPases have been shown to utilize ATP too as the phosphoryl donor. Both GTP and ATP are important molecules in the cellular environments and play multiple and discrete functional role within the cells. In our present study, we showed that one of the purine metabolic enzymes Adenylosuccinate synthetase from Leishmania donovani (LdAdSS) which belongs to the BioD-superfamily of GTPases can also carry out the catalysis by hydrolysing ATP instead of its cognate substrate GTP albeit with less efficiency. Biochemical and biophysical studies indicated its ability to bind to ATP too but at a higher concentration of ATP compared to that of GTP. Sequence analysis and molecular dynamic simulations suggested that residues of the switch loop and the G4-G5 (593SAXD596) connected motif of LdAdSS plays a role in determining the nucleotide specificity. Though the crucial interaction between Asp596 and the nucleotide is broken when ATP is bound, interactions between the Ala594 and the adenine ring of ATP could still hold ATP in the GTP binding site. The results of the present study suggested that though LdAdSS is GTP specific, it still shows ATP hydrolysing activity.

An in-silico approach to unravel the structure of 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase (DAHPS): a critical enzyme for sennoside biosynthesis in Cassia angustifolia Vahl

The laxative properties of senna are attributed to the presence of sennosides produced in the plant. The low production level of sennosides in the plant is an important impediment to their growing demand and utilization. Understanding biosynthetic pathways helps to engineer them in terms of enhanced production. The biosynthetic pathways of sennoside production in plants are not completely known yet. However, attempts to get information on genes and proteins engaged in it have been made which decode involvement of various pathways including shikimate pathway. 3-deoxy-D-arabino-heptulosonate 7-phosphate synthase (DAHPS) is a key enzyme involved in sennosides production through the shikimate pathway. Unfortunately, there is no information available on proteomic characterization of DAHPS enzyme of senna (caDAHPS) resulting in lack of knowledge about its role. We for the first time characterized DAHPS enzyme of senna using in-silico analysis. To the best of our knowledge this is the first attempt to identify the coding sequence of caDAHPS by cloning and sequencing. We found Gln179, Arg175, Glu462, Glu302, Lys357 and His420 amino acids in the active site of caDAHPS through molecular docking. followed by molecular dynamic simulation. The amino acid residues, Lys182, Cys136, His460, Leu304, Gly333, Glu334, Pro183, Asp492 and Arg433 at the surface interact with PEP by van der Waals bonds imparting stability to the enzyme-substrate complex. Docking results were further validated by molecular dynamics. The presented in-silico analysis of caDAHPS will generate opportunities to engineer the sennoside biosynthesis in plants.Communicated by Ramaswamy H. Sarma.

The Application of MD Simulation to Lead Identification, Vaccine Design, and Structural Studies in Combat against Leishmaniasis - A Review

Drug discovery, vaccine design, and protein interaction studies are rapidly moving toward the routine use of molecular dynamics simulations (MDS) and related methods. As a result of MDS, it is possible to gain insights into the dynamics and function of identified drug targets, antibody-antigen interactions, potential vaccine candidates, intrinsically disordered proteins, and essential proteins. The MDS appears to be used in all possible ways in combating diseases such as cancer, however, it has not been well documented as to how effectively it is applied to infectious diseases such as Leishmaniasis. As a result, this review aims to survey the application of MDS in combating leishmaniasis. We have systematically collected articles that illustrate the implementation of MDS in drug discovery, vaccine development, and structural studies related to Leishmaniasis. Of all the articles reviewed, we identified that only a limited number of studies focused on the development of vaccines against Leishmaniasis through MDS. Also, the PCA and FEL studies were not carried out in most of the studies. These two were globally accepted utilities to understand the conformational changes and hence it is recommended that this analysis should be taken up in similar approaches in the future.

Rohitukine content across the geographical distribution of Dysoxylum binectariferum Hook F. and its natural derivatives as potential sources of CDK inhibitors

Dysoxylum binectariferum is an important medicinal plant distributed in the Western Ghats of India. The species has gained international importance for its anticancer component, rohitukine, a chromone alkaloid. Flavopiridol, P-276-00 and IIIM-290 are the derivatives of rohitukine in clinical trials against a wide range of cancers. Flavopiridol was recently approved as an orphan drug for chronic lymphocytic leukemia treatment. In this study, we report the isolation and characterization of rohitukine from the bark of D. binectariferum. Further, rohitukine was estimated across the Western-Ghats and the North-East regions of India. Additionally, D. binectariferum is also reported (∼45 compounds) to produce many natural derivatives of rohitukine and terpenoids, which were investigated in-silico to reveal promising CDK inhibitors. The metabolite fingerprinting of tissues of D. binectariferum was studied using HPTLC and FTIR. The distribution of major chromone alkaloid rohitukine was estimated by HPLC. Further, the pharmacological potential of D. binectariferum compounds was evaluated in-silico by discovering the potential protein targets, molecular docking, ADMET analysis and MD simulation. The isolation of rohitukine has yielded 0.6% from the bark of D. binectariferum. A higher percent of rohitukine was found in the Jog populations (0.58% & 1.28%: leaf & bark), whereas least was observed in the Phasighat population (∼0.06%: both leaf & bark). Across the geographic regions, a higher percent of rohitukine was found in the Central-southern Western Ghats, whereas lower in the northern parts of the Western Ghats and Northeast regions. The leaves produce a considerably higher percent of rohitukine and could be used as a sustainable source of rohitukine. The rohitukine analogues, along with other chromone alkaloids of D. binecatariferum were found to be more interactive with the "kinases" family of proteins, majorly "Serine/threonine-protein kinase PFTAIRE-2" (CDK15) with high confidence level (0.94-0.98). The molecular docking of these chromone alkaloids found a strong binding energy with six CDKs (-3.1 to -10.6 kcal/mol) along with a promising ADMET profile. In addition, molecular dynamic simulation found that the rohitukine complexes are virtually constant with CDK-1, 2, 9 and 15, which is substantiated with MM-PBSA free energy calculations. The chromone alkaloids, majorly rohitukine and its analogues were closely clustered with flavopiridol, P-276-00 and IIIM-290 along with other chrotacumines in the chemical phylogeny. In conclusion, D. binectariferum is a rich source of chromone alkaloids, which could lead to the discovery of more potential scaffolding for CDK inhibitors as anticancer drugs.

Structural differences in 3C-like protease (Mpro) from SARS-CoV and SARS-CoV-2: molecular insights revealed by Molecular Dynamics Simulations

Novel coronavirus SARS-CoV-2 has infected millions of people with thousands of mortalities globally. The main protease (Mpro) is vital in processing replicase polyproteins. Both the CoV's Mpro shares 97% identity, with 12 mutations, but none are present in the active site. Although many therapeutics and vaccines are available to combat SARS-CoV-2, these treatments may not be practical due to their high mutational rate. On the other hand, Mpro has a high degree of conservation throughout variants, making Mpro a stout drug target. Here, we report a detailed comparison of both the monomeric Mpro and the biologically active dimeric Mpro using MD simulation to understand the impact of the 12 divergent residues (T35V, A46S, S65N, L86V, R88K, S94A, H134F, K180N, L202V, A267S, T285A and I286L) on the molecular microenvironment and the interaction between crucial residues. The present study concluded that the change in the microenvironment of residues at the entrance (T25, T26, M49 and Q189), near the catalytic site (F140, H163, H164, M165 and H172) and in the substrate-binding site (V35, N65, K88 and N180) is due to 12 mutations in the SARS-CoV-2 Mpro. Furthermore, the involvement of F140, E166 and H172 residues in dimerization stabilizes the Mpro dimer, which should be considered. We anticipate that networks and microenvironment changes identified here might guide repurposing attempts and optimization of new Mpro inhibitors.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02089-6.

Designing multi-epitope based peptide vaccine targeting spike protein SARS-CoV-2 B1.1.529 (Omicron) variant using computational approaches

Millions of lives have been infected since the SARS-CoV-2 outbreak in 2019. The high human-to-human transmission rate has warranted a need for a vaccine to protect people. Although some vaccines are in use, due to the high mutation rate in the SARS-CoV-2 multiple variants, the current vaccines may not be sufficient to immunize people against new variant threats. One of the emerging concern variants is B1.1.529 (Omicron), which carries ~ 30 mutations in the Spike protein (S) of SARS-CoV-2 and is predicted to evade antibody recognition even from vaccinated people. We used a structure-based approach and an epitope prediction server to develop a Multi-Epitope based Subunit Vaccine (MESV) involving SARS-CoV-2 B1.1.529 variant spike glycoprotein. The predicted epitope with better antigenicity and non-toxicity was used for designing and predicting vaccine construct features and structure models. In addition, the MESV construct In silico cloning in the pET28a expression vector predicted the construct to be highly translational. The proposed MESV vaccine construct was also subjected to immune simulation prediction and was found to be highly antigenic and elicit a cell-mediated immune response. Therefore, the proposed MESV in the present study has the potential to be evaluated further for vaccine production against the newly identified B1.1.529 (Omicron) variant of concern.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02027-6.

Development of Novel Anti-Leishmanials: The Case for Structure-Based Approaches

The neglected tropical disease (NTD) leishmaniasis is the collective name given to a diverse group of illnesses caused by ~20 species belonging to the genus Leishmania, a majority of which are vector borne and associated with complex life cycles that cause immense health, social, and economic burdens locally, but individually are not a major global health priority. Therapeutic approaches against leishmaniasis have various inadequacies including drug resistance and a lack of effective control and eradication of the disease spread. Therefore, the development of a rationale-driven, target based approaches towards novel therapeutics against leishmaniasis is an emergent need. The utilization of Artificial Intelligence/Machine Learning methods, which have made significant advances in drug discovery applications, would benefit the discovery process. In this review, following a summary of the disease epidemiology and available therapies, we consider three important leishmanial metabolic pathways that can be attractive targets for a structure-based drug discovery approach towards the development of novel anti-leishmanials. The folate biosynthesis pathway is critical, as Leishmania is auxotrophic for folates that are essential in many metabolic pathways. Leishmania can not synthesize purines de novo, and salvage them from the host, making the purine salvage pathway an attractive target for novel therapeutics. Leishmania also possesses an organelle glycosome, evolutionarily related to peroxisomes of higher eukaryotes, which is essential for the survival of the parasite. Research towards therapeutics is underway against enzymes from the first two pathways, while the third is as yet unexplored.

Identification of natural inhibitors against prime targets of SARS-CoV-2 using molecular docking, molecular dynamics simulation and MM-PBSA approaches

The recently emerged COVID-19 has been declared a pandemic by the World Health Organization as to date; no therapeutic drug/vaccine is available for the treatment. Due to the lack of time and the urgency to contain the pandemic, computational screening appears to be the best tool to find a therapeutic solution. Accumulated evidence suggests that many phyto-compounds possess anti-viral activity. Therefore, we identified possible phyto-compounds that could be developed and used for COVID-19 treatment. In particular, molecular docking was used to prioritize the possible active phyto-compounds against two key targets namely RNA dependent RNA polymerase (RdRp) and main protease (Mpro) of SARS-CoV-2. In this study, an antiviral drug- Remdesivir (RdRp inhibitor) and Darunavir (Mpro inhibitor) are used as reference drugs. This study revealed that phyto-molecules- Mulberroside-A/C/E/F, Emblicanin A, Nimbolide, and Punigluconin showed high binding affinity against RdRp while Andrographolides, Mulberrosides, Anolignans, Chebulic acid, Mimusopic acid, and Punigluconin showed better binding affinity against Mpro as compared with the reference drug. Furthermore, ADME profiles validated the drug-likeness properties of prioritized phyto-compounds. Besides, to assess the stability, MD simulations studies were performed along with reference inhibitors for Mpro (Darunavir) and RdRp (Remdesivir). Binding free energy calculations (MM-PBSA) revealed the estimated value (ΔG) of Mpro_Darunavir; Mpro_Mulberroside E; RdRp_Remdesivir and RdRp_Emblicanin A were -111.62 ± 6.788, -141.443 ± 9.313, 30.782 ± 5.85 and -89.424 ± 3.130 kJmol-1, respectively. Taken together, the study revealed the potential of these phyto-compounds as inhibitors of RdRp and Mpro inhibitor that could be further validated against SARS-CoV-2 for clinical benefits.Communicated by Ramaswamy H. Sarma.

Ile209 of Leishmania donovani xanthine phosphoribosyltransferase plays a key role in determining its purine base specificity

Xanthine phosphoribosyltransferase (XPRT) and hypoxanthine-guanine phosphoribosyltransferase (HGPRT) are purine salvaging enzymes of Leishmania donovani with distinct 6-oxopurine specificities. LdXPRT phosphoribosylates xanthine, hypoxanthine, and guanine, with preference toward xanthine, whereas LdHGPRT phosphoribosylates only hypoxanthine and guanine. In our study, LdXPRT was used as a model to understand these purine base specificities. Mutating I209 to V, the conserved residue found in HGPRTs, reduced the affinity of LdXPRT for xanthine, converting it to an HGXPRT-like enzyme. The Y208F mutation in the active site indicated that aromatic residue interactions with the purine ring are limited to pi-pi binding forces and do not impart purine base specificity. Deleting the unique motif (L55-Y82) of LdXPRT affected enzyme activity. Our studies established I209 as a key residue determining the 6-oxopurine specificity of LdXPRT.

Computational investigation of binding of chloroquinone and hydroxychloroquinone against PLPro of SARS-CoV-2

Novel coronavirus SARS-CoV-2 has infected 18 million people with 700,000+ mortalities worldwide and this deadly numeric figure is rapidly rising. With very few success stories, the therapeutic targeting of this epidemic has been mainly attributed to main protease (Mpro), whilst Papain-like proteases (PLpro) also plays a vital role in the processing of replicase polyprotein. Multifunctional roles of PLpro such as viral polypeptide cleavage, de-ISGlyation and immune suppression have made it a promising drug target for therapeutic interventions. Whilst there have been a number of studies and others are on-going on repurposing and new-small molecule screening, albeit previously FDA approved drugs viz. Chloroquine (CQ) and Hydroxychloroquine (HCQ) have only been found effective against this pandemic. Inspired by this fact, we have carried out molecular docking and dynamics simulation studies of FDA approved CQ and HCQ against SARS-CoV-2 PLpro. The end aim is to characterise the binding mode of CQ and HCQ and identify the key amino acid residues involved in the mechanism of action. Further, molecular dynamics simulations (MDS) were carried out with the docked complex to search for the conformational space and for understanding the integrity of binding mode. We showed that the CQ and HCQ can bind with better binding affinity with PLpro as compared to reference known PLpro inhibitor. Based on the presented findings, it can be anticipated that the SARS-CoV-2 PLpro may act as molecular target of CQ and HCQ, and can be projected for further exploration to design potent inhibitors of SARS-CoV-2 PLpro in the near future.

Identification of natural compound inhibitors against PfDXR: A hybrid structure-based molecular modeling approach and molecular dynamics simulation studies

In the present contribution, multicomplex-based pharmacophore studies were carried out on the structural proteome of Plasmodium falciparum 1-deoxy-D-xylulose-5-phosphate reductoisomerase. Among the constructed models, a representative model with complementary features, accountable for the inhibition was used as a primary filter for the screening of database molecules. Auxiliary evaluations of the screened molecules were performed via drug-likeness and molecular docking studies. Subsequently, the stability of the docked inhibitors was envisioned by molecular dynamics simulations, principle component analysis, and molecular mechanics-Poisson-Boltzmann surface area-based free binding energy calculations. The stability assessment of the hits was done by comparing with the reference (beta-substituted fosmidomycin analog, LC5) to prioritize more potent candidates. All the complexes showed stable dynamic behavior while three of them displayed higher binding free energy compared with the reference. The work resulted in the identification of the compounds with diverse scaffolds, which could be used as initial leads for the design of novel PfDXR inhibitors.