Flow-induced channel formation in the cytoplasm of motile cells

Simulations of dynamically cross-linked actin networks: Morphology, rheology, and hydrodynamic interactions

Cross-linked actin networks are the primary component of the cell cytoskeleton and have been the subject of numerous experimental and modeling studies. While these studies have demonstrated that the networks are viscoelastic materials, evolving from elastic solids on short timescales to viscous fluids on long ones, questions remain about the duration of each asymptotic regime, the role of the surrounding fluid, and the behavior of the networks on intermediate timescales. Here we perform detailed simulations of passively cross-linked non-Brownian actin networks to quantify the principal timescales involved in the elastoviscous behavior, study the role of nonlocal hydrodynamic interactions, and parameterize continuum models from discrete stochastic simulations. To do this, we extend our recent computational framework for semiflexible filament suspensions, which is based on nonlocal slender body theory, to actin networks with dynamic cross linkers and finite filament lifetime. We introduce a model where the cross linkers are elastic springs with sticky ends stochastically binding to and unbinding from the elastic filaments, which randomly turn over at a characteristic rate. We show that, depending on the parameters, the network evolves to a steady state morphology that is either an isotropic actin mesh or a mesh with embedded actin bundles. For different degrees of bundling, we numerically apply small-amplitude oscillatory shear deformation to extract three timescales from networks of hundreds of filaments and cross linkers. We analyze the dependence of these timescales, which range from the order of hundredths of a second to the actin turnover time of several seconds, on the dynamic nature of the links, solvent viscosity, and filament bending stiffness. We show that the network is mostly elastic on the short time scale, with the elasticity coming mainly from the cross links, and viscous on the long time scale, with the effective viscosity originating primarily from stretching and breaking of the cross links. We show that the influence of nonlocal hydrodynamic interactions depends on the network morphology: for homogeneous meshworks, nonlocal hydrodynamics gives only a small correction to the viscous behavior, but for bundled networks it both hinders the formation of bundles and significantly lowers the resistance to shear once bundles are formed. We use our results to construct three-timescale generalized Maxwell models of the networks.

Uni-cellular integration of complex spatial information in slime moulds and ciliates

Single-celled organisms show a fascinating faculty for integrating spatial information and adapting their behaviour accordingly. As such they are of potential interest for elucidating fundamental mechanisms of developmental biology. In this mini-review we highlight current research on two organisms, the true slime mould Physarum polycephalum and the ciliates Paramecium and Tetrahymena. For each of these, we present a case study how applying physical principles to explain behaviour can lead to the understanding of general principles possibly relevant to developmental biology.

A lumped parameter model of endoplasm flow in Physarum polycephalum explains migration and polarization-induced asymmetry during the onset of locomotion

The plasmodial slime mold Physarum polycephalum exhibits strong, periodic flow of cytoplasm through the veins of its network. In the special case of mesoplasmodia, a newly described starvation-induced, shape-constant morphotype, this periodic endoplasm streaming is the basis of locomotion. Furthermore, we presume that cytoplasm flow is also involved in signal transmission and signal processing. Mesoplasmodia motility resembles amoeboid locomotion. In contrast to other amoebae, however, mesoplasmodia move without extending pseudopods and retain a coherent, fan-shaped morphology throughout their steady locomotion. Attaining sizes of up to 2 mm², mesoplasmodia are also much bigger than other amoebae. We characterize this particular type of locomotion and identify patterns of movement. By using the analogy between pulsatile fluid flow through a network of elastic tubes and electrical circuits, we build a lumped model that explains observed fluid flow patterns. Essentially, the mesoplasmodium acts as a low-pass filter, permitting only low-frequency oscillations to propagate from back to front. This frequency selection serves to optimize flow and reduces power dissipation. Furthermore, we introduce a distributed element into the lumped model to explain cell polarization during the onset of chemotaxis: Biochemical cues (internal or external) lead to a local softening of the actin cortex, which in turn causes an increased flow of cytoplasm into that area and, thus, a net forward movement. We conclude that the internal actin-enclosed vein network gives the slime mold a high measure of control over fluid transport, especially by softening or hardening, which in turn leads to polarization and net movement.

Image-based model of the spectrin cytoskeleton for red blood cell simulation

We simulate deformable red blood cells in the microcirculation using the immersed boundary method with a cytoskeletal model that incorporates structural details revealed by tomographic images. The elasticity of red blood cells is known to be supplied by both their lipid bilayer membranes, which resist bending and local changes in area, and their cytoskeletons, which resist in-plane shear. The cytoskeleton consists of spectrin tetramers that are tethered to the lipid bilayer by ankyrin and by actin-based junctional complexes. We model the cytoskeleton as a random geometric graph, with nodes corresponding to junctional complexes and with edges corresponding to spectrin tetramers such that the edge lengths are given by the end-to-end distances between nodes. The statistical properties of this graph are based on distributions gathered from three-dimensional tomographic images of the cytoskeleton by a segmentation algorithm. We show that the elastic response of our model cytoskeleton, in which the spectrin polymers are treated as entropic springs, is in good agreement with the experimentally measured shear modulus. By simulating red blood cells in flow with the immersed boundary method, we compare this discrete cytoskeletal model to an existing continuum model and predict the extent to which dynamic spectrin network connectivity can protect against failure in the case of a red cell subjected to an applied strain. The methods presented here could form the basis of disease- and patient-specific computational studies of hereditary diseases affecting the red cell cytoskeleton.

Red blood cells are responsible for delivering oxygen to tissues throughout the body. These terminally differentiated cells have developed a fascinating flexibility and resiliency that is critical to navigating the circulatory system. Far from being rigid bodies, red blood cells adopt biconcave disk shapes at equilibrium, parachute-like shapes as they move between large vessels and small capillaries, and more extreme shapes as they traverse the endothelial slits of the spleen. Understanding the remarkable mechanical properties that allow red cells to experience such large deformations while maintaining structural integrity is a fundamental question in physiology that may help advance treatments of genetic disorders such as hereditary spherocytosis and elliptocytosis that affect red cell flexibility and can lead to severe anemia. In this work, we present a model of the red blood cell cytoskeleton based on cryoelectron tomography data. We develop an image processing technique to gather statistics from these data and use these statistics to generate a random entropic network to model the cytoskeleton. We then simulate the behavior of the resulting red blood cells in flow. As we demonstrate through simulations, this method makes it possible to examine the consequences of changes in microstructural properties such as the rate of cytoskeletal remodeling.

Current reinforcement model reproduces center-in-center vein trajectory of Physarum polycephalum

Vein networks span the whole body of the amoeboid organism in the plasmodial slime mould Physarum polycephalum, and the network topology is rearranged within an hour in response to spatio-temporal variations of the environment. It has been reported that this tube morphogenesis is capable of solving mazes, and a mathematical model, named the 'current reinforcement rule', was proposed based on the adaptability of the veins. Although it is known that this model works well for reproducing some key characters of the organism's maze-solving behaviour, one important issue is still open: In the real organism, the thick veins tend to trace the shortest possible route by cutting the corners at the turn of corridors, following a center-in-center trajectory, but it has not yet been examined whether this feature also appears in the mathematical model, using corridors of finite width. In this report, we confirm that the mathematical model reproduces the center-in-center trajectory of veins around corners observed in the maze-solving experiment.

Patterns of cell thickness oscillations during directional migration of Physarum polycephalum

The functional relationship between the velocity of cell locomotion and intracellular spatial patterns of thickness oscillations in the acellular slime mould Physarum polycephalum was studied. The freely migrating plasmodial cells of 300-800 µm length were tadpole-shaped and displayed thickness oscillations along their longitudinal (body) axis. Two distinct patterns of intracellular thickness oscillations were observed in dependence on the locomotive velocity. The first mode consisted of a single travelling wave that propagated from the rear to the front of the cell. This pattern occurred when the plasmodium migrated slowly. The second mode was a multinodal standing wave that was found during events of fast propagation. Transitions between these two types of cell thickness oscillation patterns took place in narrow propagation velocity intervals. We discuss the possible mechanism leading to these patterns, which are conjectured to modulate both the intracellular pressure and the velocity of free locomotion of the cell.

Nonlinearity in cytoplasm viscosity can generate an essential symmetry breaking in cellular behaviors

The cytoplasms of ameboid cells are nonlinearly viscous. The cell controls this viscosity by modulating the amount, localization and interactions of bio-polymers. Here we investigated how the nonlinearity infers the cellular behaviors and whether nonlinearity-specific behaviors exist. We modeled the developed plasmodium of the slime mold Physarum polycephalum as a network of branching tubes and examined the linear and nonlinear viscous cytoplasm flows in the tubes. We found that the nonlinearity in the cytoplasm׳s viscosity induces a novel type of symmetry breaking in the protoplasmic flow. We also show that symmetry breaking can play an important role in adaptive behaviors, namely, connection of behavioral modes implemented on different time scales and transportation of molecular signals from the front to the rear of the cell during cellular locomotion.

An active poroelastic model for mechanochemical patterns in protoplasmic droplets of Physarum polycephalum

Motivated by recent experimental studies, we derive and analyze a two-dimensional model for the contraction patterns observed in protoplasmic droplets of Physarum polycephalum. The model couples a description of an active poroelastic two-phase medium with equations describing the spatiotemporal dynamics of the intracellular free calcium concentration. The poroelastic medium is assumed to consist of an active viscoelastic solid representing the cytoskeleton and a viscous fluid describing the cytosol. The equations for the poroelastic medium are obtained from continuum force balance and include the relevant mechanical fields and an incompressibility condition for the two-phase medium. The reaction-diffusion equations for the calcium dynamics in the protoplasm of Physarum are extended by advective transport due to the flow of the cytosol generated by mechanical stress. Moreover, we assume that the active tension in the solid cytoskeleton is regulated by the calcium concentration in the fluid phase at the same location, which introduces a mechanochemical coupling. A linear stability analysis of the homogeneous state without deformation and cytosolic flows exhibits an oscillatory Turing instability for a large enough mechanochemical coupling strength. Numerical simulations of the model equations reproduce a large variety of wave patterns, including traveling and standing waves, turbulent patterns, rotating spirals and antiphase oscillations in line with experimental observations of contraction patterns in the protoplasmic droplets.