• Introduction
• Defining The Best Topology
• Sensitivity analysis

Model

Introduction

A regular reporter (for example, fluorescent proteins or pigment) can only detect a narrow range of concentration of input signals, (Figure 1) because they mainly possess a Hill-function type dose-response curve. The linear proportion of such a dose-response curve is relatively narrow，causing it sensitive to only a narrow range of input intensity. Thus the regular reporters are not proper devices for quantitative measurement. In order to design a quantitative, economical and convenient aromatics detector, we decided to build a band pass filter.

Selecting Network Topologies

There are a lot of potential band-pass filter topologies with an incoherent feed-forward loop as its core topology. The simplest one is a two-node topology with a positive loop and a negative loop (Fig 1.a). But it’s too difficult to make a transcriptional factor both activator and repressor at the same time.

So we set out to look for topologies with three or more nodes based on two criteria: first, considering functioning mechanisms of real protein regulators, the input node must regulate other nodes in a uniform manner, either all inhibiting or all promoting; second, to lessen the manual work required to construct the circuit, we require both the positive and negative feed-forward loop to contain no more than one internode. In the end, we decided to analyze three topologies in detail, two of them have three nodes and one has four nodes (Fig 1.b, c, d).

Figure 1.Four kinds of circuit that might form a band-pass filter. (a)Two node circuit. A isthe inducible transcriptional factor, B is the reporter GFP. (b)Three node circuit. . A is the inducible repressor, B is a repressor regulated by A, C is the reporter GFP. (c)Three node circuit. . A is the inducible activator, B is a repressor regulated by A, C is the reporter GFP. (d)Four node circuit. . A is the inducible activator, B is a repressor regulated by A, C is the activator regulated by A, D is the reporter GFP.

ODE Analysis of Circuit Networks

After determining three typical circuit networks that might function as band-pass filters, we used Ordinary Differential Equations (ODE) to analyze the three topology and their parameter preferences.

Our ODE model of the circuit networks was established based on the following assumptions:

1. Only transcription and translation process is taken into consideration. The equations only contain accumulation and degradation of the proteins represented by the nodes.
2. The interactions between nodes are described using Hill functions. We multiplied the Hill functions describing the individual regulating effects of activator and repressor in order to simulate the synergic regulating effect of the two transcription factors when bound to the same promoter.
3. The protein represented by the input node is constitutively expressed at a fixed level.
4. Leakage is present only in the expression of reporter gene (output node).

The ODE equations we used to simulate the three circuit networks are listed as following:

Circuit Network 1:

Note: Since A is expressed at a fixed level, its influence on the network can be absorbed into other constants in the model, so it doesn't appear in the equations written above.

Letting the differential term on the left handed side equal to zero, the steady state solution of the equations can be solved as following:

Circuit Network 2:

Note: For reason same as that described in Circuit Network 1, expression level of A doesn't appear in the equations written above.

The steady state solution of the equations can be solved as following:

Circuit Network 3:

Note: For reason same as that described in Circuit Network 1, expression level of A doesn't appear in the equations written above.

The steady state solution of the equations can be solved as following:

In the equations written above, k in lowercase are used to denote maximum expression rate constant of individual regulation loops. k in uppercase are used to denote the transition point of different Hill functions. A, B, C, D are used to denote expression level of proteins represented by corresponding nodes. n in lowercase denote Hill coefficient of different Hill functions. l in lowercase denote leakage terms.
The next step is to use ODE equations given above to find out the most robust circuit network. In order to do this, we applied 10,000 sets of randomly sampled network parameters to each of the three topologies and evaluated their behavior according to two important features: induction fold and peak width (Fig. ). Parameter sets resulting in a high induction fold (>32) and a reasonable peak width (<4&>0.5) are recorded separately for each topology.

Figure 1.Dose-response curve of a band-pass filter. We chosen the peak’s fold and peak’s width as output and the concentration of the ligand as the input. Peak fold is calculated as divided by and the peak width is define as the width of the peak’s half height.

Figure 2.Three kinds of circuit and the set of parameters that meet our need. (fold>32 and 4>peak width ). The blue line means repression and the blue arrow means activation. The number Q in the circuit shows the set of parameters falls into the area we choose.

Results of robustness analysis are shown in Figure 2. 20 out 10000 sets of parameters applied to circuit network 3 generated a desirable performance. The number of effective parameter sets for circuit network 1 and 2 are 5 and 13 respectively. This indicates that circuit network 3 may function in a larger proportion of the parameter space, meaning it would be more likely than the other two to work well in noisy biological context and might require less effort on fine-tuning individual components of the circuit network. So we chose the third circuit network as the blueprint for our band-pass.

Parameter Sensitivity Analysis

Having decided which circuit network to use, our next move is to determine which proteins to choose for individual nodes in the circuit network to construct a band-pass filter. First we decided to use parameter sensitivity analysis to figure out which parameters in the circuit network are most sensitive to perturbation. These parameters are crucial for own design because they must be carefully controlled within a narrow range to make sure our circuit would function properly. Once we have figured out all sensitive parameters, we may facilitate our design by choosing proteins with desirable sensitive parameter value.

Figure 3.The Circuit we choose and the equations we use to stimulate. And the equations we choose to stimulate. The pink arrow is the promoter. The green ellipse means the RBS. The blue block represents the protein. The red octagon is the terminator.

In order to be more specific in our parameter sensitivity analysis, we re-wrote the ODEs for the band-pass filter as following:

And steady state solution can be solved as following:

Symbol usages are identical to that in ODE analysis of network topologies.

We chose one set of parameters meeting our needs (peak fold>32 and 4>peak width>0.5) as our initial parameter set. We then enlarged each parameter ten times and observed the difference between pre-change behavior and post-change behavior of the band-pass filter circuit. We used log(f’/f) and log(p/p’) to characterize the change in band-pass performance, where f and f’ denote the fold change before and after parameter perturbation, p and p’ denote the peak width before and after parameter perturbation. A log(f'/f) greater than zero indicates that the fold change has increased after the parameter perturbation, making the perturbation desirable. A log(p/p') greater than zero indicates that the peak width has reduced after the parameter perturbation, also making the perturbation desirable.

The quantitative results of parameter sensitivity analysis are shown in Figure 4 and the qualitative results are summarized in Table 1. To further help understanding we drew a series of pictures to show the how the network behaves as some of the parameters change (Fig.5).

Figure 4. Fold’ (f’) and p’ indicate Fold and peak width after we changed the parameters and Fold (f) and p indicate Fold and peak width before we changed the parameters. An upward-pointing column indicates changes that are conducive to our design, and a downward-pointed column indicates changes detrimental to our design. This figure shows how the band changes with the variation of each parameter.

As can be concluded from figure 4 and table 1, parameters are sensitive to enlargement in four different ways. Some parameters, such as ksa, nb, na', Kag nb' and Kbg, once enlarged, will enhance or deteriorate the performance of induction ration and band width at the same time; another group of parameters, such as ksa, na and ksb, once enlarged, will enhance one important feature while deteriorating the other; a third group of parameters, such as kag and kbg, have little effect on one of the two important features; a final group of parameter, consisting of Ksb only, will totally eliminate the band when enlarged.

In order to tune our band pass filter most efficiently, we decided to focus on the first group of parameters because by increasing or decreasing the value of such parameters will improve the induction ratio and peak width at the same time. And since in real applications we need to change the sensor (input node) from one aromatics-responsive transcription factor to another, we don't want to set too tight a restriction on sensor-related parameters such as ksa and nb. So we limited the tunable parameters to na', nb' Kag and Kbg. As can be easily observed from figure 4, for a outstanding band-pass filter circuit, higher values for na' and nb' as well as lower values for Kag and Kbg are desirable.

Figure 5.The dose-response curve of different sets of parameters. na’ is the hill coefficient of repressor A. KAG is the ligand concentration producing half occupation. The green line means na’=2, the yellow one means na’=1.5, the red one means na’=1 and the blue one means na’=0.5. (a)KAG=1.5(b) KAG=1(c) KAG=1.From the pictures above, we can see that a small KAG and a large na’ is needed to build a good performance band-pass filter.

Table 1. The change of two important characters as the parameters rise. Green arrow has the meaning that the character value rises as the parameter rises. Red arrow has the reverse meaning. Black line means that the charater value keep relative stability as the parameter changes.

Based on all the analysis we have done above, we chose the activator from phage phiR73 as node B for its small dissociation constant (Kbg) and repressor form phage λ as node A for its large hill coefficient (na') to construct our final band-pass filter circuit (Fig. 5).

Figure 6.The final construct of our band-pass filter.