Membrane-MEDYAN: Simulating Deformable Vesicles Containing Complex Cytoskeletal Networks

On the generation of force required for actin-based motility

The fundamental question of how forces are generated in a motile cell, a lamellipodium, and a comet tail is the subject of this note. It is now well established that cellular motility results from the polymerization of actin, the most abundant protein in eukaryotic cells, into an interconnected set of filaments. We portray this process in a continuum mechanics framework, claiming that polymerization promotes a mechanical swelling in a narrow zone around the nucleation loci, which ultimately results in cellular or bacterial motility. To this aim, a new paradigm in continuum multi-physics has been designed, departing from the well-known theory of Larché-Cahn chemo-transport-mechanics. In this note, we set up the theory of network growth and compare the outcomes of numerical simulations with experimental evidence.

Biophysical Modeling of Synaptic Plasticity

Dendritic spines are small, bulbous compartments that function as postsynaptic sites and undergo intense biochemical and biophysical activity. The role of the myriad signaling pathways that are implicated in synaptic plasticity is well studied. A recent abundance of quantitative experimental data has made the events associated with synaptic plasticity amenable to quantitative biophysical modeling. Spines are also fascinating biophysical computational units because spine geometry, signal transduction, and mechanics work in a complex feedback loop to tune synaptic plasticity. In this sense, ideas from modeling cell motility can inspire us to develop multiscale approaches for predictive modeling of synaptic plasticity. In this article, we review the key steps in postsynaptic plasticity with a specific focus on the impact of spine geometry on signaling, cytoskeleton rearrangement, and membrane mechanics. We summarize the main experimental observations and highlight how theory and computation can aid our understanding of these complex processes.

Membrane curvature catalyzes actin nucleation through nano-scale condensation of N-WASP-FBP17

Actin remodeling is spatiotemporally regulated by surface topographical cues on the membrane for signaling across diverse biological processes. Yet, the mechanism dynamic membrane curvature prompts quick actin cytoskeletal changes in signaling remain elusive. Leveraging the precision of nanolithography to control membrane curvature, we reconstructed catalytic reactions from the detection of nano-scale curvature by sensing molecules to the initiation of actin polymerization, which is challenging to study quantitatively in living cells. We show that this process occurs via topographical signal-triggered condensation and activation of the actin nucleation-promoting factor (NPF), Neuronal Wiskott-Aldrich Syndrome protein (N-WASP), which is orchestrated by curvature-sensing BAR-domain protein FBP17. Such N-WASP activation is fine-tuned by optimizing FBP17 to N-WASP stoichiometry over different curvature radii, allowing a curvature-guided macromolecular assembly pattern for polymerizing actin network locally. Our findings shed light on the intricate relationship between changes in curvature and actin remodeling via spatiotemporal regulation of NPF/BAR complex condensation.

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.

Building the next generation of virtual cells to understand cellular biology

Cell science has made significant progress by focusing on understanding individual cellular processes through reductionist approaches. However, the sheer volume of knowledge collected presents challenges in integrating this information across different scales of space and time to comprehend cellular behaviors, as well as making the data and methods more accessible for the community to tackle complex biological questions. This perspective proposes the creation of next-generation virtual cells, which are dynamic 3D models that integrate information from diverse sources, including simulations, biophysical models, image-based models, and evidence-based knowledge graphs. These virtual cells would provide statistically accurate and holistic views of real cells, bridging the gap between theoretical concepts and experimental data, and facilitating productive new collaborations among researchers across related fields.

Myosin and [Formula: see text]-actinin regulation of stress fiber contractility under tensile stress

Stress fibers are actomyosin bundles that regulate cellular mechanosensation and force transduction. Interacting with the extracellular matrix through focal adhesion complexes, stress fibers are highly dynamic structures regulated by myosin motors and crosslinking proteins. Under external mechanical stimuli such as tensile forces, the stress fiber remodels its architecture to adapt to external cues, displaying properties of viscoelastic materials. How the structural remodeling of stress fibers is related to the generation of contractile force is not well understood. In this work, we simulate mechanochemical dynamics and force generation of stress fibers using the molecular simulation platform MEDYAN. We model stress fiber as two connecting bipolar bundles attached at the ends to focal adhesion complexes. The simulated stress fibers generate contractile force that is regulated by myosin motors and [Formula: see text]-actinin crosslinkers. We find that stress fibers enhance contractility by reducing the distance between actin filaments to increase crosslinker binding, and this structural remodeling ability depends on the crosslinker turnover rate. Under tensile pulling force, the stress fiber shows an instantaneous increase of the contractile forces followed by a slow relaxation into a new steady state. While the new steady state contractility after pulling depends only on the overlap between actin bundles, the short-term contractility enhancement is sensitive to the tensile pulling distance. We further show that this mechanical response is also sensitive to the crosslinker turnover rate. Our results provide new insights into the stress fiber mechanics that have significant implications for understanding cellular adaptation to mechanical signaling.

Modeling membrane reshaping driven by dynamic protein assemblies

Remodeling of membranes in living systems is almost universally coupled to self-assembly of soluble proteins. Proteins assemble into semi-rigid shells that reshape attached membranes, and into filaments that protrude membranes. These assemblies are temporary, building from reversible protein and membrane interactions that must nucleate in the proper location. The interactions are strongly influenced by the nonequilibrium environment of the cell, such as gradients of components or active modifications by kinases. From a modeling perspective, understanding mechanisms and control thus requires 1. time-dependent approaches that ideally incorporate 2. macromolecular structure, 3. out-of-equilibrium processes, and 4. deformable membranes over microns and seconds. Realistically, tradeoffs must be made with these last three features. However, we see recent developments from the highly coarsened molecule-based scale, the continuum reaction-diffusion scale, and hybrid approaches as stimulating efforts in diverse applications. We discuss here methodological advances and progress towards simulating these processes as they occur physiologically.

Biophysical models of early mammalian embryogenesis

Embryo development is a critical and fascinating stage in the life cycle of many organisms. Despite decades of research, the earliest stages of mammalian embryogenesis are still poorly understood, caused by a scarcity of high-resolution spatial and temporal data, the use of only a few model organisms, and a paucity of truly multidisciplinary approaches that combine biological research with biophysical modeling and computational simulation. Here, we explain the theoretical frameworks and biophysical processes that are best suited to modeling the early mammalian embryo, review a comprehensive list of previous models, and discuss the most promising avenues for future work.

Non-equilibrium shapes and dynamics of active vesicles

Active vesicles, constructed through the confinement of self-propelled particles (SPPs) inside a lipid membrane shell, exhibit a large variety of non-equilibrium shapes, ranging from the formation of local tethers and dendritic conformations, to prolate and bola-like structures. To better understand the behavior of active vesicles, we perform simulations of membranes modelled as dynamically triangulated surfaces enclosing active Brownian particles. A systematic analysis of membrane deformations and SPP clustering, as a function of SPP activity and volume fraction inside the vesicle is carried out. Distributions of membrane local curvature, and the clustering and mobility of SPPs obtained from simulations of active vesicles are analysed. There exists a feedback mechanism between the enhancement of membrane curvature, the formation of clusters of active particles, and local or global changes in vesicle shape. The emergence of active tension due to the activity of SPPs can well be captured by the Young-Laplace equation. Furthermore, a simple numerical method for tether detection is presented and used to determine correlations between the number of tethers, their length, and local curvature. We also provide several geometrical arguments to explain different tether characteristics for various conditions. These results contribute to the future development of steerable active vesicles or soft micro-robots whose behaviour can be controlled and used for potential applications.

Compression-induced buckling of a semiflexible filament in two and three dimensions

The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. For filaments in the low-temperature limit and under a constant compressive force, Euler buckling theory predicts a sudden transition from a compressed to a bent state in these slender rods. In this paper, we use a mean-field theory to show that if a semiflexible chain is compressed at a finite temperature with a fixed end-to-end distance (permitting fluctuations in the compressive forces), it exhibits a continuous phase transition to a buckled state at a critical level of compression. We determine a quantitatively accurate prediction of the transverse position distribution function of the midpoint of the chain that indicates this transition. We find the mean compressive forces are non-monotonic as the extension of the filament varies, consistent with the observation that strongly buckled filaments are less able to bear an external load. We also find that for the fixed extension (isometric) ensemble, the buckling transition does not coincide with the local minimum of the mean force (in contrast to Euler buckling). We also show the theory is highly sensitive to fluctuations in length in two dimensions, and that the buckling transition can still be accurately recovered by accounting for those fluctuations. These predictions may be useful in understanding the behavior of filamentous biomolecules compressed by fluctuating forces, relevant in a variety of biological contexts.

Continuum modeling of non-conservative fluid membrane for simulating long-term cell dynamics

Living cells actively deform and move by their force generations in three-dimensional (3D) space. These 3D cell dynamics occur over a long-term time scale, ranging from tens of minutes to days. On such a time scale, turnover of cell membrane constituents due to endocytosis and exocytosis cannot be ignored, i.e., the surface membrane dynamically deforms without mass conservation. Although membrane turnover is essential for large deformation of cells, there is no computational framework yet to simulate long-term cell dynamics with a non-conservative fluidic membrane. In this paper, we proposed a computational framework for simulating the long-term dynamics of a cell membrane in 3D space. For this purpose, in the proposed framework, the cell surface membrane is treated as a viscous fluid membrane without mass conservation. Cell shape is discretized by a triangular mesh, and its dynamics are expressed by effective energy and dissipation function. The mesh structure, distorted by membrane motion, is dynamically optimized by introducing a modified dynamic remeshing method. To validate the proposed framework, numerical simulations were performed, showing that the membrane flow is reproduced in a physically consistent manner and that the artificial effects of the remeshing method were negligible. To further demonstrate the applicability of the proposed framework, numerical simulations of cell migration induced by a mechanism similar to the Marangoni effect, i.e., the polarized surface tension actively generated by the cell, were performed. The observed cell behaviors agreed with existing analytical solutions, indicating that the proposed computational framework can quantitatively reproduce long-term active cell dynamics with membrane turnover. Based on the simple description of cell membrane dynamics, this framework provides a useful basis for analyzing various cell shaping and movement.

Thermo-Statistical Effects of Inclusions on Vesicles: Division into Multispheres and Polyhedral Deformation

The construction of simple cellular models has attracted much attention as a way to explore the origin of life or elucidate the mechanisms of cell division. In the absence of complex regulatory systems, some bacteria spontaneously divide through thermostatistically elucidated mechanisms, and incorporating these simple physical principles could help to construct primitive or artificial cells. Because thermodynamic interactions play an essential role in such mechanisms, this review discusses the thermodynamic aspects of spontaneous division models of vesicles that contain a high density of inclusions, with their membrane serving as a boundary. Vesicles with highly dense inclusions are deformed according to the volume-to-area ratio. The phase separation of beads at specific intermediate volume fractions and the associated polyhedral deformation of the membrane are considered in relation to the Alder transition. Current advances in the development of a membrane-growth vesicular model are summarized. The thermostatistical understanding of these mechanisms could become a cornerstone for the construction of vesicular models that display spontaneous cell division.

Integrating cellular and molecular structures and dynamics into whole-cell models

A complete description of the state of the cell requires knowledge of its size, shape, components, intracellular reactions, and interactions with its environment-all of these as a function of time and cell growth. Adding to this list is the need for theoretical models and simulations that integrate and help to interpret this daunting amount of experimental data. It seems like an overwhelming list of requirements, but progress is being made on many fronts. In this review, we discuss the current challenges and problems in obtaining sufficient information about each aspect of a dynamical whole-cell model (DWCM) for simple and well-studied bacterial systems.

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.