Table. 1 Data from experiment
*Note (1): “time” is the number of hours elapsed since the first measurement; “B1, B2, B3” is the OD measured for each sample; “OD(A)” is the average of them; “OD(A)*5” is five times of OD(A); “cell/ml” indicates how much cells per ml according to OD(A)*5, using the relationship between them:P_((cells/ml))=(133.16×OD(A)*5-1.43)×10^5; “P0” equals to 10^7 cells/ml, “P’” equals to P/P0 ; “t0” equals to 24 h, “t’” equals to t/24.
*Note (2): The reason why we set P_0 and t_0 is that when unifying units, we chose cells/ml; and for convenience, the data inserted into the computer model should not be too large, otherwise errors or running timeout may occur. Therefore, we chose to dimensionless the data. That is, P=P^'×P_0, where〖 P〗_0=10^7 cells/ml, which is a quantity with the unit. Thus, P^' is only a number, without a unit. For “time”, it is similar.
Since the data is dimensionless, the function will change a little bit into:
Therefore, the data that came back from the model was actually r' and K'. To find exact value for r and K,
Then, we used data of P’ and t’ to fit a logistic growth model, based on codes downward,
Fig.1 Codes, Results, and the Plot of logistic growth model (P’~t’)
By entering (1), (2) using codes, the regression equation that we got is:
Then, we could use the 3 parameters a_(1,) a_(2,) and a_3 to express carrying capacity K, and the constant of proportionality r, that is,
In the third picture “alpha3” line, “4.91” is the carrying capacity K'; the constant of proportionality r'=(Named num3)/(Named num1)=8.18/4.91≈1.666.
To conclude, by using data collected from experiment, we fit them with a logistic growth model and found the carrying capacity K=4.91×10^7 cells/ml; and the constant of proportionality r≈6.94×10^(-9) ml/ (cells∙h). These parameters will be used later to test the algicidal efficiency of our drug.
Fick's Second Law for the peptide diffusion:
Fick's second law states that in an unsteady diffusion process, the rate of change of concentration with time at distance X is equal to the negative rate of change of diffusion flux with distance at that point:
The codes are:
On this basis, we construct a partial differential equation with the boundary condition that the value of C is 1 when x is 0 or 1, so as to predict the diffusion of a certain water area along the axis after the input of a specific concentration of polypeptide.
To solve this function, We Pick a step h=1/N, where N is a positive integer. Then we use the approximations
We replace our function u by a discretizationu_k=u k/Nand the derivative u'' becomes -2u_i+u_(i+1)+u_(i-1.) Then translate it into a system of equations.
Then we introduced a time axis to describe the change of polypeptide concentration along the axis over time：
In Neuman boundary condition, which is
According to this, we can get
Then we just modify the system of equations, which is formula (3).
We can then use these equations to plot how the concentration of peptides changes over time
Fig.2 the concentration of peptides changes over time
To conclude, in this section, we predict the diffusion of polypeptide in a certain water area along the axis after the input of a specific concentration of it. From the result, we can learn the dispersal, and determine the concentration and the time interval to make sure that the concentration is above the minimum concentration that can kill algae within a certain range.
Algicidal efficiency model for the drug:
In this model, we want to find the algicidal efficiency. During the wet lab, we added different doses of the drug to three groups of identical samples. Then, we measured the OD of every sample about every 24 hours. Here is the graph that indicates the OD measured of algae after different doses of drug is added:
*Note (1): At Date “8/17” (the 144.5th hour), the drug is added into three samples (C3 500(10μM), C3 950(20 μM), C2 50(1 μM)).
*Note (2): “C2 B(no drug added)” is the OD of sample that no drug is added, which is a control sample; “C3 500(10μM)” is the OD of sample where the final concentration of the drug is 10 μM; “C3 950(20 μM)” is the OD of sample where the final concentration of the drug is 20 μM; “C2 50(1 μM)” is the OD of sample where the final concentration of the drug is 1 μM.
*Note (3): The record showed here closed on that day when C3 950(20 μM) reached its lowest point.
Fig.3 OD of different samples.
It is obvious that when the concentration of the drug in the solution is 20 μM, the algicidal efficiency is highest. The killing efficiency is about 51.57%. Therefore, we concentrated on this group of data and continued with our modelling to find the algicidal efficiency.
We firstly set linear differential equations below:
The first function indicates the growth of algae after the drug is delivered. The second one indicates the decrease of drug after it is put into the water area. P is the concentration of algae when t>0; r is the constant of proportionality for logistic growth of algae, which was found to be 1.666; M is the carrying capacity for logistic growth of algae, which was found to be 4.91; β is the algicidal efficiency of the drug; D is the concentration of drug when t>0; γ is a constant that indicates the durability of a drug.
Since r and M are known quantities, and the initial data for P and D (P_i and D_i) can also be recorded from wet lab, then if we enter a value for β and γ, we will get a value for C, which can be compared with the value we recorded through wet lab. Then we can find the magnitude of their difference, that is, |P_measure-P_modelling |. If we can find the exact value for β and γ that could minimize |P_measure-P_modelling |, then we could correctly find β and γ.
Again, we dimensionless the data, and the function turns to:
So, for the real value of β and γ, it should be:
We choose the last four points of C3 950(20 μM) to do the fit, that is, to be used as P_measure at certain time points, and after we unify the units and dimensionless the data, we have:
Table.3 The last four points of C3 950(20 μM)
The initial values are:
Here are the codes:
Fig.4 Codes for linear differential equations, parameters, and loss function.)
Fig.5 Plots for P (left) and D (right).
And then, we plotted the function of P and D, and get:
Based on the results of the run, P_measure and P_modelling fit well, so the results are reliable. From this model, we got:
Therefore, we got:
To conclude, in this section, we predicted the algicidal efficacy and duration of the drug. The algicidal efficiency is found to be around 51.57%, to be more specifically, 2.54×10^(-15) cells∙h/ml, and durability of a drug is found to be 3.12×10^(-9) cells∙h/ml. Combining the results of the peptide diffusion model and the algicidal efficiency model, we can obtain predictive guidance for the delivery of drugs in a certain water domain, such as the concentration of delivery, delivery interval, delivery area, etc., to achieve the best algicidal effect.
© Copyright DKU iGEM. All Rights Reserved