Getting the most out of maths: How to coordinate mathematical modelling research to support a pandemic, lessons learnt from three initiatives that were part of the COVID-19 response in the UK.

In March 2020 mathematics became a key part of the scientific advice to the UK government on the pandemic response to COVID-19. Mathematical and statistical modelling provided critical information on the spread of the virus and the potential impact of different interventions. The unprecedented scale of the challenge led the epidemiological modelling community in the UK to be pushed to its limits. At the same time, mathematical modellers across the country were keen to use their knowledge and skills to support the COVID-19 modelling effort. However, this sudden great interest in epidemiological modelling needed to be coordinated to provide much-needed support, and to limit the burden on epidemiological modellers already very stretched for time. In this paper we describe three initiatives set up in the UK in spring 2020 to coordinate the mathematical sciences research community in supporting mathematical modelling of COVID-19. Each initiative had different primary aims and worked to maximise synergies between the various projects. We reflect on the lessons learnt, highlighting the key roles of pre-existing research collaborations and focal centres of coordination in contributing to the success of these initiatives. We conclude with recommendations about important ways in which the scientific research community could be better prepared for future pandemics. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Challenges of integrating economics into epidemiological analysis of and policy responses to emerging infectious diseases.

COVID-19 has shown that the consequences of a pandemic are wider-reaching than cases and deaths. Morbidity and mortality are important direct costs, but infectious diseases generate other direct and indirect benefits and costs as the economy responds to these shocks: some people lose, others gain and people modify their behaviours in ways that redistribute these benefits and costs. These additional effects feedback on health outcomes to create a complicated interdependent system of health and non-health outcomes. As a result, interventions primarily intended to reduce the burden of disease can have wider societal and economic effects and more complicated and unintended, but possibly not anticipable, system-level influences on the epidemiological dynamics themselves. Capturing these effects requires a systems approach that encompasses more direct health outcomes. Towards this end, in this article we discuss the importance of integrating epidemiology and economic models, setting out the key challenges which such a merging of epidemiology and economics presents. We conclude that understanding people's behaviour in the context of interventions is key to developing a more complete and integrated economic-epidemiological approach; and a wider perspective on the benefits and costs of interventions (and who these fall upon) will help society better understand how to respond to future pandemics.

Commentary on the use of the reproduction number R during the COVID-19 pandemic.

Since the beginning of the COVID-19 pandemic, the reproduction number [Formula: see text] has become a popular epidemiological metric used to communicate the state of the epidemic. At its most basic, [Formula: see text] is defined as the average number of secondary infections caused by one primary infected individual. [Formula: see text] seems convenient, because the epidemic is expanding if [Formula: see text] and contracting if [Formula: see text]. The magnitude of [Formula: see text] indicates by how much transmission needs to be reduced to control the epidemic. Using [Formula: see text] in a naïve way can cause new problems. The reasons for this are threefold: (1) There is not just one definition of [Formula: see text] but many, and the precise definition of [Formula: see text] affects both its estimated value and how it should be interpreted. (2) Even with a particular clearly defined [Formula: see text], there may be different statistical methods used to estimate its value, and the choice of method will affect the estimate. (3) The availability and type of data used to estimate [Formula: see text] vary, and it is not always clear what data should be included in the estimation. In this review, we discuss when [Formula: see text] is useful, when it may be of use but needs to be interpreted with care, and when it may be an inappropriate indicator of the progress of the epidemic. We also argue that careful definition of [Formula: see text], and the data and methods used to estimate it, can make [Formula: see text] a more useful metric for future management of the epidemic.

Resource Allocation for Epidemic Control Across Multiple Sub-populations.

The number of pathogenic threats to plant, animal and human health is increasing. Controlling the spread of such threats is costly and often resources are limited. A key challenge facing decision makers is how to allocate resources to control the different threats in order to achieve the least amount of damage from the collective impact. In this paper we consider the allocation of limited resources across n independent target populations to treat pathogens whose spread is modelled using the susceptible-infected-susceptible model. Using mathematical analysis of the systems dynamics, we show that for effective disease control, with a limited budget, treatment should be focused on a subset of populations, rather than attempting to treat all populations less intensively. The choice of populations to treat can be approximated by a knapsack-type problem. We show that the knapsack closely approximates the exact optimum and greatly outperforms a number of simpler strategies. A key advantage of the knapsack approximation is that it provides insight into the way in which the economic and epidemiological dynamics affect the optimal allocation of resources. In particular using the knapsack approximation to apportion control takes into account two important aspects of the dynamics: the indirect interaction between the populations due to the shared pool of limited resources and the dependence on the initial conditions.

Filling in the gaps: estimating numbers of chlamydia tests and diagnoses by age group and sex before and during the implementation of the English National Screening Programme, 2000 to 2012.

To inform mathematical modelling of the impact of chlamydia screening in England since 2000, a complete picture of chlamydia testing is needed. Monitoring and surveillance systems evolved between 2000 and 2012. Since 2012, data on publicly funded chlamydia tests and diagnoses have been collected nationally. However, gaps exist for earlier years. We collated available data on chlamydia testing and diagnosis rates among 15-44-year-olds by sex and age group for 2000-2012. Where data were unavailable, we applied data- and evidence-based assumptions to construct plausible minimum and maximum estimates and set bounds on uncertainty. There was a large range between estimates in years when datasets were less comprehensive (2000-2008); smaller ranges were seen hereafter. In 15-19-year-old women in 2000, the estimated diagnosis rate ranged between 891 and 2,489 diagnoses per 100,000 persons. Testing and diagnosis rates increased between 2000 and 2012 in women and men across all age groups using minimum or maximum estimates, with greatest increases seen among 15-24-year-olds. Our dataset can be used to parameterise and validate mathematical models and serve as a reference dataset to which trends in chlamydia-related complications can be compared. Our analysis highlights the complexities of combining monitoring and surveillance datasets.

Modeling ion channel dynamics through reflected stochastic differential equations.

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks.

Comparison of continuous and discrete stochastic ion channel models.

The stochastic behaviour of ion channels can be described by a discrete model or by an approximate continuous approach. While the discrete approach is exact, it is also less computationally efficient, and so the continuous model is often the method of choice since it allows for incorporation into a multiscale environment. However, in recent years the accuracy of the stochastic continuous approach for calculating statistics of certain quantities in the Hodgkin-Huxley model has come into question. In this paper, we show that by correct formulation of the continuous model, the first two moments in the number of open sodium and potassium channels in the Hodgkin-Huxley model, calculated under voltage clamp conditions using the continuous approach are in good agreement with those obtained from the discrete model.