Genetic Circuits CI Dimer CI Dimer Dimerization Degradation CII Protein CI Protein Translation Repression mRNA Transcription Activation Pre Pr DNA RNAP RNAP RNAP RNAP RNAP O E O R cI cII Operator Sites Promoters Genes Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic Circuits CI Dimer CI Dimer Dimerization Degradation CII Protein CI Protein Translation Repression mRNA Transcription Activation Pre Pr DNA RNAP RNAP RNAP RNAP RNAP O E O R cI cII Operator Sites Promoters Genes Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic Circuits CI Dimer CI Dimer Dimerization Degradation CII Protein CI Protein Translation Repression mRNA Transcription Activation Pre Pr DNA RNAP RNAP RNAP RNAP RNAP O E O R cI cII Operator Sites Promoters Genes Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic Circuits CI Dimer CI Dimer Dimerization Degradation CII Protein CI Protein Translation Repression mRNA Transcription Activation Pre Pr DNA RNAP RNAP RNAP RNAP RNAP O E O R cI cII Operator Sites Promoters Genes Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic Circuits CI Dimer CI Dimer Dimerization Degradation CII Protein CI Protein Translation Repression mRNA Transcription Activation Pre Pr DNA RNAP RNAP RNAP RNAP RNAP O E O R cI cII Operator Sites Promoters Genes Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Francis Crick DNA makes RNA, RNA makes protein, and proteins make us. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Johann Von Neumann The sciences do not try to explain, they hardly even try to interpret, they mainly make models. By a model is meant a mathematical construct which, with the addition of certain verbal interpretations, describes observed phenomena. The justification of such a mathematical construct is solely and precisely that it is expected to work. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Scott Adams There are many methods for predicting the future. For example, you can read horoscopes, tea leaves, tarot cards, or crystal balls. Collectively, these methods are known as “nutty methods.” Or you can put well-researched facts into sophisticated computer models, more commonly referred to as “a complete waste of time.” Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic Circuit Model (GCM) Provides a higher level of abstraction than SBML. Includes only important species and their influences upon each other. GCMs also include structural constructs that allow us to connect GCMs for separate modules through species ports. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
A Genetic Not Gate A A C P1 C P1 c A C Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
A Genetic Nor Gate A B C A B P1 P1 C P1 c A C B Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
A Genetic Nand Gate A C A B P1 P2 P1 c B C C P2 c A C B Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
A Genetic Oscillator CI Dimer CI Dimer CI CII Protein CI Protein Pre Pr CII Pre Pr O E O R cI cII CI CII Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SBML: Main Elements Species Global parameters (ex. k1=0.1) Reactions Reactants Products Modifiers Stoichiometry Reversible Kinetic laws Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Synthesizing SBML from a GCM Representation Create degradation reactions Create open complex formation reactions Create dimerization reactions Create repression reactions Create activation reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
GCM Example CI Pre Pr CII Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Degradation Reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Open Complex Formation Reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Dimerization Reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Repression Reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Activation Reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Complete SBML Model Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Classical Chemical Kinetics Uses ordinary differential equations (ODE) to represent the system to be analyzed, and it assumes: Molecule counts are high, so concentrations can be continuous variables. Reactions occur continuously and deterministically. Genetic circuits have: Small molecule counts which must be considered as discrete variables. Gene expression reactions that occur sporadically. ODEs do not capture non-deterministic behavior. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Richard Feynmann A philosopher once said “It is necessary for the very existence of science that the same conditions always produce the same results.” Well, they do not. You set up the circumstances, with the same conditions every time, and you cannot predict behind which hole you will see the electron. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
NYTimes: Expressing Our Individuality, the Way E. Coli Do Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Rainbow and CC Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Stochastic Chemical Kinetics To more accurately predict the temporal behavior of genetic circuits, stochastic chemical kinetics formalism can be used. Use Gillespie’s Stochastic Simulation Algorithm which 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. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
� � � � � � � Automatic Abstraction Markov � Abstracted Reaction Reaction-based State-based � SAC Reaction Chain Model Abstraction Abstraction Model Analysis Model �� �� Stochastic � Results Simulation Begins with a reaction-based model in SBML. Automatically abstracts this model leveraging the quasi-steady state assumption, whenever possible. Encodes chemical species concentrations into Boolean (or n-ary) levels to produce a stochastic asynchronous circuit (SAC) model. Can now utilize Markov chain analysis. Kuwahara et al., Trans. on Comp. Sys. Bio. (2006) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Dimerization Reduction Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Dimerization Reduction Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Operator Site Reduction (PR) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Operator Site Reduction (PR) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Operator Site Reduction (PRE) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Operator Site Reduction (PRE) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Similar Reaction Combination Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Modifier Constant Propagation Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Final SBML Model 10 species and 10 reactions reduced to 2 species and 4 reactions Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
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. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
iBioSim: The Intelligent Biological Simulator Project management support. GCM Editor - creates Genetic Circuit Models (GCM). SBML Editor - creates models using the Systems Biology Markup Language (SBML). reb2sac - abstraction-based ODE, Monte Carlo, and Markov analysis. TSD Graph Editor - visualizes time series data (TSD). Probability Graph Editor - visualizes probability data. GeneNet - learns GCMs from TSD. Myers et al., Bioinformatics (2009) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
iBioSim: Genetic Circuit Editor Myers et al., Bioinformatics (2009) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
iBioSim: SBML Editor Myers et al., Bioinformatics (2009) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
iBioSim: Analysis Engine Myers et al., Bioinformatics (2009) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
ODE Results for the Simple Genetic Oscillator Comparison of ODE to SSA Results 6 5 6 0 5 5 5 0 Number of molecules 4 5 4 0 3 5 3 0 2 5 2 0 1 5 1 0 5 0 0 500 1,000 1,500 2,000 2,500 3,000 3,500 4,000 4,500 5,000 Time (s) CI_total (ODE) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
SSA Results for the Simple Genetic Oscillator Comparison of ODE to SSA Results 150 140 130 120 110 Number of molecules 100 9 0 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 0 0 0 500 1,000 1,500 2,000 2,500 3,000 3,500 4,000 4,500 5,000 Time (s) CI_total (ODE) CI_total (SSA) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic Muller C-Element A A B C 0 0 0 C 0 1 C C 1 0 C B 1 1 1 Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Toggle Switch C-Element (Genetic Circuit) A D X P1 d x B E X Y A C X Y S Q B P2 P3 e x y D R F A D Z F P7 f E F Z B E P8 P4 f z C Z Y P5 P6 c y z Nguyen et al., 13th Symposium on Async. Ckts. & Sys., 2007 ( best paper ) Nguyen et al., to appear in the Journal of Theoretical Biology Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Toggle Switch C-Element (GCM) A D X P1 d x B E X Y P2 P3 e x y D F P7 f E F Z P8 P4 f z C Z Y P5 P6 c y z Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Toggle Switch C-Element (SBML) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Toggle Switch C-Element (Abstracted) Reduced from 34 species and 31 reactions to 9 species and 15 reactions. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Toggle Switch C-Element (Simulation) Simulation time improved from 312 seconds to 20 seconds. Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Majority Gate C-Element (Genetic Circuit) A X B Y E C D Z X A Y D P1 P5 x y d Z X B E D C P2 P4 P7 P8 z x d e c Y Z D D P3 P6 y z d Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Speed-Independent C-Element (Genetic Circuit) A S1 S2 B S3 S4 C S4 A X S1 S2 S3 P1 P4 P5 P6 s4 x s1 s2 s3 B S4 Y S3 S1 Z S2 P2 P7 P8 s4 y s1 s2 z S2 S4 S3 Z C S4 P3 P9 P10 s3 z c s4 Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Genetic C-element Failures Toggle, Inputs Mixed dY=0 120 dZ=0 100 80 Stable Y 60 ? 40 20 Unstable Stable 0 0 20 40 60 80 100 120 Z Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Comparison of Failure Rates for the C-element Designs Failure Rate for Each C-Element Design 0.02 Majority Gate (High to Low) Speed-Independent (High to Low) Toggle Switch (High to Low) 0.015 Failure Rate 0.01 0.005 0 0 500 1000 1500 2000 Time (s) Chris Myers (U. of Utah) Synthetic Biology Carnegie Mellon University
Recommend
More recommend