Design Principles in Synthetic Biology Chris Myers 1 , Nathan Barker - - PowerPoint PPT Presentation

design principles in synthetic biology
SMART_READER_LITE
LIVE PREVIEW

Design Principles in Synthetic Biology Chris Myers 1 , Nathan Barker - - PowerPoint PPT Presentation

Design Principles in Synthetic Biology Chris Myers 1 , Nathan Barker 2 , Hiroyuki Kuwahara 3 , Curtis Madsen 1 , Nam Nguyen 1 , Michael Samoilov 4 , and Adam Arkin 4 1 University of Utah 2 Southern Utah University 3 Microsoft Research, Trento, Italy


slide-1
SLIDE 1

Design Principles in Synthetic Biology

Chris Myers1, Nathan Barker2, Hiroyuki Kuwahara3, Curtis Madsen1, Nam Nguyen1, Michael Samoilov4, and Adam Arkin4

1University of Utah 2Southern Utah University 3Microsoft Research, Trento, Italy 4University of California, Berkeley

Design Principles in Biological Systems April 24, 2008

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-2
SLIDE 2

Synthetic Biology

Increasing number of labs are designing more ambitious and mission critical synthetic biology projects. These projects construct synthetic genetic circuits from DNA. These synthetic genetic circuits can potentially result in:

A better understanding of how microorganisms function by examining differences in vivo compared to in silico (Sprinzak/Elowitz). More efficient pathways for the production of antimalarial drugs (Dae et al.). Bacteria that can metabolize toxic chemicals (Brazil et al.). Bacteria that can hunt and kill tumors (Anderson et al.).

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-3
SLIDE 3

Genetic Design Automation (GDA)

Electronic Design Automation (EDA) tools have facilitated the design of ever more complex integrated circuits each year. Crucial to the success of synthetic biology is an improvement in methods and tools for Genetic Design Automation (GDA). Existing GDA tools require biologists to design at the molecular level. Roughly equivalent to designing electronic circuits at the layout level. Analysis of genetic circuits is also performed at this very low level. A GDA tool that supports higher levels of abstraction is essential.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-4
SLIDE 4

Overview

This talk describes our research to develop such a GDA tool. This tool has helped us examine design principles for synthetic biology. As a case study, will describe the design of a genetic Muller C-element.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-5
SLIDE 5

Current State of GDA Tools

MIT has created a registry of standard biological parts used to design synthetic genetic circuits (http://parts.mit.edu). Methods and tools are needed to assist in the design and analysis of synthetic genetic circuits using these parts. BioJADE provides a schematic capture interface to the MIT parts registry. Systems Biology Markup Language (SBML) has been proposed as a standard representation for the simulation of biological systems. Many simulation tools have been developed that accept models in the SBML format (BioPathwise, BioSPICE, CellDesigner, SimBiology, etc.).

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-6
SLIDE 6

Systems Biology Markup Language (SBML)

SBML models biological systems at the molecular level. A typical SBML model is composed of a number of chemical species (i.e., proteins, genes, etc.) and reactions that transform these species. This is a very low level representation which is roughly equivalent to the layout level for electronic circuits. Designing and simulating genetic circuits at this level of detail is extremely tedious and time-consuming. Therefore, there is a need for higher-level abstractions for modeling, analysis, and design of genetic circuits.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-7
SLIDE 7

BioSim

Genetic Circuit

  • Insert into

Host

  • Plasmid
  • Perform

Experiments

  • Construct

Plasmid

  • Experimental

Data

  • Biological

Knowledge

  • DNA

Sequence

  • Learn Model

GCM Synthesis

  • Simulation

Data Abstraction/ Simulation

  • SBML Model
  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-8
SLIDE 8

BioSim: Analysis

Genetic Circuit

  • Insert into

Host

  • Plasmid
  • Perform

Experiments

  • Construct

Plasmid

  • Experimental

Data

  • Biological

Knowledge

  • DNA

Sequence

  • Learn Model

GCM Synthesis

  • Simulation

Data Abstraction/ Simulation

  • SBML Model
  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-9
SLIDE 9

BioSim: Design

Genetic Circuit

  • Insert into

Host

  • Plasmid
  • Perform

Experiments

  • Construct

Plasmid

  • Experimental

Data

  • Biological

Knowledge

  • DNA

Sequence

  • Learn Model

GCM Synthesis

  • Simulation

Data Abstraction/ Simulation

  • SBML Model
  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-10
SLIDE 10

BioSim: Design

Genetic Circuit

  • Insert into

Host

  • Plasmid
  • Perform

Experiments

  • Construct

Plasmid

  • Experimental

Data

  • Biological

Knowledge

  • DNA

Sequence

  • Learn Model

GCM Synthesis

  • Simulation

Data Abstraction/ Simulation

  • SBML Model
  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-11
SLIDE 11

BioSim: Genetic Circuit Model

Genetic Circuit

  • Insert into

Host

  • Plasmid
  • Perform

Experiments

  • Construct

Plasmid

  • Experimental

Data

  • Biological

Knowledge

  • DNA

Sequence

  • Learn Model

GCM Synthesis

  • Simulation

Data Abstraction/ Simulation

  • SBML Model
  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-12
SLIDE 12

Phage λ Virus

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-13
SLIDE 13

Phage λ Decision Circuit

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-14
SLIDE 14

Phage λ Decision Circuit

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-15
SLIDE 15

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites Promoters Genes Transcription cI cII CI Dimer DNA

Pre OR OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-16
SLIDE 16

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites Promoters Transcription CI Dimer DNA

Pre

Genes

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-17
SLIDE 17

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites Transcription CI Dimer DNA

Pre

Promoters Genes

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-18
SLIDE 18

Genetic Circuits

RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

RNAP CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites CI Dimer DNA

Pre

Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-19
SLIDE 19

Genetic Circuits

RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

RNAP CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites CI Dimer DNA

Pre

Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-20
SLIDE 20

Genetic Circuits

RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

RNAP CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites CI Dimer DNA

Pre

Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-21
SLIDE 21

Genetic Circuits

RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

RNAP

Pr

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites CI Dimer DNA

Pre

Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-22
SLIDE 22

Genetic Circuits

RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

RNAP

Pr

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites CI Dimer DNA

Pre

Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-23
SLIDE 23

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein Translation CII Protein Operator Sites Transcription CI Dimer DNA

Pre

mRNA Promoters Genes

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-24
SLIDE 24

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein Operator Sites CI Dimer DNA

Pre

mRNA Translation CII Protein Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-25
SLIDE 25

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein CI Dimer DNA

Pre

mRNA Translation CII Protein Operator Sites Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-26
SLIDE 26

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer CI Protein CI Dimer DNA

Pre

Activation mRNA Translation CII Protein Operator Sites Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-27
SLIDE 27

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer CI Dimer DNA

Pre

Activation CI Protein mRNA Translation CII Protein Operator Sites Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-28
SLIDE 28

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer DNA

Pre

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-29
SLIDE 29

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Degradation Repression Dimerization

Pr

CI Dimer DNA

Pre

CI Dimer Activation CI Protein mRNA Translation CII Protein Operator Sites Promoters Genes Transcription

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-30
SLIDE 30

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

mRNA Translation Transcription CI Dimer DNA

Pre

CI Dimer Activation CI Protein CII Protein Operator Sites Promoters Genes

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-31
SLIDE 31

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Degradation Dimerization Repression

Pr

Activation CI Protein mRNA Translation CII Protein Transcription CI Dimer DNA

Pre

CI Dimer Operator Sites Promoters Genes

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-32
SLIDE 32

Genetic Circuits

RNAP RNAP RNAP RNAP RNAP

Repression Degradation Dimerization

Pr

CI Dimer Activation CI Protein Transcription CI Dimer DNA

Pre

mRNA Translation CII Protein Operator Sites Promoters Genes

OR

cI cII

OE

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-33
SLIDE 33

Logical Representation

CI CII

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-34
SLIDE 34

Graphical Representation

CI CII Pre Pr

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-35
SLIDE 35

Genetic Circuit Model (GCM)

Provides a higher level of abstraction than SBML. Includes only important species and their influences upon each other. A GCM is a tuple S,P,G,I,Sd where:

S is a finite set of species; P is a finite set of promoters; G : P → 2S maps promoters to sets of species; I ⊆ S × P ×{a,r} is a finite set of influences; Sd ⊆ S is a set of species that influence as dimers.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-36
SLIDE 36

GCM Graphical Representation

A bipartite graph with species and promoters as the two types of nodes. Species are connected to promoters using influences I, and promoters are connected to species using function G. To simplify presentation, graphs shown using only species as nodes, edges are inferred using I and G, and edges are labeled with the promoter that links the species.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-37
SLIDE 37

Influences on the Same Promoter

B C

P1 P1

A

C A B P1

c

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-38
SLIDE 38

Influences on the Same Promoter

B C

P1 P1

A

C A B P1

c

B A C

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-39
SLIDE 39

Influences on Different Promoters

A B C

P1 P2

C B C A P1 P2

c c

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-40
SLIDE 40

Influences on Different Promoters

A B C

P1 P2

C B C A P1 P2

c c

C A B

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-41
SLIDE 41

GCM Parameters

Parameter Sym Structure Value Units Initial species count ns species molecule Dimerization equilibrium Kd species .05

1 molecule

Degradation rate kd species .0075

1 sec

Initial promoter count ng promoter 2 molecule Stoichiometry of production np promoter 10 molecule Degree of cooperativity nc promoter 2 molecule RNAP binding equilibrium Ko promoter .033

1 molecule

Open complex production rate ko promoter .05

1 sec

Basal production rate kb promoter .0001

1 sec

Activated production rate ka promoter .25

1 sec

Repression binding equilibrium Kr influence .5

1 moleculenc

Activation binding equilibrium Ka influence .0033

1 molecule(nc+1)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-42
SLIDE 42

GCM versus SBML Representation

CI CII Pre Pr

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-43
SLIDE 43

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-44
SLIDE 44

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-45
SLIDE 45

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-46
SLIDE 46

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-47
SLIDE 47

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-48
SLIDE 48

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-49
SLIDE 49

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-50
SLIDE 50

SBML Example

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-51
SLIDE 51

Synthesizing SBML from a GCM Representation

Create degradation reactions Create open complex formation reactions Create dimerization reactions Create repression reactions Create activation reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-52
SLIDE 52

Degradation Reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-53
SLIDE 53

Open Complex Formation Reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-54
SLIDE 54

Dimerization Reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-55
SLIDE 55

Repression Reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-56
SLIDE 56

Activation Reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-57
SLIDE 57

Complete SBML Model

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-58
SLIDE 58

Classical Chemical Kinetics

Uses ordinary differential equations (ODE) to represent the system to be analyzed, and it assumes:

A system is well-stirred. Number of molecules in a cell is high. Concentrations can be viewed as continuous variables. Reactions occur continuously and deterministically.

Genetic circuits involve small molecule counts. Gene expression can have substantial fluctuations. ODEs do not capture non-deterministic behavior.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-59
SLIDE 59

Stochastic Chemical Kinetics

To more accurately predict the temporal behavior of genetic circuits, stochastic chemical kinetics formalism can be used. Probabilistically predicts the dynamics of biochemical systems. Describes the time evolution of a system as a discrete-state jump Markov process governed by the chemical master equation (CME). Can simulate it using Gillespie’s Stochastic Simulation Algorithm (SSA). It exactly tracks the quantities of each molecular species, and treats each reaction as a separate random event. Only practical for small systems with no major time-scale separations. Abstraction is essential for efficient analysis of any realistic system.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-60
SLIDE 60

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-61
SLIDE 61

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-62
SLIDE 62

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-63
SLIDE 63

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-64
SLIDE 64

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-65
SLIDE 65

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-66
SLIDE 66

Automatic Abstraction

Reaction Model

  • Reaction-based

Abstraction Abstracted Reaction Model

  • State-based

Abstraction SAC Model

  • Markov

Chain Analysis

  • Stochastic

Simulation Results

Begins with a reaction-based model in SBML. Next, it automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Finally, it encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit model. It can now be analyzed using Markov chain analysis.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-67
SLIDE 67

Dimerization Reduction

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-68
SLIDE 68

Dimerization Reduction

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-69
SLIDE 69

Operator Site Reduction (PR)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-70
SLIDE 70

Operator Site Reduction (PR)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-71
SLIDE 71

Operator Site Reduction (PRE)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-72
SLIDE 72

Operator Site Reduction (PRE)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-73
SLIDE 73

Similar Reaction Combination

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-74
SLIDE 74

Modifier Constant Propagation

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-75
SLIDE 75

Final SBML Model

10 species and 10 reactions reduced to 2 species and 4 reactions

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-76
SLIDE 76

BioSim: Genetic Circuit Editor

http://www.async.ece.utah.edu/BioSim/

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-77
SLIDE 77

BioSim: SBML Editor

http://www.async.ece.utah.edu/BioSim/

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-78
SLIDE 78

BioSim: Simulator

http://www.async.ece.utah.edu/BioSim/

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-79
SLIDE 79

BioSim: Parameter Editor

http://www.async.ece.utah.edu/BioSim/

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-80
SLIDE 80

BioSim: Graph Editor

http://www.async.ece.utah.edu/BioSim/

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-81
SLIDE 81

GCM Advantages

Greatly increases the speed of model development and reduces the number of errors in the resulting models. Allows efficient exploration of the effects of parameter variation. Constrains SBML model such that it can be more easily abstracted resulting in substantial improvement in simulation time.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-82
SLIDE 82

Genetic Muller C-Element C

B A C’ A B C’ 1 C 1 C 1 1 1

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-83
SLIDE 83

Toggle Switch C-Element (Genetic Circuit)

B A E D F B A X Y Z C

Q S R P1 P2 P3 P7 P8 P4 P5 P6 X X Y A B E D F Z F D C Y Z E x d e x y f f z c y z

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-84
SLIDE 84

Toggle Switch C-Element (GCM)

P1 P2 P3 P7 P8 P4 P5 P6 X X Y A B E D F Z F D C Y Z E x d e x y f f z c y z

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-85
SLIDE 85

Toggle Switch C-Element (SBML)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-86
SLIDE 86

Toggle Switch C-Element (Abstracted)

Reduced from 34 species and 31 reactions to 9 species and 15 reactions.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-87
SLIDE 87

Toggle Switch C-Element (Simulation)

Simulation time improved from 312 seconds to 20 seconds.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-88
SLIDE 88

Majority Gate C-Element (Genetic Circuit)

E C X Y Z D B A

P8 P7 P5 P6 P4 P3 P2 P1 A B X D D Y C E D D Y Z Z X x y d d e c d y z z x

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-89
SLIDE 89

Majority Gate C-Element (GCM)

P8 P7 P5 P6 P4 P3 P2 P1 A B X D D Y C E D D Y Z Z X x y d d e c d y z z x

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-90
SLIDE 90

Majority Gate C-Element (Simulation)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-91
SLIDE 91

Speed-Independent C-Element (Genetic Circuit)

A B S1 S2 S3 S4 C

Z P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 S4 S4 X S1 S3 S2 S4 A C S3 S2 S4 B S2 Z S3 Y S1 x s4 y s4 z s3 s1 s2 s3 s1 s2 z c s4

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-92
SLIDE 92

Speed-Independent C-Element (GCM)

Z P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 S4 S4 X S1 S3 S2 S4 A C S3 S2 S4 B S2 Z S3 Y S1 x s4 y s4 z s3 s1 s2 s3 s1 s2 z c s4

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-93
SLIDE 93

Speed-Independent C-Element (Simulation)

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-94
SLIDE 94

Ordinary Differential Equation Analysis

Use Law of Mass Action to derive an ODE model. Study behavior of our model at steady state. Analyze nullclines to characterize the gate.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-95
SLIDE 95

ODE Analysis: Nullclines for Toggle C-Element

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs low Z Y dY=0 dZ=0

Stable

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-96
SLIDE 96

ODE Analysis: Nullclines for Toggle C-Element

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0 20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0

Stable Unstable Stable

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-97
SLIDE 97

ODE Analysis: Nullclines for Toggle C-Element

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs High Z Y dY=0 dZ=0

Stable

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-98
SLIDE 98

ODE Analysis: Nullclines for Toggle C-Element

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0 20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0

Stable Unstable Stable

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-99
SLIDE 99

Stochastic Simulation: State Change from Low to High

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0

? Stable Unstable Stable

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-100
SLIDE 100

Stochastic Simulation: State Change from Low to High

500 1000 1500 2000 0.005 0.01 0.015 0.02 0.025 0.03 Low to High Time (s) Failure Rate maj−heat−high maj−light−high tog−heat−high tog−light−high si−heat−high si−light−high

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-101
SLIDE 101

Stochastic Simulation: State Change from High to Low

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0

Stable Unstable Stable ?

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-102
SLIDE 102

Stochastic Simulation: State Change from High to Low

500 1000 1500 2000 0.005 0.01 0.015 0.02 0.025 0.03 High to Low Time (s) Failure Rate maj−heat−low maj−light−low tog−heat−low tog−light−low si−heat−low si−light−low

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-103
SLIDE 103

Effect of Gene Count

1 1.5 2 2.5 3 3.5 4 4.5 5 0.05 0.1 0.15 0.2 0.25 Low to High Number of Genes Failure Rate maj−heat−high maj−light−high tog−heat−high tog−light−high si−heat−high si−light−high

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-104
SLIDE 104

Effect of Cooperativity

1 1.5 2 2.5 3 3.5 4 4.5 5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Low to High Cooperativity Failure Rate maj−heat−high maj−light−high tog−heat−high tog−light−high si−heat−high si−light−high

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-105
SLIDE 105

Effect of Repression Strength

10

−1

10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Low to High Repression Failure Rate maj−heat−high maj−light−high tog−heat−high tog−light−high si−heat−high si−light−high

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-106
SLIDE 106

Effect of Decay Rates

0.005 0.01 0.015 0.02 0.025 0.03 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Low to High Decay Rate Failure Rate maj−heat−high maj−light−high tog−heat−high tog−light−high si−heat−high si−light−high

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-107
SLIDE 107

Effect of Dual Rail

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0

? Stable Unstable Stable

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-108
SLIDE 108

Effect of Dual Rail

500 1000 1500 2000 0.005 0.01 0.015 0.02 0.025 0.03 Low to High Time (s) Failure Rate single−tog−heat−high single−tog−light−high dual−tog−heat−high dual−tog−light−high

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-109
SLIDE 109

Effect of Dual Rail

20 40 60 80 100 120 20 40 60 80 100 120 Toggle, Inputs Mixed Z Y dY=0 dZ=0

Stable Unstable Stable ?

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-110
SLIDE 110

Effect of Dual Rail

500 1000 1500 2000 0.05 0.1 0.15 0.2 0.25 High to Low Time (s) Failure Rate single−tog−heat−low single−tog−light−low dual−tog−heat−low dual−tog−light−low

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-111
SLIDE 111

Design Principles in Synthetic Biology

Speed-independence does not necessarily imply better robustness. Higher gene counts improve production rates, higher equilibrium values, and more robust operation. Cooperativity of at least two is required to produce the necessary non-linearity for state-holding. Repressors should bind efficiently. Decay rates cannot be too high. Dual-rail outputs are essential.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-112
SLIDE 112

Future Work: Modular Design

More levels of hierarchy are needed in the GCM format. We plan to create structural constructs that allow us to connect GCM’s for separate modules through species ports. Allow design at the logical and higher levels of abstraction.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-113
SLIDE 113

Biologically Inspired Circuit Design

Human inner ear performs the equivalent of one billion floating point

  • perations per second and consumes only 14 µW while a game console

with similar performance burns about 50 W (Sarpeshkar, 2006). We believe this difference is due to over designing components in order to achieve an extremely low probability of failure in every device. Future silicon and nano-devices will be much less reliable. For Moore’s law to continue, future design methods should support the design of reliable systems using unreliable components. Biological systems constructed from very noisy and unreliable devices. GDA tools may be useful for future integrated circuit technologies. Biological systems tend to be more asynchronous and analog in nature, so future engineered circuits will likely need to be also.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008

slide-114
SLIDE 114

Acknowledgments

Nathan Barker Hiroyuki Kuwahara Nam Nguyen Curtis Madsen Michael Samoilov Adam Arkin

This work is supported by the National Science Foundation under Grants No. 0331270 and CCF07377655.

  • C. Myers et al. (U. of Utah)

Design Principles in Synthetic Biology IMA Workshop / April 24, 2008