Coinfections of the Respiratory Tract: Viral Competition for Resources

Studies have shown that simultaneous infection of the respiratory tract with at least two viruses is common in hospitalized patients, although it is not clear whether these infections are more or less severe than single virus infections. We use a mathematical model to study the dynamics of viral coinfection of the respiratory tract in an effort to understand the kinetics of these infections. Specifically, we use our model to investigate coinfections of influenza, respiratory syncytial virus, rhinovirus, parainfluenza virus, and human metapneumovirus. Our study shows that during coinfections, one virus can block another simply by being the first to infect the available host cells; there is no need for viral interference through immune response interactions. We use the model to calculate the duration of detectable coinfection and examine how it varies as initial viral dose and time of infection are varied. We find that rhinovirus, the fastest-growing virus, reduces replication of the remaining viruses during a coinfection, while parainfluenza virus, the slowest-growing virus is suppressed in the presence of other viruses.

The Restitution Portrait:. A New Method for Investigating Rate-Dependent Restitution

INTRODUCTION

Electrical restitution, relating action potential duration (APD) to diastolic interval (DI), was believed to determine the stability of heart rhythm. However, recent studies demonstrate that stability also depends on long-term APD changes caused by memory. This study presents a new method for investigation of rate- and memory-dependent aspects of restitution and for assessment of mapping models of APD.

METHODS AND RESULTS

Bullfrog ventricular myocardium was paced with a "perturbed downsweep protocol." Starting from a basic cycle length (BCL) of 1,000 ms, the tissue was paced until steady state was achieved, followed by single beats of longer and shorter cycle lengths. BCL was decreased by 50 to 100 ms and the process repeated. All APDs were plotted as a function of the preceding DI, which allowed simultaneous observation of dynamic, S1-S2, and two constant-BCL restitution curves in a "restitution portrait." Responses were classified as 1:1 (stimulus:response), transient 2:2, or persistent 2:2 (alternans) and were related to the slopes of the restitution curves. None of these slopes approached unity for the persistent 2:2 response, demonstrating that the traditional restitution condition does not predict alternans. The restitution portrait was used to evaluate three mapping models of APD. The models with no memory and with one-beat memory did not produce restitution portraits similar to the experimental one. A model with two-beat memory produced a qualitatively similar portrait.

CONCLUSION

The restitution portrait allows a more comprehensive assessment of cardiac dynamics than methods used to date. Further study of models with memory may result in a clinical criterion for electrical instability.

Differences in predictions of ODE models of tumor growth: a cautionary example

BACKGROUND

While mathematical models are often used to predict progression of cancer and treatment outcomes, there is still uncertainty over how to best model tumor growth. Seven ordinary differential equation (ODE) models of tumor growth (exponential, Mendelsohn, logistic, linear, surface, Gompertz, and Bertalanffy) have been proposed, but there is no clear guidance on how to choose the most appropriate model for a particular cancer.

METHODS

We examined all seven of the previously proposed ODE models in the presence and absence of chemotherapy. We derived equations for the maximum tumor size, doubling time, and the minimum amount of chemotherapy needed to suppress the tumor and used a sample data set to compare how these quantities differ based on choice of growth model.

RESULTS

We find that there is a 12-fold difference in predicting doubling times and a 6-fold difference in the predicted amount of chemotherapy needed for suppression depending on which growth model was used.

CONCLUSION

Our results highlight the need for careful consideration of model assumptions when developing mathematical models for use in cancer treatment planning.

Assessing mathematical models of influenza infections using features of the immune response

The role of the host immune response in determining the severity and duration of an influenza infection is still unclear. In order to identify severity factors and more accurately predict the course of an influenza infection within a human host, an understanding of the impact of host factors on the infection process is required. Despite the lack of sufficiently diverse experimental data describing the time course of the various immune response components, published mathematical models were constructed from limited human or animal data using various strategies and simplifying assumptions. To assess the validity of these models, we assemble previously published experimental data of the dynamics and role of cytotoxic T lymphocytes, antibodies, and interferon and determined qualitative key features of their effect that should be captured by mathematical models. We test these existing models by confronting them with experimental data and find that no single model agrees completely with the variety of influenza viral kinetics responses observed experimentally when various immune response components are suppressed. Our analysis highlights the strong and weak points of each mathematical model and highlights areas where additional experimental data could elucidate specific mechanisms, constrain model design, and complete our understanding of the immune response to influenza.

SARS-CoV-2 coinfections: Could influenza and the common cold be beneficial?

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly spread around the world, causing serious illness and death and creating a heavy burden on the healthcare systems of many countries. Since the virus first emerged in late November 2019, its spread has coincided with peak circulation of several seasonal respiratory viruses, yet some studies have noted limited coinfections between SARS-CoV-2 and other viruses. We use a mathematical model of viral coinfection to study SARS-CoV-2 coinfections, finding that SARS-CoV-2 replication is easily suppressed by many common respiratory viruses. According to our model, this suppression is because SARS-CoV-2 has a lower growth rate (1.8/d) than the other viruses examined in this study. The suppression of SARS-CoV-2 by other pathogens could have implications for the timing and severity of a second wave.

Neuraminidase inhibitors for treatment of human and avian strain influenza: A comparative modeling study

Treatment of seasonal influenza viral infections using antivirals such as neuraminidase inhibitors (NAIs) has been proven effective if administered within 48h post-infection. However, there is growing evidence that antiviral treatment of infections with avian-derived strains even as late as 6 days post-infection (dpi) can significantly reduce infection severity and duration. Using a mathematical model of in-host influenza viral infections which can capture the kinetics of both a short-lived, typical, seasonal infection and a severe infection exhibiting sustained viral titer, we explore differences in the effects of NAI treatment on both types of influenza viral infections. Comparison of our model's behavior against experimental data from patients naturally infected with avian strains yields estimates for the times at which patients were infected that are consistent with those reported by the patients, and estimates of drug efficacies that are lower for patients who died than for those who recovered. In addition, our model suggests that the sustained, high, viral titers often seen in more severe influenza virus infections are the reason why antiviral treatment delayed by as much as 6 dpi will still lead to reduced viral titers and shortened illness. We conclude that NAIs may be an effective and beneficial treatment strategy against more severe strains of influenza virus characterized by high, sustained, viral titers. We believe that our mathematical model will be an effective tool in guiding treatment of severe influenza viral infections with antivirals.

Exploring cell tropism as a possible contributor to influenza infection severity

Several mechanisms have been proposed to account for the marked increase in severity of human infections with avian compared to human influenza strains, including increased cytokine expression, poor immune response, and differences in target cell receptor affinity. Here, the potential effect of target cell tropism on disease severity is studied using a mathematical model for in-host influenza viral infection in a cell population consisting of two different cell types. The two cell types differ only in their susceptibility to infection and rate of virus production. We show the existence of a parameter regime which is characterized by high viral loads sustained long after the onset of infection. This finding suggests that differences in cell tropism between influenza strains could be sufficient to cause significant differences in viral titer profiles, similar to those observed in infections with certain strains of influenza A virus. The two target cell mathematical model offers good agreement with experimental data from severe influenza infections, as does the usual, single target cell model albeit with biologically unrealistic parameters. Both models predict that while neuraminidase inhibitors and adamantanes are only effective when administered early to treat an uncomplicated seasonal infection, they can be effective against more severe influenza infections even when administered late.

Modeling the role of asymptomatics in infection spread with application to SARS-CoV-2

SARS-CoV-2 started causing infections in humans in late 2019 and has spread rapidly around the world. While the number of symptomatically infected and severely ill people is high and has overwhelmed the medical systems of many countries, there is mounting evidence that some of the rapid spread of this virus has been driven by asymptomatic infections. In this study, we use a compartmental mathematical model of a viral epidemic that includes asymptomatic infection to examine the role of asymptomatic individuals in the spread of the infection. We apply the model to epidemics in California, Florida, New York, and Texas, finding that asymptomatic infections far outnumber reported symptomatic infections at the peak of the epidemic in all four states. The model suggests that relaxing of social distancing measures too quickly could lead to a rapid rise in the number of cases, driven in part by asymptomatic infections.

Quantifying the effect of remdesivir in rhesus macaques infected with SARS-CoV-2

The world is in the midst of a pandemic caused by a novel coronavirus and is desperately searching for possible treatments. The antiviral remdesivir has shown some effectiveness against SARS-CoV-2 in vitro and in a recent animal study. We use data from a study of remdesivir in rhesus macaques to fit a viral kinetics model in an effort to determine the most appropriate mathematical descripton of the effect of remdesivir. We find statistically significant differences in the viral decay rate and use this to inform a possible mathematical formulation of the effect of remdesivir. Unfortunately, this model formulation suggests that the application of remdesivir will lengthen SARS-CoV-2 infections, putting into question its potential clinical benefit.

A comparison of RSV and influenza in vitro kinetic parameters reveals differences in infecting time

Influenza and respiratory syncytial virus (RSV) cause acute infections of the respiratory tract. Since the viruses both cause illnesses with similar symptoms, researchers often try to apply knowledge gleaned from study of one virus to the other virus. This can be an effective and efficient strategy for understanding viral dynamics or developing treatment strategies, but only if we have a full understanding of the similarities and differences between the two viruses. This study used mathematical modeling to quantitatively compare the viral kinetics of in vitro RSV and influenza virus infections. Specifically, we determined the viral kinetics parameters for RSV A2 and three strains of influenza virus, A/WSN/33 (H1N1), A/Puerto Rico/8/1934 (H1N1), and pandemic H1N1 influenza virus. We found that RSV viral titer increases at a slower rate and reaches its peak value later than influenza virus. Our analysis indicated that the slower increase of RSV viral titer is caused by slower spreading of the virus from one cell to another. These results provide estimates of dynamical differences between influenza virus and RSV and help provide insight into the virus-host interactions that cause observed differences in the time courses of the two illnesses in patients.

Modelling the emergence of influenza drug resistance: The roles of surface proteins, the immune response and antiviral mechanisms

The emergence of influenza drug resistance has become of particular interest as current planning for an influenza pandemic involves using massive amounts of antiviral drugs. We use semi-stochastic simulations to examine the emergence of drug resistant mutants during the course of a single infection within a patient in the presence and absence of antiviral therapy. We specifically examine three factors and their effect on the emergence of drug-resistant mutants: antiviral mechanism, the immune response, and surface proteins. We find that adamantanes, because they act at the start of the replication cycle to prevent infection, are less likely to produce drug-resistant mutants than NAIs, which act at the end of the replication cycle. A mismatch between surface proteins and internal RNA results in drug-resistant mutants being less likely to emerge, and emerging later in the infection because the mismatch gives antivirals a second chance to prevent propagation of the mutation. The immune response subdues slow growing infections, further reducing the probability that a drug resistant mutant will emerge and yield a drug-resistant infection. These findings improve our understanding of the factors that contribute to the emergence of drug resistance during the course of a single influenza infection.

Determining drug efficacy parameters for mathematical models of influenza

Antivirals are the first line of defence against influenza, so drug efficacy should be re-evaluated for each new strain. However, due to the time and expense involved in assessing the efficacy of drug treatments both in vitro and in vivo, treatment regimens are largely not re-evaluated even when strains are found to be resistant to antivirals. Mathematical models of the infection process can help in this assessment, but for accurate model predictions, we need to measure model parameters characterizing the efficacy of antivirals. We use computer simulations to explore whether in vitro experiments can be used to extract drug efficacy parameters for use in viral kinetics models. We find that the efficacy of neuraminidase inhibitors can be determined by measuring viral load during a single cycle assay, while the efficacy of adamantanes can be determined by measuring infected cells during the preparation stage for the single cycle assay.

Initial Inoculum and the Severity of COVID-19: A Mathematical Modeling Study of the Dose-Response of SARS-CoV-2 Infections

SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) causes a variety of responses in those who contract the virus, ranging from asymptomatic infections to acute respiratory failure and death. While there are likely multiple mechanisms triggering severe disease, one potential cause of severe disease is the size of the initial inoculum. For other respiratory diseases, larger initial doses lead to more severe outcomes. We investigate whether there is a similar link for SARS-CoV-2 infections using the combination of an agent-based model (ABM) and a partial differential equation model (PDM). We use the model to examine the viral time course for different sizes of initial inocula, generating dose-response curves for peak viral load, time of viral peak, viral growth rate, infection duration, and area under the viral titer curve. We find that large initial inocula lead to short infections, but with higher viral titer peaks; and that smaller initial inocula lower the viral titer peak, but make the infection last longer.

Investigating Different Mechanisms of Action in Combination Therapy for Influenza

Combination therapy for influenza can have several benefits, from reducing the emergence of drug resistant virus strains to decreasing the cost of antivirals. However, there are currently only two classes of antivirals approved for use against influenza, limiting the possible combinations that can be considered for treatment. However, new antivirals are being developed that target different parts of the viral replication cycle, and their potential for use in combination therapy should be considered. The role of antiviral mechanism of action in the effectiveness of combination therapy has not yet been systematically investigated to determine whether certain antiviral mechanisms of action pair well in combination. Here, we use a mathematical model of influenza to model combination treatment with antivirals having different mechanisms of action to measure peak viral load, infection duration, and synergy of different drug combinations. We find that antivirals that lower the infection rate and antivirals that increase the duration of the eclipse phase perform poorly in combination with other antivirals.

Quantifying rotavirus kinetics in the REH tumor cell line using in vitro data

Globally, rotavirus is the most common cause of diarrhea in children younger than 5 years of age, however, a quantitative understanding of the infection dynamics is still lacking. In this paper, we present the first study to extract viral kinetic parameters for in vitro rotavirus infections in the REH cell tumor line. We use a mathematical model of viral kinetics to extract parameter values by fitting the model to data from rotavirus infection of REH cells. While accurate results for some of the parameters of the mathematical model were not achievable due to its global non-identifiability, we are able to quantify approximately the time course of the infection for the first time. We also find that the basic reproductive number of rotavirus, which gives the number of secondary infections from a single infected cell, is much greater than one. Quantifying the kinetics of rotavirus leads not only to a better understanding of the infection process, but also provides a method for quantitative comparison of kinetics of different strains or for quantifying the effectiveness of antiviral treatment.

Epidemiological Consequences of Viral Interference: A Mathematical Modeling Study of Two Interacting Viruses

Some viruses have the ability to block or suppress growth of other viruses when simultaneously present in the same host. This type of viral interference or viral block has been suggested as a potential interaction between some respiratory viruses including SARS-CoV-2 and other co-circulating respiratory viruses. We explore how one virus' ability to block infection with another within a single host affects spread of the viruses within a susceptible population using a compartmental epidemiological model. We find that population-level effect of viral block is a decrease in the number of people infected with the suppressed virus. This effect is most pronounced when the viruses have similar epidemiological parameters. We use the model to simulate co-circulating epidemics of SARS-CoV-2 and influenza, respiratory syncytial virus (RSV), and rhinovirus, finding that co-circulation of SARS-CoV-2 and RSV causes the most suppression of SARS-CoV-2. Paradoxically, co-circulation of SARS-CoV-2 and influenza or rhinovirus results in almost no change in the SARS-CoV-2 epidemic, but causes a shift in the timing of the influenza and rhinovirus epidemics.

The in vivo efficacy of neuraminidase inhibitors cannot be determined from the decay rates of influenza viral titers observed in treated patients

Antiviral therapy is a first line of defence against new influenza strains. Current pandemic preparations involve stock- piling oseltamivir, an oral neuraminidase inhibitor (NAI), so rapidly determining the effectiveness of NAIs against new viral strains is vital for deciding how to use the stockpile. Previous studies have shown that it is possible to extract the drug efficacy of antivirals from the viral decay rate of chronic infections. In the present work, we use a nonlinear mathematical model representing the course of an influenza infection to explore the possibility of extracting NAI drug efficacy using only the observed viral titer decay rates seen in patients. We first show that the effect of a time-varying antiviral concentration can be accurately approximated by a constant efficacy. We derive a relationship relating the true treatment dose and time elapsed between doses to the constant drug dose required to approximate the time- varying dose. Unfortunately, even with the simplification of a constant drug efficacy, we show that the viral decay rate depends not just on drug efficacy, but also on several viral infection parameters, such as infection and production rate, so that it is not possible to extract drug efficacy from viral decay rate alone.

The impact of cell regeneration on the dynamics of viral coinfection

Many mathematical models of respiratory viral infections do not include regeneration of cells within the respiratory tract, arguing that the infection is resolved before there is significant cellular regeneration. However, recent studies have found that ∼40% of patients hospitalized with influenza-like illness are infected with at least two different viruses, which could potentially lead to longer-lasting infections. In these longer infections, cell regeneration might affect the infection dynamics, in particular, allowing for the possibility of chronic coinfections. Several mathematical models have been used to describe cell regeneration in infection models, though the effect of model choice on the predicted time course of viral coinfections is not clear. We investigate four mathematical models incorporating different mechanisms of cell regeneration during respiratory viral coinfection to determine the effect of cell regeneration on infection dynamics. We perform linear stability analysis for each of the models and find the steady states analytically. The analysis suggests that chronic illness is possible but only with one viral species; chronic coexistence of two different viral species is not possible with the regeneration models considered here.

The rate of viral transfer between upper and lower respiratory tracts determines RSV illness duration

Respiratory syncytial virus can lead to serious lower respiratory infection (LRI), particularly in children and the elderly. LRI can cause longer infections, lingering respiratory problems, and higher incidence of hospitalization. In this paper, we use a simplified ordinary differential equation model of viral dynamics to study the role of transport mechanisms in the occurrence of LRI. Our model uses two compartments to simulate the upper respiratory tract and the lower respiratory tract (LRT) and assumes two distinct types of viral transfer between the two compartments: diffusion and advection. We find that a range of diffusion and advection values lead to long-lasting infections in the LRT, elucidating a possible mechanism for the severe LRI infections observed in humans.

Modeling of fusion inhibitor treatment of RSV in African green monkeys

Respiratory syncytial virus (RSV) is a respiratory infection that can cause serious illness, particularly in infants. In this study, we test four different model implementations for the effect of a fusion inhibitor, including one model that combines different drug effects, by fitting the models to data from a study of TMC353121 in African green monkeys. We use mathematical modeling to estimate the drug efficacy parameters, ε_max, the maximum efficacy of the drug, and EC₅₀, the drug concentration needed to achieve half the maximum effect. We find that if TMC353121 is having multiple effects on viral kinetics, more detailed data, using different treatment delays, is needed to detect this effect.

High-resolution high-speed panoramic cardiac imaging system

A panoramic cardiac imaging system consisting of three high-speed CCD cameras has been developed to image the surface electrophysiology of a rabbit heart via fluorescence imaging using a voltage-sensitive fluorescent dye. A robust, unique mechanical system was designed to accommodate the three cameras and to adapt to the requirements of future experiments. A unified computer interface was created for this application - a single workstation controls all three CCD cameras, illumination, stimulation, and a stepping motor that rotates the heart. The geometric reconstruction algorithms were adapted from a previous cardiac imaging system. We demonstrate the system by imaging a polymorphic cardiac tachycardia.

Estimation of viral kinetics model parameters in young and aged SARS-CoV-2 infected macaques

The SARS-CoV-2 virus disproportionately causes serious illness and death in older individuals. In order to have the greatest impact in decreasing the human toll caused by the virus, antiviral treatment should be targeted to older patients. For this, we need a better understanding of the differences in viral dynamics between SARS-CoV-2 infection in younger and older adults. In this study, we use previously published averaged viral titre measurements from the nose and throat of SARS-CoV-2 infection in young and aged cynomolgus macaques to parametrize a viral kinetics model. We find that all viral kinetics parameters differ between young and aged macaques in the nasal passages, but that there are fewer differences in parameter estimates from the throat. We further use our parametrized model to study the antiviral treatment of young and aged animals, finding that early antiviral treatment is more likely to lead to a lengthening of the infection in aged animals, but not in young animals.

A quantitative assessment of dynamical differences of RSV infections in vitro and in vivo

Experimental results in vitro and in animal models are used to guide researchers in testing vaccines or treatment in humans. However, viral kinetics are different in vitro, in animals, and in humans, so it is sometimes difficult to translate results from one system to another. In this study, we use a mathematical model to fit experimental data from multiple cycle respiratory syncytial virus (RSV) infections in vitro, in african green monkey (AGM), and in humans in order to quantitatively compare viral kinetics in the different systems. We find that there are differences in viral clearance rate, productively infectious cell lifespan, and eclipse phase duration between in vitro and in vivo systems and among different in vivo systems. We show that these differences in viral kinetics lead to different estimates of drug effectiveness of fusion inhibitors in vitro and in AGM than in humans.

Superinfection and cell regeneration can lead to chronic viral coinfections

•

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.

Model Integration in Computational Biology: The Role of Reproducibility, Credibility and Utility

During the COVID-19 pandemic, mathematical modeling of disease transmission has become a cornerstone of key state decisions. To advance the state-of-the-art host viral modeling to handle future pandemics, many scientists working on related issues assembled to discuss the topics. These discussions exposed the reproducibility crisis that leads to inability to reuse and integrate models. This document summarizes these discussions, presents difficulties, and mentions existing efforts towards future solutions that will allow future model utility and integration. We argue that without addressing these challenges, scientists will have diminished ability to build, disseminate, and implement high-impact multi-scale modeling that is needed to understand the health crises we face.