A Bayesian approach to estimate parameters of ordinary differential equation

Parameter Estimation for Kinetic Models of Chemical Reaction Networks from Partial Experimental Data of Species' Concentrations

The current manuscript addresses the problem of parameter estimation for kinetic models of chemical reaction networks from observed time series partial experimental data of species concentrations. It is demonstrated how the Kron reduction method of kinetic models, in conjunction with the (weighted) least squares optimization technique, can be used as a tool to solve the above-mentioned ill-posed parameter estimation problem. First, a new trajectory-independent measure is introduced to quantify the dynamical difference between the original mathematical model and the corresponding Kron-reduced model. This measure is then crucially used to estimate the parameters contained in the kinetic model so that the corresponding values of the species' concentrations predicted by the model fit the available experimental data. The new parameter estimation method is tested on two real-life examples of chemical reaction networks: nicotinic acetylcholine receptors and Trypanosoma brucei trypanothione synthetase. Both weighted and unweighted least squares techniques, combined with Kron reduction, are used to find the best-fitting parameter values. The method of leave-one-out cross-validation is utilized to determine the preferred technique. For nicotinic receptors, the training errors due to the application of unweighted and weighted least squares are 3.22 and 3.61 respectively, while for Trypanosoma synthetase, the application of unweighted and weighted least squares result in training errors of 0.82 and 0.70 respectively. Furthermore, the problem of identifiability of dynamical systems, i.e., the possibility of uniquely determining the parameters from certain types of output, has also been addressed.

Analysis of the COVID-19 pandemic using a compartmental model with time-varying parameters fitted by a genetic algorithm

Analyzing the COVID-19 pandemic is a critical factor in developing effective policies to deal with similar challenges in the future. However, many parameters (e.g., the actual number of infected people, the effectiveness of vaccination) are still subject to considerable debate because they are unobservable. To model a pandemic and estimate unobserved parameters, researchers use compartmental models. Most often, in such models, the transition rates are considered as constants, which allows simulating only one epidemiological wave. However, multiple waves have been reported for COVID-19 caused by different strains of the virus. This paper presents an approach based on the reconstruction of real distributions of transition rates using genetic algorithms, which makes it possible to create a model that describes several pandemic peaks. The model is fitted on registered COVID-19 cases in four countries with different pandemic control strategies (Germany, Sweden, UK, and US). Mean absolute percentage error (MAPE) was chosen as the objective function, the MAPE values of 2.168%, 2.096%, 1.208% and 1.703% were achieved for the listed countries, respectively. Simulation results are consistent with the empirical statistics of medical studies, which confirms the quality of the model. In addition to observables such as registered infected, the output of the model contains variables that cannot be measured directly. Among them are the proportion of the population protected by vaccines, the size of the exposed compartment, and the number of unregistered cases of COVID-19. According to the results, at the peak of the pandemic, between 14% (Sweden) and 25% (the UK) of the population were infected. At the same time, the number of unregistered cases exceeds the number of registered cases by 17 and 3.4 times, respectively. The average duration of the vaccine induced immune period is shorter than claimed by vaccine manufacturers, and the effectiveness of vaccination has declined sharply since the appearance of the Delta and Omicron strains. However, on average, vaccination reduces the risk of infection by about 65-70%.

Estimation of Ordinary Differential Equation Models for Gene Regulatory Networks Through Data Cloning

Ordinary differential equations (ODEs) are widely used for elucidating dynamic processes in various fields. One of the applications of ODEs is to describe dynamics of gene regulatory networks (GRNs), which is a critical step in understanding disease mechanisms. However, estimation of ODE models for GRNs is challenging because of inflexibility of the model and noisy data with complex error structures such as heteroscedasticity, correlations between genes, and time dependency. In addition, either a likelihood or Bayesian approach is commonly used for estimation of ODE models, but both approaches have benefits and drawbacks in their own right. Data cloning is a maximum likelihood (ML) estimation method through the Bayesian framework. Since it works in the Bayesian framework, it is free from local optimum problems that are common drawbacks of ML methods. Also, its inference is invariant for the selection of prior distributions, which is a major issue in Bayesian methods. This study proposes an estimation method of ODE models for GRNs through data cloning. The proposed method is demonstrated through simulation and it is applied to real gene expression time-course data.

Bayesian parameter estimation for phosphate dynamics during hemodialysis

Hyperphosphatemia in patients with renal failure is associated with increased vascular calcification and mortality. Hemodialysis is a conventional treatment for patients with hyperphosphatemia. Phosphate kinetics during hemodialysis may be described by a diffusion process and modeled by ordinary differential equations. We propose a Bayesian model approach for estimating patient-specific parameters for phosphate kinetics during hemodialysis. The Bayesian approach allows us to both analyze the full parameter space using uncertainty quantification and to compare two types of hemodialysis treatments, the conventional single-pass and the novel multiple-pass treatment. We validate and test our models on synthetic and real data. The results show limited identifiability of the model parameters when only single-pass data are available, and that the Bayesian model greatly reduces the relative standard deviation compared to existing estimates. Moreover, the analysis of the Bayesian models reveal improved estimates with reduced uncertainty when considering consecutive sessions and multiple-pass treatment compared to single-pass treatment.