Low Wall Shear Stress Is Associated with Saphenous Vein Graft Stenosis in Patients with Coronary Artery Bypass Grafting

Personalized and uncertainty-aware coronary hemodynamics simulations: From Bayesian estimation to improved multi-fidelity uncertainty quantification

Background

Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline.

Objective

We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data.

Methods

We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs.

Results

Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance.

Conclusions

The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

Superior Measurement Accuracy of Digital Thickness Gauge Versus Digital Vernier Caliper in Determining Venous Tissue Thickness

Background In determining mechanical characteristics, the accuracy of the thickness of the specimens can influence the biomechanical behavior, especially in the case of human tissues, where there is an important variability. This study aims to compare the accuracy of two routine measuring instruments, i.e., the digital vernier caliper and the digital thickness gauge, when measuring the thickness of venous specimens multiple times. Methodology In this study, we used 12 tubular vena cava specimens obtained from common breed pigs aged 18-24 weeks at the time of sacrifice from a local slaughterhouse. The measurements were performed using a digital vernier caliper (Multicomp PRO MP012475) for the first four protocols and a digital thickness gauge (Mitutoyo 547-500S) for the fifth protocol. In the first protocol, three measurements were taken on the same side, and their average was recorded as the sample thickness. The second protocol involved taking measurements on two opposite sides, and the average of these measurements was recorded as the sample thickness. In the third protocol, the thickness of each side was measured at its midpoint, and the average of the four measurements was recorded as the sample thickness. In the last protocol using a digital vernier caliper, the thickness of the vernier specimens was calculated as the average of the measurements taken at each corner of the square sample. Finally, for the fifth protocol, three consecutive measurements were taken using the digital thickness gauge, and their average was recorded as the final thickness of the sample. Results In the first protocol, we observed lower values during the first measurement in comparison to the second (0.409 ± 0.063 vs. 0.536 ± 0.064, p < 0.0001) and the third (0.409 ± 0.063 vs. 0.528 ± 0.055, p = 0.0001). Moreover, with the second protocol, we observed lower values during the first two measurements in comparison to the third measurement (p = 0.0279 and p = 0.0054). Regarding protocols three and four, we recorded higher values for the second and third measurements than the first one, with higher values for the third measurement than the second one. In the fifth protocol, there were no significant statistical differences between the three consecutive measurements (p = 0.953, p = 0.742, and p = 0.897). Further, we examined the variations in sample thickness determined using each of the protocols proposed for the digital vernier caliper, as well as the values obtained with the digital thickness gauge protocol. As a result, during the first and second measurements, we observed lower thickness values for the venous wall samples using the first four protocols compared to the fifth protocol (for all p < 0.05). However, no differences were noted between the five protocols during the third measurement. Conclusions The digital thickness gauge Mitutoyo 547-500S provided superior accuracy with no difference between three successive measurements of venous wall thickness, regardless of the examiner's experience. Accurately determining the thickness of venous specimens is crucial for calculating the tissue's biomechanical properties.

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary anatomies combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. We propose an optimization-based method to personalize multiscale coronary flow simulations by assimilating clinical CT myocardial perfusion imaging and cardiac function measurements to yield patient-specific flow distributions and model parameters. Using this proof-of-concept study on a cohort of six patients, we reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based purely on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

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.

Impact of Graft Velocity on Saphenous Vein Graft Atherosclerosis after Coronary Artery Bypass Grafting

PURPOSE

Saphenous vein grafts (SVGs) sometimes occur as vein graft stenosis or failure in coronary artery bypass grafting. The purpose of this study was to detect the factors affecting vein graft atherosclerosis.

METHODS

We performed two analysis. In the first analysis, we enrolled 120 grafts using conventionally harvested saphenous vein graft (C-SVG) and followed-up with multiple coronary computed tomography angiography (CCTA). We examined the factors that contribute to the graft atherosclerosis defined by graft failure at subsequent CCTA or substantial progression of graft stenosis (a decrease of ≥0.6 mm in diameter). In the second analysis, 66 grafts using no-touch harvested saphenous vein graft (N-SVG) were compared with those in the first analysis using C-SVG, focusing on the differences in intraoperative factors using propensity score-matched analysis.

RESULTS

In the first analysis, graft atherosclerosis+ group comprised 27 grafts, which had a larger SVG diameter, lower graft velocity, and higher graft/native ratio in diameter than the graft atherosclerosis- group. In the multivariable analysis, slow graft velocity and graft/native ≥2 in diameter were independently associated with the graft atherosclerosis. In the second analysis, the N-SVG group had a much greater graft velocity than the C-SVG group.

CONCLUSION

Lower graft velocity and higher graft/native ratio in diameter were associated with the graft atherosclerosis. The N-SVG group had increased graft velocity, which may contribute to prevent the graft atherosclerosis.(Trial registration: UMIN Clinical Trial Registry no. UMIN000050482. Registered 3 March 2023, retrospectively registered.).

Abnormal Wall Shear Stress Area is Correlated to Coronary Artery Bypass Graft Remodeling 1 Year After Surgery

Coronary artery bypass graft surgery is a common intervention for coronary artery disease; however, it suffers from graft failure, and the underlying mechanisms are not fully understood. To better understand the relation between graft hemodynamics and surgical outcomes, we performed computational fluid dynamics simulations with deformable vessel walls in 10 study participants (24 bypass grafts) based on CT and 4D flow MRI one month after surgery to quantify lumen diameter, wall shear stress (WSS), and related hemodynamic measures. A second CT acquisition was performed one year after surgery to quantify lumen remodeling. Compared to venous grafts, left internal mammary artery grafts experienced lower abnormal WSS (< 1 Pa) area one month after surgery (13.8 vs. 70.1%, p = 0.001) and less inward lumen remodeling one year after surgery (- 2.4% vs. - 16.1%, p = 0.027). Abnormal WSS area one month post surgery correlated with percent change in graft lumen diameter one year post surgery (p = 0.030). This study shows for the first time prospectively a correlation between abnormal WSS area one month post surgery and graft lumen remodeling 1 year post surgery, suggesting that shear-related mechanisms may play a role in post-operative graft remodeling and might help explain differences in failure rates between arterial and venous grafts.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Oscillatory shear stress promotes vein graft intimal hyperplasia via NADPH oxidase-related pathways

Background

Uncontrolled intimal hyperplasia (IH) after autologous saphenous vein grafting triggers a high restenosis rate; however, its association with the activation of NADPH oxidase (NOX)-related pathways is unclear. Here, we investigated the effects and mechanism of oscillatory shear stress (OSS) on grafted vein IH.

Methods

Thirty male New Zealand rabbits were randomly divided into control, high-OSS (HOSS), and low-OSS (LOSS) groups, and the vein grafts were harvested after 4 weeks. Hematoxylin and eosin staining and Masson staining assays were used to observe morphological and structural changes. Immunohistochemical staining was used to detect α-SMA, PCNA, MMP-2, and MMP-9 expression. Immunofluorescence staining was used to observe reactive oxygen species (ROS) production in the tissues. Western blotting was used to determine the expression levels of pathway-related proteins (NOX1, NOX2, AKT, p-AKT, and BIRC5), PCNA, BCL-2, BAX, and caspase-3/cleaved caspase-3 in tissues.

Results

Blood flow velocity was lower in the LOSS group than in the HOSS group, while vessel diameter did not change significantly. Shear rate was elevated in both HOSS and LOSS groups but was higher in the HOSS group. Additionally, vessel diameter increased with time in the HOSS and LOSS groups, whereas flow velocity did not. Intimal hyperplasia was significantly lower in the LOSS group than in the HOSS group. IH was dominated by smooth muscle fibers in the grafted veins and collagen fibers in the media. OSS restriction significantly reduced the α-SMA, PCNA, MMP-2, and MMP-9 levels. Moreover, ROS production and the expression of NOX1, NOX2, p-AKT, BIRC5, PCNA, BCL-2, BAX, and cleaved caspase-3 were phase-reduced in LOSS compared to the levels in the HOSS group. Total AKT was not differentially expressed among the three groups.

Conclusion

OSS promotes the proliferation, migration, and survival of subendothelial vascular smooth muscle cells in grafted veins, which may be related to the regulation of downstream p-AKT/BIRC5 levels through the increased production of ROS by NOX. Drugs inhibiting this pathway might be used to prolong vein graft survival time.

Patient-specific computational simulation of coronary artery bypass grafting

INTRODUCTION

Coronary artery bypass graft surgery (CABG) is an intervention in patients with extensive obstructive coronary artery disease diagnosed with invasive coronary angiography. Here we present and test a novel application of non-invasive computational assessment of coronary hemodynamics before and after bypass grafting.

METHODS AND RESULTS

We tested the computational CABG platform in n = 2 post-CABG patients. The computationally calculated fractional flow reserve showed high agreement with the angiography-based fractional flow reserve. Furthermore, we performed multiscale computational fluid dynamics simulations of pre- and post-CABG under simulated resting and hyperemic conditions in n = 2 patient-specific anatomies 3D reconstructed from coronary computed tomography angiography. We computationally created different degrees of stenosis in the left anterior descending artery, and we showed that increasing severity of native artery stenosis resulted in augmented flow through the graft and improvement of resting and hyperemic flow in the distal part of the grafted native artery.

CONCLUSIONS

We presented a comprehensive patient-specific computational platform that can simulate the hemodynamic conditions before and after CABG and faithfully reproduce the hemodynamic effects of bypass grafting on the native coronary artery flow. Further clinical studies are warranted to validate this preliminary data.

Wall Shear Stress Differences Between Arterial and Venous Coronary Artery Bypass Grafts One Month After Surgery

Although coronary artery bypass graft (CABG) surgery is a well-established intervention, graft failure can occur, and the underlying mechanisms remain incompletely understood. The purpose of this prospective study is to utilize computational fluid dynamics (CFD) to investigate how graft hemodynamics one month post surgery may vary among graft types, which have different long-term patency rates. Twenty-four grafts from 10 participants (64.6 ± 8.5 years, 9 men) were scanned with coronary CT angiography and 4D flow MRI one month after CABG surgery. Grafts included 10 left internal mammary arteries (LIMA), 3 radial arteries (RA), and 11 saphenous vein grafts (SVG). Image-guided CFD was used to quantify blood flow rate and wall area exposed to abnormal wall shear stress (WSS). Arterial grafts had a lower abnormal WSS area than venous grafts (17.9% vs. 70.1%; p = 0.001), and a similar trend was observed for LIMA vs. SVG (13.8% vs. 70.1%; p = 0.001). Abnormal WSS area correlated positively to lumen diameter (p < 0.001) and negatively to flow rate (p = 0.001). This CFD study is the first of its kind to prospectively reveal differences in abnormal WSS area 1 month post surgery among CABG types, suggesting that WSS may influence the differential long-term graft failure rates observed among these groups.

Colocalization of Coronary Plaque with Wall Shear Stress in Myocardial Bridge Patients

PURPOSE

Patients with myocardial bridges (MBs) have a higher prevalence of atherosclerosis. Wall shear stress (WSS) has previously been correlated with plaque in coronary artery disease patients, but such correlations have not been investigated in symptomatic MB patients. The aim of this paper was to use a multi-scale computational fluid dynamics (CFD) framework to simulate hemodynamics in MB patient, and investigate the co-localization of WSS and plaque.

METHODS

We identified N = 10 patients from a previously reported cohort of 50 symptomatic MB patients, all of whom had plaque in the proximal vessel. Dynamic 3D models were reconstructed from coronary computed tomography angiography (CCTA), intravascular ultrasound (IVUS) and catheter angiograms. CFD simulations were performed to compute WSS proximal to, within and distal to the MB. Plaque was quantified from IVUS images in 2 mm segments and registered to CFD model. Plaque area was compared to absolute and patient-normalized WSS.

RESULTS

WSS was lower in the proximal segment compared to the bridge segment (6.1 ± 2.9 vs. 16.0 ± 7.1 dynes/cm², p value < 0.01). Plaque area and plaque burden measured from IVUS peaked at 1-3 cm proximal to the MB entrance, coinciding with the first diagonal branch. Normalized WSS showed a statistically significant moderate correlation with plaque area (r = 0.41, p < 0.01).

CONCLUSION

WSS may be obtained non-invasively in MB patients and provides a surrogate marker of plaque area. Using CFD, it may be possible to non-invasively assess the extent of plaque area, and identify patients who could benefit from frequent monitoring or medical management.

Small Diameter Cell-Free Tissue-Engineered Vascular Grafts: Biomaterials and Manufacture Techniques to Reach Suitable Mechanical Properties

Vascular grafts (VGs) are medical devices intended to replace the function of a blood vessel. Available VGs in the market present low patency rates for small diameter applications setting the VG failure. This event arises from the inadequate response of the cells interacting with the biomaterial in the context of operative conditions generating chronic inflammation and a lack of regenerative signals where stenosis or aneurysms can occur. Tissue Engineered Vascular grafts (TEVGs) aim to induce the regeneration of the native vessel to overcome these limitations. Besides the biochemical stimuli, the biomaterial and the particular micro and macrostructure of the graft will determine the specific behavior under pulsatile pressure. The TEVG must support blood flow withstanding the exerted pressure, allowing the proper compliance required for the biomechanical stimulation needed for regeneration. Although the international standards outline the specific requirements to evaluate vascular grafts, the challenge remains in choosing the proper biomaterial and manufacturing TEVGs with good quality features to perform satisfactorily. In this review, we aim to recognize the best strategies to reach suitable mechanical properties in cell-free TEVGs according to the reported success of different approaches in clinical trials and pre-clinical trials.

Patient-Specific Fluid-Structure Simulations of Anomalous Aortic Origin of Right Coronary Arteries

Objectives

Anomalous aortic origin of the right coronary artery (AAORCA) may cause ischemia and sudden death. However, the specific anatomic indications for surgery are unclear, so dobutamine-stress instantaneous wave-free ratio (iFR) is increasingly used. Meanwhile, advances in fluid–structure interaction (FSI) modeling can simulate the pulsatile hemodynamics and tissue deformation. We sought to evaluate the feasibility of simulating the resting and dobutamine-stress iFR in AAORCA using patient-specific FSI models and to visualize the mechanism of ischemia within the intramural geometry and associated lumen narrowing.

Methods

We developed 6 patient-specific FSI models of AAORCA using SimVascular software. Three-dimensional geometries were segmented from coronary computed tomography angiography. Vascular outlets were coupled to lumped-parameter networks that included dynamic compression of the coronary microvasculature and were tuned to each patient's vitals and cardiac output.

Results

All cases were interarterial, and 5 of 6 had an intramural course. Measured iFRs ranged from 0.95 to 0.98 at rest and 0.80 to 0.95 under dobutamine stress. After we tuned the distal coronary resistances to achieve a stress flow rate triple that at rest, the simulations adequately matched the measured iFRs (r = 0.85, root-mean-square error = 0.04). The intramural lumen remained narrowed with simulated stress and resulted in lower iFRs without needing external compression from the pulmonary root.

Conclusions

Patient-specific FSI modeling of AAORCA is a promising, noninvasive method to assess the iFR reduction caused by intramural geometries and inform surgical intervention. However, the models’ sensitivity to distal coronary resistance suggests that quantitative stress-perfusion imaging may augment virtual and invasive iFR studies.

Creating patient-specific vein models to characterize wall shear stress in hemodialysis population

End-Stage Renal Disease (ESRD) patients require arteriovenous fistulas (AVF) that allow a mature vein to withstand hemodialysis. Unfortunately, venous thrombosis and stenosis in the cephalic vein arch after AVF placement is common and heavily influenced by hemodynamics. To better assess forces and flow behavior in the cephalic arch, we have built patient-specific millifluidic models that allow us to explore the complex interplay between patient-specific vein geometry and fluctuating hemodynamics. These 3D models were created from patient-specific intravascular ultrasound and venogram images obtained three- and twelve-months post AVF creation and fabricated into soft elastomer-based millifluidic devices. Geometric validation of fabricated phantom millifluidic device shows successful replication of original computational 3D model. Millifluidic devices were perfused with a blood-mimicking fluid containing fluorescent tracer beads under steady-state physiologic cephalic vein flow conditions (20 mL/min). Particle image velocimetry was employed to calculate wall shear stress (WSS) across the cephalic arches. Experimental WSS profile evaluation reveals that the physiologic cephalic arch model yields WSS values within physiologic range [76–760 mPa]. Moreover, upon comparing WSS profiles across all models, it is noticeable that WSS values increase as vein diameter decreases, which further supports employed experimental and analysis strategy. The presented millifluidic devices show promise for experimental WSS characterization under pathologic flow conditions to contrast from calculated physiologic hemodynamics and better understand WSS influence on thrombosis and stenosis in hemodialysis patients.

Morphological and hemodynamic analysis of the patient-specific renal cell carcinoma models

Although the morbidity of renal cell carcinoma (RCC) has been increasing as the seventh most common tumours, to our knowledge, there is few studies foucsing on the hemodynamics of the renal artery (RA) with RCC. The objective of this study is to perform morphological and hemodynamic analysis of the RA and abdominal aorta artery (AAA) in the control healthy and RCC patient groups. Three-dimensional (3D) geometries are reconstructed from 18 control healthy subjects and 15 RCC patients based on Computed Tomography Angiography (CTA) images. There is higer in the lumen diameter of the RA (6.21 ± 0.89 mm) and curvature of the RA (1.2 ± 0.07) in the RCC patient group compared with the control healthy group (4.29 ± 1.08 mm, 1.1 ± 0.1), respectively. In the hemodynamic analysis, the surface area ratio (%) of low time-averaged wall shear stress (SAR-TAWSS) at the RA (10.65 ± 11.65) and AAA (48.49 ± 12.79) in the RCC patient group is significantly higher than that in the control healthy group (0.23 ± 0.22, 21.57 ± 20.5), respectively. It is found that RCC altered the morphology of the RA in the RCC patient group, which could deteriorate the hemodynamic environment of the RA and AAA. The finding in this study could enhance us to understand the progression of vascular disease caused by RCC.

Computational Evaluation of Venous Graft Geometries in Coronary Artery Bypass Surgery

Cardiothoracic surgeons are faced with a choice of different revascularization techniques and diameters for saphenous vein grafts (SVG) in coronary artery bypass graft surgery . Using computational simulations, we virtually investigate the effect of SVG geometry on hemodynamics of both venous grafts and the target coronary arteries. We generated patient-specific 3-dimensional anatomic models of coronary artery bypass graft surgery patients and quantified mechanical stimuli. We performed virtual surgery on 3 patient-specific models by modifying the geometry vein grafts to reflect single, Y, and sequential surgical configurations with SVG diameters ranging from 2 mm to 5 mm. Our study demonstrates that the coronary artery runoffs are relatively insensitive to the choice of SVG revascularization geometry. We observe a 10% increase in runoff when the SVG diameter is changed from 2 mm to 5 mm. The wall shear stress of SVG increases dramatically when the diameter drops, following an inverse power scaling with diameter. For a fixed diameter, the average wall shear stress on the vein graft varies in ascending order as single, Y, and sequential graft in the patient cohort. The runoff to the target coronary arteries changes marginally due to the choice of graft configuration or diameter. The shear stress on the vein graft depends on both flow rate and diameter and follows an inverse power scaling consistent with a Poiseuille flow assumption. Given the similarity in runoff with different surgical configurations, choices of SVG geometries can be informed by propensity for graft failure using shear stress evaluations.