Crosstalk Circumvention Modeling

1. Introduction

In order to achieve Signal-dependent state-change circuit with crosstalk circumvention, we designed ‘’Ninja circuit’’ which is shown in Fig. 4-1-1.

Fig. 4-1-1. Our designed ‘’Ninja circuit’’

We analyzed the dynamic characteristics of the whole circuit through the modeling of “Ninja” circuit. In addition, we verified how each parameter in the model formula affects the behavior of the circuit.

To understand the difference between our Ninja model with circumvention system and the original circuit, which is not able to prevent from signal crosstalk, and, we compared with two gene circuits on equal terms.

2. Analytical method of dynamic characteristics of the circuit

2-1. Building a mathematical model

On the basis of the Hill equation, we built a mathematical model of the Ninja circuit. The mathematical model which was used this time was as below. Value for each parameter will be discussed on Chapter 2.3.

Fig. 4-1-2. Equations for our modeling

We will divide these equations into several groups to explain the mathematical model.

The first group describes change of the signaling molecules (Fig. 4-1-3). At this time, because of the constant expression of LasR and LuxR, only degradation of C6 and C12 are considered. They are added as input in a pulse manner in this scetion.

Fig. 4-1-3. Equation of degradation

The second, in Fig. 4-1-4, is differential equations of LacI. Since we use Promoter TetR that acts as a repressor, the first term shows the suppression of TetR. Although hybrid promoters which express LacI are activated by C6, we have to consider the influence of the crosstalk, which is represented as addition. Because our wet experiments shows that expression by C12 input, which cause cross talk via LasR activation, is about half of expression by C6 input, we set a3 as half of a2. About equation (3), not only from the hybrid promoter, LacI is also expressed from one side of the toggle switch. For this expression, we add a Hill equation to show suppression by CI. In the end we add expression amount of the leak, then minus the decomposition item of LacI.

Fig. 4-1-4. Equation of LacI

As the third, CI is expressed not only from Plas promoter, but also from one side of the toggle switch. As a result, we add a Hill equation to show suppression by LacI.

Fig. 4-1-5 Equation of CI

Finaly, we derived those following equations based on the circuit as well.

Fig. 4-1-6. Equation of CI434, TetR, GFP, RFP, g2p

2-2. Description of the parameters

Each parameter is defined as the following table.

Fig. 4-1-7. Description of each parameter

2-3. Parameter optimization and simulation

We carried out the simulation of this time by using MatLab. Since the amount of the parameters of the differential equation is very large, we tried to minimize the amount of the evaluation function. (Fig. 4-1-7)

Fig. 4-1-8. Evaluation function

As a result of exploring the a4, a5, a6, and a7 to minimize the evaluation function in equation (10), we found out that it would be the best selection if we take the parameters as follows. (Fig. 4-1-9)

Fig. 4-1-9. Value of each parameter

2-4. Simulation result and parametric evaluation

Bistability by the toggle switch in Ninja circuit

We firstly confirmed bistability by toggle switch. In initial state, by pulse addition of 20 nm C12, the CI side get induced and maintained even after decay of C12. After 300 min, with pulse addition of 20 nm C6, the toggle switch starts to switch to LacI “mimic” state. In addition, the state is maintained after decay of C12. Thus we could get the conclusion that the LacI state stays stable after the switching.

Fig. 4-1-10. Switching from CI to LacI

The following graph has been obtained by using the above parameters. In initial state, by the addition of 20nm C12, the CI side get induced、after 300min, with the addition of 20nm C6, the toggle switch starts to switch. We could get the conclusion that the toggle stays stable after the switching.

Fig. 4-1-11. Switching from LacI to CI

Crosstalk circumvention by TetR

After we confirmed the bistability by the toggle switch, we analyzed the switching dynamics.

From the following analysis, we modified equations (1) and (2) to describe presence of either E. samurai or E. civilian. Each of them accumulates a signaling molecule during his presense.

Firstly, we consider the situation that when E. civilian comes to attack state. Fig. 4-1-12 shows the switching from attack state to civilian state which is influenced by C6 came from E.civilian. We can know from this simulation result that when E.civilian comes on the time of 300 min, CI attack state will switch to LacI civilian state.

Fig. 4-1-12. E.civilian comes (With TetR)

During the switching from the attack state to the mimic state, absence of TetR allows activation of Pluxtet hybrid promoter. Repression of TetR production, by CI434, is indeed important in the circuit. During the mimic state, TetR accumulates to plateau level. This presence of TetR is important to prevent crosstalk, in the next figure (Fig 4-1-13), by LasR activated by C12.

Secondary, we consider the situation that when E. samurai comes to the mimic state. When switching occurs from mimic state to attack state, CI is expressed. Note that LacI expression form the Pluxtet hybrid promoter is prohibited, due to the presence of TetR, even in the presence of C12-LasR complex which can bind to the hybrid promoter for its activation. Abundant presence of C12AHL expressed from E. samurai causes CI overexpression, which stimulates expression of gp2 required for the programed phage release. Though degradation of TetR during the presence E. samurai causes LacI expression from the hybrid promoter, CI expression from Plas promoter overcomes the effect of LacI. After C12 degradation, the model shows CI-concentration decrease whih stops “shuriken” throwing in our scenario.

Fig. 4-1-13 E. samurai comes (With TetR)

Interestingly, CI expression oscillates and converges by C12AHL induction (see time point 400-1000) (Fig. 4-1-13, which is a magnified graph of Fig 4-1-12).. This is because not only the toggle switch, but also there is a repressilator by combination among TetR, LacI and CI434 which suppress with each other. As adding C12 from outside, TetR oscillates and decays. According to the effect, CI also vibrates and be settled down to a constant value. Moreover, the amount of CI decides the expression level of g2p, so we can see g2p vibration, too.

Fig. 4-1-14. Switching from caution state to attack state

3. Analytical method of effectiveness of crosstalk prevention circuit

3-1. Model without crosstalk prevention circuit

Since we already have the simulation result of, to verify the difference of “Ninja” circuit and the naive “Signal-dependent state change circuit” which cannot prevent from the crosstalk of signal, we established the model for the naive circuit. It was defined as follows.

Fig. 4-1-15. Equation of naive circuit

Value of the parameter used here is the same as the value used in Fig. 4-1-10.

3-2. Simulation results

For the dynamics of changing from LacI state to CI state, we compared the characteristics of the both circuit. The reason why we won’t compare characteristics for the opposite changing, from CI state to LacI state, is that crosstalk won’t occur in this changing. We analyzed the following case; the initial state with enough LacI changes to the state with enough CI, by the increase of C12 expression from E. samurai who came at 300min. Fig. 4-1-16 shows the changing of LacI and CI. The solid line presents the case with crosstalk prevention circuit, and the dotted line stands for the case without crosstalk prevention circuit. With C12 production, LacI is produced in a certain amount in the naïve circuit without crosstalk prevention. On contrary, Ninja circuit actually circumvent the crosstalk to reduce LacI concentration severely. . When E. samurai has gone at 1500 min, the concentration of LacI expressed by “Ninja” circuit is about 25nM. In contrast, the concentration of LacI expressed by the naive circuit, which cannot prevent from the crosstalk, is 725 nM. Regarding the convergence to 5 nM LacI, we can conclude that “Ninja” circuit will be converged faster than the naive circuit.

Fig. 4-1-16. Difference between original circuit and “Ninja” circuit