Cardiology Research and Practice

On this page

Abstract Introduction Results Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Advances in Hemodynamic Analysis in Cardiovascular Diseases

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 5903790 | https://doi.org/10.1155/2020/5903790

The Hemodynamic Effect of Enhanced External Counterpulsation Treatment on Atherosclerotic Plaque in the Carotid Artery: A Framework of Patient-Specific Computational Fluid Dynamics Analysis

Jianhang Du,^1,2,3Guangyao Wu,⁴Bokai Wu,⁵Chang Liu,⁵Zhouming Mai,^1,2Yumeng Liu,⁵Yawei Wang,⁶Pandeng Zhang,⁵Guifu Wu,^1,2,3and Jia Liu^5,7

Academic Editor: Panagiotis Korantzopoulos

Received07 Jun 2019

Revised18 Sept 2019

Accepted21 Jan 2020

Published30 Apr 2020

Abstract

Long-term enhanced external counterpulsation (EECP) therapy has been recommended for antiatherogenesis in recent clinical observations and trials. However, the precise mechanism underlying the benefits has not been fully clarified. To quantify the effect of EECP intervention on arterial hemodynamic environment, a framework of numerical assessment was introduced using a parallel computing algorithm. A 3D endothelial surface of the carotid artery with mild atherosclerotic plaque was constructed from images of magnetic resonance angiography (MRA). Physiologic boundary conditions were derived from images of the ultrasound flow velocity spectrum measured at the common carotid artery and before and during EECP intervention. Hemodynamic factors relating to wall shear stress (WSS) and its spatial and temporal fluctuations were calculated and analyzed, which included AWSS, OSI, and AWSSG. Measuring and computational results showed that diastole blood pressure, perfusion, and WSS level in carotid bifurcation were significantly increased during EECP intervention. Mean AWSS level throughout the model increased by 16.9%, while OSI level did not show a significant change during EECP. We thus suggested that long-term EECP treatment might inhibit the initiation and development of atherosclerotic plaque via improving the hemodynamic environment in the carotid artery. Meanwhile, EECP performance induced a 19.6% increase in AWSSG level, and whether it would influence the endothelial functions may need a further study. Moreover, the numerical method proposed in this study was expected to be useful for the instant assessment of clinical application of EECP .

1. Introduction

As a kind of noninvasive and atraumatic assisted circulation procedure, enhanced external counterpulsation (EECP) has exhibited itself to be an effective, safe, and economical therapy in clinics for the management of ischemic cardiovascular and cerebrovascular diseases in the recent decades [1–5] and has been thought providing a better choice for patients with chronic stable angina who failed to respond to standard revascularization procedures and aggressive pharmacotherapy [6].

The treatment of EECP (see Figure 1) involves the use of an EECP device to inflate and deflate a series of compressive cuffs wrapped around the patient’s calves, lower thighs, and upper thighs. As a result, the enhanced flow perfusion is achieved from the device’s propelling blood from veins of the lower body to arteries of the upper body and increases the blood supply for the important organs and the brain [7].

Long-term EECP intervention has been demonstrated in recent studies to be able to improve the endothelial functions and in turn may inhibit the generation and development of atherosclerosis lesion [8–11]. The hemodynamic effects, especially the wall shear stress variations, induced by EECP have been thought contributing the most important part of its benefits. Michaels et al. [12] confirmed that EECP treatment could significantly increase coronary artery flow determined by both Doppler and angiographic techniques. Braith et al. [10] suggested that EECP had a beneficial effect on peripheral artery flow-mediated dilation and endothelial-derived vasoactive agents. Our previous study [8] experimentally confirmed that EECP inhibits intimal hyperplasia and atherogenesis by modifying biomechanical stress-responsive gene expression. However, the actual influence of EECP intervention on wall shear stress (WSS) and its spatial and temporary fluctuations remained elusive.

It has been widely accepted that biomechanical stresses of large and medium arteries play an important role in maintaining the functions of endothelium and vascular remodeling progression [13]. Low and/or oscillating WSS has been commonly believed to be correlated positively with initiation and development of atherosclerosis [14–16]. Several hemodynamic factors have been proposed by different research groups to represent the biomechanical indicators connected to arterial functions, such as average wall shear stress (AWSS), oscillatory shear index (OSI), particle resident time (PRT), and wall shear stress gradient (WSSG).

This paper was aimed to conduct a pilot study on how the EECP treatment affects the hemodynamic environment and the important factors in carotid arterial bifurcation where the atherosclerotic lesion localizes characteristically. A numerical method-combined finite element method with in vivo medical imaging measurement was introduced to assess the local hemodynamic factors during EECP intervention.

2. Medical Image Acquisition and Processing

A 55-year-old coronary heart disease patient with mild carotid atheromatous plaque diagnosed (severity of stenosis was less than 20%) was enrolled to the measurement. The subject underwent the clinical protocol for carotid plaque MRI on a 3T MRI (General Electric Company, Discovery MR750). A two-element bilateral 8-channel carotid surface coil (Wk401 Jiangyin Wankang Medical Technology Co., Ltd.) was used for image acquisition.

A three-dimensional phase-contrast magnetic resonance angiography (3D Phase Contrast) sequence was performed: repetition time/echo time (TR/TE) 15.0/3.7 ms; flip angle 8°, field of view 32 cm × 32 cm, slice thickness 1.8 mm, and matrix 384 × 256. A High Resolution Three-Dimensional CUBE (computer use by engineers) T1 Weighted Imaging(HR 3D CUBE T1WI) sequence was performed as follows: repetition time/echo time (TR/TE) 575/15 ms; echo chain length (ETL) 24, slice thickness 0.8 mm, field of view 28 cm × 28 cm, and matrix 256 × 256).

3. EECP Intervention Protocol and Color Doppler Ultrasound Measurement

A short-term EECP intervention was performed using Pushikang P-ECP/TM Oxygen Saturation Monitoring Enhanced External Counterpulsation Instrument (made in Chongqing, China). The subject received a single, 45-minute session EECP treatment with the working pressure set to 0.033 MPa.

The blood velocity measurements of before EECP (rest state) and during EECP (15-25min after EECP initiated) were performed based on a Color Doppler Ultrasound System (Philip EPIQ7). (see Figure 2). The left common carotid arteries (CCA) were examined with 1.5 cm proximal to the bifurcation of the vessels. The blood velocity waveforms (see Figure 3) in cardiac cycles and before and during EECP intervention were extracted from the images of ultrasound flow velocity spectrum. Meanwhile, the diameter changes of the lumen section (see Figure 4) in cardiac cycles were extracted from images of ultrasound. The blood flow rate in cardiac cycle and perfusion in CCA could be calculated based on velocity waveforms and diameter changes.

(a)

(b)

4. 3D Reconstruction for the Endothelial Surface of the Carotid Artery

We propose a method to virtually reconstruct the endothelial surface of the carotid artery so as to visualize the carotid atheromatous plaque in 3D. The pipeline of our work is shown in Figure 5, where I_t denotes the input MR image with index t, t = 1, 2, 3, …; P_t represents the artery endothelial boundary extracted from I_t. This pipeline consists of three main steps: image preprocessing, endothelial boundary extraction, and 3D reconstruction.

5. Image Preprocessing

There exists serious inherent noise in the MR image, as shown in Figure 6(a). The noise has a detrimental influence on the accuracy of the artery extraction and 3D reconstruction. In order to reduce the noise without affecting the shape of the carotid artery, we make use of the morphological technique called open-by-reconstruction and close-by-reconstruction [17]. The processed result is illustrated in Figure 6(b). It can be seen that most of the clutters in Figure 6(a) have been removed, and the shape of the carotid artery is not influenced.

(a)

(b)

5.1. Endothelial Boundary Extraction

The extraction of the endothelial boundary of the carotid artery in each MRA image is implemented based on the result of image preprocessing. For t = 1, we manually extract the boundary of the carotid artery in I₁, so as to initialize the subsequent automatic extraction. If t > 1, given the extraction result P_t−1 (the yellow curve in Figure 7(a)),we first binarize the I_t preprocessed and then obtain the edge map E from this binary image (Figure 7(b)). Because of the great similarity between I_t and I_t−1, we eventually optimize a closed curve c = [x, y], which satisfies the following equation with calculus of variations using P_t−1 as the initial value.where s is the arc parameter of c, and α, β, as well as γ are the three scalars to balance the three terms in equation (1). The optimized c is the result P_t of I_t, as shown by the yellow curve in Figure 7(c). Equation (1) can be calculated by the method proposed in [18].

(a)

(b)

(c)

5.2. Texture Flattening and 3D Reconstruction

The carotid artery can be approximately considered as a surface of revolution (SOR). Providing that the endothelial boundary of the carotid artery in each MRI slice has been extracted, the algorithm proposed in [19] can be adopted to flatten the texture of the endothelial surface of the carotid artery and then generate a 3D texture reconstruction for the endothelial surface. The carotid artery will bifurcate at its end, so we can generate the 3D texture reconstruction for each bifurcation by the same way and finally combine all reconstructions together. The final reconstruction of the carotid artery is demonstrated in Figure 8. In this figure, the dark and concave region is the carotid atherosclerotic plaque.

6. The CFD Method and the Boundary Conditions

6.1. Geometry and Boundary Conditions

To simplify simulations, the elasticity of vessel wall is not considered in the present study (i.e., computational domain was fixed). Original geometry consists of four boundaries: inlet, outlet, and wall. An artificial extension geometry, with a length of 5 times the averaged radius of the inlet is added outward along the normal direction of the inlet boundary to acquire a fully developed velocity waveform [15], as depicted in Figure 9. Inflow velocity (V_in) measured by the carotid Doppler (see Figure 3) is specified at the inlet_ex. Opening condition is set at the outlet, and wall is assumed to be no-slip.

6.2. Mesh Generation

Most part of the geometry is meshed with tetrahedral cells by employing the commercial software ANSYS ICEM (ANSYS, Inc., USA). To capture the flow behavior where high velocity gradient exists, inflation layers are created near the wall [16] (as shown in Figure 10). To optimize the mesh size, a specific mesh-independent study is carried out for reliable results, while keeping computational loads as low as possible. As indicated in Table 1, change in AWSS is around 3% with refinement from Mesh 1 to Mesh 3, while less than 1% from Mesh 3 to Mesh 4 for both states with and without EECP. Mesh used in this study is of quantity 720085 (i.e., Mesh 3) and with a quality of ∼0.4 (measured by its orthogonality and warpage), which is at an acceptable level.

6.3. Rheology and Governing Equations

The blood fluid used in the present study is assumed to be impressible and isoviscous (i.e., Newtonian type). Therefore, the governing transport equations for this study are continuity and momentum equations which can be written in their general forms [20], as follows:where is the velocity vector, is the pressure, is the fluid density, and is the external force (assumed to be 0 here).

6.4. The solver

The governing equations are discretized by CFX code based on a finite-volume method. To meet the requirement on both robustness and accuracy, a so-called “High Resolution Advection Scheme” is implemented in this study [21]. Numerical solutions are acquired while root mean square (RMS) of both mass and momentum residuals are below 10⁻⁵. In fact, however, at most time steps, even lower RMS residual values are generally reached. We solve the unknowns with this configuration for four cardiac cycles. Results of the last cardiac cycle are presented in the following section.

7. Results

Several important hemodynamic factors such as AWSS, OSI, RRT, and WSSG were calculated in this paper, which were introduced by different research groups to represent the WSS level in a cardiac cycle and its spatial and temporary fluctuation. These factors are defined as follows [22, 23]:where is the magnitude of the instantaneous WSS vector and T is the cardiac cycle.

The numerical results are shown in Figures 11–15. Focus of previous studies were mainly put on how EECP affected the blood flow and pressure in diastole [1, 12]. In the current paper, the velocity and pressure at time points of t = 0.54T and t = 0.58T were chosen to represent the blood velocity and pressure in diastole before and during EECP intervention respectively, considering that peak blood velocity in diastole occurred at these time points based on Figure 3.

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

8. Discussions

The computational results are summarized in Table 2, which include AWSS, OSI, and AWSSG over the cardiac cycle and before and during EECP intervention. All calculations and statistics were performed based on the whole model.

The results showed that EECP performance significantly increased the blood velocity in diastole, as well as the blood pressure. Peak relative pressure during EECP increased about 2260 Pa comparing to pre-EECP state. The elevation of diastolic pressure was thought playing a key role for the increasement of perfusion [24]. A calculation based on Figures 3 and 4 showed that EECP performance induced a 12.6% elevation of the perfusion over a cardiac cycle in CCA.

WSS has been widely recognized to be an important hemodynamic factor affecting the initiation and development of atherosclerotic plaque. It is now well accepted that low and oscillating WSS correlate positively with atherosclerosis progression [25]. Our calculating results showed that EECP performance significantly increased the WSS level in carotid bifurcation and especially in bulb and ICA, and the mean AWSS level throughout the model increased by 16.9% (7.90 Pa versus 6.76 Pa).

As a factor proposed to represent the temporary fluctuation of WSS, OSI has been found to be associated with early atherosclerosis in some studies, and region of high OSI coincides with a high probability of occurrence of early atherosclerosis lesions [26, 27]. Our results showed that EECP performance didn’t induce a significant change in OSI level in the carotid artery, although this kind of intervention greatly changed the blood flow pattern.

As a factor proposed to represent the spatial fluctuation of WSS, WSSG has been suggested in some studies that might correlate with intima-medial thickness and endothelial dysfunction [27, 28]. Our current study showed that EECP performance induced a significant increase in WSSG level in the carotid artery. The mean and peak AWSSG level throughout the model, respectively, increased by 19.6% (2.01 × 10³ Pa/m versus 1.68 × 10³ Pa/m) and 22.6% (2.60 × 10⁴ Pa/m versus 2.12 × 10⁴ Pa/m).

One of the main limitations in the current study was that we did not enroll healthy subjects as comparison. Because the aim of this paper was to introduce a medical imaging-based numerical method to assess the instant hemodynamic response during EECP treatment and to conduct a pilot study of the influence of EECP on WSS and its fluctuations in carotid bifurcation with mild plaque.

9. Conclusions

We suggest that the framework of patient-specific numerical approach developed in the current paper can be potentially used in clinics for the assessment of instant hemodynamic response in the carotid artery during EECP treatment and in turn may play a role on improvement of the treatment strategies for better clinical outcome. Meanwhile, findings of this paper show that EECP treatment induced a significant augmentation of blood perfusion and WSS level in the carotid artery, which may be the main hemodynamic mechanism underlying its good clinical effect for treatment of the ischemic cerebrovascular diseases and the long-term effect for inhibition of the atherosclerosis lesion.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Jianhang Du and Guangyao Wu contributed equally to this work.

Acknowledgments

This work was supported by the National Key R&D Program of the Ministry of Science and Technology of China (Grant no. 2016YFC1301602), National Natural Science Foundation of China (Grant no. 81661168015), and Shenzhen Science and Technology Innovation Commission (Grant nos. JCYJ20160608153506088, YJ20180306174831458, and ZDSYS201703031711426).

References

W. Lin, L. Xiong, J. Han et al., “External counterpulsation augments blood pressure and cerebral flow velocities in ischemic stroke patients with cerebral intracranial large artery occlusive disease,” Stroke, vol. 43, no. 11, pp. 3007–3011, 2012.
View at: Publisher Site | Google Scholar
W. E. Lawson, J. C. K. Hui, E. D. Kennard, and G. Linnemeier, “Enhanced external counterpulsation is cost-effective in reducing hospital costs in refractory angina patients,” Clinical Cardiology, vol. 38, no. 6, pp. 344–349, 2015.
View at: Publisher Site | Google Scholar
Task Force Members, G. Montalescot, U. Sechtem et al., “2013 ESC guidelines on the management of stable coronary artery disease,” European Heart Journal, vol. 34, no. 38, pp. 2949–3003, 2013.
View at: Publisher Site | Google Scholar
E. C. Jauch, J. L. Saver, H. P. Adams et al., “Guidelines for the early management of patients with acute ischemic stroke,” Stroke, vol. 44, no. 3, pp. 870–947, 2013.
View at: Publisher Site | Google Scholar
S. D. Fihn, J. C. Blankenship, K. P. Alexander et al., “2014 ACC/AHA/AATS/PCNA/SCAI/STS focused update of the guideline for the diagnosis and management of patients with stable ischemic heart disease,” The Journal of Thoracic and Cardiovascular Surgery, vol. 149, no. 3, pp. e5–e23, 2015.
View at: Publisher Site | Google Scholar
A. Raza, K. Steinberg, J. Tartaglia et al., “Enhanced external counterpulsation therapy: past, present, and the future,” Cardiology in Review, vol. 25, no. 2, p. 59, 2016.
View at: Google Scholar
J. Du and L. Wang, “Enhanced external counterpulsation treatment may intervene the advanced atherosclerotic plaque progression by inducing the variations of mechanical factors: a 3D FSI study based on in vivo animal experiment,” Molecular & Cellular Biomechanics, vol. 12, no. 4, pp. 249–263, 2015.
View at: Google Scholar
Y. Zhang, X. He, X. Chen et al., “Enhanced external counterpulsation inhibits intimal hyperplasia by modifying shear stress-responsive gene expression in hypercholesterolemic pigs,” Circulation, vol. 116, no. 5, pp. 526–534, 2007.
View at: Publisher Site | Google Scholar
D. P. Casey, C. R. Conti, W. W. Nichols, C. Y. Choi, M. A. Khuddus, and R. W. Braith, “Effect of enhanced external counterpulsation on inflammatory cytokines and adhesion molecules in patients with angina pectoris and angiographic coronary artery disease,” The American Journal of Cardiology, vol. 101, no. 3, pp. 300–302, 2008.
View at: Publisher Site | Google Scholar
R. W. Braith, C. R. Conti, W. W. Nichols et al., “Enhanced external counterpulsation improves peripheral artery flow-mediated dilation in patients with chronic angina,” Circulation, vol. 122, no. 16, pp. 1612–1620, 2010.
View at: Publisher Site | Google Scholar
D.-Y. Yang and G.-F. Wu, “Vasculoprotective properties of enhanced external counterpulsation for coronary artery disease: beyond the hemodynamics,” International Journal of Cardiology, vol. 166, no. 1, pp. 38–43, 2013.
View at: Publisher Site | Google Scholar
A. D. Michaels, M. Accad, T. A. Ports, and W. Grossman, “Left ventricular systolic unloading and augmentation of intracoronary pressure and Doppler flow during enhanced external counterpulsation,” Circulation, vol. 106, no. 10, pp. 1237–1242, 2002.
View at: Publisher Site | Google Scholar
A. J. Brown, Z. Teng, P. C. Evans, J. H. Gillard, H. Samady, and M. R. Bennett, “Role of biomechanical forces in the natural history of coronary atherosclerosis,” Nature Reviews Cardiology, vol. 13, no. 4, pp. 210–220, 2016.
View at: Publisher Site | Google Scholar
Y. Mohamied, E. M. Rowland, E. L. Bailey, S. J. Sherwin, M. A. Schwartz, and P. D. Weinberg, “Change of direction in the biomechanics of atherosclerosis,” Annals of Biomedical Engineering, vol. 43, no. 1, pp. 16–25, 2015.
View at: Publisher Site | Google Scholar
X. Li, X. Liu, X. Li, L. Xu, X. Chen, and F. Liang, “Tortuosity of the superficial femoral artery and its influence on blood flow patterns and risk of atherosclerosis,” Biomechanics and Modeling in Mechanobiology, vol. 18, no. 4, pp. 883–896, 2019.
View at: Publisher Site | Google Scholar
L. Xu, F. Liang, B. Zhao, J. Wan, and H. Liu, “Influence of aging-induced flow waveform variation on hemodynamics in aneurysms present at the internal carotid artery: a computational model-based study,” Computers in Biology and Medicine, vol. 101, pp. 51–60, 2018.
View at: Publisher Site | Google Scholar
R. Gonzalez, R. Woods, and S. Eddins, Digital Image Processing Using Matlab, Publishing House of Electronics Industry, Beijing, China, 2004.
P. Rosin, Y. Lai, C. Liu et al., “Virtual recovery of content from X-ray micro-tomography scans of damaged historic scrolls,” Scientific Reports, vol. 8, no. 1, Article ID 11901, 2018.
View at: Publisher Site | Google Scholar
C. Liu and W. Hu, “Real-time geometric fitting and pose estimation for surface of revolution,” Pattern Recognition, vol. 85, pp. 90–108, 2019.
View at: Publisher Site | Google Scholar
R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids: Fluid Mechanics, Wiley, New York, USA, 1987.
S. Tian, “Wall effects for spherical particle in confined shear-thickening fluids,” Journal of Non-Newtonian Fluid Mechanics, vol. 257, pp. 13–21, 2018.
View at: Publisher Site | Google Scholar
F. Rikhtegar, J. A. Knight, U. Olgac et al., “Choosing the optimal wall shear parameter for the prediction of plaque location—a patient-specific computational study in human left coronary arteries,” Atherosclerosis, vol. 221, no. 2, pp. 432–437, 2012.
View at: Publisher Site | Google Scholar
C. V. Cunnane, E. M. Cunnane, and M. T. Walsh, “A review of the hemodynamic factors believed to contribute to vascular access dysfunction,” Cardiovascular Engineering and Technology, vol. 8, no. 3, pp. 280–294, 2017.
View at: Publisher Site | Google Scholar
H. Mueller, S. M. Ayres, and W. J. Grace, “Hemodynamic and myocardial metabolic response to external counterpulsation in acute myocardial infarction in man,” The American Journal of Cardiology, vol. 31, no. 1, p. 149, 1973.
View at: Publisher Site | Google Scholar
D. Tang, R. D. Kamm, C. Zheng et al., “Image-based modeling for better understanding and assessment of atherosclerotic plaque progression and vulnerability: data, modeling, validation, uncertainty and predictions,” Journal of Biomechanics, vol. 47, no. 4, pp. 834–846, 2014.
View at: Publisher Site | Google Scholar
D. N. Ku, C. K. D. P. Giddens, and S. Glagov, “Pulsatile flow and atherosclerosis in the human carotid bifurcation. Positive correlation between plaque location and low oscillating shear stress,” Arteriosclerosis: An Official Journal of the American Heart Association, Inc, , vol. 5, no. 3, pp. 293–302, 1985.
View at: Publisher Site | Google Scholar
J. R. Buchanan, C. Kleinstreuer, G. A. Truskey, and M. Lei, “Relation between non-uniform hemodynamics and sites of altered permeability and lesion growth at the rabbit aorto-celiac junction,” Atherosclerosis, vol. 143, no. 1, pp. 27–40, 1999.
View at: Publisher Site | Google Scholar
J. Murphy and F. Boyle, “Predicting neointimal hyperplasia in stented arteries using time-dependant computational fluid dynamics: a review,” Computers in Biology and Medicine, vol. 40, no. 4, pp. 408–418, 2010.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Jianhang Du et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1709

Downloads

1364

Citations