Simulation and analysis

In this modeling step, we will develop a study about different chaperones arrangement efficiency in the proteins folding, using as example our protein of interest, Limonene Synthase. As it was not possible to obtain the complete system of this toolkit, there are no available experimental data to reinforce our findings. Nevertheless, the modeling of this section is extremely significant and could be useful to another teams with some similar problems. In this sense, we decided to model our system starting from a simulated dataset where we use fake data, and in the future, real data could simply replace them.

First of all, the experiment to measure the efficiency of each chaperone and arrangement of them was conceived in the following form:

All chaperones would be tested alone (IbpA+B, ClpB and DnaK), in pairs and all together. The proteins would be encoded in plasmid pSB1C3 under control of a constitutive promoter, allied to our gene of interest under the same promoter. Protein of interest would be produced by some hours, and the cells would be harvested and lysed by ultrasound. Finally, the cell debris would be precipitated by centrifugation and the soluble and insoluble portions would be recovered. After, a SDS-PAGE the band of interest would be quantified with optical densitometry. Chaperone efficiency coefficient would be given as the reatio between soluble to insoluble bands.

Simulated data following this procedures were provided in the table bellow. A control group without any chaperones was included to comparision, and all measurements were made in triplicates. This table shows portion of well folded proteins (soluble) and misfolded ones (insoluble).

Table 1: Simulated data using different chaperone sets and the solubility after the treatment.

Chaperones efficiency coefficient still means the proportion of correct and incorrectly folded proteins. Considering the case of 100% of solubility, this coefficient tends to infinity, since the insoluble portion tends to zero. In this sense, a case of this type would be impossible to work, to prevent this, we decided work with protein yield. Protein yield coefficient would be given as percent (0-100) and calculated as follows:

(Eq. 1) $$Y = 100*[Soluble/(Soluble+Insoluble)]$$

Table 2: Chaperone arrangements related to chaperones efficiency and protein yield coefficients.

In order to observe this large set of combinations relate to each other in terms of yield, we used a hierarchical analysis model of clustering available on a free webserver, DENDROUPGMA. It has done the calculations used in the creation of dendogram shown in Figure 1. This analysis basically shows together those data that are mathematically closer and separates those who are further away. As we can see, the values obtained with all chaperones acting together, here named "all", are more isolated from other combinations for being the most efficient treatment of the others and so it is statistically more 'distant' to the other groups. Similarly we see the double "Ibp+clpB" and "Ibp+DnaK" grouped in the same node, so we conclude that their values are closer enought to not differentiate them and therefore cannot be separated.

Figure 1: Dendogram of distances using simulated dataset generating a clustering complex behaviour.

Cophenetic correlation coefficient, which brings how feaseable the dendogram is, was estimated as 0.84134 and values higher than 0.7 are considered as a good-fit. It can be seen that the combination of "all" chaperones together is farther away from the others, which was expected since its yield reached the highest value. Arrangements between "Ibp+DnaK" and "Ibp+clpB" are together on a branch, indicating that their yield values are closer than others. In addition, they are nearby branch of combination of all chaperones, indicating its yield not as high as the performance of all chaperones together, although close enough. This may indicate combinations "Ibp+DnaK" and "Ibp+clpB" as good substitutes for all three chaperones, when it would not be possible the complete arrangement.

As possible to observe, the protein yield using chaperones arrangement are not trivial, since there will be cross interactions hard to predict its behaviour jointly taking in account the single values. Predicting the interference between chaperones system, we equate all possible combinations of yield corresponding to each combination and multiply them by a mathematical factor \ (x_i \). The goal is to interpret the values of these constants, to make it possible to compare the influence that each chaperone component isolated influences on the set of combinations.

Each index \(x_i\) represents an arrangement of all studied chaperones, in this way: \(x_1 \rightarrow Ibp\); \(x_2\rightarrow Clbp\); \(x_3\rightarrow DnaK\); \(x_4\rightarrow Ibp+Clbp\);\(x_5\rightarrow Ibp+DnaK\);\(x_6 \rightarrow Clbp+DnaK\); \(x_7 \rightarrow Ibp+DnaK+Clpb\).

Results obtained for this system were:

(Exp. 1) $$x_1=(471/124)x_7-(251/248)x_6 $$

(Exp. 2) $$x_2= -(314/71)x_7 +(502/213)x_6$$

(Exp. 3) $$x_3=(471/109)x_7 - (251/218)x_6$$

(Exp. 4) $$x_4=(251/612)x_6 $$

(Exp. 5) $$x_5=(1884/503)x_7-(502/503)x_6 $$

(Exp. 6) $$x_6=x_6 $$

(Exp. 7) $$x_7=x_7 $$

In this way, it will be possible to describe all parameters taking in account just two values \(x_6\) and \(x_7\). There is a relationship between all chaperones and the two combinations (all of them and “DnaK+Clpb”. Fixing one of them, it will be possible to study the behaviour of others using graphical tools (Figure 2 and 3).

Figure 2: KKK.

Figure 3: KKK.

It was possible to see that there is no single solution to our problem, and probably even with actual values perhaps there is no single solution. Although, the graphics can nonetheless take some dependency relationships between the coefficients, such as if \(x_7 \) is fixed \(x_2 \) component increases with slope greater than \(x_1 \), which decreases with the increasing of \(x_6 \) value. This indicates that ClpB influence in the chaperone combinations are higher than Ibp in these arrangements. In this sense, this result is reversed to from the moment when increase the amounts related to \(x_7 \). As it was not possible to obtain experimental data to validate our discussions and predict outcomes, we are pleased to leave this modeling as a theoretical model that can be deployed and tested in future projects of our or another iGEM team that wants to employ it.

References

S. Garcia-Vallve, J. Palau and A. Romeu (1999) Horizontal gene transfer in glycosyl hydrolases inferred from codon usage in Escherichia coli and Bacillus subtilis. Molecular Biology and Evololution 16(9):1125-1134.