Team:GZHS-United/Model
Modeling
Background
B.s and B.ti are both used as bio-pesticide. Some researches had shown that Cry4Ba and Mtx1 have synergism when they are working in the same system, which means that if we mix or fuse the two toxins together, the pesticide will be more effective than using them separately. It’s important to find out whether the synergism exist and then estimate the best ratio of CRY4Ba and Mtx1 in the mixed system. In order to achieve this goal, we introduced probit model to analyze the data we’ve got in toxicity test and then calculated the best ratio of these two toxin by curve fitting.
Probit model
In statistics, probit model is a type of regression that the dependent variable can take only two values like “dead or not dead”. The purpose of the model is to estimate the probability of an observation with particular characteristics falling into a specific one of the categories.
Let’s suppose a response variable Y is binary, which means that it can only have two possible outcomes, 1 and 0. We also have a vector of regressors X, which are assumed to influence the outcome Y. Specifically, we assume that the model takes the form 
, where Pr denotes probability, and Φ is the Cumulative Distribution Function (CDF) of the standard normal distribution. The parameters β are typically estimated by maximum likelihood.
Estimation of toxicity
There are several index to demonstrate the toxin toxicity. Here we calculated KT_50 to measure it. KT_50 means the time a kind of toxin cost to kill half of the total number of insects in an independent system. Lower KT_50 indicates higher toxicity.
After succeeded in expressing cry4Ba in E.coli, we designed our experiment for modeling. In our experiment, every independent system contains 50mL ddH2O, 2.5mL bacterial(induced cry4Ba E.coli and Bs ) and 20 subjects of 2nd-3rd instar larvae. The ratio of cry4Ba in the mixed toxin ranged from 0 to 100%. Each group had three repetitions, and all the groups were observed for 48 hours. We recorded the number of death in 2h, 4h, 6h, 8h, 12h, 24h, and 48h. Based on the data we got in the experiment, we used probit model to estimate KT_50.
| Group | Volume of cry4Ba E.coli (mL) | Volume of Bs (mL) | Ratio of cry4Ba E.coli(%) |
|---|---|---|---|
| 1 | 0 | 2.5 | 0 |
| 2 | 0.4 | 2.1 | 16 |
| 3 | 0.8 | 1.7 | 32 |
| 4 | 1.2 | 1.3 | 48 |
| 5 | 1.6 | 0.9 | 64 |
| 6 | 2.0 | 0.5 | 80 |
| 7 | 2.5 | 0 | 100 |
| Time | 2 hr | 4 hr | 6 hr | 8 hr | 10 hr | 12 hr | 24 hr | 48 hr |
|---|---|---|---|---|---|---|---|---|
| 1 | 0 | 0 | 0 | 1.33 | 1.33 | 1.33 | 2.67 | 2.67 |
| 2 | 0 | 0 | 0 | 1.67 | 3 | 4.33 | 7 | 8.33 |
| 3 | 0 | 0 | 0 | 0.33 | 2.33 | 3.67 | 9.67 | 16 |
| 4 | 0 | 0.33 | 1 | 1.67 | 4 | 5.67 | 11.67 | 16.33 |
| 5 | 0 | 0 | 0 | 2.33 | 2.67 | 5 | 14.67 | 17.67 |
| 6 | 0 | 0 | 0 | 0.33 | 1.33 | 3 | 13.33 | 19.67 |
| 7 | 0 | 0 | 0 | 0 | 1.33 | 1.33 | 16 | 20(ALL) |
KT_50
We used SPSS to calculate KT_50. Response frequency is number of death. Covariate is time. Total observed is 15. When probability=0.5, the corresponding “Time” is KT_50. Data is shown in figures below.
Figure 1. Probit regression model and KT_50 data from the model.
Curve fitting
We’ve got KT_50 data in probit model, so the next step is to figure out whether there is a best ratio of the two toxins when they’re used to kill Aedes albopictus larvae .
| Percent of CRY4Ba(%) | KT50(h) |
|---|---|
| 16.00 | 45.53 |
| 32.00 | 25.22 |
| 48.00 | 21.03 |
| 64.00 | 18.86 |
| 80.00 | 19.54 |
| 100 | 20.43 |
Because the KT_50 of group 1 (ratio of cry4Ba is 0% and ratio of Bs is 100%) is too high and may interfere our fitting curve, we decided to abandon this data. Data we used to do the curve fitting is shown in table 3.
According to the result of curve fitting, if percent of Cry4Ba is independent variable, and KT_50 is dependent variable, we can get the equation:
y=0.008x^2 - 1.17x + 59.19
it’s easy to calculate that when y=min, x=73.13. So the best ratio of Cry4Ba is 73.13%
Discussion
Our model illustrated two significant points:
- Bs and cry4Ba E.coli actually have synergism, so it is possible that protein Cry4Ba and MTX1 also will exhibit synergism.
- The toxin mixture could kill mosquito species Aedes albopictus. Also the ideal ratio of cry4Ba E.coli in the mixed toxin with Bs is 73.13 % (target to Aedes albopictus).
Contribution to our project
This model has two important contributions to our project:
- According to the modeling results, we decided to add protein MTX1 into our design.
- We optimized the product of our entrepreneurship: pesti-jelly based on this result.
