Rohit Shukla

Distant Phe345 mutation compromises the stability and activity of Mycobacterium tuberculosis isocitrate lyase by modulating its structural flexibility

Isocitrate lyase (ICL), a potential anti-tubercular drug target, catalyzes the first step of the glyoxylate shunt. In the present investigation, we studied the conformational flexibility of MtbICL to better understand its stability and catalytic activity. Our biochemical results showed that a point mutation at Phe345, which is topologically distant (>10 Å) to the active site signature sequence (189KKCGH193), completely abolishes the activity of the enzyme. In depth computational analyses were carried out for understanding the structural alterations using molecular dynamics, time-dependent secondary structure and principal component analysis. The results showed that the mutated residue increased the structural flexibility and induced conformational changes near the active site (residues 170–210) and in the C-terminal lid region (residues 411–428). Both these regions are involved in the catalytic activity of MtbICL. Upon mutation, the residual mobility of the enzyme increased, resulting in a decrease in the stability, which was confirmed by the lower free energy of stabilization in the mutant enzyme suggesting the destabilization in the structure. Our results have both biological importance and chemical novelty. It reveals internal dynamics of the enzyme structure and also suggests that regions other than the active site should be exploited for targeting MtbICL inhibition and development of novel anti-tuberculosis compounds.

Structure-Based Screening And Molecular Dynamics Simulations Offer Novel Natural Compounds As Potential Inhibitors Of Mycobacterium Tuberculosis Isocitrate Lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with −96.462, −143.549, and −122.526 kJ mol−1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.

A combined biochemical and computational studies of the rho-class glutathione s-transferase sll1545 of Synechocystis PCC 6803.

Peroxides are one of the most important radicals that cause oxidative stress. Certain Glutathione S-transferases (GSTs) have been reported to show peroxidase activity. We report a novel peroxidase activity of Synechocystis GST- sll1545. The recombinant protein was purified to homogeneity and characterized. Low Km (0.109μM) and high Vmax (0.663μmolmin-1) values suggest a high preference of sll1545 for cumenehydroperoxide. Disc inhibition assay confirmed the ability of the enzyme to protect cells against peroxide-induced damage. sll1545 has very low sequence and structural similarity with theta and alpha class GSTs that exhibit glutathione-dependent peroxidase activity. Recent data from our laboratory shows that sll1545 is also strongly active against dichloroacetate (DCA), which is a characteristic of zeta class GST. Interestingly, sll1545 shows less than 20% sequence identity with zeta class GST. Molecular dynamic simulation results show that sll1545 was much more structurally different from alpha/theta classes. Our results suggest that sll1545 shows structural variation from zeta, theta/alpha classes of GSTs but have related enzymatic activity. Phylogenetic analysis reveal that sll1545 is evolutionally very distinct from the known GSTs. Overall, the data suggest that Synechocystis sll1545 does not belong to any known GST class and represent a novel GST class, which we have named rho.

Alterations in conformational topology and interaction dynamics caused by L418A mutation leads to activity loss of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis isocitrate lyase (MtbICL) is a key enzyme of the glyoxylate cycle that catalyzes the cleavage of isocitrate to succinate and glyoxylate and is a potential antituberculosis drug target. The aim of this research was to explore the structural alterations induced by L418A point mutation that caused the loss of enzyme activity. In-depth structural analyses were carried out for understanding the influence of L418A mutation using techniques, viz. molecular dynamics, principal component analysis, time-dependent secondary structure, residue interaction network and molecular docking. Since L418A mutation site is structurally far from the active site, it cannot influence the binding of the substrate directly. Our results showed that collective motions, residual mobility, and flexibility of the enzyme increased upon mutation. The mutated residue changed the global conformational dynamics of the system along with the residue-residue interaction network, leading to a loss of the enzyme activity. The docking results suggest that L418A mutation influenced the binding interactions of the substrate with several residues in the active site of MtbICL. This study provides information on the structural dynamics of MtbICL and highlights the importance of residue level interactions in the protein. Thus, our results may provide significant guidance to the scientific community engaged in designing potent inhibitors targeting MtbICL.

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening, molecular dynamics simulation, and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol−1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.

UDP-N-Acetylglucosamine enolpyruvyl transferase (MurA) of Acinetobacter baumannii (AbMurA): Structural and functional properties.

Peptidoglycan (PG) is the key component of the bacterial cell wall. The enzyme UDP-N-Acetylglucosamine Enolpyruvyl Transferase (MurA) catalyzes the transfer of enolpyruvate from phosphoenolpyruvate (PEP) to uridinediphospho-N-acetylglucosamine (UNAG), which is the first committed step of PG biosynthesis. Here, we present the biochemical and structural features of the MurA enzyme of the opportunistic pathogen Acinetobacter baumannii (AbMurA). The recombinant AbMurA exists as a monomer in solution and shows optimal activity at pH 7.5 and 37°C. The Km for UDP-N-acetylglucosamine was 1.062±0.09mM and for PEP was 1.806±0.23mM. The relative enzymatic activity was inhibited ∼3 fold in the presence of 50mM fosfomycin (FFQ). Superimposition of the AbMurA model with E. coli demonstrated key structural similarity in the FFQ-binding site. AbMurA also has a surface loop that contains the active site Cys116 that interact with FFQ. Sequence analysis indicates the presence of the five conserved amino acids, i.e., K22, C116, D306, D370 and L371, required for the functional activity like other MurA enzymes from different bacteria. MurA enzymes are indispensable for cell integrity and their lack of counterparts in eukaryotes suggests them to be a promising drug target.

In Silico Identification of Mimicking Molecules as Defense Inducers Triggering Jasmonic Acid Mediated Immunity against Alternaria Blight Disease in Brassica Species

Alternaria brassicae and Alternaria brassicicola are two major phytopathogenic fungi which cause Alternaria blight, a recalcitrant disease on Brassica crops throughout the world, which is highly destructive and responsible for significant yield losses. Since no resistant source is available against Alternaria blight, therefore, efforts have been made in the present study to identify defense inducer molecules which can induce jasmonic acid (JA) mediated defense against the disease. It is believed that JA triggered defense response will prevent necrotrophic mode of colonization of Alternaria brassicae fungus. The JA receptor, COI1 is one of the potential targets for triggering JA mediated immunity through interaction with JA signal. In the present study, few mimicking compounds more efficient than naturally occurring JA in terms of interaction with COI1 were identified through virtual screening and molecular dynamics simulation studies. A high quality structural model of COI1 was developed using the protein sequence of Brassica rapa. This was followed by virtual screening of 767 analogs of JA from ZINC database for interaction with COI1. Two analogs viz. ZINC27640214 and ZINC43772052 showed more binding affinity with COI1 as compared to naturally occurring JA. Molecular dynamics simulation of COI1 and COI1-JA complex, as well as best screened interacting structural analogs of JA with COI1 was done for 50 ns to validate the stability of system. It was found that ZINC27640214 possesses efficient, stable, and good cell permeability properties. Based on the obtained results and its physicochemical properties, it is capable of mimicking JA signaling and may be used as defense inducers for triggering JA mediated resistance against Alternaria blight, only after further validation through field trials.

Ameliorated de novo transcriptome assembly using Illumina paired end sequence data with Trinity Assembler.

Advent of Next Generation Sequencing has led to possibilities of de novo transcriptome assembly of organisms without availability of complete genome sequence. Among various sequencing platforms available, Illumina is the most widely used platform based on data quality, quantity and cost. Various de novo transcriptome assemblers are also available today for construction of de novo transcriptome. In this study, we aimed at obtaining an ameliorated de novo transcriptome assembly with sequence reads obtained from Illumina platform and assembled using Trinity Assembler. We found that, primary transcriptome assembly obtained as a result of Trinity can be ameliorated on the basis of transcript length, coverage, and depth and protein homology. Our approach to ameliorate is reproducible and could enhance the sensitivity and specificity of the assembled transcriptome which could be critical for validation of the assembled transcripts and for planning various downstream biological assays.

Structural insights into natural compounds as inhibitors of Fasciola gigantica thioredoxin glutathione reductase

Fascioliasis is caused by the helminth parasites of genus Fasciola. Thioredoxin glutathione reductase (TGR) is an important enzyme in parasitic helminths and plays an indispensable role in its redox biology. In the present study, we conducted a structure‐based virtual screening of natural compounds against the Fasciola gigantica TGR (FgTGR). The compounds were docked against FgTGR in four sequential docking modes. The screened ligands were further assessed for Lipinski and ADMET prediction so as to evaluate drug proficiency and likeness property. After refinement, three potential inhibitors were identified that were subjected to 50 ns molecular dynamics simulation and free energy binding analyses to evaluate the dynamics of protein‐ligand interaction and the stability of the complexes. Key residues involved in the interaction of the selected ligands were also determined. The results suggested that three top hits had a negative binding energy greater than GSSG (−91.479 KJ · mol−1), having −152.657, −141.219, and −92.931 kJ · mol−1 for compounds with IDs ZINC85878789, ZINC85879991, and ZINC36369921, respectively. Further analysis showed that the compound ZINC85878789 and ZINC85879991 displayed substantial pharmacological and structural properties to be a drug candidate. Thus, the present study might prove useful for the future design of new derivatives with higher potency and specificity.

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with −259.41, −110.09, −188.25, −163.30, −202.10, and −158.79 kJ mol−1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.