-
Notifications
You must be signed in to change notification settings - Fork 0
/
NottinghamPhageN3APIIVIODE.m
89 lines (75 loc) · 3.55 KB
/
NottinghamPhageN3APIIVIODE.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
function dy = NottinghamPhageN3APIIVIODE(t, y, sP)
% ODEs for a prey and two predators system
%
% Equations represent a batch system with an abiotic resource, prey species
% demonstrating monod kinetics, some of which show resistance vs predation,
% others which develop or lose persistance to predation and two predators
% one exhibiting a Holling type I functional response, the other a Holling
% type II response and mortality. Each predator has a seperate
% bdelloplast / infected cell stage. Development of persistence is %
% proportional to dead prey killed by bdellovibrio.
% Mode 23
%
% function dy = NottinghamPhageN3APIIVIODE(t, y, sP)
%
% dy - the rate equations
%
% t - time series
% y - initial species values
% sP - simulation parameters
% Version Author Date Affiliation
% 1.00 J K Summers 14/12/17 Kreft Lab - School of Biosciences -
% University of Birmingham
% 1.01 J K Summers 08/10/18 Ongoing developement as well as initial
% phage resistance
%
% Initial numbers of individuals in species being modeled
sub = y(1); % substrate value in fg/ml
senPrey = y(2); % Sensitive cells in cells / ml
bdPersistPrey = y(3); % Bdellovibrio persistor cells in cells / ml
phageResPrey = y(4); % phage resistant cells in cells / ml
bd = y(5); % Bdellovibrio cells in cells / ml
bdplast = y(6); % Bdelloplasts in cells / ml
phage = y(7); % Phage virions in virions / ml
infCell = y(8); % Infected cells in cells / ml
deadPrey = y(9); % Dead prey cells in cells / ml
% Rates for processes
dSub = -(senPrey + phageResPrey + bdPersistPrey * sP.persGrowth) * ...
sP.muMaxPrey * sub / ((sP.Ksn + sub) * sP.yieldNPerS) + ...
sP.yieldSPerV * sP.kV * infCell + ...
sP.yieldSPerP * sP.kP * bdplast;
dSenPrey = senPrey * sP.muMaxPrey * sub / (sP.Ksn + sub) - ...
bd * sP.muMaxPred1 * senPrey / ...
((sP.Knp + senPrey + bdPersistPrey + phageResPrey) * ...
sP.yieldBPerN) - ...
phage * sP.muMaxPred2 * senPrey / sP.yieldIPerN - ...
senPrey * sP.rateDevPersist * deadPrey + ...
bdPersistPrey * sP.rateRevert - ...
senPrey * sP.rateDevResistance;
dBdPersistPrey = bdPersistPrey * sP.muMaxPrey * sP.persGrowth * sub / ...
(sP.Ksn + sub) - ...
phage * sP.muMaxPred2 * bdPersistPrey / sP.yieldIPerN + ...
senPrey * sP.rateDevPersist * deadPrey - ...
bdPersistPrey * sP.rateRevert;
dPhageResPrey = phageResPrey * sP.muMaxPrey * sub / (sP.Ksn + sub) - ...
bd * sP.muMaxPred1 * phageResPrey / ...
((sP.Knp + senPrey + bdPersistPrey + phageResPrey) * sP.yieldBPerN) + ...
senPrey * sP.rateDevResistance;
dPred1 = sP.kP * bdplast - ...
bd * sP.muMaxPred1 * (senPrey + phageResPrey) / ...
((sP.Knp + (senPrey + bdPersistPrey + phageResPrey)) * ...
sP.yieldBPerP) - ...
sP.mortality * bd;
dBdelloplast1 = bd * sP.muMaxPred1 * (senPrey + phageResPrey) / ...
(sP.Knp + senPrey + bdPersistPrey + phageResPrey) - ...
sP.kP * bdplast / sP.yieldPPerB;
dPred2 = sP.kV * infCell - ...
phage * sP.muMaxPred2 * (senPrey + bdPersistPrey) / sP.yieldIPerV;
dBdelloplast2 = phage * sP.muMaxPred2 * (senPrey + bdPersistPrey) - ...
sP.kV * infCell / sP.yieldVPerI;
dDeadPrey = sP.kP * bdplast / sP.yieldPPerB + ...
sP.kV * infCell / sP.yieldVPerI;
%write results
dy = [dSub; dSenPrey; dBdPersistPrey; dPhageResPrey; dPred1; ...
dBdelloplast1; dPred2; dBdelloplast2; dDeadPrey];
end