Computational modeling of Fontan physiology: at the crossroads of pediatric cardiology and biomedical engineering

Virtual treatment planning in three patients with univentricular physiology using computational fluid dynamics—Pitfalls and strategies

Background

Uneven hepatic venous blood flow distribution (HFD) to the pulmonary arteries is hypothesized to be responsible for the development of intrapulmonary arteriovenous malformations (PAVM) in patients with univentricular physiology. Thus, achieving uniform distribution of hepatic blood flow is considered favorable. However, no established method for the prediction of the post-interventional hemodynamics currently exists. Computational fluid dynamics (CFD) offers the possibility to quantify HFD in patient-specific anatomies before and after virtual treatment. In this study, we evaluated the potential benefit of CFD-assisted treatment planning.

Materials and methods

Three patients with total cavopulmonary connection (TCPC) and PAVM underwent cardiovascular magnetic resonance imaging (CMR) and computed tomography imaging (CT). Based on this imaging data, the patient-specific anatomy was reconstructed. These patients were considered for surgery or catheter-based intervention aiming at hepatic blood flow re-routing. CFD simulations were then performed for the untreated state as well as for different surgical and interventional treatment options. These treatment options were applied as suggested by treating cardiologists and congenital heart surgeons with longstanding experience in interventional and surgical treatment of patients with univentricular physiology. HFD was quantified for all simulations to identify the most viable treatment decision regarding redistribution of hepatic blood flow.

Results

For all three patients, the complex TCPC anatomy could be reconstructed. However, due to the presence of metallic stent implants, hybrid models generated from CT as well as CMR data were required. Numerical simulation of pre-interventional HFD agreed well with angiographic assessment and physiologic considerations. One treatment option resulting in improvement of HFD was identified for each patient. In one patient follow-up data after treatment was available. Here, the virtual treatment simulation and the CMR flow measurements differed by 15%.

Conclusion

The combination of modern computational methods as well as imaging methods for assessment of patient-specific anatomy and flow might allow to optimize patient-specific therapy planning in patients with pronounced hepatic flow mismatch and PAVM. In this study, we demonstrate that these methods can also be applied in patients with complex univentricular physiology and extensive prior interventions. However, in those cases, hybrid approaches utilizing information of different image modalities may be required.

Non-uniform mixing of hepatic venous flow and inferior vena cava flow in the Fontan conduit

Fontan patients require a balanced hepatic blood flow distribution (HFD) to prevent pulmonary arteriovenous malformations. Currently, HFD is quantified by tracking Fontan conduit flow, assuming hepatic venous (HV) flow to be uniformly distributed within the Fontan conduit. However, this assumption may be unvalid leading to inaccuracies in HFD quantification with potential clinical impact. The aim of this study was to (i) assess the mixing of HV flow and inferior vena caval (IVC) flow within the Fontan conduit and (ii) quantify HFD by directly tracking HV flow and quantitatively comparing results with the conventional approach. Patient-specific, time-resolved computational fluid dynamic models of 15 total cavopulmonary connections were generated, including the HV and subhepatic IVC. Mixing of HV and IVC flow, on a scale between 0 (no mixing) and 1 (perfect mixing), was assessed at the caudal and cranial Fontan conduit. HFD was quantified by tracking particles from the caudal (HFD_{caudal conduit}) and cranial (HFD_{cranial conduit}) conduit and from the hepatic veins (HFD_HV). HV flow was non-uniformly distributed at both the caudal (mean mixing 0.66 ± 0.13) and cranial (mean 0.79 ± 0.11) level within the Fontan conduit. On a cohort level, differences in HFD between methods were significant but small; HFD_HV (51.0 ± 20.6%) versus HFD_{caudal conduit} (48.2 ± 21.9%, p = 0.033) or HFD_{cranial conduit} (48.0 ± 21.9%, p = 0.044). However, individual absolute differences of 8.2–14.9% in HFD were observed in 4/15 patients. HV flow is non-uniformly distributed within the Fontan conduit. Substantial individual inaccuracies in HFD quantification were observed in a subset of patients with potential clinical impact.

Hypoplastic left heart syndrome: From bedside to bench and back

Hypoplastic Left Heart Syndrome (HLHS) is a complex Congenital Heart Disease (CHD) that was almost universally fatal until the advent of the Norwood operation in 1981. Children with HLHS who largely succumbed to the disease within the first year of life, are now surviving to adulthood. However, this survival is associated with multiple comorbidities and HLHS infants have a higher mortality rate as compared to other non-HLHS single ventricle patients. In this review we (a) discuss current clinical challenges associated in the care of HLHS patients, (b) explore the use of systems biology in understanding the molecular framework of this disease, (c) evaluate induced pluripotent stem cells as a translational model to understand molecular mechanisms and manipulate them to improve outcomes, and (d) investigate cell therapy, gene therapy, and tissue engineering as a potential tool to regenerate hypoplastic cardiac structures and improve outcomes.

Outcomes of Single Ventricle Patients Undergoing the Kawashima Procedure: Can We Do Better?

Objectives:

Current technology advances in virtual surgery modeling and computational flow dynamics allow preoperative individualized computer-based design of Fontan operation. To determine potential role of those innovations in patients undergoing hepatic vein incorporation (HVI) following Kawashima operation, we retrospectively examined historic cohort of patients who underwent HVI following Kawashima with focus on regression of pulmonary arteriovenous malformations (PAVMs).

Methods:

Twenty-two children with single ventricle and interrupted inferior vena cava underwent Kawashima operation (2002-12). Twenty-one (96%) patients had left atrial isomerism and 21 (96%) had undergone prior first-stage palliation. Clinical outcomes were examined.

Results:

Mean O2 saturation (SaO2) increased from 77% ± 8% to 85% ± 6% ( P = .002) after Kawashima. Fifteen (68%) patients developed PAVMs. Eighteen patients underwent HVI (median age and interval from Kawashima: 4.4 and 3.7 years, respectively). Mean SaO2 prior to HVI was 77% ± 8% and increased to 81% ± 10% at the time of hospital discharge ( P = .250), with five patients requiring home oxygen. On follow-up, mean SaO2 increased to 95% ± 4% ( P < .001). Overall ten-year survival following Kawashima was 94%.

Conclusions:

A large number of patients develop PAVMs and subsequent cyanosis after Kawashima operation. Early following HVI, SaO2 is commonly low and insignificantly different from that prior to HVI. Although SaO2 will improve on follow-up in most patients, a number of patients continue to have low saturations, indicating incomplete resolution of PAVMs. Given the heterogeneity of those patients and lack of preoperative predictors for complete PAVM regression, our findings suggest a role for virtual surgery to determine optimal individual procedure design that would provide even distribution of hepatic blood flow to both pulmonary arteries.

The first cohort of prospective Fontan surgical planning patients with follow-up data: How accurate is surgical planning?

OBJECTIVE

Fontan surgical planning is an image-based, collaborative effort, which is hypothesized to result in improved patient outcomes. A common motivation for Fontan surgical planning is the progression (or concern for progression) of pulmonary arteriovenous malformations. The purpose of this study was to evaluate the accuracy of surgical planning predictions, specifically hepatic flow distribution (HFD), a known factor in pulmonary arteriovenous malformation progression, and identify methodological improvements needed to increase prediction accuracy.

METHODS

Twelve single-ventricle patients who were enrolled in a surgical planning protocol for Fontan surgery with pre- and postoperative cardiac imaging were included in this study. Computational fluid dynamics were used to compare HFD in the surgical planning prediction and actual postoperative conditions.

RESULTS

Overall, HFD prediction error was 17 ± 13%. This error was similar between surgery types (15 ± 18% and 18 ± 10% for revisions vs Fontan completions respectively; P = .73), but was significantly lower (6 ± 7%; P = .05) for hepatic to azygous shunts. Y-grafts and extracardiac conduits showed a strong correlation between prediction error and discrepancies in graft insertion points (r = 0.99; P < .001). Improving postoperative anatomy prediction significantly reduced overall HFD prediction error to 9 ± 6% (P = .03).

CONCLUSIONS

Although Fontan surgical planning can offer accurate HFD predictions for specific graft types, methodological improvements are needed to increase overall accuracy. Specifically, improving postoperative anatomy prediction was shown to be an important target for future work. Future efforts and refinements to the surgical planning process will benefit from an improved understanding of the current state and will rely heavily on increased follow-up data.

Aortic Arch Homograft Reconstruction of Nonconfluent Pulmonary Arteries During Extracardiac Fontan

Reconstruction of nonconfluent pulmonary arteries during Fontan completion is a challenging technical issue. In this case report, we describe the use of an aortic homograft, including the aortic arch, to complete a Fontan and reconstruct the pulmonary artery confluence in a child with discontinuous pulmonary arteries and bilateral superior caval veins who had undergone bilateral unidirectional Glenn palliation. The configuration of the aortic homograft was ideal to ensure laminar flow from the inferior vena cava to both pulmonary arteries and in maintaining durable elastance posterior to the native aorta.

Fontan Surgical Planning: Previous Accomplishments, Current Challenges, and Future Directions

The ultimate goal of Fontan surgical planning is to provide additional insights into the clinical decision-making process. In its current state, surgical planning offers an accurate hemodynamic assessment of the pre-operative condition, provides anatomical constraints for potential surgical options, and produces decent post-operative predictions if boundary conditions are similar enough between the pre-operative and post-operative states. Moving forward, validation with post-operative data is a necessary step in order to assess the accuracy of surgical planning and determine which methodological improvements are needed. Future efforts to automate the surgical planning process will reduce the individual expertise needed and encourage use in the clinic by clinicians. As post-operative physiologic predictions improve, Fontan surgical planning will become an more effective tool to accurately model patient-specific hemodynamics.

A Method for In Vitro TCPC Compliance Verification

The Fontan procedure is a common palliative intervention for sufferers of single ventricle congenital heart defects that results in an anastomosis of the venous return to the pulmonary arteries called the total cavopulmonary connection (TCPC). Local TCPC and global Fontan circulation hemodynamics are studied with in vitro circulatory models because of hemodynamic ties to Fontan patient long-term complications. The majority of in vitro studies, to date, employ a rigid TCPC model. Recently, a few studies have incorporated flexible TCPC models, but provide no justification for the model material properties. The method set forth in this study successfully utilizes patient-specific flow and pressure data from phase contrast magnetic resonance images (PCMRI) (n = 1) and retrospective pulse-pressure data from an age-matched patient cohort (n = 10) to verify the compliance of an in vitro TCPC model. These data were analyzed, and the target compliance was determined as 1.36 ± 0.78 mL/mm Hg. A method of in vitro compliance testing and computational simulations was employed to determine the in vitro flexible TCPC model material properties and then use those material properties to estimate the wall thickness necessary to match the patient-specific target compliance. The resulting in vitro TCPC model compliance was 1.37 ± 0.1 mL/mm Hg—a value within 1% of the patient-specific compliance. The presented method is useful to verify in vitro model accuracy of patient-specific TCPC compliance and thus improve patient-specific hemodynamic modeling.

Evaluation of blood flow distribution asymmetry and vascular geometry in patients with Fontan circulation using 4-D flow MRI

BACKGROUND

Asymmetrical caval to pulmonary blood flow is suspected to cause complications in patients with Fontan circulation. The aim of this study was to test the feasibility of 4-D flow MRI for characterizing the relationship between 3-D blood flow distribution and vascular geometry.

OBJECTIVE

We hypothesized that both flow distribution and geometry can be calculated with low interobserver variability and will detect a direct relationship between flow distribution and Fontan geometry.

MATERIALS AND METHODS

Four-dimensional flow MRI was acquired in 10 Fontan patients (age: 16 ± 4 years [mean ± standard deviation], range: 9-21 years). The Fontan connection was isolated by 3-D segmentation to evaluate flow distribution from the inferior vena cava (IVC) and superior vena cava (SVC) to the left and right pulmonary arteries (LPA, RPA) and to characterize geometry (cross-sectional area, caval offset, vessel angle).

RESULTS

Flow distribution results indicated SVC flow tended toward the RPA while IVC flow was more evenly distributed (SVC to RPA: 78% ± 28 [9-100], IVC to LPA: 54% ± 28 [4-98]). There was a significant relationship between pulmonary artery cross-sectional area and flow distribution (IVC to RPA: R(2)=0.50, P=0.02; SVC to LPA: R(2)=0.81, P=0.0004). Good agreement was found between observers and for flow distribution when compared to net flow values.

CONCLUSION

Four-dimensional flow MRI was able to detect relationships between flow distribution and vessel geometry. Future studies are warranted to investigate the potential of patient specific hemodynamic analysis to improve diagnostic capability.