Team:Peking/Modeling/Analysis

From 2010.igem.org

Better visual effects via FireFox ~~~




   Analyses and Results


         Modeling > Analysis and Results
We get the value of the characters (r mainly) after calculation under each set of parameters, then we analyze them trying to find one simple topology that can function well.

Identifying Minimal IOA Networks

Here we define Q value as the number of IOA function circuit among 10000 sets of network parameters, and it indicates the robustness of one topology to finish the IOA function—the larger Q value is, the more robust the topology is. We sort in reverse sequence all the network topologies according to their Q value and the x axis is their rank( Figure 6).We can observe that most network topologies have 0 or low Q value while there’s only a small part of the topologies having large Q value.

We firstly analyze the first 160 network topologies ( Q >=705 ) and list in Figure 7 all the simplest topologies that have only 3 or less direct links between the three nodes.( Figure 7) Among the 14 topologies, there are 12 three-link networks and 2 two-node networks so that we can see that the minimal number of links for the topologies to be functional is two but the most usual number is 13. The common features of the networks capable of IOA are either one negative control loop (NCL) or one negative feedback loop (NFL). Here we define NCL as a topology that has one negative control on the input-receiving node (A) from the intermediate node (B) and one positive regulation on the output node (C) from A node.(as the 1st topology in Figure 7 and the 1st one in Figure 8) Similarly we define NFL as a topology that share the positive-regulation link from A to C, while uniquely have one negative feedback from C to A.(as the 2nd topology in Figure 7 and the 2nd one in Figure 8) And the NCL topology seemingly more robust than the NFL topology, as there are 9 topologies out of the 14 simplest networks contain NCL topology compared to 6 of NFL, so that B node appears to be important, which will be discussed in details in the following part.

Qlinear.jpg


Figure 6 The Q value network topologies SequenceWe sort in reverse sequence all the network topologies according to their Q value. X Axis is their ranks, while Y Axis is their corresponding Q values. The figure indicates that most network topologies have 0 or low Q value while there’s only a small part of the topologies having large Q value.

LINEAR1.png
LINEAR2.png


Figure 7 All Functional NetworksThe numbers below each network are their ranks. The two bigger topologies have 2 links while the other 12 networks have 3 links. In each network, the green arc with one short straight line at one end stands for repression from the start node to the end node and the red arc with one arrow at one end stands for activation.

LINEAR REPRE.png


Figure 8 The NCL and NFL topologyThe green and red line have the same meaning as in Figure 7.
[TOP]

Mechanisms of Minimal IOA Networks and Key Parameters Analysis

Aimed at answering the question why the two topologies defined above (Figure 8) is functional in IOA, we unravel their mechanisms using the ODE equations in this part, also getting the parameter restrictions of each topology.

NCL Topology

When the network has built steady state:

NFunction11.png

Solve the equations:

NFunction31.png

If there is NFunction41.png, then NFunction51.png

Because 0001.png is constant, the coefficient of I is also constant. So there comes the linear correlation between 0000.png(output concentration) and I.

When NFunction41.png ,
NFunction61.png

And the concentration of B node is actually part of the slope coefficient the linear equation, we can imagine that the function of B node is to lower the concentration of A straightly without the interference of other factors as well as control the concentration of A more precisely and more freely to make the parameter restriction easier to achieve and at the same time the output range is not too small. As we know, the stochastic error may make vague the linear relationship when the values of y axis are too near. The lower concentration 0002.png is, the steeper the line is, and so the bigger the range is, which in biology means that the bioreporter is more sensitive to certain environmental signal. Through modifying the parameters of node B, we can get a proper concentration of A node to achieve a good r. In all, the node B is a proportion node.

NFL Topology

NFunction21.png

Also solve the equations:

NFunction71.png

If0NFunction81.png,thenNFunction91.pngis constant;

And ifNFunction41.png,thenNFunction101.png,which correlation is the same as the NCL Topology.
Compare X1 here with the concentration of node A in NCL topology, one can easily discover that X1 here is higher and more dull, which makes it more hard to get a balance between satisfying the second parameter restriction and the need for the range to be rational. Then we can understand why NFL topology is less than NCL in high-Q-value topologies.

Key Parameter Analysis

Aimed at unraveling the relative importance among parameters, we got all functional sets of parameters of NCL&NFL and selected randomly one set, respectively, from the two topologies and processed them with Matlab to compare the differences once the parameters change. The results are in Figure 9. From changes of the two important characters – the output range and r with the change of parameters (Table 2), we got the answer to previous questions.

Finaledition.png


Figure 9 Analysis for key parametersWe got all functional sets of parameters of NCL&NFL and selected randomly one set, respectively, from the two topologies and processed them with Matlab to compare the differences once the parameters change. When analyzing one parameter, we only change this very parameter and keep others the same, and when we change the parameter to a lower level, we get the blue line, when to a higher level, we get the red line and the black line is for the unchanged parameter set. Each line has its Pearson Correlation Coefficient r marked in the figure. The X Axis is the concentration of Hg ion as INPUT whose range is 1 to 10000 nM, and the Y Axis is the concentration of node C with the unit of nM.(A) We analyze for both topologies. (B) for NCL and for NFL (C) The degradation rate of A ( )for both topologies (D) The degradation rate of B( )for both (E)The degradation rate of C( ) for both (F) The dissociation rate of A to the binding site on the gene of C( ) for both (G) for NCL and for NFL.

For NCL, the degradation rate of A (002.png) and the value of Kd from node A to node C ( 005.png) are relatively more important to r, and the degradation rate of C (004.png) along with Kd from node A to node C ( 005.png) are essential to the output range, while the constant indicating repression from B to A(01.png) , the degradation rate of B (003.png) as well as Kd from node B to node A (006.png) influence little both characters. From Table 2 and the analysis above, we can draw the conclusion that in practice, the parameter 00.png should be larger and 004.png should be as small as possible. While for NFL, the condition is similar. It is still 002.png and 005.png that are important to r, but the key parameters on the output range are 005.png002.png00.png004.png, more than the key parameters for output range for NCL. So the suggest to practice is still smaller 004.png and larger 00.png, what’s more, there should be smaller 001.png for larger output range according to the character changes in Table 2.

              Table 2 The change of two important characters as the parameters rise

Table2222.png

[TOP]

Parameters Tendency of the IOA Networks

In order to test the correctness of our analysis about the key parameters, we look more closely at the parameters of the simplest IOA networks. For the two simplest networks in Figure 6, we examine the distribution of the parameter values which can realize linear response between input and output out of the total 10000 sets of parameters used in the search.(Figure10, Figure11)

Kbc.jpg


Figure 10 Distribution of NCL topology parameters which can establish linear response curves

Kac.jpg


Figure 11 Distribution of NFL topology parameters which can establish linear response curves.

It can be seen that the constant K of the repression (006.png and 007.png) are generally lower compared to 005.png, which means that there should be an intense repression on node A. Such a result also fits well with the conditionNFunction41.png used in the derivation above.

We also notice that the degradation rates, α, of the nodes also show apparent tendency. The degradation rate of the node A is generally higher, while the repressing nodes (B or C) have lower degradation rates. Such a tendency also shows that the repression on node A should be intense in order to establish a more linear response curve between input and output.
[TOP]

Analysis of All Possible Three-Node Networks

The above analyses focused on minimal( less than or equal to 3 links) three-node networks and identified simple topologies that are sufficient for IOA function, also unraveling the mechanism that the topologies work. But whether the topologies are necessary for the IOA function is not understood yet. In other words, are the identified minimal topologies the foundation of all possible networks, or are there more complex higher-order solutions that do not contain these minimal topologies? To answer the questions above, we analyze the first 160 topologies (Q>705) that are well capable of the IOA function.(Figure 12)

WenshiLINEAR.png


Figure 12 Analysis of the first 160 networks We count all the NCLs and NFLs and also the IOA networks, and discover that all of the IOA networks can be classified into NCL/NFL/the combination of the two. And there are more IOA functional network featured in NCL than those characterized with NFL.

Analysis of these robust topologies shows that they are overrepresented with NCL and NFL. All 160 topologies contain at least one NCL or NFL motif ( or both ). These results indicate that at least one of these motifs is necessary for IOA function.

Supplementarily, the NCL average Q value(AQV)of all 19683 topologies is 17.14 while the NFL AQV is only 9.36, which again indicate that the NCL is more robust than NFL that have drew conclusion in the minimal topology analysis.
[TOP]

Motif Combinations that Improve IOA

To investigate what additional features can improve the functional performance in some more complex and more robust networks than minimal topologies, we clustered the first 160 networks and then cluster them respectively in three categories: NCL, NFL and the combination of the two. (Figure 13)

JuleiLINEAR.png


Figure 13 The clustergrams of the networksWe use the clustergram command in matlab to get the additional features of the functional networks. The nine vertical rectangle bar stand for nine links in Figure 2 which are, respectively, from A to A, from A to B, from A to C, from B to A, from B to B, from B to C, from C to A, from C to B, from C to C. And red stands for activation, green for repression and black for no regulation. The topologies on the right are corresponding minimal topologies that is shown in the clustergrams on the left.

The results clearly indicate that apart from the link from A to C, there should be no positive regulation, and the NFL topology hates the link from A to A, while the combination topologies show no additional tendency.
[TOP]

Summary

In all, by exhaustively searching all network and analyzing the results, we draw the conclusion that NCL meets our need in application well, and we get the proper parameter range for practice.

[TOP]

go back to top