Low-density lipoprotein transport within a multi-layered arterial wall—Effect of the atherosclerotic plaque/stenosis

Impact of geometric and hemodynamic changes on a mechanobiological model of atherosclerosis

BACKGROUND AND OBJECTIVE

In this work, the analysis of the importance of hemodynamic updates on a mechanobiological model of atheroma plaque formation is proposed.

METHODS

For that, we use an idealized and axisymmetric model of carotid artery. In addition, the behavior of endothelial cells depending on hemodynamical changes is analyzed too. A total of three computational simulations are carried out and their results are compared: an uncoupled model and two models that consider the opposite behavior of endothelial cells caused by hemodynamic changes. The model considers transient blood flow using the Navier-Stokes equation. Plasma flow across the endothelium is determined with Darcy's law and the Kedem-Katchalsky equations, considering the three-pore model, which is also employed for the flow of substances across the endothelium. The behavior of the considered substances in the arterial wall is modeled with convection-diffusion-reaction equations, and the arterial wall is modeled as a hyperelastic Yeoh's material.

RESULTS

Significant variations are noted in both the morphology and stenosis ratio of the plaques when comparing the uncoupled model to the two models incorporating updates for geometry and hemodynamic stimuli. Besides, the phenomenon of double-stenosis is naturally reproduced in the models that consider both geometric and hemodynamical changes due to plaque growth, whereas it cannot be predicted in the uncoupled model.

CONCLUSIONS

The findings indicate that integrating the plaque growth model with geometric and hemodynamic settings is essential in determining the ultimate shape and dimensions of the carotid plaque.

Spatiotemporal changes of local hemodynamics and plaque components during atherosclerotic progression in rabbit

BACKGROUND AND OBJECTIVE

Recent evidence demonstrates that the atherogenic process is discontinuous. Our goal is to study changes of plaque components and local hemodynamics during atherosclerotic progression.

METHODS

The histological and immunohistochemical staining of high-fat diet rabbit aorta were evaluated at 0, 8, 10 and 12 weeks, respectively. In addition, the blood flow and LDL transport were simulated at the above four time points.

RESULTS

The plaque thickness at different characteristic regions increased at different rates. The collagen continued to increase, while the elastin, fibronectin, macrophages and smooth muscle cells increased first and then decreased. The relative surface LDL concentration decreased at 8 weeks, and then it increased first and decreased slightly. Meanwhile, the hemodynamic environment became better firstly at 8 weeks, then got slightly worse and lastly improved again.

CONCLUSIONS

The local hemodynamics and plaque components vary nonlinearly during atherosclerotic progression in rabbit aorta.

EFFECTS OF COATING PROPERTIES ON CONTROLLED DELIVERY FROM AN EMBEDDED DRUG-ELUTING STENT: A SIMULATION STUDY

The present investigation deals with the effects of biodegradable, biodurable and polymer-free coating of a stent on the release mechanism of the drug in a porous medium. The Brinkman equations for the interstitial fluid, the unsteady convection-diffusion-reaction equation for the transport of free drug in the tissue and the unsteady reaction equations for the bound as well as the internalized drug have been considered. In the coating, the transport of drug has been modeled as a diffusion process. Effects of different percentages of the embedment, convection and various coating properties of the stent on the transport of free drug, its retention and the internalization of the bound drug have been studied. Immersed Boundary Method (IBM) in the staggered grid formulation (IBM-MAC) has been used to tackle numerically the system of nonlinear governing equations. Simulated results predict the fastest release of drug from a biodegradable coating, but the averaged concentrations of all drug forms do reach a quasi-steady state in case of a biodurable coating irrespective of the degrees of embedment. Moreover, for all embedment levels of the stent, a biodegradable coating is superior to that of biodurable and polymer-free coating in the presence/absence of convection for larger times, but this superiority is lost for smaller times. Unlike biodurable coating, it is also predicted that the more the embedment level does not necessarily imply the more the effectiveness of delivery for biodegradable and polymer-free coatings of a stent.

Integrating blood cell mechanics, platelet adhesive dynamics and coagulation cascade for modelling thrombus formation in normal and diabetic blood

Normal haemostasis is an important physiological mechanism that prevents excessive bleeding during trauma, whereas the pathological thrombosis especially in diabetics leads to increased incidence of heart attacks and strokes as well as peripheral vascular events. In this work, we propose a new multiscale framework that integrates seamlessly four key components of blood clotting, namely transport of coagulation factors, coagulation kinetics, blood cell mechanics and platelet adhesive dynamics, to model the development of thrombi under physiological and pathological conditions. We implement this framework to simulate platelet adhesion due to the exposure of tissue factor in a three-dimensional microchannel. Our results show that our model can simulate thrombin-mediated platelet activation in the flowing blood, resulting in platelet adhesion to the injury site of the channel wall. Furthermore, we simulate platelet adhesion in diabetic blood, and our results show that both the pathological alterations in the biomechanics of blood cells and changes in the amount of coagulation factors contribute to the excessive platelet adhesion and aggregation in diabetic blood. Taken together, this new framework can be used to probe synergistic mechanisms of thrombus formation under physiological and pathological conditions, and open new directions in modelling complex biological problems that involve several multiscale processes.

Prediction of the development of coronary atherosclerotic plaques using computational modeling in 3D reconstructed coronary arteries

In this work we present a novel method for the prediction and generation of atherosclerotic plaques. This is performed in a two-step approach, by employing first a multilevel computational plaque growth model and second a correlation between the model's results and the 3D reconstructed follow-up plaques. In particular, computer tomography coronary angiography (CTCA) data and blood tests were collected from patients at two time points. Using the baseline data, the plaque growth is simulated using a multi-level computational model which includes: i) modeling of the blood flow dynamics, ii) modeling of low and high density lipoproteins and monocytes' infiltration in the arterial wall, and the species reactions during the atherosclerotic process, and iii) modeling of the arterial wall thickening. The correlation between the followup plaques and the simulated plaque density distribution resulted to the extraction of a threshold of the plaque density, that can be used to identify plaque areas.Clinical Relevance- The methodology presented in this work is a first step to the prediction of the plaque shape and location of patients with atherosclerosis and could be used as an additional tool for patient-specific risk stratification.

Modeling tau transport in the axon initial segment

By assuming that tau protein can be in seven kinetic states, we developed a model of tau protein transport in the axon and in the axon initial segment (AIS). Two separate sets of kinetic constants were determined, one in the axon and the other in the AIS. This was done by fitting the model predictions in the axon with experimental results and by fitting the model predictions in the AIS with the assumed linear increase of the total tau concentration in the AIS. The calibrated model was used to make predictions about tau transport in the axon and in the AIS. To the best of our knowledge, this is the first paper that presents a mathematical model of tau transport in the AIS. Our modeling results suggest that binding of free tau to microtubules creates a negative gradient of free tau in the AIS. This leads to diffusion-driven tau transport from the soma into the AIS. The model further suggests that slow axonal transport and diffusion-driven transport of tau work together in the AIS, moving tau anterogradely. Our numerical results predict an interplay between these two mechanisms: as the distance from the soma increases, the diffusion-driven transport decreases, while motor-driven transport becomes larger. Thus, the machinery in the AIS works as a pump, moving tau into the axon.

Pre-atherosclerotic flow and oncotically active solute transport across the arterial endothelium

Atherosclerosis starts with transmural (transwall) pressure-driven advective transport of blood-borne low-density lipoprotein (LDL) cholesterol across rare endothelial cell (EC) monolayer leaks into the arterial subendothelial intima (SI) wall layer where they can spread, bind to extracellular matrix and seed lesions. The local SI LDL concentration, which governs LDL's binding kinetics, depends on the overall diluting transmural liquid flow. Transmural pressures typically compress the SI at physiological pressures, which keeps this flow low. Nguyen et al. (2015) showed that aortic ECs express the water channel protein aquaporin-1 (AQP1) and the transEC (δP) portion of the transmural (ΔP) pressure difference drives, in parallel, water across AQP1s and plasma across interEC junctions. Since the lumen is isotonic, selective AQP1-mediated water flow should quickly render the ECs' lumen side hypertonic and the SI hypotonic; resulting transEC oncotic pressure differences, δπ, should oppose δP and quickly halt transEC flow. Yet Nguyen et al.'s (2015) transAQP1 flows persist for hours. To resolve this paradox, we extend our fluid filtration theory Joshi et al. (2015) to include mass transfer for oncotically active solutes like albumin. This addition nonlinearly couples mass transfer, fluid flow and wall mechanics. We simultaneously solve these problems at steady state. Surprisingly it finds that media layer filtration causes steady SI to exceed EC glycocalyx albumin concentration. Thus δπ reinforces rather than opposes δP, i.e., it sucks water from, rather than pushing water into the lumen from the SI. Endothelial AQP1s raise the overall driving force for flow across the EC above δP, most significantly at pressures too low to compress the SI, and they increase the ΔP needed for SI compression. This suggests the intriguing possibility that increasing EC AQP1 expression can raise this requisite compression pressure to physiological values. That is, increasing EC AQP1 may decompress the SI at physiological pressures, which would significantly increase SI thickness, flow and subsequently SI LDL dilution. This could retard LDL binding and delay preatherosclerotic lesion onset. The model also predicts that glycocalyx-degrading enzymes decrease overall transEC driving forces and thus lower, not raise, transmural water flux.

Computational model of stent-based delivery from a half-embedded two-layered coating

An attempt is made in the present investigation to develop a computational model for the purpose of studying the effect of interstitial flow in the porous media on the distribution of drug eluted from a half-embedded drug-eluting stent and its retention in the presence of two-layered coating of the stent. The transport of free drug inside the coatings is considered as an unsteady diffusion process while that in the tissue as an unsteady convection-diffusion-reaction process. The bound drug is governed by an unsteady reaction process only. Immersed boundary method (IBM) in the staggered grid formulation, popularly known as marker and cell (MAC) method, has been leveraged to tackle numerically the governing equations. This model highlights the benefits of consideration of two-layered coating and does predict underlying mechanism for better efficacy by tweaking the kinetics parameters. Comparisons are also made with the results available for stent-based delivery.

Effect of Interstitial Fluid Flow on Drug-Coated Balloon Delivery in a Patient-Specific Arterial Vessel with Heterogeneous Tissue Composition: A Simulation Study

Angioplasty with drug-coated balloons (DCBs) using excipients as drug carriers is emerging as a potentially viable strategy demonstrating clinical efficacy and proposing additional compliance for the treatment of obstructive vascular diseases. An attempt is made to develop an improved computational model where attention has been paid to the effect of interstitial flow, that is, plasma convection and internalization of bound drug. The present model is capable of capturing the phenomena of the transport of free drug and its retention, and also the internalization of drug in the process of endocytosis to atherosclerotic vessel of heterogeneous tissue composition comprising of healthy tissue, as well as regions of fibrous cap, fibro-fatty, calcified and necrotic core lesions. Image processing based on an unsupervised clustering technique is used for color-based segmentation of a patient-specific longitudinal image of atherosclerotic vessel obtained from intravascular ultrasound derived virtual histology. As the residence time of drug in a stent-based delivery within the arterial tissue is strongly influenced by convective forces, effect of interstitial fluid flow in case of DCB delivery can not be ruled out, and has been investigated by modeling it through unsteady Navier-Stokes equations. Transport of free drug is modeled by considering unsteady advection-reaction-diffusion process, while the bound drug, assuming completely immobilized in the tissue, by unsteady reaction process. The model also takes into account the internalization of drug through the process of endocytosis which gets degraded by the lysosomes and finally recycled into the extracellular fluid. All the governing equations representing the flow of interstitial fluid, the transport of free drug, the metabolization of free drug into bound phase and the process of internalization along with the physiologically realistic boundary and initial conditions are solved numerically using marker and cell method satisfying necessary stability criteria. Simulated results obtained predict that faster drug transfer promotes rapid saturation of binding sites despite perivascular wash out and the concentrations of all drug forms are modulated by the presence of interstitial flow. Such premier attempt of its kind would certainly be of great help in the optimization of therapeutic efficacy of drug.

Numerical simulation of haemodynamics and low-density lipoprotein transport in the rabbit aorta and their correlation with atherosclerotic plaque thickness

Two mechanisms of shear stress and mass transport have been recognized to play an important role in the development of localized atherosclerosis. However, their relationship and roles in atherogenesis are still obscure. It is necessary to investigate quantitatively the correlation among low-density lipoproteins (LDL) transport, haemodynamic parameters and plaque thickness. We simulated blood flow and LDL transport in rabbit aorta using computational fluid dynamics and evaluated plaque thickness in the aorta of a high-fat-diet rabbit. The numerical results show that regions with high luminal LDL concentration tend to have severely negative haemodynamic environments (HEs). However, for regions with moderately and slightly high luminal LDL concentration, the relationship between LDL concentration and the above haemodynamic indicators is not clear cut. Point-by-point correlation with experimental results indicates that severe atherosclerotic plaque corresponds to high LDL concentration and seriously negative HEs, less severe atherosclerotic plaque is related to either moderately high LDL concentration or moderately negative HEs, and there is almost no atherosclerotic plaque in regions with both low LDL concentration and positive HEs. In conclusion, LDL distribution is closely linked to blood flow transport, and the synergetic effects of luminal surface LDL concentration and wall shear stress-based haemodynamic indicators may determine plaque thickness.

Low-density lipoprotein transport through an arterial wall under hypertension - A model with time and pressure dependent fraction of leaky junction consistent with experiments

The influence of hypertension on low-density lipoproteins intake into the arterial wall is an important factor for understanding mechanisms of atherosclerosis. It has been experimentally observed that the increased pressure leads to the higher level of the LDL inside the wall. In this paper we attempt to construct a model of the LDL transport which reproduces quantitatively experimental outcomes. We supplement the well-known four-layer arterial wall model to include two pressure induced effects: the compression of the intima tissue and the increase of the fraction of leaky junctions. We demonstrate that such model can reach the very good agreement with experimental data.

Prediction of Atherosclerotic Plaque Development in an In Vivo Coronary Arterial Segment Based on a Multilevel Modeling Approach

OBJECTIVE

The aim of this study is to explore major mechanisms of atherosclerotic plaque growth, presenting a proof-of-concept numerical model.

METHODS

To this aim, a human reconstructed left circumflex coronary artery is utilized for a multilevel modeling approach. More specifically, the first level consists of the modeling of blood flow and endothelial shear stress (ESS) computation. The second level includes the modeling of low-density lipoprotein (LDL) and high-density lipoprotein and monocytes transport through the endothelial membrane to vessel wall. The third level comprises of the modeling of LDL oxidation, macrophages differentiation, and foam cells formation. All modeling levels integrate experimental findings to describe the major mechanisms that occur in the arterial physiology. In order to validate the proposed approach, we utilize a patient specific scenario by comparing the baseline computational results with the changes in arterial wall thickness, lumen diameter, and plaque components using follow-up data.

RESULTS

The results of this model show that ESS and LDL concentration have a good correlation with the changes in plaque area [R² = 0.365 (P = 0.029, adjusted R² = 0.307) and R² = 0.368 (P = 0.015, adjusted R² = 0.342), respectively], whereas the introduction of the variables of oxidized LDL, macrophages, and foam cells as independent predictors improves the accuracy in predicting regions potential for atherosclerotic plaque development [R² = 0.847 (P = 0.009, adjusted R² = 0.738)].

CONCLUSION

Advanced computational models can be used to increase the accuracy to predict regions which are prone to plaque development.

SIGNIFICANCE

Atherosclerosis is one of leading causes of death worldwide. For this purpose computational models have to be implemented to predict disease progression.

Analysis of non-Newtonian effects on Low-Density Lipoprotein accumulation in an artery

In this work, non-Newtonian effects on Low-Density Lipoprotein (LDL) transport across an artery are analyzed with a multi-layer model. Four rheological models (Carreau, Carreau-Yasuda, power-law and Newtonian) are used for the blood flow through the lumen. For the non-Newtonian cases, the arterial wall is modeled with a generalized momentum equation. Convection-diffusion equation is used for the LDL transport through the lumen, while Staverman-Kedem-Katchalsky, combined with porous media equations, are used for the LDL transport through the wall. Results are presented in terms of filtration velocity, Wall Shear Stresses (WSS) and concentration profiles. It is shown that non-Newtonian effects on mass transport are negligible for a healthy intramural pressure value. Non-Newtonian effects increase slightly with intramural pressure, but Newtonian assumption can still be considered reliable. Effects of arterial size are also analyzed, showing that Newtonian assumption can be considered valid for both medium and large arteries, in predicting LDL deposition. Finally, non-Newtonian effects are also analyzed for an aorta-common iliac bifurcation, showing that Newtonian assumption is valid for mass transport at low Reynolds numbers. At a high Reynolds number, it has been shown that a non-Newtonian fluid model can have more impact due to the presence of flow recirculation.

Relative Effects of Fluid Oscillations and Nutrient Transport in the In Vitro Growth of Valvular Tissues

Engineered valvular tissues are cultured dynamically, and involve specimen movement. We previously demonstrated that oscillatory shear stresses (OSS) under combined steady flow and specimen cyclic flexure (flex-flow) promote tissue formation. However, localized efficiency of specimen mass transport is also important in the context of cell viability within the growing tissues. Here, we investigated the delivery of two essential species for cell survival, glucose and oxygen, to 3-dimensional (3D) engineered valvular tissues. We applied a convective-diffusive model to characterize glucose and oxygen mass transport with and without valve-like specimen flexural movement. We found the mass transport effects for glucose and oxygen to be negligible for scaffold porosities typically present during in vitro experiments and non-essential unless the porosity was unusually low (<40%). For more typical scaffold porosities (75%) however, we found negligible variation in the specimen mass fraction of glucose and oxygen in both non-moving and moving constructs (p > 0.05). Based on this result, we conducted an experiment using bone marrow stem cell (BMSC)-seeded scaffolds under Pulsatile flow-alone states to permit OSS without any specimen movement. BMSC-seeded specimen collagen from the pulsatile flow and flex-flow environments were subsequently found to be comparable (p > 0.05) and exhibited some gene expression similarities. We conclude that a critical magnitude of fluid-induced, OSS created by either pulsatile flow or flex-flow conditions, particularly when the oscillations are physiologically-relevant, is the direct, principal stimulus that promotes engineered valvular tissues and its phenotype, whereas mass transport benefits derived from specimen movement are minimal.

Low-density lipoprotein accumulation within a carotid artery with multilayer elastic porous wall: fluid-structure interaction and non-Newtonian considerations

Low-density lipoprotein (LDL), which is recognized as bad cholesterol, typically has been regarded as a main cause of atherosclerosis. LDL infiltration across arterial wall and subsequent formation of Ox-LDL could lead to atherogenesis. In the present study, combined effects of non-Newtonian fluid behavior and fluid-structure interaction (FSI) on LDL mass transfer inside an artery and through its multilayer arterial wall are examined numerically. Navier-Stokes equations for the blood flow inside the lumen and modified Darcy's model for the power-law fluid through the porous arterial wall are coupled with the equations of mass transfer to describe LDL distributions in various segments of the artery. In addition, the arterial wall is considered as a heterogeneous permeable elastic medium. Thus, elastodynamics equation is invoked to examine effects of different wall elasticity on LDL distribution in the artery. Findings suggest that non-Newtonian behavior of filtrated plasma within the wall enhances LDL accumulation meaningfully. Moreover, results demonstrate that at high blood pressure and due to the wall elasticity, endothelium pores expand, which cause significant variations on endothelium physiological properties in a way that lead to higher LDL accumulation. Additionally, results describe that under hypertension, by increasing angular strain, endothelial junctions especially at leaky sites expand more dramatic for the high elastic model, which in turn causes higher LDL accumulation across the intima layer and elevates atherogenesis risk.

Image-based computational simulation of sub-endothelial LDL accumulation in a human right coronary artery

Accumulation of low density lipoproteins (LDL) in the vessel wall is suggested as the initiator of atherosclerosis and coronary stenosis. This process is associated with the performance of endothelium layer that regulates entering of macromolecules to the vessel wall. Therefore, the present study aims to investigate sub-endothelial accumulation of LDL molecules in a coronary tree and predict atherosclerosis prone sites. Non-Newtonian blood flow is simulated for normal and hypertensive conditions through the lumen of a right coronary artery reconstructed from computed tomography (CT) images. A three-pore model is implemented as the endothelium boundary condition and hence, plasma flow and LDL transport are simulated within the arterial wall. Based on the pore model, endothelium pathways divide into normal junctions, vesicles and leaky junctions. Most of LDL molecules pass through the leaky junctions that arise at locations with low wall shear stress (WSS). Results indicate that increase in the number of leaky junctions at branch points with low WSS can lead to both elevated levels of sub-endothelial LDL accumulation and atherosclerosis risk. Findings reveal that at the branch points with disturbed flow, sub-endothelial concentration of LDL for the hypertensive condition is higher than the normal condition, however for the rest of regions with uniform geometry and unidirectional flow, this is reversed. Comparisons of non-Newtonian and Newtonian flows show mean increases of 34% and 13% in the sub-endothelial concentrations of Newtonian flows during the normal and hypertensive conditions, respectively.

Effect of Transmural Transport Properties on Atheroma Plaque Formation and Development

We propose a mathematical model of atheroma plaque initiation and early development in coronary arteries using anisotropic transmural diffusion properties. Our current approach is on the process on plaque initiation and intimal thickening rather than in severe plaque progression and rupture phenomena. The effect of transport properties, in particular the anisotropy of diffusion properties of the artery, on plaque formation and development is investigated using the proposed mathematical model. There is not a strong influence of the anisotropic transmural properties on LDL, SMCs and collagen distribution and concentrations along the artery. On the contrary, foam cells distribution strongly depends on the value of the radial diffusion coefficient of the substances [Formula: see text] and the ratio [Formula: see text]. Decreasing [Formula: see text] or diffusion coefficients ratio means a higher concentration of the foam cells close to the intima. Due to the fact that foam cells concentration is associated to the necrotic core formation, the final distribution of foam cells is critical to evolve into a vulnerable or fibrotic plaque.

Modelling and simulation of low-density lipoprotein transport through multi-layered wall of an anatomically realistic carotid artery bifurcation

A high concentration of low-density lipoprotein (LDL) is recognized as one of the principal risk factors for development of atherosclerosis. This paper reports on modelling and simulations of the coupled mass (LDL concentration) and momentum transport through the arterial lumen and the multi-layered arterial wall of an anatomically realistic carotid bifurcation. The mathematical model includes equations for conservation of mass, momentum and concentration, which take into account a porous layer structure, the biological membranes and reactive source/sink terms in different layers of the arterial wall, as proposed in Yang & Vafai (2006). A four-layer wall model of an arterial wall with constant thickness is introduced and initially tested on a simple cylinder geometry where realistic layer properties are specified. Comparative assessment with previously published results demonstrated proper implementation of the mathematical model. Excellent agreement for the velocity and LDL concentration distributions in the arterial lumen and in the artery wall are obtained. Then, an anatomically realistic carotid artery bifurcation is studied. This is the main novelty of the presented research. We find a strong dependence between underlying blood flow pattern (and consequently the wall shear stress distributions) and the uptake of the LDL concentration in the artery wall. The radial dependency of interactions between the diffusion, convection and chemical reactions within the multi-layered artery wall is crucial for accurate predictions of the LDL concentration in the media. It is shown that a four-layer wall model produced qualitatively good agreement with the experimental results of Meyer et al. (1996) in predicting levels of LDL within the media of a rabbit aorta under identical transmural pressure conditions. Finally, it is demonstrated that the adopted model represents a good initial platform for future numerical investigations of the initial stage of atherosclerosis for patient-specific geometries.

Mathematical modelling of atheroma plaque formation and development in coronary arteries

Atherosclerosis is a vascular disease caused by inflammation of the arterial wall, which results in the accumulation of low-density lipoprotein (LDL) cholesterol, monocytes, macrophages and fat-laden foam cells at the place of the inflammation. This process is commonly referred to as plaque formation. The evolution of the atherosclerosis disease, and in particular the influence of wall shear stress on the growth of atherosclerotic plaques, is still a poorly understood phenomenon. This work presents a mathematical model to reproduce atheroma plaque growth in coronary arteries. This model uses the Navier-Stokes equations and Darcy's law for fluid dynamics, convection-diffusion-reaction equations for modelling the mass balance in the lumen and intima, and the Kedem-Katchalsky equations for the interfacial coupling at membranes, i.e. endothelium. The volume flux and the solute flux across the interface between the fluid and the porous domains are governed by a three-pore model. The main species and substances which play a role in early atherosclerosis development have been considered in the model, i.e. LDL, oxidized LDL, monocytes, macrophages, foam cells, smooth muscle cells, cytokines and collagen. Furthermore, experimental data taken from the literature have been used in order to physiologically determine model parameters. The mathematical model has been implemented in a representative axisymmetric geometrical coronary artery model. The results show that the mathematical model is able to qualitatively capture the atheroma plaque development observed in the intima layer.