MadonnaModel {Conjugative Plasmid Model} {Transfer of the plasmid converts N to C} METHOD EULER STARTTIME = 0 STOPTIME= 2000 DT = 0.005 DTOUT = 10 init N= 0 {Non-CRISPR} init C =1e8 {CRISPR} init NP=1e2 {Non Crisper Plasmid } init CP =0 {CRISPER PLASMID} init X =0 {Immune CRISPR} init R=0.065 {Resource} vn=1 {Max growth sensitive N} vnp= .95 {Max growth plasmind bearing N} vc = .9 {Max growth C} vcp=.9 {Maximum growrth CP} vx=0.9 {Maximum growth X} A=50 {Input resource} w=0.2 {flow rate} e=5E-7 {Conversion efficiency} k=0.25 {Monod constant} gnn= 1E-9 {Rate constant of plasmid transfer Np to N} gnc= 1E-9{Rate constant of plasmid transfer Np to C} gcn = 1E-9 {Rate constant of plasmid transfer Cp to N} gcc = 1e-9 {Rate constant of plasmid transfer Cp to C} tn =1e-4 {Segregation rate NP} tc = 1e=4 {Segregation rate CP} m = 1e-3 {CRISPR conversion N--> C} v =1e-8 {loss of CRISPR > X, C-->N and CP --> NP} d/dt (R) = w*(A-R) - (vn*N +vnp*NP +vc*C +vcp*CP +vx*X)* (R/(R+k))*e d/dt (N) = (vn*R/(R+k))*N - gnn*N*NP - gcn*N*CP +tn*NP +v*C +v*X- w*N d/dt (NP) = (vnp*R/(R+k))*NP+gnn*N*NP + gcn*N*CP -tn*NP +v*CP- w*NP d/dt (C) = (vc*R/(R+k))*C - gnc*NP*C*(1-m) - gcc*CP*C*(1-m) +tc*CP -v*C - w*C d/dt (CP) = (vcp*R/(R+k))*CP +gnc*NP*C +gcc*CP*C - tc*CP -v*CP - w*CP d/dt (X) =(vx*R/(R+k))*X + gnc*NP*C*m + gcc*CP*C*m- v*X - w*X Geneva5 !AMU e ]Q 9 a@ eulerMbP?MbP?MbP? $bf^ZH & ,6j.q  Q= y$& (*,e. !