SPABBATS: A pathway-discovery method based on Boolean satisfiability that facilitates the characterization of suppressor mutants

BackgroundSeveral computational methods exist to suggest rational genetic interventions that improve the productivity of industrial strains. Nonetheless, these methods are less effective to predict possible genetic responses of the strain after the intervention. This problem requires a better understanding of potential alternative metabolic and regulatory pathways able to counteract the targeted intervention.ResultsHere we present SPABBATS, an algorithm based on Boolean satisfiability (SAT) that computes alternative metabolic pathways between input and output species in a reconstructed network. The pathways can be constructed iteratively in order of increasing complexity. SPABBATS allows the accumulation of intermediates in the pathways, which permits discovering pathways missed by most traditional pathway analysis methods. In addition, we provide a proof of concept experiment for the validity of the algorithm. We deleted the genes for the glutamate dehydrogenases of the Gram-positive bacterium Bacillus subtilis and isolated suppressor mutant strains able to grow on glutamate as single carbon source. Our SAT approach proposed candidate alternative pathways which were decisive to pinpoint the exact mutation of the suppressor strain.ConclusionsSPABBATS is the first application of SAT techniques to metabolic problems. It is particularly useful for the characterization of metabolic suppressor mutants and can be used in a synthetic biology setting to design new pathways with specific input-output requirements.

[1]  J A Chudek,et al.  The effects of osmotic upshock on the intracellular solute pools of Bacillus subtilis. , 1990, Journal of general microbiology.

[2]  A. Moir,et al.  The regulation of the fumarase (citG) gene of Bacillus subtilis 168 , 1988, Molecular and General Genetics MGG.

[3]  A. Sonenshein,et al.  Role and Regulation of Bacillus subtilisGlutamate Dehydrogenase Genes , 1998, Journal of bacteriology.

[4]  M. Hecker,et al.  Transcription of glycolytic genes and operons in Bacillus subtilis: evidence for the presence of multiple levels of control of the gapA operon , 2001, Molecular microbiology.

[5]  M. Page,et al.  Search for Steady States of Piecewise-Linear Differential Equation Models of Genetic Regulatory Networks , 2008, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[6]  C. Daub,et al.  BMC Systems Biology , 2007 .

[7]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[8]  B. Palsson,et al.  Genome-scale Reconstruction of Metabolic Network in Bacillus subtilis Based on High-throughput Phenotyping and Gene Essentiality Data* , 2007, Journal of Biological Chemistry.

[9]  Fabian M. Commichau,et al.  Characterization of Bacillus subtilis Mutants with Carbon Source-Independent Glutamate Biosynthesis , 2006, Journal of Molecular Microbiology and Biotechnology.

[10]  A. Sonenshein,et al.  Transcriptional regulation of Bacillus subtilis citrate synthase genes , 1994, Journal of bacteriology.

[11]  P. Stragier,et al.  Antibiotic-resistance cassettes for Bacillus subtilis. , 1995, Gene.

[12]  J. Stelling,et al.  Combinatorial Complexity of Pathway Analysis in Metabolic Networks , 2004, Molecular Biology Reports.

[13]  Jörg Stülke,et al.  Transcriptional profiling of gene expression in response to glucose in Bacillus subtilis: regulation of the central metabolic pathways. , 2003, Metabolic engineering.

[14]  Jörg Stülke,et al.  A regulatory protein–protein interaction governs glutamate biosynthesis in Bacillus subtilis: the glutamate dehydrogenase RocG moonlights in controlling the transcription factor GltC , 2007, Molecular microbiology.

[15]  Ashish Tiwari,et al.  Analyzing Pathways Using SAT-Based Approaches , 2007, AB.

[16]  Francisco J. Planes,et al.  A critical examination of stoichiometric and path-finding approaches to metabolic pathways , 2008, Briefings Bioinform..

[17]  P Setlow,et al.  Cloning and nucleotide sequence of the Bacillus subtilis ansR gene, which encodes a repressor of the ans operon coding for L-asparaginase and L-aspartase , 1993, Journal of bacteriology.

[18]  G. Rapoport,et al.  Salt stress is an environmental signal affecting degradative enzyme synthesis in Bacillus subtilis , 1995, Journal of bacteriology.

[19]  George Chin,et al.  BioGraphE: high-performance bionetwork analysis using the Biological Graph Environment , 2008, BMC Bioinformatics.

[20]  A. Wach PCR‐synthesis of marker cassettes with long flanking homology regions for gene disruptions in S. cerevisiae , 1996, Yeast.

[21]  M. Débarbouillé,et al.  Interactions of wild-type and truncated LevR of Bacillus subtilis with the upstream activating sequence of the levanase operon. , 1994, Journal of molecular biology.

[22]  Jörg Stülke,et al.  Regulatory links between carbon and nitrogen metabolism. , 2006, Current opinion in microbiology.

[23]  Angel Rubio,et al.  Computing the shortest elementary flux modes in genome-scale metabolic networks , 2009, Bioinform..

[24]  Jörg Stülke,et al.  Connecting parts with processes: SubtiWiki and SubtiPathways integrate gene and pathway annotation for Bacillus subtilis. , 2010, Microbiology.

[25]  Haisu Ma,et al.  Synthesizing a novel genetic sequential logic circuit: a push-on push-off switch , 2010, Molecular systems biology.

[26]  Thorsten Mascher,et al.  Bacitracin sensing in Bacillus subtilis , 2008, Molecular microbiology.

[27]  L. Wray,et al.  Bacillus subtilis 168 Contains Two Differentially Regulated Genes Encoding l-Asparaginase , 2002, Journal of bacteriology.

[28]  Fabian M. Commichau,et al.  Glutamate Metabolism in Bacillus subtilis: Gene Expression and Enzyme Activities Evolved To Avoid Futile Cycles and To Allow Rapid Responses to Perturbations of the System , 2008, Journal of bacteriology.

[29]  A. Sonenshein,et al.  Modulation of Activity of Bacillus subtilis Regulatory Proteins GltC and TnrA by Glutamate Dehydrogenase , 2004, Journal of bacteriology.

[30]  Martin Fränzle,et al.  Efficient Solving of Large Non-linear Arithmetic Constraint Systems with Complex Boolean Structure , 2007, J. Satisf. Boolean Model. Comput..

[31]  Ron Piran,et al.  Algorithm of myogenic differentiation in higher-order organisms , 2009, Development.

[32]  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.

[33]  Jörg Stülke,et al.  The regulatory link between carbon and nitrogen metabolism in Bacillus subtilis: regulation of the gltAB operon by the catabolite control protein CcpA. , 2003, Microbiology.

[34]  Xiao-Jiang Feng,et al.  Metabolomics-driven quantitative analysis of ammonia assimilation in E. coli , 2009, Molecular systems biology.

[35]  Inês Lynce,et al.  Efficient Haplotype Inference with Pseudo-boolean Optimization , 2007, AB.

[36]  D. Sun,et al.  Cloning, nucleotide sequence, and expression of the Bacillus subtilis ans operon, which codes for L-asparaginase and L-aspartase , 1991, Journal of bacteriology.

[37]  J. Sambrook,et al.  Molecular Cloning: A Laboratory Manual , 2001 .

[38]  Adam M. Feist,et al.  The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli , 2008, Nature Biotechnology.

[39]  Mary Sheeran,et al.  SAT-Solving in Practice, with a Tutorial Example from Supervisory Control , 2009, Discret. Event Dyn. Syst..