Development of a novel mathematical model that explains SARS-CoV-2 infection dynamics in Caco-2 cells

Mathematical model explains differences in Omicron and Delta SARS-CoV-2 dynamics in Caco-2 and Calu-3 cells

Within-host infection dynamics of Omicron dramatically differs from previous variants of SARS-CoV-2. However, little is still known about which parameters of virus-cell interplay contribute to the observed attenuated replication and pathogenicity of Omicron. Mathematical models, often expressed as systems of differential equations, are frequently employed to study the infection dynamics of various viruses. Adopting such models for results of in vitro experiments can be beneficial in a number of aspects, such as model simplification (e.g., the absence of adaptive immune response and innate immunity cells), better measurement accuracy, and the possibility to measure additional data types in comparison with in vivo case. In this study, we consider a refinement of our previously developed and validated model based on a system of integro-differential equations. We fit the model to the experimental data of Omicron and Delta infections in Caco-2 (human intestinal epithelium model) and Calu-3 (lung epithelium model) cell lines. The data include known information on initial conditions, infectious virus titers, and intracellular viral RNA measurements at several time points post-infection. The model accurately explains the experimental data for both variants in both cell lines using only three variant- and cell-line-specific parameters. Namely, the cell entry rate is significantly lower for Omicron, and Omicron triggers a stronger cytokine production rate (i.e., innate immune response) in infected cells, ultimately making uninfected cells resistant to the virus. Notably, differences in only a single parameter (e.g., cell entry rate) are insufficient to obtain a reliable model fit for the experimental data.

Within-host delay differential model for SARS-CoV-2 kinetics with saturated antiviral responses

The present study discussed a model to describe the SARS-CoV-2 viral kinetics in the presence of saturated antiviral responses. A discrete-time delay was introduced due to the time required for uninfected epithelial cells to activate a suitable antiviral response by generating immune cytokines and chemokines. We examined the system's stability at each equilibrium point. A threshold value was obtained for which the system switched from stability to instability via a Hopf bifurcation. The length of the time delay has been computed, for which the system has preserved its stability. Numerical results show that the system was stable for the faster antiviral responses of epithelial cells to the virus concentration, i.e., quick antiviral responses stabilized patients' bodies by neutralizing the virus. However, if the antiviral response of epithelial cells to the virus increased, the system became unstable, and the virus occupied the whole body, which caused patients' deaths.