Revision as of 00:34, 18 October 2018

Project

Background Design Experiments Results

Modeling

Overview Modeling 101 Light Operon Turing Patterns CRISPR-Cas13a

Human Practices

Human Practices Public Engagement

Team

Team Members Collaborations Attributions

Top

Light Operon

1. Introduction

A deterministic model of differential equations over a continuous time interval is generally easy to implement using known numerical methods. However, the equations alone do not account for the variability in expression that we often observe at low molecular count. In other words, deterministic simulations fail to capture the actual physical basis of the reaction (Gillespie 1977; Wilkinson 2012). This is due to the intrinsic and extrinsic noise that occurs at low intracellular molecular counts; interactions do not follow traditionally defined constant parameters which limits the accurate characterization of a genetic system (Elowitz et al., 2002).

It then becomes essential when characterizing a genetic system to have a deterministic approach (to capture the large average overtime) in addition to a stochastic approach. A stochastic approach allows one to model the probabilistic trajectory based on the initial conditions. One commonly used stochastic algorithm is the Gillespie Algorithm— a simple but powerful approach to simulation that takes into account the initial state of the reactants, the reaction rate and the number of molecules present for a reaction with each timed step of the reaction drawn from a probability distribution (Gillespie 1977).

As a working example, we will be working with the pDawn and pDusk operon (BBa K161609 and BBa K1075044). Many teams have made either the pDawn/pDusk system as their inducer (IONIS Paris 2015, Wageningen 2016, NUS Singapore 2017, Cornell 2017, Kingsborough 2017). In 2017, NUS Singapore 2017 and Kingsborough 2017 used the pDawn/pDusk inducer to induce expression of MazF to induce cell death. To our knowledge, there is no literature or past iGEM project exploring a stochastic model of the pDawn/pDusk system at low molecule count.

2. Methods

In order to improve reproducibility and to build on past iGEM teams’ efforts, we used the differential equations and parameters described in the 2016 Wageningen iGEM Team’s kill switch design. Using the “smfsb” package in R for stochastic simulations and the NDSolve function in Mathematica, we modeled the 3 light-sensitive states of the Yf1 homodimer over a period of 20 hours. Yll frequency simulation was done through the "smfsb" package in R and simulated a hundred times. Histograms were made through Excel.

Equations, Constants, & Parameters

Post-Reaction Pre-Reaction Matrix for Gillespie

3. Results

The different states of Yf1 are being hit with different intensities of light, producing alternating results unique to the state. Click on any of the images to view full size.

In the stochastic results, we observe a much greater variation in expression at 95 \(\mu mol\) compared to 948 \(\mu mol\). Contrast this to the deterministic model in which the low/high concentration plots are identical with respect to a scaling factor. Gillespie causes step-behavior in the concentration. Although the plot of 948 \(\mu mol\) appears to be much smoother, we would see the same “cliffs” that we observe at 95 \(\mu mol\) by stretching the y-axis enough. The impact of a fixed gain or loss in concentration is much more “felt” in a small system than it would be in a large system (here size is relative to the initial concentration of Yf1). If we only cared about modeling a large, stable concentration however, the deterministic model clearly suffices to capture the general behavior and is less computationally expensive.

1. Starting concentration of 95 Yf1 molecules (Gillespie)


yDD (Dark-Dark state); yDL/LD (Dark-Light state); yLL (Light-Light state)

2. Starting concentration of 948 Yf1 molecules (Gillespie)


yDD (Dark-Dark state); yDL/LD (Dark-Light state); yLL (Light-Light state)


3. Starting concentration of 95 and 948 Yf1 molecules (deterministic)



(Left) Starting concentration of 95 Yf1 molecules; (Right) Starting concentration of 948 Yf1 molecules
yDD (Dark-Dark state); yDL/LD (Dark-Light state); yLL (Light-Light state)

5. Frequency of Yll at Starting Concentration of 948 Ydd at 1 umol of Light


6. Frequency of Yll at Starting Concentration of 948 Ydd at 10 umol of Light

7. Frequency of Yll at Starting Concentration of 95 Ydd at 1 umol of Light

8. Frequency of Yll at Starting Concentration of 95 Ydd at 10 umol of Light