Mathematical model to simulate SARS-CoV-2 infection within-host

Lead Research Organisation: University of Bath
Department Name: Mathematical Sciences

Abstract

Mathematical models are vital in advising strategy when dealing with pandemics, from helping to develop individual treatment strategies to guiding the national public health approach. Current modelling efforts concentrate on transmission and do not focus on variation within the host. There is clear evidence that some subgroups of patients are likely to have more severe disease and poorer outcomes. The reasons for these associations are not clear.

We have developed a within-host mathematical and computational model of SARS-CoV-2 infection that is capable of simulating viral spread in lung cells. Preliminary results illustrate how our model is able to study, in isolation, particular immune dysfunctions associated with severe COVID-19. This is difficult to achieve with biological experiments. Our results have suggested that impairing the function of Natural Killer (NK) cells, important for combatting viral infections, skews the immune response in ways that cause severe disease. Additionally, our model shows that manipulating the levels of defence molecules that immune and infected cells produce to try and fight the infection can lead to severe viral infection, similar to that observed in severe COVID-19 patients.

We have laid important groundwork for future code development in this project; parameterisation and validation, and application. In terms of application, we intend to investigate the influence of initial (and continual) viral load deposition (amount of virus that initiates infection) on the spread of infection. Additionally, we aim to consider more in-depth models of the production of defence molecules, known as cytokines (in particular a cytokine known as type I interferon). We will investigate their dysfunction in their regulatory pathway and their impact on the spread of infection. The model will enhance our understanding of COVID-19 pathophysiology. In this project we will integrate mathematical models that simulate drug distribution in the body. This will allow us to test alternative treatment strategies, such as various drug scheduling and dosing intervals, and refine therapy for specific subsets of patients. Many people remain vulnerable to this infection, but greater knowledge of how to deliver successful treatment strategies will provide hope for those who become critically unwell. It will also diminish the suffering of those who experience non-critical, but nonetheless unpleasant, disease. Understanding gained from our model simulations may also lead to improved management of the long-term effects of COVID-19 (long COVID).

The nature of the model will allow us to investigate why different subgroups are at greater risk, and why they are perhaps most likely to become infected. This can inform public health strategy to protect the most vulnerable members of society. On completion of this project, there is scope to link our within-host model with population-level and environmental models. This could help us to understand more about the course of infectiousness in individuals, aiding guidance around self-isolation and ultimately helping to reduce transmission. It will also help to understand how and why there is heterogeneity in different subsets of patients' transmissibility.

Using our mathematical framework, we will also create a mathematical tool that will allow other infectious disease researchers to model the within-host dynamics of newly emerging pathogens in the future.
 
Description In our mathematical model, we found that an increase in infection was observed with a higher number of infectious, apoptotic and removed epithelial cells being seen for the highest Multiplicity of Infection value considered. The increase in infection correlated with an increase in interferon, cytokine and chemokine levels (a similar correlation has been observed from longitudinal studies).

Increasing the IFN-I secretion delay from epithelial cells increases the spread of infection, with higher infectious, apoptotic and removed epithelial cells being observed for the longest IFN-I secretion delay. For the longest IFN-I secretion delay, reduced levels of IFN-I were observed (as expected) accompanied by higher levels of IL-6. We observed a clear decrease in the spread of infection as the IFN-I secretion value increased.

For the largest macrophage activation half-max value we observe (a) a significant influx of resting macrophages, (b) reduced infection, (c) lower IL-6 levels. We see a slight increase in generic mononuclear phagocyte chemokine, which could (in combination with IL-6) be enhancing the recruitment of resting macrophages. This suggests a possible feedback between apoptotic cells (epithelial and macrophage) and the recruitment of macrophages, driving tissue injury and inflammation. This requires further investigation.
Exploitation Route When responding to a viral infection, no component of the immune system is expected to work by itself. However, mathematical and computational studies, as well as experimental studies, necessarily require reductions and simplifications, as well as iterative development. Indeed, in this work, we chose to focus only a single type and phenotype of macrophages, as well as a limited set of cytokines and chemokines. Therefore, extending the model to include multiple immune cell populations, as well as subpopulations, and a larger set of cytokines and chemokines is a subject of future work. Furthermore, we chose to focus on a single (classical) activation route, whereas in reality macrophages can alternatively activate; therefore, extending the model to allow for multiple activation routes is also a subject of future work. Multiscale, individual-based models are inherently complex and contain many parameter values which, in general, cannot be expected to directly fit to other modelling work, let alone experimental data. Therefore, improving the parameter estimation and fitting to experimental data, either directly or indirectly via ODE models is a subject of future work. The presented model has the capability to include more complex intracellular pathways. Therefore, in the future, we would like to study, using more detailed ODE models, the crosstalk between virus replication, IFN-I induction and signalling, as well as pro-inflammatory cytokine induction and signalling, and how the crosstalk impacts the spread of infection. We would also like to consider additional models for virus entry and replication that are specific to different tissue environments. Finally, we would like to take the model to larger scales, either through computational improvements (such as massive parallelisation) or through mathematical models which can bridge-the-gap between the scales.
Sectors Healthcare

 
Description Collaboration with Dr David Hughes, University of St Andrews 
Organisation University of St Andrews
Department School of Biology
Country United Kingdom 
Sector Academic/University 
PI Contribution Results from our model simulations have inspired new work in David's lab.
Collaborator Contribution David is a virologist and he has helped us to develop our mathematical model of SARS-CoV-2.
Impact Article submitted to Royal Society Open Science. Preprint: https://www.biorxiv.org/content/10.1101/2022.05.06.490883v1
Start Year 2020
 
Description Collaboration with Ingo Johannessen, NHS Lothian 
Organisation Royal Infirmary of Edinburgh
Country United Kingdom 
Sector Hospitals 
PI Contribution Results from our model have brought new insights.
Collaborator Contribution Ingo is a clinical virologist and he has helped us to develop the SARS-CoV-2 mathematical model.
Impact Journal article submitted to Royal Society Open Science. Preprint: https://www.biorxiv.org/content/10.1101/2022.05.06.490883v1
Start Year 2020
 
Description Collaboration with Professor David Dockrell 
Organisation University of Edinburgh
Department MRC Centre for Inflammation Research
Country United Kingdom 
Sector Academic/University 
PI Contribution Results from our model have inspired work in David's lab.
Collaborator Contribution David is an immunologist who has helped us to develop the SARS-CoV-2 mathematical model.
Impact Paper submitted to Royal Society Open Science. Preprint: https://www.biorxiv.org/content/10.1101/2022.05.06.490883v1
Start Year 2020
 
Description RAMP Innovation Outreach Award 
Form Of Engagement Activity Engagement focused website, blog or social media channel
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Public/other audiences
Results and Impact I was awarded a RAMP (Rapid Assistance in Modelling the Pandemic) Innovation Outreach Award. For this, I created an animation: https://www.youtube.com/watch?v=eA-dKLmnosY
Year(s) Of Engagement Activity 2022