Tissue-engineered collagenous fibrous cap models to systematically elucidate atherosclerotic plaque rupture.

A significant amount of vascular thrombotic events are associated with rupture of the fibrous cap that overlie atherosclerotic plaques. Cap rupture is however difficult to predict due to the heterogenous composition of the plaque, unknown material properties, and the stochastic nature of the event. Here, we aim to create tissue engineered human fibrous cap models with a variable but controllable collagen composition, suitable for mechanical testing, to scrutinize the reciprocal relationships between composition and mechanical properties. Myofibroblasts were cultured in 1 × 1.5 cm-sized fibrin-based constrained gels for 21 days according to established (dynamic) culture protocols (i.e. static, intermittent or continuous loading) to vary collagen composition (e.g. amount, type and organization). At day 7, a soft 2 mm ∅ fibrin inclusion was introduced in the centre of each tissue to mimic the soft lipid core, simulating the heterogeneity of a plaque. Results demonstrate reproducible collagenous tissues, that mimic the bulk mechanical properties of human caps and vary in collagen composition due to the presence of a successfully integrated soft inclusion and the culture protocol applied. The models can be deployed to assess tissue mechanics, evolution and failure of fibrous caps or complex heterogeneous tissues in general.

SPECT/CT imaging of inflammation and calcification in human carotid atherosclerosis to identify the plaque at risk of rupture.

BACKGROUND: Calcification and inflammation are atherosclerotic plaque compositional biomarkers that have both been linked to stroke risk. The aim of this study was to evaluate their co-existing prevalence in human carotid plaques with respect to plaque phenotype to determine the value of hybrid imaging for the detection of these biomarkers. METHODS: Human carotid plaque segments, obtained from endarterectomy, were incubated in [111In]In-DOTA-butylamino-NorBIRT ([111In]In-Danbirt), targeting Leukocyte Function-associated Antigen-1 (LFA-1) on leukocytes. By performing SPECT/CT, both inflammation from DANBIRT uptake and calcification from CT imaging were assessed. Plaque phenotype was classified using histology. RESULTS: On a total plaque level, comparable levels of calcification volume existed with different degrees of inflammation and vice versa. On a segment level, an inverse relationship between calcification volume and inflammation was evident in highly calcified segments, which classify as fibrocalcific, stable plaque segments. In contrast, segments with little or no calcification presented with a moderate to high degree of inflammation, often coinciding with the more dangerous fibrous cap atheroma phenotype. CONCLUSION: Calcification imaging alone can only accurately identify highly calcified, stable, fibrocalcific plaques. To identify high-risk plaques, with little or no calcification, hybrid imaging of calcification and inflammation could provide diagnostic benefit.

The Correlation Between Wall Shear Stress and Plaque Composition in Advanced Human Carotid Atherosclerosis.

The role of wall shear stress (WSS) in atherosclerotic plaque development is evident, but the relation between WSS and plaque composition in advanced atherosclerosis, potentially resulting in plaque destabilization, is a topic of discussion. Using our previously developed image registration pipeline, we investigated the relation between two WSS metrics, time-averaged WSS (TAWSS) and the oscillatory shear index (OSI), and the local histologically determined plaque composition in a set of advanced human carotid plaques. Our dataset of 11 carotid endarterectomy samples yielded 87 histological cross-sections, which yielded 511 radial bins for analysis. Both TAWSS and OSI values were subdivided into patient-specific low, mid, and high tertiles. This cross-sectional study shows that necrotic core (NC) size and macrophage area are significantly larger in areas exposed to high TAWSS or low OSI. Local TAWSS and OSI tertile values were generally inversely related, as described in the literature, but other combinations were also found. Investigating the relation between plaque vulnerability features and different combinations of TAWSS and OSI tertile values revealed a significantly larger cap thickness in areas exposed to both low TAWSS and low OSI. In conclusion, our study confirmed previous findings, correlating high TAWSS to larger macrophage areas and necrotic core sizes. In addition, our study demonstrated new relations, correlating low OSI to larger macrophage areas, and a combination of low TAWSS and low OSI to larger cap thickness.

An MRI-based method to register patient-specific wall shear stress data to histology.

Wall shear stress (WSS), the frictional force exerted on endothelial cells by blood flow, is hypothesised to influence atherosclerotic plaque growth and composition. We developed a methodology for image registration of MR and histology images of advanced human carotid plaques and corresponding WSS data, obtained by MRI and computational fluid dynamics. The image registration method requires four types of input images, in vivo MRI, ex vivo MRI, photographs of transversally sectioned plaque tissue and histology images. These images are transformed to a shared 3D image domain by applying a combination of rigid and non-rigid registration algorithms. Transformation matrices obtained from registration of these images are used to transform subject-specific WSS data to the shared 3D image domain as well. WSS values originating from the 3D WSS map are visualised in 2D on the corresponding lumen locations in the histological sections and divided into eight radial segments. In each radial segment, the correlation between WSS values and plaque composition based on histological parameters can be assessed. The registration method was successfully applied to two carotid endarterectomy specimens. The resulting matched contours from the imaging modalities had Hausdorff distances between 0.57 and 0.70 mm, which is in the order of magnitude of the in vivo MRI resolution. We simulated the effect of a mismatch in the rigid registration of imaging modalities on WSS results by relocating the WSS data with respect to the stack of histology images. A 0.6 mm relocation altered the mean WSS values projected on radial bins on average by 0.59 Pa, compared to the output of original registration. This mismatch of one image slice did not change the correlation between WSS and plaque thickness. In conclusion, we created a method to investigate correlations between WSS and plaque composition.

Flow Patterns in Carotid Webs: A Patient-Based Computational Fluid Dynamics Study.

BACKGROUND AND PURPOSE: Carotid webs are increasingly recognized as an important cause of (recurrent) ischemic stroke in patients without other cardiovascular risk factors. Hemodynamic flow patterns induced by these lesions might be associated with thrombus formation. The aim of our study was to evaluate flow patterns of carotid webs using computational fluid dynamics. MATERIALS AND METHODS: Patients with a carotid web in the Multicenter Randomized Clinical Trial of Endovascular Treatment of Acute Ischemic Stroke in the Netherlands (MR CLEAN) were selected for hemodynamic evaluation with computational fluid dynamics models based on lumen segmentations obtained from CT angiography scans. Hemodynamic parameters, including the area of recirculation zone, time-averaged wall shear stress, transverse wall shear stress, and the oscillatory shear index, were assessed and compared with the contralateral carotid bifurcation. RESULTS: In our study, 9 patients were evaluated. Distal to the carotid webs, recirculation zones were significantly larger compared with the contralateral bifurcation (63 versus 43 mm2, P = .02). In the recirculation zones of the carotid webs and the contralateral carotid bifurcation, time-averaged wall shear stress values were comparable (both: median, 0.27 Pa; P = .30), while transverse wall shear stress and oscillatory shear index values were significantly higher in the recirculation zone of carotid webs (median, 0.25 versus 0.21 Pa; P = .02 and 0.39 versus 0.30 Pa; P = .04). At the minimal lumen area, simulations showed a significantly higher time-averaged wall shear stress in the web compared with the contralateral bifurcation (median, 0.58 versus 0.45 Pa; P = .01). CONCLUSIONS: Carotid webs are associated with increased recirculation zones and regional increased wall shear stress metrics that are associated with disturbed flow. These findings suggest that a carotid web might stimulate thrombus formation, which increases the risk of acute ischemic stroke.

The effect of the heart rate lowering drug Ivabradine on hemodynamics in atherosclerotic mice.

The heart rate lowering drug Ivabradine was shown to improve cardiac outcome in patients with previous heart failure. However, in patients without heart failure, no beneficial effect of Ivabradine was observed. Animal studies suggested a preventive effect of Ivabradine on atherosclerosis which was due to an increase in wall shear stress (WSS), the blood flow-induced frictional force exerted on the endothelium, triggering anti-inflammatory responses. However, data on the effect of Ivabradine on WSS is sparse. We aim to study the effect of Ivabradine on (i) the 3D WSS distribution over a growing plaque and (ii) plaque composition. We induced atherosclerosis in ApoE-/- mice by placing a tapered cast around the right common carotid artery (RCCA). Five weeks after cast placement, Ivabradine was administered via drinking water (15 mg/kg/day) for 2 weeks, after which the RCCA was excised for histology analyses. Before and after Ivabradine treatment, animals were imaged with Doppler Ultrasound to measure blood velocity. Vessel geometry was obtained using contrast-enhanced micro-CT. Time-averaged WSS during systole, diastole and peak WSS was subsequently computed. Ivabradine significantly decreased heart rate (459 ± 28 bpm vs. 567 ± 32 bpm, p < 0.001). Normalized peak flow significantly increased in the Ivabradine group (124.2% ± 40.5% vs. 87.3% ± 25.4%, p < 0.05), reflected by an increased normalized WSS level during systole (110.7% ± 18.4% vs. 75.4% ± 24.6%, p < 0.05). However, plaque size or composition including plaque area, relative necrotic core area and macrophage content were not altered in mice treated with Ivabradine compared to controls. We conclude that increased WSS in response to Ivabradine treatment did not affect plaque progression in a murine model.

Assessment of human left ventricle flow using statistical shape modelling and computational fluid dynamics.

Blood flow patterns in the human left ventricle (LV) have shown relation to cardiac health. However, most studies in the literature are limited to a few patients and results are hard to generalize. This study aims to provide a new framework to generate more generalized insights into LV blood flow as a function of changes in anatomy and wall motion. In this framework, we studied the four-dimensional blood flow in LV via computational fluid dynamics (CFD) in conjunction with a statistical shape model (SSM), built from segmented LV shapes of 150 subjects. We validated results in an in-vitro dynamic phantom via time-resolved optical particle image velocimetry (PIV) measurements. This combination of CFD and the SSM may be useful for systematically assessing blood flow patterns in the LV as a function of varying anatomy and has the potential to provide valuable data for diagnosis of LV functionality.

Temporal and spatial changes in wall shear stress during atherosclerotic plaque progression in mice.

Wall shear stress (WSS) is involved in atherosclerotic plaque initiation, yet its role in plaque progression remains unclear. We aimed to study (i) the temporal and spatial changes in WSS over a growing plaque and (ii) the correlation between WSS and plaque composition, using animal-specific data in an atherosclerotic mouse model. Tapered casts were placed around the right common carotid arteries (RCCA) of ApoE-/- mice. At 5, 7 and 9 weeks after cast placement, RCCA geometry was reconstructed using contrast-enhanced micro-CT. Lumen narrowing was observed in all mice, indicating the progression of a lumen intruding plaque. Next, we determined the flow rate in the RCCA of each mouse using Doppler Ultrasound and computed WSS at all time points. Over time, as the plaque developed and further intruded into the lumen, absolute WSS significantly decreased. Finally at week 9, plaque composition was histologically characterized. The proximal part of the plaque was small and eccentric, exposed to relatively lower WSS. Close to the cast a larger and concentric plaque was present, exposed to relatively higher WSS. Lower WSS was significantly correlated to the accumulation of macrophages in the eccentric plaque. When pooling data of all animals, correlation between WSS and plaque composition was weak and no longer statistically significant. In conclusion, our data showed that in our mouse model absolute WSS strikingly decreased during disease progression, which was significantly correlated to plaque area and macrophage content. Besides, our study demonstrates the necessity to analyse individual animals and plaques when studying correlations between WSS and plaque composition.

Carotid plaque elasticity estimation using ultrasound elastography, MRI, and inverse FEA - A numerical feasibility study.

The material properties of atherosclerotic plaques govern the biomechanical environment, which is associated with rupture-risk. We investigated the feasibility of noninvasively estimating carotid plaque component material properties through simulating ultrasound (US) elastography and in vivo magnetic resonance imaging (MRI), and solving the inverse problem with finite element analysis. 2D plaque models were derived from endarterectomy specimens of nine patients. Nonlinear neo-Hookean models (tissue elasticity C1) were assigned to fibrous intima, wall (i.e., media/adventitia), and lipid-rich necrotic core. Finite element analysis was used to simulate clinical cross-sectional US strain imaging. Computer-simulated, single-slice in vivo MR images were segmented by two MR readers. We investigated multiple scenarios for plaque model elasticity, and consistently found clear separations between estimated tissue elasticity values. The intima C1 (160 kPa scenario) was estimated as 125.8 ± 19.4 kPa (reader 1) and 128.9 ± 24.8 kPa (reader 2). The lipid-rich necrotic core C1 (5 kPa) was estimated as 5.6 ± 2.0 kPa (reader 1) and 8.5 ± 4.5 kPa (reader 2). A scenario with a stiffer wall yielded similar results, while realistic US strain noise and rotating the models had little influence, thus demonstrating robustness of the procedure. The promising findings of this computer-simulation study stimulate applying the proposed methodology in a clinical setting.

Geometry-based pressure drop prediction in mildly diseased human coronary arteries.

Pressure drop (△p) estimations in human coronary arteries have several important applications, including determination of appropriate boundary conditions for CFD and estimation of fractional flow reserve (FFR). In this study a △p prediction was made based on geometrical features derived from patient-specific imaging data. Twenty-two mildly diseased human coronary arteries were imaged with computed tomography and intravascular ultrasound. Each artery was modelled in three consecutive steps: from straight to tapered, to stenosed, to curved model. CFD was performed to compute the additional △p in each model under steady flow for a wide range of Reynolds numbers. The correlations between the added geometrical complexity and additional △p were used to compute a predicted △p. This predicted △p based on geometry was compared to CFD results. The mean △p calculated with CFD was 855±666Pa. Tapering and curvature added significantly to the total △p, accounting for 31.4±19.0% and 18.0±10.9% respectively at Re=250. Using tapering angle, maximum area stenosis and angularity of the centerline, we were able to generate a good estimate for the predicted △p with a low mean but high standard deviation: average error of 41.1±287.8Pa at Re=250. Furthermore, the predicted △p was used to accurately estimate FFR (r=0.93). The effect of the geometric features was determined and the pressure drop in mildly diseased human coronary arteries was predicted quickly based solely on geometry. This pressure drop estimation could serve as a boundary condition in CFD to model the impact of distal epicardial vessels.

Uniaxial tensile testing approaches for characterisation of atherosclerotic plaques.

The pathological changes associated with the development of atherosclerotic plaques within arterial vessels result in significant alterations to the mechanical properties of the diseased arterial wall. There are several methods available to characterise the mechanical behaviour of atherosclerotic plaque tissue, and it is the aim of this paper to review the use of uniaxial mechanical testing. In the case of atherosclerotic plaques, there are nine studies that employ uniaxial testing to characterise mechanical behaviour. A primary concern regarding this limited cohort of published studies is the wide range of testing techniques that are employed. These differing techniques have resulted in a large variance in the reported data making comparison of the mechanical behaviour of plaques from different vasculatures, and even the same vasculature, difficult and sometimes impossible. In order to address this issue, this paper proposes a more standardised protocol for uniaxial testing of diseased arterial tissue that allows for better comparisons and firmer conclusions to be drawn between studies. To develop such a protocol, this paper reviews the acquisition and storage of the tissue, the testing approaches, the post-processing techniques and the stress-strain measures employed by each of the nine studies. Future trends are also outlined to establish the role that uniaxial testing can play in the future of arterial plaque mechanical characterisation.

The influence of axial image resolution on atherosclerotic plaque stress computations.

Biomechanical models are used extensively to study risk factors, such as peak stresses, for vulnerable atherosclerotic plaque rupture. Typically, 3D patient-specific arterial models are reconstructed by interpolating between cross sectional contour data which have a certain axial sampling, or image, resolution. The influence of the axial sampling resolution on computed stresses, as well as the comparison of 3D with 2D simulations, is quantified in this study. A set of histological data of four atherosclerotic human coronary arteries was used which were reconstructed in 3D with a high sampling (HS) and low sampling (LS) axial resolution, and 4 slices were treated separately for 2D simulations. Stresses were calculated using finite element analysis (FEA). High stresses were found in thin cap regions and regions of thin vessel walls, low stresses were found inside the necrotic cores and media and adventitia layers. Axial sampling resolution was found to have a minor effect on general stress distributions, peak plaque/cap stress locations and the relationship between peak cap stress and minimum cap thickness. Axial sampling resolution did have a profound influence on the error in computed magnitude of peak plaque/cap stresses (±15.5% for HS vs. LS geometries and ±24.0% for HS vs. 2D geometries for cap stresses). The findings of this study show that axial under sampling does not influence the qualitative stress distribution significantly but that high axially sampled 3D models are needed when accurate computation of peak stress magnitudes is required.

Initial stress in biomechanical models of atherosclerotic plaques.

Rupture of atherosclerotic plaques is the underlying cause for the majority of acute strokes and myocardial infarctions. Rupture of the plaque occurs when the stress in the plaque exceeds the strength of the material locally. Biomechanical stress analyses are commonly based on pressurized geometries, in most cases measured by in-vivo MRI. The geometry is therefore not stress-free. The aim of this study is to identify the effect of neglecting the initial stress state on the plaque stress distribution. Fifty 2D histological sections (7 patients, 9 diseased coronary artery segments), perfusion fixed at 100 mmHg, were segmented and finite element models were created. The Backward Incremental method was applied to determine the initial stress state and the zero-pressure state. Peak plaque and cap stresses were compared with and without initial stress. The effect of initial stress on the peak stress was related to the minimum cap thickness, maximum necrotic core thickness, and necrotic core angle. When accounting for initial stress, the general relations between geometrical features and peak cap stress remain intact. However, on a patient-specific basis, accounting for initial stress has a different effect on the absolute cap stress for each plaque. Incorporating initial stress may therefore improve the accuracy of future stress based rupture risk analyses for atherosclerotic plaques.

High shear stress influences plaque vulnerability Part of the data presented in this paper were published in Stroke 2007;38:2379-81.

Shear stress of the blood at the vessel wall plays an important role in many processes in the cardiovascular system primarily focused on the regulation of vessel lumen and wall dimensions. There is ample evidence that atherosclerotic plaques are generated at low shear stress regions in the cardiovascular system, while high shear stress regions are protected. In the course of plaque progression, advanced plaques start to encroach into the lumen, and thereby start to experience high shear stress at the endothelium. Until now the consequences of high shear stress working at the endothelium of an advanced plaque are unknown. As high shear stress influences tissue regression, we hypothesised that high shear stress can destabilise the plaque by cap weakening leading to ulceration. We investigated this hypothesis in a magnetic resonance imaging (MRI) dataset of a 67-year-old woman with a plaque in the carotid artery at baseline and an ulcer at ten-month follow-up. The lumen, plaque components (lipid/necrotic core, intraplaque haemorrhage) and ulcer were reconstructed three dimensionally and the geometry at baseline was used for shear stress calculation using computational fluid dynamics. Correlation of the change in plaque composition with the shear stress at baseline showed that the ulcer was generated exclusively at the high shear stress location. In this serial MRI study we found plaque ulceration at the high shear stress location of a protruding plaque in the carotid artery. Our data suggest that high shear stress influences plaque vulnerability and therefore may become a potential parameter for predicting future events. (Neth Heart J 2008;16:280-3.).

A numerical study on the influence of vulnerable plaque composition on intravascular thermography measurements.

Intracoronary thermography is a technique that measures lumen wall temperatures for vulnerable plaque detection. In this paper the influence of vulnerable plaque composition on lumen wall temperatures was studied numerically. Concerning the vulnerable plaque heat generation, the location of the heat source and its heat production were varied. Concerning the heat transfer, the thermal properties of the lipid core and the location of the vasa vasorum were studied. The heat source location was the main determinant of the lumen wall temperature distribution. The strongest effect was noted when the heat producing macrophages were located in the shoulder region leading to focal spots of higher temperature. The maximal lumen wall temperature was mainly determined by the heat production of the macrophages and the cooling effect of blood. The insulating properties of the lipid core increased lumen wall temperatures when the heat source was located in the cap and the presence of vasa vasorum lowered the temperatures. These results show that the lumen wall temperature distribution is influenced by vulnerable plaque composition and that intracoronary thermography techniques require a high spatial resolution. To be able to couple temperature measurements to plaque vulnerability, intracoronary thermography needs to be combined with an imaging modality.

The role of shear stress in the destabilization of vulnerable plaques and related therapeutic implications.

American Heart Association type IV plaques consist of a lipid core covered by a fibrous cap, and develop at locations of eccentric low shear stress. Vascular remodeling initially preserves the lumen diameter while maintaining the low shear stress conditions that encourage plaque growth. When these plaques eventually start to intrude into the lumen, the shear stress in the area surrounding the plaque changes substantially, increasing tensile stress at the plaque shoulders and exacerbating fissuring and thrombosis. Local biologic effects induced by high shear stress can destabilize the cap, particularly on its upstream side, and turn it into a rupture-prone, vulnerable plaque. Tensile stress is the ultimate mechanical factor that precipitates rupture and atherothrombotic complications. The shear-stress-oriented view of plaque rupture has important therapeutic implications. In this review, we discuss the varying mechanobiologic mechanisms in the areas surrounding the plaque that might explain the otherwise paradoxical observations and unexpected outcomes of experimental therapies.

The role of shear stress in the generation of rupture-prone vulnerable plaques.

Blood-flow-induced shear stress acting on the arterial wall is of paramount importance in vascular biology. Endothelial cells sense shear stress and largely control its value in a feedback-control loop by adapting the arterial dimensions to blood flow. Nevertheless, to allow for variations in arterial geometry, such as bifurcations, shear stress control is modified at certain eccentrically located sites to let it remain at near-zero levels. In the presence of risk factors for atherosclerosis, low shear stress contributes to local endothelial dysfunction and eccentric plaque build up, but normal-to-high shear stress is atheroprotective. Initially, lumen narrowing is prevented by outward vessel remodeling. Maintenance of a normal lumen and, by consequence, a normal shear stress distribution, however, prolongs local unfavorable low shear stress conditions and aggravates eccentric plaque growth. While undergoing such growth, eccentric plaques at preserved lumen locations experience increased tensile stress at their shoulders making them prone to fissuring and thrombosis. Consequent loss of the plaque-free wall by coverage with thrombus and new tissue may bring shear-stress-controlled lumen preservation to an end. This change causes shear stress to increase, which as a new condition may transform the lesion into a rupture-prone vulnerable plaque. We present a discussion of the role of shear stress, in setting the stage for the generation of rupture-prone, vulnerable plaques, and how this may be prevented.

Temperature distribution in atherosclerotic coronary arteries: influence of plaque geometry and flow (a numerical study).

Intravascular coronary thermography is a method that may detect vulnerable, atherosclerotic plaques and is currently evaluated in a clinical setting. Active macrophages or enzymatic heat releasing processes in vulnerable plaques may act as heat sources. To better understand the parameters of influence on thermographic measurements, numerical simulations have been performed on a model of a coronary artery segment containing a heat source. Heat source parameters and flow were varied to study their influence on temperatures at the lumen wall. Maximal temperature differences at the lumen wall increased when the source volume increased and they differ with the source geometry. The simulations showed that blood flow acts as a coolant to the lumen wall. Blood flow decreased maximal temperatures depending on the source geometry, source volume and the maximal flow velocity. Influence of flow was highest for circumferentially extended sources, up to a factor 3.7, and lowest for longitudinally extended sources, down to a factor 1.9. When cap thickness increased, maximal temperatures decreased and the influence of flow increased. This study shows that correct interpretation of intravascular thermographic measurements requires data on the flow and on the morphologic characteristics of the atherosclerotic plaque.