3D reconstruction techniques of human coronary bifurcations for shear stress computations

Background

Heterogeneity in plaque composition in human coronary artery bifurcations is associated with blood flow induced shear stress. Shear stress is generally determined by combing 3D lumen data and computational fluid dynamics (CFD). We investigated two new procedures to generate 3D lumen reconstructions of coronary artery bifurcations for shear stress computations.

Methods

We imaged 10 patients with multislice computer tomography (MSCT) and intravascular ultrasound (IVUS). The 3D reconstruction of the main branch was based on the fusion of MSCT and IVUS. The proximal part of side branch was reconstructed using IVUS data or MSCT data, resulting in two different reconstructions of the bifurcation region. The distal part of the side branch was based on MSCT data alone. The reconstructed lumen was combined with CFD to determine the shear stress. Low and high shear stress regions were defined and shear stress patterns in the bifurcation regions were investigated.

Results

The 3D coronary bifurcations were successfully generated with both reconstruction procedures. The geometrical features of the bifurcation region for the two reconstruction procedures did not reveal appreciable differences. The shear stress maps showed a qualitative agreement, and the low and high shear stress regions were similar in size and average shear stress values were identical. The low and high shear stress regions showed an overlap of approximately 75%.

Conclusion

Reconstruction of the side branch with MSCT data alone is an adequate technique to study shear stress and wall thickness in the bifurcation region. The reconstruction procedure can be applied to further investigate the effect of shear stress on atherosclerosis in coronary bifurcations.     

Coronary stent implantation changes 3-D vessel geometry and 3-D shear stress distribution

Mechanisms of in-stent restenosis are not fully understood. Shear stress is known to play a role in plaque and thrombus formation and is sensitive to changes in regional vessel geometry. Hence, we evaluated the regional changes in 3-D geometry and shear stress induced by stent placement in coronary arteries of pigs.Methods. 3-D reconstruction was performed, applying a combined angiographic and IVUS technique (ANGUS), from seven Wallstents (diameter 3.5 (n=3) and 5mm (n=4)), which were implanted in seven coronary arteries of five pigs. This 3-D geometry was used to calculate locally the curvature, while the shear stress distribution was obtained by computational fluid dynamics. Local changes in shear stress were obtained at the entrance and exit of the stent for baseline (0. 65+/-0.22 ml/s) and hyperemic flow (2.60+/-0.86 ml/s) conditions. Results. After stent implantation, the curvature increased by 121% at the entrance and by 100% at the exit of the stent, resulting in local changes in shear stress. In general, at the entrance of the stent local maxima in shear stress were generated, while at the exit both local maxima and minima in shear stress were observed (p<0.05). Additionally, the shear stress at the entrance and exit of the stent were correlated with the local curvature (r: 0.30-0.84).Conclusion. Stent implantation changes 3-D vessel geometry in such a way that regions with decreased and increased shear stress occur close to the stent edges. These changes might be related to the asymmetric patterns of in-stent restenosis.     

Strain distribution over plaques in human coronary arteries relates to shear stress

Once plaques intrude into the lumen, the shear stress they are exposed to alters with hitherto unknown consequences for plaque composition. We investigated the relationship between shear stress and strain, a marker for plaque composition, in human coronary arteries. We imaged 31 plaques in coronary arteries with angiography and intravascular ultrasound. Computational fluid dynamics was used to obtain shear stress. Palpography was applied to measure strain. Each plaque was divided into four regions: upstream, throat, shoulder, and downstream. Average shear stress and strain were determined in each region. Shear stress in the upstream, shoulder, throat, and downstream region was 2.55+/-0.89, 2.07+/-0.98, 2.32+/-1.11, and 0.67+/-0.35 Pa, respectively. Shear stress in the downstream region was significantly lower. Strain in the downstream region was also significantly lower than the values in the other regions (0.23+/-0.08% vs. 0.48+/-0.15%, 0.43+/-0.17%, and 0.47+/-0.12%, for the upstream, shoulder, and throat regions, respectively). Pooling all regions, dividing shear stress per plaque into tertiles, and computing average strain showed a positive correlation; for low, medium, and high shear stress, strain was 0.23+/-0.10%, 0.40+/-0.15%, and 0.60+/-0.18%, respectively. Low strain colocalizes with low shear stress downstream of plaques. Higher strain can be found in all other plaque regions, with the highest strain found in regions exposed to the highest shear stresses. This indicates that high shear stress might destabilize plaques, which could lead to plaque rupture.     

Velocity and wall shear stress patterns in the human right coronary artery.

Blood flow dynamics in the human right coronary artery have not been adequately quantified despite the clinical significance of coronary atherosclerosis. In this study, a technique was developed to construct a rigid flow model from a cast of a human right coronary artery. A laser photochromic method was used to characterize the velocity and wall shear stress patterns. The flow conditions include steady flow at Reynolds numbers of 500 and 1000 as well as unsteady flow with Womersley parameter and peak Reynolds number of 1.82 and 750, respectively. Characterization of the three-dimensional geometry of the artery revealed that the largest spatial variation in curvature occurred within the almost branch-free proximal region, with the greatest curvature existing along the acute margin of the heart. In the proximal segment, high shear stresses were observed on the outer wall and lower, but not negative, stresses along the inner wall. Low shear stress on the inner wall may be related to the preferential localization of atherosclerosis in the proximal segment of the right coronary artery. However, it is possible that the large difference between the outer and inner wall shear stresses may also be involved.     

Breaking symmetry in non-planar bifurcations: distribution of flow and wall shear stress

Non-planarity in blood vessels is known to influence arterial flows and wall shear stress. To gain insight, computational fluid dynamics (CFD) has been used to investigate effects of curvature and out-of-plane geometry on the distribution of fluid flows and wall shear stresses in a hypothetical non-planar bifurcation. Three-dimensional Navier-Stokes equations for a steady state Newtonian fluid were solved numerically using a finite element method. Non-planarity in one of the two daughter vessels is found to deflect flow from the inner wall of the vessel to the outer wall and to cause changes in the distribution of wall shear stresses. Results from this study agree to experimental observations and CFD simulations in the literature, and support the view that non-planarity in blood vessels is a factor with important haemodynamic significance and may play a key role in vascular biology and pathophysiology.     

Variations in geometry and shear rate distribution in casts of human aortic bifurcations

Shear rates measured at sites along the walls of ten casts of human aortic bifurcations are presented to provide an indication of the variability in the shear rate distribution that can occur among individuals. The geometric characteristics and flow parameters for each of the casts are also given. Much of the variability in wall shear distribution seen here results from variations in arterial geometry; variations in flow parameters can widen further the normal range of shear distributions.     

Wall shear stress in normal left coronary artery tree

Despite the fact that the role of wall shear stress (WSS) as a local mechanical factor in atherogenesis is well established, its distribution over the entire normal human left coronary artery (LCA) tree has not yet been studied. A three-dimensional computer generated model of the epicardial LCA tree, based on averaged human data set extracted from angiographies, was adopted for finite-element analysis of the Navier-Stokes flow equations treating blood as non-Newtonian fluid. The LCA tree includes the left main coronary artery (LMCA), the left anterior descending (LAD), the left circumflex artery (LCxA) and their major branches. In proximal LCA tree regions where atherosclerosis frequently occurs, low WSS appears. Low WSS regions occur at bifurcations in regions opposite the flow dividers, which are anatomic sites predisposed for atherosclerotic development. On the LMCA bifurcation, at regions opposite to the flow divider, dominant low WSS values occur ranging from 0.75 to 2.25 N/m2. High WSS values are encountered at all flow dividers. This work determines, probably for the first time, the topography of the WSS in the entire normal human LCA epicardial tree. It is also the first work determining the spatial WSS differentiation between proximal and distal normal human LCA parts. The haemodynamic analysis of the entire epicardial LCA tree further verifies the implications of the WSS in atherosclerosis mechanisms.     

Wall shear stress estimates in coronary artery constrictions.

Wall shear stress estimates from laminar boundary layer theory were found to agree fairly well with the magnitude of shear stress levels along coronary artery constrictions obtained from solutions of the Navier Stokes equations for both steady and pulsatile flow. The relatively simple method can be used for in vivo estimates of wall shear stress in constrictions by using a vessel shape function determined from a coronary angiogram, along with a knowledge of the flow rate.     

Patterns of living beta-actin movement in wounded human coronary artery endothelial cells exposed to shear stress.

We previously demonstrated that physiologic levels of shear stress enhance endothelial repair. Cell spreading and migration, but not proliferation, were the major mechanisms accounting for the increases in wound closure rate (Albuquerque et al., 2000, Am. J. Physiol. Heart Circ. Physiol. 279, H293-H302). However, the patterns and movements of beta-actin filaments responsible for cell motility and translocation in human coronary artery endothelial cells (HCAECs) have not been previously investigated under physiologic flow. HCAECs transfected with beta-actin-GFP were cultured on type I collagen-coated coverslips. Confluent cell monolayers were subjected to laminar shear stress of 12 dynes/cm(2) for 18 h in a parallel-plate flow chamber to attain cellular alignment and then wounded by scraping with a metal spatula and subsequently exposed to a laminar shear stress of 20 dynes/cm(2) (S-W-sH) or static (S-W-sT) conditions. Time-lapse imaging and deconvolution microscopy was performed during the first 3 h after imposition of S-W-sH or S-W-sT conditions. The spatial and temporal dynamics of beta-actin-GFP motility and translocation during wound closure in HCAEC monolayers were analyzed under both conditions. Compared with HCAEC under S-W-sT conditions, our data show that HCAEC under S-W-sH conditions demonstrated greater beta-actin-GFP motility, filament and clumping patterns, and filament arcs used during cellular attachment and detachment. These findings demonstrate intriguing patterns of beta-actin organization and movement during wound closure in HCAEC exposed to physiological flow.     

Advanced human carotid plaque progression correlates positively with flow shear stress using follow-up scan data: An in vivo MRI multi-patient 3D FSI study

Although it has been well-accepted that atherosclerosis initiation and early progression correlate negatively with flow wall shear stresses (FSS), increasing evidence suggests mechanisms governing advanced plaque progression are not well understood. Fourteen patients were scanned 2–4 times at 18 month intervals using a histologically validated multi-contrast magnetic resonance imaging (MRI) protocol to acquire carotid plaque progression data. Thirty-two scan pairs (baseline and follow-up scans) were formed with slices matched for model construction and analysis. 3D fluid–structure interaction (FSI) models were constructed and plaque wall stress (PWS) and flow shear stress (FSS) were obtained from all matching lumen data points (400–1000 per plaque; 100 points per matched slice) to quantify correlations with plaque progression measured by vessel wall thickness increase (WTI). Using FSS and PWS data from follow-up scan, 21 out of 32 scan pairs showed a significant positive correlation between WTI and FSS (positive/negative/no significance ratio=21/8/3), and 26 out of 32 scan pairs showed a significant negative correlation between WTI and PWS (positive/negative/no significance ratio=2/26/4). The mean FSS value of lipid core nodes (n=5294) from all 47 plaque models was 63.5 dyn/cm², which was 45% higher than that from all normal vessel nodes (n=27553, p<0.00001). The results from this intensive FSI study indicate that flow shear stress from follow-up scan correlates positively with advanced plaque progression which is different from what has been observed in plaque initiation and early-stage progression. It should be noted that the correlation results do not automatically lead to any causality conclusions.     

Geometry guided data averaging enables the interpretation of shear stress related plaque development in human coronary arteries

The average low shear stress (SS) is known to determine predilection sites of atherosclerotic plaques. However, as plaques encroach into the lumen and thereby increase SS, interpretation of patient-specific data obtained at one moment in time regarding the influence of SS in the generation of atherosclerosis is not straightforward. This study aims to compare two methods of data analysis for the aid of data interpretation: (a) point-wise analysis of the raw data, (b) global analysis: to assess the history related natural SS distribution in coronary arteries by averaging the data in the axial vessel direction. Normal to mildly diseased human coronary arteries were investigated applying a combination of 3-D reconstruction technique and computational fluid dynamics (CFD). Point-wise analysis relating local wall thickness to local SS showed in only 4% of the cases an inverse relationship. In contrast, averaging the data in the axial vessel direction, showed in 38% a significant inverse relation between wall thickness and SS, resulting in an average negative slope of -0.70+/-0.46 mm/Pa. These data suggest that using a geometry guided way of data averaging may reveal history related effects of SS, which in part explains localization of atherosclerotic plaques.     

A study on the compliance of a right coronary artery and its impact on wall shear stress

A computational model incorporating physiological motion and uniform transient wall deformation of a branchless right coronary artery (RCA) was developed to assess the influence of artery compliance on wall shear stress (WSS). Arterial geometry and deformation were derived from modern medical imaging techniques, whereas the blood flow was solved numerically employing a moving-grid approach using a well-validated in-house finite element code. The simulation results indicate that artery compliance affects the WSS in the RCA heterogeneously, with the distal region mostly experiencing these effects. Under physiological inflow conditions, coronary compliance contributed to phase changes in the WSS time history, without affecting the temporal gradient of the local WSS nor the bounds of the WSS magnitude. Compliance does not cause considerable changes to the topology of WSS vector patterns nor to the localization of WSS minima along the RCA. We conclude that compliance is not an important factor affecting local hemodynamics in the proximal region of the RCA while the influence of compliance in the distal region needs to be evaluated in conjunction with the outflow to the myocardium through the major branches of the RCA.     

A new technique of coronary artery ligation: Experimental myocardial infarction in rats in vivo with reduced mortality

In vivo models of myocardial infarction following coronary artery ligation in the rat still suffer from high early mortality and a low rate of success of myocardial infarction. This study investigated the possibility of reducing early mortality and increasing the rate of myocardial infarction by modifications of surgical techniques. Eighteen rats were divided into two groups: normal control (3 rats) and ligation (15 rats). The major modifications of surgical techniques used in this study include: (1) no exteriorization of the heart, (2) ligation of the origins of the branches rather than the main trunk of the left coronary artery, (3) removal of air from the chest after closure, (4) supplying oxygen immediately after extubation. Following surgery, the rats recovered uneventfully and 11 rats were alive after 16 weeks. One rat, with a large myocardial infarction, died 2 h after surgery. Early mortality (during surgery and 1 week after surgery) was 6.7% with a success rate of myocardial infarction of 85%. The left ventricle in the ligation group showed significant dilation relative to normal and shamoperated control hearts (317% of control hearts, p < 0.001). However, myocardial mass did not increase. The average infarct size was 33%. These results demonstrate that a reduction in early mortality and an increased success rate of myocardial infarction can be achieved by modifications of surgical techniques.     

A novel workflow combining plaque imaging,plaque and plasma proteomics identifies biomarkers of human coronary atherosclerotic plaque disruption

Background

Atherosclerotic plaque rupture is the culprit event which underpins most acute vascular syndromes such as acute myocardial infarction. Novel biomarkers of plaque rupture could improve biological understanding and clinical management of patients presenting with possible acute vascular syndromes but such biomarker(s) remain elusive. Investigation of biomarkers in the context of de novo plaque rupture in humans is confounded by the inability to attribute the plaque rupture as the source of biomarker release, as plaque ruptures are typically associated with prompt down-stream events of myocardial necrosis and systemic inflammation.

Methods

We developed a novel approach to identify potential biomarkers of plaque rupture by integrating plaque imaging, using optical coherence tomography, with both plaque and plasma proteomic analysis in a human model of angioplasty-induced plaque disruption.

Results

We compared two pairs of coronary plaque debris, captured by a FilterWire Device, and their corresponding control samples and found matrix metalloproteinase 9 (MMP9) to be significantly enriched in plaque. Plaque contents, as defined by optical coherence tomography, affect the systemic changes of MMP9. Disruption of lipid-rich plaque led to prompt elevation of plasma MMP9, whereas disruption of non-lipid-rich plaque resulted in delayed elevation of plasma MMP9. Systemic MMP9 elevation is independent of the associated myocardial necrosis and systemic inflammation (measured by Troponin I and C-reactive protein, respectively). This information guided the selection of a subset of subjects of for further label free proteomics analysis by liquid chromatography tandem mass spectrometry (LC–MS/MS). We discovered five novel, plaque-enriched proteins (lipopolysaccharide binding protein, Annexin A5, eukaryotic translocation initiation factor, syntaxin 11, cytochrome B5 reductase 3) to be significantly elevated in systemic circulation at 5 min after plaque disruption.

Conclusion

This novel approach for biomarker discovery in human coronary artery plaque disruption can identify new biomarkers related to human coronary artery plaque composition and disruption.

     

A negative correlation between human carotid atherosclerotic plaque progression and plaque wall stress: in vivo MRI-based 2D/3D FSI models

It is well accepted that atherosclerosis initiation and progression correlate positively with low and oscillating flow wall shear stresses (FSS). However, this mechanism cannot explain why advanced plaques continue to grow under elevated FSS conditions. In vivo magnetic resonance imaging (MRI)-based 2D/3D multi-component models with fluid-structure interactions (FSI, 3D only) for human carotid atherosclerotic plaques were introduced to quantify correlations between plaque progression measured by wall thickness increase (WTI) and plaque wall (structure) stress (PWS) conditions. A histologically validated multi-contrast MRI protocol was used to acquire multi-year in vivo MRI images. Our results using 2D models (200-700 data points/patient) indicated that 18 out of 21 patients studied showed significant negative correlation between WTI and PWS at time 2 (T2). The 95% confidence interval for the Pearson correlation coefficient is (-0.443,-0.246), p<0.0001. Our 3D FSI model supported the 2D correlation results and further indicated that combining both plaque structure stress and flow shear stress gave better approximation results (PWS, T2: R(2)=0.279; FSS, T1: R(2)=0.276; combining both: R(2)=0.637). These pilot studies suggest that both lower PWS and lower FSS may contribute to continued plaque progression and should be taken into consideration in future investigations of diseases related to atherosclerosis.     

Wall shear stress distribution in a model human aortic arch: assessment by an electrochemical technique

Wall shear stress (WSS) distribution in a human aortic arch model is studied using 130 cathode electrodes flush-mounted on the model walls. Flow visualizations are made in a transparent geometry model to identify the regions of fluid mechanical interests, e.g. regions of flow separation, eddy formation and flow stagnancy. The 130 electrodes are strategically positioned in the arch based on information obtained from the flow visualizations. The measured data indicate that the aortic arch may be categorized into eight regions: three along the inner wall of the arch (A,B,C); and five near the outer wall (D,E,F,G,H). (1) The regions of low WSS are distributed along the inner wall of the ascending aorta A; the inner wall of the descending aorta C; and the upstream inner wall of the innominate and the common carotid branchings F. (2) The high WSS regions are distributed along the outer wall of the arch E; and the inner wall in the arch opposite to the left subclavian branching B. (3) In certain regions, high and low WSS may be found next to each other (e.g. G and H) without a definable boundary in between; and (4) as the Reynolds number increases, the areas of low WSS decrease, while the high WSS areas increase with no obvious change in magnitude of the stress along the inner wall of the arch. At the branchings, the WSS distribution is not affected by the Reynolds number within the range of observations. The measured WSS distribution is compared with Rodkiewicz's map of early atherosclerotic lesions in the aortic arch of cholesterol fed rabbits.