Microstructure-based modeling of primary cilia mechanics.

A primary cilium, made of nine microtubule doublets enclosed in a cilium membrane, is a mechanosensing organelle that bends under an external mechanical load and sends an intracellular signal through transmembrane proteins activated by cilium bending. The nine microtubule doublets are the main load-bearing structural component, while the transmembrane proteins on the cilium membrane are the main sensing component. No distinction was made between these two components in all existing models, where the stress calculated from the structural component (nine microtubule doublets) was used to explain the sensing location, which may be totally misleading. For the first time, we developed a microstructure-based primary cilium model by considering these two components separately. First, we refined the analytical solution of bending an orthotropic cylindrical shell for individual microtubule, and obtained excellent agreement between finite element simulations and the theoretical predictions of a microtubule bending as a validation of the structural component in the model. Second, by integrating the cilium membrane with nine microtubule doublets and simulating the tip-anchored optical tweezer experiment on our computational model, we found that the microtubule doublets may twist significantly as the whole cilium bends. Third, besides being cilium-length-dependent, we found the mechanical properties of the cilium are also highly deformation-dependent. More important, we found that the cilium membrane near the base is not under pure in-plane tension or compression as previously thought, but has significant local bending stress. This challenges the traditional model of cilium mechanosensing, indicating that transmembrane proteins may be activated more by membrane curvature than membrane stretching. Finally, we incorporated imaging data of primary cilia into our microstructure-based cilium model, and found that comparing to the ideal model with uniform microtubule length, the imaging-informed model shows the nine microtubule doublets interact more evenly with the cilium membrane, and their contact locations can cause even higher bending curvature in the cilium membrane than near the base.

Hydrodynamics of a multicomponent vesicle under strong confinement.

We numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

Microstructure-Based Modeling of Primary Cilia Mechanics.

A primary cilium, made of nine microtubule doublets enclosed in a cilium membrane, is a mechanosensing organelle that bends under an external mechanical load and sends an intracellular signal through transmembrane proteins activated by cilium bending. The nine microtubule doublets are the main load-bearing structural component, while the transmembrane proteins on the cilium membrane are the main sensing component. No distinction was made between these two components in all existing models, where the stress calculated from the structural component (nine microtubule doublets) was used to explain the sensing location, which may be totally misleading. For the first time, we developed a microstructure-based primary cilium model by considering these two components separately. First, we refined the analytical solution of bending an orthotropic cylindrical shell for individual microtubule, and obtained excellent agreement between finite element simulations and the theoretical predictions of a microtubule bending as a validation of the structural component in the model. Second, by integrating the cilium membrane with nine microtubule doublets, we found that the microtubule doublets may twist significantly as the whole cilium bends. Third, besides being cilium-length-dependent, we found the mechanical properties of the cilium are also highly deformation-dependent. More important, we found that the cilium membrane near the base is not under pure in-plane tension or compression as previously thought, but has significant local bending stress. This challenges the traditional model of cilium mechanosensing, indicating that transmembrane proteins may be activated more by membrane curvature than membrane stretching. Finally, we incorporated imaging data of primary cilia into our microstructure-based cilium model, and found that comparing to the ideal model with uniform microtubule length, the imaging-informed model shows the nine microtubule doublets interact more evenly with the cilium membrane, and their contact locations can cause even higher bending curvature in the cilium membrane than near the base. SIGNIFICANCE: Factors regulating the mechanical response of a primary cilium to fluid flow remain unclear. Modeling the microtubule doublet as a composite of two orthotropic shells and the ciliary axoneme as an elastic shell enclosing nine such microtubule doublets, we found that the length distribution of microtubule doublets (inferred from cryogenic electron tomography images) is the primary determining factor in the bending stiffness of primary cilia, rather than just the ciliary length. This implies ciliary-associated transmembrane proteins may be activated by membrane curvature changes rather than just membrane stretching. These insights challenge the traditional view of ciliary mechanosensation and expands our understanding of the different ways in which cells perceive and respond to mechanical stimuli.

Influence of surface viscosities on the electrodeformation of a prolate viscous drop.

Contaminants and other agents are often present at the interface between two fluids, giving rise to rheological properties such as surface shear and dilatational viscosities. The dynamics of viscous drops with interfacial viscosities has attracted greater interest in recent years, due to the influence of surface rheology on deformation and the surrounding flows. We investigate the effects of shear and dilatational viscosities on the electro-deformation of a viscous drop using the Taylor-Melcher leaky dielectric model. We use a large deformation analysis to derive an ordinary differential equation for the drop shape. Our model elucidates the contributions of each force to the overall deformation of the drop and reveals a rich range of dynamic behaviors that show the effects of surface viscosities and their dependence on rheological and electrical properties of the system. We also examine the physical mechanisms underlying the observed behaviors by analyzing the surface dilatation and surface deformation.

The effects of surface hydration on capillary adhesion under nanoscale confinement.

Nanoscale phenomena such as surface hydration and the molecular layering of liquids under strong nanoscale confinement play a critical role in liquid-mediated surface adhesion that is not accounted for by available models, which assume a uniform liquid density with or without considering surface forces and associated disjoining pressure effects. This work introduces an alternative theoretical description that via the potential of mean force (PMF) considers the strong spatial variation of the liquid number density under nanoscale confinement. This alternative description based on the PMF predicts a dual effect of surface hydration by producing: (i) strong spatial oscillations of the local liquid density and pressure and, more importantly, (ii) a configuration-dependent liquid-solid surface energy under nanoscale confinement. Theoretical analysis and molecular dynamics simulations for the case of an axisymmetric water bridge with nanoscale heights show that the latter hydration effect is critical for the accurate prediction of the surface energy and adhesion forces when a small volume of liquid is nanoscopically confined by two surfaces approaching contact.

The effect of rigid cells on blood viscosity: linking rheology and sickle cell anemia.

Sickle cell anemia (SCA) is a disease that affects red blood cells (RBCs). Healthy RBCs are highly deformable objects that under flow can penetrate blood capillaries smaller than their typical size. In SCA there is an impaired deformability of some cells, which are much stiffer and with a different shape than healthy cells, and thereby affect regular blood flow. It is known that blood from patients with SCA has a higher viscosity than normal blood. However, it is unclear how the rigidity of cells is related to the viscosity of blood, in part because SCA patients are often treated with transfusions of variable amounts of normal RBCs and only a fraction of cells will be stiff. Here, we report systematic experimental measurements of the viscosity of a suspension varying the fraction of rigid particles within a suspension of healthy cells. We also perform systematic numerical simulations of a similar mixed suspension of soft RBCs, rigid particles, and their hydrodynamic interactions. Our results show that there is a rheological signature within blood viscosity to clearly identify the fraction of rigidified cells among healthy deformable cells down to a 5% volume fraction of rigidified cells. Although aggregation of RBCs is known to affect blood rheology at low shear rates, and our simulations mimic this effect via an adhesion potential, we show that such adhesion, or aggregation, is unlikely to provide a physical rationalization for the viscosity increase observed in the experiments at moderate shear rates due to rigidified cells. Through numerical simulations, we also highlight that most of the viscosity increase of the suspension is due to the rigidity of the particles rather than their sickled or spherical shape. Our results are relevant to better characterize SCA, provide useful insights relevant to rheological consequences of blood transfusions, and, more generally, extend to the rheology of mixed suspensions having particles with different rigidities, as well as offering possibilities for developments in the field of soft material composites.

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

Primary cilia are non-motile, solitary (one per cell) microtubule-based organelles that emerge from the mother centriole after cells have exited the mitotic cycle. Identified as a mechanosensing organelle that responds to both mechanical and chemical stimuli, the primary cilium provides a fertile ground for integrative investigations of mathematical modeling, numerical simulations, and experiments. Recent experimental findings revealed considerable complexity to the underlying mechanosensory mechanisms that transmit extracellular stimuli to intracellular signaling many of which include primary cilia. In this invited review, we provide a brief survey of experimental findings on primary cilia and how these results lead to various mathematical models of the mechanics of the primary cilium bent under an external forcing such as a fluid flow or a trap. Mathematical modeling of the primary cilium as a fluid-structure interaction problem highlights the importance of basal anchorage and the anisotropic moduli of the microtubules. As theoretical modeling and numerical simulations progress, along with improved state-of-the-art experiments on primary cilia, we hope that details of ciliary regulated mechano-chemical signaling dynamics in cellular physiology will be understood in the near future.

Primary cilia have a length-dependent persistence length.

The fluctuating position of an optically trapped cilium tip under untreated and Taxol-treated conditions was used to characterize mechanical properties of the cilium axoneme and its basal body by combining experimental, analytical, and computational tools. We provide, for the first time, evidence that the persistence length of a ciliary axoneme is length-dependent; longer cilia are stiffer than shorter cilia. We demonstrate that this apparent length dependence can be understood by a combination of modeling axonemal microtubules as anisotropic elastic shells and including actomyosin-driven stochastic basal body motion. Our results also demonstrate the possibility of using observable ciliary dynamics to probe interior cytoskeletal dynamics. It is hoped that our improved characterization of cilia will result in deeper understanding of the biological function of cellular flow sensing by this organelle.

Effects of surfactant transport on electrodeformation of a viscous drop.

In this work we quantify the effects of surfactant transport on the deformation of a viscous drop under a DC electric field. We study how convective and diffusive transport of surfactants at drop surfaces influence the equilibrium and dynamic deformation of a leaky dielectric drop and a conducting drop. Focusing on the prolate drop shape (elongates along the electric field), we show the differences in equilibrium deformation and flow circulation between a leaky dielectric drop and a conducting drop. We quantify the drop electrodeformation via its dependence on the interior flow circulation and the dominant surfactant transport regime (characterized by the surface Péclet number Pe_{s}). For a leaky dielectric drop with dominant surfactant diffusion (Pe_{s}≪1), equator-to-pole (pole-to-equator) circulation yields smaller (larger) equilibrium deformation with increasing surfactant coverage, compared to a clean drop. However, when convection dominates (Pe_{s}≫1), the equilibrium drop deformation increases (decreases) with larger surfactant coverage for equator-to-pole (pole-to-equator) circulation. Larger equilibrium drop deformation is found for a leaky dielectric drop than a conducting drop when the interior flow is from equator to pole. For an interior flow from pole to equator, we identify cases where larger deformation is found for a conducting interior fluid. Finally, we study the effect of the surfactant transport on the dynamic evolution of drop shape. We found the drop undergoes an overshoot in the early deformation phase, before settling to its equilibrium shape-similar to the overshoot observed for unsteady Stokes flow.

Electrohydrodynamics of a viscous drop with inertia.

Most of the existing numerical and theoretical investigations on the electrohydrodynamics of a viscous drop have focused on the creeping Stokes flow regime, where nonlinear inertia effects are neglected. In this work we study the inertia effects on the electrodeformation of a viscous drop under a DC electric field using a novel second-order immersed interface method. The inertia effects are quantified by the Ohnesorge number Oh, and the electric field is characterized by an electric capillary number Ca_{E}. Below the critical Ca_{E}, small to moderate electric field strength gives rise to steady equilibrium drop shapes. We found that, at a fixed Ca_{E}, inertia effects induce larger deformation for an oblate drop than a prolate drop, consistent with previous results in the literature. Moreover, our simulations results indicate that inertia effects on the equilibrium drop deformation are dictated by the direction of normal electric stress on the drop interface: Larger drop deformation is found when the normal electric stress points outward, and smaller drop deformation is found otherwise. To our knowledge, such inertia effects on the equilibrium drop deformation has not been reported in the literature. Above the critical Ca_{E}, no steady equilibrium drop deformation can be found, and often the drop breaks up into a number of daughter droplets. In particular, our Navier-Stokes simulations show that, for the parameters we use, (1) daughter droplets are larger in the presence of inertia, (2) the drop deformation evolves more rapidly compared to creeping flow, and (3) complex distribution of electric stresses for drops with inertia effects. Our results suggest that normal electric pressure may be a useful tool in predicting drop pinch-off in oblate deformations.

The primary cilium is a self-adaptable, integrating nexus for mechanical stimuli and cellular signaling.

Mechanosensation is crucial for cells to sense and respond to mechanical signals within their local environment. While adaptation allows a sensor to be conditioned by stimuli within the environment and enables its operation in a wide range of stimuli intensities, the mechanisms behind adaptation remain controversial in even the most extensively studied mechanosensor, bacterial mechanosensitive channels. Primary cilia are ubiquitous sensory organelles. They have emerged as mechanosensors across diverse tissues, including kidney, liver and the embryonic node, and deflect with mechanical stimuli. Here, we show that both mechanical and chemical stimuli can alter cilium stiffness. We found that exposure to flow stiffens the cilium, which deflects less in response to subsequent exposures to flow. We also found that through a process involving acetylation, the cell can biochemically regulate cilium stiffness. Finally, we show that this altered stiffness directly affects the responsiveness of the cell to mechanical signals. These results demonstrate a potential mechanism through which the cell can regulate its mechanosensing apparatus.

Near-wall dynamics of concentrated hard-sphere suspensions: comparison of evanescent wave DLS experiments, virial approximation and simulations.

In this article we report on a study of the near-wall dynamics of suspended colloidal hard spheres over a broad range of volume fractions. We present a thorough comparison of experimental data with predictions based on a virial approximation and simulation results. We find that the virial approach describes the experimental data reasonably well up to a volume fraction of ϕ≈ 0.25 which provides us with a fast and non-costly tool for the analysis and prediction of evanescent wave DLS data. Based on this we propose a new method to assess the near-wall self-diffusion at elevated density. Here, we qualitatively confirm earlier results [Michailidou, et al., Phys. Rev. Lett., 2009, 102, 068302], which indicate that many-particle hydrodynamic interactions are diminished by the presence of the wall at increasing volume fractions as compared to bulk dynamics. Beyond this finding we show that this diminishment is different for the particle motion normal and parallel to the wall.

Efficient Brownian dynamics simulation of DNA molecules with hydrodynamic interactions in linear flows.

The coarse-grained molecular dynamics (MD) or Brownian dynamics (BD) simulation is a particle-based approach that has been applied to a wide range of biological problems that involve interactions with surrounding fluid molecules or the so-called hydrodynamic interactions (HIs). In this paper, an efficient algorithm is proposed to simulate the motion of a single DNA molecule in linear flows. The algorithm utilizes the integrating factor to cope with the effect of the linear flow of the surrounding fluid and applies the Metropolis method (MM) by Bou-Rabee, Donev, and Vanden-Eijnden [Multiscale Model. Simul. 12, 781 (2014)] to achieve more efficient BD simulation. Thus our method permits much larger time step size than previous methods while still maintaining the stability of the BD simulation, which is advantageous for long-time BD simulation. Our numerical results on λ-DNA agree very well with both experimental data and previous simulation results. Finally, when combined with fast algorithms such as the fast multipole method which has nearly optimal complexity in the total number of beads, the resulting method is parallelizable, scalable to large systems, and stable for large time step size, thus making the long-time large-scale BD simulation within practical reach. This will be useful for the study of membranes, long-chain molecules, and a large collection of molecules in the fluids.

Gating of a mechanosensitive channel due to cellular flows.

A multiscale continuum model is constructed for a mechanosensitive (MS) channel gated by tension in a lipid bilayer membrane under stresses due to fluid flows. We illustrate that for typical physiological conditions vesicle hydrodynamics driven by a fluid flow may render the membrane tension sufficiently large to gate a MS channel open. In particular, we focus on the dynamic opening/closing of a MS channel in a vesicle membrane under a planar shear flow and a pressure-driven flow across a constriction channel. Our modeling and numerical simulation results quantify the critical flow strength or flow channel geometry for intracellular transport through a MS channel. In particular, we determine the percentage of MS channels that are open or closed as a function of the relevant measure of flow strength. The modeling and simulation results imply that for fluid flows that are physiologically relevant and realizable in microfluidic configurations stress-induced intracellular transport across the lipid membrane can be achieved by the gating of reconstituted MS channels, which can be useful for designing drug delivery in medical therapy and understanding complicated mechanotransduction.

Equilibrium electrodeformation of a spheroidal vesicle in an ac electric field.

In this work, we develop a theoretical model to explain the equilibrium spheroidal deformation of a giant unilamellar vesicle (GUV) under an alternating (ac) electric field. Suspended in a leaky dielectric fluid, the vesicle membrane is modeled as a thin capacitive spheroidal shell. The equilibrium vesicle shape results from the balance between mechanical forces from the viscous fluid, the restoring elastic membrane forces, and the externally imposed electric forces. Our spheroidal model predicts a deformation-dependent transmembrane potential, and is able to capture large deformation of a vesicle under an electric field. A detailed comparison against both experiments and small-deformation (quasispherical) theory showed that the spheroidal model gives better agreement with experiments in terms of the dependence on fluid conductivity ratio, permittivity ratio, vesicle size, electric field strength, and frequency. The spheroidal model also allows for an asymptotic analysis on the crossover frequency where the equilibrium vesicle shape crosses over between prolate and oblate shapes. Comparisons show that the spheroidal model gives better agreement with experimental observations.

Dynamics of the primary cilium in shear flow.

In this work, the equilibrium shape and dynamics of a primary cilium under flow are investigated by using both theoretical modeling and experiment. The cilium is modeled as an elastic beam that may undergo large deflection due to the hydrodynamic load. Equilibrium results show that the anchoring effects of the basal body on the cilium axoneme behave as a nonlinear rotational spring. Details of the rotational spring are elucidated by coupling the elastic beam with an elastic shell. We further study the dynamics of cilium under shear flow with the cilium base angle determined from the nonlinear rotational spring, and obtain good agreement in cilium bending and relaxing dynamics when comparing between modeling and experimental results. These results potentially shed light on the physics underlying the mechanosensitive ion channel transport through the ciliary membrane.

Dynamics of a compound vesicle in shear flow.

The dynamics of a compound vesicle (a lipid bilayer membrane enclosing a fluid with a suspended particle) in shear flow is investigated by using both numerical simulations and theoretical analysis. We find that the nonlinear hydrodynamic interaction between the inclusion and the confining membrane gives rise to new features of the vesicle dynamics: The transition from tank treading to tumbling can occur in the absence of any viscosity mismatch, and a vesicle can swing if the enclosed particle is nonspherical. Our results highlight the complex effects of internal cellular structures have on cell dynamics in microcirculatory flows. For example, parasites in malaria-infected erythrocytes increase cytoplasmic viscosity, which leads to increase in blood viscosity.

Dynamics of a semiflexible polar filament in Stokes flow.

In this work, the dynamics and transport of a polarly driven filament is examined using a continuum slender-body model. Immersed in a viscous fluid, the filament gains polar propulsion from the motor proteins (anchored on the motility assay) while experiencing a viscous drag from the bottom wall. Results from the linear analysis on a straight polar filament illustrate the necessity of spatial inhomogeneity in the polar forcing for the buckling instability. The ensuing buckling leads to filament deformation, undulation, and change of its direction of motion in the numerical simulations. Repeated filament buckling in two types of motor protein concentration landscape results in diffusive transport of a polar filament on scales much larger than the mean-free path and the average duration between filament buckling events.

Hydrodynamic interactions between two semiflexible inextensible filaments in Stokes flow.

Hydrodynamic interactions between two semiflexible inextensible filaments are shown to have a significant impact on filament buckling and their subsequent motion in Stokesian fluids. In linear shear flow, hydrodynamic interactions lead to filament shear dispersion that depends on the filament aspect ratio and the initial filament separation. In linear extensional flow, hydrodynamic interactions lead to complex filament dynamics around the stagnation point. These results suggest that hydrodynamic interactions need to be taken into account to determine the self-diffusion of non-Brownian semiflexible filaments in a cellular flow [Y.-N. Young and M. J. Shelley, Phys. Rev. Lett. 99, 058303 (2007)].

Stretch-coil transition and transport of fibers in cellular flows.

It is shown that a slender elastic fiber moving in a Stokesian fluid can be susceptible to a buckling instability--termed the "stretch-coil" instability--when moving in the neighborhood of a hyperbolic stagnation point of the flow. When the stagnation point is part of an extended cellular flow, it is found that immersed fibers can move as random walkers across time-independent closed-streamline flow. It is also found that the flow is segregated into transport regions around hyperbolic stagnation points and their manifolds, and closed entrapment regions around elliptic points.