Atherosclerotic plaques form in artery walls due to a chronic inflammatory response driven by lipid accumulation. A key component of the inflammatory response is the interaction between monocyte-derived macrophages and extracellular lipid. Although concentrations of low-density lipoprotein (LDL) and high-density lipoprotein (HDL) particles in the blood are known to affect plaque progression, their impact on the lipid load of plaque macrophages remains unexplored. In this paper, we develop a lipid-structured mathematical model to investigate the impact of blood LDL/HDL levels on plaque composition, and lipid distribution in plaque macrophages. A reduced subsystem, derived by summing the equations of the full model, describes the dynamics of biophysical quantities relating to plaque composition (e.g. total number of macrophages, total amount of intracellular lipid). We also derive a continuum approximation of the model to facilitate analysis of the macrophage lipid distribution. The results, which include time-dependent numerical solutions and asymptotic analysis of the unique steady state solution, indicate that plaque lipid content is sensitive to the influx of LDL relative to HDL capacity. The macrophage lipid distribution evolves in a wave-like manner towards an equilibrium profile which may be monotone decreasing, quasi-uniform or unimodal, attaining its maximum value at a non-zero lipid level. Our model also reveals that macrophage uptake may be severely impaired by lipid accumulation. We conclude that lipid accumulation in plaque macrophages may serve as a partial explanation for the defective uptake of apoptotic cells (efferocytosis) often reported in atherosclerotic plaques.

The retinal vascular network supplies perfusion to vital visual structures, including retinal ganglion cells responsible for vision. Impairments in retinal blood flow and oxygenation are involved in the progression of many ocular diseases, including glaucoma. In this study, an established theoretical hybrid model of a retinal microvascular network is extended to include the effects of local blood flow regulation on oxygenation. A heterogeneous representation of the arterioles based on confocal microscopy images is combined with a compartmental description of the downstream capillaries and venules. A Green's function method is used to simulate oxygen transport in the arterioles, and a Krogh cylinder model is applied to the capillary and venular compartments. Acute blood flow regulation is simulated in response to changes in pressure, shear stress, and metabolism. Model results predict that both increased intraocular pressure and impairment of blood flow regulation can cause decreased tissue oxygenation, indicating that both mechanisms represent factors that could lead to impaired oxygenation characteristic of ocular disease. Results also indicate that the metabolic response mechanism reduces the fraction of poorly oxygenated tissue but that the pressure- and shear stress-dependent response mechanisms may hinder the vascular response to changes in oxygenation. Importantly, the heterogeneity of the vascular network demonstrates that traditionally reported average values of tissue oxygen levels hide significant localized defects in tissue oxygenation that may be involved in disease processes, including glaucoma. Ultimately, the model framework presented in this study will facilitate future comparisons to sectorial-specific clinical data to better assess the role of impaired blood flow regulation in ocular disease.