function dndt = mc_simple(t,n,param) % CONSTITUTIVE EXPRESSION SIMPLE MODEL: % Establishes the equations of the designed constitutive model. It % gives us the variation of the concentrations of the biochemical % species of the reactions in each instant of time. %Abbreviations: %PoI: protein of interest %mRNA: messenger RNA % param = [Cr,p,dm,dp,mu,Kmax] Vector of parameters % Cr: constitutive transcription rate(molec.min-1) % p: translation rate (min-1) % dm: mRNA degradation rate(min-1) % dp: PoI degradation rate(min-1) % mu: dilution rate (min-1) % Kmax: maximum growth capacity (cells) Cr = param(1); p = param(2); dm = param(3); dp = param(4); mu = param(5); Kmax = param(6); % Differential equations: %n(1): mRNA concentration (molec) %n(2): PoI concentration (molec) %n(3): number of cells dn1 = Cr - (dm+mu)*n(1); dn2 = p*n(1) - (dp+mu)*n(2); dn3 = mu*n(3)*(1 - n(3)/Kmax); dndt = [dn1; dn2; dn3]; return