학술논문

A Mathematical Model for Simulation of Vasoplegic Shock and Vasopressor Therapy
Document Type
Periodical
Source
IEEE Transactions on Biomedical Engineering IEEE Trans. Biomed. Eng. Biomedical Engineering, IEEE Transactions on. 70(5):1565-1574 May, 2023
Subject
Bioengineering
Computing and Processing
Components, Circuits, Devices and Systems
Communication, Networking and Broadcast Technologies
Mathematical models
Baroreflex
Modulation
Computational modeling
Physiology
Biomedical monitoring
Drugs
Vasopressor
vasoplegia
phenylephrine
digital twin
autonomous critical care
physiological closed-loop control
Language
ISSN
0018-9294
1558-2531
Abstract
Objective: To develop a high-fidelity mathematical model intended to replicate the cardiovascular (CV) responses of a critically ill patient to vasoplegic shock-induced hypotension and vasopressor therapy. Methods: The mathematical model consists of a lumped-parameter CV physiology model with baroreflex modulation feedback and a phenomenological dynamic dose-response model of a vasopressor. The adequacy of the proposed mathematical model was investigated using an experimental dataset acquired from 10 pigs receiving phenylephrine (PHP) therapy after vasoplegic shock induced via sodium nitroprusside (SNP). Results: Upon calibration, the mathematical model could (i) faithfully replicate the effects of PHP on dynamic changes in blood pressure (BP), cardiac output (CO), and systemic vascular resistance (SVR) (root-mean-squared errors between measured and calibrated mathematical responses: mean arterial BP 2.5+/−1.0 mmHg, CO 0.2+/−0.1 lpm, SVR 2.4+/−1.5 mmHg/lpm; r value: mean arterial BP 0.96+/−0.01, CO 0.65+/−0.45, TPR 0.92+/−0.10) and (ii) predict physiologically plausible behaviors of unmeasured internal CV variables as well as secondary baroreflex modulation effects. Conclusion: This mathematical model is perhaps the first of its kind that can comprehensively replicate both primary (i.e., direct) and secondary (i.e., baroreflex modulation) effects of a vasopressor drug on an array of CV variables, rendering it ideally suited to pre-clinical virtual evaluation of the safety and efficacy of closed-loop control algorithms for autonomous vasopressor administration once it is extensively validated. Significance: This mathematical model architecture incorporating both direct and baroreflex modulation effects may generalize to serve as part of an effective platform for high-fidelity in silico simulation of CV responses to vasopressors during vasoplegic shock.