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

Modeling the bystander effect during viral coinfection

Viral coinfections are responsible for a significant portion of cases of patients hospitalized with influenza-like illness. As our awareness of viral coinfections has increased, researchers have started to experimentally examine some of the virus-virus interactions underlying these infections. One mechanism of interaction between viruses is through the innate immune response. This seems to occur primarily through the interferon response, which generates an antiviral state in nearby uninfected cells, a phenomenon know as the bystander effect. Here, we develop a mathematical model of two viruses interacting through the bystander effect. We find that when the rate of removal of cells to the protected state is high, growth of the first virus is suppressed, while the second virus enjoys sole access to the protected cells, enhancing its growth. Conversely, growth of the second virus can be fully suppressed if its ability to infect the protected cells is limited.

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.

Flying under the radar - impact and factors influencing asymptomatic DENV infections

The clinical outcome of DENV and other Flaviviruses infections represents a spectrum of severity that ranges from mild manifestations to severe disease, which can ultimately lead to death. Nonetheless, most of these infections result in an asymptomatic outcome that may play an important role in the persistent circulation of these viruses. Also, although little is known about the mechanisms that lead to these asymptomatic infections, they are likely the result of a complex interplay between viral and host factors. Specific characteristics of the infecting viral strain, such as its replicating efficiency, coupled with host factors, like gene expression of key molecules involved in the immune response or in the protection against disease, are among crucial factors to study. This review revisits recent data on factors that may contribute to the asymptomatic outcome of the world's widespread DENV, highlighting the importance of silent infections in the transmission of this pathogen and the immune status of the host.

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.

A Trendline and Predictive Analysis of the First-Wave COVID-19 Infections in Malta

Following the first COVID-19 infected cases, Malta rapidly imposed strict lockdown measures, including restrictions on international travel, together with national social distancing measures, such as prohibition of public gatherings and closure of workplaces. The study aimed to elucidate the effect of the intervention and relaxation of the social distancing measures upon the infection rate by means of a trendline analysis of the daily case data. In addition, the study derived a predictive model by fitting historical data of the SARS-CoV-2 positive cases within a two-parameter Weibull distribution, whilst incorporating swab-testing rates, to forecast the infection rate at minute computational expense. The trendline analysis portrayed the wave of infection to fit within a tri-phasic pattern, where the primary phase was imposed with social measure interventions. Following the relaxation of public measures, the two latter phases transpired, where the two peaks resolved without further escalation of national measures. The derived forecasting model attained accurate predictions of the daily infected cases, attaining a high goodness-of-fit, utilising uncensored government-official infection-rate and swabbing-rate data within the first COVID-19 wave in Malta.

Mathematical modeling to study the impact of immigration on the dynamics of the COVID-19 pandemic: A case study for Venezuela

We propose two different mathematical models to study the effect of immigration on the COVID-19 pandemic. The first model does not consider immigration, whereas the second one does. Both mathematical models consider five different subpopulations: susceptible, exposed, infected, asymptomatic carriers, and recovered. We find the basic reproduction number R₀ using the next-generation matrix method for the mathematical model without immigration. This threshold parameter is paramount because it allows us to characterize the evolution of the disease and identify what parameters substantially affect the COVID-19 pandemic outcome. We focus on the Venezuelan scenario, where immigration and emigration have been important over recent years, particularly during the pandemic. We show that the estimation of the transmission rates of the SARS-CoV-2 are affected when the immigration of infected people is considered. This has an important consequence from a public health perspective because if the basic reproduction number is less than unity, we can expect that the SARS-CoV-2 would disappear. Thus, if the basic reproduction number is slightly above one, we can predict that some mild non-pharmaceutical interventions would be enough to decrease the number of infected people. The results show that the dynamics of the spread of SARS-CoV-2 through the population must consider immigration to obtain better insight into the outcomes and create awareness in the population regarding the population flow.

Theoretical Explanation for the Rarity of Antibody-Dependent Enhancement of Infection (ADE) in COVID-19

Global vaccination against the SARS-CoV-2 virus has proved to be highly effective. However, the possibility of antibody-dependent enhancement of infection (ADE) upon vaccination remains underinvestigated. Here, we aimed to theoretically determine conditions for the occurrence of ADE in COVID-19. We developed a series of mathematical models of antibody response: model Ab-a model of antibody formation; model Cv-a model of infection spread in the body; and a complete model, which combines the two others. The models describe experimental data on SARS-CoV and SARS-CoV-2 infections in humans and cell cultures, including viral load dynamics, seroconversion times and antibody concentration kinetics. The modelling revealed that a significant proportion of macrophages can become infected only if they bind antibodies with high probability. Thus, a high probability of macrophage infection and a sufficient amount of pre-existing antibodies are necessary for the development of ADE in SARS-CoV-2 infection. However, from the point of view of the dynamics of pneumocyte infection, the two cases where the body has a high concentration of preexisting antibodies and a high probability of macrophage infection and where there is a low concentration of antibodies in the body and no macrophage infection are indistinguishable. This conclusion could explain the lack of confirmed ADE cases for COVID-19.

SARS-CoV-2 Droplet and Airborne Transmission Heterogeneity

The spread dynamics of the SARS-CoV-2 virus have not yet been fully understood after two years of the pandemic. The virus's global spread represented a unique scenario for advancing infectious disease research. Consequently, mechanistic epidemiological theories were quickly dismissed, and more attention was paid to other approaches that considered heterogeneity in the spread. One of the most critical advances in aerial pathogens transmission was the global acceptance of the airborne model, where the airway is presented as the epicenter of the spread of the disease. Although the aerodynamics and persistence of the SARS-CoV-2 virus in the air have been extensively studied, the actual probability of contagion is still unknown. In this work, the individual heterogeneity in the transmission of 22 patients infected with COVID-19 was analyzed by close contact (cough samples) and air (environmental samples). Viral RNA was detected in 2/19 cough samples from patient subgroups, with a mean Ct (Cycle Threshold in Quantitative Polymerase Chain Reaction analysis) of 25.7 ± 7.0. Nevertheless, viral RNA was only detected in air samples from 1/8 patients, with an average Ct of 25.0 ± 4.0. Viral load in cough samples ranged from 7.3 × 105 to 8.7 × 108 copies/mL among patients, while concentrations between 1.1-4.8 copies/m3 were found in air, consistent with other reports in the literature. In patients undergoing follow-up, no viral load was found (neither in coughs nor in the air) after the third day of symptoms, which could help define quarantine periods in infected individuals. In addition, it was found that the patient's Ct should not be considered an indicator of infectiousness, since it could not be correlated with the viral load disseminated. The results of this work are in line with proposed hypotheses of superspreaders, which can attribute part of the heterogeneity of the spread to the oversized emission of a small percentage of infected people.

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.

Complications Post-COVID-19 and Risk Factors among Patients after Six Months of a SARS-CoV-2 Infection: A Population-Based Prospective Cohort Study

In October 2020, we conducted a population-based prospective cohort study to determine post-COVID-19 complications, recovery, return to usual health, and associated risk factors in 536 cases of COVID-19 outbreak in Borriana (Spain) by administering an epidemiological questionnaire via phone interviews. A total of 484 patients participated (90.3%), age mean 37.2 ± 17.1 years, and 301 females (62.2%). Mild illness was the most common COVID-19 manifestation. After six months, 160 patients (33.1%) suffered at least one complication post-COVID-19, and 47 (29.4%) of them sought medical assistance. The most frequent persistent symptoms were hair loss, fatigue, loss of smell or taste, and headache. Risk factors associated with a complication were female sex (adjusted relative risk, [aRR] = 1.93 95% confidence interval [CI] 1.41-2.65), age 35 years and above (aRR = 1.50 95% CI 1.14-1.99), B blood group (aRR = 1.51 95% CI 1.04-2.16), current smoker (RR = 1.61 95% CI 1.02-2.54), and at least a COVID-19 exposure (aRR = 2.13 95% CI 1.11-4.09). Male sex, age younger than 35 years, and low COVID-19 exposures were associated with better recovery and return to usual health. A third of patients presented persistent symptoms compatible with the long-COVID-19 syndrome. In conclusion, an active medical follow-up of post-COVID-19 patients must be implemented.

Mathematical modelling of SARS-CoV-2 infection of human and animal host cells reveals differences in the infection rates and delays in viral particle production by infected cells

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV -2), a causative agent of COVID-19 disease, poses a significant threat to public health. Since its outbreak in December 2019, Wuhan, China, extensive collection of diverse data from cell culture and animal infections as well as population level data from an ongoing pandemic, has been vital in assessing strategies to battle its spread. Mathematical modelling plays a key role in quantifying determinants that drive virus infection dynamics, especially those relevant for epidemiological investigations and predictions as well as for proposing efficient mitigation strategies. We utilized a simple mathematical model to describe and explain experimental results on viral replication cycle kinetics during SARS-CoV-2 infection of animal and human derived cell lines, green monkey kidney cells, Vero-E6, and human lung epithelium cells, A549-ACE2, respectively. We conducted cell infections using two distinct initial viral concentrations and quantified viral loads over time. We then fitted the model to our experimental data and quantified the viral parameters. We showed that such cellular tropism generates significant differences in the infection rates and incubation times of SARS-CoV-2, that is, the times to the first release of newly synthesised viral progeny by SARS-CoV-2-infected cells. Specifically, the rate at which A549-ACE2 cells were infected by SARS-CoV-2 was 15 times lower than that in the case of Vero-E6 cell infection and the duration of latent phase of A549-ACE2 cells was 1.6 times longer than that of Vero-E6 cells. On the other hand, we found no statistically significant differences in other viral parameters, such as viral production rate or infected cell death rate. Since in vitro infection assays represent the first stage in the development of antiviral treatments against SARS-CoV-2, discrepancies in the viral parameter values across different cell hosts have to be identified and quantified to better target vaccine and antiviral research.

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.

Molecular Communication Theoretical Modeling and Analysis of SARS-CoV2 Transmission in Human Respiratory System

Severe Acute Respiratory Syndrome-CoronaVirus 2 (SARS-CoV2) caused the ongoing pandemic. This pandemic devastated the world by killing more than a million people, as of October 2020. It is imperative to understand the transmission dynamics of SARS-CoV2 so that novel and interdisciplinary prevention, diagnostic, and therapeutic techniques could be developed. In this work, we model and analyze the transmission of SARS-CoV2 through the human respiratory tract from a molecular communication perspective. We consider that virus diffusion occurs in the mucus layer so that the shape of the tract does not have a significant effect on the transmission. Hence, this model reduces the inherent complexity of the human respiratory system. We further provide the impulse response of SARS-CoV2-ACE2 receptor binding event to determine the proportion of the virus population reaching different regions of the respiratory tract. Our findings confirm the results in the experimental literature on higher mucus flow rate causing virus migration to the lower respiratory tract. These results are especially important to understand the effect of SARS-CoV2 on the different human populations at different ages who have different mucus flow rates and ACE2 receptor concentrations in the different regions of the respiratory tract.

Analysis of Key Factors of a SARS-CoV-2 Vaccination Program: A Mathematical Modeling Approach

The administration of vaccines against the coronavirus disease 2019 (COVID-19) started in early December of 2020. Currently, there are only a few approved vaccines, each with different efficacies and mechanisms of action. Moreover, vaccination programs in different regions may vary due to differences in implementation, for instance, simply the availability of the vaccine. In this article, we study the impact of the pace of vaccination and the intrinsic efficacy of the vaccine on prevalence, hospitalizations, and deaths related to the SARS-CoV-2 virus. Then we study different potential scenarios regarding the burden of the COVID-19 pandemic in the near future. We construct a compartmental mathematical model and use computational methodologies to study these different scenarios. Thus, we are able to identify some key factors to reach the aims of the vaccination programs. We use some metrics related to the outcomes of the COVID-19 pandemic in order to assess the impact of the efficacy of the vaccine and the pace of the vaccine inoculation. We found that both factors have a high impact on the outcomes. However, the rate of vaccine administration has a higher impact in reducing the burden of the COVID-19 pandemic. This result shows that health institutions need to focus on increasing the vaccine inoculation pace and create awareness in the population about the importance of COVID-19 vaccines.

Current and prospective computational approaches and challenges for developing COVID-19 vaccines

SARS-CoV-2, which causes COVID-19, was first identified in humans in late 2019 and is a coronavirus which is zoonotic in origin. As it spread around the world there has been an unprecedented effort in developing effective vaccines. Computational methods can be used to speed up the long and costly process of vaccine development. Antigen selection, epitope prediction, and toxicity and allergenicity prediction are areas in which computational tools have already been applied as part of reverse vaccinology for SARS-CoV-2 vaccine development. However, there is potential for computational methods to assist further. We review approaches which have been used and highlight additional bioinformatic approaches and PK modelling as in silico methods which may be useful for SARS-CoV-2 vaccine design but remain currently unexplored. As more novel viruses with pandemic potential are expected to arise in future, these techniques are not limited to application to SARS-CoV-2 but also useful to rapidly respond to novel emerging viruses.

Impact of a New SARS-CoV-2 Variant on the Population: A Mathematical Modeling Approach

Several SARS-CoV-2 variants have emerged around the world, and the appearance of other variants depends on many factors. These new variants might have different characteristics that can affect the transmissibility and death rate. The administration of vaccines against the coronavirus disease 2019 (COVID-19) started in early December of 2020 and in some countries the vaccines will not soon be widely available. For this article, we studied the impact of a new more transmissible SARS-CoV-2 strain on prevalence, hospitalizations, and deaths related to the SARS-CoV-2 virus. We studied different scenarios regarding the transmissibility in order to provide a scientific support for public health policies and bring awareness of potential future situations related to the COVID-19 pandemic. We constructed a compartmental mathematical model based on differential equations to study these different scenarios. In this way, we are able to understand how a new, more infectious strain of the virus can impact the dynamics of the COVID-19 pandemic. We studied several metrics related to the possible outcomes of the COVID-19 pandemic in order to assess the impact of a higher transmissibility of a new SARS-CoV-2 strain on these metrics. We found that, even if the new variant has the same death rate, its high transmissibility can increase the number of infected people, those hospitalized, and deaths. The simulation results show that health institutions need to focus on increasing non-pharmaceutical interventions and the pace of vaccine inoculation since a new variant with higher transmissibility, such as, for example, VOC-202012/01 of lineage B.1.1.7, may cause more devastating outcomes in the population.