Modeling leucine's metabolic pathway and knockout prediction improving the production of surfactin, a biosurfactant from Bacillus subtilis.

A Bacillus subtilis mutant strain overexpressing surfactin biosynthetic genes was previously constructed. In order to further increase the production of this biosurfactant, our hypothesis is that the surfactin precursors, especially leucine, must be overproduced. We present a three step approach for leucine overproduction directed by methods from computational biology. Firstly, we develop a new algorithm for gene knockout prediction based on abstract interpretation, which applies to a recent modeling language for reaction networks with partial kinetic information. Secondly, we model the leucine metabolic pathway as a reaction network in this language, and apply the knockout prediction algorithm with the target of leucine overproduction. Out of the 21 reactions corresponding to potential gene knockouts, the prediction algorithm selects 12 reactions. Six knockouts were introduced in B. subtilis 168 derivatives strains to verify their effects on surfactin production. For all generated mutants, the specific surfactin production is increased from 1.6- to 20.9-fold during the exponential growth phase, depending on the medium composition. These results show the effectiveness of the knockout prediction approach based on formal models for metabolic reaction networks with partial kinetic information, and confirms our hypothesis that precursors supply is one of the main parameters to optimize surfactin overproduction.

[1]  S. Fisher,et al.  Regulation of nitrogen metabolism in Bacillus subtilis: vive la différence! , 1999, Molecular microbiology.

[2]  Jörg Stülke,et al.  SubtiWiki—a comprehensive community resource for the model organism Bacillus subtilis , 2011, Nucleic Acids Res..

[3]  M. Landy,et al.  Bacillomycin; an antibiotic from Bacillus subtilis active against pathogenic fungi. , 1948, Proceedings of the Society for Experimental Biology and Medicine. Society for Experimental Biology and Medicine.

[4]  H. Furuta,et al.  A Two-Piece Derivative of a Group I Intron RNA as a Platform for Designing Self-Assembling RNA Templates to Promote Peptide Ligation , 2012, Journal of nucleic acids.

[5]  B. Palsson,et al.  Genome-scale models of microbial cells: evaluating the consequences of constraints , 2004, Nature Reviews Microbiology.

[6]  Fangfang Xia,et al.  Building the repertoire of dispensable chromosome regions in Bacillus subtilis entails major refinement of cognate large-scale metabolic model , 2012, Nucleic acids research.

[7]  R. Losick,et al.  Additional Targets of the Bacillus subtilis Global Regulator CodY Identified by Chromatin Immunoprecipitation and Genome-Wide Transcript Analysis , 2003, Journal of bacteriology.

[8]  J. Nielsen,et al.  Industrial systems biology. , 2010, Biotechnology and bioengineering.

[9]  A. Barabasi,et al.  Global organization of metabolic fluxes in the bacterium Escherichia coli , 2004, Nature.

[10]  A. Sonenshein,et al.  Activation of the Bacillus subtilis global regulator CodY by direct interaction with branched‐chain amino acids , 2004, Molecular microbiology.

[11]  Jason A. Papin,et al.  Comparison of network-based pathway analysis methods. , 2004, Trends in biotechnology.

[12]  Jennifer L. Reed,et al.  OptORF: Optimal metabolic and regulatory perturbations for metabolic engineering of microbial strains , 2010, BMC Systems Biology.

[13]  Anthony J. Wilkinson,et al.  Genetic and Biochemical Analysis of the Interaction of Bacillus subtilis CodY with Branched-Chain Amino Acids , 2009, Journal of bacteriology.

[14]  Jörg Stülke,et al.  Essential genes in Bacillus subtilis: a re-evaluation after ten years. , 2013, Molecular bioSystems.

[15]  T. Kaneda,et al.  Iso- and anteiso-fatty acids in bacteria: biosynthesis, function, and taxonomic significance. , 1991, Microbiological reviews.

[16]  Costas D. Maranas,et al.  OptForce: An Optimization Procedure for Identifying All Genetic Manipulations Leading to Targeted Overproductions , 2010, PLoS Comput. Biol..

[17]  Sang Yup Lee,et al.  In silico genome-scale metabolic analysis of Pseudomonas putida KT2440 for polyhydroxyalkanoate synthesis, degradation of aromatics and anaerobic survival. , 2010, Biotechnology journal.

[18]  P. Jacques Surfactin and Other Lipopeptides from Bacillus spp. , 2011 .

[19]  Susanne Hennig,et al.  Transcriptional Organization and Posttranscriptional Regulation of the Bacillus subtilis Branched-Chain Amino Acid Biosynthesis Genes , 2004, Journal of bacteriology.

[20]  Joachim Niehren,et al.  Qualitative Reasoning for Reaction Networks with Partial Kinetic Information , 2015, CMSB.

[21]  E. Duitman Nonribosomal peptide synthesis in Bacillus subtilis , 2003 .

[22]  A. Sonenshein,et al.  Bacillus subtilis CodY represses early-stationary-phase genes by sensing GTP levels. , 2001, Genes & development.

[23]  M. Arnaud,et al.  Role of BkdR, a Transcriptional Activator of the SigL-Dependent Isoleucine and Valine Degradation Pathway inBacillus subtilis , 1999, Journal of bacteriology.

[24]  J. Errington,et al.  Essential Bacterial Functions Encoded by Gene Pairs , 2006, Journal of bacteriology.

[25]  Uwe Sauer,et al.  Regulation of CodY Activity through Modulation of Intracellular Branched-Chain Amino Acid Pools , 2010, Journal of bacteriology.

[26]  J. Guez,et al.  Effect of pps disruption and constitutive expression of srfA on surfactin productivity, spreading and antagonistic properties of Bacillus subtilis 168 derivatives , 2010, Journal of applied microbiology.

[27]  B. Schwikowski,et al.  Condition-Dependent Transcriptome Reveals High-Level Regulatory Architecture in Bacillus subtilis , 2012, Science.

[28]  Jörg Stülke,et al.  SPABBATS: A pathway-discovery method based on Boolean satisfiability that facilitates the characterization of suppressor mutants , 2011, BMC Systems Biology.

[29]  H. E. Umbarger,et al.  Threonine deaminase from Bacillus subtilis. II. The steady state kinetic properties. , 1970, Journal of Biological Chemistry.

[30]  A. Sonenshein,et al.  Contributions of Multiple Binding Sites and Effector-Independent Binding to CodY-Mediated Regulation in Bacillus subtilis , 2010, Journal of bacteriology.

[31]  Joachim Niehren,et al.  Knockout Prediction for Reaction Networks with Partial Kinetic Information , 2012, VMCAI.

[32]  B. Palsson,et al.  Constraints-based models: regulation of gene expression reduces the steady-state solution space. , 2003, Journal of theoretical biology.

[33]  P. Jauregi,et al.  Micelle size characterization of lipopeptides produced by B. subtilis and their recovery by the two-step ultrafiltration process. , 2013 .

[34]  Y. Fujita,et al.  CcpA-Mediated Catabolite Activation of the Bacillus subtilis ilv-leu Operon and Its Negation by Either CodY- or TnrA-Mediated Negative Regulation , 2014, Journal of bacteriology.

[35]  J. Calvo,et al.  Transcriptional regulation of the ilv-leu operon of Bacillus subtilis , 1992, Journal of bacteriology.

[36]  Bozhong Mu,et al.  Effects of Different Amino Acids in Culture Media on Surfactin Variants Produced by Bacillus subtilis TD7 , 2012, Applied Biochemistry and Biotechnology.

[37]  W. Wiechert 13C metabolic flux analysis. , 2001, Metabolic engineering.

[38]  A. Sonenshein,et al.  Dissecting Complex Metabolic Integration Provides Direct Genetic Evidence for CodY Activation by Guanine Nucleotides , 2011, Journal of bacteriology.

[39]  Christian Jungreuthmayer,et al.  Designing optimal cell factories: integer programming couples elementary mode analysis with regulation , 2012, BMC Systems Biology.

[40]  Zhaoxin Lu,et al.  Enhancement of surfactin production of Bacillus subtilis fmbR by replacement of the native promoter with the Pspac promoter. , 2009, Canadian journal of microbiology.

[41]  Steffen Klamt,et al.  In silico profiling of Escherichia coli and Saccharomyces cerevisiae as terpenoid factories , 2013, Microbial Cell Factories.

[42]  Yasutaro Fujita,et al.  Elaborate transcription regulation of the Bacillus subtilis ilv‐leu operon involved in the biosynthesis of branched‐chain amino acids through global regulators of CcpA, CodY and TnrA , 2005, Molecular microbiology.

[43]  J. Guez,et al.  Production of a novel mixture of mycosubtilins by mutants of Bacillus subtilis. , 2013, Bioresource technology.

[44]  Vincent Fromion,et al.  Reconstruction and analysis of the genetic and metabolic regulatory networks of the central metabolism of Bacillus subtilis , 2008, BMC Systems Biology.

[45]  A. Burgard,et al.  Optknock: A bilevel programming framework for identifying gene knockout strategies for microbial strain optimization , 2003, Biotechnology and bioengineering.