Xi Zhuo Jiang
Kai H. Luo
Yiannis Ventikos- 1School of Mechanical Engineering and Automation, Northeastern University, Shenyang, China
- 2Department of Mechanical Engineering, University College London, London, United Kingdom
Endothelial glycocalyx (EG) is a forest-like structure, covering the lumen side of blood vessel walls. EG is exposed to the mechanical forces of blood flow, mainly shear, and closely associated with vascular regulation, health, diseases, and therapies. One hallmark function of the EG is mechanotransduction, which means the EG senses the mechanical signals from the blood flow and then transmits the signals into the cells. Using numerical modelling methods or in silico experiments to investigate EG-related topics has gained increasing momentum in recent years, thanks to tremendous progress in supercomputing. Numerical modelling and simulation allows certain very specific or even extreme conditions to be fulfilled, which provides new insights and complements experimental observations. This mini review examines the application of numerical methods in EG-related studies, focusing on how computer simulation contributes to the understanding of EG as a mechanotransducer. The numerical methods covered in this review include macroscopic (i.e., continuum-based), mesoscopic [e.g., lattice Boltzmann method (LBM) and dissipative particle dynamics (DPD)] and microscopic [e.g., molecular dynamics (MD) and Monte Carlo (MC) methods]. Accounting for the emerging trends in artificial intelligence and the advent of exascale computing, the future of numerical simulation for EG-related problems is also contemplated.
Introduction
The inner surface of blood vessel walls is covered by a layer of dendritic structures termed endothelial glycocalyx (EG). A typical glycocalyx unit consists of a few glycosaminoglycan (GAG) chains and a core protein which anchors the GAGs to cell membrane. The EG is exposed to the blood flow and is the first barrier in direct contact with blood. Such a unique location allows EG to coordinate microvascular mass transport (Kang et al., 2021), regulate cell adhesion (Robert et al., 2006) and participate in mechanotransduction (Tarbell and Pahakis, 2006). Mechanotransduction describes a process of bio-structures perceiving and transferring their ambient mechanical signals. Specifically for EG, mechanotransduction means the EG senses the mechanical stimuli from blood flow and transmits such signals to the cytoplasm. The response of EG to the haemodynamic forces is critical to maintaining vascular health and function (Goligorsky and Sun, 2020; Weinbaum et al., 2021). Malfunction of the EG in sensing and transmitting mechanical signals can cause severe vascular diseases such as atherosclerosis, hypertension, stroke, and sepsis (Tarbell and Cancel, 2016; Goligorsky and Sun, 2020), just to name a few.
Using wet-lab experiments, especially enzymatic degradation, is an effective way to explore the role of EG and its components in regulating mechanotransduction. Experiments have probed the glycocalyx as a mechanotrasnducer, revealed the potential manner in which the shear (mechanical stimuli) is felt and uncovered the activities involved in the message delivery after sensing the mechanical signal. Experimental evidences supporting the relationship between EG and the mechanotransduction are reviewed in Tarbell and Pahakis (2006).
Complementary to wet-lab experimental methods, in silico experiments, also named numerical simulations or computational modelling, approach the problems of EG mechanotransduction from a different angle. As the name itself indicates, in silico experiments are implemented on computers, which allow precise control of EG properties in order to delineate how EG and its components play their roles in mechanotransduction. Initially, numerical studies of EG adopted simple models only with a limited number of glycocalyx features, and the aim was mainly to reproduce wet-lab results by simulations. The recent progresses in advanced numerical methods and supercomputers empower in silico experiments to tackle far more challenging EG-related problems with high spatial and temporal resolutions. Breakthroughs in crystallographic structure determination lend numerical simulations the possibility to offer accurate results with a resolution at the atomic/molecular level.
The aim of this mini review is to critically summarise recent developments in numerical methods and their applications in the study of the EG mechanotransduction, demonstrating knowledge generation through in silico experiments. Principles of the relevant numerical methods are individually described briefly, together with their contributions to the understanding of the EG functionality as a mechanotransducer. The future of numerical studies in EG functionality is also considered, accounting for trends in artificial intelligence, big data and the advent of exascale computing.
Numerical Methods
Numerical methods can be classified into macroscopic, mesoscopic, and microscopic approaches in terms of length and time scales resolved, as illustrated in Figure 1. Macroscopic simulations have the longest history in solving EG-related problems and are still a popular and powerful tool. Macroscopic simulations deal with problems at cell- or tissue-scales, like endothelial cell deformation, force distribution on glycocalyx ultrastructures and flow shear stress over the EG layer (Tarbell and Shi, 2013; Dabagh et al., 2014). Mesoscopic methods can further explore detailed fluid-structure interactions (O’connor and Revell, 2019) between blood flow and endothelium surface structures but at slightly higher computational cost. Microscopic in silico experiments have been successfully applied to capture the atomic/molecular interactions in a single proteoglycan unit (Jiang et al., 2017). Such microscopic simulations, however, confined to a computation domain measured in nanometers and a physical time in nanoseconds, due to high computational cost.
Figure 1. Numerical methods used in the investigation of the endothelial glycocalyx (EG) mechanotransduction. Computational fluid dynamics (CFD) and the finite element method (FEM) are two established macroscopic methods; lattice Boltzmann method (LBM) and dissipative particle dynamics (DPD) are typical mesoscopic methods; molecular dynamics (MD) and Monte Carlo (MC) are deterministic and stochastic microscopic methods, respectively. Sources: MC is from Xu et al. (2016), MD from Jiang et al. (2017), LBM from O’connor and Revell (2019), DPD from Jiang et al. (2021), CFD from Lee et al. (2020), and FEM from Dabagh et al. (2014). Reuse of figures are permitted by the Cambridge University Press, Royal Society of Chemistry and the Royal Society (United Kingdom).
In this section, the basic principles of these methods are described, and their applications are demonstrated. It should be noted that the main purpose of this mini review is to show advances in using numerical methods to understand the EG mechanotransduction while the algorithmic details of individual simulation methods are described in the relevant references.
Macroscopic Simulations
In macroscopic methods, the EG is generally treated as a continuum layer with a limited number of bio-features, and the simulation outcomes are the macroscopic properties of the endothelial dendritic structures. Computational fluid dynamics (CFD; Tu et al., 2019) and the finite element method (FEM; Zienkiewicz et al., 2005) are two commonly used techniques in macroscopic simulations, which originated from the fluid mechanics and solid mechanics communities, respectively. Generally, continuum-based conservation laws and macroscopic physical models for bio-features are solved to obtain flow profile over and/or through the EG layer (Wang, 2007; Lee et al., 2016, 2020), deformation/force of the anchoring cells or membranes (Tarbell and Shi, 2013; Dabagh et al., 2014), order of magnitude of force transmitted via the glycocalyx (Saez et al., 2019).
In numerical modelling, simplifications of the problems under concern are inevitable due to the complexity of EG structures. With a focus on different aspects of the EG features, mechanotransduction can be approached from different perspectives. For example, the endothelial surface is covered by a layer of interweaving, dynamic and brush-like GAG chains. Such a layer can be simplified as an array of side chains or a layer of porous structures. By assuming the EG layer as an organised array of elastic chains, with individual chains having one end fixed to the cell membrane, Secomb et al. (2001) were able to identify the role of EG in transmitting shear stress to the endothelial cell surface and conclude that the EG layer would mediate rapid flow fluctuations and attenuate shear stress. By treating the EG layer as a porous medium, Tarbell and Shi reported that the solid stress of the endothelial cell surface is one to two orders of magnitude higher than the surface fluid stress in the presence of EG layer, suggesting that the interstitial flow shear stress can be sensed by the EG layer (Tarbell and Shi, 2013). As such, the role of EG is to strengthen the mechanical signals. Combining these results, a possible pathway for force transmission via the EG layer is proposed, i.e., flow shear stress → EG → cell membrane, as discussed in Tarbell and Pahakis (2006).
Simplification based on experimental observation is imperative for achieving numerical results of acceptable accuracy in reasonable timeframes. Take geometric periodicity as an example: Researchers (Lee et al., 2020) have set up different periodic geometries for core proteins and GAGs and analysed the pressure gradient (an indicator used to quantify mechanotransduction in the research) differences among different periodic geometries. Results suggest that the periodicity type determines the pressure gradient within the EG layer: hexagonal periodicity of the core proteins leads to greater mechanotransduction, and rectangular periodicity of aggregated GAGs leads to the greatest conversion of pressure gradients to wall tractions. Indeed, a hexagonal periodicity is favoured by previous studies (Weinbaum et al., 2003; Zeng et al., 2013), and conclusion based on the hexagonal periodicity is sensible.
For macroscopic methods, a “complete” model with all the features of the system is definitely beneficial for accurate results from numerical simulations. On the other hand, the computational cost would increase as more details are involved in computational models. Such a complete model is intimately related with access to accurate structural information which in turn relies on the availability (and possibly invention) of cutting-edge microscopic observation equipment.
Mesoscopic Simulations
The lattice Boltzmann method (LBM) and the dissipative particle dynamics (DPD) are two widely used mesoscopic methodologies in mechnotransduction studies. Both methods are based on the kinetic theory of particle dynamics.
In the framework of LBM, the fluid is modelled as a set of fictitious particles, and such particles interact with each other on a regular discrete lattice mesh through a two-stage “collision and streaming” process. A systematic review of the LBM in flow studies can be found in Chen and Doolen (1998) and Li et al. (2016). In DPD models, the system of interest is simplified to a set of interacting beads, each bead representing a cluster of constitutive molecules or atoms. The governing equation in the DPD method is Newton’s Law of Motion, i.e., F = ma, where F is the force on the bead, m is the mass of the bead, and a is the acceleration of motion of the bead. The force term comprises three elements: a conservative force, a random force and a dissipative force (Liu et al., 2014). A detailed review of DPD can be found in Espanol and Warren (2017). The particle-based nature of both mesoscale methods facilitates the construction of complex structures at the fluid-structure interface. Thus, recent mesoscale simulation studies are mostly about interpreting the fluid-structure interactions between the blood flow and the complex endothelium surface structures (Yin and Zhang, 2013).
At the mesoscales, more features of the glycocalyx can be included, which enables detailed comparison with experiments. By coupling the LBM with immersed boundary methods, alongside structural solvers, researchers revealed a broad range of behaviours for slender structures, from a single flap in a periodic array, to a small finite array of flaps, and finally to a large finite array (Yin and Zhang, 2013; O’connor et al., 2016; O’connor and Revell, 2019). The findings suggest that the flow instability over the slender array depends on the natural frequencies of the flow and the surface structures (O’connor and Revell, 2019). As a simple extension, one may expect that the shedding or impairment of EG chains would change flow patterns, cause flow instability and further disturb mechanotransduction. The role of EG in regulating flow field was also demonstrated in a recent DPD study (Jiang et al., 2021).
Beyond the fluid-structure interactions between the blood flow and EG layer, mesoscopic methods have been applied to a wide spectrum of haemodynamic problems, like cell suspension and leukocyte adhesion over endothelium (Yan et al., 2019), the glycocalyx–endothelium–erythrocyte interaction in the microcirculation (Pontrelli et al., 2015; Jiang et al., 2021) and the formation of the cell-free layer (Deng et al., 2012).
Microscopic Simulations
Molecular dynamics (MD) is a typical microscopic simulation method that describes how atoms and molecules move and interact during a period of time. The governing equation of MD is also Newton’s Law of Motion as in the DPD method, but a wide range of interatomic and intermolecular forces can be included and often parameterised as force fields. For the theoretical and numerical details of the MD method, the readers are referred to Frenkel and Smit (2002) and Rapaport (2011). When the structural configurations of biomolecules are known and a set of force field parameters are properly chosen, the trajectories of the constitutive atoms can be calculated by MD and the behaviour of the biomolecular system is then fully characterised. In principle, all biophysical problems could be tackled by MD if the structural information at the atomic scale were known. The bottleneck preventing the broad application of MD is its high computational cost, which limits not only the number of observed individual molecules to be included but also the duration of observations (i.e., simulated physical time). Despite such constraints, MD has been increasingly employed in studying dynamics of EG and surrounding molecules, thanks to its unique ability to gain atomic level insight. In the classic vasoregulation signalling problem – in what manner does glycocalyx transmit force to the cytoskeleton from the flowing blood?, the glycocalyx was regarded as a rigid body and the force was transmitted by the local torque induced by the movement of the glycocalyx core protein. In recent MD simulations, Jiang et al. (2020b) elucidated the dynamics of Syndecan-4 (a typical type of transmembrane protein of the glycocalyx), and concluded that a scissor-like motion of the Syndecan-4 transmembrane dimer, together with its bending, is responsible for the force transmission and the force transmitted into the cytoplasm is of an order 10–100 pN. The MD model and force transmission mode is summarised in Figure 2. The same all-atom model for MD has also been applied to understand how the local lipid membrane around glycocalyx core proteins responds to the flow shear stress (Jiang et al., 2020a). By artificially switching off the charges of the core protein, the contributions of electrostatic and van der Waals interactions to cell membrane deformation induced by flow were quantified for the first time (Jiang et al., 2020a). The capability of isolating specific effects in MD demonstrates an important advantage of in silico experiments in exploring the mechanisms of mechanotransduction.
Figure 2. Force transmission mode of Syndecan-4 revealed by MD simulations. (Left) An MD model for EG. (Middle) Secondary structure of Syndecan-4. (Right) Force transmission mode via the transmembrane part of Syndecan-4 (Jiang et al., 2020b; reuse of figure under CC BY 4.0 permission).
It is worthwhile to mention that fluid-structure interactions between blood flow and EG can also be interrogated by MD. When the current supercomputing platforms are employed, large-scale MD simulations are able to tackle flow problems in the proximity of the EG layer. As reported in recent studies, the blood flow changes the conformation of EG (Cruz-Chu et al., 2014; Pikoula et al., 2018), and the EG spatial configuration in turn modifies the blood flow profile (Jiang et al., 2017, 2018). This is in accordance with the clearly detrimental cycle between the endothelial dysfunction and the degradation of glycocalyx (Zhang et al., 2018).
Together with MD, the Monte Carlo (MC) method, based on the stochastic theory, is also able to capture microscopic features of bio-systems. A recent MC study has demonstrated that the glycocalyx promotes cooperative binding and clustering of adhesion receptors (Xu et al., 2016). The role of glycocalyx in nanocarrier-cell adhesion is also investigated by the MC method (Agrawal and Radhakrishnan, 2007).
Discussion
In the previous sections, the main features of macroscopic, meoscopic and microscopic simulation methods and their applications in the EG studies are reviewed. Each of these methods (i.e., CFD, FEM; LBM, DPD; MD, and MC) is designed to cover a specific range of length and time scales, and plays a different role in revealing the multiscale multiphysics nature of EG-related phenomena. These methods, however, are not mutually exclusive and can be used in combination. For example, MD, DPD, and LBM have all been used to study the fluid-structure interactions between blood flow and EG, which provide different levels of understanding. MC and LBM have been both used to explore the acceptor-receptor binding, revealing the stochastic and deterministic features, respectively. Moreover, the cell membrane deformation phenomena have been investigated by continuum models and MD, which can give interesting comparisons between macroscopic and microscopic properties of the same process. In short, a hierarchy of numerical simulations have been developed and applied to EG studies, which complement experiments and provide time-dependent and 3D information.
In this mini review, we have focused on the biophysical aspects of mechanotransduction, as the vast majority of numerical studies have been devoted to so far. The biochemical aspects of mechanotransduction, on the other hand, have been rarely studied using numerical simulations, due to the added complexity of biochemical reactions, like the generation of NO (Ignarro, 1989; Bartosch et al., 2017). In principle, the hierarchy of numerical methods can all be developed to reproduce the biochemical processes in mechanotransduction. For example, quantum mechanics/molecular mechanics (QM/MM) computations, which fall under microscopic methods (subatomic scale, to be precise) and advances in probing the electron transfer process, can be applied to find reaction mechanisms of enzymes, like the activation of eNOS (endothelial nitric oxide synthase) as reviewed in Bignon et al. (2019). Also, reactive force field MD (Senftle et al., 2016) simulations, which uses QM-trained force fields to mimic bond breaking and formation and have been extensively proved to be effective in capturing reaction pathways, can be applied to work out the chain reactions after the release of Ca2+ initiated by mechanotransduction.
Although mechanotransduction has been investigated using a multitude of methods, in most case, these methods are applied individually or sequentially. A truly multiscale numerical method requires the dynamic coupling of different simulation methods in one simulation (Fedosov and Karniadakis, 2009). For example, a coupling between DPD and MD via a special multi-scale interface server (Wang et al., 2019) could solve the fluid problem and force transmission simultaneously. The challenge, however, is to find a suitable mapping between high- and low-resolution simulations in both spatial and temporal senses. Despite such inherent difficulties, dynamic coupling of simulations is expected to play an ever increasing role in studying multiscale phenomena such as mechanotranduction, as more powerful and cheaper computing platforms become as available.
Finally, we anticipate increasingly coupled use of artificial intelligence and numerical simulations in the coming years. Already, machine learning techniques have been used to empower MD to simulate a molecular system of 100 million atoms with ab initio accuracy (Jia et al., 2020), which is unprecedented. Together with the advent of exascale supercomputers, the integration of artificial intelligence and big data strategies into numerical simulations will revolutionise the way to conduct in silico experiments, leading to better understanding of the detailed mechanisms of mechanotransduction. Another significant development based on numerical simulations is the construction of a digital twin for a bio-system such as an EG layer, which will enable real-time diagnosis of diseases and, better still, instant repair of diseased components.
Author Contributions
XJ drafted the manuscript. KL and YV offered the direction, proposed ideas, and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was in part supported by the UK Engineering and Physical Sciences Research Council under the project “UK Consortium on Mesoscale Engineering Sciences (UKCOMES)” (Grant No. EP/R029598/1). XJ acknowledges the financial support from the Starting Fund of Northeastern University, China.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Agrawal, N. J., and Radhakrishnan, R. (2007). The Role of Glycocalyx in Nanocarrier-Cell Adhesion Investigated Using a Thermodynamic Model and Monte Carlo Simulations. J. Phys. Chem. C Nanomater Interfaces 111, 15848–15856. doi: 10.1021/jp074514x
Bartosch, A. M. W., Mathews, R., and Tarbell, J. M. (2017). Endothelial Glycocalyx-Mediated Nitric Oxide Production in Response to Selective AFM Pulling. Biophys. J. 113, 101–108. doi: 10.1016/j.bpj.2017.05.033
Bignon, E., Rizza, S., Filomeni, G., and Papaleo, E. (2019). Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase. Comput. Struct. Biotechnol. J. 17, 415–429. doi: 10.1016/j.csbj.2019.03.011
Chen, S., and Doolen, G. D. (1998). Lattice Boltzmann Method for Fluid Flows. Annu. Rev. Fluid Mech. 30, 329–364.
Cruz-Chu, E. R., Malafeev, A., Pajarskas, T., Pivkin, I. V., and Koumoutsakos, P. (2014). Structure and response to flow of the glycocalyx layer. Biophys. J. 106, 232–243. doi: 10.1016/j.bpj.2013.09.060
Dabagh, M., Jalali, P., Butler, P. J., and Tarbell, J. M. (2014). Shear-induced force transmission in a multicomponent, multicell model of the endothelium. J. R. Soc. Interface. 11:20140431. doi: 10.1098/rsif.2014.0431
Deng, M., Li, X., Liang, H., Caswell, B., and Karniadakis, G. E. (2012). Simulation and Modelling of Slip Flow Over Surfaces Grafted with Polymer Brushes and Glycocalyx Fibres. United Kingdom: Cambridge University Press.
Espanol, P., and Warren, P. B. (2017). Perspective: dissipative particle dynamics. J. Chem. Phys. 146:150901. doi: 10.1063/1.4979514
Fedosov, D. A., and Karniadakis, G. E. (2009). Triple-decker: Interfacing atomistic–mesoscopic–continuum flow regimes. J. Comput. Phys. 228, 1157–1171. doi: 10.1016/j.jcp.2008.10.024
Goligorsky, M. S., and Sun, D. (2020). Glycocalyx in Endotoxemia and Sepsis. Am. J. Pathol. 190, 791–798. doi: 10.1016/j.ajpath.2019.06.017
Ignarro, L. J. (1989). Biological actions and properties of endothelium-derived nitric oxide formed and released from artery and vein. Circ. Res. 65, 1–21. doi: 10.1161/01.res.65.1.1
Jia, W., Wang, H., Chen, M., Lu, D., Lin, L., Car, R., et al. (2020). “Pushing the Limit of Molecular Dynamics with Ab Initio Accuracy to 100 Million Atoms with Machine Learning” in SC20: international Conference for High Performance Computing, Networking, Storage and Analysis. (United States: ACM). 1–14.
Jiang, X. Z., Feng, M., Ventikos, Y., and Luo, K. H. (2018). Regimes of Flow over Complex Structures of Endothelial Glycocalyx: a Molecular Dynamics Simulation Study. Sci. Rep. 8:5732.
Jiang, X. Z., Goligorsky, M. S., and Luo, K. H. (2021). Cross Talk between Endothelial and Red Blood Cell Glycocalyces via Near-field Flow. Biophys. J. 120, 1–12. doi: 10.1016/j.bpj.2021.06.002
Jiang, X. Z., Gong, H., Luo, K. H., and Ventikos, Y. (2017). Large-scale molecular dynamics simulation of coupled dynamics of flow and glycocalyx: towards understanding atomic events on an endothelial cell surface. J. R. Soc. Interface 14:20170780. doi: 10.1098/rsif.2017.0780
Jiang, X. Z., Guo, L., Luo, K. H., and Ventikos, Y. (2020a). Membrane Deformation of Endothelial Surface Layer Interspersed with Syndecan-4: a Molecular Dynamics Study. Ann. Biomed. Eng. 48, 357–366. doi: 10.1007/s10439-019-02353-7
Jiang, X. Z., Luo, K. H., and Ventikos, Y. (2020b). Principal mode of Syndecan-4 mechanotransduction for the endothelial glycocalyx is a scissor-like dimer motion. Acta Physiol. 228:e13376.
Kang, H., Yang, J., Zhang, W., Lu, J., Ma, X., Sun, A., et al. (2021). Effect of endothelial glycocalyx on water and LDL transport through the rat abdominal aorta. Am. J. Physiol. Heart Circ. Physiol. 320, H1724–H1737.
Lee, T. C., Long, D. S., and Clarke, R. J. (2016). Effect of endothelial glycocalyx layer redistribution upon microvessel poroelastohydrodynamics. J. Fluid Mech. 798, 812–852. doi: 10.1017/jfm.2016.337
Lee, T. C., Suresh, V., and Clarke, R. J. (2020). Influence of Endothelial Glycocalyx Layer Microstructure upon its Role as a Mechanotransducer. United Kingdom: Cambridge University Press.
Li, Q., Luo, K. H., Kang, Q. J., He, Y. L., Chen, Q., and Liu, Q. (2016). Lattice Boltzmann methods for multiphase flow and phase-change heat transfer. Prog. Energy Combust. Sci. 52, 62–105.
Liu, M. B., Liu, G. R., Zhou, L. W., and Chang, J. Z. (2014). Dissipative Particle Dynamics (DPD): an Overview and Recent Developments. Arch. Comput. Methods Eng. 22, 529–556. doi: 10.1007/s11831-014-9124-x
O’connor, J., and Revell, A. (2019). Dynamic interactions of multiple wall-mounted flexible flaps. J. Fluid Mech. 870, 189–216. doi: 10.1017/jfm.2019.266
O’connor, J., Revell, A., Mandal, P., and Day, P. (2016). Application of a lattice Boltzmann-immersed boundary method for fluid-filament dynamics and flow sensing. J. Biomech. 49, 2143–2151. doi: 10.1016/j.jbiomech.2015.11.057
Pikoula, M., Tessier, M. B., Woods, R. J., and Ventikos, Y. (2018). Oligosaccharide model of the vascular endothelial glycocalyx in physiological flow. Microfluid Nanofluidics 22:21.
Pontrelli, G., Halliday, I., Spencer, T. J., Konig, C. S., and Collins, M. W. (2015). Modelling the glycocalyx-endothelium-erythrocyte interaction in the microcirculation: a computational study. Comput. Methods Biomech. Biomed. Eng. 18, 351–361. doi: 10.1080/10255842.2013.799146
Rapaport, D. C. (2011). The Art of Molecular Dynamics Simulation. United Kingdom: Cambridge University Press.
Robert, P., Limozin, L., Benoliel, A.-M., Pierres, A., and Bongrand, P. (2006). “7 - GYCOCALYX REGULATION OF CELL ADHESION” in Principles of Cellular Engineering. ed. M. R. King (Burlington: Academic Press). 143–169. doi: 10.1016/b978-012369392-1/50008-5
Saez, P., Gallo, D., and Morbiducci, U. (2019). Mechanotransmission of haemodynamic forces by the endothelial glycocalyx in a full-scale arterial model. R. Soc. Open Sci. 6:190607. doi: 10.1098/rsos.190607
Secomb, T. W., Hsu, R., and Pries, A. R. (2001). Effect of the endothelial surface layer on transmission of fluid shear stress to endothelial cells. Biorheology 38, 143–150.
Senftle, T. P., Hong, S., Islam, M. M., Kylasa, S. B., Zheng, Y., Shin, Y. K., et al. (2016). The ReaxFF reactive force-field: development, applications and future directions. NPJ Comput. Mater. 2:15011.
Tarbell, J. M., and Cancel, L. M. (2016). The glycocalyx and its significance in human medicine. J. Intern. Med. 280, 97–113. doi: 10.1111/joim.12465
Tarbell, J. M., and Pahakis, M. Y. (2006). Mechanotransduction and the glycocalyx. J. Intern. Med. 259, 339–350. doi: 10.1111/j.1365-2796.2006.01620.x
Tarbell, J. M., and Shi, Z. D. (2013). Effect of the glycocalyx layer on transmission of interstitial flow shear stress to embedded cells. Biomech. Model Mechanobiol. 12, 111–121. doi: 10.1007/s10237-012-0385-8
Tu, J., Yeoh, G.-H., and Liu, C. (2019). Computational Fluid Dynamics (Third Edition), Elsevier Butterworth Heinemann.
Wang, W. (2007). Change in properties of the glycocalyx affects the shear rate and stress distribution on endothelial cells. J. Biomech. Eng. 129, 324–329. doi: 10.1115/1.2720909
Wang, Y., Li, Z., Xu, J., Yang, C., and Karniadakis, G. E. (2019). Concurrent coupling of atomistic simulation and mesoscopic hydrodynamics for flows over soft multi-functional surfaces. Soft Matter 15, 1747–1757. doi: 10.1039/c8sm02170h
Weinbaum, S., Cancel, L. M., Fu, B. M., and Tarbell, J. M. (2021). The Glycocalyx and Its Role in Vascular Physiology and Vascular Related Diseases. Cardiovasc. Eng. Technol. 12, 37–71. doi: 10.1007/s13239-020-00485-9
Weinbaum, S., Zhang, X., Han, Y., Vink, H., and Cowin, S. C. (2003). Mechanotransduction and flow across the endothelial glycocalyx. Proc. Natl. Acad. Sci. U. S. A. 100, 7988–7995. doi: 10.1073/pnas.1332808100
Xu, G. K., Qian, J., and Hu, J. (2016). The glycocalyx promotes cooperative binding and clustering of adhesion receptors. Soft Matter 12, 4572–4583. doi: 10.1039/c5sm03139g
Yan, W., Liu, Y., and Fu, B. (2019). LBM simulations on the influence of endothelial SGL structure on cell adhesion in the micro-vessels. Comput. Math. Appl. 78, 1182–1193. doi: 10.1016/j.camwa.2016.07.005
Yin, X., and Zhang, J. (2013). Modeling the dynamic flow-fiber interaction for microscopic biofluid systems. J. Biomech. 46, 314–318. doi: 10.1016/j.jbiomech.2012.11.001
Zeng, Y., Waters, M., Andrews, A., Honarmandi, P., Ebong, E. E., Rizzo, V., et al. (2013). Fluid shear stress induces the clustering of heparan sulfate via mobility of glypican-1 in lipid rafts. Am. J. Physiol. Heart Circ. Physiol. 305, H811–H820.
Zhang, X., Sun, D., Song, J. W., Zullo, J., Lipphardt, M., Coneh-Gould, L., et al. (2018). Endothelial cell dysfunction and glycocalyx - A vicious circle. Matrix Biol. 7, 421–431. doi: 10.1016/j.matbio.2018.01.026
Keywords: microscale, mesoscale, macroscale, computational simulation, mechanotransduction, endothelial glycocalyx, numerical modelling
Citation: Jiang XZ, Luo KH and Ventikos Y (2021) Understanding the Role of Endothelial Glycocalyx in Mechanotransduction via Computational Simulation: A Mini Review. Front. Cell Dev. Biol. 9:732815. doi: 10.3389/fcell.2021.732815
Received: 29 June 2021; Accepted: 27 July 2021;
Published: 17 August 2021.
Edited by:
Bingmei M. Fu, City College of New York, United StatesReviewed by:
Yang Liu, Hong Kong Polytechnic University, ChinaMiguel Vicente-Manzanares, Consejo Superior de Investigaciones Científicas, Spain
Copyright © 2021 Jiang, Luo and Ventikos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xi Zhuo Jiang, jiangxz@mail.neu.edu.cn; Kai H. Luo, k.luo@ucl.ac.uk; Yiannis Ventikos, y.ventikos@ucl.ac.uk