Theory and algorithms to compute Helfrich bending forces: a review

Shaping membrane vesicles by tuning the activity of confined active polymer chains

Semi-flexible polymers, such as actin filaments, can deform the shape of membrane when confined in a membrane vesicle, playing an important role in biological processes. Here, we use dynamic Monte Carlo simulations to study an active polymer chain confined in a membrane vesicle. For flexible polymer chains, the membrane shape is governed by the competition between membrane bending rigidity and polymer activity. Stiff membrane is unaffected by small active forces, but moderate forces cause the polymer to alternate between stretched and disordered configurations, increasing the asphericity of both the polymer and the vesicle. For semi-flexible polymer chains, their stiffness can significantly impact both the vesicle and polymer shapes. We identify distinct classes of configurations that emerge as a function of polymer stiffness, membrane bending rigidity, and polymer activity. A weak polymer activity can cause the polymer to align along its contour, effectively increasing its stiffness. However, a moderate polymer activity softens the polymer chain. For membranes with low bending rigidities κ, large-scale deformations, such as wormlike or tadpole-shaped vesicles, appear at a weak polymer activity and high polymer stiffness. In the wormlike configuration, the polymer chain adopts a hairpin configuration to minimize the polymer bending energy. As the polymer stiffness increases, a tadpole-like vesicle forms, with part of the polymer deforming the membrane into a protrusion while the rest remaining confined in a bud-like structure. For stiffer membranes, we observe oblate vesicles containing toroidal polymer chains, resulting from the high cost of membrane bending energy. A moderate polymer activity causes the softening of the polymer chain, leading to a nearly spherical vesicle with slight shape fluctuation. We further characterize the order parameter of toroidal polymer chains in oblate vesicles and reveal that a slight increase in polymer activity leads to a more ordered helical structure of polymer chains.

Investigating Cell-Induced Mixing Dynamics in Microfluidic Droplets Using the Lattice Boltzmann Method

This study investigates the impact of cell dynamics on mixing efficiency within a microfluidic droplet, emphasizing the relationship between cell motion, deformability, and resultant asymmetry in velocity and concentration fields. Simulations were conducted for droplets containing encapsulated cells at varying Peclet numbers (Pe = 100-800) and coupling constants (kt = 0.0025, 0.005, 0.0075). The mixing index was significantly enhanced by the presence of the encapsulated cell, particularly at high Peclet numbers, where cell-induced disturbances in the velocity field disrupted symmetrical flow patterns, improving mixing. An asymmetry index quantified deviations in velocity and concentration fields caused by cell motion. Results revealed a complex interplay between cell deformability and fluid-cell interactions. Lower coupling constants corresponded to weaker velocity field asymmetry but, paradoxically, higher mixing indices. This counterintuitive finding was examined by analyzing the asymmetry in the x-component of the velocity field, aligned with the primary concentration gradient. Disturbances in this direction enhanced convective transport across the diffusion interface, crucial for efficient mixing. The study also examined cell trajectory and membrane deformation during droplet generation. Cells with moderate deformability exhibited greater off-center movement, leading to increased rotational dynamics and chaotic flow patterns conducive to enhanced mixing. In contrast, cells with higher or lower deformability followed more constrained paths, resulting in less effective mixing. These findings suggest optimal mixing within microfluidic droplets occurs when cells exhibit moderate deformability, balancing fluid-cell coupling and the ability to induce flow field asymmetry. These insights could inform the design of microfluidic systems for applications requiring precise mixing control, such as biomedical diagnostics and chemical synthesis.

Mitochondria and MICOS - function and modeling

An extensive review is presented on mitochondrial structure and function, mitochondrial proteins, the outer and inner membranes, cristae, the role of F1FO-ATP synthase, the mitochondrial contact site and cristae organizing system (MICOS), the sorting and assembly machinery morphology and function, and phospholipids, in particular cardiolipin. Aspects of mitochondrial regulation under physiological and pathological conditions are outlined, in particular the role of dysregulated MICOS protein subunit Mic60 in Parkinson's disease, the relations between mitochondrial quality control and proteins, and mitochondria as signaling organelles. A mathematical modeling approach of cristae and MICOS using mechanical beam theory is introduced and outlined. The proposed modeling is based on the premise that an optimization framework can be used for a better understanding of critical mitochondrial function and also to better map certain experiments and clinical interventions.

Mesoscopic elasticity controls dynamin-driven fission of lipid tubules

Mesoscale physics bridges the gap between the microscopic degrees of freedom of a system and its large-scale continuous behavior and highlights the role of a few key quantities in complex and multiscale phenomena, like dynamin-driven fission of lipid membranes. The dynamin protein wraps the neck formed during clathrin-mediated endocytosis, for instance, and constricts it until severing occurs. Although ubiquitous and fundamental for life, the cooperation between the GTP-consuming conformational changes within the protein and the full-scale response of the underlying lipid substrate is yet to be unraveled. In this work, we build an effective mesoscopic model from constriction to fission of lipid tubules based on continuum membrane elasticity and implicitly accounting for ratchet-like power strokes of dynamins. Localization of the fission event, the overall geometry, and the energy expenditure we predict comply with the major experimental findings. This bolsters the idea that a continuous picture emerges soon enough to relate dynamin polymerization length and membrane rigidity and tension with the optimal pathway to fission. We therefore suggest that dynamins found in in vivo processes may optimize their structure accordingly. Ultimately, we shed light on real-time conductance measurements available in literature and predict the fission time dependency on elastic parameters.

Mesoscale simulation of biomembranes with FreeDTS

We present FreeDTS software for performing computational research on biomembranes at the mesoscale. In this software, a membrane is represented by a dynamically triangulated surface equipped with vertex-based inclusions to integrate the effects of integral and peripheral membrane proteins. Several algorithms are included in the software to simulate complex membranes at different conditions such as framed membranes with constant tension, vesicles and high-genus membranes with various fixed volumes or constant pressure differences and applying external forces to membrane regions. Furthermore, the software allows the user to turn off the shape evolution of the membrane and focus solely on the organization of proteins. As a result, we can take realistic membrane shapes obtained from, for example, cryo-electron tomography and backmap them into a finer simulation model. In addition to many biomembrane applications, this software brings us a step closer to simulating realistic biomembranes with molecular resolution. Here we provide several interesting showcases of the power of the software but leave a wide range of potential applications for interested users.

The local variation of the Gaussian modulus enables different pathways for fluid lipid vesicle fusion

Viral infections, fertilization, neurotransmission, and many other fundamental biological processes rely on membrane fusion. Straightforward calculations based on the celebrated Canham-Helfrich elastic model predict a large topological energy barrier that prevents the fusion process from being thermally activated. While such high energy is in accordance with the physical barrier function of lipid membranes, it is difficult to reconcile with the biological mechanisms involved in fusion processes. In this work, we use a Ginzburg-Landau type of free energy that recovers the Canham-Helfrich model in the limit of small width-to-vesicle-extension ratio, with the additional ability to handle topological transitions. We show that a local modification of the Gaussian modulus in the merging region both dramatically lowers the elastic energy barrier and substantially changes the minimal energy pathway for fusion, in accordance with experimental evidence. Therefore, we discuss biological examples in which such a modification might play a crucial role.

Forceful patterning: theoretical principles of mechanochemical pattern formation

Biological pattern formation is essential for generating and maintaining spatial structures from the scale of a single cell to tissues and even collections of organisms. Besides biochemical interactions, there is an important role for mechanical and geometrical features in the generation of patterns. We review the theoretical principles underlying different types of mechanochemical pattern formation across spatial scales and levels of biological organization.

Effects of membrane viscoelasticity on the red blood cell dynamics in a microcapillary

The mechanical properties of red blood cells (RBCs) play key roles in their biological functions in microcirculation. In particular, RBCs must deform significantly to travel through microcapillaries with sizes comparable with or even smaller than their own. Although the dynamics of RBCs in microcapillaries have received considerable attention, the effect of membrane viscoelasticity has been largely overlooked. In this work, we present a computational study based on the boundary integral method and thin-shell mechanics to examine how membrane viscoelasticity influences the dynamics of RBCs flowing through straight and constricted microcapillaries. Our results reveal that the cell with a viscoelastic membrane undergoes substantially different motion and deformation compared with results based on a purely elastic membrane model. Comparisons with experimental data also suggest the importance of accounting for membrane viscoelasticity to properly capture the transient dynamics of an RBC flowing through a microcapillary. Taken together, these findings demonstrate the significant effects of membrane viscoelasticity on RBC dynamics in different microcapillary environments. The computational framework also lays the groundwork for more accurate quantitative modeling of the mechanical response of RBCs in their mechanotransduction process in subsequent investigations.

The effects of internal forces and membrane heterogeneity on three-dimensional cell shapes

The shape of cells and the control thereof plays a central role in a variety of cellular processes, including endo- and exocytosis, cell division and cell movement. Intra- and extracellular forces control the shapes, and while the shape changes in some processes such as exocytosis are intracellularly-controlled and localized in the cell, movement requires force transmission to the environment, and the feedback from it can affect the cell shape and mode of movement used. The shape of a cell is determined by its cytoskeleton (CSK), and thus shape changes involved in various processes involve controlled remodeling of the CSK. While much is known about individual components involved in these processes, an integrated understanding of how intra- and extracellular signals are coupled to the control of the mechanical changes involved is not at hand for any of them. As a first step toward understanding the interaction between intracellular forces imposed on the membrane and cell shape, we investigate the role of distributed surrogates for cortical forces in producing the observed three-dimensional shapes. We show how different balances of applied forces lead to such shapes, that there are different routes to the same end state, and that state transitions between axisymmetric shapes need not all be axisymmetric. Examples of the force distributions that lead to protrusions are given, and the shape changes induced by adhesion of a cell to a surface are studied. The results provide a reference framework for developing detailed models of intracellular force distributions observed experimentally, and provide a basis for studying how movement of a cell in a tissue or fluid is influenced by its shape.

Activation energy and force fields during topological transitions of fluid lipid vesicles

Topological transitions of fluid lipid membranes are fundamental processes for cell life. For example, they are required for endo- and exocytosis or to enable neurotransmitters to cross the neural synapses. Here, inspired by the idea that fusion and fission proteins could have evolved in Nature in order to carry out a minimal work expenditure, we evaluate the minimal free energy pathway for the transition between two spherical large unilamellar vesicles and a dumbbell-shaped one. To address the problem, we propose and successfully use a Ginzburg-Landau type of free energy, which allows us to uniquely describe without interruption the whole, full-scale topological change. We also compute the force fields needed to overcome the involved energy barriers. The obtained forces are in excellent agreement, in terms of intensity, scale, and spatial localization with experimental data on typical fission protein systems, whereas they suggest the presence of additional features in fusion proteins.

Mem3DG: Modeling membrane mechanochemical dynamics in 3D using discrete differential geometry

Biomembranes adopt varying morphologies that are vital to cellular functions. Many studies use computational modeling to understand how various mechanochemical factors contribute to membrane shape transformations. Compared with approximation-based methods (e.g., finite element method [FEM]), the class of discrete mesh models offers greater flexibility to simulate complex physics and shapes in three dimensions; its formulation produces an efficient algorithm while maintaining coordinate-free geometric descriptions. However, ambiguities in geometric definitions in the discrete context have led to a lack of consensus on which discrete mesh model is theoretically and numerically optimal; a bijective relationship between the terms contributing to both the energy and forces from the discrete and smooth geometric theories remains to be established. We address this and present an extensible framework, Mem3DG, for modeling 3D mechanochemical dynamics of membranes based on discrete differential geometry (DDG) on triangulated meshes. The formalism of DDG resolves the inconsistency and provides a unifying perspective on how to relate the smooth and discrete energy and forces. To demonstrate, Mem3DG is used to model a sequence of examples with increasing mechanochemical complexity: recovering classical shape transformations such as 1) biconcave disk, dumbbell, and unduloid; and 2) spherical bud on spherical, flat-patch membrane; investigating how the coupling of membrane mechanics with protein mobility jointly affects phase and shape transformation. As high-resolution 3D imaging of membrane ultrastructure becomes more readily available, we envision Mem3DG to be applied as an end-to-end tool to simulate realistic cell geometry under user-specified mechanochemical conditions.

A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flows

When red blood cells (RBCs) experience non-physiologically high stresses, e.g., in medical devices, they can rupture in a process called hemolysis. Directly simulating this process is computationally unaffordable given that the length scales of a medical device are several orders of magnitude larger than that of a RBC. To overcome this separation of scales, the present work introduces an affordable computational framework that accurately resolves the stress and deformation of a RBC in a spatially and temporally varying macroscale flow field such as those found in a typical medical device. The underlying idea of the present framework is to treat RBCs as one-way coupled tracers in the macroscale flow by capturing the effect of the flow on their dynamics but neglecting their effect on the flow at the macroscale. As a result, the RBC dynamics are simulated after those of the flow in a postprocessing step by receiving the fluid velocity gradient tensor measured along the RBC trajectory as the input. To resolve the fluid velocity in the immediate vicinity of the RBC as well as the motion of the membrane, we employ the boundary integral method coupled to a structural solver. The governing equations are discretized in space using spherical harmonics, yielding spectral integration accuracy. The predictions produced by this formulation are in good agreement with those obtained from simulations of spherical capsules in shear flows and optical tweezers experiments. The accuracy of the present method is evaluated using unbounded shear flow as a benchmark. Its computational cost grows proportional to p 5, where p is the degree of the spherical harmonic. It also exhibits a fast convergence rate that is approximately O ( p 6 ) for p ⪅ 20.

Accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit, and customized 16-bit number formats

Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here we evaluate the possibility to use even FP16 and posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32. For this, we first show that the commonly occurring number range in the LBM is a lot smaller than the FP16 number range. Based on this observation, we develop customized 16-bit formats-based on a modified IEEE-754 and on a modified posit standard-that are specifically tailored to the needs of the LBM. We then carry out an in-depth characterization of LBM accuracy for six different test systems with increasing complexity: Poiseuille flow, Taylor-Green vortices, Karman vortex streets, lid-driven cavity, a microcapsule in shear flow (utilizing the immersed-boundary method), and, finally, the impact of a raindrop (based on a volume-of-fluid approach). We find that the difference in accuracy between FP64 and FP32 is negligible in almost all cases, and that for a large number of cases even 16-bit is sufficient. Finally, we provide a detailed performance analysis of all precision levels on a large number of hardware microarchitectures and show that significant speedup is achieved with mixed FP32/16-bit.

Experimental Investigations on the Conductance of Lipid Membranes under Differential Hydrostatic Pressure

The unassisted transport of inorganic ions through lipid membranes has become increasingly relevant to an expansive range of biological phenomena. Recent simulations indicate a strong influence of a lipid membrane's curvature on its permeability, which may be part of the overall cell sensitivity to mechanical stimulation. However, most ionic permeability experiments employ a flat, uncurved lipid membrane, which disregards the physiological relevance of curvature on such investigations. To fill this gap in our knowledge, we adapted a traditional experimental system consisting of a planar lipid membrane, which we exposed to a controlled, differential hydrostatic pressure. Our electrophysiology experiments indicate a strong correlation between the changes in membrane geometry elicited by the application of pressure, as inferred from capacitance measurements, and the resulting conductance. Our experiments also confirmed the well-established influence of cholesterol addition to lipid membranes in adjusting their mechanical properties and overall permeability. Therefore, the proposed experimental system may prove useful for a better understanding of the intricate connections between membrane mechanics and adjustments of cellular functionalities upon mechanical stimulation, as well as for confirmation of predictions made by simulations and theoretical modeling.

Squeezing multiple soft particles into a constriction: Transition to clogging

We study numerically how multiple deformable capsules squeeze into a constriction. This situation is largely encountered in microfluidic chips designed to manipulate living cells, which are soft entities. We use fully three-dimensional simulations based on the lattice Boltzmann method to compute the flow of the suspending fluid and on the immersed boundary method to achieve the two-way fluid-structure interaction. The mechanics of the capsule membrane elasticity is computed with the finite-element method. We obtain two main states: continuous passage of the particles and their blockage that leads to clogging the constriction. The transition from one state to another is dictated by the ratio between the size of the capsules and the constriction width and by the capsule membrane deformability. The latter is found to enhance particle passage through narrower constrictions, where rigid particles with similar diameter are blocked and lead to clogging.

Red blood cell shape transitions and dynamics in time-dependent capillary flows

The dynamics of single red blood cells (RBCs) determine microvascular blood flow by adapting their shape to the flow conditions in the narrow vessels. In this study, we explore the dynamics and shape transitions of RBCs on the cellular scale under confined and unsteady flow conditions using a combination of microfluidic experiments and numerical simulations. Tracking RBCs in a comoving frame in time-dependent flows reveals that the mean transition time from the symmetric croissant to the off-centered, nonsymmetric slipper shape is significantly faster than the opposite shape transition, which exhibits pronounced cell rotations. Complementary simulations indicate that these dynamics depend on the orientation of the RBC membrane in the channel during the time-dependent flow. Moreover, we show how the tank-treading movement of slipper-shaped RBCs in combination with the narrow channel leads to oscillations of the cell's center of mass. The frequency of these oscillations depends on the cell velocity, the viscosity of the surrounding fluid, and the cytosol viscosity. These results provide a potential framework to identify and study pathological changes in RBC properties.

Biomembranes undergo complex, non-axisymmetric deformations governed by Kirchhoff-Love kinematicsand revealed by a three-dimensional computational framework

Biomembranes play a central role in various phenomena like locomotion of cells, cell-cell interactions, packaging and transport of nutrients, transmission of nerve impulses, and in maintaining organelle morphology and functionality. During these processes, the membranes undergo significant morphological changes through deformation, scission, and fusion. Modelling the underlying mechanics of such morphological changes has traditionally relied on reduced order axisymmetric representations of membrane geometry and deformation. Axisymmetric representations, while robust and extensively deployed, suffer from their inability to model-symmetry breaking deformations and structural bifurcations. To address this limitation, a three-dimensional computational mechanics framework for high fidelity modelling of biomembrane deformation is presented. The proposed framework brings together Kirchhoff–Love thin-shell kinematics, Helfrich-energy-based mechanics, and state-of-the-art numerical techniques for modelling deformation of surface geometries. Lipid bilayers are represented as spline-based surface discretizations immersed in a three-dimensional space; this enables modelling of a wide spectrum of membrane geometries, boundary conditions, and deformations that are physically admissible in a three-dimensional space. The mathematical basis of the framework and its numerical machinery are presented, and their utility is demonstrated by modelling three classical, yet non-trivial, membrane deformation problems: formation of tubular shapes and their lateral constriction, Piezo1-induced membrane footprint generation and gating response, and the budding of membranes by protein coats during endocytosis. For each problem, the full three-dimensional membrane deformation is captured, potential symmetry-breaking deformation paths identified, and various case studies of boundary and load conditions are presented. Using the endocytic vesicle budding as a case study, we also present a ‘phase diagram’ for its symmetric and broken-symmetry states.

Exploring new roles for actin upon LTP induction in dendritic spines

Dendritic spines, small protrusions of the dendrites, enlarge upon LTP induction, linking morphological and functional properties. Although the role of actin in spine enlargement has been well studied, little is known about its relationship with mechanical membrane properties, such as membrane tension, which is involved in many cell processes, like exocytosis. Here, we use a 3D model of the dendritic spine to investigate how polymerization of actin filaments can effectively elevate the membrane tension to trigger exocytosis in a domain close to the tip of the spine. Moreover, we show that the same pool of actin promotes full membrane fusion after exocytosis and spine stabilization.

Numerical simulations of capsule deformation using a dual time-stepping lattice Boltzmann method

In this work a quasisteady, dual time-stepping lattice Boltzmann method is proposed for simulation of capsule deformation. At each time step the steady-state lattice Boltzmann equation is solved using the full approximation storage multigrid scheme for nonlinear equations. The capsule membrane is modeled as an infinitely thin shell suspended in an ambient fluid domain with the fluid structure interaction computed using the immersed boundary method. A finite element method is used to compute the elastic forces exerted by the capsule membrane. Results for a wide range of parameters and initial configurations are presented. The proposed method is found to reduce the computational time by a factor of ten.

Vortical flow structures induced by red blood cells in capillaries

OBJECTIVE

Knowledge about the flow field of the plasma around the red blood cells in capillary flow is important for a physical understanding of blood flow and the transport of micro- and nanoparticles and molecules in the flowing plasma. We conducted an experimental study on the flow field around red blood cells in capillary flow that is complemented by simulations of vortical flow between red blood cells.

METHODS

Red blood cells were injected in a 10 × 12 µm rectangular microchannel at a low hematocrit, and the flow field around one or two cells was captured by a high-speed camera that tracked 250 nm nanoparticles in the flow field, acting as tracers.

RESULTS

While the flow field around a steady "croissant" shape is found to be similar to that of a rigid sphere, the flow field around a "slipper" shape exhibits a small vortex at the rear of the red blood cell. Even more pronounced are vortex-like structures observed in the central region between two neighboring croissants.

CONCLUSIONS

The rotation frequency of the vortices is to a good approximation, inversely proportional to the distance between the cells. Our experimental data are complemented by numerical simulations.

Static interface profiles for contact lines on an elastic membrane with the Willmore energy

We consider a fluid interface in contact with an elastic membrane and study the static profiles of the interface and the membrane. Equilibrium conditions are derived by minimizing the total energy of the system with volume constraints. The total energy consists of surface energies and the Willmore energy; the latter penalizes the bending of the membrane. It is found that, while the membrane is locally flat at the contact line with the contact angle satisfying the Young-Dupré equation, the gradient of the mean curvature of the membrane exhibits a jump across the contact line. This jump balances the surface tension of the fluid interface in the normal direction of the membrane. Asymptotic solutions are obtained for two-dimensional systems in the limits as the reduced bending modulus ν tends to +∞ and 0, respectively. In the stiff limit as ν→+∞, the leading-order solution is given by that of a droplet sitting on a rigid substrate with the contact angle satisfying the Young-Dupré equation; in contrast, in the soft limit as ν→0, a transition layer appears near the contact line and the interfaces have constant curvatures in the outer region with apparent contact angles obeying Neumann's law. These solutions are validated by numerical experiments.

Interfacial Curvature as a Potential Index for Prognosis of Colon Adenocarcinoma

Tumor invasion and metastasis are complex interfacial mechanical processes between the tumor and its surrounding tissue, with the interfacial curvature of tumor playing an important role in cancer progression. In this study, the potential role of interfacial curvature in the prognosis of patients with colon adenocarcinoma is investigated. The front edge interfacial curvature of adenocarcinoma from biopsies of patients in different tumor, lymph node, and metastasis (TNM) stages are calculated and compared, and prognosis assessment is conducted using Kaplan-Meier and Cox proportional hazards regression analyses. Results reveal that patients with larger interfacial curvature of adenocarcinoma are more likely to belong to higher TNM stages. Concomitantly, in the same TNM stage, patients with increased adenocarcinoma interfacial curvature show worse prognosis with higher recurrence and lower survival rates. Besides, interfacial curvature is an independent prognostic factor for cause-specific survival and relapse-free survival among all selected patients. Mechanical models of colon adenocarcinoma invasion and metastasis are established to better understand the close association between interfacial curvature and tumor progression. The results together with hematoxylin and eosin staining indicate that metastasis in stages T3N0M0 and T3N1M0 may be linked to large interfacial curvatures. Therefore, interfacial curvature may serve as a potential index for predicting prognosis in patients with colon adenocarcinoma.

Quantifying the roles of space and stochasticity in computer simulations for cell biology and cellular biochemistry

Most of the fascinating phenomena studied in cell biology emerge from interactions among highly organized multimolecular structures embedded into complex and frequently dynamic cellular morphologies. For the exploration of such systems, computer simulation has proved to be an invaluable tool, and many researchers in this field have developed sophisticated computational models for application to specific cell biological questions. However, it is often difficult to reconcile conflicting computational results that use different approaches to describe the same phenomenon. To address this issue systematically, we have defined a series of computational test cases ranging from very simple to moderately complex, varying key features of dimensionality, reaction type, reaction speed, crowding, and cell size. We then quantified how explicit spatial and/or stochastic implementations alter outcomes, even when all methods use the same reaction network, rates, and concentrations. For simple cases, we generally find minor differences in solutions of the same problem. However, we observe increasing discordance as the effects of localization, dimensionality reduction, and irreversible enzymatic reactions are combined. We discuss the strengths and limitations of commonly used computational approaches for exploring cell biological questions and provide a framework for decision making by researchers developing new models. As computational power and speed continue to increase at a remarkable rate, the dream of a fully comprehensive computational model of a living cell may be drawing closer to reality, but our analysis demonstrates that it will be crucial to evaluate the accuracy of such models critically and systematically.

Tools for computational analysis of moving boundary problems in cellular mechanobiology

A cell's ability to change shape is one of the most fundamental biological processes and is essential for maintaining healthy organisms. When the ability to control shape goes awry, it often results in a diseased system. As such, it is important to understand the mechanisms that allow a cell to sense and respond to its environment so as to maintain cellular shape homeostasis. Because of the inherent complexity of the system, computational models that are based on sound theoretical understanding of the biochemistry and biomechanics and that use experimentally measured parameters are an essential tool. These models involve an inherent feedback, whereby shape is determined by the action of regulatory signals whose spatial distribution depends on the shape. To carry out computational simulations of these moving boundary problems requires special computational techniques. A variety of alternative approaches, depending on the type and scale of question being asked, have been used to simulate various biological processes, including cell motility, division, mechanosensation, and cell engulfment. In general, these models consider the forces that act on the system (both internally generated, or externally imposed) and the mechanical properties of the cell that resist these forces. Moving forward, making these techniques more accessible to the non-expert will help improve interdisciplinary research thereby providing new insight into important biological processes that affect human health. This article is categorized under: Cancer > Cancer>Computational Models Cancer > Cancer>Molecular and Cellular Physiology.

Single bubble and drop techniques for characterizing foams and emulsions

The physics of foams and emulsions has traditionally been studied using bulk foam/emulsion tests and single film platforms such as the Scheludko cell. Recently there has been a renewed interest in a third class of techniques that we term as single bubble/drop tests, which employ isolated whole bubbles and drops to probe the characteristics of foams and emulsions. Single bubble and drop techniques provide a convenient framework for investigating a number of important characteristics of foams and emulsions, including the rheology, stabilization mechanisms, and rupture dynamics. In this review we provide a comprehensive discussion of the various single bubble/drop platforms and the associated experimental measurement protocols including the construction of coalescence time distributions, visualization of the thin film profiles and characterization of the interfacial rheological properties. Subsequently, we summarize the recent developments in foam and emulsion science with a focus on the results obtained through single bubble/drop techniques. We conclude the review by presenting important venues for future research.

Mechanical design of apertures and the infolding of pollen grain

When pollen grains become exposed to the environment, they rapidly desiccate. To protect themselves until rehydration, the grains undergo characteristic infolding with the help of special structures in the grain wall-apertures-where the otherwise thick exine shell is absent or reduced in thickness. Recent theoretical studies have highlighted the importance of apertures for the elastic response and the folding of the grain. Experimental observations show that different pollen grains sharing the same number and type of apertures can nonetheless fold in quite diverse fashions. Using the thin-shell theory of elasticity, we show how both the absolute elastic properties of the pollen wall and the relative elastic differences between the exine wall and the apertures play an important role in determining pollen folding upon desiccation. Focusing primarily on colpate pollen, we delineate the regions of pollen elastic parameters where desiccation leads to a regular, complete closing of all apertures and thus to an infolding which protects the grain against water loss. Phase diagrams of pollen folding pathways indicate that an increase in the number of apertures leads to a reduction of the region of elastic parameters where the apertures close in a regular fashion. The infolding also depends on the details of the aperture shape and size, and our study explains how the features of the mechanical design of apertures influence the pollen folding patterns. Understanding the mechanical principles behind pollen folding pathways should also prove useful for the design of the elastic response of artificial inhomogeneous shells.

Active particles induce large shape deformations in giant lipid vesicles

Biological cells generate intricate structures by sculpting their membrane from within to actively sense and respond to external stimuli or to explore their environment^1-4. Several pathogenic bacteria also provide examples of how localized forces strongly deform cell membranes from inside, leading to the invasion of neighbouring healthy mammalian cells⁵. Giant unilamellar vesicles have been successfully used as a minimal model system with which to mimic biological cells^6-11, but the realization of a minimal system with localized active internal forces that can strongly deform lipid membranes from within and lead to dramatic shape changes remains challenging. Here we present a combined experimental and simulation study that demonstrates how self-propelled particles enclosed in giant unilamellar vesicles can induce a plethora of non-equilibrium shapes and active membrane fluctuations. Using confocal microscopy, in the experiments we explore the membrane response to local forces exerted by self-phoretic Janus microswimmers. To quantify dynamic membrane changes, we perform Langevin dynamics simulations of active Brownian particles enclosed in thin membrane shells modelled by dynamically triangulated surfaces. The most pronounced shape changes are observed at low and moderate particle loadings, with the formation of tether-like protrusions and highly branched, dendritic structures, whereas at high volume fractions globally deformed vesicle shapes are observed. The resulting state diagram predicts the conditions under which local internal forces generate various membrane shapes. A controlled realization of such distorted vesicle morphologies could improve the design of artificial systems such as small-scale soft robots and synthetic cells.

Flow-accelerated platelet biogenesis is due to an elasto-hydrodynamic instability

Blood platelets are formed by fragmentation of long membrane extensions from bone marrow megakaryocytes in the blood flow. Using lattice-Boltzmann/immersed boundary simulations we propose a biological Rayleigh-Plateau instability as the biophysical mechanism behind this fragmentation process. This instability is akin to the surface tension-induced breakup of a liquid jet but is driven by active cortical processes including actomyosin contractility and microtubule sliding. Our fully three-dimensional simulations highlight the crucial role of actomyosin contractility, which is required to trigger the instability, and illustrate how the wavelength of the instability determines the size of the final platelets. The elasto-hydrodynamic origin of the fragmentation explains the strong acceleration of platelet biogenesis in the presence of an external flow, which we observe in agreement with experiments. Our simulations then allow us to disentangle the influence of specific flow conditions: While a homogeneous flow with uniform velocity leads to the strongest acceleration, a shear flow with a linear velocity gradient can cause fusion events of two developing platelet-sized swellings during fragmentation. A fusion event may lead to the release of larger structures which are observable as preplatelets in experiments. Together, our findings strongly indicate a mainly physical origin of fragmentation and regulation of platelet size in flow-accelerated platelet biogenesis.

Newton-Cartan submanifolds and fluid membranes

We develop the geometric description of submanifolds in Newton-Cartan spacetime. This provides the necessary starting point for a covariant spacetime formulation of Galilean-invariant hydrodynamics on curved surfaces. We argue that this is the natural geometrical framework to study fluid membranes in thermal equilibrium and their dynamics out of equilibrium. A simple model of fluid membranes that only depends on the surface tension is presented and, extracting the resulting stresses, we show that perturbations away from equilibrium yield the standard result for the dispersion of elastic waves. We also find a generalization of the Canham-Helfrich bending energy for lipid vesicles that takes into account the requirements of thermal equilibrium.

Modeling the Shape of Synaptic Spines by Their Actin Dynamics

Dendritic spines are the morphological basis of excitatory synapses in the cortex and their size and shape correlates with functional synaptic properties. Recent experiments show that spines exhibit large shape fluctuations that are not related to activity-dependent plasticity but nonetheless might influence memory storage at their synapses. To investigate the determinants of such spontaneous fluctuations, we propose a mathematical model for the dynamics of the spine shape and analyze it in 2D-related to experimental microscopic imagery-and in 3D. We show that the spine shape is governed by a local imbalance between membrane tension and the expansive force from actin bundles that originates from discrete actin polymerization foci. Experiments have shown that only few such polymerization foci co-exist at any time in a spine, each having limited life time. The model shows that the momentarily existing set of such foci pushes the membrane along certain directions until foci are replaced and other directions may now be affected. We explore these relations in depth and use our model to predict shape and temporal characteristics of spines from the different biophysical parameters involved in actin polymerization. Approximating the model by a single recursive equation we finally demonstrate that the temporal evolution of the number of active foci is sufficient to predict the size of the model-spines. Thus, our model provides the first platform to study the relation between molecular and morphological properties of the spine with a high degree of biophysical detail.

Protein Recruitment through Indirect Mechanochemical Interactions

Some of the key proteins essential for important cellular processes are capable of recruiting other proteins from the cytosol to phospholipid membranes. The physical basis for this cooperativity of binding is, surprisingly, still unclear. Here, we suggest a general feedback mechanism that explains cooperativity through mechanochemical coupling mediated by the mechanical properties of phospholipid membranes. Our theory predicts that protein recruitment, and therefore also protein pattern formation, involves membrane deformation and is strongly affected by membrane composition.

Antimargination of Microparticles and Platelets in the Vicinity of Branching Vessels

We investigate the margination of microparticles/platelets in blood flow through complex geometries typical for in vivo vessel networks: a vessel confluence and a bifurcation. Using three-dimensional lattice Boltzmann simulations, we confirm that behind the confluence of two vessels, a cell-free layer devoid of red blood cells develops in the channel center. Despite its small size of roughly 1 μm, this central cell-free layer persists for up to 100 μm after the confluence. Most importantly, we show from simulations that this layer also contains a significant amount of microparticles/platelets and validate this result by in vivo microscopy in mouse venules. At bifurcations, however, a similar effect does not appear, and margination is largely unaffected by the geometry. This antimargination toward the vessel center after a confluence may explain earlier in vivo observations, which found that platelet concentrations near the vessel wall are seen to be much higher on the arteriolar side (containing bifurcations) than on the venular side (containing confluences) of the vascular system.

Frequency-dependent higher-order Stokes singularities near a planar elastic boundary: Implications for the hydrodynamics of an active microswimmer near an elastic interface

The emerging field of self-driven active particles in fluid environments has recently created significant interest in the biophysics and bioengineering communities owing to their promising future for biomedical and technological applications. These microswimmers move autonomously through aqueous media, where under realistic situations they encounter a plethora of external stimuli and confining surfaces with peculiar elastic properties. Based on a far-field hydrodynamic model, we present an analytical theory to describe the physical interaction and hydrodynamic couplings between a self-propelled active microswimmer and an elastic interface that features resistance toward shear and bending. We model the active agent as a superposition of higher-order Stokes singularities and elucidate the associated translational and rotational velocities induced by the nearby elastic boundary. Our results show that the velocities can be decomposed in shear and bending related contributions which approach the velocities of active agents close to a no-slip rigid wall in the steady limit. The transient dynamics predict that contributions to the velocities of the microswimmer due to bending resistance are generally more pronounced than those due to shear resistance. Bending can enhance (suppress) the velocities resulting from higher-order singularities whereas the shear related contribution decreases (increases) the velocities. Most prominently, we find that near an elastic interface of only energetic resistance toward shear deformation, such as that of an elastic capsule designed for drug delivery, a swimming bacterium undergoes rotation of the same sense as observed near a no-slip wall. In contrast to that, near an interface of only energetic resistance toward bending, such as that of a fluid vesicle or liposome, we find a reversed sense of rotation. Our results provide insight into the control and guidance of artificial and synthetic self-propelling active microswimmers near elastic confinements.

Mesoscale simulation of soft particles with tunable contact angle in multicomponent fluids

Soft particles at fluid interfaces play an important role in many aspects of our daily life, such as the food industry, paints and coatings, and medical applications. Analytical methods are not capable of describing the emergent effects of the complex dynamics of suspensions of many soft particles, whereas experiments typically either only capture bulk properties or require invasive methods. Computational methods are therefore a great tool to complement experimental work. However, an efficient and versatile numerical method is needed to model dense suspensions of many soft particles. In this article we propose a method to simulate soft particles in a multicomponent fluid, both at and near fluid-fluid interfaces, based on the lattice Boltzmann method, and characterize the error stemming from the fluid-structure coupling for the particle equilibrium shape when adsorbed onto a fluid-fluid interface. Furthermore, we characterize the influence of the preferential contact angle of the particle surface and the particle softness on the vertical displacement of the center of mass relative to the fluid interface. Finally, we demonstrate the capability of our model by simulating a soft capsule adsorbing onto a fluid-fluid interface with a shear flow parallel to the interface, and the covering of a droplet suspended in another fluid by soft particles with different wettability.

Computational modeling of active deformable membranes embedded in three-dimensional flows

Active gel theory has recently been very successful in describing biologically active materials such as actin filaments or moving bacteria in temporally fixed and simple geometries such as cubes or spheres. Here we develop a computational algorithm to compute the dynamic evolution of an arbitrarily shaped, deformable thin membrane of active material embedded in a three-dimensional flowing liquid. For this, our algorithm combines active gel theory with the classical theory of thin elastic shells. To compute the actual forces resulting from active stresses, we apply a parabolic fitting procedure to the triangulated membrane surface. Active forces are then dynamically coupled via an immersed-boundary method to the surrounding fluid whose dynamics can be solved by any standard, e.g., Lattice-Boltzmann, flow solver. We validate our algorithm using the Green's functions of Berthoumieux et al. [New J. Phys. 16, 065005 (2014)10.1088/1367-2630/16/6/065005] for an active cylindrical membrane subjected (i) to a locally increased active stress and (ii) to a homogeneous active stress. For the latter scenario, we predict in addition a nonaxisymmetric instability. We highlight the versatility of our method by analyzing the flow field inside an actively deforming cell embedded in external shear flow. Further applications may be cytoplasmic streaming or active membranes in blood flows.

Glycan-decorated protocells: novel features for rebuilding cellular processes

In synthetic biology approaches, lipid vesicles are widely used as protocell models. While many compounds have been encapsulated in vesicles (e.g. DNA, cytoskeleton and enzymes), the incorporation of glycocalyx components in the lipid bilayer has attracted much less attention so far. In recent years, glycoconjugates have been integrated in the membrane of giant unilamellar vesicles (GUVs). These minimal membrane systems have largely contributed to shed light on the molecular mechanisms of cellular processes. In this review, we first introduce several preparation and biophysical characterization methods of GUVs. Then, we highlight specific applications of protocells investigating glycolipid-mediated endocytosis of toxins, viruses and bacteria. In addition, we delineate how prototissues have been assembled from glycan-decorated protocells by using lectin-mediated cross-linking of opposed glycoreceptors (e.g. glycolipids and glycopeptides). In future applications, glycan-decorated protocells might be useful for investigating cell-cell interactions (e.g. adhesion and communication). We also speculate about the implication of lectin-glycoreceptor interactions in membrane fusion processes.

Effect of particles with repulsive interactions enclosed in both rigid spherical shells and flexible fluid vesicles studied by Monte Carlo simulation

Experimental observations indicate that the repulsion of particles is a factor that induces the transformation of vesicles containing multiple particles. Metropolis Monte Carlo simulations are performed with two models in which repulsive particles are enclosed inside a vesicle. The distribution of the particles and the effective bending coefficient and surface tension of the membrane are analyzed. The shape and internal structure of the vesicle containing the particles are investigated as the vesicle volume is decreased. It is revealed that the repulsive interaction between particles produces a layered structure and stiffens the membrane. When particles repulsively interact over a long range, the membrane takes on a dumbbell form.

Hydrodynamic mobility of a solid particle near a spherical elastic membrane. II. Asymmetric motion

In this paper, we derive analytical expressions for the leading-order hydrodynamic mobility of a small solid particle undergoing motion tangential to a nearby large spherical capsule whose membrane possesses resistance toward shearing and bending. Together with the results obtained in the first part [Daddi-Moussa-Ider and Gekle, Phys. Rev. E 95, 013108 (2017)2470-004510.1103/PhysRevE.95.013108], where the axisymmetric motion perpendicular to the capsule membrane is considered, the solution of the general mobility problem is thus determined. We find that shearing resistance induces a low-frequency peak in the particle self-mobility, resulting from the membrane normal displacement in the same way, although less pronounced, to what has been observed for the axisymmetric motion. In the zero-frequency limit, the self-mobility correction near a hard sphere is recovered only if the membrane has a nonvanishing resistance toward shearing. We further compute the in-plane mean-square displacement of a nearby diffusing particle, finding that the membrane induces a long-lasting subdiffusive regime. Considering capsule motion, we find that the correction to the pair-mobility function is solely determined by membrane shearing properties. Our analytical calculations are compared and validated with fully resolved boundary integral simulations where a very good agreement is obtained.