Skip Navigation Links
Journal of Applied Nonlinear Dynamics
Miguel A. F. Sanjuan (editor), Albert C.J. Luo (editor)
Miguel A. F. Sanjuan (editor)

Department of Physics, Universidad Rey Juan Carlos, 28933 Mostoles, Madrid, Spain


Albert C.J. Luo (editor)

Department of Mechanical and Industrial Engineering, Southern Illinois University Ed-wardsville, IL 62026-1805, USA

Fax: +1 618 650 2555 Email:

Stability Analysis of Mathematical Modeling of Atherosclerotic Plaque Formation

Journal of Applied Nonlinear Dynamics 9(3) (2020) 361--389 | DOI:10.5890/JAND.2020.09.003

Debasmita Mukherjee, Lakshmi Narayan Guin, Santabrata Chakravarty

Department of Mathematics, Visva-Bharati, Santiniketan, 731235, West Bengal, India

Download Full Text PDF



The present theoretical investigation is dealt with stability of the model system characterizing the formation of atherosclerotic plaques emerging from the interactions between several cellular species in the media of blood stream and its surrounding vascular region. The biochemical processes behind the formation of atherosclerotic plaque involve the interactions between several cellular species like low density lipoprotein, free radicals, chemoattractants, monocytes, macrophages, T-cells, smooth muscle cells, foam cells and collagen. Taking all these complex events into account, an appropriate mathematical model depicting the onset of atherosclerotic plaques in the arterial lumen is constructed through the system of ten nonlinear ordinary differential equations (ODEs) for the concentrations of most pertinent components of the atherosclerotic constraints. Besides conducting an in-depth study of the model including several sub-models for their stability criteria, special emphasis is also paid on a reduced model system having adequate relevance to the present system following Quasi Steady State Approximation (QSSA) theory. Both the local and the global stability together with the bifurcation analysis for the reduced model system are carried out analytically. Results of numerical simulation based on the model parameter values reveal the time-series representations for the concentrations of all the interacting species, the local and the global stability and bifurcations with respect to some parameters of significance for the reduced model system under study. The complex features of several subsystems are also examined through several phase portraits in order to explore the clinical implications of atherosclerotic lesion in the arterial lumen. The model validation is also performed through comparison of the present results with those of previous ones [1].


The authors gratefully acknowledge the financial support by Special Assistance Programme (SAP-III) sponsored by the University Grants Commission (UGC), New Delhi, India [Grants no. F.510/3/DRSIII/2015(SAP-I)].We would like to express thank the anonymous referee and the editor for supportive remarks and ideas.


  1. [1]  Bulelzai, M.A. and Dubbeldam, J.L. (2012), Long time evolution of atherosclerotic plaques, Journal of Theoretical Biology, 297, 1-10.
  2. [2]  Crowther, M.A. (2005), Pathogenesis of Atherosclerosis, ASH Education Program Book, 436-441.
  3. [3]  Malek, A.M., Alper, S.L., and Izumo, S. (1999), Hemodynamic shear stress and its role in atherosclerosis, Journal of the American Medical Association, 282, 2035-2042.
  4. [4]  Li, Z.Y., Howarth, S.P., Tang, T., and Gillard, J.H. (2006), How critical is fibrous cap thickness to carotid plaque stability?, Stroke, 37, 1195-1199.
  5. [5]  Gijsen, F.J., Wentzel, J.J., Thury, A., Mastik, F., Schaar, J.A., Schuurbiers, J.C., Slager, C.J., van der Giessen, W.J., de Feyter, A.F., Pim J van der Steen, P.W., and Serruys (2008), Strain distribution over plaques in human coronary arteries relates to shear stress, American Journal of Physiology-Heart and Circu- latory Physiology, 295, H1608-H1614.
  6. [6]  Hansson, G.K. and Libby, P. (2006), The immune response in atherosclerosis: a double-edged sword, Nature Reviews Immunology, 6, 508-519.
  7. [7]  Newby, A.C. and Zaltsman, A.B. (1999), Fibrous cap formation or destruction-the critical importance of vascular smooth muscle cell proliferation, migration and matrix formation, Cardiovascular Research, 41, 345-360.
  8. [8]  Han, K.H., Hong, K.H., Park, J.H., Ko, J., Kang, D.H., Choi, K.J., Hong, M.K., Park, S.W., and Park, S.J. (2004), C-reactive protein promotes monocyte chemoattractant protein-1-mediated chemotaxis through upregulating CC chemokine receptor 2 expression in human monocytes, Circulation, 109, 2566-2571.
  9. [9]  Cohen, A., Myerscough,M.R., and Thompson, R.S. (2014), Athero-protective effects of High Density Lipoproteins (HDL): An ODE model of the early stages of atherosclerosis, Bulletin of Mathematical Biology, 76, 1117–1142.
  10. [10]  Libby, P., Ridker, P.M., andMaseri, A. (2002), Inflammation and atherosclerosis, Circulation, 105, 1135-1143.
  11. [11]  Koskinas, K.C., Chatzizisis, Y.S., Baker, A.B., Edelman, E.R., Stone, P.H., and Feldman, C.L. (2009), The role of low endothelial shear stress in the conversion of atherosclerotic lesions from stable to unstable plaque, Current Opinion in Cardiology, 24, 580-590.
  12. [12]  McKay, C., McKee, S., Mottram, N., Mulholland, T., Wilson, S., Kennedy, S., and Wadsworth, R. (2005), Towards a model of atherosclerosis, University of Strathclyde, 1-29.
  13. [13]  Ougrinovskaia, A., Thompson, R.S., and Myerscough, M.R. (2010), An ODE model of early stages of atherosclerosis: mechanisms of the inflammatory response, Bulletin of Mathematical Biology, 72, 1534-1561.
  14. [14]  El Khatib, N., Genieys, S., and Volpert, V. (2007), Atherosclerosis initiation modeled as an inflammatory process, Mathematical Modelling of Natural Phenomena, 2, 126-141.
  15. [15]  Fok, P.W. (2012), Mathematical model of intimal thickening in atherosclerosis: vessel stenosis as a free boundary problem, Journal of Theoretical Biology, 314, 23-33.
  16. [16]  Chalmers, A.D., Cohen, A., Bursill, C.A., and Myerscough, M.R. (2015), Bifurcation and dynamics in a mathematical model of early atherosclerosis, Journal of Mathematical Biology, 711451–1480.
  17. [17]  Hao, W. and Friedman, A. (2014), The ldl-hdl profile determines the risk of atherosclerosis: a mathematical model, PloS One, 9, e90497.
  18. [18]  Friedman, A. and Hao, W. (2015), A mathematical model of atherosclerosis with reverse cholesterol transport and associated risk factors, Bulletin of Mathematical Biology, 77, 758-781.
  19. [19]  Wang, H.H. (2001), Analytical models of atherosclerosis, Atherosclerosis, 159, 1-7.
  20. [20]  Zohdi, T., Holzapfel, G., and Berger, S. (2004), A phenomenological model for atherosclerotic plaque growth and rupture, Journal of Theoretical Biology, 227, 437-443.
  21. [21]  Calvez, V., Ebde, A., Meunier, N., and Raoult, A. (2009), Mathematical modelling of the atherosclerotic plaque formation, 28 1-12.
  22. [22]  Di Tomaso, G., Díaz-Zuccarini, V., and Pichardo-Almarza, C. (2011), A multiscale model of atherosclerotic plaque formation at its early stage, IEEE Transactions on Biomedical Engineering, 58, 3460-3463.
  23. [23]  Cilla, M., Peña, E., and Martínez, M.A. (2014), Mathematical modelling of atheroma plaque formation and development in coronary arteries, Journal of The Royal Society Interface, 11, 20130866.
  24. [24]  Cilla, M., Martínez, M., and Peña, E. (2015), Effect of transmural transport properties on atheroma plaque formation and development, Annals of Biomedical Engineering, 43, 1516–1530.
  25. [25]  Anlamlert, W., Lenbury, Y., and Bell, J. (2017), Modeling fibrous cap formation in atherosclerotic plaque development: stability and oscillatory behavior, Advances in Difference Equations, 2017, 195.
  26. [26]  Li, X. and Xu, J.X. (2016), A mathematical prognosis model for pancreatic cancer patients receiving immunotherapy, Journal of Theoretical Biology, 406, 42-51.
  27. [27]  Zhang, S., Ritter, L., and Ibragimov, A. (2013), Foam cell formation in atherosclerosis: Hdl and macrophage reverse cholesterol transport, Discrete Contin Dyn Syst Supplements, 825-835.
  28. [28]  Flach, E.H. and Schnell, S. (2006), Use and abuse of the quasi-steady-state approximation, IEE Proceedings- Systems Biology, 153, 187-191.
  29. [29]  Nise, N.S. (2011), (Control system engineering), John Wiley & Sons, Inc, New York.
  30. [30]  Hale, J. (1969), Ordinary differential equations, (Pure and Applied Mathematics), Wiley-Interscience.
  31. [31]  Perko, L. (2008), Differential Equations and Dynamical Systems, Texts in Applied Mathematics, Springer, New York.
  32. [32]  Agarwal, A., Adams, R., Castellani, G.C., and Shouval, H.Z. (2012), On the precision of quasi steady state assumptions in stochastic dynamics, The Journal of chemical physics, 137, 044105.
  33. [33]  Cobbold, C., Sherratt, J., and Maxwell, S. (2002), Lipoprotein oxidation and its significance for atherosclerosis: a mathematical approach, Bulletin of Mathematical Biology, 64, 65-95.