A Bayesian approach for the identification of patient-specific parameters in a dialysis kinetic model

Hemodialysis is the most common therapy to treat renal insufficiency. However, notwithstanding the recent improvements, hemodialysis is still associated with a non-negligible rate of comorbidities, which could be reduced by customizing the treatment. Many differential compartment models have been developed to describe the mass balance of blood electrolytes and catabolites during hemodialysis, with the goal of improving and controlling hemodialysis sessions. However, these models often refer to an average uremic patient, while on the contrary the clinical need for customization requires patient-specific models. In this work, we assume that the customization can be obtained by means of patient-specific model parameters. We propose and validate a Bayesian approach to estimate the patient-specific parameters of a multi-compartment model, and to predict the single patient’s response to the treatment, in order to prevent intra-dialysis complications. The likelihood function is obtained by means of a discretized version of the multi-compartment model, where the discretization is in terms of a Runge–Kutta method to guarantee convergence, and the posterior densities of model parameters are obtained through Markov Chain Monte Carlo simulation. Results show fair estimations and the applicability in the clinical practice.

[1]  J. Kooman,et al.  On-line monitoring of electrolytes in hemodialysis: on the road towards individualizing treatment , 2016, Expert review of medical devices.

[2]  Michele Conti,et al.  A Clinically Applicable Stochastic Approach for Noninvasive Estimation of Aortic Stiffness Using Computed Tomography Data , 2015, IEEE Transactions on Biomedical Engineering.

[3]  F. Locatelli,et al.  Optimization of sodium removal in paired filtration dialysis by single pool sodium and conductivity kinetic models. , 1997, Blood purification.

[4]  J T Daugirdas,et al.  Compartment Effects in Hemodialysis , 2001, Seminars in dialysis.

[5]  A. Peixoto,et al.  Methods of Assessment of Volume Status and Intercompartmental Fluid Shifts in Hemodialysis Patients: Implications in Clinical Practice , 2004, Seminars in dialysis.

[6]  Fabrizio Ruggeri,et al.  Inertance estimation in a lumped-parameter hydraulic simulator of human circulation. , 2013, Journal of biomechanical engineering.

[7]  W. Clark,et al.  Quantifying the effect of changes in the hemodialysis prescription on effective solute removal with a mathematical model. , 1999, Journal of the American Society of Nephrology : JASN.

[8]  F. Locatelli,et al.  The role of blood volume reduction in the genesis of intradialytic hypotension. , 2002, American journal of kidney diseases : the official journal of the National Kidney Foundation.

[9]  M J Lysaght,et al.  Kinetics of hemodiafiltration. I. In vitro transport characteristics of a hollow-fiber blood ultrafilter. , 1975, The Journal of laboratory and clinical medicine.

[10]  F. Ruggeri,et al.  Bayesian Estimation of Thermal Conductivity and Temperature Profile in a Homogeneous Mass , 2014 .

[11]  M. Gabriela M. Gomes,et al.  A Bayesian Framework for Parameter Estimation in Dynamical Models , 2011, PloS one.

[12]  A. Cheung,et al.  A Pseudo‐One Compartment Model of Phosphorus Kinetics During Hemodialysis: Further Supporting Evidence , 2017, Artificial organs.

[13]  F. Gotch,et al.  Mathematic modeling of dialysis therapy. , 1980, Kidney international. Supplement.

[14]  Domenico Vito,et al.  Patient-Specific Modeling of Multicompartmental Fluid and Mass Exchange during Dialysis , 2016, The International journal of artificial organs.

[15]  R. Fluck,et al.  Potassium kinetics during hemodialysis , 2015, Hemodialysis international. International Symposium on Home Hemodialysis.

[16]  J T Daugirdas,et al.  Dialysis hypotension: a hemodynamic analysis. , 1991, Kidney international.

[17]  J. Leypoldt Kinetics of β2‐Microglobulin and Phosphate during Hemodialysis: Effects of Treatment Frequency and Duration , 2005, Seminars in dialysis.

[18]  J. Michael Textbook of Medical Physiology , 2005 .

[19]  S. Prakash,et al.  Central, Peripheral, and Other Blood Volume Changes During Hemodialysis , 2002, ASAIO journal.

[20]  M. Ursino,et al.  Prediction of Solute Kinetics, Acid-Base Status, and Blood Volume Changes During Profiled Hemodialysis , 2000, Annals of Biomedical Engineering.

[21]  Julia L. Blanchard,et al.  Parameter uncertainty of a dynamic multispecies size spectrum model1 , 2016 .

[22]  T. Faniran Numerical Solution of Stochastic Differential Equations , 2015 .

[23]  M. Tapolyai,et al.  Which fluid space is affected by ultrafiltration during hemodiafiltration? , 2014, Hemodialysis international. International Symposium on Home Hemodialysis.

[24]  F. Locatelli,et al.  Assessment of intradialysis calcium mass balance by a single pool variable‐volume calcium kinetic model , 2017, Hemodialysis international. International Symposium on Home Hemodialysis.

[25]  S Stefoni,et al.  A Mathematical Model for the Prediction of Solute Kinetics, Osmolarity and Fluid Volume Changes during Hemodiafiltration with On-Line Regeneration of Ultrafiltrate (HFR) , 2006, The International journal of artificial organs.

[26]  F. Ruggeri,et al.  Nonlinear functional response parameter estimation in a stochastic predator-prey model. , 2011, Mathematical biosciences and engineering : MBE.

[27]  Ettore Lanzarone,et al.  A Bayesian estimation approach for the mortality in a stage-structured demographic model , 2017, Journal of mathematical biology.

[28]  S. Santos,et al.  Sodium Balance in Maintenance Hemodialysis , 2010, Seminars in dialysis.

[29]  A. Cheung,et al.  Prescribing hemodialysis using a weekly urea mass balance model. , 1991, Blood purification.

[30]  B. Øksendal Stochastic differential equations : an introduction with applications , 1987 .