SARS-CoV-2 infection represents a global threat to human health. Various approaches were employed to reveal the pathogenetic mechanisms of COVID-19. Mathematical and computational modelling is a powerful tool to describe and analyze the infection dynamics in relation to a plethora of processes contributing to the observed disease phenotypes. In our study here, we formulate and calibrate a deterministic model of the SARS-CoV-2 life cycle. It provides a kinetic description of the major replication stages of SARS-CoV-2. Sensitivity analysis of the net viral progeny with respect to model parameters enables the identification of the life cycle stages that have the strongest impact on viral replication. These three most influential parameters are (i) degradation rate of positive sense vRNAs in cytoplasm (negative effect), (ii) threshold number of non-structural proteins enhancing vRNA transcription (negative effect), and (iii) translation rate of non-structural proteins (positive effect). The results of our analysis could be used for guiding the search for antiviral drug targets to combat SARS-CoV-2 infection.

There are many studies that model the within-host population dynamics of Human Immunodeficiency Virus Type 1 (HIV-1) infection. However, the within-infected-cell replication of HIV-1 remains to be not comprehensively addressed. There exist rather few quantitative models describing the regulation of the HIV-1 life cycle at the intracellular level. In treatment of HIV-1 infection, there remain issues related to side-effects and drug-resistance that require further search “...for new and better drugs, ideally targeting multiple independent steps in the HIV-1 replication cycle” (as highlighted recently by Tedbury & Freed, The Future of HIV-1 Therapeutics, 2015). High-resolution mathematical models of HIV-1 growth in infected cells provide an additional analytical tool in identifying novel drug targets. We formulate a high-dimensional model describing the biochemical reactions underlying the replication of HIV-1 in target cells. The model considers a nonlinear regulation of the transcription of HIV-1 mediated by Tat and the Rev-dependent transport of fully spliced and singly spliced transcripts from the nucleus to the cytoplasm. The model is calibrated using available information on the kinetics of various stages of HIV-1 replication. The sensitivity analysis of the model is performed to rank the biochemical processes of HIV-1 replication with respect to their impact on the net production of virions by one actively infected cell. The ranking of the sensitivity factors provides a quantitative basis for identifying novel targets for antiviral therapy. Our analysis suggests that HIV-1 assembly depending on Gag and Tat-Rev regulation of transcription and mRNA distribution present two most critical stages in HIV-1 replication that can be targeted to effectively control virus production. These processes are not covered by current antiretroviral treatments.

A calibrated mathematical model of antiviral immune response to SARS-CoV-2 infection is developed. The model considers the innate and antigen-specific responses to SARS-CoV-2 infection. Recently published data sets from human challenge studies with SARS-CoV-2 were used for parameter evaluation. The calibration of the mathematical model of SARS-CoV-2 infection is based on combining the parameter guesses from our earlier study of influenza A virus infection, some recent quantitative models of SARS-CoV-2 infection and clinical data-based parameter estimation of a subset of the model parameters. Hence, the calibrated mathematical model represents a theoretical exploration type of study, i.e., ‘in silico patient’ with mild-to-moderate severity phenotype, rather than a completely validated quantitative model of COVID-19 with respect to all its state-space variables. Understanding the regulation of multiple intertwined reaction components of the immune system is necessary for linking the kinetics of immune responses with the clinical phenotypes of COVID-19. Consideration of multiple immune reaction components in a single calibrated mathematical model allowed us to address some fundamental issues related to the pathogenesis of COVID-19, i.e., the sensitivity of the peak viral load to the parameters characterizing the antiviral specific response components, the kinetic coordination of the individual innate and adaptive immune responses, and the factors favoring a prolonged viral persistence. The model provides a tool for predicting the infectivity of patients, i.e., the amount of virus which is transmitted via droplets from the person infected with SARS-CoV-2, depending on the time of infection. The thresholds for variations of the innate and adaptive response parameters which lead to a prolonged persistence of SARS-CoV-2 due to the loss of a kinetic response synchrony/coordination between them were identified.

The human lymphatic system (HLS) is a complex network of lymphatic organs linked through the lymphatic vessels. We present a graph theory-based approach to model and analyze the human lymphatic network. Two different methods of building a graph are considered: the method using anatomical data directly and the method based on a system of rules derived from structural analysis of HLS. A simple anatomical data-based graph is converted to an oriented graph by quantifying the steady-state fluid balance in the lymphatic network with the use of the Poiseuille equation in vessels and the mass conservation at vessel junctions. A computational algorithm for the generation of the rule-based random graph is developed and implemented. Some fundamental characteristics of the two types of HLS graph models are analyzed using different metrics such as graph energy, clustering, robustness, etc.

Animal models of diseases, particularly mice, are considered to be the cornerstone for translational research in immunology. The aim of the present study is to model the geometry and analyze the network structure of the murine lymphatic system (LS). The algorithm for building the graph model of the LS makes use of anatomical data. To identify the edge directions of the graph model, a mass balance approach to lymph dynamics based on the Hagen–Poiseuille equation is applied. It is the first study in which a geometric model of the murine LS has been developed and characterized in terms of its structural organization and the lymph transfer function. Our study meets the demand for quantitative mechanistic approaches in the growing field of immunoengineering to utilize or exploit the lymphatic system for immunotherapy.

Replication of Human Immunodeficiency Virus type 1 (HIV) in infected CD4${}^{+}$ T cells represents a key driver of HIV infection. The HIV life cycle is characterised by the heterogeneity of infected cells with respect to multiplicity of infection and the variability in viral progeny. This heterogeneity can result from the phenotypic diversity of infected cells as well as from random effects and fluctuations in the kinetics of biochemical reactions underlying the virus replication cycle. To quantify the contribution of stochastic effects to the variability of HIV life cycle kinetics, we propose a high-resolution mathematical model formulated as a Markov chain jump process. The model is applied to generate the statistical characteristics of the (i) cell infection multiplicity, (ii) cooperative nature of viral replication, and (iii) variability in virus secretion by phenotypically identical cells. We show that the infection with a fixed number of viruses per CD4${}^{+}$ T cell leads to some heterogeneity of infected cells with respect to the number of integrated proviral genomes. The bottleneck factors in the virus production are identified, including the Gag-Pol proteins. Sensitivity analysis enables ranking of the model parameters with respect to the strength of their impact on the size of viral progeny. The first three globally influential parameters are the transport of genomic mRNA to membrane, the tolerance of transcription activation to Tat-mediated regulation, and the degradation of free and mature virions. These can be considered as potential therapeutical targets.