On the Parametrization of Epidemiologic Models: Lessons from Modelling COVID-19 Epidemic

Numerous prediction models of SARS-CoV-2 pandemic were proposed in the past. Unknown parameters of these models are often estimated based on observational data. However, lag in case-reporting, changing testing policy or incompleteness of data lead to biased estimates. Moreover, parametrization is time-dependent due to changing age-structures, emerging virus variants, non-pharmaceutical interventions, and vaccination programs. To cover these aspects, we propose a principled approach to parametrize a SIR-type epidemiologic model by embedding it as a hidden layer into an input-output non-linear dynamical system (IO-NLDS). Observable data are coupled to hidden states of the model by appropriate data models considering possible biases of the data. This includes data issues such as known delays or biases in reporting. We estimate model parameters including their time-dependence by a Bayesian knowledge synthesis process considering parameter ranges derived from external studies as prior information. We applied this approach on a specific SIR-type model and data of Germany and Saxony demonstrating good prediction performances. Our approach can estimate and compare the relative effectiveness of non-pharmaceutical interventions and provide scenarios of the future course of the epidemic under specified conditions. It can be translated to other data sets, i.e., other countries and other SIR-type models

Fluid-Retention Side-Effects of the Chemotherapy-Supportive Treatment Interleukin-11: Mathematical Modelling of Putative Underlying Mechanisms

Interleukin-11 (IL-11) is a pleiotropic thrombopoietic cytokine and immune modulator, clinically approved for alleviation of chemotherapy-induced thrombocytopenia in non-myeloid malignancies. IL-11 therapy exerts fluid accumulation-associated adverse effects, complicating its administration and limiting its use. Implementation of standard biomathematical techniques to assess these effects is not possible, due to incomplete knowledge of the underlying mechanisms. This study investigates IL-11-induced blood volume expansion (BVE) by a new mathematical modelling methodology. Alternative models for BVE following IL-11 therapy were constructed, calibrated with clinical information and simulated in a number of treatment scenarios. The models demonstrated high compliance and were equally capable of reliably predicting BVE in a wide range of treatments, provided sufficient data. Model simulations indicate that frequent and low dose IL-11 regimens are favored for ensuring minimal fluid retention, upon the current IL-11 therapy

National and subnational short-term forecasting of COVID-19 in Germany and Poland during early 2021

We compare forecasts of weekly case and death numbers for COVID-19 in Germany and Poland based on 15 different modelling approaches. These cover the period from January to April 2021 and address numbers of cases and deaths one and two weeks into the future, along with the respective uncertainties. We find that combining different forecasts into one forecast can enable better predictions. However, case numbers over longer periods were challenging to predict. Additional data sources, such as information about different versions of the SARS-CoV-2 virus present in the population, might improve forecasts in the future

An Integrated Disease/Pharmacokinetic/Pharmacodynamic Model Suggests Improved Interleukin-21 Regimens Validated Prospectively for Mouse Solid Cancers

Interleukin (IL)-21 is an attractive antitumor agent with potent immunomodulatory functions. Yet thus far, the cytokine has yielded only partial responses in solid cancer patients, and conditions for beneficial IL-21 immunotherapy remain elusive. The current work aims to identify clinically-relevant IL-21 regimens with enhanced efficacy, based on mathematical modeling of long-term antitumor responses. For this purpose, pharmacokinetic (PK) and pharmacodynamic (PD) data were acquired from a preclinical study applying systemic IL-21 therapy in murine solid cancers. We developed an integrated disease/PK/PD model for the IL-21 anticancer response, and calibrated it using selected “training” data. The accuracy of the model was verified retrospectively under diverse IL-21 treatment settings, by comparing its predictions to independent “validation” data in melanoma and renal cell carcinoma-challenged mice (R2>0.90). Simulations of the verified model surfaced important therapeutic insights: (1) Fractionating the standard daily regimen (50 µg/dose) into a twice daily schedule (25 µg/dose) is advantageous, yielding a significantly lower tumor mass (45% decrease); (2) A low-dose (12 µg/day) regimen exerts a response similar to that obtained under the 50 µg/day treatment, suggestive of an equally efficacious dose with potentially reduced toxicity. Subsequent experiments in melanoma-bearing mice corroborated both of these predictions with high precision (R2>0.89), thus validating the model also prospectively in vivo. Thus, the confirmed PK/PD model rationalizes IL-21 therapy, and pinpoints improved clinically-feasible treatment schedules. Our analysis demonstrates the value of employing mathematical modeling and in silico-guided design of solid tumor immunotherapy in the clinic

Modeling individual time courses of thrombopoiesis during multi-cyclic chemotherapy.

BACKGROUND:Thrombocytopenia is a major side-effect of cytotoxic cancer therapies. The aim of precision medicine is to develop therapy modifications accounting for the individual's risk. METHODOLOGY/PRINCIPLE FINDINGS:To solve this task, we develop an individualized bio-mechanistic model of the dynamics of bone marrow thrombopoiesis, circulating platelets and therapy effects thereon. Comprehensive biological knowledge regarding cell differentiation, amplification, apoptosis rates, transition times and corresponding regulations are translated into ordinary differential equations. A model of osteoblast/osteoclast interactions was incorporated to mechanistically describe bone marrow support of quiescent cell stages. Thrombopoietin (TPO) as a major regulator is explicitly modelled including pharmacokinetics and-dynamics of TPO injections. Effects of cytotoxic drugs are modelled by transient depletions of proliferating cells. To calibrate the model, we used population data from the literature and close-meshed individual data of N = 135 high-grade non-Hodgkin's lymphoma patients treated with CHOP-like chemotherapies. To limit the number of free parameters, several parsimony assumptions were derived from biological data and tested via Likelihood methods. Heterogeneity of patients was explained by a few model parameters. The over-fitting issue of individual parameter estimation was successfully dealt with a virtual participation of each patient in population-based experiments. The model qualitatively and quantitatively explains a number of biological observations such as the role of osteoblasts in explaining long-term toxic effects, megakaryocyte-mediated feedback on stem cells, bi-phasic stimulation of thrombopoiesis by TPO, dynamics of megakaryocyte ploidies and non-exponential platelet degradation. Almost all individual time series could be described with high precision. We demonstrated how the model can be used to provide predictions regarding individual therapy adaptations. CONCLUSIONS:We propose a mechanistic thrombopoiesis model of unprecedented comprehensiveness in both, biological mechanisms considered and experimental data sets explained. Our innovative method of parameter estimation allows robust determinations of individual parameter settings facilitating the development of individual treatment adaptations during chemotherapy

Individual modelling of haematotoxicity with NARX neural networks: A knowledge transfer approach

Cytotoxic cancer therapy often results in dose-limiting haematotoxic side effects. Predicting an individual's risk is a major objective in precision medicine of cancer treatment. In this regard, patient heterogeneity presents a significant challenge. In this paper, we explore the use of hypothesis-free machine learning models based on recurrent nonlinear auto-regressive networks with exogenous inputs (NARX) as an approach to achieve this goal. Also, we propose a knowledge transfer approach to ameliorate the issue of sparse individual data, which typically hampers learning of individual networks. We demonstrate the feasibility of our approach based on a virtual patient population generated using a semi-mechanistic model of haematopoiesis and imposing different cytotoxic therapy scenarios on it. Employing different techniques of model optimisation, we derive robust and parsimonious individual networks with good generalisation performances. Moreover, we analyse in detail possible factors influencing the generalisation performance. Results suggest that our transfer learning approach using NARX networks can provide robust predictions of individual patient's response to treatment. As a practical perspective, we apply our approach to individual time series data of two patients with non-Hodgkin's lymphoma

National and subnational short-term forecasting of COVID-19 in Germany and Poland during early 2021

We compare forecasts of weekly case and death numbers for COVID-19 in Germany and Poland based on 15 different modelling approaches. These cover the period from January to April 2021 and address numbers of cases and deaths one and two weeks into the future, along with the respective uncertainties. We find that combining different forecasts into one forecast can enable better predictions. However, case numbers over longer periods were challenging to predict. Additional data sources, such as information about different versions of the SARS-CoV-2 virus present in the population, might improve forecasts in the future

European Covid-19 Forecast Hub

European Covid-19 Forecast Hub