Modeling the Excess Cell Surface Stored in a Complex Morphology of Bleb-Like Protrusions
[['Cells maintain their structural integrity while being flexible enough to adopt a variety of shapes .', 'In general , it is the cytoskeleton of eukaryotic cells that drives shape transformation leading to cell movement and provides the structural support to the cytoplasm and the means to resist external forces .', 'The periphery of cells , consisting of the plasma membrane ( PM ) and the acto-myosin cortex , is highly dynamic to accommodate shape change .', 'The plasma membrane ( PM ) consists of a high density of proteins [1] embedded in a phospholipid bilayer of 5–10 nm thickness , with a very limited ability to extend without rupture [2 , 3] but highly amenable to bending [4 , 5 , 6] .', 'The thin ( 50–500 nm ) layer of cytoskeleton structure immediately subjacent to the plasma membrane , known as the cell cortex , consists of a dense F-actin network that is cross-linked by actin binding proteins and is amenable to contractility mediated by myosin motors .', 'Interposed between the cortex and the PM is a thin spectrin-actin network , forming a ‘fishnet’ with a mesh size of ~100 nm [7 , 8] .', 'This structure is anchored both to the PM and cortex by adaptor proteins .', 'In the following , we term the plasma membrane and spectrin-actin network as the “cell surface” .', 'Previously we [9] suggested that most dynamical shape changes exhibited by non-spread ( rounded ) cells originate from a membrane-cortex folding-unfolding process and an excess of cell surface area is a necessary requirement for such changes .', 'We investigated the dynamics of periodically protruding cells and hypothesized that the plasma membrane and thin cortical layer remain coupled during all stages of shape transformation .', 'We also assumed that densely compressed cell surface folds and small protrusions could be kept intact by the underlying actin-myosin network residing in the cortex proper .', 'While this notion may be applicable to many shape transitions occurring in non-spread cells , in this paper we reconsider this hypothesis in context of one of the most drastic changes of cell shape: the transformation from a fully spread to a rounded state .', 'If a cell transitions from a spread to rounded state while maintaining a constant volume , it will experience an excess of surface area over the minimum needed to cover the enclosed volume .', 'Because this process typically happens rapidly ( ~30s- ) , there is insufficient time for excess membrane to be internalized by endocytosis .', 'Thus , another mechanism for storing surface area at the plasma membrane must exist .', 'Indeed , there is significant evidence from both electron and fluorescence microscopy that during the rounding process the cell surface adopts a tightly folded morphology [9 , 10 , 11 , 12] .', 'While there are a number of models for cell shape , most of them treat the cell surface as smooth [13 , 14 , 15 , 16 , 17] and do not take into account the possibility that rounded cells store excess surface area in a dense distribution of bleb-like protrusions ( BLiPs ) .', 'Thus , new modeling approaches are needed to understand the dynamics of cell shape changes that involve active use of this surface storage .', 'We introduce three complementary modeling approaches , each incorporating the concept of excess surface area .', 'The first approach is a 2D model based on a thin cell surface structure that is coupled to a thicker , contractile actomyosin layer .', 'This model allows us to investigate the folding of the excess surface and to estimate the bending energy in different configurations .', 'The second approach is a random “seed and growth” model that produces 3D morphologies consistent with the distributions of BLiP size and number estimated from scanning electron micrographs .', 'This model yields insight into how large numbers of BLiPs are efficiently packed on the cell periphery .', 'The third approach is a multi-compartment phase field model .', 'By faithfully capturing the physical properties of the cortex , cytosol , and cell surface , the model predicts the stress and pressure distributions associated with a highly folded 2D morphology and a dense distribution of 3D BLIPs .', 'Phase field models have been widely used to study complex systems comprised of distinct material phases and their adjacent interfaces .', 'When the separate material phases are immiscible , the phase field approach is to prescribe a finite thickness of a “diffuse” interface within which there is a mixture of the two materials [18 , 19] .', 'The phase field method is an alternative to sharp interface methods; in both methods the shape and evolution of the sharp versus diffuse interface are part of the solution .', 'For every pair of adjacent material components , a phase field variable is introduced that interpolates from one material phase to the other through the finite thickness boundary .', 'Phase field models have been employed to describe shapes of lipid bilayer vesicles in which the surface tension and Helfrich bending energy are approximated using a bulk energy defined within the diffuse interfacial layer [19] .', 'Phase field models have been applied to many interfacial problems including liquid drops , multiphase complex fluids [20] , and fractures in solid-state materials [21] .', 'The phase field model simulations achieve separate goals .', 'From either a 2D transmission electron micrograph or a 3D image reconstruction of the cell morphology , the model “learns” the spontaneous curvature functional of the rounded , BLiP-rich , morphology .', 'Since the phase field model faithfully captures material properties of each cellular compartment , the model converges to the cell target morphology while constructing self-consistent stress and isotropic pressure distributions for the cell surface , cortex and cytoplasm , as well as estimating the nematic orientation within the cortex .', 'Storage of excess cell surface in folds or bleb-like protrusions at the periphery is likely to be important for a variety of rapid cell shape changes , taking place over a time scale of a few minutes or less , such as those that occur in forms of amoeboid migration or either within or in the transitions between the phases of mitosis .', 'It seems likely that rapid cell shape changes can be accomplished more quickly by calling upon a reserve of excess membrane stored in the BLiP distributions rather than relying on extensive membrane-cortex remodeling and exocytosis .', 'Thus , the theoretical approaches presented here should be applicable in a number of different biological contexts .'], ['When spread cells ( Fig 1A ) are chemically detached from an underlying substrate , they rapidly transition to a rounded state on a characteristic time scale of ~30-60s ( Fig 1B ) .', 'Numerous studies suggest that in media with constant osmolarity , cell volume is stable [3 , 22] .', 'We estimated cell volume by reconstructing 3D geometries from Z-stacks of spinning disc fluorescence images of cells undergoing rounding .', 'The mean volume for Chinese hamster ovary ( CHO ) cells in the spread state is 6 . 5±2 . 82*103 μm3 , while the mean volume in the rounded state is 5 . 7±2 . 30*103 μm3 , indicating a slight decrease in cell volume after rounding .', 'Because this slight decrease in cell volume would only increase excess surface area , in all our models , we assume that cell volume remains constant during rounding .', 'The surface area of a spread cell is estimated as twice the area measured from images to account for dorsal and ventral surfaces .', 'In reality cells are not completely flat and have more surface area due to finite thickness , particularly around the nucleus .', 'Therefore , we are underestimating the surface area of a spread cell .', 'In the rounded state , the minimal surface area needed to enclose the measured cell volume can be found by assuming the cell is spherical and calculating the radius .', 'For example , the cell 1 in Fig 1A has a surface area of ~ 31000 μm2 while the surface area needed to enclose the rounded state is only ~ 2200 μm2 .', 'Therefore after rounding , this cell has ~ 14 times more surface area than is required to enclose its volume .', 'This image presents an extreme case of surface area excess .', 'For cell 2 in Fig 1A , which is less spread before detachment , an excess surface area of about five times the required amount is accumulated following rounding .', 'It is important to note that the amount of excess surface area that is accumulated during rounding depends on cell type and characteristics of the spreading and detachment for individual cells .', 'Using DIC and fluorescence microscopy , we studied populations of cells before and after detachment , and individual cells rounding during trypsinization .', 'The histogram in Fig 1C presents the distribution of surface areas for spread CHO cells ( blue bars; population mean = 4310±3600 μm2 , N = 199 ) and cells immediately after rounding ( red bars; population mean = 892±284 μm2 , N = 1646 ) .', 'The distribution of rounded cell sizes is narrow with majority of cell radii ( Fig 1E ) being between 7 and 9 . 5 μm ( mean = 8 . 36±1 . 24 μm; N = 1646 ) .', 'Separate experiments , where we followed the change in morphology of individual cells during rounding , demonstrated that for CHO cells the average excess surface , defined as the ratio of spread cell area to that required to smoothly cover a sphere with radius corresponding to that of the rounded cell , accumulated due to detachment and rounding is 3 . 8± 2 . 06 with maximum value of 12 ( N = 99 ) .', 'To gain insight into how much excess surface area can be stored in BLiPs , we first consider the case of a rounded cell uniformly covered with equally sized sphererical BLiPs .', 'It is easy to show ( S1 Appendix ) that as spherical BLiPs become smaller , more excess surface area can be accommodated .', 'The maximum possible surface excess that can be stored in the equally sized sphererical BLiPs is 5 ( in the limit of BliP radius r→0 ) .', 'The fact that we observed rounded cells with the surface excess as high as 14 , means that cells utilize a more efficient packing strategy .', 'Also no limits on the surface excess ratio would be imposed if we did not require BLiPs to be spherical , but rather allow for an arbitrarily high curvature of the surface , as occurs , for example , in tubules .', 'Yet , the majority of BLiPs appear to be rounded immediately after detachment .', 'These considerations suggest that the actual morphology of the folded cell surface is dictated by a balance between the necessity to pack tightly a very large number of BLiPs and the necessity to generate , regulate , and maintain significant surface curvature .', 'To better understand the process of packing cell surface excess into a convoluted surface morphology , we constructed a 2D geometric model designed to produce BLiPs .', 'We hypothesize that the cell surface and underlying contractile cortex form a two-layer structure that is coupled at certain fixed points .', 'The first layer , which we term the cell surface , is passive and consists of the plasma membrane and membrane associated cytoskeleton .', 'This layer is assumed to be similar to the spectrin-actin network that is coupled to the plasma membrane of red cells [23] .', 'Such structures have been shown by Kusumi and co-workers to exist in many other cell types [8] .', 'The membrane associated cytoskeleton has been termed the membrane skeleton fence [8] .', 'It is basically a very thin filamentous meshwork that provides a “fishnet” with a mesh size of approximately 100 nm immediately underlying the PM .', 'This layer is coupled to the plasma membrane via adaptor proteins including the ankyrin and ERM families as well as by interactions of the membrane skeleton fence with lipids in the inner monolayer of the PM .', 'The layer is thought to be passive undergoing only thermal motions , serving to anchor some transmembrane proteins and restrict the free diffusion of others .', 'We assume that this thin layer is coupled via adaptor proteins to a thicker , active contractile layer containing actin and myosin .', 'This view of the cortex-cell surface couple is consistent with that advanced by Charras et al ( 2006 ) [24] in the context of spontaneously blebbing cells .', 'Additional evidence for this structure comes by imaging employing confocal and electron microscopy .', 'Fig 1F shows a confocal image of the actin-myosin cortex in rounded cells as visualized with GFP-lifeact and RFP-myosin merged with a DIC image of the same cell .', 'The green signal for Lifeact marks F-actin filaments associated with the folded morphology of the cell periphery and this fluorescent signal originates from both the thin layer immediately subjacent to the membrane and the thicker contractile layer .', 'The fluorescent signal from myosin ( red ) shows that this protein is localized mainly to a thin circle located below the BLiPs and more toward the cell interior .', 'S1 Fig presents Z-stack images of the same cell .', 'From this image , it is clearly visible that the convoluted morphology covers the whole cell .', 'Fig 1G and 1H shows immunogold TEM images of GFP-lifeact where F-actin is seen underlying BLiPs ( arrow ) and also in a layer closer toward the center of the cell ( arrowhead ) .', 'To construct the 2D geometric model introduced qualitatively above , we implement a two-layer architecture in a 2D bead-spring model of the cell membrane and cortex ( model description in Methods ) .', 'The bead-spring model consists of two-layers ( Fig 2A and S2 Fig ) , where one layer ( outer layer ) represents the membrane and underlying actin mesh ( i . e . the cell surface ) and the other layer ( inner layer ) represents the actomyosin-rich contractile cortex .', 'In each layer , beads are connected pairwise by springs and contact points serve to connect the two layers .', 'By minimizing the bending energy we explored the steady-state shapes of BLiPs generated during cell rounding when the cell is rapidly presented with a substantial excess surface .', 'For simplicity we define excess surface ratio ( ER ) as a ratio of perimeters of the surface layer and contracted cortex .', 'Here we define a normalized total bending energy ( E ) as for', '( 1 )', 'where κ is the local curvature measured between two neighboring beads , L is the perimeter , and N is the number of beads in the outer layer [25 , 26 , 27 , 28 , 29] .', 'The number of contact points determines the number of folds ( M ) .', 'In the simulation the total Hamiltonian of the two-layer system is minimized with the result that a folded configuration of outer layer is produced .', 'Fig 2B shows the resulting shapes as a function of both M and ER .', 'Fig 2C shows a portion of a model cell with BLiPs at steady state where the gold line represents the contractile part of the cortex with contact points .', 'While the appearance of folds is expected , the shape of folds and the bending energy stored in each configuration is of particular interest .', 'Fig 2D gives a comparison of the fold configuration for several different sets of parameters with the calculated bending energy for each shape .', 'Inspection of the fold shapes shows that in order to accommodate more surface , the folds tend to develop long necks .', '( Note that in case where there is heterogeneity in the size of folds , this effect would allow small folds to grow under larger ones , an effect that permits accommodation of more excess surface . )', 'The smallest possible bending energy will be achieved when ER = 1 and M = 0 ( no surface excess and no BLiPs ) .', 'For a given value of ER , the energy increases with the number of BLiPs ( Fig 2D and S2 Fig ) .', 'However , the bending energy is decreasing while the excess surface ratio is increasing .', 'Although this result looks counterintuitive , it can be explained .', 'The local curvature is the inverse of local radius .', 'Folds with a longer perimeter have bigger inner radii which substantially decreases bending energy with the square of local radius ( Fig 2D and S3 Fig ) .', 'Thus , morphologies with longer perimeters corresponding to larger ERs will have lower energy compared to the shapes with the same number of folds but with shorter perimeters ( i . e . smaller ERs ) .', 'The analysis of the area which is stored inside the folds ( i . e . , volume in 3D ) shows that for the same surface excess , more area is stored in folds when the number of folds used for accommodation of this surface surplus is smaller ( S3D Fig ) .', 'At the resolution achievable by standard fluorescence microscopy , the convoluted cell surface often appears as a thickening of membrane and cortical stains ( Fig 3A and 3B ) .', 'However , surface morphology can be imaged at higher resolution using both scanning ( SEM ) and thin-section transmission ( TEM ) electron micrographs ( Figs 3C and 1H , respectively ) .', 'At this scale , bleb-like protrusions ( BLiPs ) and other cell surface protuberances are clearly visible .', 'Note that only fully spread cells have a smooth surface essentially devoid of protrusions ( S4 Fig ) .', 'To determine length , area and volume metrics of BLiPs , we manually segmented SEM images of cells that were fixed after rounding ( S5 Fig ) .', 'Each protrusion was approximated as a sphere and the area of the protrusion visible on the image was interpreted as a two dimensional projection of that sphere .', 'We calculated the radius that corresponds to a projection of that size , and consider it as the radius of the BLiP .', 'The distributions of BLiP radii derived from 10 SEM images ( 25 cells ) that include 7096 BLiPs is presented in Fig 3D .', 'We find that the distribution of radii is skewed with a preponderance of small BLiPs and a decreasing frequency of larger BLiPs .', 'The mean BLiP radius is R = 0 . 25 μm with a median of 0 . 22 μm and mode of 0 . 19 μm .', 'It is important to mention that during the processes of detachment and rounding some part of the cell surface can be lost due to incomplete detachment from the substrate or because it remains in retraction fibers .', 'However , the area remaining in retraction fibers is quite small .', 'Using SEM images from cells that spread for 24 h and rounded 5 minutes before fixation , we estimated that the surface area that might be stored in retraction fibers represents between ~0 . 5–5% of the cell surface area in the spread state .', 'To investigate how the large number of BLiPs required to accommodate the excess cell surface are packed on the cell periphery , we constructed two models .', 'As a plausible starting point , we employed a Voronoi approach , in which a spherical ball of radius R=S/4π contains the surface area , S of the spread cell before detachment and rounding; n seed point locations are randomly sampled from a uniform spatial distribution on the ball .', 'The ball is then partitioned by a Voronoi tessellation according to the n seed points ( Fig 4A ) such that any point in each Voronoi cell is closer to the parent seed point than any other seed point .', 'The area of each Voronoi cell is then determined .', 'In this configuration , we assumed that , upon the cell rounding to its final state with the “BLiPed” morphology , each Voronoi cell of area v morphs into a spherical BLiP with radius r=ν/4π .', 'Fig 4D demonstrates that , constructed in this way , the distribution of BLiP radii has a well-defined length scale with the bell-shape distribution , which is not consistent with the skewed distribution of experimental data .', 'This result arises from the fact that Voronoi cells corresponding to two very closely positioned seeds are not necessarily small themselves , as might be expected from two closely positioned BLiPs .', 'In order to mitigate this effect , we introduce an alternative 3D “seed and growth” model ( Fig 4B and 4C ) , in which BLiP radii are proportional to spacing between randomly distributed seeds .', 'In this model , spheres are generated from each seed point by increasing their radii at a uniform rate .', 'Simultaneously , the locations of the seed points are moved outward radially at the same rate , so that the spheres always remain tangent to the cell .', 'When one sphere encounters another , it stops growing .', 'When all spheres have stopped growth , a spherical cell coated by different sized BLiPs is produced ( Fig 4B ) .', 'The resulting BLiP radius distribution in Fig 4D is more consistent with the experimental distribution than that produced by the Voronoi model .', 'The generated structure is also consistent with SEM image data , which show approximately spherical BLiPs largely covering the cell but with some areas devoid of BLiPs .', 'We reproduced 2D cross-sectional views from the simulated 3D geometries ( Fig 4F ) ; these show similarity to the thin section TEM images of rounded cells ( Fig 4G ) where some of the BLiPs appear to be detached from the cell body because BLiPs are not always sectioned through their centers .', 'In the “seed and growth” model , a larger number of BLiPs results in a higher surface area excess ratio ( Fig 4E ) and a smaller percentage of the cell volume stored within the BLiPs , which is consistent with our simplified estimations based on an 2D equal-sized BLiP distribution .', 'In principle , BLiPs that are not of equal size could allow a more efficient packing of excessive cell surface ( with smaller BLiPs filling the space between larger ones ) , which is important for accommodating very high excess ratios ( >5 ) .', 'However , in this model the packing is still inefficient because it always generates areas devoid of BLiPs .', 'A potential improvement to our model might be to incorporate stochastic seeding of new BLiPs and occasional “shrinking” BLiPs that are in contact , so that BLiPs are dynamic and continue to adjust themselves toward the most efficient filling of the available space .', 'Such a “seed , growth , and shrinking” model would be consistent with the BLiP dynamics observed in our experiments , but is beyond the scope of the current paper .', 'Another potential improvement would be to make final BLiP size proportional to the rate of expansion of the BLiP; this has been found to be the case in an earlier study of blebbing cells [30] .', 'The preceding models help build mechanistic intuition , yet , while predictive , they do not capture all of the essential physics of the rounded phenotype .', 'In order to approach this goal , we formulated 2D and 3D phase field models for a cell immersed in the aqueous extracellular environment .', 'The model is formulated in 3 space dimensions ( 3D ) , but it also restricts to 2D for purposes of modeling a cell cross-section .', 'In our case , we have three phases ( exterior aqueous medium , cortex , interior cytosol ) and two diffuse interfaces .', 'The external aqueous medium and interior cytosol are modeled as viscous fluids with specified viscosities and the cortex is modeled as a nematic ( liquid crystal ) gel [31] .', 'The first diffuse ( i . e . , finite thickness ) interface is what we have termed the cell surface , consisting of the plasma membrane and very thin underlying filamentous “fishnet” , that separates the aqueous medium and the cortex proper .', 'As described below in Methods , a particular level set function in the phase field formulation will afford our definition of the “cell surface” .', 'The second diffuse interface is the cortex-cytosol transition layer .', 'Fig 5 is a 2D schematic of a 3D cell cross-section with individual components and diffuse interfaces labeled , along with the phase variables defined below .', '( We do not explicitly model the nucleus within the cytoplasm for this paper since we are primarily concerned with the stationary rounded morphology . )', 'A more complete mechanical formulation giving the total system free energy in terms of its components is found in Methods .', 'We summarize the key numerical results of the phase field modeling of a 2D cell surface morphology due to an imposed excess arc length enclosing the 2D area .', 'We require a 2D image of the membrane morphology , taken from 2D transmission electron micrographs .', 'From the image file , we posit an initial smooth membrane boundary , and then evolve the phase field model while adjusting the spontaneous curvature function C1 until the model converges to the image dataset .', 'We first illustrate the ability of the phase field model to match an arbitrary specified 2D boundary by “learning” the spontaneous curvature function; the results are shown in Fig 6A for an illustrative benchmark in which the cell perimeter contains 25 regularly spaced , uniform “BLiPs” .', '( In 2D , this is achieved by superimposing the appropriate Fourier mode on a circle . )', 'Next , we used as input the actual periphery of a rounded cell from a 2D transmission electron micrograph ( TEM ) image .', 'The results in Fig 6B show the convergence of the phase field membrane morphology to the TEM image , where the nematic phase ordering ( representing F-actin orientation ) in the cortex is depicted .', 'This result assumes tangential anchoring condition of F-actin at both cortical diffuse interfaces .', 'Fig 6C–6E shows the phase field predictions of the pressure distribution ( C ) and the first invariant of the dominant stored stress , the Ericksen stress , for the converged stationary morphology ( D , E ) shown in Fig 6B .', 'These results reveal the orders of magnitude as well as spatial localization of pressure and stored stress for the target 2D morphology .', 'A 3D simulation is depicted in Fig 7 to demonstrate the capability of our phase field model to converge to a target 3D cell surface morphology .', 'The excess surface area ratio for this illustration is s0 = 3 .', 'Because it is impossible to reproduce a full 3D morphology of a rounded cell from a single scanning electron micrograph , we use the “seed and growth” model to simulate a 3D cell surface target .', 'This model ( Fig 4 ) provides a 3D surface morphology consistent with the measured BLiP size distribution data; therefore , we posit the output image from this model as the target morphology for the 3D phase field simulation .', 'As shown earlier in 2D and here in Fig 7 in 3D , the phase field model converges to the target 3D morphology from an initial posited surface , while satisfying the volume and excess surface area constraints .', 'The model does so by iterating the spontaneous curvature function until all constraints are satisfied; once converged , the model then yields the pressure and stresses within the cell surface and cortex that are self-consistent with the 3D surface morphology and constitutive properties of the exterior and cell compartments .', 'In Fig 7A , the target cell morphology is shown .', 'The evolution of cell morphology , from an initial rounded cell guess to the target cell shape , is provided in ( Fig 7B and 7C ) .', 'Fig 7D depicts 2D projections in three mutually orthogonal planes of the cell surface morphology as well as the cortical layer and interior cytosol domains , displaying the values of the phase field variables for each domain .', 'In Fig 7E and 7F the model predictions for hydrostatic pressure distributions and stored stress in the same orthogonal planar sections for the stationary morphology are given correspondingly .', 'The pressure values are not unreasonable ( e . g . a 1 mm depth of water at atmospheric pressure exerts a hydrostatic pressure of 9 . 8 Pa ) .', 'The pressure is low and positive in the external aqueous medium and cytosol; therefore , an inward pressure is exerted from the external medium to the cell surface and an outward pressure from the cell interior ( cytoplasm ) to the cortex .', 'The pressure is negative in the plasma membrane and cortical layer meaning that this layer experiences an inward compressive pressure from the external aqueous medium and cytosol .', 'The highest ( compressive ) pressures arise in the cell surface “interphase” nearby high curvature BLiPs , with about an order of magnitude lower values within the cortical layer itself .', 'These stationary pressure gradients suggest a propensity for fluid absorption from the exterior aqueous medium into the cortex phase and cell surface interface .', 'It is important to note that our simulations assume a stationary morphology , and the pressure and stress distributions are a consequence of the stationary assumption .', 'In reality , these morphologies are non-stationary , and in particular there is flow of the cytoplasm that fills the BLiPs .', 'The new balance of pressure and stress will then dictate the directional flux and flow of the external aqueous medium and cytosol through the various cell compartments .', 'The 3D orientational distributions of cortical actin-filaments , comprising the nematic cortical phase , are given in Fig 8A; planar 2D slices are shown at the specified positions in Fig 8B–8D .', 'These figures convey the degree and direction of order within the F-actin filaments of the cortex , and their strong correlation with the cell surface morphology .', 'Note that cortical F-actin is assumed to prefer parallel alignment at both cortical interfaces for this illustration of the 3D phase field model .', 'The anchoring energy at the interfaces together with the presumed strength of the nematic potential are responsible for the relatively high degree of alignment of the F-actin; these parameters are tunable to match experimental data on nematic order within the cortex .', 'The nematic order was changed by varying Flory order parameter from 1 to values approaching 0 .', 'The results of these simulations are shown in S6 Fig which shows the nematic order superimposed on the 3D morphology for order parameters that range from 1 . 0 to 0 . 01 .', 'Cross-sections taken at three orthogonal planes are shown in S7 Fig . Note that the values assumed for K ( the Franck elastic constant ) , and h1 and h2 will not significantly affect the pressure or the Ericksen stress of the stationary morphology; this is because the dominant contribution to the Ericksen stress is from spatial gradients of the level set function ϕ1 = 1/2 that defines the cell surface .', 'As shown in S6 and S7 Figs , h1 and h2 only affect the nematic order of the stationary state and K does not affect any of the predictions for the stationary morphology .', 'However , for dynamic processes , these parameter values will strongly dictate the results of the simulation .'], ['Using comparative measurements of many individual cells in two distinct configurations , spread on a substrate versus in a rounded state detached from the substrate , we showed that the cell surface area in the rounded state is highly convoluted and far exceeds the surface area of a sphere that would enclose the volume of the cortex , cytosol and nucleus .', 'We analyzed the size and distribution of bleb-like protrusions ( BLiPs ) on the cell periphery that served as storage for the excess of surface area .', 'We then developed three complementary modeling approaches that incorporate the concept of excess surface area on rounded cells in different ways .', 'Employing a 2D discrete geometric model , we tested whether the two-layer composition of the cortex , with the outer layer termed the cell surface , giving the local shape of the BLiPs and the inner layer responsible for contraction , is sufficient to reproduce the highly folded surface observed after rounding .', 'This model demonstrated that morphologies with longer perimeters corresponding to larger ERs will have lower energy compared to the shapes with the same number of folds but with shorter perimeters ( i . e . smaller ERs ) .', 'This result predicts that during cell rounding , larger folds , which are energetically more favorable , should appear early in the process .', 'On the other hand , smaller folds will appear later because it can require time to build structures that will support the higher curvature folds .', 'Indeed , several preliminary experiments in which cells are detached and imaged appear to support this notion .', 'In the early stages of rounding , big folds often appear on the cell surface but , as time goes on , the cell breaks the large folds into the smaller ones .', 'This is supported by the fact that we see a prevalence of small BLiPs in the distribution of fold sizes from the SEM and fluorescence imaging .', 'It is also possible that smaller BLiPs are required as the cell approaches a rounded , steady state because the large folds stored too m
['Introduction', 'Results', 'Discussion', 'Methods']
['Cells transition from spread to rounded morphologies in diverse physiological contexts including mitosis and mesenchymal-to-amoeboid transitions .', 'When these drastic shape changes occur rapidly , cell volume and surface area are approximately conserved .', 'Consequently , the rounded cells are suddenly presented with a several-fold excess of cell surface whose area far exceeds that of a smooth sphere enclosing the cell volume .', 'This excess is stored in a population of bleb-like protrusions ( BLiPs ) , whose size distribution is shown by electron micrographs to be skewed .', 'We introduce three complementary models of rounded cell morphologies with a prescribed excess surface area .', 'A 2D Hamiltonian model provides a mechanistic description of how discrete attachment points between the cell surface and cortex together with surface bending energy can generate a morphology that satisfies a prescribed excess area and BLiP number density .', 'A 3D random seed-and-growth model simulates efficient packing of BLiPs over a primary rounded shape , demonstrating a pathway for skewed BLiP size distributions that recapitulate 3D morphologies .', 'Finally , a phase field model ( 2D and 3D ) posits energy-based constitutive laws for the cell membrane , nematic F-actin cortex , interior cytosol , and external aqueous medium .', 'The cell surface is equipped with a spontaneous curvature function , a proxy for the cell surface-cortex couple , that is a priori unknown , which the model “learns” from the thin section transmission electron micrograph image ( 2D ) or the “seed and growth” model image ( 3D ) .', 'Converged phase field simulations predict self-consistent amplitudes and spatial localization of pressure and stress throughout the cell for any posited stationary morphology target and cell compartment constitutive properties .', 'The models form a general framework for future studies of cell morphological dynamics in a variety of biological contexts .']
['Individual cells must have the capability for rapid morphological transformations under various physiological conditions .', 'One of the most drastic shape transformations occurs during the transition from spread to rounded morphologies .', 'When this transition occurs rapidly , there is insufficient time for significant changes in surface area to occur , although the final size of the rounded cell indicates a significant reduction in apparent cell surface area at light microscope resolution .', 'By contrast , high-resolution scanning electron micrographs of rapidly rounded cells reveal that a large amount of surface area is stored in a highly convoluted surface morphology consisting of bleb-like protrusions ( BLiPs ) and other small structures that are unrecognizable at lower resolution .', 'This surface reserve is an important part of the mechanism that allows rapid and efficient large scale transformations of cell shape .', 'Remarkably , although this convoluted morphology has been observed for decades , there has been very little effort recognizing and including this surface surplus in mathematical modeling of cell morphology and physiology .', 'In this paper , we develop three complementary models to fill this void and lay the foundation for future investigations of the mechanisms that drive cellular morphological dynamics .']
['crystallization techniques', 'perimeters', 'biological cultures', 'radii', 'geometry', 'interface diffusion', 'mathematics', 'cellular structures and organelles', 'thermodynamics', 'research and analysis methods', 'separation processes', 'cho cells', 'cell lines', 'cell membranes', 'free energy', 'physics', 'cell biology', 'curvature', 'biology and life sciences', 'physical sciences', 'cytosol']