Virus-mediated cell fusion of SARS-CoV-2 variants.

SARS-CoV-2 has the ability to form large multi-nucleated cells known as syncytia. Little is known about how syncytia affect the dynamics of the infection or severity of the disease. In this manuscript, we extend a mathematical model of cell-cell fusion assays to estimate both the syncytia formation rate and the average duration of the fusion phase for five strains of SARS-CoV-2. We find that the original Wuhan strain has the slowest rate of syncytia formation (6.4×10-4/h), but takes only 4.0 h to complete the fusion process, while the Alpha strain has the fastest rate of syncytia formation (0.36 /h), but takes 7.6 h to complete the fusion process. The Beta strain also has a fairly fast syncytia formation rate (9.7×10-2/h), and takes the longest to complete fusion (8.4 h). The D614G strain has a fairly slow syncytia formation rate (2.8×10-3/h), but completes fusion in 4.0 h. Finally, the Delta strain is in the middle with a syncytia formation rate of 3.2×10-2/h and a fusing time of 6.1 h. We note that for these SARS-CoV-2 strains, there appears to be a tradeoff between the ease of forming syncytia and the speed at which they complete the fusion process.

Different routes of infection of H5N1 lead to changes in infecting time.

Influenza virus infection can result in a wide range of clinical outcomes from asymptomatic infection to severe disease and death. While there are undoubtedly many factors that contribute to the severity of disease, one possible contributing factor that needs more investigation is the route of infection. In this study, we use previously published data from cynomolgus macaques infected with A/Vietnam/1203/04 (H5N1) via either aerosol (with and without bronchoalveolar lavages (BAL)) or a combined intrabronchial, oral, and intranasal route. We fit a mathematical model of within host viral kinetics to the data and find that when the macaques are infected via the aerosol route with subsequent BAL, the infecting time is significantly lower than for the other two groups. A lower infecting time indicates that the virus spreads from cell to cell more rapidly for aerosol infection with BAL than for the combined deposition or aerosol deposition alone. This study helps elucidate the mechanism behind different infection outcomes caused by differences in routes of infection.

Mathematical Modeling of Oncolytic Virus Therapy Reveals Role of the Immune Response.

Oncolytic adenoviruses (OAds) present a promising path for cancer treatment due to their selectivity in infecting and lysing tumor cells and their ability to stimulate the immune response. In this study, we use an ordinary differential equation (ODE) model of tumor growth inhibited by oncolytic virus activity to parameterize previous research on the effect of genetically re-engineered OAds in A549 lung cancer tumors in murine models. We find that the data are best fit by a model that accounts for an immune response, and that the immune response provides a mechanism for elimination of the tumor. We also find that parameter estimates for the most effective OAds share characteristics, most notably a high infection rate and low viral clearance rate, that might be potential reasons for these viruses' efficacy in delaying tumor growth. Further studies observing E1A and P19 recombined viruses in different tumor environments may further illuminate the extent of the effects of these genetic modifications.

Modeling of oncolytic viruses in a heterogeneous cell population to predict spread into non-cancerous cells.

New cancer treatment modalities that limit patient discomfort need to be developed. One possible new therapy is the use of oncolytic (cancer-killing) viruses. It is only recently that our ability to manipulate viral genomes has allowed us to consider deliberately infecting cancer patients with viruses. One key consideration is to ensure that the virus exclusively targets cancer cells and does not harm nearby non-cancerous cells. Here, we use a mathematical model of viral infection to determine the characteristics a virus would need to have in order to eradicate a tumor, but leave non-cancerous cells untouched. We conclude that the virus must differ in its ability to infect the two different cell types, with the infection rate of non-cancerous cells needing to be less than one hundredth of the infection rate of cancer cells. Differences in viral production rate or infectious cell death rate alone are not sufficient to protect non-cancerous cells.

Quantifying the effect of defective viral genomes in respiratory syncytial virus infections.

Defective viral genomes (DVGs) are viral genomes that contain only a partial viral RNA and so cannot replicate within cells on their own. If a cell containing DVGs is subsequently infected with a complete viral genome, the DVG can then use the missing proteins expressed by the full genome in order to replicate itself. Since the cell is producing defective genomes, it has less resources to produce fully functional virions and thus release of complete virions is often suppressed. Here, we use data from challenge studies of respiratory syncytial virus (RSV) in healthy adults to quantify the effect of DVGs. We use a mathematical model to fit the data, finding that late onset of DVGs and prolonged DVG detection are associated with lower infection rates and higher clearance rates. This result could have implications for the use of DVGs as a therapeutic.

Effect of cellular regeneration and viral transmission mode on viral spread.

Illness negatively affects all aspects of life and one major cause of illness is viral infections. Some viral infections can last for weeks; others, like influenza (the flu), can resolve quickly. During infections, uninfected cells can replicate in order to replenish the cells that have died due to the virus. Many viral models, especially those for short-lived infections like influenza, tend to ignore cellular regeneration since many think that uncomplicated influenza resolves much faster than cells regenerate. This research accounts for cellular regeneration, using an agent-based framework, and varies the regeneration rate in order to understand how cell regeneration affects viral infection dynamics under assumptions of different modes of transmission. We find that although the general trends in peak viral load, time of viral peak, and chronic viral load as regeneration rate changes are the same for cell-free or cell-to-cell transmission, the changes are more extreme for cell-to-cell transmission due to limited access of infected cells to newly generated cells.

Treatment of Respiratory Viral Coinfections.

With the advent of rapid multiplex PCR, physicians have been able to test for multiple viral pathogens when a patient presents with influenza-like illness. This has led to the discovery that many respiratory infections are caused by more than one virus. Antiviral treatment of viral coinfections can be complex because treatment of one virus will affect the time course of the other virus. Since effective antivirals are only available for some respiratory viruses, careful consideration needs to be given on the effect treating one virus will have on the dynamics of the other virus, which might not have available antiviral treatment. In this study, we use mathematical models of viral coinfections to assess the effect of antiviral treatment on coinfections. We examine the effect of the mechanism of action, relative growth rates of the viruses, and the assumptions underlying the interaction of the viruses. We find that high antiviral efficacy is needed to suppress both infections. If high doses of both antivirals are not achieved, then we run the risk of lengthening the duration of coinfection or even of allowing a suppressed virus to replicate to higher viral titers.

Estimation of virus-mediated cell fusion rate of SARS-CoV-2.

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.

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.

Using mathematical modeling to estimate time-independent cancer chemotherapy efficacy parameters.

One of the primary cancer treatment modalities is chemotherapy. Unfortunately, traditional anti-cancer drugs are often not selective and cause damage to healthy cells, leading to serious side effects for patients. For this reason more targeted therapeutics and drug delivery methods are being developed. The effectiveness of new treatments is initially determined via in vitro cell viability assays, which determine the IC 50  of the drug. However, these assays are known to result in estimates of IC 50  that depend on the measurement time, possibly resulting in over- or under-estimation of the IC 50 . Here, we test the possibility of using cell growth curves and fitting of mathematical models to determine the IC 50  as well as the maximum efficacy of a drug ( ε max ). We measured cell growth of MCF-7 and HeLa cells in the presence of different concentrations of doxorubicin and fit the data with a logistic growth model that incorporates the effect of the drug. This method leads to measurement time-independent estimates of IC 50  and ε max , but we find that ε max  is not identifiable. Further refinement of this methodology is needed to produce uniquely identifiable parameter estimates.

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.

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 study of the effects of age on the dynamics of RSV in animal models.

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.

The role of syncytia during viral infections.

Respiratory syncytial virus (RSV) is a common, contagious infection of the lungs and the respiratory tract. RSV is characterized by syncytia, which are multinuclear cells created by cells that have fused together. We use a mathematical model to study how different assumptions about the viral production and lifespan of syncytia change the resulting infection time course. We find that the effect of syncytia on viral titer is only apparent when the basic reproduction number for infection via syncytia formation is similar to the reproduction number for cell free viral transmission. When syncytia fusion rate is high, we find the presence of syncytia can lead to slowly growing infections if viral production is suppressed in syncytia. Our model provides insight into how the presence of syncytia can affect the time course of a viral infection.

Quantifying the effect of trypsin and elastase on in vitro SARS-CoV infections.

The SARS coronavirus (SARS-CoV) has the potential to cause serious disease that can spread rapidly around the world. Much of our understanding of SARS-CoV pathogenesis comes from in vitro experiments. Unfortunately, in vitro experiments cannot replicate all the complexity of the in vivo infection. For example, proteases in the respiratory tract cleave the SARS-CoV surface protein to facilitate viral entry, but these proteases are not present in vitro. Unfortunately, proteases might also have an effect on other parts of the replication cycle. Here, we use mathematical modeling to estimate parameters characterizing viral replication for SARS-CoV in the presence of trypsin or elastase, and in the absence of either. In addition to increasing the infection rate, the addition of trypsin and elastase causes lengthening of the eclipse phase duration and the infectious cell lifespan.

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.

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.

Energy Requirements for Loss of Viral Infectivity.

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.

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.

Understanding the effect of measurement time on drug characterization.

In order to determine correct dosage of chemotherapy drugs, the effect of the drug must be properly quantified. There are two important values that characterize the effect of the drug: εmax is the maximum possible effect of a drug, and IC50 is the drug concentration where the effect diminishes by half. There is currently a problem with the way these values are measured because they are time-dependent measurements. We use mathematical models to determine how the εmax and IC50 values depend on measurement time and model choice. Seven ordinary differential equation models (ODE) are used for the mathematical analysis; the exponential, Mendelsohn, logistic, linear, surface, Bertalanffy, and Gompertz models. We use the models to simulate tumor growth in the presence and absence of treatment with a known IC50 and εmax. Using traditional methods, we then calculate the IC50 and εmax values over fifty days to show the time-dependence of these values for all seven mathematical models. The general trend found is that the measured IC50 value decreases and the measured εmax increases with increasing measurement day for most mathematical models. Unfortunately, the measured values of IC50 and εmax rarely matched the values used to generate the data. Our results show that there is no optimal measurement time since models predict that IC50 estimates become more accurate at later measurement times while εmax is more accurate at early measurement times.