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

Building Digital Twins for Cardiovascular Health: From Principles to Clinical Impact

The past several decades have seen rapid advances in diagnosis and treatment of cardiovascular diseases and stroke, enabled by technological breakthroughs in imaging, genomics, and physiological monitoring, coupled with therapeutic interventions. We now face the challenge of how to (1) rapidly process large, complex multimodal and multiscale medical measurements; (2) map all available data streams to the trajectories of disease states over the patient's lifetime; and (3) apply this information for optimal clinical interventions and outcomes. Here we review new advances that may address these challenges using digital twin technology to fulfill the promise of personalized cardiovascular medical practice. Rooted in engineering mechanics and manufacturing, the digital twin is a virtual representation engineered to model and simulate its physical counterpart. Recent breakthroughs in scientific computation, artificial intelligence, and sensor technology have enabled rapid bidirectional interactions between the virtual-physical counterparts with measurements of the physical twin that inform and improve its virtual twin, which in turn provide updated virtual projections of disease trajectories and anticipated clinical outcomes. Verification, validation, and uncertainty quantification builds confidence and trust by clinicians and patients in the digital twin and establishes boundaries for the use of simulations in cardiovascular medicine. Mechanistic physiological models form the fundamental building blocks of the personalized digital twin that continuously forecast optimal management of cardiovascular health using individualized data streams. We present exemplars from the existing body of literature pertaining to mechanistic model development for cardiovascular dynamics and summarize existing technical challenges and opportunities pertaining to the foundation of a digital twin.

Guidelines for mechanistic modeling and analysis in cardiovascular research

Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.

Numerical study of hemodynamic flow in the aortic vessel of Williams syndrome patient with congenital heart disease

Congenital arterial stenosis such as supravalvar aortic stenosis (SVAS) are highly prevalent in Williams syndrome (WS) and other arteriopathies pose a substantial health risk. Conventional tools for severity assessment, including clinical findings and pressure gradient estimations, often fall short due to their susceptibility to transient physiological changes and disease stage influences. Moreover, in the pediatric population, the severity of these and other congenital heart defects (CHDs) often restricts the applicability of invasive techniques for obtaining crucial physiological data. Conversely, evaluating CHDs and their progression requires a comprehensive understanding of intracardiac blood flow. Current imaging modalities, such as blood speckle imaging (BSI) and four-dimensional magnetic resonance imaging (4D MRI) face limitations in resolving flow data, especially in cases of elevated flow velocities. To address these challenges, we devised a computational framework employing zero-dimensional (0D) lumped parameter models coupled with patient-specific reconstructed geometries pre- and post-surgical intervention to execute computational fluid dynamic (CFD) simulations. This framework facilitates the analysis and visualization of intricate blood flow patterns, offering insights into geometry and flow dynamics alterations impacting cardiac function. In this study, we aim to assess the efficacy of surgical intervention in correcting an extreme aortic defect in a patient with WS, leading to reductions in wall shear stress (WSS), maximum velocity magnitude, pressure drop, and ultimately a decrease in cardiac workload.

Stress Load and Ascending Aortic Aneurysms: An Observational, Longitudinal, Single-Center Study Using Computational Fluid Dynamics

Ascending aortic aneurysm (AAoA) is a silent disease with high mortality; however, the factors associated with a worse prognosis are not completely understood. The objective of this observational, longitudinal, single-center study was to identify the hemodynamic patterns and their influence on AAoA growth using computational fluid dynamics (CFD), focusing on the effects of geometrical variations on aortic hemodynamics. Personalized anatomic models were obtained from angiotomography scans of 30 patients in two different years (with intervals of one to three years between them), of which 16 (53%) showed aneurysm growth (defined as an increase in the ascending aorta volume by 5% or more). Numerically determined velocity and pressure fields were compared with the outcome of aneurysm growth. Through a statistical analysis, hemodynamic characteristics were found to be associated with aneurysm growth: average and maximum high pressure (superior to 100 Pa); average and maximum high wall shear stress (superior to 7 Pa) combined with high pressure (>100 Pa); and stress load over time (maximum pressure multiplied by the time interval between the exams). This study provides insights into a worse prognosis of this serious disease and may collaborate for the expansion of knowledge about mechanobiology in the progression of AAoA.

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.

Haemodynamic effects of hypertension and type 2 diabetes: Insights from a 4D flow MRI-based personalized cardiovascular mathematical model

Type 2 diabetes (T2D) and hypertension increase the risk of cardiovascular diseases mediated by whole-body changes to metabolism, cardiovascular structure and haemodynamics. The haemodynamic changes related to hypertension and T2D are complex and subject-specific, however, and not fully understood. We aimed to investigate the haemodynamic mechanisms in T2D and hypertension by comparing the haemodynamics between healthy controls and subjects with T2D, hypertension, or both. For all subjects, we combined 4D flow magnetic resonance imaging data, brachial blood pressure and a cardiovascular mathematical model to create a comprehensive subject-specific analysis of central haemodynamics. When comparing the subject-specific haemodynamic parameters between the four groups, the predominant haemodynamic difference is impaired left ventricular relaxation in subjects with both T2D and hypertension compared to subjects with only T2D, only hypertension and controls. The impaired relaxation indicates that, in this cohort, the long-term changes in haemodynamic load of co-existing T2D and hypertension cause diastolic dysfunction demonstrable at rest, whereas either disease on its own does not. However, through subject-specific predictions of impaired relaxation, we show that altered relaxation alone is not enough to explain the subject-specific and group-related differences; instead, a combination of parameters is affected in T2D and hypertension. These results confirm previous studies that reported more adverse effects from the combination of T2D and hypertension compared to either disease on its own. Furthermore, this shows the potential of personalized cardiovascular models in providing haemodynamic mechanistic insights and subject-specific predictions that could aid in the understanding and treatment planning of patients with T2D and hypertension. KEY POINTS: The combination of 4D flow magnetic resonance imaging data and a cardiovascular mathematical model allows for a comprehensive analysis of subject-specific haemodynamic parameters that otherwise cannot be derived non-invasively. Using this combination, we show that diastolic dysfunction in subjects with both type 2 diabetes (T2D) and hypertension is the main group-level difference between controls, subjects with T2D, subjects with hypertension, and subjects with both T2D and hypertension. These results suggest that, in this relatively healthy population, the additional load of both hypertension and T2D affects the haemodynamic function of the left ventricle, whereas each disease on its own is not enough to cause significant effects under resting conditions. Finally, using the subject-specific model, we show that the haemodynamic effects of diastolic dysfunction alone are not sufficient to explain all the observed haemodynamic differences. Instead, additional subject-specific variations in cardiac and vascular function combine to explain the complex haemodynamics of subjects affected by hypertension and/or T2D.

Evaluation of Different Cannulation Strategies for Aortic Arch Surgery Using a Cardiovascular Numerical Simulator

Aortic disease has a significant impact on quality of life. The involvement of the aortic arch requires the preservation of blood supply to the brain during surgery. Deep hypothermic circulatory arrest is an established technique for this purpose, although neurological injury remains high. Additional techniques have been used to reduce risk, although controversy still remains. A three-way cannulation approach, including both carotid arteries and the femoral artery or the ascending aorta, has been used successfully for aortic arch replacement and redo procedures. We developed circuits of the circulation to simulate blood flow during this type of cannulation set up. The CARDIOSIM© cardiovascular simulation platform was used to analyse the effect on haemodynamic and energetic parameters and the benefit derived in terms of organ perfusion pressure and flow. Our simulation approach based on lumped-parameter modelling, pressure-volume analysis and modified time-varying elastance provides a theoretical background to a three-way cannulation strategy for aortic arch surgery with correlation to the observed clinical practice.

Structure (Epicardial Stenosis) and Function (Microvascular Dysfunction) That Influence Coronary Fractional Flow Reserve Estimation

Background. The treatment of coronary stenosis is decided by performing high risk invasive surgery to generate the fractional flow reserve diagnostics index, a ratio of distal to proximal pressures in respect of coronary atherosclerotic plaques. Non-invasive methods are a need of the times that necessitate the use of mathematical models of coronary hemodynamic physiology. This study proposes an extensible mathematical description of the coronary vasculature that provides an estimate of coronary fractional flow reserve. Methods. By adapting an existing computational model of human coronary blood flow, the effects of large vessel stenosis and microvascular disease on fractional flow reserve were quantified. Several simulations generated flow and pressure information, which was used to compute fractional flow reserve under several conditions including focal stenosis, diffuse stenosis, and microvascular disease. Sensitivity analysis was used to uncover the influence of model parameters on fractional flow reserve. The model was simulated as coupled non-linear ordinary differential equations and numerically solved using our implicit higher order method. Results. Large vessel stenosis affected fractional flow reserve. The model predicts that the presence, rather than severity, of microvascular disease affects coronary flow deleteriously. Conclusions. The model provides a computationally inexpensive instrument for future in silico coronary blood flow investigations as well as clinical-imaging decision making. A combination of focal and diffuse stenosis appears to be essential to limit coronary flow. In addition to pressure measurements in the large epicardial vessels, diagnosis of microvascular disease is essential. The independence of the index with respect to heart rate suggests that computationally inexpensive steady state simulations may provide sufficient information to reliably compute the index.

Personalized intervention cardiology with transcatheter aortic valve replacement made possible with a non-invasive monitoring and diagnostic framework

One of the most common acute and chronic cardiovascular disease conditions is aortic stenosis, a disease in which the aortic valve is damaged and can no longer function properly. Moreover, aortic stenosis commonly exists in combination with other conditions causing so many patients suffer from the most general and fundamentally challenging condition: complex valvular, ventricular and vascular disease (C3VD). Transcatheter aortic valve replacement (TAVR) is a new less invasive intervention and is a growing alternative for patients with aortic stenosis. Although blood flow quantification is critical for accurate and early diagnosis of C3VD in both pre and post-TAVR, proper diagnostic methods are still lacking because the fluid-dynamics methods that can be used as engines of new diagnostic tools are not well developed yet. Despite remarkable advances in medical imaging, imaging on its own is not enough to quantify the blood flow effectively. Moreover, understanding of C3VD in both pre and post-TAVR and its progression has been hindered by the absence of a proper non-invasive tool for the assessment of the cardiovascular function. To enable the development of new non-invasive diagnostic methods, we developed an innovative image-based patient-specific computational fluid dynamics framework for patients with C3VD who undergo TAVR to quantify metrics of: (1) global circulatory function; (2) global cardiac function as well as (3) local cardiac fluid dynamics. This framework is based on an innovative non-invasive Doppler-based patient-specific lumped-parameter algorithm and a 3-D strongly-coupled fluid-solid interaction. We validated the framework against clinical cardiac catheterization and Doppler echocardiographic measurements and demonstrated its diagnostic utility by providing novel analyses and interpretations of clinical data in eleven C3VD patients in pre and post-TAVR status. Our findings position this framework as a promising new non-invasive diagnostic tool that can provide blood flow metrics while posing no risk to the patient. The diagnostic information, that the framework can provide, is vitally needed to improve clinical outcomes, to assess patient risk and to plan treatment.

An integrated lumped-parameter model of the cardiovascular system for the simulation of acute ischemic stroke: description of instantaneous changes in hemodynamics

Acute Ischemic Stroke (AIS) is defined as the acute condition of occlusion of a cerebral artery and is often caused by a Hypertensive Condition (HC). Due to its sudden occurrence, AIS is not observable the right moment it occurs, thus information about instantaneous changes in hemodynamics is limited. This study aimed to propose an integrated Lumped Parameter (LP) model of the cardiovascular system to simulate an AIS and describe instantaneous changes in hemodynamics. In the integrated LP model of the cardiovascular system, heart chambers have been modelled with elastance systems with controlled pressure inputs; heart valves have been modelled with static open/closed pressure-controlled valves; eventually, the vasculature has been modelled with resistor-inductor-capacitor (RLC) direct circuits and have been linked to the rest of the system through a series connection. After simulating physiological conditions, HC has been simulated by changing pressure inputs and constant RLC parameters. Then, AIS occurring in arteries of different sizes have been simulated by considering time-dependent RLC parameters due to the elimination from the model of the occluding artery; instantaneous changes in hemodynamics have been evaluated by Systemic Arteriolar Flow (Qa) and Systemic Arteriolar Pressure (Pa) drop with respect to those measured in HC. Occlusion of arteries of different sizes leaded to an average Qa drop of 0.38 ml/s per cardiac cycle (with minimum and maximum values of 0.04 ml/s and 1.93 ml/s) and average Pa drop of 0.39 mmHg, (with minimum and maximum values of 0.04 mmHg and 1.98 mmHg). In conclusion, hemodynamic variations due to AIS are very small with respect to HC. A direct relation between the inverse of the length of the artery in which the occlusion occurs and the hemodynamic variations has been highlighted; this may allow to link the severity of AIS to the length of the interested artery.

Hybrid modelling for stroke care: Review and suggestions of new approaches for risk assessment and simulation of scenarios

Stroke is an example of a complex and multi-factorial disease involving multiple organs, timescales, and disease mechanisms. To deal with this complexity, and to realize Precision Medicine of stroke, mathematical models are needed. Such approaches include: 1) machine learning, 2) bioinformatic network models, and 3) mechanistic models. Since these three approaches have complementary strengths and weaknesses, a hybrid modelling approach combining them would be the most beneficial. However, no concrete approach ready to be implemented for a specific disease has been presented to date. In this paper, we both review the strengths and weaknesses of the three approaches, and propose a roadmap for hybrid modelling in the case of stroke care. We focus on two main tasks needed for the clinical setting: a) For stroke risk calculation, we propose a new two-step approach, where non-linear mixed effects models and bioinformatic network models yield biomarkers which are used as input to a machine learning model and b) For simulation of care scenarios, we propose a new four-step approach, which revolves around iterations between simulations of the mechanistic models and imputations of non-modelled or non-measured variables. We illustrate and discuss the different approaches in the context of Precision Medicine for stroke.

Impact of prosthetic mitral valve orientation on the ventricular flow field: Comparison using patient-specific computational fluid dynamics

Significant mitral valve regurgitation creates progressive adverse remodeling of the left ventricle (LV). Replacement of the failing valve with a prosthesis generally improves patient outcomes but leaves the patient with non-physiological intracardiac flow patterns that might contribute to their future risk of thrombus formation and embolism. It has been suggested that the angular orientation of the implanted valve might modify the postoperative distortion of the intraventricular flow field. In this study, we investigated the effect of prosthetic valve orientation on LV flow patterns by using heart geometry from a patient with LV dysfunction and a competent native mitral valve to calculate intracardiac flow fields with computational fluid dynamics (CFD). Results were validated using in vivo 4D Flow MRI. The computed flow fields were compared to calculations following virtual implantation of a mechanical heart valve oriented in four different angles to assess the effect of leaflet position. Flow patterns were visualized in long- and short-axes and quantified with flow component analysis. In comparison to a native valve, valve implantation increased the proportion of the mitral inflow remaining in the basal region and further increased the residual volume in the apical area. Only slight changes due to valve orientation were observed. Using our numerical framework, we demonstrated quantitative changes in left ventricular blood flow due to prosthetic mitral replacement. This framework may be used to improve design of prosthetic heart valves and implantation procedures to minimize the potential for apical flow stasis, and potentially assist personalized treatment planning.

From CT to artificial intelligence for complex assessment of plaque-associated risk

The recent technological developments in the field of cardiac imaging have established coronary computed tomography angiography (CCTA) as a first-line diagnostic tool in patients with suspected coronary artery disease (CAD). CCTA offers robust information on the overall coronary circulation and luminal stenosis, also providing the ability to assess the composition, morphology, and vulnerability of atherosclerotic plaques. In addition, the perivascular adipose tissue (PVAT) has recently emerged as a marker of increased cardiovascular risk. The addition of PVAT quantification to standard CCTA imaging may provide the ability to extract information on local inflammation, for an individualized approach in coronary risk stratification. The development of image post-processing tools over the past several years allowed CCTA to provide a significant amount of data that can be incorporated into machine learning (ML) applications. ML algorithms that use radiomic features extracted from CCTA are still at an early stage. However, the recent development of artificial intelligence will probably bring major changes in the way we integrate clinical, biological, and imaging information, for a complex risk stratification and individualized therapeutic decision making in patients with CAD. This review aims to present the current evidence on the complex role of CCTA in the detection and quantification of vulnerable plaques and the associated coronary inflammation, also describing the most recent developments in the radiomics-based machine learning approach for complex assessment of plaque-associated risk.

Monitoring of cardiovascular physiology augmented by a patient-specific biomechanical model during general anesthesia. A proof of concept study

During general anesthesia (GA), direct analysis of arterial pressure or aortic flow waveforms may be inconclusive in complex situations. Patient-specific biomechanical models, based on data obtained during GA and capable to perform fast simulations of cardiac cycles, have the potential to augment hemodynamic monitoring. Such models allow to simulate Pressure-Volume (PV) loops and estimate functional indicators of cardiovascular (CV) system, e.g. ventricular-arterial coupling (Vva), cardiac efficiency (CE) or myocardial contractility, evolving throughout GA. In this prospective observational study, we created patient-specific biomechanical models of heart and vasculature of a reduced geometric complexity for n = 45 patients undergoing GA, while using transthoracic echocardiography and aortic pressure and flow signals acquired in the beginning of GA (baseline condition). If intraoperative hypotension (IOH) appeared, diluted norepinephrine (NOR) was administered and the model readjusted according to the measured aortic pressure and flow signals. Such patients were a posteriori assigned into a so-called hypotensive group. The accuracy of simulated mean aortic pressure (MAP) and stroke volume (SV) at baseline were in accordance with the guidelines for the validation of new devices or reference measurement methods in all patients. After NOR administration in the hypotensive group, the percentage of concordance with 10% exclusion zone between measurement and simulation was >95% for both MAP and SV. The modeling results showed a decreased Vva (0.64±0.37 vs 0.88±0.43; p = 0.039) and an increased CE (0.8±0.1 vs 0.73±0.11; p = 0.042) in hypotensive vs normotensive patients. Furthermore, Vva increased by 92±101%, CE decreased by 13±11% (p < 0.001 for both) and contractility increased by 14±11% (p = 0.002) in the hypotensive group post-NOR administration. In this work we demonstrated the application of fast-running patient-specific biophysical models to estimate PV loops and functional indicators of CV system using clinical data available during GA. The work paves the way for model-augmented hemodynamic monitoring at operating theatres or intensive care units to enhance the information on patient-specific physiology.

A diagnostic, monitoring, and predictive tool for patients with complex valvular, vascular and ventricular diseases

Hemodynamics quantification is critically useful for accurate and early diagnosis, but we still lack proper diagnosticmethods for many cardiovascular diseases. Furthermore, as most interventions intend to recover the healthy condition, the ability to monitor and predict hemodynamics following interventions can have significant impacts on saving lives. Predictive methods are rare, enabling prediction of effects of interventions, allowing timely and personalized interventions and helping critical clinical decision making about life-threatening risks based on quantitative data. In this study, an innovative non-invasive imaged-based patient-specific diagnostic, monitoring and predictive tool (called C3VI-CMF) was developed, enabling quantifying (1) details of physiological flow and pressures through the heart and circulatory system; (2) heart function metrics. C3VI-CMF also predicts the breakdown of the effects of each disease constituents on the heart function. Presently, neither of these can be obtained noninvasively in patients and when invasive procedures are undertaken, the collected metrics cannot be by any means as complete as the ones C3VI-CMF provides. C3VI-CMF purposefully uses a limited number of noninvasive input parameters all of which can be measured using Doppler echocardiography and sphygmomanometer. Validation of C3VI-CMF, against cardiac catheterization in forty-nine patients with complex cardiovascular diseases, showed very good agreement with the measurements.

Dynamic Modulation of Device-Arterial Coupling to Determine Cardiac Output and Vascular Resistance

Clinical adoption of mechanical circulatory support for shock is rapidly expanding. Achieving optimal therapeutic benefit requires metrics of state to guide titration and weaning of support. Using the transvalvular positioning of a percutaneous ventricular assist device (pVAD), device:heart interactions are leveraged to determine cardiac output (CO) and systemic vascular resistance (SVR) near-continuously without disrupting therapeutic function. An automated algorithm rapidly alternates between device support levels to dynamically modulate physiological response. Employing a two-element lumped parameter model of the vasculature, SVR and CO are quantified directly from measurements obtained by the pVAD without external calibration or invasive catheters. The approach was validated in an acute porcine model across a range of cardiac (CO = 3-10.6 L/min) and vascular (SVR = 501-1897 dyn s/cm⁵) states. Cardiac output calculations closely correlated (r = 0.82) to measurements obtained by the pulmonary artery catheter-based thermodilution method with a mean bias of 0.109 L/min and limits of agreement from - 1.67 to 1.89 L/min. SVR was also closely correlated (r = 0.86) to traditional catheter-based measurements with a mean bias of 62.1 dyn s/cm⁵ and limits of agreement from - 260 to 384 dyn s/cm⁵. Use of diagnostics integrated into therapeutic device function enables the potential for optimizing support to improve outcomes for cardiogenic shock.

Synergistic Model of Cardiac Function with a Heart Assist Device

The breakdown of cardiac self-organization leads to heart diseases and failure, the number one cause of death worldwide. The left ventricular pressure-volume relation plays a key role in the diagnosis and treatment of heart diseases. Lumped-parameter models combined with pressure-volume loop analysis are very effective in simulating clinical scenarios with a view to treatment optimization and outcome prediction. Unfortunately, often invoked in this analysis is the traditional, time-varying elastance concept, in which the ratio of the ventricular pressure to its volume is prescribed by a periodic function of time, instead of being calculated consistently according to the change in feedback mechanisms (e.g., the lack or breakdown of self-organization) in heart diseases. Therefore, the application of the time-varying elastance for the analysis of left ventricular assist device (LVAD)-heart interactions has been questioned. We propose a paradigm shift from the time-varying elastance concept to a synergistic model of cardiac function by integrating the mechanical, electric, and chemical activity on microscale sarcomere and macroscale heart levels and investigating the effect of an axial rotary pump on a failing heart. We show that our synergistic model works better than the time-varying elastance model in reproducing LVAD-heart interactions with sufficient accuracy to describe the left ventricular pressure-volume relation.

Leveraging Device-Arterial Coupling to Determine Cardiac and Vascular State

OBJECTIVE

Limitations in available diagnostic metrics restrict the efficacy of managing therapies for cardiogenic shock. In current clinical practice, cardiovascular state is inferred through measurement of pulmonary capillary wedge pressure and reliance on linear approximations between pressure and flow to estimate peripheral vascular resistance. Mechanical circulatory support devices residing within the left ventricle and aorta provide an opportunity for both determining cardiac and vascular state and offering therapeutic benefit. We leverage the controllable mode of operation and transvalvular position of an indwelling percutaneous ventricular assist device to assess vascular and, in turn, cardiac state through the effects of device-arterial coupling across different levels of device support.

METHODS

Vascular state is determined by measuring changes in the pressure waveforms induced through intentional variation in the device generated blood flow. We evaluate this impact by applying a lumped parameter model to quantify state-specific vascular resistance and compliance and calculate beat-to-beat stroke volume and cardiac output in both animal models and retrospective patient data without external calibration.

RESULTS

Vascular state was accurately predicted in patients and animals in both baseline and experimental conditions. In the animal, stroke volume was predicted within a total root mean square error of 3.71 mL (n = 482).

CONCLUSION

We demonstrate that device-arterial coupling is a powerful tool for evaluating patient and state specific parameters of cardiovascular function.

SIGNIFICANCE

These insights may yield improved clinical care and support the development of next generation mechanical circulatory support devices that determine and operate in tandem with the supported organ.

Pulse transit time estimation of aortic pulse wave velocity and blood pressure using machine learning and simulated training data

Recent developments in cardiovascular modelling allow us to simulate blood flow in an entire human body. Such model can also be used to create databases of virtual subjects, with sizes limited only by computational resources. In this work, we study if it is possible to estimate cardiovascular health indices using machine learning approaches. In particular, we carry out theoretical assessment of estimating aortic pulse wave velocity, diastolic and systolic blood pressure and stroke volume using pulse transit/arrival timings derived from photopletyshmography signals. For predictions, we train Gaussian process regression using a database of virtual subjects generated with a cardiovascular simulator. Simulated results provides theoretical assessment of accuracy for predictions of the health indices. For instance, aortic pulse wave velocity can be estimated with a high accuracy (r > 0.9) when photopletyshmography is measured from left carotid artery using a combination of foot-to-foot pulse transmit time and peak location derived for the predictions. Similar accuracy can be reached for diastolic blood pressure, but predictions of systolic blood pressure are less accurate (r > 0.75) and the stroke volume predictions are mostly contributed by heart rate.

Recently there has been a strong trend for self-monitoring of your cardiovascular health and new wearable sport trackers and mobile applications are coming to the market everyday. However, such solutions are mostly taking advantage of heart rate measurement. Other health indices such as blood pressure and pulse wave velocity reflecting to the condition of cardiovascular system would also be of great interest, but such solutions for continuous monitoring are barely existing or are at least unreliable. In this paper, we use computational modelling to assess theoretical capabilities of such measurements. We concentrate on predicting health indices using on pulse transmit time type of measurements. Such measurements could be carried out, for example, with photopletyshmography sensor or an optical sensor already found from several wearable sport trackers. We use cardiovascular modelling to create a database of “virtual subjects”, which is applied with machine learning to construct predictors for health indices. Our findings suggest that aortic pulse wave velocity and diastolic blood pressured could be predicted with a high accuracy, but predictions of systolic blood pressure are less accurate.

In silico Optimization of Left Atrial Appendage Occluder Implantation Using Interactive and Modeling Tools

According to clinical studies, around one third of patients with atrial fibrillation (AF) will suffer a stroke during their lifetime. Between 70 and 90% of these strokes are caused by thrombus formed in the left atrial appendage. In patients with contraindications to oral anticoagulants, a left atrial appendage occluder (LAAO) is often implanted to prevent blood flow entering in the LAA. A limited range of LAAO devices is available, with different designs and sizes. Together with the heterogeneity of LAA morphology, these factors make LAAO success dependent on clinician's experience. A sub-optimal LAAO implantation can generate thrombi outside the device, eventually leading to stroke if not treated. The aim of this study was to develop clinician-friendly tools based on biophysical models to optimize LAAO device therapies. A web-based 3D interactive virtual implantation platform, so-called VIDAA, was created to select the most appropriate LAAO configurations (type of device, size, landing zone) for a given patient-specific LAA morphology. An initial LAAO configuration is proposed in VIDAA, automatically computed from LAA shape features (centreline, diameters). The most promising LAAO settings and LAA geometries were exported from VIDAA to build volumetric meshes and run Computational Fluid Dynamics (CFD) simulations to assess blood flow patterns after implantation. Risk of thrombus formation was estimated from the simulated hemodynamics with an index combining information from blood flow velocity and complexity. The combination of the VIDAA platform with in silico indices allowed to identify the LAAO configurations associated to a lower risk of thrombus formation; device positioning was key to the creation of regions with turbulent flows after implantation. Our results demonstrate the potential for optimizing LAAO therapy settings during pre-implant planning based on modeling tools and contribute to reduce the risk of thrombus formation after treatment.

Non-invasive Assessment of Systolic and Diastolic Cardiac Function During Rest and Stress Conditions Using an Integrated Image-Modeling Approach

Background: The possibility of non-invasively assessing load-independent parameters characterizing cardiac function is of high clinical value. Typically, these parameters are assessed during resting conditions. However, for diagnostic purposes, the parameter behavior across a physiologically relevant range of heart rate and loads is more relevant than the isolated measurements performed at rest. This study sought to evaluate changes in non-invasive estimations of load-independent parameters of left-ventricular contraction and relaxation patterns at rest and during dobutamine stress. Methods: We applied a previously developed approach that combines non-invasive measurements with a physiologically-based, reduced-order model of the cardiovascular system to provide subject-specific estimates of parameters characterizing left ventricular function. In this model, the contractile state of the heart at each time point along the cardiac cycle is modeled using a time-varying elastance curve. Non-invasive data, including four-dimensional magnetic resonance imaging (4D Flow MRI) measurements, were acquired in nine subjects without a known heart disease at rest and during dobutamine stress. For each of the study subjects, we constructed two personalized models corresponding to the resting and the stress state. Results: Applying the modeling framework, we identified significant increases in the left ventricular contraction rate constant [from 1.5 ± 0.3 to 2 ± 0.5 (p = 0.038)] and relaxation constant [from 37.2 ± 6.9 to 46.1 ± 12 (p = 0.028)]. In addition, we found a significant decrease in the elastance diastolic time constant from 0.4 ± 0.04 s to 0.3 ± 0.03 s (p = 0.008). Conclusions: The integrated image-modeling approach allows the assessment of cardiovascular function given as model-based parameters. The agreement between the estimated parameter values and previously reported effects of dobutamine demonstrates the potential of the approach to assess advanced metrics of pathophysiology that are otherwise difficult to obtain non-invasively in clinical practice.

A Novel Hierarchical Template Matching Model for Cardiac Motion Estimation

Cardiovascular disease diagnosis and prognosis can be improved by measuring patient-specific in-vivo local myocardial strain using Magnetic Resonance Imaging. Local myocardial strain can be determined by tracking the movement of sample muscles points during cardiac cycle using cardiac motion estimation model. The tracking accuracy of the benchmark Free Form Deformation (FFD) model is greatly affected due to its dependency on tunable parameters and regularisation function. Therefore, Hierarchical Template Matching (HTM) model, which is independent of tunable parameters, regularisation function, and image-specific features, is proposed in this article. HTM has dense and uniform points correspondence that provides HTM with the ability to estimate local muscular deformation with a promising accuracy of less than half a millimetre of cardiac wall muscle. As a result, the muscles tracking accuracy has been significantly (p < 0.001) improved (30%) compared to the benchmark model. Such merits of HTM provide reliably calculated clinical measures which can be incorporated into the decision-making process of cardiac disease diagnosis and prognosis.