Coronary artery wall shear stress is associated with progression and transformation of atherosclerotic plaque and arterial remodeling in patients with coronary artery disease

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.

Computation in the rabbit aorta of a new metric - the transverse wall shear stress - to quantify the multidirectional character of disturbed blood flow

Spatial variation of the haemodynamic stresses acting on the arterial wall is commonly assumed to explain the focal development of atherosclerosis. Disturbed flow in particular is thought to play a key role. However, widely-used metrics developed to quantify its extent are unable to distinguish between uniaxial and multidirectional flows. We analysed pulsatile flow fields obtained in idealised and anatomically-realistic arterial geometries using computational fluid dynamics techniques, and in particular investigated the multidirectionality of the flow fields, capturing this aspect of near-wall blood flow with a new metric - the transverse wall shear stress (transWSS) - calculated as the time-average of wall shear stress components perpendicular to the mean flow direction. In the idealised branching geometry, multidirectional flow was observed downstream of the branch ostium, a region of flow stagnation, and to the sides of the ostium. The distribution of the transWSS was different from the pattern of traditional haemodynamic metrics and more dependent on the velocity waveform imposed at the branch outlet. In rabbit aortas, transWSS patterns were again different from patterns of traditional metrics. The near-branch pattern varied between intercostal ostia, as is the case for lesion distribution; for some branches there were striking resemblances to the age-dependent patterns of disease seen in rabbit and human aortas. The new metric may lead to improved understanding of atherogenesis.

Correlations between carotid plaque progression and mechanical stresses change sign over time: a patient follow up study using MRI and 3D FSI models

BACKGROUND

Increasing evidence suggests that mechanisms governing advanced plaque progression may be different from those for early progression and require further investigation. Serial MRI data and 3D fluid-structure interaction (FSI) models were employed to identify possible correlations between mechanical stresses and advanced plaque progression measured by vessel wall thickness increase (WTI). Long-term patient follow up was used to gather data and investigate if the correlations identified above were reproducible.

METHODS

In vivo MRI data were acquired from 16 patients in a follow-up study with 2 to 4 scans for each patient (scan interval: average 18 months and standard deviation 6.8 months). A total of 38 scan pairs (baseline and follow-up) were formed for analysis using the carotid bifurcation as the registration point. 3D FSI models were constructed to obtain plaque wall stress (PWS) and flow shear stress (FSS) to quantify their correlations with plaque progression. The Linear Mixed-Effects models were used to study possible correlations between WTI and baseline PWS and FSS with nodal dependence taken into consideration.

RESULTS

Of the 38 scan pairs, 22 pairs showed positive correlation between baseline PWS and WTI, 1 pair showed negative correlation, and 15 pairs showed no correlation. Thirteen patients changed their correlation sign (81.25%). Between baseline FSS and WTI, 16 pairs showed negative correlation, 1 pair showed positive correlation. Twelve patients changed correlation sign (75%).

CONCLUSION

Our results showed that advanced plaque progression had an overall positive correlation with plaque wall stress and a negative correlation with flow shear stress at baseline. However, long-term follow up showed that correlations between plaque progress and mechanical stresses (FSS and PWS) identified for one time period were not re-producible for most cases (>80%). Further investigations are needed to identify the reasons causing the correlation sign changes.

Platelet biology: the role of shear

Platelets are anucleated fragments produced by megakaryocytes that circulate in the blood. Platelets are involved in the initial cellular response to damaged endothelium and migrate to this area to prevent excessive bleeding. What is becoming more acknowledged over the last few decades is that blood flow (hemodynamics) plays a critical role in platelet function. The purpose of this review is to summarize the current understanding of platelet biology with particular focus on the role of hemodynamics. The emerging concept of shear microgradients, which are challenging the traditional model of platelet function, will also be introduced in the review.

Computational Fluid Dynamics Simulations of Hemodynamics in Plaque Erosion

PURPOSE

We investigated whether local hemodynamics were associated with sites of plaque erosion and hypothesized that patients with plaque erosion have locally elevated WSS magnitude in regions where erosion has occurred.

METHODS

We generated 3D, patient-specific models of coronary arteries from biplane angiographic images in 3 human patients with plaque erosion diagnosed by optical coherence tomography (OCT). Using computational fluid dynamics, we simulated pulsatile blood flow and calculated both wall shear stress (WSS) and oscillatory shear index (OSI). We also investigated anatomic features of plaque erosion sites by examining branching and local curvature in x-ray angiograms of barium-perfused autopsy hearts.

RESULTS

Neither high nor low magnitudes of mean WSS were associated with sites of plaque erosion. OSI and local curvature were also not associated with erosion. Anatomically, 8 of 13 hearts had a nearby bifurcation upstream of the site of plaque erosion.

CONCLUSIONS

This study provides preliminary evidence that neither hemodynamics nor anatomy are predictors of plaque erosion, based upon a very unique dataset. Our sample sizes are small, but this dataset suggests that high magnitudes of wall shear stress, one potential mechanism for inducing plaque erosion, are not necessary for erosion to occur.

Comparison of lipid deposition at coronary bifurcations versus at nonbifurcation portions of coronary arteries as determined by near-infrared spectroscopy

Atherosclerosis has been shown to develop preferentially at sites of coronary bifurcation, yet culprit lesions resulting in ST-elevation myocardial infarction do not occur more frequently at these sites. We hypothesized that these findings can be explained by similarities in intracoronary lipid and that lipid and lipid core plaque would be found with similar frequency in coronary bifurcation and nonbifurcation segments. One hundred seventy bifurcations were identified, 156 of which had comparative nonbifurcation segments proximal and/or distal to the bifurcation. We compared lipid deposition at bifurcation and nonbifurcation segments in coronary arteries using near-infrared spectroscopy (NIRS), a novel method for the in vivo detection of coronary lipid. Any NIRS signal for the presence of lipid was found with similar frequency in bifurcation and nonbifurcation segments (79% vs 74%, p = NS). Lipid core burden index, a measure of total lipid quantity indexed to segment length, was similar across bifurcation segments as well as their proximal and distal controls (lipid core burden index 66.3 ± 106, 67.1 ± 116, and 66.6 ± 104, p = NS). Lipid core plaque, identified as a high-intensity focal NIRS signal, was found in 21% of bifurcation segments, and 20% of distal nonbifurcation segments (p = NS). In conclusion, coronary bifurcations do not appear to have higher levels of intracoronary lipid or lipid core plaque than their comparative nonbifurcation regions.

Does low and oscillatory wall shear stress correlate spatially with early atherosclerosis? A systematic review

Low and oscillatory wall shear stress is widely assumed to play a key role in the initiation and development of atherosclerosis. Indeed, some studies have relied on the low shear theory when developing diagnostic and treatment strategies for cardiovascular disease. We wished to ascertain if this consensus is justified by published data. We performed a systematic review of papers that compare the localization of atherosclerotic lesions with the distribution of haemodynamic indicators calculated using computational fluid dynamics. The review showed that although many articles claim their results conform to the theory, it has been interpreted in different ways: a range of metrics has been used to characterize the distribution of disease, and they have been compared with a range of haemodynamic factors. Several studies, including all of those making systematic point-by-point comparisons of shear and disease, failed to find the expected relation. The various pre- and post-processing techniques used by different groups have reduced the range of shears over which correlations were sought, and in some cases are mutually incompatible. Finally, only a subset of the known patterns of disease has been investigated. The evidence for the low/oscillatory shear theory is less robust than commonly assumed. Longitudinal studies starting from the healthy state, or the collection of average flow metrics derived from large numbers of healthy vessels, both in conjunction with point-by-point comparisons using appropriate statistical techniques, will be necessary to improve our understanding of the relation between blood flow and atherogenesis.

Thin-capped atheromata with reduced collagen content in pigs develop in coronary arterial regions exposed to persistently low endothelial shear stress

OBJECTIVE

The mechanisms promoting the focal formation of rupture-prone coronary plaques in vivo remain incompletely understood. This study tested the hypothesis that coronary regions exposed to low endothelial shear stress (ESS) favor subsequent development of collagen-poor, thin-capped plaques.

APPROACH AND RESULTS

Coronary angiography and 3-vessel intravascular ultrasound were serially performed at 5 consecutive time points in vivo in 5 diabetic, hypercholesterolemic pigs. ESS was calculated along the course of each artery with computational fluid dynamics at all 5 time points. At follow-up, 184 arterial segments with previously identified in vivo ESS underwent histopathologic analysis. Compared with other plaque types, eccentric thin-capped atheromata developed more in segments that experienced lower ESS during their evolution. Compared with lesions with higher preceding ESS, segments persistently exposed to low ESS (<1.2 Pa) exhibited reduced intimal smooth muscle cell content; marked intimal smooth muscle cell phenotypic modulation; attenuated procollagen-I gene expression; increased gene and protein expression of the interstitial collagenases matrix-metalloproteinase-1, -8, -13, and -14; increased collagenolytic activity; reduced collagen content; and marked thinning of the fibrous cap.

CONCLUSIONS

Eccentric thin-capped atheromata, lesions particularly prone to rupture, form more frequently in coronary regions exposed to low ESS throughout their evolution. By promoting an imbalance of attenuated synthesis and augmented collagen breakdown, low ESS favors the focal evolution of early lesions toward plaques with reduced collagen content and thin fibrous caps-2 critical determinants of coronary plaque vulnerability.

Early detection and invasive passivation of future culprit lesions: a future potential or an unrealistic pursuit of chimeras?

New advances in image and signal processing have allowed the development of numerous invasive and noninvasive imaging modalities that have revealed details of plaque pathology and allowed us to study in vivo the atherosclerotic evolution. Recent natural history of atherosclerosis studies permitted us to evaluate changes in the compositional and morphological characteristics of the plaque and identify predictors of future events. The idea of being able to identify future culprit lesions and passivate these plaques has gradually matured, and small scale studies have provided proofs about the feasibility of this concept. This review article summarizes the recent advances in the study of atherosclerosis, cites the current evidence, highlights our limitations in understanding the evolution of the plaque and in predicting plaque destabilization, and discusses the potentiality of an early invasive sealing of future culprit lesions.

Fractional flow reserve for the assessment of complex multivessel disease in a patient after hybrid coronary revascularization

We present a case of a 43-year-old woman with history of hybrid coronary revascularization [endoscopic atraumatic coronary artery bypass (ACAB)] of left internal mammary artery (LIMA) to the left anterior descending artery (LAD) and stent implantation in right coronary artery (RCA), who presented 6 years later with recurrent atypical angina. Coronary angiography revealed patent LIMA to LAD and RCA stent, with a new lesion in an obtuse marginal artery and significant progression of disease in the proximal/mid LAD proximal to LIMA touchdown. To further evaluate the hemodynamic significance of these new disease segments, the patient underwent fractional flow reserve (FFR) assessment of the left coronary system with subsequent stent implantation in the proximal/mid LAD. This case illustrates (1) the critical value of FFR assessment in determining the ischemia provoking lesions in this post ACAB patient with complex multivessel coronary artery disease; and (2) the accelerated progression of atherosclerosis in bypassed segments as compared to segments proximal to stents.

Wall shear stress variations and unsteadiness of pulsatile blood-like flows in 90-degree bifurcations

Complex and slow interaction of different mechanical and biochemical processes in hemodynamics is believed to govern atherogenesis. Over the last decades studies have shown that fluid mechanical factors such as the Wall Shear Stress (WSS) and WSS gradients can play an important role in the pathological changes of the endothelium. This study provides further indications that the effects of fluid mechanical aspects are correlated with the diseased regions of the larger arteries. Unsteady high temporal WSS gradients (TWSSG), a function of the shear-thinning property of the non-Newtonian viscosity, move with the separation bubble. Red Blood Cell (RBC) dilution due to the secondary flows determines the magnitudes of the WSS and TWSSG. The results indicate that the focal nature of the TWSSG may have implications on the response of the endothelium.

A finite element study on variations in mass transport in stented porcine coronary arteries based on location in the coronary arterial tree

Drug-eluting stents have a significant clinical advantage in late-stage restenosis due to the antiproliferative drug release. Understanding how drug transport occurs between coronary arterial locations can better help guide localized drug treatment options. Finite element models with properties from specific porcine coronary artery sections (left anterior descending (LAD), right (RCA); proximal, middle, distal regions) were created for stent deployment and drug delivery simulations. Stress, strain, pore fluid velocity, and drug concentrations were exported at different time points of simulation (0-180 days). Tests indicated that the highest stresses occurred in LAD sections. Higher-than-resting homeostatic levels of stress and strain existed at upwards of 3.0 mm away from the stented region, whereas concentration of species only reached 2.7 mm away from the stented region. Region-specific concentration showed 2.2 times higher concentrations in RCA artery sections at times corresponding to vascular remodeling (peak in the middle segment) compared to all other segments. These results suggest that wall transport can occur differently based on coronary artery location. Awareness of peak growth stimulators and where drug accumulation occurs in the vasculature can better help guide local drug delivery therapies.

Exploring coronary atherosclerosis with intravascular imaging

Coronary angiography has been widely used for five decades to evaluate a range of vascular pathologies and triage patients to therapeutic interventions. The inability to directly visualize the artery wall with conventional angiographic techniques has stimulated development of a number of intravascular imaging modalities. These approaches have the potential to provide a more comprehensive characterization of the burden, composition and functionality of atherosclerotic plaque, neointimal hyperplasia and allograft vasculopathy that develop within coronary arteries. The ability to use these modalities in vivo and in a serial fashion has provided a greater insight into the factors that underlie the disease process and guide therapeutic interventions.

Compound ex vivo and in silico method for hemodynamic analysis of stented arteries

Hemodynamic factors such as low wall shear stress have been shown to influence endothelial healing and atherogenesis in stent-free vessels. However, in stented vessels, a reliable quantitative analysis of such relations has not been possible due to the lack of a suitable method for the accurate acquisition of blood flow. The objective of this work was to develop a method for the precise reconstruction of hemodynamics and quantification of wall shear stress in stented vessels. We have developed such a method that can be applied to vessels stented in or ex vivo and processed ex vivo. Here we stented the coronary arteries of ex vivo porcine hearts, performed vascular corrosion casting, acquired the vessel geometry using micro-computed tomography and reconstructed blood flow and shear stress using computational fluid dynamics. The method yields accurate local flow information through anatomic fidelity, capturing in detail the stent geometry, arterial tissue prolapse, radial and axial arterial deformation as well as strut malapposition. This novel compound method may serve as a unique tool for spatially resolved analysis of the relationship between hemodynamic factors and vascular biology. It can further be employed to optimize stent design and stenting strategies.

Hemodynamic impacts of various types of stenosis in the left coronary artery bifurcation: a patient-specific analysis

This study investigates the hemodynamic changes to various types of coronary stenosis in the left coronary artery bifurcation, based on a patient-specific analysis. Twenty two patients with left coronary artery disease were included in this study. All stenoses involving the left coronary artery bifurcation were classified into four types, according to their locations: A) left circumflex (LCx) and left anterior descending (LAD), B) LCx only, C) left main stem only, and D) LAD only. Computational fluid dynamics (CFD) was performed to analyze the flow and wall shear stress (WSS) changes in all reconstructed left coronary geometries. Our results showed that the flow velocity and WSS were significantly increased at stenotic locations. High WSS was found at >70% lumen stenosis, which ranged from 2.5 Pa to 3.5 Pa. This study demonstrates that in patients with more than 50% stenosis in the left coronary artery bifurcation, WSS plays an important role in providing information about the extent of coronary atherosclerosis in the left coronary artery branch.

Flow recirculation zone length and shear rate are differentially affected by stenosis severity in human coronary arteries

Flow recirculation zones and shear rate are associated with distinct pathogenic biological pathways relevant to thrombosis and atherogenesis. The interaction between stenosis severity and lesion eccentricity in determining the length of flow recirculation zones and peak shear rate in human coronary arteries in vivo is unclear. Computational fluid dynamic simulations were performed under resting and hyperemic conditions on computer-generated models and three-dimensional (3-D) reconstructions of coronary arteriograms of 25 patients. Boundary conditions for 3-D reconstructions simulations were obtained by direct measurements using a pressure-temperature sensor guidewire. In the computer-generated models, stenosis severity and lesion eccentricity were strongly associated with recirculation zone length and maximum shear rate. In the 3-D reconstructions, eccentricity increased recirculation zone length and shear rate when lesions of the same stenosis severity were compared. However, across the whole population of coronary lesions, eccentricity did not correlate with recirculation zone length or shear rate ( P = not signficant for both), whereas stenosis severity correlated strongly with both parameters ( r = 0.97, P < 0.001, and r = 0.96, P < 0.001, respectively). Nonlinear regression analyses demonstrated that the relationship between stenosis severity and peak shear was exponential, whereas the relationship between stenosis severity and recirculation zone length was sigmoidal, with an apparent threshold effect, demonstrating a steep increase in recirculation zone length between 40% and 60% diameter stenosis. Increasing stenosis severity and lesion eccentricity can both increase flow recirculation and shear rate in human coronary arteries. Flow recirculation is much more sensitive to mild changes in the severity of intermediate stenoses than is peak shear.

A novel method for quantifying spatial correlations between patterns of atherosclerosis and hemodynamic factors

Studies investigating the relation between the focal nature of atherosclerosis and hemodynamic factors are employing increasingly rigorous approaches to map the disease and calculate hemodynamic metrics. However, no standardized methodology exists to quantitatively compare these distributions. We developed a statistical technique that can be used to determine if hemodynamic and lesion maps are significantly correlated. The technique, which is based on a surrogate data analysis, does not require any assumptions (such as linearity) on the nature of the correlation. Randomized sampling was used to ensure the independence of data points, another basic assumption of commonly-used statistical methods that is often disregarded. The novel technique was used to compare previously-obtained maps of lesion prevalence in aortas of immature and mature cholesterol-fed rabbits to corresponding maps of wall shear stress, averaged across several animals in each age group. A significant spatial correlation was found in the proximal descending thoracic aorta, but not further downstream. Around intercostal branch openings the correlation was borderline significant in immature but not in mature animals. The results confirm the need for further investigation of the relation between the localization of atherosclerosis and blood flow, in conjunction with appropriate statistical techniques such as the method proposed here.

Image-based modeling and precision medicine: patient-specific carotid and coronary plaque assessment and predictions

Atherosclerotic plaques may rupture without warning and cause acute cardiovascular events such as heart attack and stroke. Current clinical screening tools are insufficient to identify those patients with risks early and prevent the adverse events from happening. Medical imaging and image-based modeling have made considerable progress in recent years in identifying plaque morphological and mechanical risk factors which may be used in developing improved patient screening strategies. The key steps and factors in image-based models for human carotid and coronary plaques were illustrated, as well as grand challenges facing the researchers in the field to develop more accurate screening tools.

Model generation of coronary artery bifurcations from CTA and single plane angiography

PURPOSE

To generate accurate and realistic models of coronary artery bifurcations before and after percutaneous coronary intervention (PCI), using information from two image modalities. Because bifurcations are regions where atherosclerotic plaque appears frequently and intervention is more challenging, generation of such realistic models could be of high value to predict the risk of restenosis or thrombosis after stent implantation, and to study geometrical and hemodynamical changes.

METHODS

Two image modalities have been employed to generate the bifurcation models: computer tomography angiography (CTA) to obtain the 3D trajectory of vessels, and 2D conventional coronary angiography (CCA) to obtain radius information of the vessel lumen, due to its better contrast and image resolution. In addition, CCA can be acquired right before and after the intervention in the operation room; therefore, the combination of CTA and CCA allows the generation of realistic preprocedure and postprocedure models of coronary bifurcations. The method proposed is semiautomatic, based on landmarks manually placed on both image modalities.

RESULTS

A comparative study of the models obtained with the proposed method with models manually obtained using only CTA, shows more reliable results when both modalities are used together. The authors show that using preprocedure CTA and postprocedure CCA, realistic postprocedure models can be obtained. Analysis carried out of the Murray's law in all patient bifurcations shows the geometric improvement of PCI in our models, better than using manual models from CTA alone. An experiment using a cardiac phantom also shows the feasibility of the proposed method.

CONCLUSIONS

The authors have shown that fusion of CTA and CCA is feasible for realistic generation of coronary bifurcation models before and after PCI. The method proposed is efficient, and relies on minimal user interaction, and therefore is of high value to study geometric and hemodynamic changes of treated patients.

Elevated electrochemical impedance in the endoluminal regions with high shear stress: implication for assessing lipid-rich atherosclerotic lesions

BACKGROUND

Identifying metabolically active atherosclerotic lesions remains an unmet clinical challenge during coronary intervention. Electrochemical impedance (EIS) increased in response to oxidized low density lipoprotein (oxLDL)-laden lesions. We hereby assessed whether integrating EIS with intravascular ultrasound (IVUS) and shear stress (ISS) provided a new strategy to assess oxLDL-laden lesions in the fat-fed New Zealand White (NZW) rabbits.

METHODS AND RESULTS

A micro-heat transfer sensor was deployed to acquire the ISS profiles at baseline and post high-fat diet (HD) in the NZW rabbits (n=8). After 9 weeks of HD, serum oxLDL levels (mg/dL) increased by 140 fold, accompanied by a 1.5-fold increase in kinematic viscosity (cP) in the HD group. Time-averaged ISS (ISSave) in the thoracic aorta also increased in the HD group (baseline: 17.61±0.24 vs. 9 weeks: 25.22±0.95dyne/cm(2), n=4), but remained unchanged in the normal diet group (baseline: 22.85±0.53dyn/cm(2) vs. 9 weeks: 22.37±0.57dyne/cm(2), n=4). High-frequency intravascular ultrasound (IVUS) revealed atherosclerotic lesions in the regions with augmented ISSave, and concentric bipolar microelectrodes demonstrated elevated EIS signals, which were correlated with prominent anti-oxLDL immuno-staining (oxLDL-free regions: 497±55Ω, n=8 vs. oxLDL-rich lesions: 679±125Ω, n=12, P<0.05). The equivalent circuit model for tissue resistance between the lesion-free and ox-LDL-rich lesions further validated the experimental EIS signals.

CONCLUSIONS

By applying electrochemical impedance in conjunction with shear stress and high-frequency ultrasound sensors, we provided a new strategy to identify oxLDL-laden lesions. The study demonstrated the feasibility of integrating EIS, ISS, and IVUS for a catheter-based approach to assess mechanically unstable plaque.

High wall shear stress and spatial gradients in vascular pathology: a review

Cardiovascular pathologies such as intracranial aneurysms (IAs) and atherosclerosis preferentially localize to bifurcations and curvatures where hemodynamics are complex. While extensive knowledge about low wall shear stress (WSS) has been generated in the past, due to its strong relevance to atherogenesis, high WSS (typically >3 Pa) has emerged as a key regulator of vascular biology and pathology as well, receiving renewed interests. As reviewed here, chronic high WSS not only stimulates adaptive outward remodeling, but also contributes to saccular IA formation (at bifurcation apices or outer curves) and atherosclerotic plaque destabilization (in stenosed vessels). Recent advances in understanding IA pathogenesis have shed new light on the role of high WSS in pathological vascular remodeling. In complex geometries, high WSS can couple with significant spatial WSS gradient (WSSG). A combination of high WSS and positive WSSG has been shown to trigger aneurysm initiation. Since endothelial cells (ECs) are sensors of WSS, we have begun to elucidate EC responses to high WSS alone and in combination with WSSG. Understanding such responses will provide insight into not only aneurysm formation, but also plaque destabilization and other vascular pathologies and potentially lead to improved strategies for disease management and novel targets for pharmacological intervention.

Association of coronary wall shear stress with atherosclerotic plaque burden, composition, and distribution in patients with coronary artery disease

Background

Extremes of wall shear stress (WSS) have been associated with plaque progression and transformation, which has raised interest in the clinical assessment of WSS. We hypothesized that calculated coronary WSS is predicted only partially by luminal geometry and that WSS is related to plaque composition.

Methods and Results

Twenty‐seven patients with coronary artery disease underwent virtual histology intravascular ultrasound and Doppler velocity measurement for computational fluid dynamics modeling for WSS calculation in each virtual histology intravascular ultrasound segment (N=3581 segments). We assessed the association of WSS with plaque burden and distribution and with plaque composition. WSS remained relatively constant across the lower 3 quartiles of plaque burden (P=0.08) but increased in the highest quartile of plaque burden (P<0.001). Segments distal to lesions or within bifurcations were more likely to have low WSS (P<0.001). However, the majority of segments distal to lesions (80%) and within bifurcations (89%) did not exhibit low WSS. After adjustment for plaque burden, there was a negative association between WSS and percent necrotic core and calcium. For every 10 dynes/cm² increase in WSS, percent necrotic core decreased by 17% (P=0.01), and percent dense calcium decreased by 17% (P<0.001). There was no significant association between WSS and percent of fibrous or fibrofatty plaque components (P=NS).

Conclusions

In patients with coronary artery disease: (1) Luminal geometry predicts calculated WSS only partially, which suggests that detailed computational techniques must be used to calculate WSS. (2) Low WSS is associated with plaque necrotic core and calcium, independent of plaque burden, which suggests a link between WSS and coronary plaque phenotype. (J Am Heart Assoc. 2012;1:e002543 doi: 10.1161/JAHA.112.002543.)

Correlations of coronary plaque wall thickness with wall pressure and wall pressure gradient: a representative case study

BACKGROUND

There are two major hemodynamic stresses imposed at the blood arterial wall interface by flowing blood: the wall shear stress (WSS) acting tangentially to the wall, and the wall pressure (WP) acting normally to the wall. The role of flow wall shear stress in atherosclerosis progression has been under intensive investigation, while the impact of blood pressure on plaque progression has been under-studied.

METHOD

The correlations of wall thickness (WT) with wall pressure (WP, blood pressure on the lumen wall) and spatial wall pressure gradient (WPG) in a human atherosclerotic right coronary artery were studied. The pulsatile blood flow was simulated using a three dimensional mathematical model. The blood was treated as an incompressible viscous non-Newtonian fluid. The geometry of the artery was re-constructed using an in vivo intravascular ultrasound (IVUS) 44-slice dataset obtained from a patient with consent obtained. The WT, the WP and the WPG were averaged on each slice, respectively, and Pearson correlation analysis was performed on slice averaged base. Each slice was then divided into 8 segments and averaged vessel WT, WP and WPG were collected from all 352 segments for correlation analysis. Each slice was also divided into 2 segments (inner semi-wall of bend and outer semi-wall of bend) and the correlation analysis was performed on the 88 segments.

RESULTS

Under mean pressure, the Pearson coefficient for correlation between WT and WP was r = - 0.52 (p < 0.0001) by 2-segment analysis and r = - 0.81 (p < 0.0001) by slice averaged analysis, respectively. The Pearson coefficient for correlation between WT and WPG was r = 0.30 (p = 0.004) by 2-segment analysis and r = 0.45 (p = 0.002) by slice averaged analysis, respectively. The r-values corresponding to systole and diastole pressure conditions were similar.

CONCLUSIONS

Results from this representative case report indicated that plaque wall thickness correlated negatively with wall pressure (r = -0.81 by slice) and positively with wall pressure gradient (r = 0.45). The slice averaged WT has a strong linear relationship with the slice averaged WP. Large-scale patient studies are needed to further confirm our findings.

Prediction of progression of coronary artery disease and clinical outcomes using vascular profiling of endothelial shear stress and arterial plaque characteristics: the PREDICTION Study

BACKGROUND

Atherosclerotic plaques progress in a highly individual manner. The purposes of the Prediction of Progression of Coronary Artery Disease and Clinical Outcome Using Vascular Profiling of Shear Stress and Wall Morphology (PREDICTION) Study were to determine the role of local hemodynamic and vascular characteristics in coronary plaque progression and to relate plaque changes to clinical events.

METHODS AND RESULTS

Vascular profiling, using coronary angiography and intravascular ultrasound, was used to reconstruct each artery and calculate endothelial shear stress and plaque/remodeling characteristics in vivo. Three-vessel vascular profiling (2.7 arteries per patient) was performed at baseline in 506 patients with an acute coronary syndrome treated with a percutaneous coronary intervention and in a subset of 374 (74%) consecutive patients 6 to 10 months later to assess plaque natural history. Each reconstructed artery was divided into sequential 3-mm segments for serial analysis. One-year clinical follow-up was completed in 99.2%. Symptomatic clinical events were infrequent: only 1 (0.2%) cardiac death; 4 (0.8%) patients with new acute coronary syndrome in nonstented segments; and 15 (3.0%) patients hospitalized for stable angina. Increase in plaque area (primary end point) was predicted by baseline large plaque burden; decrease in lumen area (secondary end point) was independently predicted by baseline large plaque burden and low endothelial shear stress. Large plaque size and low endothelial shear stress independently predicted the exploratory end points of increased plaque burden and worsening of clinically relevant luminal obstructions treated with a percutaneous coronary intervention at follow-up. The combination of independent baseline predictors had a 41% positive and 92% negative predictive value to predict progression of an obstruction treated with a percutaneous coronary intervention.

CONCLUSIONS

Large plaque burden and low local endothelial shear stress provide independent and additive prediction to identify plaques that develop progressive enlargement and lumen narrowing.

CLINICAL TRIAL REGISTRATION

URL: http:www.//clinicaltrials.gov. Unique Identifier: NCT01316159.

Endothelial shear stress in the evolution of coronary atherosclerotic plaque and vascular remodelling: current understanding and remaining questions

The heterogeneity of plaque formation, the vascular remodelling response to plaque formation, and the consequent phenotype of plaque instability attest to the extraordinarily complex pathobiology of plaque development and progression, culminating in different clinical coronary syndromes. Atherosclerotic plaques predominantly form in regions of low endothelial shear stress (ESS), whereas regions of moderate/physiological and high ESS are generally protected. Low ESS-induced compensatory expansive remodelling plays an important role in preserving lumen dimensions during plaque progression, but when the expansive remodelling becomes excessive promotes continued influx of lipids into the vessel wall, vulnerable plaque formation and potential precipitation of an acute coronary syndrome. Advanced plaques which start to encroach into the lumen experience high ESS at their most stenotic region, which appears to promote plaque destabilization. This review describes the role of ESS from early atherogenesis to early plaque formation, plaque progression to advanced high-risk stenotic or non-stenotic plaque, and plaque destabilization. The critical implication of the vascular remodelling response to plaque growth is also discussed. Current developments in technology to characterize local ESS and vascular remodelling in vivo may provide a rationale for innovative diagnostic and therapeutic strategies for coronary patients that aim to prevent clinical coronary syndromes.

Role of endothelial shear stress in stent restenosis and thrombosis: pathophysiologic mechanisms and implications for clinical translation

Restenosis and thrombosis are potentially fatal complications of coronary stenting with a recognized multifactorial etiology. The effect of documented risk factors, however, cannot explain the preponderance of certain lesion types, stent designs, and implantation configurations for the development of these complications. Local hemodynamic factors, low endothelial shear stress (ESS) in particular, are long known to critically affect the natural history of atherosclerosis. Increasing evidence now suggests that ESS may also contribute to the development of restenosis and thrombosis upon stenting of atherosclerotic plaques, in conjunction with well-appreciated risk factors. In this review, we present in vivo and mechanistic evidence associating ESS with the localization and progression of neointimal hyperplasia and in-stent clotting. Clinical studies have associated stent design features with the risk of restenosis. Importantly, computational simulations extend these observations by directly linking specific stent geometry and positioning characteristics with the post-stenting hemodynamic milieu and with the stent's thrombogenicity and pro-restenotic potential, thereby indicating ways to clinical translation. An enhanced understanding of the pathophysiologic role of ESS in restenosis and thrombosis might dictate hemodynamically favorable stent designs and deployment configurations to reduce the potential for late lumen loss and thrombotic obstruction. Recent methodologies for in vivo ESS profiling at a clinical level might allow for early identification of patients at high risk for the development of restenosis or thrombosis and might thereby guide individualized, risk-tailored treatment strategies to prevent devastating complications of endovascular interventions.

Plaque characteristics and arterial remodeling in coronary and peripheral arterial systems

BACKGROUND

Few studies have examined plaque characteristics among multiple arterial beds in vivo. The purpose of this study was to compare the plaque morphology and arterial remodeling between coronary and peripheral arteries using gray-scale and radiofrequency intravascular ultrasound (IVUS) at clinical presentation.

METHODS AND RESULTS

IVUS imaging was performed in 68 patients with coronary and 93 with peripheral artery lesions (29 carotid, 50 renal, and 14 iliac arteries). Plaques were classified as fibroatheroma (VH-FA) (further subclassified as thin-capped [VH-TCFA] and thick-capped [VH-ThCFA]), fibrocalcific plaque (VH-FC) and pathological intimal thickening (VH-PIT). Plaque rupture (13% of coronary, 7% of carotid, 6% of renal, and 7% of iliac arteries; P = NS) and VH-TCFA (37% of coronary, 24% of carotid, 16% of renal, and 7% of iliac arteries; P = 0.02) were observed in all arteries. Compared with coronary arteries, VH-FA was less frequently observed in renal (P < 0.001) and iliac arteries (P < 0.006). Lesions with positive remodeling demonstrated more characteristics of VH-FA in coronary (84% vs. 25%, P < 0.001), carotid (72% vs. 20%, P = 0.001), and renal arteries (42% vs. 4%, P = 0.001) compared with those with intermediate/negative remodeling. There was positive relationship between remodeling index and percent necrotic area in all four arteries.

CONCLUSIONS

Atherosclerotic plaque phenotypes were heterogeneous among four different arteries; renal and iliac arteries had more stable phenotypes compared with coronary artery. In contrast, the associations of remodeling pattern with plaque phenotype and composition were similar among the various arterial beds.

Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery

We present a computational framework for multiscale modeling and simulation of blood flow in coronary artery bypass graft (CABG) patients. Using this framework, only CT and non-invasive clinical measurements are required without the need to assume pressure and/or flow waveforms in the coronaries and we can capture global circulatory dynamics. We demonstrate this methodology in a case study of a patient with multiple CABGs. A patient-specific model of the blood vessels is constructed from CT image data to include the aorta, aortic branch vessels (brachiocephalic artery and carotids), the coronary arteries and multiple bypass grafts. The rest of the circulatory system is modeled using a lumped parameter network (LPN) 0 dimensional (0D) system comprised of resistances, capacitors (compliance), inductors (inertance), elastance and diodes (valves) that are tuned to match patient-specific clinical data. A finite element solver is used to compute blood flow and pressure in the 3D (3 dimensional) model, and this solver is implicitly coupled to the 0D LPN code at all inlets and outlets. By systematically parameterizing the graft geometry, we evaluate the influence of graft shape on the local hemodynamics, and global circulatory dynamics. Virtual manipulation of graft geometry is automated using Bezier splines and control points along the pathlines. Using this framework, we quantify wall shear stress, wall shear stress gradients and oscillatory shear index for different surgical geometries. We also compare pressures, flow rates and ventricular pressure-volume loops pre- and post-bypass graft surgery. We observe that PV loops do not change significantly after CABG but that both coronary perfusion and local hemodynamic parameters near the anastomosis region change substantially. Implications for future patient-specific optimization of CABG are discussed.

Impact of plaques in the left coronary artery on wall shear stress and pressure gradient in coronary side branches

In this study, we investigate plaques located at the left coronary bifurcation. We focus on the effect that the resulting changes in wall shear stress (WSS) and wall pressure stress gradient (WPSG) have on atherosclerotic progress in coronary artery disease. Coronary plaques were simulated and placed at the left main stem and the left anterior descending to produce >50% narrowing of the coronary lumen. Computational fluid dynamics analysis was carried out, simulating realistic physiological conditions that show the in vivo cardiac haemodynamic. WSS and WPSG in the left coronary artery were calculated and compared in the left coronary models, with and without the presence of plaques during cardiac cycles. Our results showed that WSS decreased while WPSG was increased in coronary side branches due to the presence of plaques. There is a direct correlation between coronary plaques and subsequent WSS and WPSG variations based on the bifurcation plaques simulated in the realistic coronary models.

Microcalcifications increase coronary vulnerable plaque rupture potential: a patient-based micro-CT fluid-structure interaction study

Asymptomatic vulnerable plaques (VP) in coronary arteries accounts for significant level of morbidity. Their main risk is associated with their rupture which may prompt fatal heart attacks and strokes. The role of microcalcifications (micro-Ca), embedded in the VP fibrous cap, in the plaque rupture mechanics has been recently established. However, their diminutive size offers a major challenge for studying the VP rupture biomechanics on a patient specific basis. In this study, a highly detailed model was reconstructed from a post-mortem coronary specimen of a patient with observed VP, using high resolution micro-CT which captured the microcalcifications embedded in the fibrous cap. Fluid-structure interaction (FSI) simulations were conducted in the reconstructed model to examine the combined effects of micro-Ca, flow phase lag and plaque material properties on plaque burden and vulnerability. This dynamic fibrous cap stress mapping elucidates the contribution of micro-Ca and flow phase lag VP vulnerability independently. Micro-Ca embedded in the fibrous cap produced increased stresses predicted by previously published analytical model, and corroborated our previous studies. The 'micro-CT to FSI' methodology may offer better diagnostic tools for clinicians, while reducing morbidity and mortality rates for patients with vulnerable plaques and ameliorating the ensuing healthcare costs.