Automatic Design of Synthetic Gene Circuits through Mixed Integer Non-linear Programming

Automatic design of synthetic gene circuits poses a significant challenge to synthetic biology, primarily due to the complexity of biological systems, and the lack of rigorous optimization methods that can cope with the combinatorial explosion as the number of biological parts increases. Current optimization methods for synthetic gene design rely on heuristic algorithms that are usually not deterministic, deliver sub-optimal solutions, and provide no guaranties on convergence or error bounds. Here, we introduce an optimization framework for the problem of part selection in synthetic gene circuits that is based on mixed integer non-linear programming (MINLP), which is a deterministic method that finds the globally optimal solution and guarantees convergence in finite time. Given a synthetic gene circuit, a library of characterized parts, and user-defined constraints, our method can find the optimal selection of parts that satisfy the constraints and best approximates the objective function given by the user. We evaluated the proposed method in the design of three synthetic circuits (a toggle switch, a transcriptional cascade, and a band detector), with both experimentally constructed and synthetic promoter libraries. Scalability and robustness analysis shows that the proposed framework scales well with the library size and the solution space. The work described here is a step towards a unifying, realistic framework for the automated design of biological circuits.

[1]  Madhukar S. Dasika,et al.  OptCircuit: An optimization based method for computational design of genetic circuits , 2008, BMC Systems Biology.

[2]  Bor-Sen Chen,et al.  Robust synthetic gene network design via library-based search method , 2011, Bioinform..

[3]  Andrew Phillips,et al.  Towards programming languages for genetic engineering of living cells , 2009, Journal of The Royal Society Interface.

[4]  Jay D Keasling,et al.  Importance of systems biology in engineering microbes for biofuel production. , 2008, Current opinion in biotechnology.

[5]  Ron Weiss,et al.  Engineering signal processing in cells: Towards molecular concentration band detection , 2002, Natural Computing.

[6]  M. Elowitz,et al.  Modeling a synthetic multicellular clock: repressilators coupled by quorum sensing. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Jörg Stelling,et al.  Automatic Design of Digital Synthetic Gene Circuits , 2011, PLoS Comput. Biol..

[8]  Tom E. Bishop,et al.  Blind Image Restoration Using a Block-Stationary Signal Model , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[9]  Ron Weiss,et al.  Parameter estimation for two synthetic gene networks: a case study , 2005, Proceedings. (ICASSP '05). IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005..

[10]  Terence Hwa,et al.  Combinatorial transcriptional control of the lactose operon of Escherichia coli , 2007, Proceedings of the National Academy of Sciences.

[11]  J. Griffith Mathematics of cellular control processes. II. Positive feedback to one gene. , 1968, Journal of theoretical biology.

[12]  V. Hakim,et al.  Design of genetic networks with specified functions by evolution in silico. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[13]  J. Griffith,et al.  Mathematics of cellular control processes. I. Negative feedback to one gene. , 1968, Journal of theoretical biology.

[14]  A. L. Koch,et al.  Protein degradation in Escherichia coli. I. Measurement of rapidly and slowly decaying components. , 1970, The Journal of biological chemistry.

[15]  R. Weiss,et al.  Ultrasensitivity and noise propagation in a synthetic transcriptional cascade. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[16]  J. Cronan,et al.  Long-term and homogeneous regulation of the Escherichia coli araBAD promoter by use of a lactose transporter of relaxed specificity , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[17]  J. Collins,et al.  DIVERSITY-BASED, MODEL-GUIDED CONSTRUCTION OF SYNTHETIC GENE NETWORKS WITH PREDICTED FUNCTIONS , 2009, Nature Biotechnology.

[18]  R. Weiss,et al.  Optimizing genetic circuits by global sensitivity analysis. , 2004, Biophysical journal.

[19]  S. Basu,et al.  A synthetic multicellular system for programmed pattern formation , 2005, Nature.

[20]  M. Elowitz,et al.  A synthetic oscillatory network of transcriptional regulators , 2000, Nature.

[21]  R. Weiss,et al.  Automatic Compilation from High-Level Biologically-Oriented Programming Language to Genetic Regulatory Networks , 2011, PloS one.

[22]  R. Bellman Stability theory of differential equations , 1953 .

[23]  Timothy K Lu,et al.  Engineered bacteriophage targeting gene networks as adjuvants for antibiotic therapy , 2009, Proceedings of the National Academy of Sciences.

[24]  James N. Weiss The Hill equation revisited: uses and misuses , 1997, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[25]  Leo Liberti,et al.  Branching and bounds tighteningtechniques for non-convex MINLP , 2009, Optim. Methods Softw..

[26]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[27]  Douglas Densmore,et al.  Rule based constraints for the construction of genetic devices , 2010, Proceedings of 2010 IEEE International Symposium on Circuits and Systems.

[28]  J. Collins,et al.  Construction of a genetic toggle switch in Escherichia coli , 2000, Nature.

[29]  B. Goodwin Oscillatory behavior in enzymatic control processes. , 1965, Advances in enzyme regulation.

[30]  Alfonso Jaramillo,et al.  Computational design of synthetic regulatory networks from a genetic library to characterize the designability of dynamical behaviors , 2011, Nucleic acids research.

[31]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[32]  W. J. Meek,et al.  THE FEDERATION OF AMERICAN SOCIETIES FOR EXPERIMENTAL BIOLOGY. , 1938, Science.

[33]  Alfonso Jaramillo,et al.  Genetdes: automatic design of transcriptional networks , 2007, Bioinform..

[34]  Joseph H. Davis,et al.  Design, construction and characterization of a set of insulated bacterial promoters , 2010, Nucleic acids research.

[35]  Yiannis N Kaznessis,et al.  Model-driven designs of an oscillating gene network. , 2005, Biophysical journal.

[36]  Efstratios N. Pistikopoulos,et al.  New algorithms for mixed-integer dynamic optimization , 2003, Comput. Chem. Eng..