Cytoskeleton

The intracellular transport system is pivotal for cellular function and integrity, facilitated by cytoskeletal motor proteins such as dynein, which traverse along microtubules (MTs). The heterogeneity of the tubulin isotypes composing MTs introduces functional diversity, potentially affecting cytoskeletal motor proteins interactions with the MT. This in silico study investigated the influence of amino acid sequence variations in the C-terminal tails (CTTs) of six different Homo sapiens tubulin isotypes, TUBB2A, TUBB2B, TUBB2C, TUBB3, TUBB4A, and TUBB5, highly expressed in human brain tumors, and assessed the isotypes effect on the binding of motor protein dynein to MT. Among these isotypes, TUBB2A, TUBB2B, and TUBB2C were found to affect conformational motions of the dyneins microtubule-binding domain (MTBD) and stalk domain. The investigation highlighted the novel role of isotype-specific variations in lateral interactions between tubulin protofilaments (PFs) in determining the proximity of the {beta}-CTT of the adjacent PF to the MTBD, potentially affecting dyneins motility and suggesting how changes in isotype expression directly influence dyneins velocity and processivity and contribute to transport defects associated with neurological disorders and cancers. Isolating specific tubulin isotypes experimentally is challenging due to their high sequence similarity and complex interactions with other microtubule-associated proteins. This makes it challenging to distinguish between different tubulin isotypes and their effects, particularly in tissues where multiple isotypes are co-expressed. Additionally, these isotypes are heavily modified in vivo by post-translational modifications, which further complicate the isolation of a single, unmodified tubulin isotype. As a result, computational approaches have been necessary in this study for exploring these effects in a controlled, isotype-specific manner.

In multicellular organisms, various cellular structures exhibit cytoplasmic sharing, where cells remain interconnected. While essential for development and function in contexts such as germ cell formation and insect early embryos, the physical basis of cell volume regulation in these systems remains poorly defined. Germline cysts are formed by interconnected sister cells via intercellular bridges. In mice, germline cysts form during gametogenesis in fetal ovaries and testes. In mouse fetal female cysts, cells with numerous bridges preferentially differentiate into oocytes by selectively increasing their volume, a process that may be mediated through cytoplasmic flow. This volume bias may be influenced by hydrostatic pressure within the cytoplasm. Here, we theoretically investigate how the mechanical properties of cells affect cytoplasmic pressure and volume distribution within interconnected cells. Our soap-bubble model revealed that cells with more bridges exhibit increased volume when they have large cell-cell contact areas, as observed in fetal cysts. We found that incorporating cell cycle (including cell growth and cell division) significantly enhances the likelihood of volume bias in favor of cells with more bridges. These theoretical findings suggest that intrinsic mechanical properties, coupled with cell cycle, establish robust cyst development in fetal female germline cysts. Our findings also provide insights into the volume dynamics observed in adult male germline cysts, which are characterized by smaller cell-cell contact areas. Impact statementA theoretical model demonstrates how mechanical properties and cell cycle dynamics regulate volume distribution in germline cysts, providing a physical basis for oocyte differentiation.

The proper assembly, architecture, and maintenance of microtubule, actin and other cytoskeletal networks require regulation by various polymer binding proteins. Microtubules, rely on both the tubulin building blocks, but also many tubulin- and microtubule binding proteins. TOG domain-containing proteins comprise one family of tubulin-binding proteins that regulate microtubule dynamics. Here we identify two previously uncharacterized TOG domain-containing proteins (TOD-1 and TOD-2) in the nematode, C. elegans. These proteins are unique in that they are members of the XMAP215 family but contain reduced numbers of TOG domains and, in one case, a divergent TOG domain. TOD-1 and TOD-2 are expressed in and contribute to the normal function of sperm. The single TOG domain of TOD-1 and both TOG domains of TOD-2 are predicted to bind free tubulin dimers and not microtubule lattice. Deletion of either tod gene resulted in an increased laying of unfertilized oocytes. Inspection of mutant hermaphrodites revealed a premature onset of sperm migration failure. Together, these findings suggest that C. elegans requires regulation of tubulin dimers and/or microtubules for sperm localization and function. The amoeboid movement of C. elegans sperm has been considered microtubule-independent, our results open a new avenue of research into their unique motility.

Ventricular myosin light chain-1 (MLC1v) is a key structural and function-modulating component of the {beta}-cardiac myosin ({beta}M-II) motor complex. Single-point mutations in MLC1v are linked to severe forms of hypertrophic cardiomyopathy (HCM) and sudden cardiac death (SCD) at a young age. However, the molecular mechanisms underlying the motor dysfunction responsible for HCM phenotype development are not fully understood. Here, we investigated native {beta}M-II motors isolated from septal myectomy sample of an HCM patient, harboring a rare homozygous mutation in MLC1v (A57D). Using a pure population of mutant motors (MUT), and sensitive single-molecule functional analysis approach, we directly assessed the primary functional alterations in {beta}M-II bearing A57D MLC1v mutation. In optical trap single-molecules measurements, the mutant motors displayed increased actomyosin (AM) interaction duration in strongly bound state (ton), corresponding to 3-fold reduced AM detachment rate than wild type myosin (WT). The MUT myosin also generated a shorter powerstroke size ({delta}). Ensemble average analysis of AM interaction events demonstrated that both the first powerstroke ({delta}1) associated with Pi release and the second powerstroke ({delta}2) linked to ADP release were reduced in MUT myosin. Moreover, the increased actomyosin cross-bridge stiffness in the AM.ADP state was observed for MUT compared to WT motors. Consistent with slower AM detachment rate and shorter stroke size, reconstituted human mutant {beta}M-II displayed slower actin filament gliding speed. Alterations in sarcomere-level mechanics included increased Ca2+ sensitivity of force generation and prolonged relaxation time, as predicted by FiberSim modelling. Molecular dynamics simulations indicated that the substitution of alanine by aspartate altered MLC1v interactions with myosin heavy chain (MyHC) and light chain 2 (MLC2v), affecting the curvature and flexibility of the lever arm. Overall, these studies establish the molecular mechanism underlying the primary myosin dysfunction due to A57D MLC1v mutation and further highlight the crucial role of MLC1v-mediated regulation of myosin function. Understanding the precise changes in the mutant myosins biomechanical properties is directly relevant to comprehending the initial triggers for pathological cardiac remodeling in HCM patients and designing tailored therapeutic interventions.

Mammalian spermatogenesis occurs in the seminiferous tubules, which exhibit unique spatiotemporal differentiation patterns known as cellular association patterns. In mice, these patterns can be regarded as one-dimensional wavetrains that consistently propagate inward from both ends, resulting in one or more "sites of reversal." Segmented wavetrain pattern, in which the wave propagation direction spatially switches, was observed in our previous three-species reaction-diffusion model for interspecific species difference in spermatogenic waves (Kawamura et al., 2021). However, the biological mechanisms of the formation of sites of reversal and of this directional bias, as well as the principle of pattern formation, remain unknown. Here, we refined our previous model to match the actual biological spatiotemporal scale and examined its dynamics through extensive numerical simulations. The modified model frequently generated segmented wavetrain patterns, corresponding to the sites of reversal, but without directional bias. We systematically examined possible biological mechanisms for the bias and found that tubule elongation, especially near the rete testis, most effectively accounts for the bias among the tested. Extensive simulations revealed that the segmented pattern is numerically stable, emerges more frequently in longer domains, and shows an exponential segment size distribution with a lower limit for the stably existing segment length. These explorations imply that locally emerged unidirectional wavetrains serve as building blocks to generate the stable segmented wavetrains through their interactions. HighlightsO_LISegmented wavetrains reflect sites of reversal in seminiferous tubules. C_LIO_LISegmented patterns frequently emerge but show no inherent directional bias. C_LIO_LITubule elongation may contribute to inward propagation near the rete testis. C_LIO_LISegmented wavetrains are numerically stable and more frequent in longer domains. C_LIO_LIInteractions of local unidirectional wavetrains generate stable segmented structures. C_LI

Mechanosensing involves proteins detecting mechanical changes in the cytoskeleton or at cell adhesion sites. These interactions initiate signaling cascades that produce biochemical effects such as post-translational modifications or cytoskeletal rearrangements. Filamin is a ubiquitous mechanosensing protein that binds actin filaments and senses pulling forces within the cytoskeleton. Drosophila Filamin (Cheerio) is structurally similar to mammalian Filamin, with roles in egg chamber development, embryo cellularization, and integrity of muscle attachment sites and Z discs in Drosophila indirect flight muscles (IFMs). Here we report a potential novel binding partner of Drosophila Filamins: the death-associated protein kinase Drak that functions as a myosin light chain kinase. We found that Drak biochemically bound to an open mutant of Filamin that resembles the mechanically activated form partially bound to wild type Filamin and did not bind to closed mutant of Filamin. The interaction site was mapped to the intrinsically unfolded C-terminal region of Drak. To study the functional role of Drak-Filamin interaction, we studied two developmental events where Drak has been earlier shown to be expressed and where Filamin also functions: early embryonic cellularization and indirect flight muscle development at pupal stages. We found partial colocalization between Drak-GFP and Filamin-mCherry during the initiation of cellularization furrow, and at the time of myotube attachment site maturation in tendon cells. However, functionally we could not show direct correlation between Filamin and Drak. Our studies reveal interesting new expression patterns of Drak during Drosophila development and provide detailed information about Filamin localization during IFM development.

The Drosophila scaffolding protein Zasp52 is required to maintain structure at the muscle Z-disc which experiences strong forces during contraction. It is alternatively spliced into many isoforms, some of which contain a long intrinsically disordered region (IDR). We show that this region is primarily expressed in the indirect flight muscle (IFM) and is required for maintaining integrity of the Z-disc. Deleting the IDR-encoding exon15e results in flightlessness and structural IFM defects, including sarcomere bending at the Z-disc and an inability to de-contract. These defects are indicative of a lack of proper thin filament anchoring to the Z-disc. This is further supported by a genetic interaction between exon15e and actin. Fluorescence recovery after photobleaching of an isoform lacking exon15e shows that the IDR is required for maintaining Zasp52 at the Z-disc and thereby stabilizing Z-discs. Lastly, we can rescue these phenotypes by restricting IFM use. Together, these results suggest that Zasp52s IDR confers thin filament stability at the Z-disc of IFM.

Motile cells can sense and exert forces on the extracellular environment through dynamic actin networks. Increased stress against the polymerizing barbed ends of branched actin networks has been shown to lead to an increase in the density of these networks through a force feedback mechanism, though this phenomenon has not been explored through the examination of real-time responses of endogenous actin networks in cells. Here, we utilize mouse embryonic fibroblast CRISPR knock-in lines with labeled ARP2/3 complex to identify cellular and extracellular conditions that regulate branched actin density and enrichment at the leading edge of lamellipodial protrusions. A common theme shared among all branched actin density-increasing conditions is higher levels of interface stress between the plasma membrane and the barbed ends of the lamellipodial actin network. Among these conditions, we find that ARP2/3 is specifically required for robust spreading and protrusion in response to increased extracellular viscosity. Interestingly, time-lapse traction force microscopy of ARP2/3-dependent viscosity responses show significantly reduced changes in strain energy applied to the substrate when compared to spreading and motility through cell-matrix adhesion. In addition, we find that increased extracellular viscosity can bypass the need for extracellular matrix proteins to support lamellipodial protrusion driven by optogenetic Rac activation. Our studies provide strong support for in vitro models of branched actin force feedback responses and further characterize an essential role for branched actin in mediating dramatic cell shape changes in response to increased extracellular viscosity.

The actin cytoskeleton is an inherently disordered active system. the actomyosin cortex and reconstituted actomyosin systems are globally disordered, yet undergo transitions between distinct disordered states as parameters like motor and crosslinker concentration and filament length and rigidity change. In cells these changes are related to genetic mutations or differences in cell state and dictate fundamental biological processes. However, we dont have well established methods to detect and classify differences in disordered polymer networks. Image-based morphology techniques provide a non-invasive, high-throughput method of extracting information about a system. In this work we simulate biopolymer networks under varying conditions and develop and use morphological descriptors to construct trajectories in morphospace. Using statistical analysis we find that morphological descriptors are able to distinguish between different trajectories of the system, including differences not apparent to the eye. However, no single descriptor alone is able to capture all the differences in the simulated trajectories. Nematic order parameters typically perform the worst for our simulations while curvature and texture descriptors can collectively distinguish between dynamic trajectories. This work helps develop quantification of cytoskeleton dynamics for classification and data-driven modeling.

Cell contractility plays numerous essential roles in a healthy organism, while its malfunctioning can lead to disease. The ubiquitous actin-dependent motors of the nonmuscle myosin 2 (NM2) family, which includes NM2A and NM2B, are chiefly responsible for cell contractility because of their ability to polymerize into bipolar filaments. Polymerization/depolymerization of NM2 filaments allows cells to quickly reorganize their contractile system according to constantly changing shapes and positions of nonmuscle cells. Bipolar filament depolymerization is known to depend on the C-terminal features of the NM2A heavy chain. Here, we show that the motor activity of NM2A is another key component of NM2As depolymerization mechanism, which cooperates with tail-dependent mechanisms to facilitate NM2A turnover in cells and, through copolymerization with NM2B, to reorganize and dynamize NM2B in trans, thus generating a proper intracellular NM2A/NM2B distribution needed for efficient cell migration. Together, we show that NM2A motor activity is a key component of the bipolar filament depolymerization mechanism.

Pollen tubes are dynamic tip-growing cells that deliver sperm nuclei to female gametes in flowering plants, allowing for sexual reproduction and seed formation. Actin and microtubule cytoskeletons both play important roles in directional pollen tube growth and guidance. While actin dynamics are well-studied in pollen tubes, the role of microtubules and the interactions between these two cytoskeletal filaments are less well understood. To address this knowledge gap, we imaged growing Arabidopsis thaliana pollen tubes co-expressing fluorescently-labeled tubulin and actin markers and observed partial co-localization of actin and microtubule filaments. We found that treatment with microtubule disrupting drugs did not affect the actin cytoskeleton. In contrast, when actin filaments were depolymerized, microtubules in the medial region of pollen tubes were disrupted, while microtubules at the cell cortex remained intact. Thus, the microtubule cytoskeleton in A. thaliana pollen tubes relies on the actin cytoskeleton in a spatially dependent manner. Furthermore, we utilized native expression of the microtubule plus-end binding protein EB1b to track microtubule orientation in growing pollen tubes. We found the microtubule array to be largely parallel, with plus ends growing away from the tube apex. Together, these findings offer new insights into the dynamics and organization of microtubules in growing pollen tubes and the interactions between actin filaments and microtubules.

Cellular transport processes along microtubules are often facilitated by multi-motor complexes, which are connected by adapter proteins and cargoes. The nuclear pore protein Nup358, for example, interacts with the dynein adapter Bicaudal D2 (BicD2), which in turn recruits minus-end directed dynein motors and plus-end directed kinesin-1 motors for a nuclear positioning pathway that is essential for brain development. How motor recruitment is regulated by interactions of BicD2 with Nup358 is not well understood. Here, we characterize the structure of a minimal complex of kinesin-1 light chain 2 (KLC2), Nup358 and BicD2 by cryo-electron microscopy and small angle X-ray scattering. KLC2/Nup358 assumes a rod-like structure that increases in thickness, when BicD2 is bound. The addition of BicD2 also shifts the KLC2/Nup358/BicD2 complex towards a 2:2:2 stoichiometry, promoting dimerization at lower protein concentrations than without BicD2. Similarly, the presence of the Nup358/KLC2 interaction results in a shift towards a 2:2:2 stoichiometry. Based on these results, we hypothesize that KLC2 and BicD2 are recruited to Nup358 in a cooperative manner, and cooperativity may be promoted by modulation of the oligomeric state.

Centriolar satellites are dynamic pericentrosomal structures implicated in centrosomal protein homeostasis and ciliogenesis. Centriolar satellites have been identified in vertebrates and were only recently described in flies. In C. elegans similar pericentriolar structures were reported for the Sjogrens Syndrome Nuclear Antigen 1 (SSNA-1). However, whether these foci have characteristics resembling centriolar satellites of vertebrates, has not been explored. We show that Spindle Assembly-1 (SAS-1), the interaction partner of SSNA-1 forms similar satellite-like structures that localize to a pericentrosomal space in a cell cycle-dependent manner. SAS-1 satellite-like structures associate with and are dependent on the microtubule cytoskeleton. Furthermore, we demonstrate that they form in a dose dependent manner, are dynamic and sensitive to agents disrupting weak hydrophobic interactions, characteristics of biomolecular condensates. We conclude that C. elegans has bona fide centriolar satellites highlighting their evolutionary conservation and importance across species, and at the same time opening new avenues for future mechanistic studies.

Vimentin intermediate filaments play essential roles in maintaining cell integrity and regulating numerous cellular functions. In particular, vimentin cooperates with the actin cytoskeleton in key cellular processes that rely on actin dynamics, such as migration, division, and mechanosensing. While there is evidence that these two cytoskeletal components interact in cells, the underlying molecular mechanisms are only partially understood. Actin and vimentin can interact through biochemical signaling pathways or via cross-linkers, but whether they engage in a direct protein-protein interaction has remained controversial, in part because such interactions are difficult to isolate and characterize in cells. Using in vitro reconstitution coupled to theoretical modeling, and total internal reflection fluorescence microscopy to monitor the elongation of single actin filaments, we show that vimentin promotes actin elongation by stabilizing actin subunits at the barbed end in a dose-dependent manner. Strikingly, this effect depends on the nucleotide state of actin, as the acceleration is only observed for the elongation from ATP-actin, and not ADP-actin monomers. We further establish that neither the vimentin tail nor head domains are required for this effect, and both filamentous and non-filamentous vimentin enhance actin elongation. Finally, we find that vimentin promotes the nucleation of actin filaments. Consistently, magnetic pull-down assays demonstrate a direct interaction between vimentin and ATP-actin monomers. Altogether, these findings identify vimentin as an unexpected new actor in the regulation of actin dynamics at the barbed end and bring new insights into the functional role of vimentin through cytoskeletal crosstalk.

Variants in ACTB gene encoding for cytoplasmic {beta}-actin result in a group of rare disorders called non-muscle actinopathies (NMA). We investigated the cellular effects of a missense variant, G302A, and a four-amino-acid deletion, S338-I341, associated with the subgroup of NMA - ACTB pLoF (predicted loss-of-function) disorder in patient-derived fibroblast cells. We found that neither of the mutations affected the organization of actin or the width of the actin-filament bundles, while the mutation G302A reduced the stiffness of the cells as measured by using atomic force microscopy. The latter effect might be associated with the misorganization of tubulin and with the increased size and number of focal adhesions. When we challenged the cells by monolayer stretching and followed the mechanically-induced reorganization of the actin cytoskeleton, we found that G302A mutant cells showed more dense actin filament bundles within the cells compared to wild type cells. At the same time, the extent of cofilin reorganization from the cell periphery was increased upon stretch, and this correlated with an increased cofilin phosphorylation. In the case of the deletion, while the extent of cofilin phosphorylation increased, the extent of reorganization was unaltered; rather, the phosphorylation of myosin light chain, important in counteracting external force, was drastically reduced. We could partially rescue this fascinating effect by overexpressing the active form of the formin mDia. Our findings open the possibility to validate the cellular phenotype in the most affected patients cells, in neurons.

1Nuclear envelope stretch and rupture are common to cell spreading and migration in a variety of microenvironments, leading to marked changes in nucleocytoplasmic transport. Predicting cell response to different mechanochemical cues that are transmitted to the nucleus remains an open problem in the field of mechanomedicine. We developed a predictive modeling framework to examine how nuclear deformation on substrates with different nanotopographies influences nucleocytoplasmic transport and rearrangement of the nuclear lamina. Using the finite element method, we simulated nuclear compression by the perinuclear actin cap on substrates with arrays of nanopillars, modeling the nuclear envelope as a nonlinear elastic structure and coupling deformations to a biochemical model of lamin remodeling and nucleocytoplasmic transport. These simulations predicted regions of high nuclear envelope stretch adjacent to cell-nanopillar contacts, leading to maximized nuclear envelope tension on small nanopillars spaced by 4-5 microns. We then considered the effects on nuclear transport of YAP and TAZ and found that increased nuclear compression led to YAP/TAZ nuclear localization in agreement with previous experiments. Furthermore, the simulated force load per lamin was maximized on nanopillar substrates with high nuclear stretch. The magnitude of this load was modulated by the rate of actin cap assembly and the overall expression level of lamin A/C - decreasing lamin content in the nuclear envelope led to a higher likelihood of rupture. We validated this prediction in subsequent experiments with lamin-depleted U2OS cells, establishing the central importance of lamin transport and microenvironment nanotopography to nuclear mechanotransduction. 2 SignificanceCell nuclei commonly experience large strains, but existing computational models do not explain the coupling between such deformations and molecular transport. Here, we present a modeling framework that includes the mechanics of nuclear deformations and the reaction-transport of molecules within the cytoplasm, nuclear envelope, and nuclear interior. As a well-controlled setup for comparing experiments and simulations, we consider nuclear indentations exhibited by cells on nanopillar substrates. Our simulations recapitulate measurements of nuclear YAP/TAZ localization from the literature and predict that low-lamin cells experience higher force loads at the nuclear envelope. We validate this prediction experimentally, showing that lamin-depleted cells are more likely to exhibit nuclear rupture. Overall, our framework presents opportunities to predict nuclear mechanoadaptation to different microenvironments.

Embryonic development in multicellular organisms exhibits diverse morphogenetic patterns, which can generally be categorized into fundamental types such as monolayer and multilayer spheres, as well as cell masses. Furthermore, we identify two distinct processes for the formation of spherical structures. These basic patterns are thought to be governed by the microscopic properties of intercellular adhesion. However, the specific mechanisms linking the microscopic factors to the emergence of distinct macroscopic morphogenetic patterns remain poorly understood. In this study, we explore how different morphogenetic patterns arise by employing a computational model that incorporates intercellular adhesion and polarity. Our results demonstrate that all fundamental morphogenetic patterns can be generated through the interplay of two key parameters: the strength of cell polarity and the regulation of polarity via mechanical signals. Furthermore, analytical discussions reveal partial mechanisms underlying the formation of these patterns. These findings highlight the critical role of physical constraints in morphogenesis and suggest potential applications in the design principles for artificial tissues and organoids. Author summaryLiving organisms build their bodies through morphogenesis, during which cells autonomously arrange themselves into functional structures such as sheets, tubes, and spheres. From simple monolayered spheres to complex multilayered tissues organized by adhesion, it remains unclear how such diverse forms arise. Here, we mathematically modeled a population of proliferating cells governed only by two microscopic factors: the strength of polarity-dependent adhesion and the time scale at which polarity is regulated by cell-cell contact. Surprisingly, we found that this minimal model reproduces five basic morphological types observed in living embryos, including monolayer/multilayer structures and two distinct modes of cavity formation: by wrapping around or by inflating from the inside. Systematic simulations revealed that these macroscale outcomes are determined solely by two parameters controlling polarity strength and its regulation, suggesting that simple physical rules underlie diverse developmental architectures. Analysis of the model uncovers phase transitions between the five morphogenetic types and reveals how varying polarity and adhesion can recapitulate features of real embryogenesis. Our work proposes a unified framework that connects microscopic polarity mechanics to diverse developmental morphologies and provides a foundation for future applications in organoid design and tissue engineering.

Microtubule networks are major determinants of cell architecture and logistics. Microtubule organization and density are regulated by severing enzymes, which cut microtubule lattices or affect their growth and shortening. These activities can lead to microtubule amplification or disassembly, depending on the presence of microtubule stabilizers or destabilizers, but the interplay between these factors is poorly understood. Here, we reconstituted in vitro the activity of microtubule severase katanin together with microtubule minus-end stabilizers CAMSAPs, their binding partner WDR47 and microtubule depolymerase kinesin-13/MCAK. We confirmed that katanin can amplify or destroy microtubules in a concentration-dependent manner. CAMSAPs recruit katanin to microtubules and reduce katanin concentration needed for both amplification and destruction, whereas kinesin-13 completely abolishes microtubule amplification. WDR47 binds to microtubules decorated by CAMSAPs and suppresses katanin binding and severing. In addition, both katanin and WDR47 inhibit polymerization of CAMSAP-decorated microtubule minus ends. These data explain how these proteins act together to fine-tune microtubule minus-end stability without strongly increasing microtubule abundance. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=169 SRC="FIGDIR/small/714132v1_ufig1.gif" ALT="Figure 1"> View larger version (20K): org.highwire.dtl.DTLVardef@746fe3org.highwire.dtl.DTLVardef@5dd5a8org.highwire.dtl.DTLVardef@762373org.highwire.dtl.DTLVardef@1192db_HPS_FORMAT_FIGEXP M_FIG Graphical abstract C_FIG

The dynamic assembly of actin filaments underlies diverse cellular morphologies such as lamellipodia, filopodia, and reticulated networks. However, how filament-scale interactions among actin-binding proteins produce distinct actin architectures remains unclear. We developed a filament-resolved computational model of actin self-organization regulated by the Arp2/3 complex and fascin. Individual F-actin filaments are represented as elastic chains, and their stochastic polymerization, Arp2/3-mediated branching, and fascin-mediated crosslinking and bundling are explicitly modeled. The simulations reproduce three actin architectures observed in minimal reconstitution experiments, including lamellipodia-like branched networks, filopodia-like bundled protrusions, and reticulated meshworks, as a function of Arp2/3 and fascin concentrations. We quantify these regimes using actin density, orientational order, and spikiness, which robustly separate the three morphologies across conditions. To connect filament organization to shape change, we further couple the actin network to membrane deformation using a phase-field formulation. This coupling shows how localized remodeling concentrates load to drive pseudopodial protrusions, whereas highly branched networks distribute stresses and stabilize rounded shapes. The model links molecular interactions to emergent architecture and cell-scale morphodynamics.

Haptotaxis is an understudied form of directed cell migration in which movements are biased by gradients of immobilized ligands. For example, fibroblasts and other mesenchymal cells sense and respond to gradients of extracellular matrix (ECM) composition, which is relevant during tissue morphogenesis and repair. As a step towards understanding how haptotactic gradients spatially bias cell adhesion, intracellular signal transduction, and cytoskeletal dynamics, we formulated a phase field model of whole-cell migration, in which the occupancy of potential adhesion sites changes stochastically with time. With careful assignment of parameter values, the model predicts significant haptotactic bias for adhesion-site gradient steepness of a few percent across the cell. We then used the model to predict how the cells removal of surface-bound ECM ligand (as observed in experiment) and/or the presence of a competing, chemotactic gradient influence(s) haptotactic fidelity. An emergent principle is that gains in directional persistence naturally offset losses of directional bias, at the cost of greater cell-to-cell heterogeneity of the response. In the case of orthogonally oriented gradients, this offset manifests as a remarkable robustness of the multi-cue response.