Multiscale Modeling of Superior Cavopulmonary Circulation: Hemi-Fontan and Bidirectional Glenn Are Equivalent

Advancement in computational simulation and validation of congenital heart disease: a review

The improvement in congenital heart disease (CHD) treatment and management has increased the life expectancy in infants. However, the long-term efficacy is difficult to assess and thus, computational modelling has been applied for evaluating this. Here, we provide an overview of the applications of computational modelling in CHD based on three categories; CHD involving large blood vessels only, heart chambers only, and CHD that occurs at multiple heart structures. We highlight the advancement of computational simulation of CHD that uses multiscale and multiphysics modelling to ensure a complete representation of the heart and circulation. We provide a brief future direction of computational modelling of CHD such as to include growth and remodelling, detailed conduction system, and occurrence of myocardial infarction. We also proposed validation technique using advanced three-dimensional (3D) printing and particle image velocimetry (PIV) technologies to improve the model accuracy.

Hemodynamics and Wall Mechanics of Vascular Graft Failure

Blood vessels are subjected to complex biomechanical loads, primarily from pressure-driven blood flow. Abnormal loading associated with vascular grafts, arising from altered hemodynamics or wall mechanics, can cause acute and progressive vascular failure and end-organ dysfunction. Perturbations to mechanobiological stimuli experienced by vascular cells contribute to remodeling of the vascular wall via activation of mechanosensitive signaling pathways and subsequent changes in gene expression and associated turnover of cells and extracellular matrix. In this review, we outline experimental and computational tools used to quantify metrics of biomechanical loading in vascular grafts and highlight those that show potential in predicting graft failure for diverse disease contexts. We include metrics derived from both fluid and solid mechanics that drive feedback loops between mechanobiological processes and changes in the biomechanical state that govern the natural history of vascular grafts. As illustrative examples, we consider application-specific coronary artery bypass grafts, peripheral vascular grafts, and tissue-engineered vascular grafts for congenital heart surgery as each of these involves unique circulatory environments, loading magnitudes, and graft materials.

Hemi-Fontan and bidirectional Glenn operations result in flow-mediated viscous energy loss at the time of stage II palliation

Background

Superior cavopulmonary connection (SCPC) for stage II palliation of hypoplastic left heart syndrome (HLHS) is achieved most frequently by either a bidirectional Glenn (BDG) or hemi-Fontan (HF) operation. The comparison of flow hemodynamic efficiency at the region of surgical reconstruction and in proximal pulmonary arteries has been evaluated primarily using computational modeling techniques with conflicting reports. The purpose of this descriptive study was to compare flow hemodynamics following stage II (BDG vs HF) using 4-dimensional flow magnetic resonance imaging (4D-Flow MRI) with particular focus on flow-mediated viscous energy loss (EL') under matched hemodynamic conditions.

Methods

Patients with hypoplastic left heart syndrome (HLHS) who underwent either HF or BDG as part of stage II palliation underwent pre-Fontan 4D-Flow MRI. Patients were matched by the pulmonary vascular resistance index, net superior vena cava (SVC) flow, right pulmonary artery (RPA) and left pulmonary artery (LPA) size, and age. Maximum EL' throughout the cardiac cycle was calculated along the SVC-RPA and SVC-LPA tracts.

Results

Eight patients who underwent HF as part of their stage II single ventricle palliation were matched with 8 patients who underwent BDG. There were no differences between the 2 groups in median volumetric indices, including end-diastolic volume (P = .278) and end-systolic volume (P = .213). Moreover, no differences were observed in ejection fraction (P = .091) and cardiac index (P = .324). There also were no differences in peak EL' measured along the SVC-RPA tract (median, 0.05 mW for HF vs 0.04 mW for BDG; P = .365) or along the SVC-LPA tract (median, 0.05 mW vs 0.04 mW; P = .741).

Conclusions

The second stage of surgical palliation of HLHS using either HF or BDG results in similar flow-mediated viscous energy loss throughout the SCPC junction. 4D-Flow MRI and computational methods should be applied together to investigate flow hemodynamic patterns throughout the Fontan palliation and overall efficiency of the Fontan circuit.

Validation of a laser projection platform for the preparation of surgical patches used in paediatric cardiac surgery

OBJECTIVES

Reconstruction of cardiovascular anatomy with patch material is integral to the repair of congenital heart disease. We present validation of a laser projection platform for the preparation of surgical patches as a proof-of-concept for intraoperative use in patient-specific planning of paediatric cardiac surgery reconstructions.

METHODS

The MicroLASERGUIDE, a compact laser projection system that displays computer-aided designs onto 2D/3D surfaces, serves as an alternative to physical templates. A non-inferiority comparison of dimensional measurements was conducted between laser projection ('laser') and OZAKI AVNeo Template ('template') methods in creation of 51 (each group) size 13 valve leaflets from unfixed bovine pericardium. A digital version of the OZAKI AVNeo Template dimensions served as control. Feasibility testing was performed with other common patch materials (fixed bovine pericardium, PTFE and porcine main pulmonary artery as a substitute for pulmonary homograft) and sizes (13, 23) (n = 3 each group).

RESULTS

Compared to control (height 21.5, length 21.0 mm), template height and length were smaller (height and length differences of -0.3 [-0.5 to 0.0] and -0.4 [-0.8 to -0.1] mm, P < 0.01 each); whereas, both laser height and length were relatively similar (height and length differences of height 0.0 [-0.2 to 0.2], P = 0.804, and 0.2 [-0.1 to 0.4] mm, P = 0.029). Template percent error for height and length was -1.5 (-2.3 to 0.0)% and -1.9 (-3.7 to -0.6)% vs 0.2 (-1.0 to 1.1)% and 1.0 (-0.5 to 1.8)% for the laser. Similar results were found with other materials and sizes. Overall, laser sample dimensions differed by a maximum of 5% (∼1 mm) from the control.

CONCLUSIONS

The laser projection platform has demonstrated promise as an alternative methodology for the preparation of surgical patches for use in cardiac surgery. This technology has potential to revolutionize preoperative surgical planning for numerous congenital anomalies that require patient-specific patch-augmented repair.

Determination of Contrast Timing by Time-Resolved Magnetic Resonance Angiography in Patients With Bidirectional Glenn and Hemi-Fontan Anastomoses

PURPOSE

Children with single-ventricle congenital heart disease undergo a series of operations to maintain their pulmonary circulation including bidirectional Glenn (BDG) or hemi-Fontan in the second stage to create a superior cavopulmonary anastomosis. We aimed to optimize cardiovascular angiography protocols by determining optimal contrast timing of pulmonary and systemic circulation on magnetic resonance angiography (MRA) performed with the technique of time-resolved imaging with interleaved stochastic trajectories (TWIST).

METHODS AND MATERIALS

Cardiac TWIST MRA with lower extremity (LE) contrast injection was analyzed in 92 consecutive patients with a BDG or hemi-Fontan anastomosis. Contrast arrival time to inferior vena cava was set to zero to determine the relative time-to-peak (TTP) of the target vessels. Time-to-peak of each vessel was compared by age (<2 or ≥2 y), ejection fraction (<54% or ≥54%), the median values of heart rate (<111 or ≥111 beats per minute), body surface area (BSA, <0.59 or ≥0.59), cardiac index (<6.04 or ≥6.04), and indexed ascending aorta flow (AscAo_i, <5.3 or ≥5.3). The TTP of the vessels was also correlated with the volumetric parameters.

RESULTS

The mean age of 92 patients (32 female, 60 male) was 3.1 years (0.7-5.6 years). With LE injection, the first peak was depicted in AscAo. Time-to-peak of the pulmonary arteries was approximately 9 seconds later than AscAo. The TTP difference between pulmonary arteries and AscAo was shorter in high heart rate group (8.3 vs 10 seconds, P < 0.001). The TTP difference between AscAo and the mean of pulmonary arteries was significantly shorter in high cardiac index group (8.4 vs 9.9 seconds, P < 0.01) and high AscAo_i group (8.7 vs 9.7 seconds, P = 0.03). The TTP differences were not significant by age, ejection fraction, and BSA. Cardiac index and AscAo_i were negatively correlated with all TTPs except AscAo. The ejection fraction, stroke volume, and atrioventricular regurgitation fraction did not correlate with the TTP.

CONCLUSIONS

In patients with BDG or hemi-Fontan anastomosis, TTP of the pulmonary arteries on TWIST MRA via LE intravenous injection is approximately 9 seconds later than AscAo, approximately 8 and 10 seconds later in high and low heart rate groups, respectively. Cardiac index and AscAo_i have less effect on the TTP than the heart rate. There was no TTP difference of the pulmonary arteries by age, BSA, and ejection fraction and no correlation with ejection fraction, stroke volume, and atrioventricular regurgitation fraction. These data can be used to guide timing of pulmonary arterial enhancement of single-ventricle patients after BDG or hemi-Fontan anastomosis.

Hemodynamic Response to Device Titration in the Shunted Single Ventricle Circulation: A Patient Cohort Modeling Study

Clinical outcomes of ventricular assist device (VAD) support for shunted single ventricle patients trail the larger population due in part to the challenges in optimizing VAD support and balancing systemic and pulmonary circulations. We sought to understand the response to VAD titration in the shunted circulation using a lumped-parameter network modeling six patient-specific clinical cases. Hemodynamic data from six patients (mean body surface area = 0.30 m2) with a systemic-to-pulmonary shunt was used to construct simulated cases of heart failure and hemodynamic response to increasing VAD flow from 5 to 10 L/min/m2. With increasing VAD flow, the pulmonary arterial pressure stayed relatively constant in five patient cases and increased in one patient case. The mean VAD flow needed to attain an arterial-venous O2 saturation difference of 30% was 6.5 ± 1.2 L/min/m2, which is higher than that in the equivalent nonshunted scenario due to the partial diversion of flow to the pulmonary circulation. The hemodynamic responses to VAD support can vary significantly between specific patient cases; therefore hemodynamic modeling may help guide an individualized approach to perioperative VAD management in the shunted single-ventricle circulation and to understand the patients who may benefit the most from VAD support.

Engineering Perspective on Cardiovascular Simulations of Fontan Hemodynamics: Where Do We Stand with a Look Towards Clinical Application

BACKGROUND

Cardiovascular simulations for patients with single ventricles undergoing the Fontan procedure can assess patient-specific hemodynamics, explore surgical advances, and develop personalized strategies for surgery and patient care. These simulations have not yet been broadly accepted as a routine clinical tool owing to a number of limitations. Numerous approaches have been explored to seek innovative solutions for improving methodologies and eliminating these limitations.

PURPOSE

This article first reviews the current state of cardiovascular simulations of Fontan hemodynamics. Then, it will discuss the technical progress of Fontan simulations with the emphasis of its clinical impact, noting that substantial improvements have been made in the considerations of patient-specific anatomy, flow, and blood rheology. The article concludes with insights into potential future directions involving clinical validation, uncertainty quantification, and computational efficiency. The advancements in these aspects could promote the clinical usage of Fontan simulations, facilitating its integration into routine clinical practice.

Computational Modeling to Support Surgical Decision Making in Single Ventricle Physiology

Many of the advances in congenital heart surgery were built upon lessons and insights gained from model simulations. While animal and mock-circuit models have historically been the main arena to test new operative techniques and concepts, the recognition that complex cardiovascular anatomy and circulation can be modeled mathematically ushered a new era of collaboration between surgeons and engineers. In 1996, the computational age in congenital heart surgery began when investigators in London and Milan tapped the power of the computer to simulate the Fontan procedure and introduced operative improvements. Since then, computational modeling has led to numerous contributions in congenial heart surgery as continuing sophistication and advances in numerical and imaging methods furthered the ability to refine anatomic and physiologic details. Idealized generic models have given way to precise patient-specific simulations of the 3-dimensional anatomy, reconstructed circulation, affected hemodynamics, and altered physiology. Tools to perform virtual surgery, and predict flow dynamic and circulatory results, have been developed for some of the most complex defects, such as those requiring single ventricle palliation. In today's quest for personalized medicine and precision care, computational modeling's role to assist surgical planning in complex congenital heart surgery will continue to grow and evolve. With ever closer collaboration between surgeons and engineers, and clear understanding of modeling limitations, computational simulations can be a valuable adjunct to support preoperative surgical decision making.