Mathematical modeling of atherosclerotic plaque destabilization: Role of neovascularization and intraplaque hemorrhage

Observational studies have identified angiogenesis from the adventitial vasa vasorum and intraplaque hemorrhage (IPH) as critical factors in atherosclerotic plaque progression and destabilization. Here we propose a mathematical model incorporating intraplaque neovascularization and hemodynamic calculation with plaque destabilization for the quantitative evaluation of the role of neoangiogenesis and IPH in the vulnerable atherosclerotic plaque formation. An angiogenic microvasculature is generated by two-dimensional nine-point discretization of endothelial cell proliferation and migration from the vasa vasorum. Three key cells (endothelial cells, smooth muscle cells and macrophages) and three key chemicals (vascular endothelial growth factors, extracellular matrix and matrix metalloproteinase) are involved in the plaque progression model, and described by the reaction–diffusion partial differential equations. The hemodynamic calculation of the microcirculation on the generated microvessel network is carried out by coupling the intravascular, interstitial and transvascular flow. The plasma concentration in the interstitial domain is defined as the description of IPH area according to the diffusion and convection with the interstitial fluid flow, as well as the extravascular movement across the leaky vessel wall. The simulation results demonstrate a series of pathophysiological phenomena during the vulnerable progression of an atherosclerotic plaque, including the expanding necrotic core, the exacerbated inflammation, the high microvessel density (MVD) region at the shoulder areas, the transvascular flow through the capillary wall and the IPH. The important role of IPH in the plaque destabilization is evidenced by simulations with varied model parameters. It is found that the IPH can significantly speed up the plaque vulnerability by increasing necrotic core and thinning fibrous cap. In addition, the decreased MVD and vessel permeability may slow down the process of plaque destabilization by reducing the IPH dramatically. We envision that the present model and its future advances can serve as a valuable theoretical platform for studying the dynamic changes in the microenvironment during the plaque destabilization