Modeling holo-ACP:DH and holo-ACP:KR complexes of modular polyketide synthases: a docking and molecular dynamics study

BMC Structural Biology

Debasisa Mohanty¹  Swadha Anand¹ 

[1] Bioinformatics Center, National Institute of Immunology, Aruna Asaf Ali Marg, New Delhi, 110067, India

关键词: Ketoreductase domain;    Dehydratase domain;    Modular polyketide synthase;    Evolutionary conservation;    Substrate binding site;    Protein-protein interaction;    Protein-ligand docking;    Molecular dynamics;   

Others  :  1092236
DOI  :  10.1186/1472-6807-12-10

 received in 2011-11-18, accepted in 2012-05-28,  发布年份 2012

【 摘 要 】

Background

Modular polyketide synthases are multifunctional megasynthases which biosynthesize a variety of secondary metabolites using various combinations of dehydratase (DH), ketoreductase (KR) and enoyl-reductase (ER) domains. During the catalysis of various reductive steps these domains act on a substrate moiety which is covalently attached to the phosphopantetheine (P-pant) group of the holo-Acyl Carrier Protein (holo-ACP) domain, thus necessitating the formation of holo-ACP:DH and holo-ACP:KR complexes. Even though three dimensional structures are available for DH, KR and ACP domains, no structures are available for DH or KR domains in complex with ACP or substrate moieties. Since Ser of holo-ACP is covalently attached to a large phosphopantetheine group, obtaining complexes involving holo-ACP by standard protein-protein docking has been a difficult task.

Results

We have modeled the holo-ACP:DH and holo-ACP:KR complexes for identifying specific residues on DH and KR domains which are involved in interaction with ACP, phosphopantetheine and substrate moiety. A novel combination of protein-protein and protein-ligand docking has been used to first model complexes involving apo-ACP and then dock the phosphopantetheine and substrate moieties using covalent connectivity between ACP, phosphopantetheine and substrate moiety as constraints. The holo-ACP:DH and holo-ACP:KR complexes obtained from docking have been further refined by restraint free explicit solvent MD simulations to incorporate effects of ligand and receptor flexibilities. The results from 50 ns MD simulations reveal that substrate enters into a deep tunnel in DH domain while in case of KR domain the substrate binds a shallow surface exposed cavity. Interestingly, in case of DH domain the predicted binding site overlapped with the binding site in the inhibitor bound crystal structure of FabZ, the DH domain from E.Coli FAS. In case of KR domain, the substrate binding site identified by our simulations was in proximity of the known stereo-specificity determining residues.

Conclusions

We have modeled the holo-ACP:DH and holo-ACP:KR complexes and identified the specific residues on DH and KR domains which are involved in interaction with ACP, phosphopantetheine and substrate moiety. Analysis of the conservation profile of binding pocket residues in homologous sequences of DH and KR domains indicated that, these results can also be extrapolated to reductive domains of other modular PKS clusters.

【 授权许可】

   
2012 Anand and Mohanty; licensee BioMed Central Ltd.

【 图 表 】

Figure 1 .

Figure 2 .

Figure 3 .

Figure 4 .

Figure 5 .

Figure 6 .

Figure 7 .

Figure 8 .

Figure 9 .

Figure 10 .

Figure 11 .

Figure 12 .

Figure 13 .

【 参考文献 】

• [1]Fischbach MA, Walsh CT: Assembly-line enzymology for polyketide and nonribosomal peptide antibiotics: logic, machinery, and mechanisms. Chem Rev 2006, 106(8):3468-3496.
• [2]Ansari MZ, Yadav G, Gokhale RS, Mohanty D: NRPS-PKS: a knowledge-based resource for analysis of NRPS/PKS megasynthases. Nucleic Acids Res 2004, 32:405-413.
• [3]Dillon SC, Bateman A: The hotdog fold: wrapping up a superfamily of thioesterases and dehydratases. BMC Bioinformatics 2004, 5:109. BioMed Central Full Text
• [4]Starcevic A, Zucko J, Simunkovic J, Long PF, Cullum J, Hranueli D: ClustScan: an integrated program package for the semi-automatic annotation of modular biosynthetic gene clusters and in silico prediction of novel chemical structures. Nucleic Acids Res 2008, 36(21):6882-6892.
• [5]Keatinge-Clay AT: A tylosin ketoreductase reveals how chirality is determined in polyketides. Chem Biol 2007, 14(8):898-908.
• [6]Caffrey P: Conserved amino acid residues correlating with ketoreductase stereospecificity in modular polyketide synthases. Chembiochem 2003, 4(7):654-657.
• [7]Kimber MS, Martin F, Lu Y, Houston S, Vedadi M, Dharamsi A, Fiebig KM, Schmid M, Rock CO: The structure of (3R)-hydroxyacyl-acyl carrier protein dehydratase (FabZ) from Pseudomonas aeruginosa. J Biol Chem 2004, 279(50):52593-52602.
• [8]Kostrewa D, Winkler FK, Folkers G, Scapozza L, Perozzo R: The crystal structure of PfFabZ, the unique beta-hydroxyacyl-ACP dehydratase involved in fatty acid biosynthesis of plasmodium falciparum. Protein Sci 2005, 14(6):1570-1580.
• [9]Leesong M, Henderson BS, Gillig JR, Schwab JM, Smith JL: Structure of a dehydratase-isomerase from the bacterial pathway for biosynthesis of unsaturated fatty acids: two catalytic activities in one active site. Structure 1996, 4(3):253-264.
• [10]Korman TP, Tan YH, Wong J, Luo R, Tsai SC: Inhibition kinetics and emodin cocrystal structure of a type II polyketide ketoreductase. Biochemistry 2008, 47(7):1837-1847.
• [11]Korman TP, Hill JA, Vu TN, Tsai SC: Structural analysis of actinorhodin polyketide ketoreductase: cofactor binding and substrate specificity. Biochemistry 2004, 43(46):14529-14538.
• [12]Hadfield AT, Limpkin C, Teartasin W, Simpson TJ, Crosby J, Crump MP: The crystal structure of the actIII actinorhodin polyketide reductase: proposed mechanism for ACP and polyketide binding. Structure 2004, 12(10):1865-1875.
• [13]Keatinge-Clay AT, Maltby DA, Medzihradszky KF, Khosla C, Stroud RM: An antibiotic factory caught in action. Nat Struct Mol Biol 2004, 11(9):888-893.
• [14]Javidpour P, Korman TP, Shakya G, Tsai SC: Structural and biochemical analyses of regio- and stereospecificities observed in a type II polyketide ketoreductase. Biochemistry 2011, 50(21):4638-4649.
• [15]Javidpour P, Das A, Khosla C, Tsai SC: Structural and Biochemical Studies of the hedamycin type ii polyketide ketoreductase (HedKR): molecular basis of stereo- and regiospecificities. Biochemistry 2011, 50(34):7426-7439.
• [16]Zheng J, Keatinge-Clay AT: Structural and functional analysis of C2-type ketoreductases from modular polyketide synthases. J Mol Biol , 410(1):105-117.
• [17]Zheng J, Taylor CA, Piasecki SK, Keatinge-Clay AT: Structural and functional analysis of a-type ketoreductases from the amphotericin modular polyketide synthase. Structure , 18(8):913-922.
• [18]Akey DL, Razelun JR, Tehranisa J, Sherman DH, Gerwick WH, Smith JL: Crystal structures of dehydratase domains from the curacin polyketide biosynthetic pathway. Structure , 18(1):94-105.
• [19]Keatinge-Clay A: Crystal structure of the erythromycin polyketide synthase dehydratase. J Mol Biol 2008, 384(4):941-953.
• [20]Keatinge-Clay AT, Stroud RM: The structure of a ketoreductase determines the organization of the beta-carbon processing enzymes of modular polyketide synthases. Structure 2006, 14(4):737-748.
• [21]Maier T, Leibundgut M, Ban N: The crystal structure of a mammalian fatty acid synthase. Science 2008, 321(5894):1315-1322.
• [22]Gajendrarao P, Krishnamoorthy N, Sakkiah S, Lazar P, Lee KW: Molecular modeling study on orphan human protein CYP4A22 for identification of potential ligand binding site. J Mol Graph Model 2010, 28(6):524-532.
• [23]Bouaziz-Terrachet S, Toumi-Maouche A, Maouche B, Tairi-Kellou S: Modeling the binding modes of stilbene analogs to cyclooxygenase-2: a molecular docking study. J Mol Model 2010, 16(12):1919-1929.
• [24]Khare G, Gupta V, Gupta RK, Gupta R, Bhat R, Tyagi AK: Dissecting the role of critical residues and substrate preference of a fatty acyl-coa synthetase (FadD13) of Mycobacterium tuberculosis. PLoS One 2009, 4(12):e8387.
• [25]Zhao P, Liao QH, Ren CF, Wei J: Identification of ligand binding site on RXRgamma using molecular docking and dynamics methods. J Mol Model , 17(6):1259-1265.
• [26]Gabb HA, Jackson RM, Sternberg MJ: Modelling protein docking using shape complementarity, electrostatics and biochemical information. J Mol Biol 1997, 272(1):106-120.
• [27]Moont G, Gabb HA, Sternberg MJ: Use of pair potentials across protein interfaces in screening predicted docked complexes. Proteins 1999, 35(3):364-373.
• [28]Parris KD, Lin L, Tam A, Mathew R, Hixon J, Stahl M, Fritz CC, Seehra J, Somers WS: Crystal structures of substrate binding to Bacillus subtilis holo-(acyl carrier protein) synthase reveal a novel trimeric arrangement of molecules resulting in three active sites. Structure 2000, 8(8):883-895.
• [29]Anand S, Prasad MV, Yadav G, Kumar N, Shehara J, Ansari MZ, Mohanty D: SBSPKS: structure based sequence analysis of polyketide synthases. Nucleic Acids Res 2010, 38:W487-W496.
• [30]Morris GM, Goodsell DS, Halliday RS, Huey R, Hart W, Belew RK, Olson AJ: Automated docking using a lamarckian genetic algorithm and an empirical binding free energy function. JComputChem 1998, 19:1639-1662.
• [31]Alekseyev VY, Liu CW, Cane DE, Puglisi JD, Khosla C: Solution structure and proposed domain domain recognition interface of an acyl carrier protein domain from a modular polyketide synthase. Protein Sci 2007, 16(10):2093-2107.
• [32]Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol 1993, 234(3):779-815.
• [33]Zornetzer GA, Fox BG, Markley JL: Solution structures of spinach acyl carrier protein with decanoate and stearate. Biochemistry 2006, 45(16):5217-5227.
• [34]Gasteiger J, Marsili M: Iterative partial equalization of orbital electronegativity - a rapid access to atomic charges. Tetrahedron 1980, 36(22):3219-3228.
• [35]Case DA, Cheatham TE, Darden T, Gohlke H, Luo R, Merz KM, Onufriev A, Simmerling C, Wang B, Woods RJ: The amber biomolecular simulation programs. J Comput Chem 2005, 26(16):1668-1688.
• [36]Case DA, Darden TA, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo R, Merz KM, Pearlman DA, Crowley M, et al.: AMBER 9. University of California, San Francisco; 2006.
• [37]Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T, et al.: A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J Comput Chem 2003, 24(16):1999-2012.
• [38]Ryckaert JP, Ciccotti G, Berendsen HJC: Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys 1977, 23(3):327-341.
• [39]Darden TA, Pedersen LG: Molecular modeling: an experimental tool. Environ Health Perspect 1993, 101(5):410-412.