Molecular modeling studies of Fatty acyl-CoA synthetase (FadD13) from Mycobacterium tuberculosis—a potential target for the development of antitubercular drugs

Tuberculosis (TB) is a global health problem and the situation has become more precarious due to the advent of HIV infections and continuous rise in the number of multi-drug resistant strains of Mycobacterium tuberculosis (M. tb). Biochemical studies on Fatty Acyl-CoA Synthetases (FadD13), one of the gene products of mymA operon, have provided insights into the involvement of this protein in the activation of fatty acids. Due to non-availability of the crystal structure of FadD13, we have employed in silico approaches to resolve and characterize the structure of this important protein of M. tb. A three dimensional model of M. tb FadD13 was predicted by a de novo structure prediction server that integrates fragment assembly with SimFold energy function. With the aid of molecular mechanics and dynamics methods, the final model was obtained and assessed subsequently for global and local accuracy by various assessment programs. With this model, a flexible docking study with the substrates was performed. Results of ligand interactions with key amino acids in the binding site are also summarized. The molecular model for the M. tb FadD13 obtained sheds light on the topographical features of the binding pocket of the protein and provides atomic insight into the possible modes of substrate recognition. The three-dimensional model of FadD13 presented here would be helpful in guiding both enzymatic studies as well as design of specific inhibitors.

[1]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

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

[3]  Amit Singh,et al.  Requirement of the mymA Operon for Appropriate Cell Wall Ultrastructure and Persistence of Mycobacterium tuberculosis in the Spleens of Guinea Pigs , 2005, Journal of bacteriology.

[4]  Jie Liang,et al.  CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues , 2006, Nucleic Acids Res..

[5]  Richard M. Jackson,et al.  Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites , 2005, Bioinform..

[6]  A. Elofsson,et al.  Can correct protein models be identified? , 2003, Protein science : a publication of the Protein Society.

[7]  P. Groot,et al.  Fatty acid activation: specificity, localization, and function. , 1976, Advances in lipid research.

[8]  Tetsuya Hori,et al.  Structural Basis of the Substrate-specific Two-step Catalysis of Long Chain Fatty Acyl-CoA Synthetase Dimer* , 2004, Journal of Biological Chemistry.

[9]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[10]  S. Anderson,et al.  The role of lysine 529, a conserved residue of the acyl-adenylate-forming enzyme superfamily, in firefly luciferase. , 2000, Biochemistry.

[11]  B. Bloom,et al.  Tuberculosis Pathogenesis, Protection, and Control , 1994 .

[12]  J M Thornton,et al.  LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. , 1995, Protein engineering.

[13]  Alexander V. Diemand,et al.  iMolTalk: an interactive, internet-based protein structure analysis server , 2004, Nucleic Acids Res..

[14]  Shoji Takada,et al.  SimFold energy function for de novo protein structure prediction: Consensus with Rosetta , 2005, Proteins.

[15]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[16]  K. Reynolds,et al.  Antibacterial targets in fatty acid biosynthesis. , 2007, Current opinion in microbiology.

[17]  A. Horswill,et al.  Characterization of the propionyl-CoA synthetase (PrpE) enzyme of Salmonella enterica: residue Lys592 is required for propionyl-AMP synthesis. , 2002, Biochemistry.

[18]  Matthew P. Repasky,et al.  Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. , 2006, Journal of medicinal chemistry.

[19]  A W Munro,et al.  The TB structural genomics consortium: a resource for Mycobacterium tuberculosis biology. , 2003, Tuberculosis.

[20]  K. Mdluli,et al.  Novel targets for tuberculosis drug discovery. , 2006, Current opinion in pharmacology.

[21]  B. Barrell,et al.  Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence , 1998, Nature.

[22]  Michael J E Sternberg,et al.  Exploring the extremes of sequence/structure space with ensemble fold recognition in the program Phyre , 2008, Proteins.

[23]  D. Eisenberg,et al.  Assessment of protein models with three-dimensional profiles , 1992, Nature.

[24]  Ben M. Webb,et al.  Comparative Protein Structure Modeling Using Modeller , 2006, Current protocols in bioinformatics.

[25]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[26]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[27]  William L. Jorgensen,et al.  Temperature and size dependence for Monte Carlo simulations of TIP4P water , 1985 .

[28]  Christophe G. Lambert,et al.  ESyPred3D: Prediction of proteins 3D structures , 2002, Bioinform..

[29]  Lloyd N. Friedman,et al.  Diagnostic standards and classification of tuberculosis. , 1991, The American review of respiratory disease.

[30]  David Eisenberg,et al.  Structural genomics of Mycobacterium tuberculosis: a preliminary report of progress at UCLA. , 2003, Biophysical chemistry.

[31]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[32]  L. Holm,et al.  The Pfam protein families database , 2005, Nucleic Acids Res..

[33]  Anil K. Tyagi,et al.  Dissecting the Role of Critical Residues and Substrate Preference of a Fatty Acyl-CoA Synthetase (FadD13) of Mycobacterium tuberculosis , 2009, PloS one.

[34]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[35]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[36]  Manfred J. Sippl,et al.  Thirty years of environmental health research--and growing. , 1996, Nucleic Acids Res..

[37]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[38]  E. Baker Structural genomics as an approach towards understanding the biology of tuberculosis , 2007, Journal of Structural and Functional Genomics.

[39]  T. Yeates,et al.  Verification of protein structures: Patterns of nonbonded atomic interactions , 1993, Protein science : a publication of the Protein Society.

[40]  Robert D. Finn,et al.  The Pfam protein families database , 2004, Nucleic Acids Res..

[41]  P. Black,et al.  Mutational Analysis of a Fatty Acyl-Coenzyme A Synthetase Signature Motif Identifies Seven Amino Acid Residues That Modulate Fatty Acid Substrate Specificity* , 1997, The Journal of Biological Chemistry.

[42]  Samiul Hasan,et al.  Prioritizing Genomic Drug Targets in Pathogens: Application to Mycobacterium tuberculosis , 2006, PLoS Comput. Biol..

[43]  Thomas M. Shinnick,et al.  Microarray Analysis of the Mycobacterium tuberculosis Transcriptional Response to the Acidic Conditions Found in Phagosomes , 2002, Journal of bacteriology.

[44]  J. Porter,et al.  The global burden of tuberculosis. , 1999 .

[45]  Manuel C. Peitsch,et al.  SWISS-MODEL: an automated protein homology-modeling server , 2003, Nucleic Acids Res..

[46]  S. Nosé A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[47]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[48]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[49]  H. Ropers,et al.  FACL4, encoding fatty acid-CoA ligase 4, is mutated in nonspecific X-linked mental retardation , 2002, Nature Genetics.

[50]  Christopher Dye,et al.  Global burden of tuberculosis , 1999 .

[51]  M. Daffé,et al.  The envelope layers of mycobacteria with reference to their pathogenicity. , 1998, Advances in microbial physiology.