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

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 Lipid-Structured Model of Atherosclerosis with Macrophage Proliferation

Atherosclerotic plaques are fatty deposits that form in the walls of major arteries and are one of the major causes of heart attacks and strokes. Macrophages are the main immune cells in plaques and macrophage dynamics influence whether plaques grow or regress. Macrophage proliferation is a key process in atherosclerosis, particularly in the development of mid-stage plaques, but very few mathematical models include proliferation. In this paper we reframe the lipid-structured model of Ford et al. (J Theor Biol 479:48-63, 2019. https://doi.org/10.1016/j.jtbi.2019.07.003 ) to account for macrophage proliferation. Proliferation is modelled as a non-local decrease in the lipid structural variable. Steady state analysis indicates that proliferation assists in reducing eventual necrotic core lipid content and spreads the lipid load of the macrophage population amongst the cells. The contribution of plaque macrophages from proliferation relative to recruitment from the bloodstream is also examined. The model suggests that a more proliferative plaque differs from an equivalent (defined as having the same lipid content and cell numbers) recruitment-dominant plaque in the way lipid is distributed amongst the macrophages. The macrophage lipid distribution of an equivalent proliferation-dominant plaque is less skewed and exhibits a local maximum near the endogenous lipid content.

Development of physiologically-informed computational coronary artery plaques for use in virtual imaging trials

BACKGROUND

As a leading cause of death, worldwide, cardiovascular disease is of great clinical importance. Among cardiovascular diseases, coronary artery disease (CAD) is a key contributor, and it is the attributed cause of death for 10% of all deaths annually. The prevalence of CAD is commensurate with the rise in new medical imaging technologies intended to aid in its diagnosis and treatment. The necessary clinical trials required to validate and optimize these technologies require a large cohort of carefully controlled patients, considerable time to complete, and can be prohibitively expensive. A safer, faster, less expensive alternative is using virtual imaging trials (VITs), utilizing virtual patients or phantoms combined with accurate computer models of imaging devices.

PURPOSE

In this work, we develop realistic, physiologically-informed models for coronary plaques for application in cardiac imaging VITs.

METHODS

Histology images of plaques at micron-level resolution were used to train a deep convolutional generative adversarial network (DC-GAN) to create a library of anatomically variable plaque models with clinical anatomical realism. The stability of each plaque was evaluated by finite element analysis (FEA) in which plaque components and vessels were meshed as volumes, modeled as specialized tissues, and subjected to the range of normal coronary blood pressures. To demonstrate the utility of the plaque models, we combined them with the whole-body XCAT computational phantom to perform initial simulations comparing standard energy-integrating detector (EID) CT with photon-counting detector (PCD) CT.

RESULTS

Our results show the network is capable of generating realistic, anatomically variable plaques. Our simulation results provide an initial demonstration of the utility of the generated plaque models as targets to compare different imaging devices.

CONCLUSIONS

Vast, realistic, and variable CAD pathologies can be generated to incorporate into computational phantoms for VITs. There they can serve as a known truth from which to optimize and evaluate cardiac imaging technologies quantitatively.

A stochastic lipid structured model for macrophage dynamics in atherosclerotic plaques

We propose to model certain aspects of the dynamics of a macrophage that moves randomly in a one dimensional space in arterial wall tissue and grows by accumulating localized lipid particles, thus reducing its motility. This phenomenon has been observed in the context of atherosclerotic plaque formation. For this purpose, we use a system of stochastic differential equations satisfied by the position and diffusion coefficient of a Brownian particle whose diffusion coefficient is modified at each visit to the origin and with a dumping coefficient. The novelty of the model, with respect to Bénichou et al. (Phys Rev E 85(2):021137, 2012), Meunier et al. (Acta Appl Math 161:107-126, 2019), is to include offloading of lipids through the dumping term. We find explicit necessary and sufficient conditions for macrophage trapping in the locally enriched region.

Haemodynamic Effects on the Development and Stability of Atherosclerotic Plaques in Arterial Blood Vessel

Studying the formation and stability of atherosclerotic plaques in the hemodynamic field is essential for understanding the growth mechanism and preventive treatment of atherosclerotic plaques. In this paper, based on a multiplayer porous wall model, we established a two-way fluid-solid interaction with time-varying inlet flow. The lipid-rich necrotic core (LRNC) and stress in atherosclerotic plaque were described for analyzing the stability of atherosclerotic plaques during the plaque growth by solving advection-diffusion-reaction equations with finite-element method. It was found that LRNC appeared when the lipid levels of apoptotic materials (such as macrophages, foam cells) in the plaque reached a specified lower concentration, and increased with the plaque growth. LRNC was positively correlated with the blood pressure and was negatively correlated with the blood flow velocity. The maximum stress was mainly located at the necrotic core and gradually moved toward the left shoulder of the plaque with the plaque growth, which increases the plaque instability and the risk of the plaque shedding. The computational model may contribute to understanding the mechanisms of early atherosclerotic plaque growth and the risk of instability in the plaque growth.

Investigating key cell types and molecules dynamics in PyMT mice model of breast cancer through a mathematical model

The most common kind of cancer among women is breast cancer. Understanding the tumor microenvironment and the interactions between individual cells and cytokines assists us in arriving at more effective treatments. Here, we develop a data-driven mathematical model to investigate the dynamics of key cell types and cytokines involved in breast cancer development. We use time-course gene expression profiles of a mouse model to estimate the relative abundance of cells and cytokines. We then employ a least-squares optimization method to evaluate the model's parameters based on the mice data. The resulting dynamics of the cells and cytokines obtained from the optimal set of parameters exhibit a decent agreement between the data and predictions. We perform a sensitivity analysis to identify the crucial parameters of the model and then perform a local bifurcation on them. The results reveal a strong connection between adipocytes, IL6, and the cancer population, suggesting them as potential targets for therapies.

Modeling cholesterol metabolism and atherosclerosis

Atherosclerotic cardiovascular disease (ASCVD) is the leading cause of morbidity and mortality among Western populations. Many risk factors have been identified for ASCVD; however, elevated low-density lipoprotein cholesterol (LDL-C) remains the gold standard. Cholesterol metabolism at the cellular and whole-body level is maintained by an array of interacting components. These regulatory mechanisms have complex behavior. Likewise, the mechanisms which underpin atherogenesis are nontrivial and multifaceted. To help overcome the challenge of investigating these processes mathematical modeling, which is a core constituent of the systems biology paradigm has played a pivotal role in deciphering their dynamics. In so doing models have revealed new insights about the key drivers of ASCVD. The aim of this review is fourfold; to provide an overview of cholesterol metabolism and atherosclerosis, to briefly introduce mathematical approaches used in this field, to critically discuss models of cholesterol metabolism and atherosclerosis, and to highlight areas where mathematical modeling could help to investigate in the future. This article is categorized under: Cardiovascular Diseases > Computational Models.

Mathematical Modeling of Severe Acute Respiratory Syndrome Coronavirus 2 Infection Network with Cytokine Storm, Oxidative Stress, Thrombosis, Insulin Resistance, and Nitric Oxide Pathways

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is a systemic disease affecting not only the lungs but also multiple organ systems. Clinical studies implicate that SARS-CoV-2 infection causes imbalance of cellular homeostasis and immune response that trigger cytokine storm, oxidative stress, thrombosis, and insulin resistance. Mathematical modeling can offer in-depth understanding of the SARS-CoV-2 infection and illuminate how subcellular mechanisms and feedback loops underpin disease progression and multiorgan failure. We report here a mathematical model of SARS-CoV-2 infection pathway network with cytokine storm, oxidative stress, thrombosis, insulin resistance, and nitric oxide (NO) pathways. The biochemical systems theory model shows autocrine loops with positive feedback enabling excessive immune response, cytokines, transcription factors, and interferons, which can imbalance homeostasis of the system. The simulations suggest that changes in immune response led to uncontrolled release of cytokines and chemokines, including interleukin (IL)-1β, IL-6, and tumor necrosis factor α (TNFα), and affect insulin, coagulation, and NO signaling pathways. Increased production of NETs (neutrophil extracellular traps), thrombin, PAI-1 (plasminogen activator inhibitor-1), and other procoagulant factors led to thrombosis. By analyzing complex biochemical reactions, this model forecasts the key intermediates, potential biomarkers, and risk factors at different stages of COVID-19. These insights can be useful for drug discovery and development, as well as precision treatment of multiorgan implications of COVID-19 as seen in systems medicine.

A Turing mechanism in order to explain the patchy nature of Crohn's disease

Crohn's disease is an inflammatory bowel disease (IBD) that is not well understood. In particular, unlike other IBDs, the inflamed parts of the intestine compromise deep layers of the tissue and are not continuous but separated and distributed through the whole gastrointestinal tract, displaying a patchy inflammatory pattern. In the present paper, we introduce a toy-model which might explain the appearance of such patterns. We consider a reaction-diffusion system involving bacteria and phagocyte and prove that, under certain conditions, this system might reproduce an activator-inhibitor dynamic leading to the occurrence of Turing-type instabilities. In other words, we prove the existence of stable stationary solutions that are spatially periodic and do not vanish in time. We also propose a set of parameters for which the system exhibits such phenomena and compare it with realistic parameters found in the literature. This is the first time, as far as we know, that a Turing pattern is investigated in inflammatory models.

Controlling the thickness of the atherosclerotic plaque by statin medication

Atherosclerosis, a chronic inflammatory disorder of the arterial wall, is a complex process whose dynamics are affected by multiple factors. The disease control consists of restraining it by administering statins. Slowing down or halting the plaque growth depends on the patient age at which the statin treatment begins and on the thickness of the intima-media (IMT) at that time. In this paper, we propose a mathematical model to estimate the sets of atherosclerosis states, from which the use of statins can restrain the disease. Our model is control-theoretic, and the estimated sets are the viability kernels, in the parlance of viability theory. To our best knowledge, this way of modelling the atherosclerosis progression is original. We compute two viability kernels, each for a different statin-treatment dose. Each kernel is composed of the vector [age, IMT] from which the disease can be restrained. By extension, the disease can’t be restrained from the kernel complements, this being mainly because of the disease and patient-age advancement. The kernels visualise tradeoffs between early and late treatments, which helps the clinician to decide when to start the statin treatment and which statin dose may be sufficient.

An integrated approach to simulating the vulnerable atherosclerotic plaque

Analyses of individual atherosclerotic plaques are mostly descriptive, relying, for example, on histological classification by spectral analysis of ultrasound waves or staining and observing particular cellular components. Such passive methods have proved useful for characterizing the structure and vulnerability of plaques but have little quantitative predictive power. Our aim is to introduce and discuss a computational framework to provide insight to clinicians and help them visualize internal plaque dynamics. We use partial differential equations (PDEs) with macrophages, necrotic cells, oxidized lipids, oxygen concentration, and platelet-derived growth factor (PDGF) as primary variables coupled to a biomechanical model to describe vessel growth. The model is deterministic, providing mechanical, morphological, and histological characteristics of an atherosclerotic vessel at any desired future time point. We use our model to create computer-generated animations of a plaque evolution that are in qualitative agreement with published serial ultrasound images and hypothesize possible atherogenic mechanisms. A systems biology model consisting of five differential equations is able to capture the morphology of necrotic cores residing within vulnerable atherosclerotic plaque. In the context of the model, the distribution of oxidized low-density lipoprotein (Ox-LDL) particles, endothelial inflammation, plaque oxygenation (via the presence of vasa vasora), and intimal oxygenation are four important factors that drive changes in core morphology.NEW & NOTEWORTHY In this article, we propose a quantitative framework to describe the evolution of atherosclerotic plaque. We use partial differential equations (PDEs) with macrophages, necrotic cells, oxidized lipids, oxygen concentration, and PDGF as primary variables coupled to a biomechanical model to describe vessel growth. A feature of our method is that it outputs color-coded vessel sections corresponding to regions of the plaque that are necrotic and fibrous, qualitatively similar to images generated by enhanced intravascular ultrasound.

Simple model of atherosclerosis in cylindrical arteries: impact of anisotropic growth on Glagov remodeling

In 1987, Seymour Glagov observed that arteries went through a two-stage remodeling process as a result of plaque growth: first, a compensatory phase where the lumen area remains approximately constant and second, an encroachment phase where the lumen area decreases over time. In this paper, we investigate the effect of growth anisotropy on Glagov remodeling in five different cases: pure radial, pure circumferential, pure axial, isotropic and general anisotropic growth where the elements of the growth tensor are chosen to minimize the total energy. We suggest that the nature of anisotropy is inclined towards the growth direction that requires the least amount of energy. Our framework is the theory of morphoelasticity on an axisymmetric arterial domain. For each case, we explore their specific effect on the Glagov curves. For the latter two cases, we also provide the changes in collagen fiber orientation and length in the intima, media and adventitia. In addition, we compare the total energy produced by growth in radial, circumferential and axial direction and deduce that using a radially dominant anisotropic growth leads to lower strain energy than isotropic growth.

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.

A Spatially Resolved and Quantitative Model of Early Atherosclerosis

Atherosclerosis is a major burden for all societies, and there is a great need for a deeper understanding of involved key inflammatory, immunological and biomechanical processes. A decisive step for the prevention and medical treatment of atherosclerosis is to predict what conditions determine whether early atherosclerotic plaques continue to grow, stagnate or become regressive. The driving biological and mechanobiological mechanisms that determine the stability of plaques are yet not fully understood. We develop a spatially resolved and quantitative mathematical model of key contributors of early atherosclerosis. The stability of atherosclerotic model plaques is assessed to identify and classify progression-prone and progression-resistant atherosclerotic regions based on measurable or computable in vivo inputs, such as blood cholesterol concentrations and wall shear stresses. The model combines Darcy's law for the transmural flow through vessels walls, the Kedem-Katchalsky equations for endothelial fluxes of lipoproteins, a quantitative model of early plaque formation from a recent publication and a novel submodel for macrophage recruitment. The parameterization and analysis of the model suggest that the advective flux of lipoproteins through the endothelium is decisive, while the influence of the advective transport within the artery wall is negligible. Further, regions in arteries with an approximate wall shear stress exposure below 20% of the average exposure and their surroundings are potential regions where progression-prone atherosclerotic plaques develop.

Identifying important parameters in the inflammatory process with a mathematical model of immune cell influx and macrophage polarization

In an inflammatory setting, macrophages can be polarized to an inflammatory M1 phenotype or to an anti-inflammatory M2 phenotype, as well as existing on a spectrum between these two extremes. Dysfunction of this phenotypic switch can result in a population imbalance that leads to chronic wounds or disease due to unresolved inflammation. Therapeutic interventions that target macrophages have therefore been proposed and implemented in diseases that feature chronic inflammation such as diabetes mellitus and atherosclerosis. We have developed a model for the sequential influx of immune cells in the peritoneal cavity in response to a bacterial stimulus that includes macrophage polarization, with the simplifying assumption that macrophages can be classified as M1 or M2. With this model, we were able to reproduce the expected timing of sequential influx of immune cells and mediators in a general inflammatory setting. We then fit this model to in vivo experimental data obtained from a mouse peritonitis model of inflammation, which is widely used to evaluate endogenous processes in response to an inflammatory stimulus. Model robustness is explored with local structural and practical identifiability of the proposed model a posteriori. Additionally, we perform sensitivity analysis that identifies the population of apoptotic neutrophils as a key driver of the inflammatory process. Finally, we simulate a selection of proposed therapies including points of intervention in the case of delayed neutrophil apoptosis, which our model predicts will result in a sustained inflammatory response. Our model can therefore provide hypothesis testing for therapeutic interventions that target macrophage phenotype and predict outcomes to be validated by subsequent experimentation.

Using experimental data and mathematical analysis, we develop a model for the inflammatory response that includes macrophage polarization between M1 and M2 phenotypes. Dysfunction of this phenotypic switch can disrupt the timely influx and egress of immune cells during the healing process and lead to chronic wounds or disease. The modulation of macrophage population has been suggested as a strategy to dampen inflammation in diseases that feature chronic inflammation, such as diabetes and atherosclerosis. It is therefore important that we learn more about which components of the system drive the population level switch in phenotype. Our model is able to reproduce the expected timing of sequential influx of neutrophils and macrophages in response to an inflammatory stimulus. Model parameters were estimated with weighted least squares fitting to in vivo experimental data from a mouse model of peritonitis while considering identifiability of parameter sets. We perform sensitivity analysis that identifies primary drivers of the system, and predict the effects of variations in these key parameters on immune cell populations.

A lipid-structured model for macrophage populations in atherosclerotic plaques

Atherosclerosis is a chronic inflammatory disease driven by the accumulation of pro-inflammatory, lipid-loaded macrophages at sites inside artery walls. These accumulations lead to the development of atherosclerotic plaques. The rupture of plaques that contain lipid-rich necrotic cores can trigger heart attacks and strokes via occlusion of blood vessels. We construct and analyse a system of partial integro-differential equations that model lipid accumulation by macrophages, the generation of apoptotic cells and the formation of the necrotic core. The model accounts for the following cell behaviours: monocyte recruitment into the plaque and differentiation into macrophages; macrophage ingestion of low density lipoproteins (LDL) and of apoptotic cells and necrotic material; lipid offloading to high density lipoproteins (HDL); macrophage emigration; and apoptosis of macrophages and necrosis of apoptotic cells. With this model, we study how changes in parameters predict the characteristic features of plaque pathology. In particular, we find the qualitative form of lipid distribution across the macrophage population and show that high lipid loads can occur in the absence of LDL ingestion. We also demonstrate the importance of macrophage emigration in mitigating and resolving inflammation and plaque lipid accumulation.

A Control-Theoretic Model of Atherosclerosis

We propose a control-theoretic aggregate model of the progression of atherosclerosis plaque, a chronic inflammatory disease of the arterial wall, to study the basic features of this disease. In the model, we exploit the role of inflammation in the disease progression, and use statins-drugs commonly recommended in atherosclerosis-to control this progression. We use a logistic function to allow for constrained growth of plaque. In the model, both the patient's age and overall health impact the plaque growth and its sensitivity to statins. The model parameters are estimated using original data, or calibrated using published research as well as our own clinical and laboratory studies. We contend that our model helps to gauge the statins' impact on a patient's plaque thickness, hence the disease's progression and cardiovascular risk, without requiring artery scans.

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.

Nonlinear dynamics of early atherosclerotic plaque formation may determine the efficacy of high density lipoproteins (HDL) in plaque regression

We use a computational model to explore the effect of foam cell accumulation on plaque regression following an increase in high density lipoprotein (HDL) influx into the plaque. Atherosclerotic plaque formation is the outcome of cellular and cytokine responses to low density lipoproteins (LDL) that penetrate the artery wall following an injury to the endothelium and become modified. We modelled the cells and cytokines that are most important in plaque formation using partial differential equations. The model includes monocytes and macrophages, foam cells, macrophage chemoattractants, endothelium-stimulating cytokines, modified low density lipoproteins (mod LDL) and HDL. We included interactions both at the endothelium surface and inside the artery wall. The model predicts that when HDL influx into a well-established plaque with large numbers of foam cells is increased, the plaque may not regress but may continue to grow at a slower rate. If HDL influx is increased when a model plaque is recently established and has fewer foam cells, then the plaque does regress. If modLDL influx into the plaque is lowered at the same time that HDL influx increased or the capacity of the HDL to remove cholesterol from foam cells is increased, then the plaque is more likely to regress. The predictions of the model are in qualitative agreement with experimental studies in mice and rabbits. The results suggest that the intrinsic dynamics of reverse cholesterol transport by HDL are important in determining the success of HDL raising in promoting plaque regression.

A multiphysics approach for modeling early atherosclerosis

This work is devoted to the development of a mathematical model of the early stages of atherosclerosis incorporating processes of all time scales of the disease and to show their interactions. The cardiovascular mechanics is modeled by a fluid-structure interaction approach coupling a non-Newtonian fluid to a hyperelastic solid undergoing anisotropic growth and a change of its constitutive equation. Additionally, the transport of low-density lipoproteins and its penetration through the endothelium is considered by a coupled set of advection-diffusion-reaction equations. Thereby, the permeability of the endothelium is wall-shear stress modulated resulting in a locally varying accumulation of foam cells triggering a novel growth and remodeling formulation. The model is calibrated and applied to an murine-specific case study, and a qualitative validation of the computational results is performed. The model is utilized to further investigate the influence of the pulsatile blood flow and the compliance of the artery wall to the atherosclerotic process. The computational results imply that the pulsatile blood flow is crucial, whereas the compliance of the aorta has only a minor influence on atherosclerosis. Further, it is shown that the novel model is capable to produce a narrowing of the vessel lumen inducing an adaption of the endothelial permeability pattern.

Anti-atherosclerotic effect of the angiotensin 1-7 mimetic AVE0991 is mediated by inhibition of perivascular and plaque inflammation in early atherosclerosis

BACKGROUND AND PURPOSE

Inflammation plays a key role in atherosclerosis. The protective role of angiotensin 1-7 (Ang-(1-7)) in vascular pathologies suggested the therapeutic use of low MW, non-peptide Ang-(1-7) mimetics, such as AVE0991. The mechanisms underlying the vaso-protective effects of AVE0991, a Mas receptor agonist, remain to be explored.

EXPERIMENTAL APPROACH

We investigated the effects of AVE0991 on the spontaneous atherosclerosis in apolipoprotein E (ApoE)-/- mice, in the context of vascular inflammation and plaque stability.

KEY RESULTS

AVE0991 has significant anti-atherosclerotic properties in ApoE-/- mice and increases plaque stability, by reducing plaque macrophage content, without effects on collagen. Using the descending aorta of chow-fed ApoE-/- mice, before significant atherosclerotic plaque develops, we gained insight to early events in atherosclerosis. Interestingly, perivascular adipose tissue (PVAT) and adventitial infiltration with macrophages and T-cells precedes atherosclerotic plaque or the impairment of endothelium-dependent NO bioavailability (a measure of endothelial function). AVE0991 inhibited perivascular inflammation, by reducing chemokine expression in PVAT and through direct actions on monocytes/macrophages inhibiting their activation, characterized by production of IL-1β, TNF-α, CCL2 and CXCL10, and differentiation to M1 phenotype. Pretreatment with AVE0991 inhibited migration of THP-1 monocytes towards supernatants of activated adipocytes (SW872). Mas receptors were expressed in PVAT and in THP-1 cells in vitro, and the anti-inflammatory effects of AVE0991 were partly Mas dependent.

CONCLUSIONS AND IMPLICATIONS

The selective Mas receptor agonist AVE0991 exhibited anti-atherosclerotic and anti-inflammatory actions, affecting monocyte/macrophage differentiation and recruitment to the perivascular space during early stages of atherosclerosis in ApoE-/- mice.

LINKED ARTICLES

This article is part of a themed section on Targeting Inflammation to Reduce Cardiovascular Disease Risk. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v174.22/issuetoc and http://onlinelibrary.wiley.com/doi/10.1111/bcp.v82.4/issuetoc.

Ox-LDL Upregulates IL-6 Expression by Enhancing NF-κB in an IGF2-Dependent Manner in THP-1 Macrophages

Interleukin 6 (IL-6) is a pro-inflammatory cytokine that is well established as a vital factor in determining the risk of coronary heart disease and pathogenesis of atherosclerosis. Moreover, accumulating evidences have shown that oxidized low-density lipoprotein (ox-LDL) can promote IL-6 expression in macrophages. Nevertheless, the underlying mechanism of how ox-LDL upregulates IL-6 expression remains largely unexplained. We found that the expression of insulin-like growth factor 2 (IGF2), nuclear factor kappa B (NF-κB), and IL-6 was upregulated at both the messenger RNA (mRNA) and protein levels in a dose-dependent manner when treated with 0, 25, 50, or 100 μg/mL of ox-LDL for 48 h in THP-1 macrophages. Moreover, overexpression of IGF2 significantly upregulated NF-κB and IL-6 expressions in THP-1 macrophages. However, the upregulation of NF-κB and IL-6 expressions induced by ox-LDL were significantly abolished by IGF2 small interfering RNA (siRNA) in THP-1 macrophages. Further studies indicated the upregulation of IL-6 induced by ox-LDL could be abolished when treated with NF-κB siRNA in THP-1 macrophages. Ox-LDL might upregulate IL-6 in the cell and its secretion via enhancing NF-κB in an IGF2-dependent manner in THP-1 macrophages.