Abstract:An explicit immersed boundary-lattice Boltzmann method is applied to numerically investigate the dynamics of primary cilium in pulsatile blood flows with two-way fluid–structure interaction considered. To well characterize the effect of cilium basal body on cilium dynamics, the cilium base is modeled as a nonlinear rotational spring attached to the cilium's basal end as proposed by Resnick (Biophys J 109:18–25, <a href="/article/10.1007/s10237-019-01192-8?shared-article-renderer#ref-CR34">2015</a>. <a href="https://doi.org/10.1016/j.bpj.2015.05.031">https://doi.org/10.1016/j.bpj.2015.05.031</a>). After several careful validations, the fluid–cilium interaction system is investigated in detail at various pulsatile flow conditions that are characterized by peak Reynolds numbers (\(Re_{{\rm peak}}\)) and Womersley numbers (\(Wo\)). The periodic flapping of primary cilium observed in our simulations is very similar to the in vivo ciliary oscillation captured by O'Connor et al. (Cilia 2:8, <a href="/article/10.1007/s10237-019-01192-8?shared-article-renderer#ref-CR29">2013</a>. <a href="https://doi.org/10.1186/2046-2530-2-8">https://doi.org/10.1186/2046-2530-2-8</a>). The cilium's dynamics is found to be closely related to the \(Re_{{\rm peak}}\) and \(Wo\). Increase the \(Re_{{\rm peak}}\) or decrease the \(Wo\) bring to an increase in the cilium's flapping amplitude, tip angular speed, basal rotation, and maximum tensile stress. It is also demonstrated that by reducing the \(Re_{{\rm peak}}\) or enhancing the \(Wo\) to a certain level, one can shift the flapping pattern of cilium from its original two-side one to a one-side one, making the stretch only happen on one particular side. During the flapping process, the location of the maximum tensile stress is not always found at the basal region; instead, it is able to propagate from time to time within a certain distance to the base. Due to the obstruction of the primary cilium, the distribution of wall shear stress no longer remains uniform as in the absence of cilia. It oscillates in space with the minimum magnitude which is always found near where the cilium is located. The presence of cilium also reduces the overall level of wall shear stress, especially at the region near the cilium's anchor point.

Primary cilium: a paradigm for integrating mathematical modeling with experiments and numerical simulations in mechanobiology

Microstructure‐based modeling of primary cilia mechanics

The cilium as a force sensor−myth versus reality

MECHANOTRANSDUCTION OF THE CELL AND ITS PRIMARY CILIUM IN THE MICROFLUIDIC CHANNEL

Ciliary Central Apparatus Structure Reveals Mechanisms of Microtubule Patterning

Cell Cycle Regulation Through Primary Cilium: A Long-Forgotten Story

Mechano-Sensing and Shear Stress-Shielding by Endothelial Primary Cilia: Structure, Composition, and Function

Advances in Modeling Cellular Mechanical Perceptions and Responses Via the Membrane-Cytoskeleton-Nucleus Machinery

Emerging mechanistic understanding of cilia function in cellular signalling

Cellular mechanosensing of the biophysical microenvironment: A review of mathematical models of biophysical regulation of cell responses.

Numerical study on the dynamics of primary cilium in pulsatile flows by the immersed boundary-lattice Boltzmann method

Primary cilia as a novel horizon between neuron and environment

Evaluating efficiency and robustness in cilia design

Functional Aspects of Primary Cilium in Signaling, Assembly and Microenvironment in Cancer.

Review and Perspective on Soft Matter Modeling in Cellular Mechanobiology: Cell Contact, Adhesion, Mechanosensing, and Motility

Computational Model for Cell Shape Regulation Through Mechanosensing and Mechanical Feedback

The role of primary cilium in bone tissue mechanotransduction

Effect of curvature and normal forces on motor regulation of cilia

Primary Cilium-Dependent Signaling Mechanisms

The Cilium: Cellular Antenna and Central Processing Unit

Realizing the Physics of Motile Cilia Synchronization with Driven Colloids