A new correlation for inclusion of leaky junctions in macroscopic modeling of atherosclerotic lesion initiation

NUMERICAL MODELING OF ATHEROSCLEROSIS LESION EVOLUTION IN TIME I. INITIAL STAGE OF THE DISEASE

Atherosclerosis is one of the main causes of death in the developed world. The disease, which is an inflammatory disease, has been the focus of many studies. A few studies attempted to model atherosclerosis lesion development mathematically while no attention has been paid to the multistage nature of the disease. The present study provides a mathematical model for atherosclerosis evolution by focusing on the inflammatory responses of the initial stage of the disease. In the model, the inflammatory response in type I lesion, which includes endothelium dysfunction, LDL oxidation, monocytes entry, foam cell formation and intima property changes, are coupled with the transport equations of blood and LDL in lumen and arterial wall. The innovation of the model is determination of the duration of the initial stage of lesion propagation for a specific patient while the presence of leaky junction in endothelial layer and LDL oxidation in the intima layer are considered. The greatest advantage of the study in comparison with previous studies is to provide a model for the initiating stage of the atherosclerosis development so that a more precise result of the disease evolution is obtained.

Mass transport of low density lipoprotein in reconstructed hemodynamic environments of human carotid arteries: the role of volume and solute flux through the endothelium

The accumulation of low density lipoprotein (LDL) in the arterial intima is a critical step in the initiation and progression of atheromatous lesions. In this study we examine subject-specific LDL transport into the intima of carotid bifurcations in three human subjects using a three-pore model for LDL mass transfer. Subject-specific carotid artery computational models were derived using magnetic resonance imaging (MRI) to obtain the geometry and phase-contract MRI (PC-MRI) to acquire pulsatile inflow and outflow boundary conditions for each subject. The subjects were selected to represent a wide range of anatomical configurations and different stages of atherosclerotic development from mild to moderate intimal thickening. A fluid-solid interaction (FSI) model was implemented in the computational fluid dynamics (CFD) approach in order to consider the effects of a compliant vessel on wall shear stress (WSS). The WSS-dependent response of the endothelium to LDL mass transfer was modeled by multiple pathways to include the contributions of leaky junctions, normal junctions, and transcytosis to LDL solute and plasma volume flux from the lumen into the intima. Time averaged WSS (TAWSS) over the cardiac cycle was computed to represent the spatial WSS distribution, and wall thickness (WTH) was determined from black blood MRI (BBMRI) so as to visualize intimal thickening patterns in the bifurcations. The regions which are exposed to low TAWSS correspond to elevated WTH and higher mass and volume flux via the leaky junctions. In all subjects, the maximum LDL solute flux was observed to be immediately downstream of the stenosis, supporting observations that existing atherosclerotic lesions tend to progress in the downstream direction of the stenosis.