Morphological Transformation and Force Generation of Active Cytoskeletal Networks

Emergence of diverse patterns driven by molecular motors in the motility assay

Actomyosin contractility originating from interactions between F-actin and myosin motors in the actin cytoskeleton generates mechanical forces and drives a wide range of cellular processes including cell migration and cytokinesis. To probe the interactions between F-actin and myosin motors, the myosin motility assay has been popularly employed, which consists of myosin heads attached to a glass surface and F-actins gliding on the surface via interactions with the heads. Several experiments have shown that F-actins move in a collective fashion due to volume-exclusion effects between neighboring F-actins. Furthermore, Computational models have shown how changes in key parameters lead to diverse pattern formation in motility assay. However, in most of the computational models, myosin motors were implicitly considered by applying a constant propulsion force to filaments to reduce computational cost. This simplification limits the physiological relevance of the insights provided by the models and potentially leads to artifacts. In this study, we employed an agent-based computational model for the motility assay with explicit immobile motors interacting with filaments. We rigorously account for the kinetics of myosin motors including the force-velocity relationship for walking and the binding and unbinding behaviors. We probed the effects of the length, rigidity, and concentration of filaments and repulsive strength on collective movements and pattern formation. It was found that four distinct types of structures-homogeneous networks, flocks, bands, and rings-emerged as a result of collisions between gliding filaments. We further analyzed the frequency and morphology of these structures and the curvature, alignment, and rotational motions of filaments. Our study provides better insights into the origin and properties of patterns formed by gliding filaments beyond what was shown before.

Myosin-induced F-actin fragmentation facilitates contraction of actin networks

Mechanical forces play a crucial role in diverse physiological processes, such as cell migration, cytokinesis, and morphogenesis. The actin cytoskeleton generates a large fraction of the mechanical forces via molecular interactions between actin filaments (F-actins) and myosin motors. Recent studies have shown that the common tendency of actomyosin networks to contract into a smaller structure deeply involves F-actin buckling induced by motor activities, fragmentation of F-actins, and the force-dependent unbinding of cross-linkers that inter-connect F-actins. The fragmentation of F-actins was shown to originate from either buckling or tensile force from previous single-molecule experiments. While the role of buckling in network contraction has been studied extensively, to date, the role of tension-induced F-actin fragmentation in network contraction has not been investigated. In this study, we employed in vitro experiments and an agent-based computational model to illuminate when and how the tension-induced F-actin fragmentation facilitates network contraction. Our experiments demonstrated that F-actins can be fragmented due to tensile forces, immediately followed by catastrophic rupture and contraction of networks. Using the agent-based model, we showed that F-actin fragmentation by tension results in distinct rupture dynamics different from that observed in networks only with cross-linker unbinding. Moreover, we found that tension-induced F-actin fragmentation is particularly important for the contraction of networks with high connectivity. Results from our study shed light on an important regulator of the contraction of actomyosin networks which has been neglected. In addition, our results provide insights into the rupture mechanisms of polymeric network structures and bio-inspired materials.

E-cadherin adhesion dynamics as revealed by an accelerated force ramp are dependent upon the presence of α-catenin

Tissue remodeling and shape changes often rely on force-induced cell rearrangements occurring via cell-cell contact dynamics. Epithelial cell-cell contact shape changes are particularly dependent upon E-cadherin adhesion dynamics which are directly influenced by cell-generated and external forces. While both the mobility of E-cadherin adhesions and their adhesion strength have been reported before, it is not clear how these two aspects of E-cadherin adhesion dynamics are related. Here, using magnetic pulling cytometry, we applied an accelerated force ramp on the E-cadherin adhesion between an E-cadherin-coated magnetic microbead and an epithelial cell to ascertain this relationship. Our approach enables the determination of the adhesion strength and force-dependent mobility of individual adhesions, which revealed a direct correlation between these key characteristics. Since α-catenin has previously been reported to play a role in both E-cadherin mobility and adhesion strength when studied independently, we also probed epithelial cells in which α-catenin has been knocked out. We found that, in the absence of α-catenin, E-cadherin adhesions not only had lower adhesion strength, as expected, but were also more mobile. We observed that α-catenin was required for the recovery of strained cell-cell contacts and propose that the adhesion strength and force-dependent mobility of E-cadherin adhesions act in tandem to regulate cell-cell contact homeostasis. Our approach introduces a method which relates the force-dependent adhesion mobility to adhesion strength and highlights the morphological role played by α-catenin in E-cadherin adhesion dynamics.

Contractile ring composition dictates kinetics of in silico contractility

Constriction kinetics of the cytokinetic ring are expected to depend on dynamic adjustment of contractile ring composition, but the impact of ring component abundance dynamics on ring constriction is understudied. Computational models generally assume that contractile networks maintain constant total amounts of components, which is not always true. To test how compositional dynamics affect constriction kinetics, we first measured F-actin, non-muscle myosin II, septin, and anillin during Caenorhabditis elegans zygotic mitosis. A custom microfluidic device that positioned the cell with the division plane parallel to a light sheet allowed even illumination of the cytokinetic ring. Measured component abundances were implemented in a three-dimensional agent-based model of a membrane-associated contractile ring. With constant network component amounts, constriction completed with biologically unrealistic kinetics. However, imposing the measured changes in component quantities allowed this model to elicit realistic constriction kinetics. Simulated networks were more sensitive to changes in motor and filament amounts than those of crosslinkers and tethers. Our findings highlight the importance of network composition for actomyosin contraction kinetics.

α-Catenin Dependent E-cadherin Adhesion Dynamics as Revealed by an Accelerated Force Ramp

Tissue remodeling and shape changes often rely on force-induced cell rearrangements occurring via cell-cell contact dynamics. Epithelial cell-cell contact shape changes are particularly dependent upon E-cadherin adhesion dynamics which are directly influenced by cell-generated and external forces. While both the mobility of E-cadherin adhesions and their adhesion strength have been reported before, it is not clear how these two aspects of E-cadherin adhesion dynamics are related. Here, using magnetic pulling cytometry, we applied an accelerated force ramp on the E-cadherin adhesion between an E-cadherin-coated magnetic microbead and an epithelial cell to ascertain this relationship. Our approach enables the determination of the adhesion strength and force-dependent mobility of individual adhesions, which revealed a direct correlation between these key characteristics. Since α-catenin has previously been reported to play a role in both E-cadherin mobility and adhesion strength when studied independently, we also probed epithelial cells in which α-catenin has been knocked out. We found that, in the absence of α-catenin, E-cadherin adhesions not only had lower adhesion strength, as expected, but were also more mobile. We observed that α-catenin was required for the recovery of strained cell-cell contacts and propose that the adhesion strength and force-dependent mobility of E-cadherin adhesions act in tandem to regulate cell-cell contact homeostasis. Our approach introduces a method which relates the force-dependent adhesion mobility to adhesion strength and highlights the morphological role played by α-catenin in E-cadherin adhesion dynamics.

Computational model of integrin adhesion elongation under an actin fiber

Cells create physical connections with the extracellular environment through adhesions. Nascent adhesions form at the leading edge of migrating cells and either undergo cycles of disassembly and reassembly, or elongate and stabilize at the end of actin fibers. How adhesions assemble has been addressed in several studies, but the exact role of actin fibers in the elongation and stabilization of nascent adhesions remains largely elusive. To address this question, here we extended our computational model of adhesion assembly by incorporating an actin fiber that locally promotes integrin activation. The model revealed that an actin fiber promotes adhesion stabilization and elongation. Actomyosin contractility from the fiber also promotes adhesion stabilization and elongation, by strengthening integrin-ligand interactions, but only up to a force threshold. Above this force threshold, most integrin-ligand bonds fail, and the adhesion disassembles. In the absence of contraction, actin fibers still support adhesions stabilization. Collectively, our results provide a picture in which myosin activity is dispensable for adhesion stabilization and elongation under an actin fiber, offering a framework for interpreting several previous experimental observations.

Magnetic microrheometry of tumor-relevant stiffness levels and probabilistic quantification of viscoelasticity differences inside 3D cell culture matrices

The progression of breast cancer involves cancer-cell invasions of extracellular matrices. To investigate the progression, 3D cell cultures are widely used along with different types of matrices. Currently, the matrices are often characterized using parallel-plate rheometry for matrix viscoelasticity, or liquid-like viscous and stiffness-related elastic characteristics. The characterization reveals averaged information and sample-to-sample variation, yet, it neglects internal heterogeneity within matrices, experienced by cancer cells in 3D culture. Techniques using optical tweezers and magnetic microrheometry have measured heterogeneity in viscoelasticity in 3D culture. However, there is a lack of probabilistic heterogeneity quantification and cell-size-relevant, microscale-viscoelasticity measurements at breast-tumor tissue stiffness up to ≃10 kPa in Young's modulus. Here, we have advanced methods, for the purpose, which use a magnetic microrheometer that applies forces on magnetic spheres within matrices, and detects the spheres displacements. We present probabilistic heterogeneity quantification using microscale-viscoelasticity measurements in 3D culture matrices at breast-tumor-relevant stiffness levels. Bayesian multilevel modeling was employed to distinguish heterogeneity in viscoelasticity from the effects of experimental design and measurement errors. We report about the heterogeneity of breast-tumor-relevant agarose, GrowDex, GrowDex-collagen and fibrin matrices. The degree of heterogeneity differs for stiffness, and phase angle (i.e. ratio between viscous and elastic characteristics). Concerning stiffness, agarose and GrowDex show the lowest and highest heterogeneity, respectively. Concerning phase angle, fibrin and GrowDex-collagen present the lowest and the highest heterogeneity, respectively. While this heterogeneity information involves softer matrices, probed by ≃30 μm magnetic spheres, we employ larger ≃100 μm spheres to increase magnetic forces and acquire a sufficient displacement signal-to-noise ratio in stiffer matrices. Thus, we show pointwise microscale viscoelasticity measurements within agarose matrices up to Young's moduli of 10 kPa. These results establish methods that combine magnetic microrheometry and Bayesian multilevel modeling for enhanced heterogeneity analysis within 3D culture matrices.

Nanoscale Tracking Combined with Cell-Scale Microrheology Reveals Stepwise Increases in Force Generated by Cancer Cell Protrusions

In early breast cancer progression, cancer cells invade through a nanoporous basement membrane (BM) as a first key step toward metastasis. This invasion is thought to be mediated by a combination of proteases, which biochemically degrade BM matrix, and physical forces, which mechanically open up holes in the matrix. To date, techniques that quantify cellular forces of BM invasion in 3D culture have been unavailable. Here, we developed cellular-force measurements for breast cancer cell invasion in 3D culture that combine multiple-particle tracking of force-induced BM-matrix displacements at the nanoscale, and magnetic microrheometry of localized matrix mechanics. We find that cancer-cell protrusions exert forces from picoNewtons up to nanoNewtons during invasion. Strikingly, the protrusions extension involves stepwise increases in force, in steps of 0.2 to 0.5 nN exerted from every 30 s to 6 min. Thus, this technique reveals previously unreported dynamics of force generation by invasive protrusions in cancer cells.

F-actin bending facilitates net actomyosin contraction By inhibiting expansion with plus-end-located myosin motors

Contraction of actomyosin networks underpins important cellular processes including motility and division. The mechanical origin of actomyosin contraction is not fully-understood. We investigate whether contraction arises on the scale of individual filaments, without needing to invoke network-scale interactions. We derive discrete force-balance and continuum partial differential equations for two symmetric, semi-flexible actin filaments with an attached myosin motor. Assuming the system exists within a homogeneous background material, our method enables computation of the stress tensor, providing a measure of contractility. After deriving the model, we use a combination of asymptotic analysis and numerical solutions to show how F-actin bending facilitates contraction on the scale of two filaments. Rigid filaments exhibit polarity-reversal symmetry as the motor travels from the minus to plus-ends, such that contractile and expansive components cancel. Filament bending induces a geometric asymmetry that brings the filaments closer to parallel as a myosin motor approaches their plus-ends, decreasing the effective spring force opposing motor motion. The reduced spring force enables the motor to move faster close to filament plus-ends, which reduces expansive stress and gives rise to net contraction. Bending-induced geometric asymmetry provides both new understanding of actomyosin contraction mechanics, and a hypothesis that can be tested in experiments.

Rapid assembly of a polar network architecture drives efficient actomyosin contractility

Actin network architecture and dynamics play a central role in cell contractility and tissue morphogenesis. RhoA-driven pulsed contractions are a generic mode of actomyosin contractility, but the mechanisms underlying how their specific architecture emerges and how this architecture supports the contractile function of the network remain unclear. Here we show that, during pulsed contractions, the actin network is assembled by two subpopulations of formins: a functionally inactive population (recruited) and formins actively participating in actin filament elongation (elongating). We then show that elongating formins assemble a polar actin network, with barbed ends pointing out of the pulse. Numerical simulations demonstrate that this geometry favors rapid network contraction. Our results show that formins convert a local RhoA activity gradient into a polar network architecture, causing efficient network contractility, underlying the key function of kinetic controls in the assembly and mechanics of cortical network architectures.

RhoA-driven actomyosin contractility plays a key role in driving cell and tissue contractility during morphogenesis. Tracking individual formins, Costache et al. show that the network assembled downstream of RhoA displays a polar architecture, barbed ends pointing outward, a feature that supports efficient contractility and force transmission during pulsed contractions.

Interplay between Brownian motion and cross-linking kinetics controls bundling dynamics in actin networks

Morphology changes in cross-linked actin networks are important in cell motility, division, and cargo transport. Here, we study the transition from a weakly cross-linked network of actin filaments to a heavily cross-linked network of actin bundles through microscopic Brownian dynamics simulations. We show that this transition occurs in two stages: first, a composite bundle network of small and highly aligned bundles evolves from cross-linking of individual filaments and, second, small bundles coalesce into the clustered bundle state. We demonstrate that Brownian motion speeds up the first stage of this process at a faster rate than the second. We quantify the time to reach the composite bundle state and show that it strongly increases as the mesh size increases only when the concentration of cross-links is small and that it remains roughly constant if we decrease the relative ratio of cross-linkers as we increase the actin concentration. Finally, we examine the dependence of the bundling timescale on filament length, finding that shorter filaments bundle faster because they diffuse faster.

Plastin and spectrin cooperate to stabilize the actomyosin cortex during cytokinesis

Cytokinesis, the process that partitions the mother cell into two daughter cells, requires the assembly and constriction of an equatorial actomyosin network. Different types of non-motor F-actin crosslinkers localize to the network, but their functional contribution remains poorly understood. Here, we describe a synergy between the small rigid crosslinker plastin and the large flexible crosslinker spectrin in the C. elegans one-cell embryo. In contrast to single inhibitions, co-inhibition of plastin and the βH-spectrin (SMA-1) results in cytokinesis failure due to progressive disorganization and eventual collapse of the equatorial actomyosin network. Cortical localization dynamics of non-muscle myosin II in co-inhibited embryos mimic those observed after drug-induced F-actin depolymerization, suggesting that the combined action of plastin and spectrin stabilizes F-actin in the contractile ring. An in silico model predicts that spectrin is more efficient than plastin at stabilizing the ring and that ring formation is relatively insensitive to βH-spectrin length, which is confirmed in vivo with a sma-1 mutant that lacks 11 of its 29 spectrin repeats. Our findings provide the first evidence that spectrin contributes to cytokinesis and highlight the importance of crosslinker interplay for actomyosin network integrity.

Discrete mechanical model of lamellipodial actin network implements molecular clutch mechanism and generates arcs and microspikes

Mechanical forces, actin filament turnover, and adhesion to the extracellular environment regulate lamellipodial protrusions. Computational and mathematical models at the continuum level have been used to investigate the molecular clutch mechanism, calculating the stress profile through the lamellipodium and around focal adhesions. However, the forces and deformations of individual actin filaments have not been considered while interactions between actin networks and actin bundles is not easily accounted with such methods. We develop a filament-level model of a lamellipodial actin network undergoing retrograde flow using 3D Brownian dynamics. Retrograde flow is promoted in simulations by pushing forces from the leading edge (due to actin polymerization), pulling forces (due to molecular motors), and opposed by viscous drag in cytoplasm and focal adhesions. Simulated networks have densities similar to measurements in prior electron micrographs. Connectivity between individual actin segments is maintained by permanent and dynamic crosslinkers. Remodeling of the network occurs via the addition of single actin filaments near the leading edge and via filament bond severing. We investigated how several parameters affect the stress distribution, network deformation and retrograde flow speed. The model captures the decrease in retrograde flow upon increase of focal adhesion strength. The stress profile changes from compression to extension across the leading edge, with regions of filament bending around focal adhesions. The model reproduces the observed reduction in retrograde flow speed upon exposure to cytochalasin D, which halts actin polymerization. Changes in crosslinker concentration and dynamics, as well as in the orientation pattern of newly added filaments demonstrate the model's ability to generate bundles of filaments perpendicular (actin arcs) or parallel (microspikes) to the protruding direction.

Filament-guided filament assembly provides structural memory of filament alignment during cytokinesis

During cytokinesis, animal cells rapidly remodel the equatorial cortex to build an aligned array of actin filaments called the contractile ring. Local reorientation of filaments by active equatorial compression is thought to underlie the emergence of filament alignment during ring assembly. Here, combining single molecule analysis and modeling in one-cell C. elegans embryos, we show that filaments turnover is far too fast for reorientation of individual filaments by equatorial compression to explain the observed alignment, even if favorably oriented filaments are selectively stabilized. By tracking single formin/CYK-1::GFP particles to monitor local filament assembly, we identify a mechanism that we call filament-guided filament assembly (FGFA), in which existing filaments serve as templates to orient the growth of new filaments. FGFA sharply increases the effective lifetime of filament orientation, providing structural memory that allows cells to build highly aligned filament arrays in response to equatorial compression, despite rapid turnover of individual filaments.

Protein friction and filament bending facilitate contraction of disordered actomyosin networks

We use mathematical modeling and computation to investigate how protein friction facilitates contraction of disordered actomyosin networks. We simulate two-dimensional networks using an agent-based model, consisting of a system of force-balance equations for myosin motor proteins and semiflexible actin filaments. A major advantage of our approach is that it enables direct calculation of the network stress tensor, which provides a quantitative measure of contractility. Exploiting this, we use repeated simulations of disordered networks to confirm that both protein friction and actin filament bending are required for contraction. We then use simulations of elementary two-filament systems to show that filament bending flexibility can facilitate contraction on the microscopic scale. Finally, we show that actin filament turnover is necessary to sustain contraction and prevent filament aggregation. Simulations with and without turnover also exhibit contractile pulses. However, these pulses are aperiodic, suggesting that periodic pulsation can only arise because of additional regulatory mechanisms or more complex mechanical behavior.

Tools for computational analysis of moving boundary problems in cellular mechanobiology

A cell's ability to change shape is one of the most fundamental biological processes and is essential for maintaining healthy organisms. When the ability to control shape goes awry, it often results in a diseased system. As such, it is important to understand the mechanisms that allow a cell to sense and respond to its environment so as to maintain cellular shape homeostasis. Because of the inherent complexity of the system, computational models that are based on sound theoretical understanding of the biochemistry and biomechanics and that use experimentally measured parameters are an essential tool. These models involve an inherent feedback, whereby shape is determined by the action of regulatory signals whose spatial distribution depends on the shape. To carry out computational simulations of these moving boundary problems requires special computational techniques. A variety of alternative approaches, depending on the type and scale of question being asked, have been used to simulate various biological processes, including cell motility, division, mechanosensation, and cell engulfment. In general, these models consider the forces that act on the system (both internally generated, or externally imposed) and the mechanical properties of the cell that resist these forces. Moving forward, making these techniques more accessible to the non-expert will help improve interdisciplinary research thereby providing new insight into important biological processes that affect human health. This article is categorized under: Cancer > Cancer>Computational Models Cancer > Cancer>Molecular and Cellular Physiology.

Modeling multicellular dynamics regulated by extracellular-matrix-mediated mechanical communication via active particles with polarized effective attraction

Collective cell migration is crucial to many physiological and pathological processes such as embryo development, wound healing, and cancer invasion. Recent experimental studies have indicated that the active traction forces generated by migrating cells in a fibrous extracellular matrix (ECM) can mechanically remodel the ECM, giving rise to bundlelike mesostructures bridging individual cells. Such fiber bundles also enable long-range propagation of cellular forces, leading to correlated migration dynamics regulated by the mechanical communication among the cells. Motivated by these experimental discoveries, we develop an active-particle model with polarized effective attractions (APPA) to investigate emergent multicellular migration dynamics resulting from ECM-mediated mechanical communications. In particular, the APPA model generalizes the classic active-Brownian-particle (ABP) model by imposing a pairwise polarized attractive force between the particles, which depends on the instantaneous dynamic states of the particles and mimics the effective mutual pulling between the cells via the fiber bundle bridge. The APPA system exhibits enhanced aggregation behaviors compared to the classic ABP system, and the contrast is more apparent at lower particle densities and higher rotational diffusivities. Importantly, in contrast to the classic ABP system where the particle velocities are not correlated for all particle densities, the high-density phase of the APPA system exhibits strong dynamic correlations, which are characterized by the slowly decaying velocity correlation functions with a correlation length comparable to the linear size of the high-density phase domain (i.e., the cluster of particles). The strongly correlated multicellular dynamics predicted by the APPA model is subsequently verified in in vitro experiments using MCF-10A cells. Our studies indicate the importance of incorporating ECM-mediated mechanical coupling among the migrating cells for appropriately modeling emergent multicellular dynamics in complex microenvironments.

Tensile force-induced cytoskeletal remodeling: Mechanics before chemistry

Understanding cellular remodeling in response to mechanical stimuli is a critical step in elucidating mechanical activation of biochemical signaling pathways. Experimental evidence indicates that external stress-induced subcellular adaptation is accomplished through dynamic cytoskeletal reorganization. To study the interactions between subcellular structures involved in transducing mechanical signals, we combined experimental data and computational simulations to evaluate real-time mechanical adaptation of the actin cytoskeletal network. Actin cytoskeleton was imaged at the same time as an external tensile force was applied to live vascular smooth muscle cells using a fibronectin-functionalized atomic force microscope probe. Moreover, we performed computational simulations of active cytoskeletal networks under an external tensile force. The experimental data and simulation results suggest that mechanical structural adaptation occurs before chemical adaptation during filament bundle formation: actin filaments first align in the direction of the external force by initializing anisotropic filament orientations, then the chemical evolution of the network follows the anisotropic structures to further develop the bundle-like geometry. Our findings present an alternative two-step explanation for the formation of actin bundles due to mechanical stimulation and provide new insights into the mechanism of mechanotransduction.

Remodeling the cytoskeletal network in response to external force is key to cellular mechanotransduction. Despite much focus on cytoskeletal remodeling in recent years, a comprehensive understanding of actin remodeling in real-time in cells under mechanical stimuli is still lacking. We integrated tensile stress-induced 3D actin remodeling and 3D computational simulations of actin cytoskeleton to study how the actin cytoskeleton form bundles and how these bundles evolve over time upon external tensile stress. We found that actin network remodels through a two-step process in which rapid alignment of actin filaments is followed by slower actin bundling. Based on these results, we propose a “mechanics before chemistry” model of actin cytoskeleton remodeling under external tensile force.

Bond Type and Discretization of Nonmuscle Myosin II Are Critical for Simulated Contractile Dynamics

Molecular motors drive cytoskeletal rearrangements to change cell shape. Myosins are the motors that move, cross-link, and modify the actin cytoskeleton. The primary force generator in contractile actomyosin networks is nonmuscle myosin II (NMMII), a molecular motor that assembles into ensembles that bind, slide, and cross-link actin filaments (F-actin). The multivalence of NMMII ensembles and their multiple roles have confounded the resolution of crucial questions, including how the number of NMMII subunits affects dynamics and what affects the relative contribution of ensembles' cross-linking versus motoring activities. Because biophysical measurements of ensembles are sparse, modeling of actomyosin networks has aided in discovering the complex behaviors of NMMII ensembles. Myosin ensembles have been modeled via several strategies with variable discretization or coarse graining and unbinding dynamics, and although general assumptions that simplify motor ensembles result in global contractile behaviors, it remains unclear which strategies most accurately depict cellular activity. Here, we used an agent-based platform, Cytosim, to implement several models of NMMII ensembles. Comparing the effects of bond type, we found that ensembles of catch-slip and catch motors were the best force generators and binders of filaments. Slip motor ensembles were capable of generating force but unbound frequently, resulting in slower contractile rates of contractile networks. Coarse graining of these ensemble types from two sets of 16 motors on opposite ends of a stiff rod to two binders, each representing 16 motors, reduced force generation, contractility, and the total connectivity of filament networks for all ensemble types. A parallel cluster model, previously used to describe ensemble dynamics via statistical mechanics, allowed better contractility with coarse graining, though connectivity was still markedly reduced for this ensemble type with coarse graining. Together, our results reveal substantial tradeoffs associated with the process of coarse graining NMMII ensembles and highlight the robustness of discretized catch-slip ensembles in modeling actomyosin networks.

Modeling cell migration regulated by cell extracellular-matrix micromechanical coupling

Cell migration in fibrous extracellular matrix (ECM) is crucial to many physiological and pathological processes such as tissue regeneration, immune response, and cancer progression. During migration, individual cells can generate active pulling forces via actomyosin contraction, which are transmitted to the ECM fibers through focal adhesion complexes, remodel the ECM, and eventually propagate to and can be sensed by other cells in the system. The microstructure and physical properties of the ECM can also significantly influence cell migration, e.g., via durotaxis and contact guidance. Here, we develop a computational model for two-dimensional cell migration regulated by cell-ECM micromechanical coupling. Our model explicitly takes into account a variety of cellular-level processes, including focal adhesion formation and disassembly, active traction force generation and cell locomotion due to actin filament contraction, transmission and propagation of tensile forces in the ECM, as well as the resulting ECM remodeling. We validate our model by accurately reproducing single-cell dynamics of MCF-10A breast cancer cells migrating on collagen gels and show that the durotaxis and contact guidance effects naturally arise as a consequence of the cell-ECM micromechanical interactions considered in the model. Moreover, our model predicts strongly correlated multicellular migration dynamics, which are resulted from the ECM-mediated mechanical coupling among the migrating cell and are subsequently verified in in vitro experiments using MCF-10A cells. Our computational model provides a robust tool to investigate emergent collective dynamics of multicellular systems in complex in vivo microenvironment and can be utilized to design in vitro microenvironments to guide collective behaviors and self-organization of cells.

Myosin II Filament Dynamics in Actin Networks Revealed with Interferometric Scattering Microscopy

The plasma membrane and the underlying cytoskeletal cortex constitute active platforms for a variety of cellular processes. Recent work has shown that the remodeling acto-myosin network modifies local membrane organization, but the molecular details are only partly understood because of difficulties with experimentally accessing the relevant time and length scales. Here, we use interferometric scattering microscopy to investigate a minimal acto-myosin network linked to a supported lipid bilayer membrane. Using the magnitude of the interferometric contrast, which is proportional to molecular mass, and fast acquisition rates, we detect and image individual membrane-attached actin filaments diffusing within the acto-myosin network and follow individual myosin II filament dynamics. We quantify myosin II filament dwell times and processivity as functions of ATP concentration, providing experimental evidence for the predicted ensemble behavior of myosin head domains. Our results show how decreasing ATP concentrations lead to both increasing dwell times of individual myosin II filaments and a global change from a remodeling to a contractile state of the acto-myosin network.

Collective and contractile filament motions in the myosin motility assay

Cells require mechanical forces for their physiological functions. The forces are generated mainly from molecular interactions between actin filaments, cross-linking proteins, and myosin motors in the actin cytoskeleton. To better understand the molecular interactions, many studies employed myosin motility assays with actin filaments propelled by myosin heads fixed on a surface. Various interesting behaviors of actin filaments have been observed in the motility assay experiments. Despite the popularity of the motility assays, there were only a few computational models designed for simulating the motility assay systems. Most of the previous models have limitations which precluded full understanding of molecular origins for behaviors of actin filaments. In this study, we used an agent-based computational model based on Brownian dynamics for simulating the motility assay system. Our model rigorously describes the mechanics, dynamics, and interactions of actin filaments, cross-linking proteins, and molecular motors. Using the model, we first investigated how properties of actin filaments and motors affect gliding motions of actin filaments without volume-exclusion effects as a base study. We found that actin filaments can continuously glide at relative fast speed only when they are sufficiently longer than the average spacing between neighboring motors and that the gliding speed of F-actins shows a biphasic dependence on processivity of motors. Then, we showed that volume-exclusion effects between actin filaments can induce diverse collective movements and alignment of actin filaments, thus creating thick bundles and ring-like structures in the absence of cross-linking proteins. Lastly, we demonstrated that cross-linking proteins can lead to distinct contractile behaviors of actin networks depending on the density and kinetics of the cross-linking proteins. Results from our study show the ability of our model to simulate the motility assay system under various conditions and provide insights into understanding of different behaviors of actin filaments.

Classical and Emerging Regulatory Mechanisms of Cytokinesis in Animal Cells

The primary goal of cytokinesis is to produce two daughter cells, each having a full set of chromosomes. To achieve this, cells assemble a dynamic structure between segregated sister chromatids called the contractile ring, which is made up of filamentous actin, myosin-II, and other regulatory proteins. Constriction of the actomyosin ring generates a cleavage furrow that divides the cytoplasm to produce two daughter cells. Decades of research have identified key regulators and underlying molecular mechanisms; however, many fundamental questions remain unanswered and are still being actively investigated. This review summarizes the key findings, computational modeling, and recent advances in understanding of the molecular mechanisms that control the formation of the cleavage furrow and cytokinesis.

Remarkable structural transformations of actin bundles are driven by their initial polarity, motor activity, crosslinking, and filament treadmilling

Bundled actin structures play a key role in maintaining cellular shape, in aiding force transmission to and from extracellular substrates, and in affecting cellular motility. Recent studies have also brought to light new details on stress generation, force transmission and contractility of actin bundles. In this work, we are primarily interested in the question of what determines the stability of actin bundles and what network geometries do unstable bundles eventually transition to. To address this problem, we used the MEDYAN mechano-chemical force field, modeling several micron-long actin bundles in 3D, while accounting for a comprehensive set of chemical, mechanical and transport processes. We developed a hierarchical clustering algorithm for classification of the different long time scale morphologies in our study. Our main finding is that initially unipolar bundles are significantly more stable compared with an apolar initial configuration. Filaments within the latter bundles slide easily with respect to each other due to myosin activity, producing a loose network that can be subsequently severely distorted. At high myosin concentrations, a morphological transition to aster-like geometries was observed. We also investigated how actin treadmilling rates influence bundle dynamics, and found that enhanced treadmilling leads to network fragmentation and disintegration, while this process is opposed by myosin and crosslinking activities. Interestingly, treadmilling bundles with an initial apolar geometry eventually evolve to a whole gamut of network morphologies based on relative positions of filament ends, such as sarcomere-like organization. We found that apolar bundles show a remarkable sensitivity to environmental conditions, which may be important in enabling rapid cytoskeletal structural reorganization and adaptation in response to intracellular and extracellular cues.

Multiscale model of integrin adhesion assembly

The ability of adherent cells to form adhesions is critical to numerous phases of their physiology. The assembly of adhesions is mediated by several types of integrins. These integrins differ in physical properties, including rate of diffusion on the plasma membrane, rapidity of changing conformation from bent to extended, affinity for extracellular matrix ligands, and lifetimes of their ligand-bound states. However, the way in which nanoscale physical properties of integrins ensure proper adhesion assembly remains elusive. We observe experimentally that both β-1 and β-3 integrins localize in nascent adhesions at the cell leading edge. In order to understand how different nanoscale parameters of β-1 and β-3 integrins mediate proper adhesion assembly, we therefore develop a coarse-grained computational model. Results from the model demonstrate that morphology and distribution of nascent adhesions depend on ligand binding affinity and strength of pairwise interactions. Organization of nascent adhesions depends on the relative amounts of integrins with different bond kinetics. Moreover, the model shows that the architecture of an actin filament network does not perturb the total amount of integrin clustering and ligand binding; however, only bundled actin architectures favor adhesion stability and ultimately maturation. Together, our results support the view that cells can finely tune the expression of different integrin types to determine both structural and dynamic properties of adhesions.

Mobility of Molecular Motors Regulates Contractile Behaviors of Actin Networks

Cells generate mechanical forces primarily from interactions between F-actin, cross-linking proteins, myosin motors, and other actin-binding proteins in the cytoskeleton. To understand how molecular interactions between the cytoskeletal elements generate forces, a number of in vitro experiments have been performed but are limited in their ability to accurately reproduce the diversity of motor mobility. In myosin motility assays, myosin heads are fixed on a surface and glide F-actin. By contrast, in reconstituted gels, the motion of both myosin and F-actin is unrestricted. Because only these two extreme conditions have been used, the importance of mobility of motors for network behaviors has remained unclear. In this study, to illuminate the impacts of motor mobility on the contractile behaviors of the actin cytoskeleton, we employed an agent-based computational model based on Brownian dynamics. We find that if motors can bind to only one F-actin like myosin I, networks are most contractile at intermediate mobility. In this case, less motor mobility helps motors stably pull F-actins to generate tensile forces, whereas higher motor mobility allows F-actins to aggregate into larger clustering structures. The optimal intermediate motor mobility depends on the stall force and affinity of motors that are regulated by mechanochemical rates. In addition, we find that the role of motor mobility can vary drastically if motors can bind to a pair of F-actins. A network can exhibit large contraction with high motor mobility because motors bound to antiparallel pairs of F-actins can exert similar forces regardless of their mobility. Results from this study imply that the mobility of molecular motors may critically regulate contractile behaviors of actin networks in cells.

Unite to divide - how models and biological experimentation have come together to reveal mechanisms of cytokinesis

Cytokinesis is the fundamental and ancient cellular process by which one cell physically divides into two. Cytokinesis in animal and fungal cells is achieved by contraction of an actomyosin cytoskeletal ring assembled in the cell cortex, typically at the cell equator. Cytokinesis is essential for the development of fertilized eggs into multicellular organisms and for homeostatic replenishment of cells. Correct execution of cytokinesis is also necessary for genome stability and the evasion of diseases including cancer. Cytokinesis has fascinated scientists for well over a century, but its speed and dynamics make experiments challenging to perform and interpret. The presence of redundant mechanisms is also a challenge to understand cytokinesis, leaving many fundamental questions unresolved. For example, how does a disordered cytoskeletal network transform into a coherent ring? What are the long-distance effects of localized contractility? Here, we provide a general introduction to 'modeling for biologists', and review how agent-based modeling and continuum mechanics modeling have helped to address these questions.

Nonequilibrium phase diagrams for actomyosin networks

Living cells dynamically modulate the local morphologies of their actin networks to perform biological functions, including force transduction, intracellular transport, and cell division. A major challenge is to understand how diverse structures of the actin cytoskeleton are assembled from a limited set of molecular building blocks. Here we study the spontaneous self-assembly of a minimal model of cytoskeletal materials, consisting of semiflexible actin filaments, crosslinkers, and molecular motors. Using coarse-grained simulations, we demonstrate that by changing concentrations and kinetics of crosslinkers and motors, as well as filament lengths, we can generate three distinct structural phases of actomyosin assemblies: bundled, polarity-sorted, and contracted. We introduce new metrics to distinguish these structural phases and demonstrate their functional roles. We find that the binding kinetics of motors and crosslinkers can be tuned to optimize contractile force generation, motor transport, and mechanical response. By quantitatively characterizing the relationships between the modes of cytoskeletal self-assembly, the resulting structures, and their functional consequences, our work suggests new principles for the design of active materials.

Cell-cell adhesion interface: orthogonal and parallel forces from contraction, protrusion, and retraction

The epithelial lateral membrane plays a central role in the integration of intercellular signals and, by doing so, is a principal determinant in the emerging properties of epithelial tissues. Mechanical force, when applied to the lateral cell-cell interface, can modulate the strength of adhesion and influence intercellular dynamics. Yet the relationship between mechanical force and epithelial cell behavior is complex and not completely understood. This commentary aims to provide an investigative look at the usage of cellular forces at the epithelial cell-cell adhesion interface.

Cross-linkers both drive and brake cytoskeletal remodeling and furrowing in cytokinesis

Cell shape changes such as cytokinesis are driven by the actomyosin contractile cytoskeleton. The molecular rearrangements that bring about contractility in nonmuscle cells are currently debated. Specifically, both filament sliding by myosin motors, as well as cytoskeletal cross-linking by myosins and nonmotor cross-linkers, are thought to promote contractility. Here we examined how the abundance of motor and nonmotor cross-linkers affects the speed of cytokinetic furrowing. We built a minimal model to simulate contractile dynamics in the Caenorhabditis elegans zygote cytokinetic ring. This model predicted that intermediate levels of nonmotor cross-linkers are ideal for contractility; in vivo, intermediate levels of the scaffold protein anillin allowed maximal contraction speed. Our model also demonstrated a nonlinear relationship between the abundance of motor ensembles and contraction speed. In vivo, thorough depletion of nonmuscle myosin II delayed furrow initiation, slowed F-actin alignment, and reduced maximum contraction speed, but partial depletion allowed faster-than-expected kinetics. Thus, cytokinetic ring closure is promoted by moderate levels of both motor and nonmotor cross-linkers but attenuated by an over-abundance of motor and nonmotor cross-linkers. Together, our findings extend the growing appreciation for the roles of cross-linkers in cytokinesis and reveal that they not only drive but also brake cytoskeletal remodeling.

A theory that predicts behaviors of disordered cytoskeletal networks

Morphogenesis in animal tissues is largely driven by actomyosin networks, through tensions generated by an active contractile process. Although the network components and their properties are known, and networks can be reconstituted in vitro, the requirements for contractility are still poorly understood. Here, we describe a theory that predicts whether an isotropic network will contract, expand, or conserve its dimensions. This analytical theory correctly predicts the behavior of simulated networks, consisting of filaments with varying combinations of connectors, and reveals conditions under which networks of rigid filaments are either contractile or expansile. Our results suggest that pulsatility is an intrinsic behavior of contractile networks if the filaments are not stable but turn over. The theory offers a unifying framework to think about mechanisms of contractions or expansion. It provides the foundation for studying a broad range of processes involving cytoskeletal networks and a basis for designing synthetic networks.

Progress in Integrative Biomaterial Systems to Approach Three-Dimensional Cell Mechanotransduction

Mechanotransduction between cells and the extracellular matrix regulates major cellular functions in physiological and pathological situations. The effect of mechanical cues on biochemical signaling triggered by cell-matrix and cell-cell interactions on model biomimetic surfaces has been extensively investigated by a combination of fabrication, biophysical, and biological methods. To simulate the in vivo physiological microenvironment in vitro, three dimensional (3D) microstructures with tailored bio-functionality have been fabricated on substrates of various materials. However, less attention has been paid to the design of 3D biomaterial systems with geometric variances, such as the possession of precise micro-features and/or bio-sensing elements for probing the mechanical responses of cells to the external microenvironment. Such precisely engineered 3D model experimental platforms pave the way for studying the mechanotransduction of multicellular aggregates under controlled geometric and mechanical parameters. Concurrently with the progress in 3D biomaterial fabrication, cell traction force microscopy (CTFM) developed in the field of cell biophysics has emerged as a highly sensitive technique for probing the mechanical stresses exerted by cells onto the opposing deformable surface. In the current work, we first review the recent advances in the fabrication of 3D micropatterned biomaterials which enable the seamless integration with experimental cell mechanics in a controlled 3D microenvironment. Then, we discuss the role of collective cell-cell interactions in the mechanotransduction of engineered tissue equivalents determined by such integrative biomaterial systems under simulated physiological conditions.