Mathematical model of intimal thickening in atherosclerosis: vessel stenosis as a free boundary problem

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.

Blood Lipoproteins Shape the Phenotype and Lipid Content of Early Atherosclerotic Lesion Macrophages: A Dual-Structured Mathematical Model

Macrophages in atherosclerotic lesions exhibit a spectrum of behaviours or phenotypes. The phenotypic distribution of monocyte-derived macrophages (MDMs), its correlation with MDM lipid content, and relation to blood lipoprotein densities are not well understood. Of particular interest is the balance between low density lipoproteins (LDL) and high density lipoproteins (HDL), which carry bad and good cholesterol respectively. To address these issues, we have developed a mathematical model for early atherosclerosis in which the MDM population is structured by phenotype and lipid content. The model admits a simpler, closed subsystem whose analysis shows how lesion composition becomes more pathological as the blood density of LDL increases relative to the HDL capacity. We use asymptotic analysis to derive a power-law relationship between MDM phenotype and lipid content at steady-state. This relationship enables us to understand why, for example, lipid-laden MDMs have a more inflammatory phenotype than lipid-poor MDMs when blood LDL lipid density greatly exceeds HDL capacity. We show further that the MDM phenotype distribution always attains a local maximum, while the lipid content distribution may be unimodal, adopt a quasi-uniform profile or decrease monotonically. Pathological lesions exhibit a local maximum in both the phenotype and lipid content MDM distributions, with the maximum at an inflammatory phenotype and near the lipid content capacity respectively. These results illustrate how macrophage heterogeneity arises in early atherosclerosis and provide a framework for future model validation through comparison with single-cell RNA sequencing data.

A biomathematical model of atherosclerosis in mice

Atherosclerosis is one of the leading causes of death worldwide. Biomathematical modelling of the underlying disease and therapy processes might be a useful aid to develop and improve preventive and treatment concepts of atherosclerosis. We here propose a biomathematical model of murine atherosclerosis under different diet and treatment conditions including lipid modulating compound and antibiotics. The model is derived by translating known biological mechanisms into ordinary differential equations and by assuming appropriate response kinetics to the applied interventions. We explicitly describe the dynamics of relevant immune cells and lipid species in atherosclerotic lesions including the degree of blood vessel occlusion due to growing plaques. Unknown model parameters were determined by fitting the predictions of model simulations to time series data derived from mice experiments. Parameter fittings resulted in a good agreement of model and data for all 13 experimental scenarios considered. The model can be used to predict the outcome of alternative treatment schedules of combined antibiotic, immune modulating, and lipid lowering agents under high fat or normal diet. We conclude that we established a comprehensive biomathematical model of atherosclerosis in mice. We aim to validate the model on the basis of further experimental data.

WELL-POSEDNESS OF A MATHEMATICAL MODEL OF DIABETIC ATHEROSCLEROSIS WITH ADVANCED GLYCATION END-PRODUCTS

Atherosclerosis is a leading cause of death worldwide; it emerges as a result of multiple dynamical cell processes including hemodynamics, endothelial damage, innate immunity and sterol biochemistry. Making matters worse, nearly 463 million people have diabetes, which increases atherosclerosis-related inflammation, diabetic patients are twice as likely to have a heart attack or stroke. The pathophysiology of diabetic vascular disease is generally understood. Dyslipidemia with increased levels of atherogenic LDL, hyperglycemia, oxidative stress and increased inflammation are factors that increase the risk and accelerate development of atherosclerosis. In a recent paper [53], we have developed mathematical model that includes the effect of hyperglycemia and insulin resistance on plaque growth. In this paper, we propose a more comprehensive mathematical model for diabetic atherosclerosis which include more variables; in particular it includes the variable for Advanced Glycation End-Products (AGEs)concentration. Hyperglycemia trigger vascular damage by forming AGEs, which are not easily metabolized and may accelerate the progression of vascular disease in diabetic patients. The model is given by a system of partial differential equations with a free boundary. We also establish local existence and uniqueness of solution to the model. The methodology is to use Hanzawa transformation to reduce the free boundary to a fixed boundary and reduce the system of partial differential equations to an abstract evolution equation in Banach spaces, and apply the theory of analytic semigroup.

Mathematical modeling of plaque progression and associated microenvironment: How far from predicting the fate of atherosclerosis?

Mathematical modeling contributes to pathophysiological research of atherosclerosis by helping to elucidate mechanisms and by providing quantitative predictions that can be validated. In turn, the complexity of atherosclerosis is well suited to quantitative approaches as it provides challenges and opportunities for new developments of modeling. In this review, we summarize the current 'state of the art' on the mathematical modeling of the effects of biomechanical factors and microenvironmental factors on the plaque progression, and its potential help in prediction of plaque development. We begin with models that describe the biomechanical environment inside and outside the plaque and its influence on its growth and rupture. We then discuss mathematical models that describe the dynamic evolution of plaque microenvironmental factors, such as lipid deposition, inflammation, smooth muscle cells migration and intraplaque hemorrhage, followed by studies on plaque growth and progression using these modelling approaches. Moreover, we present several key questions for future research. Mathematical models can complement experimental and clinical studies, but also challenge current paradigms, redefine our understanding of mechanisms driving plaque vulnerability and propose future potential direction in therapy for cardiovascular disease.

Effects of the Haemodynamic Stimulus on the Location of Carotid Plaques Based on a Patient-Specific Mechanobiological Plaque Atheroma Formation Model

In this work, we propose a mechanobiological atheroma growth model modulated by a new haemodynamic stimulus. To test this model, we analyse the development of atheroma plaques in patient-specific bifurcations of carotid arteries for a total time of 30 years. In particular, eight geometries (left or right carotid arteries) were segmented from clinical images and compared with the solutions obtained computationally to validate the model. The influence of some haemodynamical stimuli on the location and size of plaques is also studied. Plaques predicted by the mechanobiological models using the time average wall shear stress (TAWSS), the oscillatory shear index (OSI) and a new index proposed in this work are compared. The new index predicts the shape index of the endothelial cells as a combination of TAWSS and OSI values and was fitted using data from the literature. The mechanobiological model represents an evolution of the one previously proposed by the authors. This model uses Navier-Stokes equations to simulate blood flow along the lumen in the transient mode. It also employs Darcy's law and Kedem-Katchalsky equations for plasma and substance flow across the endothelium using the three-pore model. The mass balances of all the substances that have been considered in the model are implemented by convection-diffusion-reaction equations, and finally the growth of the plaques has been computed. The results show that by using the new mechanical stimulus proposed in this study, prediction of plaques is, in most cases, better than only using TAWSS or OSI with a minimal and maximal errors on stenosis ratio of 2.77 and 32.89 %, respectively. However, there are a few geometries in which haemodynamics cannot predict the location of plaques, and other biological or genetic factors would be more relevant than haemodynamics. In particular, the model predicts correctly eleven of the fourteen plaques presented in all the geometries considered. Additionally, a healthy geometry has been computed to check that plaque is not developed with the model in this case.

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.

A multiphase model of growth factor-regulated atherosclerotic cap formation

Atherosclerosis is characterised by the growth of fatty plaques in the inner artery wall. In mature plaques, vascular smooth muscle cells (SMCs) are recruited from adjacent tissue to deposit a collagenous cap over the fatty plaque core. This cap isolates the thrombogenic plaque content from the bloodstream and prevents the clotting cascade that leads to myocardial infarction or stroke. Despite the protective role of the cap, the mechanisms that regulate cap formation and maintenance are not well understood. It remains unclear why some caps become stable, while others become vulnerable to rupture. We develop a multiphase PDE model with non-standard boundary conditions to investigate collagen cap formation by SMCs in response to diffusible growth factor signals from the endothelium. Platelet-derived growth factor stimulates SMC migration, proliferation and collagen degradation, while transforming growth factor (TGF)-[Formula: see text] stimulates SMC collagen synthesis and inhibits collagen degradation. The model SMCs respond haptotactically to gradients in the collagen phase and have reduced rates of migration and proliferation in dense collagenous tissue. The model, which is parameterised using in vivo and in vitro experimental data, reproduces several observations from plaque growth in mice. Numerical and analytical results demonstrate that a stable cap can be formed by a relatively small SMC population and emphasise the critical role of TGF-[Formula: see text] in effective cap formation. These findings provide unique insight into the mechanisms that may lead to plaque destabilisation and rupture. This work represents an important step towards the development of a comprehensive in silico plaque model.

Mathematical Modelling and Simulation of Atherosclerosis Formation and Progress: A Review

Cardiovascular disease (CVD) is a major threat to human health since it is the leading cause of death in western countries. Atherosclerosis is a type of CVD related to hypertension, diabetes, high levels of cholesterol, smoking, oxidative stress, and age. Atherosclerosis primarily occurs in medium and large arteries, such as coronary and the carotid artery and, in particular, at bifurcations and curvatures. Atherosclerosis is compared to an inflammatory disease where a thick, porous material comprising cholesterol fat, saturated sterols, proteins, fatty acids, calcium etc., is covered by an endothelial membrane and a fragile fibrous tissue which makes atheromatic plaque prone to rupture that could lead to the blockage of the artery due to the released plaque material. Despite the great progress achieved, the nature of the disease is not fully understood. This paper reviews the current state of modelling of all levels of atherosclerosis formation and progress and discusses further challenges in atherosclerosis modelling. The objective is to pave a way towards more precise computational tools to predict and eventually reengineer the fate of atherosclerosis.

Coupled Modeling of Lipid Deposition, Inflammatory Response and Intraplaque Angiogenesis in Atherosclerotic Plaque

We propose a multiphysical mathematical model by fully coupling lipid deposition, monocytes/macrophages recruitment and angiogenesis to investigate the pathophysiological responses of an atherosclerotic plaque to the dynamic changes in the microenvironment. The time evolutions of cellular (endothelial cells, macrophages, smooth muscle cells, etc.) and acellular components (low density lipoprotein, proinflammatory cytokines, extravascular plasma concentration, etc.) within the plaque microenvironment are assessed quantitatively. The thickening of the intima, the distributions of the lipid and inflammatory factors, and the intraplaque hemorrhage show a qualitative consistency with the MRI and histology data. Models with and without angiogenesis are compared to demonstrate the important role of neovasculature in the accumulation of blood-borne components in the atherosclerotic lesion by extravasation from the leaky vessel wall, leading to the formation of a lipid core and an inflammatory microenvironment, which eventually promotes plaque destabilization. This model can serve as a theoretical platform for the investigation of the pathological mechanisms of plaque progression and may contribute to the optimal design of atherosclerosis treatment strategies, such as lipid-lowering or anti-angiogenetic therapies.

A two-phase model of early fibrous cap formation in atherosclerosis

Atherosclerotic plaque growth is characterised by chronic, non-resolving inflammation that promotes the accumulation of cellular debris and extracellular fat in the inner artery wall. This material is highly thrombogenic, and plaque rupture can lead to the formation of blood clots that occlude major arteries and cause myocardial infarction or stroke. In advanced plaques, vascular smooth muscle cells (SMCs) are recruited from deeper in the artery wall to synthesise a cap of fibrous tissue that stabilises the plaque and sequesters the thrombogenic plaque content from the bloodstream. The fibrous cap provides crucial protection against the clinical consequences of atherosclerosis, but the mechanisms of cap formation are poorly understood. In particular, it is unclear why certain plaques become stable and robust while others become fragile and dangerously vulnerable to rupture. We develop a multiphase model with non-standard boundary conditions to investigate early fibrous cap formation in the atherosclerotic plaque. The model is parameterised using data from a range of in vitro and in vivo studies, and includes highly nonlinear mechanisms of SMC proliferation and migration in response to an endothelium-derived chemical signal. We demonstrate that the model SMC population naturally evolves towards a steady-state, and predict a rate of cap formation and a final plaque SMC content consistent with experimental observations in mice. Parameter sensitivity simulations show that SMC proliferation makes a limited contribution to cap formation, and demonstrate that stable cap formation relies primarily on a critical balance between the rates of SMC recruitment to the plaque, chemotactic SMC migration within the plaque and SMC loss by apoptosis or phenotype change. This model represents the first detailed in silico study of fibrous cap formation in atherosclerosis, and establishes a multiphase modelling framework that can be readily extended to investigate many other aspects of plaque development.

A biochemical and mechanical model of injury-induced intimal thickening

In this paper, we investigate an axisymmetric model of intimal thickening using hyperelasticity theory. Our model describes the growth of the arterial intima due to cell proliferation which, in turn, is driven by the release of a cytokine such as platelet-derived growth factor (PDGF). With the growth rate tied to both local stress and the local concentration of PDGF, we derive a quadruple free boundary problem with different regions of the vessel wall characterized by different homeostatic stress. We compare our model predictions to rabbit and rodent models of atherosclerosis and find that in order to achieve the growth rates reported in the experiments, growth must be mainly cytokine induced rather than stress induced. Our model is also able to reproduce Glagov remodelling where, as a vessel becomes more diseased, the lumen expands before rapidly contracting.

Patient-Specific, Multi-Scale Modeling of Neointimal Hyperplasia in Vein Grafts

Neointimal hyperplasia is amongst the major causes of failure of bypass grafts. The disease progression varies from patient to patient due to a range of different factors. In this paper, a mathematical model will be used to understand neointimal hyperplasia in individual patients, combining information from biological experiments and patient-specific data to analyze some aspects of the disease, particularly with regard to mechanical stimuli due to shear stresses on the vessel wall. By combining a biochemical model of cell growth and a patient-specific computational fluid dynamics analysis of blood flow in the lumen, remodeling of the blood vessel is studied by means of a novel computational framework. The framework was used to analyze two vein graft bypasses from one patient: a femoro-popliteal and a femoro-distal bypass. The remodeling of the vessel wall and analysis of the flow for each case was then compared to clinical data and discussed as a potential tool for a better understanding of the disease. Simulation results from this first computational approach showed an overall agreement on the locations of hyperplasia in these patients and demonstrated the potential of using new integrative modeling tools to understand disease progression.

Bifurcation and dynamics in a mathematical model of early atherosclerosis: How acute inflammation drives lesion development

We present here a mathematical model describing the primary mechanisms that drive the early stages of atherosclerosis. This involves the interactions between modified low density lipoprotein (LDL), monocytes/macrophages, cytokines and foam cells. This model suggests that there is an initial inflammatory phase associated with atherosclerotic lesion development and a longer, quasi-static process of plaque development inside the arterial wall that follows the initial transient. We will show results that show how different LDL concentrations in the blood stream and different immune responses can affect the development of a plaque. Through numerical bifurcation analysis, we show the existence of a fold bifurcation when the flux of LDL from the blood is sufficiently high. By analysing the model presented in this paper, we gain a greater insight into this inflammatory response qualitatively and quantitatively.

Using a Systems Pharmacology Approach to Study the Effect of Statins on the Early Stage of Atherosclerosis in Humans

More than 100,000 people have participated in controlled trials of statins (lowering cholesterol drugs) since the introduction of lovastatin in the 1980s. Meta-analyses of this data have shown that statins have a beneficial effect on treated groups compared to control groups, reducing cardiovascular risk. Inhibiting the HMG-CoA reductase in the liver, statins can reduce cholesterol levels, thus reducing LDL levels in circulation. Published data from intravascular ultrasound studies (IVUS) was used in this work to develop and validate a unique integrative system model; this consisted of analyzing control groups from two randomized controlled statins trials (24/97 subjects respectively), one treated group (40 subjects, simvastatin trial), and 27 male subjects (simvastatin, pharmacokinetic study). The model allows to simulate the pharmacokinetics of statins and its effect on the dynamics of lipoproteins (e.g., LDL) and the inflammatory pathway while simultaneously exploring the effect of flow-related variables (e.g., wall shear stress) on atherosclerosis progression.

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.