Mathematical model of interaction Escherichia coli and Coliphages

We propose a mathematical model based in ordinary differential equations between bacterial pathogen and Bacteriophages to describe the infection dynamics of these populations, for which we use a nonlinear function with an inhibitory effect. We study the stability of the model using the Lyapunov theory and the second additive compound matrix and perform a global sensitivity analysis to elucidate the most influential parameters in the model, besides we make a parameter estimation using growth data of Escherichia coli (E.coli) bacteria in presence of Coliphages (bacteriophages that infect E.coli) with different multiplicity of infection. We found a threshold that indicates whether the bacteriophage concentration will coexist with the bacterium (the coexistence equilibrium ) or become extinct (phages extinction equilibrium), the first equilibrium is locally asymptotically stable while the other is globally asymptotically stable depending on the magnitude of this threshold. Beside we found that the dynamics of the model is particularly affected by infection rate of bacteria and Half-saturation phages density. Parameter estimation show that all multiplicities of infection are effective in eliminating infected bacteria but the smaller one leaves a higher number of bacteriophages at the end of this elimination.