Computational study of the risk of restenosis in coronary bypasses

The Outcomes of Coronary Artery Bypass Surgery after 18 Months-Is There an Influence of the Initial Right Ventricle Diastolic Dysfunction?

BACKGROUND

This study aimed to investigate the association of preoperative right heart filling indicators with outcomes after coronary artery bypass grafting (CABG) at an 18 month follow up.

METHODS

Patients who underwent CABG at a single center were included in this study. In addition to the baseline preoperative indicators and perioperative data, initial parameters of the right ventricle (RV) systolic and diastolic function were assessed.

RESULTS

Among the 189 patients, a total of 19 (10.0%) MACE (cardiovascular death, nonfatal myocardial infarction and stroke) were recorded during an 18 month follow up. In patients with the development of MACE during the initial examination, the following changes in RV function were revealed compared with the group without MACE: a decrease in the e't index (8.2 versus 9.6 cm/s, p = 0.029), an increase in the Et/e't ratio (5.25 vs. 4.42, p = 0.049) and more frequent presence of RV pseudonormal filling (p = 0.03). In the binary logistic regression analysis, the development of MACE 18 months after CABG was associated with the nonconduction of PCI before surgery, the presence of peripheral atherosclerosis, an increase in IVST and Et/e't and a decrease in LVEF.

CONCLUSIONS

RV diastolic dysfunction in the preoperative period was associated with the development of MACE within 18 months after CABG, and the ratio Et/e't was one of the independent predictors of MACE in a multiple regression analysis. This makes it expedient to include an assessment of not only systolic but also diastolic RV function in the preoperative examination. The inclusion of an assessment of RV diastolic function in the pre-CABG evaluation of patients deserves further study.

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.

Numerical analysis of stenoses severity and aortic wall mechanics in patients with supravalvular aortic stenosis

Supravalvular aortic stenosis (SVAS) is an aortic malformation characterized by a narrowing of the ascending aorta, resulting in abnormal hemodynamics and pressure drop across the stenosed region. It has been observed that the pressure drops measured from Doppler ultrasound exams often tend to be higher than those obtained from invasive cardiac catheterization. These misleadingly elevated pressure measurements may drive the decision to refer patients for surgical treatment prematurely. Considering this strong clinical association, the purpose of this work is to develop a computational modeling approach using a two-way coupled fluid-structure interaction methodology to determine an accurate prediction of trans-stenotic pressure drop and to further highlight the discrepancy between the SVAS assessment methods. Blood is modeled using Navier-Stokes equations while the aortic wall is simulated by a composite poroelastic structure to represent the three main layers of the arterial wall. The relationship between aortic wall elasticity and the blood flow conditions is examined in varying levels of stenosis, ranging from mild to severe degrees of vessel diameter narrowing. A substantial overestimation of the traditional Doppler pressure drop measurement is observed, especially for severe stenosis levels. The simulation results indicate that elasticity of the aortic wall has a relatively little effect on trans-stenotic pressure drop for the range of mild to moderate SVAS cases, but predicted to have a profound effect for severe SVAS cases. Moreover, significant sensitivity to the pressure drop across the SVAS region from stenosis severity is observed.

On the relevance of boundary conditions and viscosity models in blood flow simulations in patient-specific aorto-coronary bypass models

Physiologically realistic results are the aim of every blood flow simulation. This is not different in aorto-coronary bypasses where the properties of the coronary circulation may significantly affect the relevance of the performed simulations. By considering three patient-specific bypass geometries, the present article focuses on two aspects of the coronary blood flow - its phasic flow pattern and its behaviour affected by blood rheology. For the phasic flow property, a multiscale modelling approach is chosen as a means to assess the ability of five different types of coronary boundary conditions (mean arterial pressure, Windkessel model and three lumped parameter models) to attain realistic coronary haemodynamics. From the analysed variants of boundary conditions, the best option in terms of physiological characteristics and its potential for use in patient-based simulations, is utilised to account for the effect of shear-dependent viscosity on the resulting haemodynamics and wall shear stress stimulation. Aside from the Newtonian model, the blood rheology is approximated by two non-Newtonian models in order to determine whether the choice of a viscosity model is important in simulations involving coronary circulation. A comprehensive analysis of obtained results demonstrated notable superiority of all lumped parameter models, especially in comparison to the constant outlet pressure, which regardless of bypass type gave overestimated and physiologically misleading results. In terms of rheology, it was noted that blood in undamaged coronary arteries behaves as a Newtonian fluid, whereas in vessels with atypical lumen geometry, such as that of anastomosis or stenosis, its shear-thinning behaviour should not be ignored.

Numerical study of hemodynamics in a complete coronary bypass with venous and arterial grafts and different degrees of stenosis

Cardiovascular diseases are among the leading causes of death in the world. The coronary blockage is one of most common types of these diseases that in the majority of cases has been treated by bypass surgery. In the bypass surgery, a graft is implemented to alter the blocked coronary and allow the blood supply process. The hemodynamic characteristics of the bypass strongly depend on the geometry and mechanical properties of the graft. In the present study, the fluid-structure interaction (FSI) analysis is conducted to investigate the bypass performance for a thoracic artery as well as a saphenous vein graft. Blood flow introduces a pressure on the walls of the graft which behaves as a hyperelastic material. A complete coronary bypass with stenosis degrees of 70% and 100% is modeled. To consider the nonlinear stress-strain behavior of the grafts, a five parameter Mooney-Rivlin hyperplastic model is implemented for the structural analysis and blood is assumed to behave as a Newtonian fluid. The simulations are performed for a structured grid to solve the governing equations using finite element method (FEM). The results show that wall shear stress (WSS) for saphenous vein is larger than that of thoracic artery while the total deformation of the thoracic artery is larger compared to the saphenous vein. Also, for the venous grafts or lower stenosis degree, the oscillatory shear index (OSI) is higher at both left and right anastomoses meaning that venous grafts as well as lower degree of stenosis are more critical in terms of restenosis.

Assessment of surface roughness and blood rheology on local coronary haemodynamics: a multi-scale computational fluid dynamics study

The surface roughness of the coronary artery is associated with the onset of atherosclerosis. The study applies, for the first time, the micro-scale variation of the artery surface to a 3D coronary model, investigating the impact on haemodynamic parameters which are indicators for atherosclerosis. The surface roughness of porcine coronary arteries have been detailed based on optical microscopy and implemented into a cylindrical section of coronary artery. Several approaches to rheology are compared to determine the benefits/limitations of both single and multiphase models for multi-scale geometry. Haemodynamic parameters averaged over the rough/smooth sections are similar; however, the rough surface experiences a much wider range, with maximum wall shear stress greater than 6 Pa compared to the approximately 3 Pa on the smooth segment. This suggests the smooth-walled assumption may neglect important near-wall haemodynamics. While rheological models lack sufficient definition to truly encompass the micro-scale effects occurring over the rough surface, single-phase models (Newtonian and non-Newtonian) provide numerically stable and comparable results to other coronary simulations. Multiphase models allow for phase interactions between plasma and red blood cells which is more suited to such multi-scale models. These models require additional physical laws to govern advection/aggregation of particulates in the near-wall region.

Blood residence time to assess significance of coronary artery stenosis

Coronary artery stenosis is a narrowing of coronary lumen space caused by an atherosclerotic lesion. Fractional flow reserve (FFR) is the gold standard metric to assess physiological significance of coronary stenosis, but requires an invasive procedure. Computational modeling in conjunction with patient-specific imaging demonstrates formation of regions of recirculatory flow distal to a stenosis, increasing mean blood residence time relative to uninhibited flow. A new computational parameter, mean blood residence time (BloodRT), was computed for 100 coronary artery segments for which FFR was known. A threshold for BloodRT was determined to assess the physiological significance of a stenosis, analogous to diagnostic threshold for FFR. Model sensitivity and specificity of BloodRT for diagnosis of hemodynamically significant coronary stenosis was 98% and 96% respectively, compared with FFR. When applied to clinical practice, this could potentially allow practicing cardiologists to accurately assess the severity of coronary stenosis without resorting to invasive techniques.

Myocardial Fluid Balance and Pathophysiology of Myocardial Edema in Coronary Artery Bypass Grafting

Myocardial edema is one of the most common complications of coronary artery bypass grafting (CABG) that is linearly related to many coronary artery diseases. Myocardial edema can cause several consequences including systolic dysfunction, diastolic dysfunction, arrhythmia, and cardiac tissue fibrosis that can increase mortality in CABG. Understanding myocardial fluid balance and tissue and systemic fluid regulation is crucial in order to ultimately link how coronary artery bypass grafting can cause myocardial edema in such a setting. The identification of susceptible patients by using imaging modalities is still challenging. Future studies about the technique of imaging modalities, examination protocols, prevention, and treatment of myocardial edema should be carried out, in order to limit myocardial edema occurrence and prevent complications.

The effect of anastomotic angle and diameter ratio on flow field in the distal end-to-side anastomosis

Flow fields in the distal end-to-side anastomosis of coronary artery bypass graft are associated with intimal hyperplasia and bypass failure. This work aims to demonstrate the effect of anastomotic angle and diameter ratio on flow field of coronary artery bypass graft. The flow fields inside polydimethylsiloxane models of coronary artery bypass graft with various anastomotic angles (α = 30°, 45°, 60° and 75°) and different diameter ratios (Φ = 0.78 and 1.11) are investigated using particle image velocimetry and computational fluid dynamics method under pulsatile flow condition. The results show that the anastomotic angle is positively correlated with the number and area of the recirculation zone, and the flow field disturbance at the anastomosis will develop in the same direction. Compared with that of Φ = 0.78, when Φ = 1.11, the flow fields at the anastomosis are relatively smoother with less turbulence.

NUMERICAL INVESTIGATION OF THE EFFECTS OF BED SHAPE ON THE END-TO-SIDE CABG HEMODYNAMICS

Disrupted flow initiates and aggravates intimal thickening in the end-to-side (ETS) coronary artery bypass grafting (CABG), which may lead to failure. To enhance the post-intervention hemodynamics, the geometry is either optimized or totally reconfigured. Majority of configurations proposed by researchers have not suited CABG surgery, for they entailed rigorous manipulation on conventional grafts in situ, which was neither swift nor straightforward. The aim of the present study is, thus, to introduce a slight, yet effective, modification to a conventional ETS CABG configuration, and numerically investigate its effects on updated hemodynamic and structural environment, anticipating the longevity of proposed configuration and CABG success. This fairly simple modification may easily be made positioning a pre-designed anastomotic device between the bed of host artery in the conventional ETS CABG and its surrounding tissues. Conducting comprehensive numerical simulations, performance of the proposed configuration was assessed using idealized and patient-specific geometries of the conventional ETS CABG. Blood flow was simulated in a conventional and an updated CABG configuration considering 2-way fluid–structure interaction. Results revealed that, although the proposed configuration may induce higher structural stresses in vessels walls, it may improve important hemodynamic metrics such as wall shear stress gradient, oscillatory shear index, and relative residence time on host artery bed reducing disruption of flow. This study may also set the stage for design engineers and regulatory officials to evolve ETS CABG toward more hemodynamics-friendly approaches. Further in vitro, preclinical, and clinical experiments are, yet, entailed to accomplish ideal designs of procedural guidelines/grafts.

Patient-Specific, Multi-Scale Modeling of Neointimal Hyperplasia in Vein Grafts

Neointimal hyperplasia is amongst the major causes of failure of bypass grafts. The disease progression varies from patient to patient due to a range of different factors. In this paper, a mathematical model will be used to understand neointimal hyperplasia in individual patients, combining information from biological experiments and patient-specific data to analyze some aspects of the disease, particularly with regard to mechanical stimuli due to shear stresses on the vessel wall. By combining a biochemical model of cell growth and a patient-specific computational fluid dynamics analysis of blood flow in the lumen, remodeling of the blood vessel is studied by means of a novel computational framework. The framework was used to analyze two vein graft bypasses from one patient: a femoro-popliteal and a femoro-distal bypass. The remodeling of the vessel wall and analysis of the flow for each case was then compared to clinical data and discussed as a potential tool for a better understanding of the disease. Simulation results from this first computational approach showed an overall agreement on the locations of hyperplasia in these patients and demonstrated the potential of using new integrative modeling tools to understand disease progression.