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.

[1]  Md. Sohail Akhtar,et al.  Guanidinium chloride- and urea-induced unfolding of the dimeric enzyme glucose oxidase. , 2002, Biochemistry.

[2]  C. Pace Determination and analysis of urea and guanidine hydrochloride denaturation curves. , 1986, Methods in enzymology.

[3]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[4]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[5]  J. Graham,et al.  Identification of Mycobacterium tuberculosis RNAs synthesized in response to phagocytosis by human macrophages by selective capture of transcribed sequences (SCOTS). , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[6]  C. Nathan,et al.  Isocitrate lyase mediates broad antibiotic tolerance in Mycobacterium tuberculosis , 2014, Nature Communications.

[7]  Md. Sohail Akhtar,et al.  Isocitrate lyase of Mycobacterium tuberculosis is inhibited by quercetin through binding at N-terminus. , 2015, International journal of biological macromolecules.

[8]  T. Tripathi,et al.  UDP-N-Acetylglucosamine enolpyruvyl transferase (MurA) of Acinetobacter baumannii (AbMurA): Structural and functional properties. , 2017, International journal of biological macromolecules.

[9]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[10]  Ranjeet Kumar,et al.  Mycobacterium tuberculosis isocitrate lyase (MtbIcl): Role of divalent cations in modulation of functional and structural properties , 2008, Proteins.

[11]  D. Boehr,et al.  An NMR perspective on enzyme dynamics. , 2006, Chemical reviews.

[12]  David G. Russell,et al.  Characterization of Activity and Expression of Isocitrate Lyase in Mycobacterium avium andMycobacterium tuberculosis , 1999, Journal of bacteriology.

[13]  R. Hartkoorn,et al.  Tuberculosis drugs: new candidates and how to find more. , 2011, Future microbiology.

[14]  Sujata Sharma,et al.  Structure of isocitrate lyase, a persistence factor of Mycobacterium tuberculosis , 2000, Nature Structural Biology.

[15]  A. Singh,et al.  Insight into the structural flexibility and function of Mycobacterium tuberculosis isocitrate lyase. , 2015, Biochimie.

[16]  L. Longo,et al.  A logical OR redundancy within the Asx-Pro-Asx-Gly type I beta-turn motif. , 2008, Journal of molecular biology.

[17]  C Ratledge,et al.  Use of carbon sources for lipid biosynthesis in Mycobacterium leprae: a comparison with other pathogenic mycobacteria. , 1988, Journal of general microbiology.

[18]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[19]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[20]  G. Krishnamoorthy Motional dynamics in proteins and nucleic acids control their function : revelation by time-domain fluorescence , 2012 .

[21]  Jihun Lee,et al.  Spackling the crack: stabilizing human fibroblast growth factor-1 by targeting the N and C terminus beta-strand interactions. , 2007, Journal of molecular biology.

[22]  Lynn Rasmussen,et al.  High-throughput screening for inhibitors of Mycobacterium tuberculosis H37Rv. , 2009, Tuberculosis.

[23]  B. McFadden,et al.  Lysine 194 Is Functional in Isocitrate Lyase from Escherichia coli , 1997, Current Microbiology.

[24]  Jihun Lee,et al.  Redesigning symmetry‐related “mini‐core” regions of FGF‐1 to increase primary structure symmetry: Thermodynamic and functional consequences of structural symmetry , 2005, Protein science : a publication of the Protein Society.

[25]  T. Creighton,et al.  Experimental studies of protein folding and unfolding. , 1978, Progress in biophysics and molecular biology.

[26]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[27]  B. McFadden,et al.  Site-directed mutagenesis of lysine 193 in Escherichia coli isocitrate lyase by use of unique restriction enzyme site elimination , 1993, Journal of bacteriology.

[28]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[29]  D Eisenberg,et al.  3D domain swapping: A mechanism for oligomer assembly , 1995, Protein science : a publication of the Protein Society.

[30]  T. Tripathi,et al.  A combined biochemical and computational studies of the rho-class glutathione s-transferase sll1545 of Synechocystis PCC 6803. , 2017, International journal of biological macromolecules.

[31]  Yang Liu,et al.  Transcriptional Adaptation of Mycobacterium tuberculosis within Macrophages , 2003, The Journal of experimental medicine.

[32]  Antje Wolf,et al.  Principal component and clustering analysis on molecular dynamics data of the ribosomal L11·23S subdomain , 2012, Journal of Molecular Modeling.

[33]  B. McFadden,et al.  The importance of four histidine residues in isocitrate lyase from Escherichia coli , 1994, Journal of bacteriology.

[34]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[35]  Peter Güntert,et al.  Spatial elucidation of motion in proteins by ensemble-based structure calculation using exact NOEs , 2012, Nature Structural &Molecular Biology.

[36]  A. Sharma,et al.  Isocitrate lyase: a potential target for anti-tubercular drugs. , 2013, Recent patents on inflammation & allergy drug discovery.

[37]  T. Tripathi CALCULATION OF THERMODYNAMIC PARAMETERS OF PROTEIN UNFOLDING USING FAR-ULTRAVIOLET CIRCULAR DICHROISM , 2013 .

[38]  B. McFadden,et al.  Comparative structure, function and regulation of isocitrate lyase, an important assimilatory enzyme. , 1990, Comparative biochemistry and physiology. B, Comparative biochemistry.

[39]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[40]  Charles C. David,et al.  Principal component analysis: a method for determining the essential dynamics of proteins. , 2014, Methods in molecular biology.