Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Exploring optimal drug targets through subtractive proteomics analysis and pangenomic insights for tailored drug design in tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis, ranks among the top causes of global human mortality, as reported by the World Health Organization's 2022 TB report. The prevalence of M. tuberculosis strains that are multiple and extensive-drug resistant represents a significant barrier to TB eradication. Fortunately, having many completely sequenced M. tuberculosis genomes available has made it possible to investigate the species pangenome, conduct a pan-phylogenetic investigation, and find potential new drug targets. The 442 complete genome dataset was used to estimate the pangenome of M. tuberculosis. This study involved phylogenomic classification and in-depth analyses. Sequential filters were applied to the conserved core genome containing 2754 proteins. These filters assessed non-human homology, virulence, essentiality, physiochemical properties, and pathway analysis. Through these intensive filtering approaches, promising broad-spectrum therapeutic targets were identified. These targets were docked with FDA-approved compounds readily available on the ZINC database. Selected highly ranked ligands with inhibitory potential include dihydroergotamine and abiraterone acetate. The effectiveness of the ligands has been supported by molecular dynamics simulation of the ligand-protein complexes, instilling optimism that the identified lead compounds may serve as a robust basis for the development of safe and efficient drugs for TB treatment, subject to further lead optimization and subsequent experimental validation.

Exploring the promising potential of noscapine for cancer and neurodegenerative disease therapy through inhibition of integrin-linked kinase-1

Integrin-linked kinase (ILK), a β1-integrin cytoplasmic domain interacting protein, supports multi-protein complex formation. ILK-1 is involved in neurodegenerative diseases by promoting neuro-inflammation. On the other hand, its overexpression induces epithelial-mesenchymal transition (EMT), which is a major hallmark of cancer and activates various factors associated with a tumorigenic phenotype. Thus, ILK-1 is considered as an attractive therapeutic target. We investigated the binding affinity and ILK-1 inhibitory potential of noscapine (NP) using spectroscopic and docking approaches followed by enzyme inhibition activity. A strong binding affinity of NP was measured for the ILK-1 with estimated Ksv (M-1) values of 1.9 × 105, 3.6 × 105, and 4.0 × 105 and ∆G0 values (kcal/mol) -6.19554, -7.8557 and -8.51976 at 298 K, 303 K, and 305 K, respectively. NP binds to ILK-1 with a docking score of -6.6 kcal/mol and forms strong interactions with active-site pocket residues (Lys220, Arg323, and Asp339). The binding constant for the interaction of NP to ILK-1 was 1.04 × 105 M-1, suggesting strong affinity and excellent ILK-1 inhibitory potential (IC50 of ∼5.23μM). Conformational dynamics of ILK-1 were also studied in the presence of NP. We propose that NP presumably inhibits ILK-1-mediated phosphorylation of various downstream signalling pathways that are involved in cancer cell survival and neuroinflammation.

Exploration of molecular interactions between scoparone and associated compounds with Constitutive androstane receptor (CAR) leading to gallstone prevention: an in silico investigation

Scoparone (6, 7 dimethylesculetin) is a biologically active compound derived from the herb Artemisia capillaris having anti-inflammatory, anti-lipemic, and anti-allergic roles. Activation of the constitutive androstane receptor (CAR) in primary hepatocytes of both wild-type and humanized CAR mice by scoparone, accelerates bilirubin and cholesterol clearance in vivo. This can prevent gallstones which is a dreaded gastrointestinal disease. To date, surgery is regarded as the gold standard for treating gallstones. The molecular interactions between scoparone and CAR leading to gallstone prevention are not yet explored. In this study, we have analyzed these interactions through an insilico approach. After extracting the CAR structures (mice and human) from the protein databank and 6, 7-dimethylesuletin from PubChem, energy minimization of both the receptors was done to make them stable followed by docking. Next, a simulation was performed to stabilize the docked complexes. Through docking, H-bonds and pi-pi interactions were found in the complexes, which imply a stable interaction, thus activating the CAR. A similarity search for scoparone was performed and the selected compounds were docked with the CAR receptors. Esculentin acetate and scopoletin acetate interacted with human CAR through pi-alkyl and H-bond respectively. While Fraxidin methyl ether, fraxinol methyl ether, and 6, 7 diethoxycoumarin interacted with mice CAR through H-bond and Pi-Pi T-shaped bonds. The selected complexes were simulated further. Our results are in accordance with the hypothesis in the literature. We have also analyzed the drug likeliness, absorption, non-carcinogenicity, and other properties of scoparone which can support further in vivo studies.Communicated by Ramaswamy H. Sarma.

Golgi Apparatus-Targeted Photodynamic Therapy for Enhancing Tumor Immunogenicity by Eliciting NLRP3 Protein-Dependent Pyroptosis

Innate and adaptive immunity is important for initiating and maintaining immune function. The nucleotide-binding oligomerization domain-like receptor family pyrin domain-containing 3 (NLRP3) inflammasome serves as a checkpoint in innate and adaptive immunity, promoting the secretion of pro-inflammatory cytokines and gasdermin D-mediated pyroptosis. As a highly inflammatory form of cell death distinct from apoptosis, pyroptosis can trigger immunogenic cell death and promote systemic immune responses in solid tumors. Previous studies proposed that NLRP3 was activated by translocation to the mitochondria. However, a recent authoritative study has challenged this model and proved that the Golgi apparatus might be a prerequisite for the activation of NLRP3. In this study, we first developed a Golgi apparatus-targeted photodynamic strategy to induce the activation of NLRP3 by precisely locating organelles. We found that Golgi apparatus-targeted photodynamic therapy could significantly upregulate NLRP3 expression to promote the subsequent release of intracellular proinflammatory contents such as IL-1β or IL-18, creating an inflammatory storm to enhance innate immunity. Moreover, this acute NLRP3 upregulation also activated its downstream classical caspase-1-dependent pyroptosis to enhance tumor immunogenicity, triggering adaptive immunity. Pyroptosis eventually led to immunogenic cell death, promoted the maturation of dendritic cells, and effectively activated antitumor immunity and long-lived immune memory. Overall, this Golgi apparatus-targeted strategy provided molecular insights into the occurrence of immunogenic pyroptosis and offered a platform to remodel the tumor microenvironment.

Designing a multi-epitope vaccine to control porcine epidemic diarrhea virus infection using immunoinformatics approaches

Porcine epidemic diarrhea virus (PEDV), a continuously evolving pathogen, causes severe diarrhea in piglets with high mortality rates. However, current vaccines cannot provide complete protection against PEDV, so vaccine development is still necessary and urgent. Here, with the help of immunoinformatics approaches, we attempted to design a multi-epitope vaccine named rPMEV to prevent and control PEDV infection. The epitopes of rPMEV were constructed by 9 cytotoxic T lymphocyte epitopes (CTLs), 11 helper T lymphocyte epitopes (HTLs), 6 linear B cell epitopes (LBEs), and 4 conformational B cell epitopes (CBEs) based on the S proteins from the four representative PEDV G2 strains. To enhance immunogenicity, porcine β-defensin-2 (PBD-2) was adjoined to the N-terminal of the vaccine as an adjuvant. All of the epitopes and PBD-2 were joined by corresponding linkers and recombined into the multivalent vaccine, which is stable, antigenic, and non-allergenic. Furthermore, we adopted molecular docking and molecular dynamics simulation methods to analyze the interaction of rPMEV with the Toll-like receptor 4 (TLR4): a stable interaction between them created by 13 hydrogen bonds. In addition, the results of the immune simulation showed that rPMEV could stimulate both cellular and humoral immune responses. Finally, to raise the expression efficiency, the sequence of the vaccine protein was cloned into the pET28a (+) vector after the codon optimization. These studies indicate that the designed multi-epitope vaccine has a potential protective effect, providing a theoretical basis for further confirmation of its protective effect against PEDV infection in vitro and in vivo studies.

Identification of a potential inhibitor for New Delhi metallo-β-lactamase 1 (NDM-1) from FDA approved chemical library- a drug repurposing approach to combat carbapenem resistance

Superbugs producing New Delhi metallo-β-lactamase 1 (NDM-1) enzyme is a growing crisis, that is adversely affecting the global health care system. NDM-1 empowers the bacteria to inactivate entire arsenal of β-lactam antibiotics including carbapenem (the last resort antibiotic) and remains ineffective to all the available β lactamase inhibitors used in the clinics. Limited therapeutic option available for rapidly disseminating NDM-1 producing bacteria makes it imperative to identify a potential inhibitor for NDM-1 enzyme. With drug repurposing approach, in this study, we used virtual screening of available Food and Drug Administration (FDA) approved chemical library (ZINC12 database) and captured 'adapalene' (FDA drug) as a potent inhibitor candidate for NDM-1 enzyme. Active site docking with NDM-1, showed adapalene with binding energy -9.21 kcal/mol and interacting with key amino acid residues (Asp124, His122, His189, His250, Cys208) in the active site of NDM-1. Further, molecular dynamic simulation of NDM-1 docked with the adapalene at 100 ns displayed a stable conformation dynamic, with relative RMSD and RMSF in the acceptable range. Subsequently, in vitro enzyme assays using recombinant NDM-1 protein demonstrated inhibition of NDM-1 by adapalene. Further, the combination of adapalene plus meropenem (carbapenem antibiotic) showed synergistic effect against the NDM-1 producing carbapenem (meropenem) resistant clinical isolates (Escherichia coli and Klebsiella pneumoniae). Overall, our data indicated that adapalene can be a potential inhibitor candidate for NDM-1 enzyme that can contribute to the development of a suitable adjuvant to save the activity of carbapenem antibiotic against infections caused by NDM-1 positive gram-negative bacteria. Communicated by Ramaswamy H. Sarma.

Computational insights into the role of structurally diverse plant secondary metabolites as inhibitors against Imidazole Glycerol Phosphate Dehydratase of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is one of the major causes of death worldwide and there is a pressing need for the development of novel drug leads. The Imidazole Glycerol Phosphate Dehydratase (IGPD) of Mtb is one of the key enzymes in the histidine biosynthesis pathway and has been recognized as the potentially underexploited drug target for anti-tuberculosis treatment. In the present study, 6063 structurally diverse plant secondary metabolites (PSM) were screened for their efficiency in inhibiting the catalytic activity of IGPD through molecular docking. The top 150 PSMs with the lowest binding energy represent the chemical classes, including Tannins (34%), Flavonoid Glycosides (14%), Terpene Glycosides (10%), Steroid Lactones (9.3%), Flavonoids (6.6%), Steroidal Glycosides (4.6%), etc. Bismahanine, Ashwagandhanolide, and Daurisoline form stable IGPD-inhibitor complexes with binding free energies of -291.3 ± 16.5, -279.0 ± 25.0, and -279.8 ± 17.6 KJ/mol, respectively, as determined by molecular dynamics simulations. These PSM demonstrated strong H-bond interactions with the amino acid residues Ile279, Arg281, and Lys276 in the catalytic region of IGPD, as revealed by structural snapshots. On the basis of our findings, these three PSM could be considered as possible leads against IGPD and should be explored in vitro and in vivo.Communicated by Ramaswamy H. Sarma.

In silico study on evaluation of corosolic acid of Lagerstroemia speciosa against Alzheimer’s disease

Purpose

Alzheimer’s disease (AD), the most common cause of dementia, is a neurodegenerative disorder caused by the aggregation of amyloid-beta (Aβ) at outside of neuron cells and also due to tau aggregation inside the cell. Corosolic acid is aimed to be selected as a main active constituent of Lagerstroemia speciosa for the study.

Design/methodology/approach

In the present study, molecular docking of corosolic acid and tau protein was examined using PyRx-v.0.8 software. Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties were described and a molecular dynamics study of the bound complex was performed using Desmond.

Findings

The docking score and interactions suggested that the corosolic acid (CID:6918774) could bind to tau protein to prevent the fibrillar network, to prevent AD. During simulation corosolic acid-bound protein root mean square deviation (RMSD) values showed more stability when compared to the Apo form of protein. Molecular dynamics study of tau protein and corosolic acid complex gave the insights to develop a drug-like candidate against AD.

Originality/value

The use of corosolic acid of Lagerstroemia speciosa to prevent AD is supported by preliminary analysis on a computational basis. This compound should explore in terms of experimental strategies for the further drug development process. However, in vitro and in vivo evaluation studies are required to suggest the use of corosolic acid against AD.

Mechanistic insights into the role of vitamin D and computational identification of potential lead compounds for Parkinson's disease

Parkinson's disease (PD) is the second most common neurodegenerative disorder that affects dopaminergic neurons in the midbrain. A recent study suggests that Orphan Nuclear Receptor 1 (NURR1) impairment may contribute to PD pathogenesis. Our study found three potent agonists for NURR1 protein based on structural and ligand-based screening methods. The pharmacophore is comprised of a hydrogen bond donor, a hydrophobic group, and two aromatic rings (DHRR). The Pharmacophore screening method screened 3142 compounds, of which 3 were screened using structure-based screening. An analysis of the molecules using Molecular Mechanics-Generalized Born Surface Area (binding free energy) revealed a range of -46.77 to -59.06 Kcal/mol. After that, chemical reactivity was investigated by density functional theory, and molecular dynamics simulation was performed (protein-ligand stability). Based on the computational studies, Lifechemical_16901310, Maybridge_2815310, and NPACT_392450 are promising agonists with respect to NURR1. To confirm the potency of the identified compounds, further validation and experiments must be conducted.

Machine learning driven drug repurposing strategy for identification of potential RET inhibitors against non-small cell lung cancer

Non-small cell lung cancer (NSCLC) remains the leading cause of mortality and morbidity worldwide accounting about 85% of total lung cancer cases. The receptor REarranged during Transfection (RET) plays an important role by ligand independent activation of kinase domain resulting in carcinogenesis. Presently, the treatment for RET driven NSCLC is limited to multiple kinase inhibitors. This situation necessitates the discovery of novel and potent RET specific inhibitors. Thus, we employed high throughput screening strategy to repurpose FDA approved compounds from DrugBank comprising of 2509 molecules. It is worth noting that the initial screening is accomplished with the aid of in-house machine learning model built using IC₅₀ values corresponding to 2854 compounds obtained from BindingDB repository. A total of 497 compounds (19%) were predicted as actives by our generated model. Subsequent in silico validation process such as molecular docking, MMGBSA and density function theory analysis resulted in identification of two lead compounds named DB09313 and DB00471. The simulation study highlights the potency of DB00471 (Montelukast) as potential RET inhibitor among the investigated compounds. In the end, the half-minimal inhibitory activity of montelukast was also predicted against RET protein expressing LC-2/ad cell lines demonstrated significant anticancer activity. Collective analysis from our study highlights that montelukast could be a promising candidate for the management of RET specific NSCLC.

Novel In Silico Insights into Rv1417 and Rv2617c as Potential Protein Targets: The Importance of the Medium on the Structural Interactions with Exported Repetitive Protein (Erp) of Mycobacterium tuberculosis

Nowadays, tuberculosis is the second leading cause of death from a monopathogenic transmitted disease, only ahead of COVID-19. The role of exported repetitive protein (Erp) in the virulence of Mycobacterium tuberculosis has been extensively demonstrated. In vitro and in vivo assays have identified that Erp interacts with Rv1417 and Rv2617c proteins, forming putative transient molecular complexes prior to localization to the cell envelope. Although new insights into the interactions and functions of Erp have emerged over the years, knowledge about its structure and protein–protein interactions at the atomistic level has not been sufficiently explored. In this work, we have combined several in silico methodologies to gain new insights into the structural relationship between these proteins. Two system conditions were evaluated by MD simulations: Rv1417 and Rv2617c embedded in a lipid membrane and another with a semi-polar solvent to mimic the electrostatic conditions on the membrane surface. The Erp protein was simulated as an unanchored structure. Stabilized structures were docked, and complexes were evaluated to recognize the main residues involved in protein–protein interactions. Our results show the influence of the medium on the structural conformation of proteins. Globular conformations were favored under high polarity conditions and showed a higher energetic affinity in complex formation. Meanwhile, disordered conformations were favored under semi-polar conditions and an increase in the number of contacts between residues was observed. In addition, the electrostatic potential analysis showed remarkable changes in protein interactions due to the polarity of the medium, demonstrating the relevance of Erp protein in heterodimer formation. On the other hand, contact analysis showed that several C-terminal residues of Erp were involved in the protein interactions, which seems to contradict experimental observations; however, these complexes could be transient forms. The findings presented in this work are intended to open new perspectives in the studies of Erp protein molecular interactions and to improve the knowledge about its function and role in the virulence of Mycobacterium tuberculosis.

Designing Novel Compounds for the Treatment and Management of RET-Positive Non-Small Cell Lung Cancer-Fragment Based Drug Design Strategy

Rearranged during transfection (RET) is an oncogenic driver receptor that is overexpressed in several cancer types, including non-small cell lung cancer. To date, only multiple kinase inhibitors are widely used to treat RET-positive cancer patients. These inhibitors exhibit high toxicity, less efficacy, and specificity against RET. The development of drug-resistant mutations in RET protein further deteriorates this situation. Hence, in the present study, we aimed to design novel drug-like compounds using a fragment-based drug designing strategy to overcome these issues. About 18 known inhibitors from diverse chemical classes were fragmented and bred to form novel compounds against RET proteins. The inhibitory activity of the resultant 115 hybrid molecules was evaluated using molecular docking and RF-Score analysis. The binding free energy and chemical reactivity of the compounds were computed using MM-GBSA and density functional theory analysis, respectively. The results from our study revealed that the developed hybrid molecules except for LF21 and LF27 showed higher reactivity and stability than Pralsetinib. Ultimately, the process resulted in three hybrid molecules namely LF1, LF2, and LF88 having potent inhibitory activity against RET proteins. The scrutinized molecules were then subjected to molecular dynamics simulation for 200 ns and MM-PBSA analysis to eliminate a false positive design. The results from our analysis hypothesized that the designed compounds exhibited significant inhibitory activity against multiple RET variants. Thus, these could be considered as potential leads for further experimental studies.

Itaconate: an antimicrobial metabolite of macrophages

Itaconate is a conjugated 1,4-dicarboxylate produced by macrophages. This small molecule has recently received increasing attention due to its role in modulating the immune response of macrophages upon exposure to pathogens. Itaconate has also been proposed to play an antimicrobial function; however, this has not been explored as intensively. Consistent with the latter, itaconate is known to show antibacterial activity in vitro and was reported to inhibit isocitrate lyase, an enzyme required for survival of bacterial pathogens in mammalian systems. Recent studies have revealed bacterial growth inhibition under biologically relevant conditions. In addition, an antimicrobial role for itaconate is substantiated by the high concentration of itaconate found in bacteria-containing vacuoles, and by the production of itaconate-degrading enzymes in pathogens such as Salmonella enterica ser. Typhimurium, Pseudomonas aeruginosa, and Yersinia pestis. This review describes the current state of literature in understanding the role of itaconate as an antimicrobial agent in host–pathogen interactions.

Pathway Driven Target Selection in Klebsiella pneumoniae: Insights Into Carbapenem Exposure

Carbapenem-resistant Klebsiella pneumoniae (CR-KP) represents an emerging threat to public health. CR-KP infections result in elevated morbidity and mortality. This fact, coupled with their global dissemination and increasingly limited number of therapeutic options, highlights the urgency of novel antimicrobials. Innovative strategies linking genome-wide interrogation with multi-layered metabolic data integration can accelerate the early steps of drug development, particularly target selection. Using the BioCyc ontology, we generated and manually refined a metabolic network for a CR-KP, K. pneumoniae Kp13. Converted into a reaction graph, we conducted topological-based analyses in this network to prioritize pathways exhibiting druggable features and fragile metabolic points likely exploitable to develop novel antimicrobials. Our results point to the aptness of previously recognized pathways, such as lipopolysaccharide and peptidoglycan synthesis, and casts light on the possibility of targeting less explored cellular functions. These functions include the production of lipoate, trehalose, glycine betaine, and flavin, as well as the salvaging of methionine. Energy metabolism pathways emerged as attractive targets in the context of carbapenem exposure, targeted either alone or in conjunction with current therapeutic options. These results prompt further experimental investigation aimed at controlling this highly relevant pathogen.

Discovery of a Potent Candidate for RET-Specific Non-Small-Cell Lung Cancer-A Combined In Silico and In Vitro Strategy

Rearranged during transfection (RET) is a tyrosine kinase oncogenic receptor, activated in several cancers including non-small-cell lung cancer (NSCLC). Multiple kinase inhibitors vandetanib and cabozantinib are commonly used in the treatment of RET-positive NSCLC. However, specificity, toxicity, and reduced efficacy limit the usage of multiple kinase inhibitors in targeting RET protein. Thus, in the present investigation, we aimed to figure out novel and potent candidates for the inhibition of RET protein using combined in silico and in vitro strategies. In the present study, screening of 11,808 compounds from the DrugBank repository was accomplished by different hypotheses such as pharmacophore, e-pharmacophore, and receptor cavity-based models in the initial stage. The results from the different hypotheses were then integrated to eliminate the false positive prediction. The inhibitory activities of the screened compounds were tested by the glide docking algorithm. Moreover, RF score, Tanimoto coefficient, prime-MM/GBSA, and density functional theory calculations were utilized to re-score the binding free energy of the docked complexes with high precision. This procedure resulted in three lead molecules, namely DB07194, DB03496, and DB11982, against the RET protein. The screened lead molecules together with reference compounds were then subjected to a long molecular dynamics simulation with a 200 ns time duration to validate the inhibitory activity. Further analysis of compounds using MM-PBSA and mutation studies resulted in the identification of potent compound DB07194. In essence, a cell viability assay with RET-specific lung cancer cell line LC-2/ad was also carried out to confirm the in vitro biological activity of the resultant compound, DB07194. Indeed, the results from our study conclude that DB07194 can be effectively translated for this new therapeutic purpose, in contrast to the properties for which it was originally designed and synthesized.

An Insilico evaluation of phytocompounds from Albizia amara and Phyla nodiflora as cyclooxygenase-2 enzyme inhibitors

PURPOSE

The enzyme Cyclooxygenases (COX-1 and COX-2) catalyze the formation of prostaglandin, a mediator of the inflammatory pathway. Inflammation related pathological conditions may be alleviated by targeting the Cox enzymes.COX-2 inhibitors that are currently available in the market causes undesirable side effects. Our present study focuses on the in-silico inhibition of COX -2 enzyme by the phytocompounds from Albizia amara and Phyla nodiflora.

METHODS

The phytochemicals present in Albizia amara and Phyla nodiflora were analyzed for their COX-2 inhibition potential. Eight compounds from Albizia amara and eleven compounds from Phyla nodiflora obtained from GC-MS analysis was used for the current study. Molecular docking was performed using AutoDock vina. The crystal structure of COX-2 (PDB ID: 5IKR) was obtained from Protein data bank. PyMol was used to remove any solvent, organic and inorganic molecules. Energy minimization of the protein was carried out using SPDBV software. Geometrical optimizations of the ligands were performed using Avogadro software. Celecoxib was used as the positive control. ADMET properties of the compounds were analyzed using SwissADME and ProtoxII online servers. Molecular mechanics/generalized born surface area (MM/GBSA) calculations were performed to evaluate the binding efficiency. Molecular dynamics of the protein and protein-ligand complex was studied for about 100 ns using Desmond package of Schrodinger suite.

RESULTS

Among the eighteen compounds, Squalene present in both the plants showed a better binding energy of -7.7 kcal/mol, when compare to other phytocompounds present in the extract. The control celecoxib showed a binding energy of about - 9.4 kcal/mol. The toxicity and ADMET properties of squalene indicated that it is non-toxic and followed Lipinski's rule. Molecular Dynamics (MD) analysis showed that the binding of squalene to the enzyme was stable.

CONCLUSION

Squalene could potentially inhibit COX2 and o wing to its properties, squalene can be formulated in gels/creams and could be possibly used for external edema and inflammation.

Study of lipase producing gene in wheat - an in silico approach

BACKGROUND

Lipases (EC 3.1.1.3) catalyze the hydrolysis of oil into free fatty acids and glycerol forming the 3rd largest group of commercialized enzymes. Plant lipases grab attention recently because of their specificity, less production and purified cost, and easy availability. In silico approach is the first step to identify different genes coding for lipase in a most common indigenous plant, wheat, to explore the possibility of this plant as an alternative source for commercial lipase production. As the hierarchy organization of genes reflects an ancient process of gene duplication and divergence, many of the theoretical and analytical tools of the phylogenetic systematics can be utilized for comparative genomic studies. Also, in addition to experimental identification and characterization of genes, for computational genomic analysis, Arabidopsis has become a popular strategy to identify crop genes which are economically important, as Arabidopsis genes had been well identified and characterized for lipase. A number of articles had been reported in which genes of wheat have shown strong homology with Arabidopsis. The complete genome sequences of rice and Arabidopsis constitute a valuable resource for comparative genome analysis as they are representatives of the two major evolutionary lineages within the angiosperms. Here, in this in silico approach, Arabidopsis and Oryza sativa serve as models for dicotyledonous and monocotyledonous species, respectively, and the genomic sequence data available was used to identify the lipase genes in wheat.

RESULTS

In this present study, Ensembl Plants database was explored for lipase producing gene present in wheat genome and 21 genes were screened down as they contain specific domain and motif for lipase (GXSXG). According to the evolutionary analysis, it was found that the gene TraesCS5B02G157100, located in 5B chromosome, has 58.35% sequence similarity with the reported lipase gene of Arabidopsis thaliana and gene TraesCS3A02G463500 located in the 3A chromosome has 51.74% sequence similarity with the reported lipase gene of Oryza sativa. Homology modeling was performed using protein sequences coded by aforementioned genes and optimized by molecular dynamic simulations. Further with the help of molecular docking of modeled structures with tributyrin, binding efficiency was checked, and the difference in energies (DE) was -9.83 kcal/mol and -6.67 kcal/mol, respectively.

CONCLUSIONS

The present work provides a basic understanding of the gene-encoding lipase in wheat, which could be easily accessible and used as a potent industrial enzyme. The study enlightens another direction which can be used further to explore plant lipases.

High-throughput screening of natural compounds and inhibition of a major therapeutic target HsGSK-3β for Alzheimer's disease using computational approaches

BACKGROUND

Alzheimer's disease is a leading neurodegenerative disease worldwide and is the 6th leading cause of death in the USA. AD is a very complex disease and the drugs available in the market cannot fully cure it. The glycogen synthase kinase 3 beta plays a major role in the hyperphosphorylation of tau protein which forms the neurofibrillary tangles which is a major hallmark of AD. In this study, we have used a series of computational approaches to find novel inhibitors against GSK-3β to reduce the TAU hyperphosphorylation.

RESULTS

We have retrieved a set of compounds (n=167,741) and screened against GSK-3β in four sequential steps. The resulting analysis of virtual screening suggested that 404 compounds show good binding affinity and can be employed for pharmacokinetic analysis. From here, we have selected 20 compounds those were good in terms of pharmacokinetic parameters. All these compounds were re-docked by using Autodock Vina followed by Autodock. Four best compounds were employed for MDS and here predicted RMSD, RMSF, Rg, hydrogen bonds, SASA, PCA, and binding-free energy. From all these analyses, we have concluded that out of 167,741 compounds, the ZINC15968620, ZINC15968622, and ZINC70707119 can act as lead compounds against HsGSK-3β to reduce the hyperphosphorylation.

CONCLUSION

The study suggested three compounds (ZINC15968620, ZINC15968622, and ZINC70707119) have great potential to be a drug candidate and can be tested using in vitro and in vivo experiments for further characterization and applications.

Repurposing of FDA approved drugs and their validation against potential drug targets for Salmonella enterica through molecular dynamics simulation

Salmonella is a widely distributed pathogen causing infection of intestinal tract, typhoid, and paratyphoid fever. Number of drugs was developed against salmonella, but in the last few decades due to the emergence of drug resistant strains, most of these drugs became dormant. As a result Salmonellosis emerges as a trivial cause of human mortality worldwide; therefore, there is an urgent need for unexploited drug targets and drugs to treat Salmonellosis. As development of new drug molecules is very time consuming and costly, drug repurposing is in consideration as a better alternative. With the aim to identify a new drug molecule against the Salmonella through repurposing approach, we utilized 14 well reported druggable targets known to play a vital role in the life cycle of pathogens. These targets were used to screen DrugBank and got 53 FDA approved drugs against them. To find the interaction between considered target proteins and screened drugs, molecular docking was performed. Fourteen docked drug-target complexes with reasonable binding affinities were subjected to Molecular Dynamics Simulation (MDS) at 150 ns, using Amber18. At the end MMPBSA and MMGBSA calculations were performed for all stable complexes and finally, got 3 precise and favourable complexes, i.e. ArcB-Cefpiramide, MrcB-Cefoperazone, and PhoQ-Carindacillin. Rigorous structural and energetic analysis for these complexes validates the potential of drug molecules to act as therapeutic drugs against Salmonella enterica. With this study we hypothesize that the drugs Cefpiramide (DB00430), Cefoperazone (DB01329) and Carindacillin (DB09319) will be the good repurposed-drugs for the treatment of Salmonellosis.Communicated by Ramaswamy H. Sarma.

Identification of potent Antigen 85C inhibitors of Mycobacterium tuberculosis via in-house lichen library and binding free energy studies Part-II

Tuberculosis remains the cause of mortality throughout the world. Currently, the available anti-tubercular drugs are not effective because of the existence of Multi-Drug resistant tuberculosis (MDR-TB) and Extensively-Drug resistant tuberculosis (XDR-TB). It has, therefore, become necessary to develop novel drugs that inhibit the activity of drug-resistant Mycobacterium tuberculosis. Due to the existence of MDR and XDR-TB, Mtb Ag85C has risen out as a propitious molecular drug target as it has importance in the synthesis of main components of the Mtb cell envelope which are essential for the virulence and survival of Mtb. In a previous paper, we studied a potential drug target by virtual high throughput screening of compounds and in continuation of the study on Mtb Ag85C, we further studied the role of lichen compounds in the inhibition of Ag85C. In the current research work, virtual screening of a lichen compounds library was performed against Ag85C. Further, ADMET analysis was employed to filter out the screened lichen compounds. Bioactivity score and toxicity prediction finalized four lichen compounds i.e. Portentol, Aspicilin, Parietinic acid and Polyporic acid as potential inhibitors of Ag85C. The stability and dynamic behavior of four compounds were analyzed by using Molecular dynamics simulation which indicated that they may be potential inhibitors of Ag85C. Therefore, based on the above results, Portentol, Aspicilin, Parietinic acid and Polyporic acid may be potential drug candidates against Mtb. We suggest that the use of these compounds can minimize the treatment time-period and the various side effects associated with the currently available anti-tubercular drugs.

Well-dispersed platelet lysate entrapped nanoparticles incorporate with injectable PDLLA-PEG-PDLLA triblock for preferable cartilage engineering application

Platelet lysate (PL) as a cost-effective cocktail of growth factors is an emerging ingredient in regenerative medicine, especially in cartilage tissue engineering. However, most studies fail to pay attention to PL's intrinsic characteristics and incorporate it directly with scaffolds or hydrogels by simple mixture. Currently, the particle size distribution of PL was determined to be scattered. Directly introducing PL into a thermosensitive poly(d,l-lactide)-poly(ethylene glycol)-poly(d,l-lactide) (PLEL) hydrogel disturbed its sol-gel transition. Electrostatic self-assembly heparin (Hep) and ε-poly-l-lysine (EPL) nanoparticles (NPs) were fabricated to improve the dispersity of PL. Such PL-NPs-incorporated PLEL gels retained the initial gelling capacity and showed a long-term PL-releasing ability. Moreover, the PL-loaded composite hydrogels inhibited the inflammatory response and dedifferentiation of IL-1β-induced chondrocytes. For in vivo applications, the PLEL@PL-NPs system ameliorated the early cartilage degeneration and promoted cartilage repair in the late stage of osteoarthritis. RNA sequencing analysis indicated that PL's protective effects might be associated with modulating hyaluronan synthase 1 (HAS-1) expression. Taken together, these results suggest that well-dispersed PL by Hep/EPL NPs is a preferable approach for its incorporation into hydrogels and the constructed PLEL@PL-NPs system is a promising cell-free and stepwise treatment option for cartilage tissue engineering.

Enhanced Production of Active Ecumicin Component with Higher Antituberculosis Activity by the Rare Actinomycete Nonomuraea sp. MJM5123 Using a Novel Promoter-Engineering Strategy

Ecumicins are potent antituberculosis natural compounds produced by the rare actinomycete Nonomuraea sp. MJM5123. Here, we report an efficient genetic manipulation platform of this rare actinomycete. CRISPR/Cas9-based genome editing was achieved based on successful sporulation. Two genes in the ecumicin gene cluster were further investigated, ecuN and ecuE, which potentially encode a pretailoring cytochrome P450 hydroxylase and the core peptide synthase, respectively. Deletion of ecuN led to an enhanced ratio of the ecumicin compound EcuH16 relative to that of EcuH14, indicating that EcuN is indeed a P450 hydroxylase, and there is catalyzed hydroxylation at the C-3 position in unit₁₂ phenylalanine to transform EcuH16 to the compound EcuH14. Furthermore, promoter engineering of ecuE by employing the strong promoter kasO*P was performed and optimized. We found that integrating the endogenous ribosome-binding site (RBS) of ecuE together with the RBS from kasO*P led to improved ecumicin production and resulted in a remarkably high EcuH16/EcuH14 ratio. Importantly, production of the more active component EcuH16 was considerably increased in the double RBSs engineered strain EPR1 compared to that in the wild-type strain, reaching 310 mg/L. At the same time, this production level was 2.3 times higher than that of the control strain EPA1 with only one RBS from kasO*P. To the best of our knowledge, this is the first report of genome editing and promoter engineering on the rare actinomycete Nonomuraea.

Demystifying the catalytic pathway of Mycobacterium tuberculosis isocitrate lyase

Pulmonary tuberculosis, caused by Mycobacterium tuberculosis, is one of the most persistent diseases leading to death in humans. As one of the key targets during the latent/dormant stage of M. tuberculosis, isocitrate lyase (ICL) has been a subject of interest for new tuberculosis therapeutics. In this work, the cleavage of the isocitrate by M. tuberculosis ICL was studied using quantum mechanics/molecular mechanics method at M06-2X/6-31+G(d,p): AMBER level of theory. The electronic embedding approach was applied to provide a better depiction of electrostatic interactions between MM and QM regions. Two possible pathways (pathway I that involves Asp108 and pathway II that involves Glu182) that could lead to the metabolism of isocitrate was studied in this study. The results suggested that the core residues involved in isocitrate catalytic cleavage mechanism are Asp108, Cys191 and Arg228. A water molecule bonded to Mg²⁺ acts as the catalytic base for the deprotonation of isocitrate C(2)–OH group, while Cys191 acts as the catalytic acid. Our observation suggests that the shuttle proton from isocitrate hydroxyl group C(2) atom is favourably transferred to Asp108 instead of Glu182 with a lower activation energy of 6.2 kcal/mol. Natural bond analysis also demonstrated that pathway I involving the transfer of proton to Asp108 has a higher intermolecular interaction and charge transfer that were associated with higher stabilization energy. The QM/MM transition state stepwise catalytic mechanism of ICL agrees with the in vitro enzymatic assay whereby Asp108Ala and Cys191Ser ICL mutants lost their isocitrate cleavage activities.

Discovery of Potential Chemical Probe as Inhibitors of CXCL12 Using Ligand-Based Virtual Screening and Molecular Dynamic Simulation

CXCL12 are small pro-inflammatory chemo-attractant cytokines that bind to a specific receptor CXCR4 with a role in angiogenesis, tumor progression, metastasis, and cell survival. Globally, cancer metastasis is a major cause of morbidity and mortality. In this study, we targeted CXCL12 rather than the chemokine receptor (CXCR4) because most of the drugs failed in clinical trials due to unmanageable toxicities. Until now, no FDA approved medication has been available against CXCL12. Therefore, we aimed to find new inhibitors for CXCL12 through virtual screening followed by molecular dynamics simulation. For virtual screening, active compounds against CXCL12 were taken as potent inhibitors and utilized in the generation of a pharmacophore model, followed by validation against different datasets. Ligand based virtual screening was performed on the ChEMBL and in-house databases, which resulted in successive elimination through the steps of pharmacophore-based and score-based screenings, and finally, sixteen compounds of various interactions with significant crucial amino acid residues were selected as virtual hits. Furthermore, the binding mode of these compounds were refined through molecular dynamic simulations. Moreover, the stability of protein complexes, Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), and radius of gyration were analyzed, which led to the identification of three potent inhibitors of CXCL12 that may be pursued in the drug discovery process against cancer metastasis.

Molecular Insights into the Interaction of Ursolic Acid and Cucurbitacin from Colocynth with Therapeutic Targets of Mycobacterium tuberculosis

Aims:

Medicinal plants like Citrullus colocynthis are a potential choice to produce helpful novel antimycobacterial drugs. The existence of a range of natural products in the plants, especially Ursolic Acid (UA) and cucurbitacin E 2-0-β-d-glucopyranoside (CEG), with promising antibacterial activity against a variety of bacteria, prompted the need to check its actions against Mycobacterium tuberculosis (Mtb).

Background:

Mycobacterium tuberculosis (Mtb), an obligate human pathogen causes tuberculosis and is one of the major causes of death worldwide. A few combinations of drugs are currently accessible for treating TB patients, but these are inadequate to tackle worldwide TB cases.

Objective:

The molecular interactions between ursolic acid and cucurbitacin E with the eight potential Mtb target proteins were investigated with the objective of finding drug-like inhibitors.

Methods:

Avogadro v.1.2.0 and Openbabel v.2.4.1 were used for creating file formats required for docking analysis. Molecular docking was performed with eight different proteins essential for Mtb metabolism and survival. AutoDock v.4.2 and AutoDock vina v.1.1.2 were used for docking and Gromacs 5.1.4 was used for simulation studies.

Results and Discussion:

Among the two ligands used in this research, cucurbitacin E showed a better docking score relative to the drugs presently available for all the target proteins. Rifampicin showed the best binding affinity (among known inhibitors) i.e. -10.8 kcal/mol with C terminal caspase recruitment domain. Moreover, ursolic acid and cucurbitacin E showed uniform binding score (above -7.5 kcal/mol) with all the target proteins, acknowledged its availability as a potential multi-target drug.

Conclusion:

Ursolic acid can be useful in the creation of novel, multi-targeted and effective anti- TB medicines since it showed stable structure with FabH.

In-silico anti-inflammatory potential of guaiane dimers from Xylopia vielana targeting COX-2

Natural products of herbal origin are prodigious to display diverse pharmacological activities. In the present study, five guaiane-type sesquiterpene dimers, xylopidimers A - E (1-5), isolated from Xylopia vielana species were tested against COX-2 protein target (PDB: 1CX2), a potent target for anti-inflammatory agents. To better understand the pharmacological properties of all these compounds, in this work, a systemic in silico study was performed on xylopidimers A-E using molecular docking, ADMET analysis and MD simulations. During ADMET predictions the two compounds xylopidimer C, D displayed best results as compared to others. The compound xylopidimer C was further evaluated for its MD simulations and its molecular interactions with COX2 complex showed clear interactions with active gorge of the enzyme through hydrogen bonding as well as hydrophobic contacts. The xylopidimer C has shown the best binding potential with -10.57Kcal/mol energy with 17.92 nano molar of predicted inhibition constant better than Ibuprofen and Felbinac. These findings provide enough significant information for designing and developing novel targeted base anti-inflammatory drugs from guaiane dimers.

In silico identification of natural fungicide from Melia azedarach against isocitrate lyase of Fusarium graminearum

Isocitrate Lyase (ICL) is a crucial enzyme involved in the Glyoxylate pathway, essential for the virulence of several fungal pathogens including Fusarium graminearum. ICL is a novel target for the discovery of antifungal compounds and F. graminearum ICL inhibitors can be used to control the growth of this fungus. Although, several inhibitors of ICL have been identified, however, most of these inhibitors are not environment-friendly. Hence there is still a need to discover natural inhibitors of ICL that can be more effective. To identify a potential antifungal compound, we performed a structure-based screening of phytochemicals of Melia azedarach against the FgICL structure by molecular docking and 104 ligands were found to have a better docking score as compared to the reference molecule. These compounds were assessed for drug-likeness and ADMET prediction. After molecular docking, drug-likeness and toxicity analysis, six potential compounds (Melianoninol (-6.6 kcal/mol), Nimbinene (-7.7 kcal/mol), Vilasinin (-8.1 kcal/mol), Fraxinellone (-6.7 kcal/mol), Gedunin (-7.8 kcal/mol), and Meldenin (-7.8 kcal/mol)) were subjected for rescoring by X-Score. The structural stability and dynamics of screened compounds at the active site of FgICL were examined using MD simulation and MM-PBSA analysis. The result of MM-PBSA revealed that four phytochemicals viz. Melianoninol, Nimbinene, Vilasinin, and Fraxinellone had binding free energy of -17.25 kcal/mol, -59.35 kcal/mol, -64.79 kcal/mol, and -29.86 kcal/mol, respectively. Molecular dynamics simulation and MM-PBSA demonstrated that these four phytochemicals displayed considerable significant structural and pharmacological properties and could be probable antifungal drug candidates against F. graminearum. These phyotchemicals of M. azedarach may be suitable candidates for further experimental analysis. [Formula: see text]Communicated by Ramaswamy H. Sarma.

Structure-based discovery of phenyl-diketo acids derivatives as Mycobacterium tuberculosis malate synthase inhibitors

Mycobacterium tuberculosis remains one of the most successful bacterial pathogens worldwide. The development of drug-resistant strains and the ability of the bacteria to persist in a latent form in the host are major problems for tuberculosis (TB) control. Glyoxylate shunt is a metabolic bypass of the Krebs cycle and is the key for M. tuberculosis to survive under latent conditions. Malate synthase (MtbMS) catalyzes the second step of the glyoxylate cycle and converts glyoxylate into malate. Phenyl-diketo acid (PDKA) is a potent inhibitor of MtbMS, and its efficacy is validated in a mouse model of TB. To identify novel PDKA analogs as anti-TB compounds, PDKA analogs that obeyed the Lipinski rules (n = 5473) were analyzed and docked with MtbMS structure in three sequential modes. These compounds were then assessed for ADMET parameters. Of the compounds examined, 19 were found to fit well for redocking studies. After optimization, four prospective inhibitors were identified, that along with the reference compound PDKA were subjected to 50 ns molecular dynamics simulation and binding-free energy analyses to evaluate the complex dynamics after ligand binding, the stability of the bound complexes, and the intermolecular interactions between the complexes. The MtbMS-PDKA complex showed the binding free energy of -57.16 kJ·mol^-1. After a thorough analysis, our results suggested that three compounds which had binding-free energy of -127.96, -97.60, and -83.98 kJ·mol^-1, with PubChem IDs 91937661, 14016246, and 126487337, respectively, have the potential to inhibit MtbMS and can be taken as lead compounds for drug discovery against TB.Communicated by Ramaswamy H. Sarma.

Applications of Molecular Simulation in the Discovery of Antituberculosis Drugs: A Review

After decades of efforts, tuberculosis has been well controlled in most places. The existing drugs are no longer sufficient for the treatment of drug-resistant Mycobacterium tuberculosis due to significant toxicity and selective pressure, especially for XDR-TB. In order to accelerate the development of high-efficiency, low-toxic antituberculosis drugs, it is particularly important to use Computer Aided Drug Design (CADD) for rational drug design. Here, we systematically reviewed the specific role of molecular simulation in the discovery of new antituberculosis drugs. The purpose of this review is to overview current applications of molecular simulation methods in the discovery of antituberculosis drugs. Furthermore, the unique advantages of molecular simulation was discussed in revealing the mechanism of drug resistance. The comprehensive use of different molecular simulation methods will help reveal the mechanism of drug resistance and improve the efficiency of rational drug design. With the help of molecular simulation methods such as QM/MM method, the mechanisms of biochemical reactions catalyzed by enzymes at atomic level in Mycobacterium tuberculosis has been deeply analyzed. QSAR and virtual screening both accelerate the development of highefficiency, low-toxic potential antituberculosis drugs. Improving the accuracy of existing algorithms and developing more efficient new methods for CADD will always be a hot topic in the future. It is of great value to utilize molecular dynamics simulation to investigate complex systems that cannot be studied in experiments, especially for drug resistance of Mycobacterium tuberculosis.

Strategic Development of a Next-Generation Multi-Epitope Vaccine To Prevent Nipah Virus Zoonotic Infection

Nipah virus (NiV) is an emerging zoonotic pathogen, reported for the recent severe outbreaks of encephalitis and respiratory illness in humans and animals, respectively. Many antiviral drugs have been discovered to inhibit this pathogen, but none of them were that much efficient. To overcome the complications associated with this severe pathogenic virus, we have designed a multi-epitope subunit vaccine using computational immunology strategies. Identification of structural and nonstructural proteins of Nipah virus assisted in the vaccine designing. The selected proteins are known to be involved in the survival of the virus. The antigenic binders (B-cell, HTL, and CTL) from the selected proteins were prognosticated. These antigenic binders will be able to generate the humoral as well as cell-mediated immunity. All the epitopes were united with the help of suitable linkers and with an adjuvant at the N-terminal of the vaccine, for the enhancement of immunogenicity. The physiological characterization, along with antigenicity and allergenicity of the designed vaccine candidates, was estimated. The 3D structure prediction and its validation were performed. The validated vaccine model was then docked and simulated with the TLR-3 receptor to check the stability of the docked complex. This next-generation approach will provide a new vision for the development of a high immunogenic vaccine against the NiV.

Identification of novel small molecules against GSK3β for Alzheimer's disease using chemoinformatics approach

Alzheimer's disease is a rapidly increasing neurodegenerative disease. It is a multifactorial disease and also a global threat. Several enzymes are implicated in the disease in which Glycogen Synthase Kinase 3 beta is a key enzyme to increase the disease progression by the hyperphosphorylation of the tau protein. We have used an integrative chemoinformatics and pharmacokinetics approach for the identification of novel small molecules. We have retrieved a subset from the ZINC database (n = 5,36,709) and screened against GSK3β in four steps. From here top 298 potent compounds were selected and employed for their pharmacokinetics analysis. We had seen that 29 compounds showed the key characteristics to be a novel drug candidate therefore, all these compounds were employed for redocking studies using Autodock Vina and Autodock. This analysis revealed that four compounds were showing good binding affinity. All these four compounds were employed for MDS analysis of 100 ns From here using a bunch of MD analyses we have found that out of four compounds GSK3β-ZINC21011059 and GSK3β-ZINC21011066 act as a stable protein-ligand complex. Therefore we proposed ZINC21011059 and ZINC21011066 can serve as a novel compounds against GSK3β and predicted scaffold can be used for further optimization towards the improvement of isoform selectivity, and warranting further investigations towards their in vitro and in vivo validation of the bioactivity.

Mycobacterium tuberculosis nucleoside diphosphate kinase shows interaction with putative ATP binding cassette (ABC) transporter, Rv1273c

Protein-protein interactions are crucial for all biological processes. Compiling this network provides many new insights into protein function and gives directions for the development of new drugs targeted to the pathogen. Mycobacterium tuberculosis Nucleoside diphosphate kinase (Mtb Ndk) has been reported to promote survival of mycobacterium within the macrophage and contribute significantly to mycobacterium virulence. Hence, the present study was aimed to identify and characterize the interacting partner for Ndk. The in vitro experiments, pull down and far western blotting have demonstrated that Mtb Ndk interacts with Rv1273c, a probable drug ABC transporter ATP-binding protein annotated to export drugs across the membrane. This observation was further confirmed by molecular docking and dynamic simulations studies. The homology model of Rv1273c was constructed and docked with Mtb Ndk for protein-protein interaction analysis. The critical residues involved at interface of Rv1273c-Ndk interaction were identified. MDS and Principal Component analysis carried out for conformational feasibility and stability concluded that the complex between the two proteins is more stable as compared to apo proteins. Our findings would be expected to improve the dissection of protein-protein interaction network and significantly advance our understanding of tuberculosis infection.Communicated by Ramaswamy H. Sarma.

Applications of Ensemble Docking in Potential Inhibitor Screening for Mycobacterium tuberculosis Isocitrate Lyase Using a Local Plant Database

Isocitrate lyase (ICL) is a persistent factor for the survival of dormant stage Mycobacterium tuberculosis (MTB), thus a potential drug target for tuberculosis treatment. In this work, ensemble docking approach was used to screen for potential inhibitors of ICL. The ensemble conformations of ICL active site were obtained from molecular dynamics simulation on three dimer form systems, namely the apo ICL, ICL in complex with metabolites (glyoxylate and succinate), and ICL in complex with substrate (isocitrate). Together with the ensemble conformations and the X-ray crystal structures, 22 structures were used for the screening against Malaysian Natural Compound Database (NADI). The top 10 compounds for each ensemble conformation were selected. The number of compounds was then further narrowed down to 22 compounds that were within the Lipinski's Rule of Five for drug-likeliness and were also docked into more than one ensemble conformation. Theses 22 compounds were furthered evaluate using whole cell assay. Some compounds were not commercially available; therefore, plant crude extracts were used for the whole cell assay. Compared to itaconate (the known inhibitor of ICL), crude extracts from Manilkara zapota, Morinda citrifolia, Vitex negundo, and Momordica charantia showed some inhibition activity. The MIC/MBC value were 12.5/25, 12.5/25, 0.78/1.6, and 0.39/1.6 mg/mL, respectively. This work could serve as a preliminary study in order to narrow the scope for high throughput screening in the future.

Virtual screening, pharmacokinetics, molecular dynamics and binding free energy analysis for small natural molecules against cyclin-dependent kinase 5 for Alzheimer's disease

Alzheimer's disease (AD) is a progressive neurodegenerative disorder and characterized by brain cell death, memory loss and is the most common form of dementia. Although AD has devastating effects, however, drugs which can treat the AD remain limited. The cyclin-dependent kinase 5 (CDK5) has been recognized as being involved in the pathological hyperphosphorylation of tau protein, which leads to the formation of neurofibrillary tangles (NFTs). We utilized the structure-based virtual screening (SBVS) approach to find the potential inhibitors against HsCDK5. The natural compound subset from the ZINC database (n = 167,741) was retrieved and screened by using SBVS method. From here, we have predicted 297 potent inhibitors. These 297 compounds were evaluated through their pharmacokinetic properties by ADMET (absorption, distribution, metabolism, elimination/excretion and toxicity) descriptors. Finally, 17 compounds were selected and used for re-docking. After the refinement by molecular docking and by using drug-likeness analysis, we have identified four potential inhibitors (ZINC85877721, ZINC96114862, ZINC96115616 and ZINC96116231). All these four ligands were employed for 100 ns MDS study. From the root mean square deviation (RMSD), root mean square fluctuation (RMSF), Rg, number of hydrogen bonds, solvent accessible surface area (SASA), principal component analysis (PCA) and binding free energy analysis we have found that out of four inhibitors ZINC85877721 and ZINC96116231 showed good binding free energy of -198.84 and -159.32 kJ.mol^-1, respectively, and also good in other structural analyses. Both compounds displayed excellent pharmacological and structural properties to be the drug candidates. Collectively, these findings recommend that two compounds have great potential to be a promising agent against AD to reduce the CDK5 induced hyperphosphorylation and could be considered as therapeutic agents for the AD.Communicated by Ramaswamy H. Sarma.

In silico identification of mimicking molecule(s) triggering von Willebrand factor in human: a molecular drug target for regulating coagulation pathway

Blood coagulation is a complex and dynamic process wherein the body activates its emergency mechanism to stop bleeding and wound healing via the interactions of prothrombotic and antithrombotic agents. von Willebrand factor (VWF) is a complex glycoprotein and initial component of the hemostasis pathway which serves a multipurpose role in blood coagulation process. There are reports of various plants that contain several bioactive compounds possessing properties of inducing blood coagulation directly or indirectly. In the present study, efforts have been made to identify bioactive compounds that may play a significant role in regulation of the coagulation cascade by accelerating VWF and thus enhance the hemostasis process. An antidiuretic peptide drug, Desmopressin, works on VWF and releases them in circulation. Forty-seven compounds from different plant sources were screened through molecular docking, out of which two compounds, Emodin and Peruvianoside II, showed more binding affinity than the reference drug Desmopressin. Emodin and Peruvianoside II showed binding energies -7.2 and -7.0 kcal/mol, respectively, when docked with VWF, whereas Desmopressin displayed less binding energy (-6.9 kcal/mol). Emodin belongs to anthraquinone from Rumex hastasus and Peruvianoside II belongs to flavanone glycosides from Thevetia peruviana. The mimicking potential of top identified molecules with respect to the drug was confirmed through simulation analysis. Besides, the molecular dynamics simulation (MDS) study (for 20 ns) showed that the Peruvianoside II protein complex was energetically more stable than Emodin protein complex. Based on the results, Peruvianoside II possesses great potential and thus may be considered for development of drugs for hemostasis.

Structural and energetic understanding of novel natural inhibitors of Mycobacterium tuberculosis malate synthase

Persistent infection by Mycobacterium tuberculosis requires the glyoxylate shunt. This is a bypass to the tricarboxylic acid cycle in which isocitrate lyase (ICL) and malate synthase (MS) catalyze the net incorporation of carbon during mycobacterial growth on acetate or fatty acids as the primary carbon source. To identify a potential antitubercular compound, we performed a structure-based screening of natural compounds from the ZINC database (n = 1 67 740) against the M tuberculosis MS (MtbMS) structure. The ligands were screened against MtbMS, and 354 ligands were found to have better docking score. These compounds were assessed for Lipinski and absorption, distribution, metabolism, excretion, and toxicity prediction where 15 compounds were found to fit well for redocking studies. After refinement by molecular docking and drug-likeness analysis, four potential inhibitors (ZINC1483899, ZINC1754310, ZINC2269664, and ZINC15729522) were identified. These four ligands with phenyl-diketo acid were further subjected to molecular dynamics simulation to compare the dynamics and stability of the protein structure after ligand binding. The binding energy analysis was calculated to determine the intermolecular interactions. Our results suggested that the four compounds had a binding free energy of -201.96, -242.02, -187.03, and -169.02 kJ·mol^-1 , for compounds with IDs ZINC1483899, ZINC1754310, ZINC2269664, and ZINC15729522, respectively. We concluded that two compounds (ZINC1483899 and ZINC1754310) displayed considerable structural and pharmacological properties and could be probable drug candidates to fight against M tuberculosis parasites.

In silico screening of deleterious single nucleotide polymorphisms (SNPs) and molecular dynamics simulation of disease associated mutations in gene responsible for oculocutaneous albinism type 6 (OCA 6) disorder

Solute carrier family 24 member 5 (SLC24A5) is a gene that is associated with oculocutaneous albinism type 6 (OCA6) disorder and is involved in skin and hair pigmentation. It is involved in the maturation of melanosomes and melanin synthesis. SLC24A5 gene is located in the chromosomal position of 15q21.1. The present study involves the use of computational techniques in order to obtain a detailed picture of the most probable mutations that are associated with SLC24A5. From the observed result it was found that the mutation S145F is most deleterious and disease associated is predicted using several bioinformatics tools. The 3-D structures of native and mutant (S145F) were modeled in order to understand protein functionality using ab initio Robetta server. The modeled structure validation was done with ERRAT, Verify-3D, Procheck and RAMPAGE Ramachandran plot analysis. The most validated structure undergoes molecular dynamics simulations (MDS) study to understand the structural and functional behaviour of the native and mutant proteins. The MDS result showed the more flexibility in the native SLC24A5 structure. Due to mutation in the SLC24A5 protein structure it became more rigid and might disturb the conformational changes and glycosylation function of protein structure and might play role in inducing the OCA6. This study provides a significant insight into the underlying molecular mechanism involved in albinism associated with OCA6. It further helps scientists to develop a drug therapy against OCA 6 disease. Communicated by Ramaswamy H. Sarma.

Structure-function studies of the asparaginyl-tRNA synthetase from Fasciola gigantica: understanding the role of catalytic and non-catalytic domains

The asparaginyl-tRNA synthetase (NRS) catalyzes the attachment of asparagine to its cognate tRNA during translation. NRS first catalyzes the binding of Asn and ATP to form the NRS-asparaginyl adenylate complex, followed by the esterification of Asn to its tRNA. We investigated the role of constituent domains in regulating the structure and activity of Fasciola gigantica NRS (FgNRS). We cloned the full-length FgNRS, along with its various truncated forms, expressed, and purified the corresponding proteins. Size exclusion chromatography indicated a role of the anticodon-binding domain (ABD) of FgNRS in protein dimerization. The N-terminal domain (NTD) was not essential for cognate tRNA binding, and the hinge region between the ABD and the C-terminal domain (CTD) was crucial for regulating the enzymatic activity. Molecular docking and fluorescence quenching experiments elucidated the binding affinities of the substrates to various domains. The molecular dynamics simulation of the modeled protein showed the presence of an unstructured region between the NTD and ABD that exhibited a large number of conformations over time, and further analysis indicated this region to be intrinsically disordered. The present study provides information on the structural and functional regulation, protein-substrate(s) interactions and dynamics, and the role of non-catalytic domains in regulating the activity of FgNRS.

Combinatorial screening algorithm to engineer multiepitope subunit vaccine targeting human T-lymphotropic virus-1 infection

Human T-lymphotropic virus (HTLV), the first human retrovirus has been discovered which is known to cause the age-old assassinating disease HTLV-1 associated myelopathy. Cancer caused by this virus is adult T cell leukemia/lymphoma which targets 10-20 million throughout the world. The effect of this virus extends to the fact that it causes chronic disease to the spinal cord resulting in loss of sensation and further causes blood cancer. So, to overcome the complications, we designed a subunit vaccine by the assimilation of B-cell, cytotoxic T-lymphocyte , and helper T-lymphocyte epitopes. The epitopes were joined together along with adjuvant and linkers and a vaccine was fabricated which was further subjected to 3D modeling. The physiochemical properties, allergenicity, and antigenicity were evaluated. Molecular docking and dynamics were performed with the obtained 3D model against toll like receptor (TLR-3) immune receptor. Lastly, in silico cloning was performed to ensure the expression of the designed vaccine in pET28a(+) expression vector. The future prospects of the study entailed the in vitro and in vivo experimental analysis for evaluating the immune response of the designed vaccine construct.

Targeting Nucleotide Binding Domain of Multidrug Resistance-associated Protein-1 (MRP1) for the Reversal of Multi Drug Resistance in Cancer

Multidrug resistance (MDR) is the major cause, by which cancer cells expel the drugs out, developing a challenge against the current chemotherapeutic drugs regime. This mechanism is attributed to the over expression of ABC transporters like MRP1 on the surface of cells. Since nucleotide binding domains (NBD) of ABC transporters are the site of ATP binding and hydrolysis, thereby in this study we have targeted NBD1 of MRP1using molecular docking and molecular dynamic simulations (MDS). The compounds present in the FDA approved library were docked against NBD1 of the human multidrug resistance associated protein 1 (PDB ID: 2CBZ). For the docking studies, Standard Precision and Extra Precision methods were employed. After the EP docking studies, ligands showed an extremely low docking score that was indicative of very high binding affinity of the ligands to the NBD. Apart from the low docking score, another short listing criterion in simulation studies was the interaction of incoming ligand with the desired conserved residues of NDB involved in ATP binding and hydrolysis. Based on these measures, potassium citrate (DB09125) and technetium Tc-99m medronate (DB09138) were chosen and subjected to 100 ns simulation studies. From the MDS study we concluded that between these two compounds, potassium citrate is a better candidate for targeting MRP1.

Drug targeted virtual screening and molecular dynamics of LipU protein of Mycobacterium tuberculosis and Mycobacterium leprae

The lipolytic protein LipU was conserved in mycobacterium sp. including M. tuberculosis (MTB LipU) and M. leprae (MLP LipU). The MTB LipU was identified in extracellular fraction and was reported to be essential for the survival of mycobacterium. Therefore to address the problem of drug resistance in pathogen, LipU was selected as a drug target and the viability of finding out some FDA approved drugs as LipU inhibitors in both the cases was explored. Three-dimensional (3D) model structures of MTB LipU and MLP LipU were generated and stabilized through molecular dynamics (MD). FDA approved drugs were screened against these proteins. The result showed that the top-scoring compounds for MTB LipU were Diosmin, Acarbose and Ouabain with the Glide XP score of -12.8, -11.9 and -11.7 kcal/mol, respectively, whereas for MLP LipU protein, Digoxin (-9.2 kcal/mol), Indinavir (-8.2 kcal/mol) and Travoprost (-8.2 kcal/mol) showed highest affinity. These drugs remained bound in the active site pocket of MTB LipU and MLP LipU structure and interaction grew stronger after dynamics. RMSD, RMSF and Rg were found to be persistent throughout the simulation period. Hydrogen bonds along with large number of hydrophobic interactions stabilized the complex structures. Binding free energies obtained through Prime/MM-GBSA were found in the significant range from -63.85 kcal/mol to -34.57 kcal/mol for MTB LipU and -71.33 kcal/mol to -23.91 kcal/mol for MLP LipU. The report suggested high probability of these drugs to demolish the LipU activity and could be probable drug candidates to combat TB and leprosy disease.

Comprehensive analysis of the catalytic and structural properties of a mu-class glutathione s-transferase from Fasciola gigantica

Glutathione S‒transferases (GSTs) play an important role in the detoxification of xenobiotics. They catalyze the nucleophilic addition of glutathione (GSH) to nonpolar compounds, rendering the products water-soluble. In the present study, we investigated the catalytic and structural properties of a mu-class GST from Fasciola gigantica (FgGST1). The purified recombinant FgGST1 formed a homodimer composed of 25 kDa subunit. Kinetic analysis revealed that FgGST1 displays broad substrate specificity and shows high GSH conjugation activity toward 1-chloro-2,4-dinitrobenzene, 4-nitroquinoline-1-oxide, and trans-4-phenyl-3-butene-2-one and peroxidase activity towards trans-2-nonenal and hexa-2,4-dienal. The FgGST1 was highly sensitive to inhibition by cibacron blue. The cofactor (GSH) and inhibitor (cibacron blue) were docked, and binding sites were identified. The molecular dynamics studies and principal component analysis indicated the stability of the systems and the collective motions, respectively. Unfolding studies suggest that FgGST1 is a highly cooperative molecule because, during GdnHCl-induced denaturation, a simultaneous unfolding of the protein without stabilization of any partially folded intermediate is observed. The protein is stabilized with a conformational free energy of about 10 ± 0.3 kcal mol-1. Additionally, the presence of conserved Pro-53 and structural motifs such as N-capping box and hydrophobic staple, further aided in the stability and proper folding of FgGST1.