RUS  ENG
Full version
JOURNALS // Computer Research and Modeling // Archive

Computer Research and Modeling, 2023 Volume 15, Issue 1, Pages 111–124 (Mi crm1047)

ANALYSIS AND MODELING OF COMPLEX LIVING SYSTEMS

Bistability and damped oscillations in the homogeneous model of viral infection

A. A. Tokarevabc, N. O. Rodinb, V. A. Vol'pertbd

a N. N. Semenov Federal Research Center for Chemical Physics RAS, 4/1 Kosygina st., 119991 Moscow, Russia
b Ðeoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya st., Moscow, 117198, Russia
c Bukhara Engineering Technological Institute, 15 Murtazoyev st., Bukhara, 200100, Uzbekistan
d Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622 Villeurbanne, France

Abstract: The development of a viral infection in the organism is a complex process which depends on the competition race between virus replication in the host cells and the immune response. To study different regimes of infection progression, we analyze the general mathematical model of immune response to viral infection. The model consists of two ODEs for virus and immune cells non-dimensionalized concentrations. The proliferation rate of immune cells in the model is represented by a bell-shaped function of the virus concentration. This function increases for small virus concentrations describing the antigen-stimulated clonal expansion of immune cells, and decreases for sufficiently high virus concentrations describing down-regulation of immune cells proliferation by the infection. Depending on the virus virulence, strength of the immune response, and the initial viral load, the model predicts several scenarios: (a) infection can be completely eliminated, (b) it can remain at a low level while the concentration of immune cells is high; (c) immune cells can be essentially exhausted, or (d) completely exhausted, which is accompanied (c, d) by high virus concentration. The analysis of the model shows that virus concentration can oscillate as it gradually converges to its equilibrium value. We show that the considered model can be obtained by the reduction of a more general model with an additional equation for the total viral load provided that this equation is fast. In the case of slow kinetics of the total viral load, this more general model should be used.

Keywords: dynamics of viral infection, immune response, bistability, damped oscillations, mathematical modeling, qualitative analysis of ordinary differential equations.

UDC: 577.31, 577.38, 51-76

Received: 25.05.2022
Revised: 06.10.2022
Accepted: 13.12.2022

Language: English

DOI: 10.20537/2076-7633-2023-15-1-111-124



© Steklov Math. Inst. of RAS, 2026