ObservableHQ notebook: Identification of phage load and injection time in simulated phage therapy occurrences

Dynamic code for case 1 reported in the paper. The scenario refers to the data reported by Payne, R.J. and Jansen, V.A. "Understanding Bacteriophage Therapy as a Density-Dependent Kinetic Process", J. Theor. Biol. 2001, 208, 37–48.

The default values for the parameters below as well as the ODEs are set according to the paper. The output of this computation is the plot of the densities of the hypothetical bacterial and phagial populations. The plot has been extended five-fold with respect to the work of Payne and Jansen (100 hours versus 20 hours).

The user can have a feeling on how phage therapy depends on case-specific conditions by modifying the parameters using the slider below. The list of parameters is as follows:

• tᵢₙ = time of administration of the phage dose (h)
• vᵢₙ = amount of the phage dose (PFU/mL)
• S₀ = initial density of susceptible bacteria (CFU/mL)
• K = carrying capacity of the environment (CFU/mL)
• μ = maximum growth rate of the bacteria (1/h)
• ω = outflow of the environment (1/mL)
• δ = infection rate of the phage (1/min)
• η = inverse of latency time (min)
• λ = degradation rate of the phage (1/day)
• β = burst size of the phage (PFU)