Simulating membrane-bound cytoskeletal dynamics

Abstract

The cell membrane defines the shape of the cell and plays an indispensable role in bridging intra- and extra-cellular environments. The membrane, consisting of a lipid bilayer and various attaching proteins, mechanochemically interacts with the active cytoskeletal network that dynamically self-organizes, playing a vital role in cellular biomechanics and mechanosensing. Comprehensive simulations of membrane-cytoskeleton dynamics can bring insight in understanding how the cell mechanochemically responds to external signals, but a computational model that captures the complex cytoskeleton-membrane with both refined details and computational efficiency is lacking. To address this, we introduce in this thesis a triangulated membrane model and incorporate it with the active biological matter simulation platform MEDYAN ("Mechanochemical Dynamics of Active Networks"). This model accurately captures the membrane physical properties, showing how the membrane rigidity, the structure of actin networks and local chemical environments regulate the membrane deformations. Then, we present a new method for simulating membrane proteins, using stochastic reaction-diffusion sampling on unstructured membrane meshes. By incorporating a surface potential energy field into the reaction-diffusion sampling, we demonstrate rich membrane protein collective behaviors such as confined diffusion, liquid-liquid phase separation and membrane curvature sensing. Finally, in order to capture stretching, bending, shearing and twisting of actin filaments which are not all available with traditional actomyosin simulations, we introduce new finite-radius filament models based off Cosserat theory of elastic rods, with efficient implementation using finite-dimensional configurational spaces. Using the new filament models, we show that the filaments' torsional compliance can induce chiral symmetry breaking in a cross-linked actin bundle. All the new models are implemented in the MEDYAN platform, shedding light on whole cell simulations, paving way for a better understanding of the membrane-cytoskeleton system and its role in cellular dynamics.