A Novel Stochastic Multi-Scale Model of Francisella tularensis Infection to Predict Risk of Infection in a Laboratory

Quantifying in vitro B. anthracis growth and PA production and decay: a mathematical modelling approach

Protective antigen (PA) is a protein produced by Bacillus anthracis. It forms part of the anthrax toxin and is a key immunogen in US and UK anthrax vaccines. In this study, we have conducted experiments to quantify PA in the supernatants of cultures of B. anthracis Sterne strain, which is the strain used in the manufacture of the UK anthrax vaccine. Then, for the first time, we quantify PA production and degradation via mathematical modelling and Bayesian statistical techniques, making use of this new experimental data as well as two other independent published data sets. We propose a single mathematical model, in terms of delay differential equations (DDEs), which can explain the in vitro dynamics of all three data sets. Since we did not heat activate the B. anthracis spores prior to inoculation, germination occurred much slower in our experiments, allowing us to calibrate two additional parameters with respect to the other data sets. Our model is able to distinguish between natural PA decay and that triggered by bacteria via proteases. There is promising consistency between the different independent data sets for most of the parameter estimates. The quantitative characterisation of B. anthracis PA production and degradation obtained here will contribute towards the ambition to include a realistic description of toxin dynamics, the host immune response, and anti-toxin treatments in future mechanistic models of anthrax infection.

Using mixture density networks to emulate a stochastic within-host model of Francisella tularensis infection

For stochastic models with large numbers of states, analytical techniques are often impractical, and simulations time-consuming and computationally demanding. This limitation can hinder the practical implementation of such models. In this study, we demonstrate how neural networks can be used to develop emulators for two outputs of a stochastic within-host model of Francisella tularensis infection: the dose-dependent probability of illness and the incubation period. Once the emulators are constructed, we employ Markov Chain Monte Carlo sampling methods to parameterize the within-host model using records of human infection. This inference is only possible through the use of a mixture density network to emulate the incubation period, providing accurate approximations of the corresponding probability distribution. Notably, these estimates improve upon previous approaches that relied on bacterial counts from the lungs of macaques. Our findings reveal a 50% infectious dose of approximately 10 colony-forming units and we estimate that the incubation period can last for up to 11 days following low dose exposure.

A Stochastic Intracellular Model of Anthrax Infection With Spore Germination Heterogeneity

We present a stochastic mathematical model of the intracellular infection dynamics of Bacillus anthracis in macrophages. Following inhalation of B. anthracis spores, these are ingested by alveolar phagocytes. Ingested spores then begin to germinate and divide intracellularly. This can lead to the eventual death of the host cell and the extracellular release of bacterial progeny. Some macrophages successfully eliminate the intracellular bacteria and will recover. Here, a stochastic birth-and-death process with catastrophe is proposed, which includes the mechanism of spore germination and maturation of B. anthracis. The resulting model is used to explore the potential for heterogeneity in the spore germination rate, with the consideration of two extreme cases for the rate distribution: continuous Gaussian and discrete Bernoulli. We make use of approximate Bayesian computation to calibrate our model using experimental measurements from in vitro infection of murine peritoneal macrophages with spores of the Sterne 34F2 strain of B. anthracis. The calibrated stochastic model allows us to compute the probability of rupture, mean time to rupture, and rupture size distribution, of a macrophage that has been infected with one spore. We also obtain the mean spore and bacterial loads over time for a population of cells, each assumed to be initially infected with a single spore. Our results support the existence of significant heterogeneity in the germination rate, with a subset of spores expected to germinate much later than the majority. Furthermore, in agreement with experimental evidence, our results suggest that most of the spores taken up by macrophages are likely to be eliminated by the host cell, but a few germinated spores may survive phagocytosis and lead to the death of the infected cell. Finally, we discuss how this stochastic modelling approach, together with dose-response data, allows us to quantify and predict individual infection risk following exposure.

Stochastic dynamics of Francisella tularensis infection and replication

We study the pathogenesis of Francisella tularensis infection with an experimental mouse model, agent-based computation and mathematical analysis. Following inhalational exposure to Francisella tularensis SCHU S4, a small initial number of bacteria enter lung host cells and proliferate inside them, eventually destroying the host cell and releasing numerous copies that infect other cells. Our analysis of disease progression is based on a stochastic model of a population of infectious agents inside one host cell, extending the birth-and-death process by the occurrence of catastrophes: cell rupture events that affect all bacteria in a cell simultaneously. Closed expressions are obtained for the survival function of an infected cell, the number of bacteria released as a function of time after infection, and the total bacterial load. We compare our mathematical analysis with the results of agent-based computation and, making use of approximate Bayesian statistical inference, with experimental measurements carried out after murine aerosol infection with the virulent SCHU S4 strain of the bacterium Francisella tularensis, that infects alveolar macrophages. The posterior distribution of the rate of replication of intracellular bacteria is consistent with the estimate that the time between rounds of bacterial division is less than 6 hours in vivo.

Infecting a host cell is required for the replication of many types of bacteria and viruses. In some cases, infected cells release new infectious agents continuously over their lifetime. In others, such as the Francisella tularensis bacterium studied here, they are released in a single burst that coincides with the cell’s death. We show how a stochastic model, the birth-and-death process with catastrophe, can be used to characterise infection in a single cell, thereby allowing us to account for burst events and quantify the kinetics of pathogenesis in the lung, the initial site of infection, as well as in other organs that the infection spreads to. We learn about the parameters of the mathematical model of Francisella tularensis infection making use of the experimental measurements of bacterial loads, together with approximate Bayesian statistical inference methods. The most important parameter describing the pathogenesis is the rate of replication of intracellular bacteria.

A Multicompartment SIS Stochastic Model with Zonal Ventilation for the Spread of Nosocomial Infections: Detection, Outbreak Management, and Infection Control

In this work, we study the environmental and operational factors that influence airborne transmission of nosocomial infections. We link a deterministic zonal ventilation model for the airborne distribution of infectious material in a hospital ward, with a Markovian multicompartment SIS model for the infection of individuals within this ward, in order to conduct a parametric study on ventilation rates and their effect on the epidemic dynamics. Our stochastic model includes arrival and discharge of patients, as well as the detection of the outbreak by screening events or due to symptoms being shown by infective patients. For each ventilation setting, we measure the infectious potential of a nosocomial outbreak in the hospital ward by means of a summary statistic: the number of infections occurred within the hospital ward until end or declaration of the outbreak. We analytically compute the distribution of this summary statistic, and carry out local and global sensitivity analysis in order to identify the particular characteristics of each ventilation regime with the largest impact on the epidemic spread. Our results show that ward ventilation can have a significant impact on the infection spread, especially under slow detection scenarios or in overoccupied wards, and that decreasing the infection risk for the whole hospital ward might increase the risk in specific areas of the health-care facility. Moreover, the location of the initial infective individual and the protocol in place for outbreak declaration both form an interplay with ventilation of the ward.

A Multicompartment SIS Stochastic Model with Zonal Ventilation for the Spread of Nosocomial Infections: Detection, Outbreak Management, and Infection Control

In this work, we study the environmental and operational factors that influence airborne transmission of nosocomial infections. We link a deterministic zonal ventilation model for the airborne distribution of infectious material in a hospital ward, with a Markovian multicompartment SIS model for the infection of individuals within this ward, in order to conduct a parametric study on ventilation rates and their effect on the epidemic dynamics. Our stochastic model includes arrival and discharge of patients, as well as the detection of the outbreak by screening events or due to symptoms being shown by infective patients. For each ventilation setting, we measure the infectious potential of a nosocomial outbreak in the hospital ward by means of a summary statistic: the number of infections occurred within the hospital ward until end or declaration of the outbreak. We analytically compute the distribution of this summary statistic, and carry out local and global sensitivity analysis in order to identify the particular characteristics of each ventilation regime with the largest impact on the epidemic spread. Our results show that ward ventilation can have a significant impact on the infection spread, especially under slow detection scenarios or in overoccupied wards, and that decreasing the infection risk for the whole hospital ward might increase the risk in specific areas of the health-care facility. Moreover, the location of the initial infective individual and the protocol in place for outbreak declaration both form an interplay with ventilation of the ward.