HPC Optimization Algorithm for Assessing Haemodynamic Parameters in Synthetic Patient Cohorts

Authors

  • Artem V. Rogov Moscow Institute of Physics and TechnologyWorld-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First State Medical University, Moscow, Russian Federation; Moscow Institute of Physics and Technology, Dolgoprudny, Russian Federation https://orcid.org/0009-0000-1610-2034
  • Timur M. Gamilov World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First State Medical University, Moscow, Russian Federation; Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow, Russian Federation https://orcid.org/0000-0002-1914-3859
  • Yaroslav Yu. Kirichenko I.M. Sechenov First Moscow State Medical University (Sechenov University), Moscow, Russian Federation https://orcid.org/0009-0001-7755-3995
  • Philipp Yu. Kopylov World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First State Medical University, Moscow, Russian Federation https://orcid.org/0000-0001-5124-6383
  • Sergey S. Simakov Moscow Institute of Physics and Technology, Dolgoprudny, Russian Federation; Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow, Russian Federation; I.M. Sechenov First Moscow State Medical University (Sechenov University), Moscow, Russian Federation https://orcid.org/0000-0003-3406-9623

DOI:

https://doi.org/10.14529/jsfi240408

Keywords:

Unscented Kalman Filter, hemodynamic model, inverse problem, synthetic data base, virtual population, parallel computing, high-performance computing

Abstract

A computational framework for the generation of a synthetic pulse wave database is developed. This framework demonstrates the feasibility of generating large-scale, high-fidelity virtual patient cohorts for biomedical research. Pulse waves are generated using a one-dimensional hemodynamic model of the systemic circulation coupled with a model of the left heart. Each virtual patient in the database is defined by a set of physiological parameters, including systolic and diastolic blood pressure, stroke volume, and heart rate. The parameters are optimized to match the desired outputs by solving an inverse problem using the Unscented Kalman Filter (UKF). The UKF is selected for its ability to accurately and efficiently estimate parameters in nonlinear systems. The generation of a single virtual patient requires between one and several hundred iterations of the UKF, depending on the complexity of the desired outputs. To meet the computational demands of generating a database with thousands of virtual patients, a computing cluster with 24 CPU nodes, each containing 52 cores, is employed. Two levels of parallelization are implemented, resulting in a speedup factor of 8.

References

Abohtyra, R.M., Vincent, T.L.: Nonlinear parameter and state estimation approach in end-stage kidney disease patients. bioRxiv (2023). https://doi.org/10.1101/2022.04.02.486844

Allen, R.J., Rieger, T.R., Musante, C.J.: Efficient generation and selection of virtual populations in quantitative systems pharmacology models. CPT Pharmacometrics Syst Pharmacol 5(3), 140–146 (2016). https://doi.org/10.1002/psp4.12063

Arthurs, C., Xiao, N., Moireau, P., et al.: A flexible framework for sequential estimation of model parameters in computational hemodynamics. Advances in Modeling and Simulation in Engineering Sciences 7, 48 (2020).