1
|
Williams B, Carruthers J, Gillard JJ, Lythe G, Perelson AS, Ribeiro RM, Molina-París C, López-García M. The reproduction number and its probability distribution for stochastic viral dynamics. J R Soc Interface 2024; 21:20230400. [PMID: 38264928 PMCID: PMC10806437 DOI: 10.1098/rsif.2023.0400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Accepted: 12/18/2023] [Indexed: 01/25/2024] Open
Abstract
We consider stochastic models of individual infected cells. The reproduction number, R, is understood as a random variable representing the number of new cells infected by one initial infected cell in an otherwise susceptible (target cell) population. Variability in R results partly from heterogeneity in the viral burst size (the number of viral progeny generated from an infected cell during its lifetime), which depends on the distribution of cellular lifetimes and on the mechanism of virion release. We analyse viral dynamics models with an eclipse phase: the period of time after a cell is infected but before it is capable of releasing virions. The duration of the eclipse, or the subsequent infectious, phase is non-exponential, but composed of stages. We derive the probability distribution of the reproduction number for these viral dynamics models, and show it is a negative binomial distribution in the case of constant viral release from infectious cells, and under the assumption of an excess of target cells. In a deterministic model, the ultimate in-host establishment or extinction of the viral infection depends entirely on whether the mean reproduction number is greater than, or less than, one, respectively. Here, the probability of extinction is determined by the probability distribution of R, not simply its mean value. In particular, we show that in some cases the probability of infection is not an increasing function of the mean reproduction number.
Collapse
Affiliation(s)
- Bevelynn Williams
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds, UK
| | | | - Joseph J. Gillard
- CBR Division, Defence Science and Technology Laboratory, Salisbury, UK
| | - Grant Lythe
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds, UK
| | - Alan S. Perelson
- T-6, Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Ruy M. Ribeiro
- T-6, Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Carmen Molina-París
- T-6, Theoretical Biology and Biophysics, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Martín López-García
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds, UK
| |
Collapse
|
2
|
Dobrovolny HM. Mathematical Modeling of Virus-Mediated Syncytia Formation: Past Successes and Future Directions. Results Probl Cell Differ 2024; 71:345-370. [PMID: 37996686 DOI: 10.1007/978-3-031-37936-9_17] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2023]
Abstract
Many viruses have the ability to cause cells to fuse into large multi-nucleated cells, known as syncytia. While the existence of syncytia has long been known and its importance in helping spread viral infection within a host has been understood, few mathematical models have incorporated syncytia formation or examined its role in viral dynamics. This review examines mathematical models that have incorporated virus-mediated cell fusion and the insights they have provided on how syncytia can change the time course of an infection. While the modeling efforts are limited, they show promise in helping us understand the consequences of syncytia formation if future modeling efforts can be coupled with appropriate experimental efforts to help validate the models.
Collapse
Affiliation(s)
- Hana M Dobrovolny
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, USA.
| |
Collapse
|
3
|
Peña-Gil N, Randazzo W, Carmona-Vicente N, Santiso-Bellón C, Cárcamo-Cálvo R, Navarro-Lleó N, Monedero V, Yebra MJ, Buesa J, Gozalbo-Rovira R, Rodríguez-Díaz J. Culture of Human Rotaviruses in Relevant Models Shows Differences in Culture-Adapted and Nonculture-Adapted Strains. Int J Mol Sci 2023; 24:17362. [PMID: 38139191 PMCID: PMC10743750 DOI: 10.3390/ijms242417362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2023] [Revised: 12/04/2023] [Accepted: 12/07/2023] [Indexed: 12/24/2023] Open
Abstract
Rotavirus (RV) is the leading cause of acute gastroenteritis (AGE) in children under 5 years old worldwide, and several studies have demonstrated that histo-blood group antigens (HBGAs) play a role in its infection process. In the present study, human stool filtrates from patients diagnosed with RV diarrhea (genotyped as P[8]) were used to infect differentiated Caco-2 cells (dCaco-2) to determine whether such viral strains of clinical origin had the ability to replicate in cell cultures displaying HBGAs. The cell culture-adapted human RV Wa model strain (P[8] genotype) was used as a control. A time-course analysis of infection was conducted in dCaco-2 at 1, 24, 48, 72, and 96 h. The replication of two selected clinical isolates and Wa was further assayed in MA104, undifferentiated Caco-2 (uCaco-2), HT29, and HT29-M6 cells, as well as in monolayers of differentiated human intestinal enteroids (HIEs). The results showed that the culture-adapted Wa strain replicated more efficiently in MA104 cells than other utilized cell types. In contrast, clinical virus isolates replicated more efficiently in dCaco-2 cells and HIEs. Furthermore, through surface plasmon resonance analysis of the interaction between the RV spike protein (VP8*) and its glycan receptor (the H antigen), the V7 RV clinical isolate showed 45 times better affinity compared to VP8* from the Wa strain. These findings support the hypothesis that the differences in virus tropism between clinical virus isolates and RV Wa could be a consequence of the different HBGA contents on the surface of the cell lines employed. dCaco-2, HT29, and HT29M6 cells and HIEs display HBGAs on their surfaces, whereas MA104 and uCaco-2 cells do not. These results indicate the relevance of using non-cell culture-adapted human RV to investigate the replication of rotavirus in relevant infection models.
Collapse
Affiliation(s)
- Nazaret Peña-Gil
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
- Instituto de Investigación INCLIVA, Hospital Clínico Universitario de Valencia, 46010 Valencia, Spain
| | - Walter Randazzo
- Department of Preservation and Food Safety Technologies, IATA-CSIC, Av. Agustín Escardino 7, 46980 Paterna, Spain;
| | - Noelia Carmona-Vicente
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
| | - Cristina Santiso-Bellón
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
- Instituto de Investigación INCLIVA, Hospital Clínico Universitario de Valencia, 46010 Valencia, Spain
| | - Roberto Cárcamo-Cálvo
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
- Instituto de Investigación INCLIVA, Hospital Clínico Universitario de Valencia, 46010 Valencia, Spain
| | - Noemi Navarro-Lleó
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
| | - Vicente Monedero
- Department of Biotechnology, IATA-CSIC, Av. Agustín Escardino 7, 46980 Paterna, Spain; (V.M.); (M.J.Y.)
| | - María J. Yebra
- Department of Biotechnology, IATA-CSIC, Av. Agustín Escardino 7, 46980 Paterna, Spain; (V.M.); (M.J.Y.)
| | - Javier Buesa
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
- Instituto de Investigación INCLIVA, Hospital Clínico Universitario de Valencia, 46010 Valencia, Spain
| | - Roberto Gozalbo-Rovira
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
- Instituto de Investigación INCLIVA, Hospital Clínico Universitario de Valencia, 46010 Valencia, Spain
| | - Jesús Rodríguez-Díaz
- Department of Microbiology, School of Medicine, University of Valencia, Av. Blasco Ibáñez 15, 46010 Valencia, Spain; (N.P.-G.); (N.C.-V.); (C.S.-B.); (R.C.-C.); (N.N.-L.); (J.B.)
- Instituto de Investigación INCLIVA, Hospital Clínico Universitario de Valencia, 46010 Valencia, Spain
| |
Collapse
|
4
|
Quirouette C, Cresta D, Li J, Wilkie KP, Liang H, Beauchemin CAA. The effect of random virus failure following cell entry on infection outcome and the success of antiviral therapy. Sci Rep 2023; 13:17243. [PMID: 37821517 PMCID: PMC10567758 DOI: 10.1038/s41598-023-44180-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 10/04/2023] [Indexed: 10/13/2023] Open
Abstract
A virus infection can be initiated with very few or even a single infectious virion, and as such can become extinct, i.e. stochastically fail to take hold or spread significantly. There are many ways that a fully competent infectious virion, having successfully entered a cell, can fail to cause a productive infection, i.e. one that yields infectious virus progeny. Though many stochastic models (SMs) have been developed and used to estimate a virus infection's establishment probability, these typically neglect infection failure post virus entry. The SM presented herein introduces parameter [Formula: see text] which corresponds to the probability that a virion's entry into a cell will result in a productive cell infection. We derive an expression for the likelihood of infection establishment in this new SM, and find that prophylactic therapy with an antiviral reducing [Formula: see text] is at least as good or better at decreasing the establishment probability, compared to antivirals reducing the rates of virus production or virus entry into cells, irrespective of the SM parameters. We investigate the difference in the fraction of cells consumed by so-called extinct versus established virus infections, and find that this distinction becomes biologically meaningless as the probability of establishment approaches zero. We explain why the release of virions continuously over an infectious cell's lifespan, rather than as a single burst at the end of the cell's lifespan, does not result in an increased risk of infection extinction. We show, instead, that the number of virus released, not the timing of the release, affects infection establishment and associated critical antiviral efficacy.
Collapse
Affiliation(s)
| | - Daniel Cresta
- Department of Physics, Toronto Metropolitan University, Toronto, Canada
| | - Jizhou Li
- Interdisciplinary Theoretical and Mathematical Sciences (iTHEMS), RIKEN, Wako, Japan
| | - Kathleen P Wilkie
- Department of Mathematics, Toronto Metropolitan University, Toronto, Canada
| | - Haozhao Liang
- Nishina Center for Accelerator-Based Science (RNC), RIKEN, Wako, Japan
- Department of Physics, University of Tokyo, Tokyo, Japan
| | - Catherine A A Beauchemin
- Department of Physics, Toronto Metropolitan University, Toronto, Canada.
- Interdisciplinary Theoretical and Mathematical Sciences (iTHEMS), RIKEN, Wako, Japan.
| |
Collapse
|
5
|
Amidei A, Dobrovolny HM. Estimation of virus-mediated cell fusion rate of SARS-CoV-2. Virology 2022; 575:91-100. [PMID: 36088794 PMCID: PMC9449781 DOI: 10.1016/j.virol.2022.08.016] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 08/27/2022] [Accepted: 08/28/2022] [Indexed: 12/22/2022]
Abstract
Several viruses have the ability to form large multinucleated cells known as syncytia. Many properties of syncytia and the role they play in the evolution of a viral infection are not well understood. One basic question that has not yet been answered is how quickly syncytia form. We use a novel mathematical model of cell-cell fusion assays and apply it to experimental data from SARS-CoV-2 fusion assays to provide the first estimates of virus-mediated cell fusion rate. We find that for SARS-CoV2, the fusion rate is in the range of 6 × 10−4–12×10−4/h. We also use our model to compare fusion rates when the protease TMPRSS2 is overexpressed (2–4 times larger fusion rate), when the protease furin is removed (one third the original fusion rate), and when the spike protein is altered (1/10th the original fusion rate). The use of mathematical models allows us to provide additional quantitative information about syncytia formation.
Collapse
Affiliation(s)
- Ava Amidei
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, TX, USA
| | - Hana M Dobrovolny
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, USA.
| |
Collapse
|
6
|
Sultana S, González-Parra G, Arenas AJ. Dynamics of toxoplasmosis in the cat's population with an exposed stage and a time delay. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:12655-12676. [PMID: 36654016 DOI: 10.3934/mbe.2022591] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
We propose a new mathematical model to investigate the effect of the introduction of an exposed stage for the cats who become infected with the T. gondii parasite, but that are not still able to produce oocysts in the environment. The model considers a time delay in order to represent the duration of the exposed stage. Besides the cat population the model also includes the oocysts related to the T. gondii in the environment. The model includes the cats since they are the only definitive host and the oocysts, since they are relevant to the dynamics of toxoplasmosis. The model considers lifelong immunity for the recovered cats and vaccinated cats. In addition, the model considers that cats can get infected through an effective contact with the oocysts in the environment. We find conditions such that the toxoplasmosis disease becomes extinct. We analyze the consequences of considering the exposed stage and the time delay on the stability of the equilibrium points. We numerically solve the constructed model and corroborated the theoretical results.
Collapse
Affiliation(s)
- Sharmin Sultana
- Department of Mathematics, New Mexico Tech, New Mexico, 87801, USA
| | | | - Abraham J Arenas
- Departamento de Matemáticas y Estadística, Universidad de Córdoba, Montería, Colombia
| |
Collapse
|
7
|
Khan S, Dobrovolny HM. A study of the effects of age on the dynamics of RSV in animal models. Virus Res 2021; 304:198524. [PMID: 34329697 DOI: 10.1016/j.virusres.2021.198524] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Revised: 06/24/2021] [Accepted: 07/17/2021] [Indexed: 01/18/2023]
Abstract
Respiratory syncytial virus can cause severe illness and even death, particularly in infants. The increased severity of disease in young children is thought to be due to a lack of previous exposure to the virus as well as the limited immune response in infants. While studies have examined the clinical differences in disease between infants and adults, there has been limited examination of how the viral dynamics differ as infants develop. In this study, we apply a mathematical model to data from cotton rats and ferrets of different ages to assess how viral kinetics parameters change as the animals age. We find no clear trend in the viral decay rate, infecting time, and basic reproduction number as the animals age. We discuss possible reasons for the null result including the limited data, lack of detail of the mathematical model, and the limitations of animal models.
Collapse
Affiliation(s)
- Shaheer Khan
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX USA
| | - Hana M Dobrovolny
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX USA.
| |
Collapse
|
8
|
Yan AWC, Zhou J, Beauchemin CAA, Russell CA, Barclay WS, Riley S. Quantifying mechanistic traits of influenza viral dynamics using in vitro data. Epidemics 2020; 33:100406. [PMID: 33096342 DOI: 10.1016/j.epidem.2020.100406] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2019] [Revised: 07/10/2020] [Accepted: 09/04/2020] [Indexed: 11/28/2022] Open
Abstract
When analysing in vitro data, growth kinetics of influenza virus strains are often compared by computing their growth rates, which are sometimes used as proxies for fitness. However, analogous to mathematical models for epidemics, the growth rate can be defined as a function of mechanistic traits: the basic reproduction number (the average number of cells each infected cell infects) and the mean generation time (the average length of a replication cycle). Fitting a model to previously published and newly generated data from experiments in human lung cells, we compared estimates of growth rate, reproduction number and generation time for six influenza A strains. Of four strains in previously published data, A/Canada/RV733/2003 (seasonal H1N1) had the lowest basic reproduction number, followed by A/Mexico/INDRE4487/2009 (pandemic H1N1), then A/Indonesia/05/2005 (spill-over H5N1) and A/Anhui/1/2013 (spill-over H7N9). This ordering of strains was preserved for both generation time and growth rate, suggesting a positive biological correlation between these quantities which have not been previously observed. We further investigated these potential correlations using data from reassortant viruses with different internal proteins (from A/England/195/2009 (pandemic H1N1) and A/Turkey/05/2005 (H5N1)), and the same surface proteins (from A/Puerto Rico/8/34 (lab-adapted H1N1)). Similar correlations between traits were observed for these viruses, confirming our initial findings and suggesting that these patterns were related to the degree of human adaptation of internal genes. Also, the model predicted that strains with a smaller basic reproduction number, shorter generation time and slower growth rate underwent more replication cycles by the time of peak viral load, potentially accumulating mutations more quickly. These results illustrate the utility of mathematical models in inferring traits driving observed differences in in vitro growth of influenza strains.
Collapse
Affiliation(s)
- Ada W C Yan
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London W2 1PG, United Kingdom
| | - Jie Zhou
- Section of Virology, Department of Medicine, Imperial College London, London W2 1PG, United Kingdom
| | - Catherine A A Beauchemin
- Department of Physics, Ryerson University, 350 Victoria Street, Toronto, Ontario, M5B 2K3, Canada; Interdisciplinary Theoretical and Mathematical Sciences (iTHEMS), RIKEN, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Colin A Russell
- Laboratory of Applied Evolutionary Biology, Department of Medical Microbiology, Academic Medical Center, University of Amsterdam, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Wendy S Barclay
- Section of Virology, Department of Medicine, Imperial College London, London W2 1PG, United Kingdom
| | - Steven Riley
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London W2 1PG, United Kingdom.
| |
Collapse
|
9
|
Rowell CER, Dobrovolny HM. Energy Requirements for Loss of Viral Infectivity. FOOD AND ENVIRONMENTAL VIROLOGY 2020; 12:281-294. [PMID: 32757142 PMCID: PMC7405386 DOI: 10.1007/s12560-020-09439-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/21/2020] [Accepted: 07/31/2020] [Indexed: 06/11/2023]
Abstract
Outside the host, viruses will eventually lose their ability to infect cells due to conformational changes that occur to proteins on the viral capsid. In order to undergo a conformational change, these proteins require energy to activate the chemical reaction that leads to the conformational change. In this study, data from the literature is used to calculate the energy required for viral inactivation for a variety of different viruses by means of the Arrhenius equation. We find that some viruses (rhinovirus, poliovirus, human immunodeficiency virus, Alkhumra hemorrhagic fever virus, and hepatitis A virus) have high inactivation energies, indicative of breaking of a chemical double bond. We also find that several viruses (respiratory syncytial virus, poliovirus, and norovirus) have nonlinear Arrhenius plots, suggesting that there is more than a single pathway for inactivation of these viruses.
Collapse
Affiliation(s)
- Caroline E R Rowell
- Department of Chemistry, Wingate University, Hendersonville, NC, USA
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX, USA
| | - Hana M Dobrovolny
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX, USA.
| |
Collapse
|
10
|
Liao LE, Carruthers J, Smither SJ, Weller SA, Williamson D, Laws TR, García-Dorival I, Hiscox J, Holder BP, Beauchemin CAA, Perelson AS, López-García M, Lythe G, Barr JN, Molina-París C. Quantification of Ebola virus replication kinetics in vitro. PLoS Comput Biol 2020; 16:e1008375. [PMID: 33137116 PMCID: PMC7660928 DOI: 10.1371/journal.pcbi.1008375] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Revised: 11/12/2020] [Accepted: 09/23/2020] [Indexed: 12/17/2022] Open
Abstract
Mathematical modelling has successfully been used to provide quantitative descriptions of many viral infections, but for the Ebola virus, which requires biosafety level 4 facilities for experimentation, modelling can play a crucial role. Ebola virus modelling efforts have primarily focused on in vivo virus kinetics, e.g., in animal models, to aid the development of antivirals and vaccines. But, thus far, these studies have not yielded a detailed specification of the infection cycle, which could provide a foundational description of the virus kinetics and thus a deeper understanding of their clinical manifestation. Here, we obtain a diverse experimental data set of the Ebola virus infection in vitro, and then make use of Bayesian inference methods to fully identify parameters in a mathematical model of the infection. Our results provide insights into the distribution of time an infected cell spends in the eclipse phase (the period between infection and the start of virus production), as well as the rate at which infectious virions lose infectivity. We suggest how these results can be used in future models to describe co-infection with defective interfering particles, which are an emerging alternative therapeutic.
Collapse
Affiliation(s)
- Laura E. Liao
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM, USA 87545
| | - Jonathan Carruthers
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK
| | | | | | - Simon A. Weller
- Defence Science and Technology Laboratory, Salisbury SP4 0JQ, UK
| | - Diane Williamson
- Defence Science and Technology Laboratory, Salisbury SP4 0JQ, UK
| | - Thomas R. Laws
- Defence Science and Technology Laboratory, Salisbury SP4 0JQ, UK
| | - Isabel García-Dorival
- Institute of Infection and Global Health, University of Liverpool, Liverpool, L69 7BE, UK
| | - Julian Hiscox
- Institute of Infection and Global Health, University of Liverpool, Liverpool, L69 7BE, UK
| | - Benjamin P. Holder
- Department of Physics, Grand Valley State University, Allendale, MI, USA 49401
| | - Catherine A. A. Beauchemin
- Department of Physics, Ryerson University, Toronto, ON, Canada M5B 2K3
- Interdisciplinary Theoretical and Mathematical Sciences (iTHEMS) Research Program at RIKEN, Wako, Saitama, Japan, 351-0198
| | - Alan S. Perelson
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM, USA 87545
| | - Martín López-García
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK
| | - Grant Lythe
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK
| | - John N. Barr
- School of Molecular and Cellular Biology, University of Leeds, Leeds LS2 9JT, UK
| | - Carmen Molina-París
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK
- * E-mail:
| |
Collapse
|
11
|
Abstract
In this paper, we study and explore two control strategies to decrease the spread of Zika virus in the human and mosquito populations. The control strategies that we consider in this study are awareness and spraying campaigns. We solve several optimal control problems relying on a mathematical epidemic model of Zika that considers both human and mosquito populations. The first control strategy is broad and includes using information campaigns, encouraging people to use bednetting, wear long-sleeve shirts, or similar protection actions. The second control is more specific and relies on spraying insecticides. The control system relies on a Zika mathematical model with control functions. To develop the optimal control problem, we use Pontryagins’ maximum principle, which is numerically solved as a boundary value problem. For the mathematical model of the Zika epidemic, we use parameter values extracted from real data from an outbreak in Colombia. We study the effect of the costs related to the controls and infected populations. These costs are important in real life since they can change the outcomes and recommendations for health authorities dramatically. Finally, we explore different options regarding which control measures are more cost-efficient for society.
Collapse
|
12
|
Schmid H, Dobrovolny HM. An approximate solution of the interferon-dependent viral kinetics model of influenza. J Theor Biol 2020; 498:110266. [PMID: 32339545 DOI: 10.1016/j.jtbi.2020.110266] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Revised: 01/10/2020] [Accepted: 04/01/2020] [Indexed: 11/25/2022]
Abstract
The analysis of viral kinetics models is mostly achieved by numerical methods. We present an approach via a Magnus expansion that allows us to give an approximate solution to the interferon-dependent viral infection model of influenza which is compared with numerical results. The time of peak viral load is calculated from the approximation and stays within 10% in the studied range of interferon (IFN) efficacy ϵ ∈ [0, 1000]. We utilize our solution to interpret the effect of varying IFN efficacy, suggesting a competition between virions and interferon that can cause an additional peak in the usually exponential increase in the viral load.
Collapse
Affiliation(s)
- Harald Schmid
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, USA
| | - Hana M Dobrovolny
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, USA.
| |
Collapse
|
13
|
Pinky L, Gonzalez-Parra G, Dobrovolny HM. Effect of stochasticity on coinfection dynamics of respiratory viruses. BMC Bioinformatics 2019; 20:191. [PMID: 30991939 PMCID: PMC6469119 DOI: 10.1186/s12859-019-2793-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2018] [Accepted: 04/03/2019] [Indexed: 12/17/2022] Open
Abstract
Background Respiratory viral infections are a leading cause of mortality worldwide. As many as 40% of patients hospitalized with influenza-like illness are reported to be infected with more than one type of virus. However, it is not clear whether these infections are more severe than single viral infections. Mathematical models can be used to help us understand the dynamics of respiratory viral coinfections and their impact on the severity of the illness. Most models of viral infections use ordinary differential equations (ODE) that reproduce the average behavior of the infection, however, they might be inaccurate in predicting certain events because of the stochastic nature of viral replication cycle. Stochastic simulations of single virus infections have shown that there is an extinction probability that depends on the size of the initial viral inoculum and parameters that describe virus-cell interactions. Thus the coinfection dynamics predicted by the ODE might be difficult to observe in reality. Results In this work, a continuous-time Markov chain (CTMC) model is formulated to investigate probabilistic outcomes of coinfections. This CTMC model is based on our previous coinfection model, expressed in terms of a system of ordinary differential equations. Using the Gillespie method for stochastic simulation, we examine whether stochastic effects early in the infection can alter which virus dominates the infection. Conclusions We derive extinction probabilities for each virus individually as well as for the infection as a whole. We find that unlike the prediction of the ODE model, for similar initial growth rates stochasticity allows for a slower growing virus to out-compete a faster growing virus. Electronic supplementary material The online version of this article (10.1186/s12859-019-2793-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Lubna Pinky
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, USA.,Department of Pediatrics, University of Tennessee Health Science Center, Memphis, TN, USA
| | | | - Hana M Dobrovolny
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, USA.
| |
Collapse
|
14
|
Pinky L, González-Parra G, Dobrovolny HM. Superinfection and cell regeneration can lead to chronic viral coinfections. J Theor Biol 2019; 466:24-38. [PMID: 30639572 PMCID: PMC7094138 DOI: 10.1016/j.jtbi.2019.01.011] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2018] [Revised: 11/14/2018] [Accepted: 01/08/2019] [Indexed: 12/15/2022]
Abstract
Clinical researchers have found that coinfection of the respiratory tract can cause distinct disease outcome, sometimes leading to long-lasting infection, compared to single viral infection. The impact of coinfections in human respiratory tract have not yet been evaluated in either theoretical or experimental studies on a large scale. A few experiments confirm that different respiratory viruses can infect the same cell (superinfection). Superinfection alone cannot cause long-lasting viral coinfections. The combined mechanism of superinfection and cell regeneration provides a plausible mechanism for chronic viral coinfections.
Molecular diagnostic techniques have revealed that approximately 43% of the patients hospitalized with influenza-like illness are infected by more than one viral pathogen, sometimes leading to long-lasting infections. It is not clear how the heterologous viruses interact within the respiratory tract of the infected host to lengthen the duration of what are usually short, self-limiting infections. We develop a mathematical model which allows for single cells to be infected simultaneously with two different respiratory viruses (superinfection) to investigate the possibility of chronic coinfections. We find that a model with superinfection and cell regeneration has a stable chronic coinfection fixed point, while superinfection without cell regeneration produces only acute infections. This analysis suggests that both superinfection and cell regeneration are required to sustain chronic coinfection via this mechanism since coinfection is maintained by superinfected cells that allow slow-growing infections a chance to infect cells and continue replicating. This model provides a possible mechanism for chronic coinfection independent of any viral interactions via the immune response.
Collapse
Affiliation(s)
- Lubna Pinky
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX, United States.
| | - Gilberto González-Parra
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX, United States; Department of Mathematics, New Mexico Tech, Socorro, NM, United States
| | - Hana M Dobrovolny
- Department of Physics and Astronomy, Texas Christian University, Fort Worth, TX, United States
| |
Collapse
|