Integrated intracellular organization and its variations in human iPS cells.

Understanding how a subset of expressed genes dictates cellular phenotype is a considerable challenge owing to the large numbers of molecules involved, their combinatorics and the plethora of cellular behaviours that they determine1,2. Here we reduced this complexity by focusing on cellular organization-a key readout and driver of cell behaviour3,4-at the level of major cellular structures that represent distinct organelles and functional machines, and generated the WTC-11 hiPSC Single-Cell Image Dataset v1, which contains more than 200,000 live cells in 3D, spanning 25 key cellular structures. The scale and quality of this dataset permitted the creation of a generalizable analysis framework to convert raw image data of cells and their structures into dimensionally reduced, quantitative measurements that can be interpreted by humans, and to facilitate data exploration. This framework embraces the vast cell-to-cell variability that is observed within a normal population, facilitates the integration of cell-by-cell structural data and allows quantitative analyses of distinct, separable aspects of organization within and across different cell populations. We found that the integrated intracellular organization of interphase cells was robust to the wide range of variation in cell shape in the population; that the average locations of some structures became polarized in cells at the edges of colonies while maintaining the 'wiring' of their interactions with other structures; and that, by contrast, changes in the location of structures during early mitotic reorganization were accompanied by changes in their wiring.

DLITE Uses Cell-Cell Interface Movement to Better Infer Cell-Cell Tensions.

Cell shapes and connectivities evolve over time as the colony changes shape or embryos develop. Shapes of intercellular interfaces are closely coupled with the forces resulting from actomyosin interactions, membrane tension, or cell-cell adhesions. Although it is possible to computationally infer cell-cell forces from a mechanical model of collective cell behavior, doing so for temporally evolving forces in a manner robust to digitization difficulties is challenging. Here, we introduce a method for dynamic local intercellular tension estimation (DLITE) that infers such evolution in temporal force with less sensitivity to digitization ambiguities or errors. This method builds upon previous work on single time points (cellular force-inference toolkit). We validate our method using synthetic geometries. DLITE's inferred cell colony tension evolutions correlate better with ground truth for these synthetic geometries as compared to tension values inferred from methods that consider each time point in isolation. We introduce cell connectivity errors, angle estimate errors, connection mislocalization, and connection topological changes to synthetic data and show that DLITE has reduced sensitivity to these conditions. Finally, we apply DLITE to time series of human-induced pluripotent stem cell colonies with endogenously expressed GFP-tagged zonulae occludentes-1. We show that DLITE offers improved stability in the inference of cell-cell tensions and supports a correlation between the dynamics of cell-cell forces and colony rearrangement.

Interpretable representation learning for 3D multi-piece intracellular structures using point clouds.

A key challenge in understanding subcellular organization is quantifying interpretable measurements of intracellular structures with complex multi-piece morphologies in an objective, robust and generalizable manner. Here we introduce a morphology-appropriate representation learning framework that uses 3D rotation invariant autoencoders and point clouds. This framework is used to learn representations of complex multi-piece morphologies that are independent of orientation, compact, and easy to interpret. We apply our framework to intracellular structures with punctate morphologies (e.g. DNA replication foci) and polymorphic morphologies (e.g. nucleoli). We systematically compare our framework to image-based autoencoders across several intracellular structure datasets, including a synthetic dataset with pre-defined rules of organization. We explore the trade-offs in the performance of different models by performing multi-metric benchmarking across efficiency, generative capability, and representation expressivity metrics. We find that our framework, which embraces the underlying morphology of multi-piece structures, facilitates the unsupervised discovery of sub-clusters for each structure. We show how our approach can also be applied to phenotypic profiling using a dataset of nucleolar images following drug perturbations. We implement and provide all representation learning models using CytoDL, a python package for flexible and configurable deep learning experiments.

Colony context and size-dependent compensation mechanisms give rise to variations in nuclear growth trajectories.

To investigate the fundamental question of how cellular variations arise across spatiotemporal scales in a population of identical healthy cells, we focused on nuclear growth in hiPS cell colonies as a model system. We generated a 3D timelapse dataset of thousands of nuclei over multiple days, and developed open-source tools for image and data analysis and an interactive timelapse viewer for exploring quantitative features of nuclear size and shape. We performed a data-driven analysis of nuclear growth variations across timescales. We found that individual nuclear volume growth trajectories arise from short timescale variations attributable to their spatiotemporal context within the colony. We identified a strikingly time-invariant volume compensation relationship between nuclear growth duration and starting volume across the population. Notably, we discovered that inheritance plays a crucial role in determining these two key nuclear growth features while other growth features are determined by their spatiotemporal context and are not inherited.

Mechanistic insights into actin force generation during vesicle formation from cryo-electron tomography.

Actin assembly provides force for a multitude of cellular processes. Compared to actin-assembly-based force production during cell migration, relatively little is understood about how actin assembly generates pulling forces for vesicle formation. Here, cryo-electron tomography identified actin filament number, organization, and orientation during clathrin-mediated endocytosis in human SK-MEL-2 cells, showing that force generation is robust despite variance in network organization. Actin dynamics simulations incorporating a measured branch angle indicate that sufficient force to drive membrane internalization is generated through polymerization and that assembly is triggered from ∼4 founding "mother" filaments, consistent with tomography data. Hip1R actin filament anchoring points are present along the entire endocytic invagination, where simulations show that it is key to pulling force generation, and along the neck, where it targets filament growth and makes internalization more robust. Actin organization described here allowed direct translation of structure to mechanism with broad implications for other actin-driven processes.

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.

Principles of self-organization and load adaptation by the actin cytoskeleton during clathrin-mediated endocytosis.

Force generation by actin assembly shapes cellular membranes. An experimentally constrained multiscale model shows that a minimal branched actin network is sufficient to internalize endocytic pits against membrane tension. Around 200 activated Arp2/3 complexes are required for robust internalization. A newly developed molecule-counting method determined that ~200 Arp2/3 complexes assemble at sites of clathrin-mediated endocytosis in human cells. Simulations predict that actin self-organizes into a radial branched array with growing ends oriented toward the base of the pit. Long actin filaments bend between attachment sites in the coat and the base of the pit. Elastic energy stored in bent filaments, whose presence was confirmed by cryo-electron tomography, contributes to endocytic internalization. Elevated membrane tension directs more growing filaments toward the base of the pit, increasing actin nucleation and bending for increased force production. Thus, spatially constrained actin filament assembly utilizes an adaptive mechanism enabling endocytosis under varying physical constraints.

The outer membrane of a cell is a tight but elastic barrier that controls what enters or leaves the cell. Large molecules typically cannot cross this membrane unaided. Instead, to enter the cell, they must be packaged into a pocket of the membrane that is then pulled inside. This process, called endocytosis, shuttles material into a cell hundreds of times a minute. Endocytosis relies on molecular machines that assemble and disassemble at the membrane as required. One component, a protein called actin, self-assembles near the membrane into long filaments with many repeated subunits. These filaments grow against the membrane, pulling it inwards. But it was not clear how actin filaments organize in such a way that allows them to pull on the membrane with enough force ­ and without a template to follow. Akamatsu et al. set about identifying how actin operates during endocytosis by using computer simulations that were informed by measurements made in living cells. The simulations included information about the location of actin and other essential molecules, along with the details of how these molecules work individually and together. Akamatsu et al. also developed a method to count the numbers of molecules of a key protein at individual sites of endocytosis. High-resolution imaging was then used to create 3D pictures of actin and endocytosis in action in human cells grown in the laboratory. The analysis showed the way actin filaments arrange themselves depends on the starting positions of a few key molecules that connect to actin. Imaging confirmed that, like a pole-vaulting pole, the flexible actin filaments bend to store energy and then release it to pull the membrane inwards during endocytosis. Finally, the simulations predicted that the collection of filaments adapts its shape and size in response to the resistance of the elastic membrane. This makes the system opportunistic and adaptable to the unpredictable environment within cells.

Applications and Challenges of Machine Learning to Enable Realistic Cellular Simulations.

In this perspective, we examine three key aspects of an end-to-end pipeline for realistic cellular simulations: reconstruction and segmentation of cellular structures; generation of cellular structures; and mesh generation, simulation, and data analysis. We highlight some of the relevant prior work in these distinct but overlapping areas, with a particular emphasis on current use of machine learning technologies, as well as on future opportunities.