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

Automated segmentation of soft X-ray tomography: native cellular structure with sub-micron resolution at high throughput for whole-cell quantitative imaging in yeast

Soft X-ray tomography (SXT) is an invaluable tool for quantitatively analyzing cellular structures at sub-optical isotropic resolution. However, it has traditionally depended on manual segmentation, limiting its scalability for large datasets. Here, we leverage a deep learning-based auto-segmentation pipeline to segment and label cellular structures in hundreds of cells across three Saccharomyces cerevisiae strains. This task-based pipeline employs manual iterative refinement to improve segmentation accuracy for key structures, including the cell body, nucleus, vacuole, and lipid droplets, enabling high-throughput and precise phenotypic analysis. Using this approach, we quantitatively compared the 3D whole-cell morphometric characteristics of wild-type, VPH1-GFP, and vac14 strains, uncovering detailed strain-specific cell and organelle size and shape variations. We show the utility of SXT data for precise 3D curvature analysis of entire organelles and cells and detection of fine morphological features using surface meshes. Our approach facilitates comparative analyses with high spatial precision and statistical throughput, uncovering subtle morphological features at the single cell and population level. This workflow significantly enhances our ability to characterize cell anatomy and supports scalable studies on the mesoscale, with applications in investigating cellular architecture, organelle biology, and genetic research across diverse biological contexts.

Significance Statement

Soft X-ray tomography offers many powerful features for whole-cell multi-organelle imaging, but, like other high resolution volumetric imaging modalities, is typically limited by low throughput due to laborious segmentation.Auto-segmentation for soft X-ray tomography overcomes this limitation, enabling statistical 3D morphometric analysis of multiple organelles in whole cells across cell populations. The combination of high 3D resolution of SXT data with statistically useful throughput represents an avenue for more thorough characterizations of cells in toto and opens new mesoscale biological questions and statistical whole-cell modeling of organelle and cell morphology, interactions, and responses to perturbations.

Comparing Multifunctional Viral and Eukaryotic Proteins for Generating Scission Necks in Membranes

Deterministic formation of membrane scission necks by protein machinery with multiplexed functions is critical in biology. A microbial example is M2 viroporin, a proton pump from the influenza A virus that is multiplexed with membrane remodeling activity to induce budding and scission in the host membrane during viral maturation. In comparison, the dynamin family constitutes a class of eukaryotic proteins implicated in mitochondrial fission, as well as various budding and endocytosis pathways. In the case of Dnm1, the mitochondrial fission protein in yeast, the membrane remodeling activity is multiplexed with mechanoenzyme activity to create fission necks. It is not clear why these functions are combined in these scission processes, which occur in drastically different compositions and solution conditions. In general, direct experimental access to changing neck sizes induced by individual proteins or peptide fragments is challenging due to the nanoscale dimensions and influence of thermal fluctuations. Here, we use a mechanical model to estimate the size of scission necks by leveraging small-angle X-ray scattering structural data of protein-lipid systems under different conditions. The influence of interfacial tension, lipid composition, and membrane budding morphology on the size of the induced scission necks is systematically investigated using our data and molecular dynamic simulations. We find that the M2 budding protein from the influenza A virus has robust pH-dependent membrane activity that induces nanoscopic necks within the range of spontaneous hemifission for a broad range of lipid compositions. In contrast, the sizes of scission necks generated by mitochondrial fission proteins strongly depend on lipid composition, which suggests a role for mechanical constriction.

Modeling the mechanochemical feedback for membrane-protein interactions using a continuum mesh model

The Helfrich free energy is widely used to model the generation of membrane curvature due to different physical and chemical components. The governing equations resulting from the energy minimization procedure are a system of coupled higher order partial differential equations. Simulations of membrane deformation for obtaining quantitative comparisons against experimental observations require computational schemes that will allow us to solve these equations without restrictions to axisymmetric coordinates. Here, we describe one such tool that we developed in our group based on discrete differential geometry to solve these equations along with examples.

Biophysical Modeling of Actin-Mediated Structural Plasticity Reveals Mechanical Adaptation in Dendritic Spines

Synaptic plasticity is important for learning and memory formation; it describes the strengthening or weakening of connections between synapses. The postsynaptic part of excitatory synapses resides in dendritic spines, which are small protrusions on the dendrites. One of the key features of synaptic plasticity is its correlation with the size of these spines. A long-lasting synaptic strength increase [long-term potentiation (LTP)] is only possible through the reconfiguration of the actin spine cytoskeleton. Here, we develop an experimentally informed three-dimensional computational model in a moving boundary framework to investigate this reconfiguration. Our model describes the reactions between actin and actin-binding proteins leading to the cytoskeleton remodeling and their effect on the spine membrane shape to examine the spine enlargement upon LTP. Moreover, we find that the incorporation of perisynaptic elements enhances spine enlargement upon LTP, exhibiting the importance of accounting for these elements when studying structural LTP. Our model shows adaptation to repeated stimuli resulting from the interactions between spine proteins and mechanical forces.

The ins and outs of membrane bending by intrinsically disordered proteins

Membrane curvature is essential to diverse cellular functions. While classically attributed to structured domains, recent work illustrates that intrinsically disordered proteins are also potent drivers of membrane bending. Specifically, repulsive interactions among disordered domains drive convex bending, while attractive interactions drive concave bending, creating membrane-bound, liquid-like condensates. How might disordered domains that contain both repulsive and attractive domains affect curvature? Here, we examined chimeras that combined attractive and repulsive interactions. When the attractive domain was closer to the membrane, its condensation amplified steric pressure among repulsive domains, leading to convex curvature. In contrast, when the repulsive domain was closer to the membrane, attractive interactions dominated, resulting in concave curvature. Further, a transition from convex to concave curvature occurred with increasing ionic strength, which reduced repulsion while enhancing condensation. In agreement with a simple mechanical model, these results illustrate a set of design rules for membrane bending by disordered proteins.

Formation of protein-mediated bilayer tubes is governed by a snapthrough transition

Plasma membrane tubes are ubiquitous in cellular membranes and in the membranes of intracellular organelles. They play crucial roles in trafficking, ion transport, and cellular motility. These tubes can be formed due to localized forces acting on the membrane or by the curvature induced by membrane-bound proteins. Here, we present a mathematical framework to model cylindrical tubular protrusions formed by proteins that induce anisotropic spontaneous curvature. Our analysis revealed that the tube radius depends on an effective tension that includes contributions from the bare membrane tension and the protein-induced curvature. We also found that the length of the tube undergoes an abrupt transition from a short, dome-shaped membrane to a long cylinder and this transition is characteristic of a snapthrough instability. Finally, we show that the snapthrough instability depends on the different parameters including coat area, bending modulus, and extent of protein-induced curvature. Our findings have implications for tube formation due to BAR-domain proteins in processes such as endocytosis, t-tubule formation in myocytes, and cristae formation in mitochondria.

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.