ABSTRACT
BACKGROUND: Dialysis is associated with a non-negligible rate of morbidity, requiring treatment customization. Many mathematical models have been developed describing solute kinetics during hemodialysis (HD) for an average uremic patient. The clinical need can be more adequately addressed by developing a patient-specific, multicompartmental model. MATERIALS AND METHODS: The data from 148 sessions (20 patients), recorded at the Regional Hospital of Lugano, Switzerland, were used to develop and validate the mathematical model. Diffusive and convective interactions among patient, dialysate and substitution fluid were considered. Three parameters, related to mass transfer efficiency at the cell membrane, at the dialyzer and at the capillary wall, were used to tune the model. The ability of the model to describe the clinical evolution of a specific HD session was evaluated by comparing model outputs with clinically acquired data on solutes and catabolite concentrations. RESULTS: The model developed in this study allows electrolyte and catabolite concentration trends during each HD session to be described. The errors obtained before the estimation of the patient-specific parameters drastically decrease after their identification. With the optimized model, plasmatic concentration trends can be described with an average percent error lower than 2.1% for Na+, Cl-, Ca2+ and HCO3-, lower than 5% for K+ and lower than 8% for urea. CONCLUSIONS: The peculiarity of the proposed model is the possibility it offers to perform a real-time simulation enabling quantitative appraisal of hematochemical quantities whose direct measurement is prohibitive. These will be beneficial to dialysis therapy planning, reducing intradialysis complications and improving patients' quality of life.
