Data assimilation and modelling of patient-specific single-ventricle physiology with and without valve regurgitation

Simulating impaired left ventricular-arterial coupling in aging and disease: a systematic review

Aortic stenosis, hypertension, and left ventricular hypertrophy often coexist in the elderly, causing a detrimental mismatch in coupling between the heart and vasculature known as ventricular-vascular (VA) coupling. Impaired left VA coupling, a critical aspect of cardiovascular dysfunction in aging and disease, poses significant challenges for optimal cardiovascular performance. This systematic review aims to assess the impact of simulating and studying this coupling through computational models. By conducting a comprehensive analysis of 34 relevant articles obtained from esteemed databases such as Web of Science, Scopus, and PubMed until July 14, 2022, we explore various modeling techniques and simulation approaches employed to unravel the complex mechanisms underlying this impairment. Our review highlights the essential role of computational models in providing detailed insights beyond clinical observations, enabling a deeper understanding of the cardiovascular system. By elucidating the existing models of the heart (3D, 2D, and 0D), cardiac valves, and blood vessels (3D, 1D, and 0D), as well as discussing mechanical boundary conditions, model parameterization and validation, coupling approaches, computer resources and diverse applications, we establish a comprehensive overview of the field. The descriptions as well as the pros and cons on the choices of different dimensionality in heart, valve, and circulation are provided. Crucially, we emphasize the significance of evaluating heart-vessel interaction in pathological conditions and propose future research directions, such as the development of fully coupled personalized multidimensional models, integration of deep learning techniques, and comprehensive assessment of confounding effects on biomarkers.

An In Silico Model for Predicting the Efficacy of Edge-to-Edge Repair for Mitral Regurgitation

In recent years, transcatheter edge-to-edge repair (TEER) has been widely adopted as an effective treatment for mitral regurgitation (MR). The aim of this study is to develop a personalized in silico model to predict the effect of edge-to-edge repair in advance to the procedure for each individual patient. For this purpose, we propose a combination of a valve deformation model for computing the mitral valve (MV) orifice area (MVOA) and a lumped parameter model for the hemodynamics, specifically mitral regurgitation volume (RVol). Although we cannot obtain detailed information on the three-dimensional flow field near the mitral valve, we can rapidly simulate the important medical parameters for the clinical decision support. In the present method, we construct the patient-specific pre-operative models by using the parameter optimization and then simulate the postoperative state by applying the additional clipping condition. The computed preclip MVOAs show good agreement with the clinical measurements, and the correlation coefficient takes 0.998. In addition, the MR grade in terms of RVol also has good correlation with the grade by ground truth MVOA. Finally, we try to investigate the applicability for the predicting the postclip state. The simulated valve shapes clearly show the well-known double orifice and the improvement of the MVOA, compared with the preclip state. Similarly, we confirmed the improved reverse flow and MR grade in terms of RVol. A total computational time is approximately 8 h by using general-purpose PC. These results obviously indicate that the present in silico model has good capability for the assessment of edge-to-edge repair.

Computational modeling of ventricular-ventricular interactions suggest a role in clinical conditions involving heart failure

Introduction: The left (LV) and right (RV) ventricles are linked biologically, hemodynamically, and mechanically, a phenomenon known as ventricular interdependence. While LV function has long been known to impact RV function, the reverse is increasingly being realized to have clinical importance. Investigating ventricular interdependence clinically is challenging given the invasive measurements required, including biventricular catheterization, and confounding factors such as comorbidities, volume status, and other aspects of subject variability. Methods: Computational modeling allows investigation of mechanical and hemodynamic interactions in the absence of these confounding factors. Here, we use a threesegment biventricular heart model and simple circulatory system to investigate ventricular interdependence under conditions of systolic and diastolic dysfunction of the LV and RV in the presence of compensatory volume loading. We use the end-diastolic pressure-volume relationship, end-systolic pressure-volume relationship, Frank Starling curves, and cardiac power output as metrics. Results: The results demonstrate that LV systolic and diastolic dysfunction lead to RV compensation as indicated by increases in RV power. Additionally, RV systolic and diastolic dysfunction lead to impaired LV filling, interpretable as LV stiffening especially with volume loading to maintain systemic pressure. Discussion: These results suggest that a subset of patients with intact LV systolic function and diagnosed to have impaired LV diastolic function, categorized as heart failure with preserved ejection fraction (HFpEF), may in fact have primary RV failure. Application of this computational approach to clinical data sets, especially for HFpEF, may lead to improved diagnosis and treatment strategies and consequently improved outcomes.

Monitoring variability in parameter estimates for lumped parameter models of the systemic circulation using longitudinal hemodynamic measurements

BACKGROUND

Physics-based cardiovascular models are only recently being considered for disease diagnosis or prognosis in clinical settings. These models depend on parameters representing the physical and physiological properties of the modeled system. Personalizing these parameters may give insight into the specific state of the individual and etiology of disease. We applied a relatively fast model optimization scheme based on common local optimization methods to two model formulations of the left ventricle and systemic circulation. One closed-loop model and one open-loop model were applied. Intermittently collected hemodynamic data from an exercise motivation study were used to personalize these models for data from 25 participants. The hemodynamic data were collected for each participant at the start, middle and end of the trial. We constructed two data sets for the participants, both consisting of systolic and diastolic brachial pressure, stroke volume, and left-ventricular outflow tract velocity traces paired with either the finger arterial pressure waveform or the carotid pressure waveform.

RESULTS

We examined the feasibility of separating parameter estimates for the individual from population estimates by assessing the variability of estimates using the interquartile range. We found that the estimated parameter values were similar for the two model formulations, but that the systemic arterial compliance was significantly different ([Formula: see text]) depending on choice of pressure waveform. The estimates of systemic arterial compliance were on average higher when using the finger artery pressure waveform as compared to the carotid waveform.

CONCLUSIONS

We found that for the majority of participants, the variability of parameter estimates for a given participant on any measurement day was lower than the variability both across all measurement days combined for one participant, and for the population. This indicates that it is possible to identify individuals from the population, and that we can distinguish different measurement days for the individual participant by parameter values using the presented optimization method.

Inverse problems in blood flow modeling: A review

Mathematical and computational modeling of the cardiovascular system is increasingly providing non-invasive alternatives to traditional invasive clinical procedures. Moreover, it has the potential for generating additional diagnostic markers. In blood flow computations, the personalization of spatially distributed (i.e., 3D) models is a key step which relies on the formulation and numerical solution of inverse problems using clinical data, typically medical images for measuring both anatomy and function of the vasculature. In the last years, the development and application of inverse methods has rapidly expanded most likely due to the increased availability of data in clinical centers and the growing interest of modelers and clinicians in collaborating. Therefore, this work aims to provide a wide and comparative overview of literature within the last decade. We review the current state of the art of inverse problems in blood flows, focusing on studies considering fully dimensional fluid and fluid-solid models. The relevant physical models and hemodynamic measurement techniques are introduced, followed by a survey of mathematical data assimilation approaches used to solve different kinds of inverse problems, namely state and parameter estimation. An exhaustive discussion of the literature of the last decade is presented, structured by types of problems, models and available data.

Multiscale modelling of Potts shunt as a potential palliative treatment for suprasystemic idiopathic pulmonary artery hypertension: a paediatric case study

Potts shunt (PS) was suggested as palliation for patients with suprasystemic pulmonary arterial hypertension (PAH) and right ventricular (RV) failure. PS, however, can result in poorly understood mortality. Here, a patient-specific geometrical multiscale model of PAH physiology and PS is developed for a paediatric PAH patient with stent-based PS. In the model, 7.6mm-diameter PS produces near-equalisation of the aortic and PA pressures and [Formula: see text] (oxygenated vs deoxygenated blood flow) ratio of 0.72 associated with a 16% decrease of left ventricular (LV) output and 18% increase of RV output. The flow from LV to aortic arch branches increases by 16%, while LV contribution to the lower body flow decreases by 29%. Total flow in the descending aorta (DAo) increases by 18% due to RV contribution through the PS with flow into the distal PA branches decreasing. PS induces 18% increase of RV work due to its larger stroke volume pumped against lower afterload. Nonetheless, larger RV work does not lead to increased RV end-diastolic volume. Three-dimensional flow assessment demonstrates the PS jet impinging with a high velocity and wall shear stress on the opposite DAo wall with the most of the shunt flow being diverted to the DAo. Increasing the PS diameter from 5mm up to 10mm results in a nearly linear increase in post-operative shunt flow and a nearly linear decrease in shunt pressure-drop. In conclusion, this model reasonably represents patient-specific haemodynamics pre- and post-creation of the PS, providing insights into physiology of this complex condition, and presents a predictive tool that could be useful for clinical decision-making regarding suitability for PS in PAH patients with drug-resistant suprasystemic PAH.

A method of parameter estimation for cardiovascular hemodynamics based on deep learning and its application to personalize a reduced-order model

Precise model personalization is a key step towards the application of cardiovascular physical models. In this manuscript, we propose to use deep learning (DL) to solve the parameter estimation problem in cardiovascular hemodynamics. Based on the convolutional neural network (CNN) and fully connected neural network (FCNN), a multi-input deep neural network (DNN) model is developed to map the nonlinear relationship between measurements and the parameters to be estimated. In this model, two separate network structures are designed to extract the features of two types of measurement data, including pressure waveforms and a vector composed of heart rate (HR) and pulse transit time (PTT), and a shared structure is used to extract their combined dependencies on the parameters. Besides, we try to use the transfer learning (TL) technology to further strengthen the personalized characteristics of a trained-well network. For assessing the proposed method, we conducted the parameter estimation using synthetic data and in vitro data respectively, and in the test with synthetic data, we evaluated the performance of the TL algorithm through two individuals with different characteristics. A series of estimation results show that the estimated parameters are in good agreement with the true values. Furthermore, it is also found that the estimation accuracy can be significantly improved by a multicycle combination strategy. Therefore, we think that the proposed method has the potential to be used for parameter estimation in cardiovascular hemodynamics, which can provide an immediate, accurate, and sustainable personalization process, and deserves more attention in the future.

Parameter estimation for closed-loop lumped parameter models of the systemic circulation using synthetic data

Physics-based models can be applied to describe mechanisms in both health and disease, which has the potential to accelerate the development of personalized medicine. The aim of this study was to investigate the feasibility of personalizing a model of systemic hemodynamics by estimating model parameters. We investigated the feasibility of estimating model parameters for a closed-loop lumped parameter model of the left heart and systemic circulation using the step-wise subset reduction method. This proceeded by first investigating the structural identifiability of the model parameters. Secondly we performed sensitivity analysis to determine which parameters were most influential on the most relevant model outputs. Finally, we constructed a sequence of progressively smaller subsets including parameters based on their ranking by model output influence. The model was then optimized to data for each set of parameters to evaluate how well the parameters could be estimated for each subset. The subsequent results allowed assessment of how different data sets, and noise affected the parameter estimates. In the noiseless case, all parameters could be calibrated to less than 10^-3% error using time series data, while errors using clinical index data could reach over 100%. With 5% normally distributed noise the accuracy was limited to be within 10% error for the five most sensitive parameters, while the four least sensitive parameters were unreliably estimated for waveform data. The three least sensitive parameters were particularly challenging to estimate so these should be prioritised for measurement. Cost functions based on time series such as pressure waveforms, were found to give better parameter estimates than cost functions based on standard indices used in clinical assessment of the cardiovascular system, for example stroke volume (SV) and pulse pressure (PP). Averaged parameter estimate errors were reduced by several orders of magnitude by choosing waveforms for noiseless synthetic data. Also when measurement data were noisy, the parameter estimation procedure based on continuous waveforms was more accurate than that based on clinical indices. By application of the step-wise subset reduction method we demonstrated that by the addition of venous pressure to the cost function, or conversely fixing the systemic venous compliance parameter at an accurate value improved all parameter estimates, especially the diastolic filling parameters which have least influence on the aortic pressure.

A physiologically realistic virtual patient database for the study of arterial haemodynamics

This study creates a physiologically realistic virtual patient database (VPD), representing the human arterial system, for the primary purpose of studying the effects of arterial disease on haemodynamics. A low dimensional representation of an anatomically detailed arterial network is outlined, and a physiologically realistic posterior distribution for its parameters constructed through the use of a Bayesian approach. This approach combines both physiological/geometrical constraints and the available measurements reported in the literature. A key contribution of this work is to present a framework for including all such available information for the creation of virtual patients (VPs). The Markov Chain Monte Carlo (MCMC) method is used to sample random VPs from this posterior distribution, and the pressure and flow-rate profiles associated with each VP computed through a physics based model of pulse wave propagation. This combination of the arterial network parameters (representing a virtual patient) and the haemodynamics waveforms of pressure and flow-rates at various locations (representing functional response and potential measurements that can be acquired in the virtual patient) makes up the VPD. While 75,000 VPs are sampled from the posterior distribution, 10,000 are discarded as the initial burn-in period of the MCMC sampler. A further 12,857 VPs are subsequently removed due to the presence of negative average flow-rate, reducing the VPD to 52,143. Due to undesirable behaviour observed in some VPs-asymmetric under- and over-damped pressure and flow-rate profiles in left and right sides of the arterial system-a filter is proposed to remove VPs showing such behaviour. Post application of the filter, the VPD has 28,868 subjects. It is shown that the methodology is appropriate by comparing the VPD statistics to those reported in literature across real populations. Generally, a good agreement between the two is found while respecting physiological/geometrical constraints.

Predictive Modeling of Secondary Pulmonary Hypertension in Left Ventricular Diastolic Dysfunction

Diastolic dysfunction is a common pathology occurring in about one third of patients affected by heart failure. This condition may not be associated with a marked decrease in cardiac output or systemic pressure and therefore is more difficult to diagnose than its systolic counterpart. Compromised relaxation or increased stiffness of the left ventricle induces an increase in the upstream pulmonary pressures, and is classified as secondary or group II pulmonary hypertension (2018 Nice classification). This may result in an increase in the right ventricular afterload leading to right ventricular failure. Elevated pulmonary pressures are therefore an important clinical indicator of diastolic heart failure (sometimes referred to as heart failure with preserved ejection fraction, HFpEF), showing significant correlation with associated mortality. However, accurate measurements of this quantity are typically obtained through invasive catheterization and after the onset of symptoms. In this study, we use the hemodynamic consistency of a differential-algebraic circulation model to predict pulmonary pressures in adult patients from other, possibly non-invasive, clinical data. We investigate several aspects of the problem, including the ability of model outputs to represent a sufficiently wide pathologic spectrum, the identifiability of the model's parameters, and the accuracy of the predicted pulmonary pressures. We also find that a classifier using the assimilated model parameters as features is free from the problem of missing data and is able to detect pulmonary hypertension with sufficiently high accuracy. For a cohort of 82 patients suffering from various degrees of heart failure severity, we show that systolic, diastolic, and wedge pulmonary pressures can be estimated on average within 8, 6, and 6 mmHg, respectively. We also show that, in general, increased data availability leads to improved predictions.

Implementation and Calibration of a Deep Neural Network to Predict Parameters of Left Ventricular Systolic Function Based on Pulmonary and Systemic Arterial Pressure Signals

The evaluation of cardiac contractility by the assessment of the ventricular systolic elastance function is clinically challenging and cannot be easily obtained at the bedside. In this work, we present a framework characterizing left ventricular systolic function from clinically readily available data, including systemic and pulmonary arterial pressure signals. We implemented and calibrated a deep neural network (DNN) consisting of a multi-layer perceptron with 4 fully connected hidden layers and with 16 neurons per layer, which was trained with data obtained from a lumped model of the cardiovascular system modeling different levels of cardiac function. The lumped model included a function of circulatory autoregulation from carotid baroreceptors in pulsatile conditions. Inputs for the DNN were systemic and pulmonary arterial pressure curves. Outputs from the DNN were parameters of the lumped model characterizing left ventricular systolic function, especially end-systolic elastance. The DNN adequately performed and accurately recovered the relevant hemodynamic parameters with a mean relative error of less than 2%. Therefore, our framework can easily provide complex physiological parameters of cardiac contractility, which could lead to the development of invaluable tools for the clinical evaluation of patients with severe cardiac dysfunction.

Uncertainty in model-based treatment decision support: Applied to aortic valve stenosis

Patient outcome in trans-aortic valve implantation (TAVI) therapy partly relies on a patient's haemodynamic properties that cannot be determined from current diagnostic methods alone. In this study, we predict changes in haemodynamic parameters (as a part of patient outcome) after valve replacement treatment in aortic stenosis patients. A framework to incorporate uncertainty in patient-specific model predictions for decision support is presented. A 0D lumped parameter model including the left ventricle, a stenotic valve and systemic circulatory system has been developed, based on models published earlier. The unscented Kalman filter (UKF) is used to optimize model input parameters to fit measured data pre-intervention. After optimization, the valve treatment is simulated by significantly reducing valve resistance. Uncertain model parameters are then propagated using a polynomial chaos expansion approach. To test the proposed framework, three in silico test cases are developed with clinically feasible measurements. Quality and availability of simulated measured patient data are decreased in each case. The UKF approach is compared to a Monte Carlo Markov Chain (MCMC) approach, a well-known approach in modelling predictions with uncertainty. Both methods show increased confidence intervals as measurement quality decreases. By considering three in silico test-cases we were able to show that the proposed framework is able to incorporate optimization uncertainty in model predictions and is faster and the MCMC approach, although it is more sensitive to noise in flow measurements. To conclude, this work shows that the proposed framework is ready to be applied to real patient data.

Cardiovascular models for personalised medicine: Where now and where next?

The aim of this position paper is to provide a brief overview of the current status of cardiovascular modelling and of the processes required and some of the challenges to be addressed to see wider exploitation in both personal health management and clinical practice. In most branches of engineering the concept of the digital twin, informed by extensive and continuous monitoring and coupled with robust data assimilation and simulation techniques, is gaining traction: the Gartner Group listed it as one of the top ten digital trends in 2018. The cardiovascular modelling community is starting to develop a much more systematic approach to the combination of physics, mathematics, control theory, artificial intelligence, machine learning, computer science and advanced engineering methodology, as well as working more closely with the clinical community to better understand and exploit physiological measurements, and indeed to develop jointly better measurement protocols informed by model-based understanding. Developments in physiological modelling, model personalisation, model outcome uncertainty, and the role of models in clinical decision support are addressed and 'where-next' steps and challenges discussed.

Data-Augmented Modeling of Intracranial Pressure

Precise management of patients with cerebral diseases often requires intracranial pressure (ICP) monitoring, which is highly invasive and requires a specialized ICU setting. The ability to noninvasively estimate ICP is highly compelling as an alternative to, or screening for, invasive ICP measurement. Most existing approaches for noninvasive ICP estimation aim to build a regression function that maps noninvasive measurements to an ICP estimate using statistical learning techniques. These data-based approaches have met limited success, likely because the amount of training data needed is onerous for this complex applications. In this work, we discuss an alternative strategy that aims to better utilize noninvasive measurement data by leveraging mechanistic understanding of physiology. Specifically, we developed a Bayesian framework that combines a multiscale model of intracranial physiology with noninvasive measurements of cerebral blood flow using transcranial Doppler. Virtual experiments with synthetic data are conducted to verify and analyze the proposed framework. A preliminary clinical application study on two patients is also performed in which we demonstrate the ability of this method to improve ICP prediction.

Virtual medicine: Utilization of the advanced cardiac imaging patient avatar for procedural planning and facilitation

Advances in imaging technology have led to a paradigm shift from planning of cardiovascular procedures and surgeries requiring the actual patient in a "brick and mortar" hospital to utilization of the digitalized patient in the virtual hospital. Cardiovascular computed tomographic angiography (CCTA) and cardiovascular magnetic resonance (CMR) digitalized 3-D patient representation of individual patient anatomy and physiology serves as an avatar allowing for virtual delineation of the most optimal approaches to cardiovascular procedures and surgeries prior to actual hospitalization. Pre-hospitalization reconstruction and analysis of anatomy and pathophysiology previously only accessible during the actual procedure could potentially limit the intrinsic risks related to time in the operating room, cardiac procedural laboratory and overall hospital environment. Although applications are specific to areas of cardiovascular specialty focus, there are unifying themes related to the utilization of technologies. The virtual patient avatar computer can also be used for procedural planning, computational modeling of anatomy, simulation of predicted therapeutic result, printing of 3-D models, and augmentation of real time procedural performance. Examples of the above techniques are at various stages of development for application to the spectrum of cardiovascular disease processes, including percutaneous, surgical and hybrid minimally invasive interventions. A multidisciplinary approach within medicine and engineering is necessary for creation of robust algorithms for maximal utilization of the virtual patient avatar in the digital medical center. Utilization of the virtual advanced cardiac imaging patient avatar will play an important role in the virtual health care system. Although there has been a rapid proliferation of early data, advanced imaging applications require further assessment and validation of accuracy, reproducibility, standardization, safety, efficacy, quality, cost effectiveness, and overall value to medical care.

Inverse problems in reduced order models of cardiovascular haemodynamics: aspects of data assimilation and heart rate variability

Inverse problems in cardiovascular modelling have become increasingly important to assess each patient individually. These problems entail estimation of patient-specific model parameters from uncertain measurements acquired in the clinic. In recent years, the method of data assimilation, especially the unscented Kalman filter, has gained popularity to address computational efficiency and uncertainty consideration in such problems. This work highlights and presents solutions to several challenges of this method pertinent to models of cardiovascular haemodynamics. These include methods to (i) avoid ill-conditioning of the covariance matrix, (ii) handle a variety of measurement types, (iii) include a variety of prior knowledge in the method, and (iv) incorporate measurements acquired at different heart rates, a common situation in the clinic where the patient state differs according to the clinical situation. Results are presented for two patient-specific cases of congenital heart disease. To illustrate and validate data assimilation with measurements at different heart rates, the results are presented on a synthetic dataset and on a patient-specific case with heart valve regurgitation. It is shown that the new method significantly improves the agreement between model predictions and measurements. The developed methods can be readily applied to other pathophysiologies and extended to dynamical systems which exhibit different responses under different sets of known parameters or different sets of inputs (such as forcing/excitation frequencies).

Bridging the gap between measurements and modelling: a cardiovascular functional avatar

Lumped parameter models of the cardiovascular system have the potential to assist researchers and clinicians to better understand cardiovascular function. The value of such models increases when they are subject specific. However, most approaches to personalize lumped parameter models have thus far required invasive measurements or fall short of being subject specific due to a lack of the necessary clinical data. Here, we propose an approach to personalize parameters in a model of the heart and the systemic circulation using exclusively non-invasive measurements. The personalized model is created using flow data from four-dimensional magnetic resonance imaging and cuff pressure measurements in the brachial artery. We term this personalized model the cardiovascular avatar. In our proof-of-concept study, we evaluated the capability of the avatar to reproduce pressures and flows in a group of eight healthy subjects. Both quantitatively and qualitatively, the model-based results agreed well with the pressure and flow measurements obtained in vivo for each subject. This non-invasive and personalized approach can synthesize medical data into clinically relevant indicators of cardiovascular function, and estimate hemodynamic variables that cannot be assessed directly from clinical measurements.