The emergence of sarcomeric, graded-polarity and spindle-like patterns in bundles of short cytoskeletal polymers and two opposite molecular motors

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.

A computational model of the early stages of acentriolar meiotic spindle assembly

The mitotic spindle is an ensemble of microtubules responsible for the repartition of the chromosomal content between the two daughter cells during division. In metazoans, spindle assembly is a gradual process involving dynamic microtubules and recruitment of numerous associated proteins and motors. During mitosis, centrosomes organize and nucleate the majority of spindle microtubules. In contrast, oocytes lack canonical centrosomes but are still able to form bipolar spindles, starting from an initial ball that self-organizes in several hours. Interfering with early steps of meiotic spindle assembly can lead to erroneous chromosome segregation. Although not fully elucidated, this process is known to rely on antagonistic activities of plus end- and minus end-directed motors. We developed a model of early meiotic spindle assembly in mouse oocytes, including key factors such as microtubule dynamics and chromosome movement. We explored how the balance between plus end- and minus end-directed motors, as well as the influence of microtubule nucleation, impacts spindle morphology. In a refined model, we added spatial regulation of microtubule stability and minus-end clustering. We could reproduce the features of early stages of spindle assembly from 12 different experimental perturbations and predict eight additional perturbations. With its ability to characterize and predict chromosome individualization, this model can help deepen our understanding of spindle assembly.

Polarity sorting of axonal microtubules: a computational study

We present a computational model to test a "polarity sorting" mechanism for microtubule (MT) organization in developing axons. We simulate the motor-based axonal transport of short MTs to test the hypothesis that immobilized cytoplasmic dynein motors transport short MTs with their plus ends leading, so "mal-oriented" MTs with minus-end-out are transported toward the cell body while "correctly" oriented MTs are transported in the anterograde direction away from the soma. We find that dynein-based transport of short MTs can explain the predominately plus-end-out polarity pattern of axonal MTs but that transient attachments of plus-end-directed motor proteins and nonmotile cross-linker proteins are needed to explain the frequent pauses and occasional reversals observed in live-cell imaging of MT transport. Static cross-linkers increase the likelihood of a stalled "tug-of-war" between retrograde and anterograde forces on the MT, providing an explanation for the frequent pauses of short MTs and the immobility of longer MTs. We predict that inhibition of the proposed static cross-linker will produce disordered transport of short MTs and increased mobility of longer MTs. We also predict that acute inhibition of cytoplasmic dynein will disrupt the polarity sorting of MTs by increasing the likelihood of "incorrect" sorting of MTs by plus-end-directed motors.

A Combination of Actin Treadmilling and Cross-Linking Drives Contraction of Random Actomyosin Arrays

We investigate computationally the self-organization and contraction of an initially random actomyosin ring. In the framework of a detailed physical model for a ring of cross-linked actin filaments and myosin-II clusters, we derive the force balance equations and solve them numerically. We find that to contract, actin filaments have to treadmill and to be sufficiently cross linked, and myosin has to be processive. The simulations reveal how contraction scales with mechanochemical parameters. For example, they show that the ring made of longer filaments generates greater force but contracts slower. The model predicts that the ring contracts with a constant rate proportional to the initial ring radius if either myosin is released from the ring during contraction and actin filaments shorten, or if myosin is retained in the ring, while the actin filament number decreases. We demonstrate that a balance of actin nucleation and compression-dependent disassembly can also sustain contraction. Finally, the model demonstrates that with time pattern formation takes place in the ring, worsening the contractile process. The more random the actin dynamics are, the higher the contractility will be.

MEDYAN: Mechanochemical Simulations of Contraction and Polarity Alignment in Actomyosin Networks

Active matter systems, and in particular the cell cytoskeleton, exhibit complex mechanochemical dynamics that are still not well understood. While prior computational models of cytoskeletal dynamics have lead to many conceptual insights, an important niche still needs to be filled with a high-resolution structural modeling framework, which includes a minimally-complete set of cytoskeletal chemistries, stochastically treats reaction and diffusion processes in three spatial dimensions, accurately and efficiently describes mechanical deformations of the filamentous network under stresses generated by molecular motors, and deeply couples mechanics and chemistry at high spatial resolution. To address this need, we propose a novel reactive coarse-grained force field, as well as a publicly available software package, named the Mechanochemical Dynamics of Active Networks (MEDYAN), for simulating active network evolution and dynamics (available at www.medyan.org). This model can be used to study the non-linear, far from equilibrium processes in active matter systems, in particular, comprised of interacting semi-flexible polymers embedded in a solution with complex reaction-diffusion processes. In this work, we applied MEDYAN to investigate a contractile actomyosin network consisting of actin filaments, alpha-actinin cross-linking proteins, and non-muscle myosin IIA mini-filaments. We found that these systems undergo a switch-like transition in simulations from a random network to ordered, bundled structures when cross-linker concentration is increased above a threshold value, inducing contraction driven by myosin II mini-filaments. Our simulations also show how myosin II mini-filaments, in tandem with cross-linkers, can produce a range of actin filament polarity distributions and alignment, which is crucially dependent on the rate of actin filament turnover and the actin filament’s resulting super-diffusive behavior in the actomyosin-cross-linker system. We discuss the biological implications of these findings for the arc formation in lamellipodium-to-lamellum architectural remodeling. Lastly, our simulations produce force-dependent accumulation of myosin II, which is thought to be responsible for their mechanosensation ability, also spontaneously generating myosin II concentration gradients in the solution phase of the simulation volume.

Active matter systems have the distinct ability to convert energy from their surroundings into mechanical work, which gives rise to them having highly dynamic properties. Modeling active matter systems and capturing their complex behavior has been a great challenge in past years due to the many coupled interactions between their constituent parts, including not only distinct chemical and mechanical properties, but also feedback between them. One of the most intriguing biological active matter systems is the cell cytoskeleton, which can dynamically respond to chemical and mechanical cues to control cell structure and shape, playing a central role in many higher-order cellular processes. To model these systems and reproduce their behavior, we present a new modeling approach which combines the chemical, mechanical, and molecular transport aspects of active matter systems, all represented with equivalent complexity, while also allowing for various forms of mechanochemical feedback. This modeling approach, named MEDYAN, and software implementation is flexible so that a wide range of active matter systems can be simulated with a high level of detail, and ultimately can help to describe active matter phenomena, and in particular, the dynamics of the cell cytoskeleton. In this work, we have used MEDYAN to simulate a cytoskeletal network consisting of actin filaments, cross-linking proteins, and myosin II molecular motors. We found that these systems show rich dynamical behaviors, undergoing alignment and bundling transitions, with an emergent contractility, as the concentrations of myosin II and cross-linking proteins, as well as actin filament turnover rates, are varied.

Architecture and Connectivity Govern Actin Network Contractility

Actomyosin contractility plays a central role in a wide range of cellular processes, including the establishment of cell polarity, cell migration, tissue integrity, and morphogenesis during development. The contractile response is variable and depends on actomyosin network architecture and biochemical composition. To determine how this coupling regulates actomyosin-driven contraction, we used a micropatterning method that enables the spatial control of actin assembly. We generated a variety of actin templates and measured how defined actin structures respond to myosin-induced forces. We found that the same actin filament crosslinkers either enhance or inhibit the contractility of a network, depending on the organization of actin within the network. Numerical simulations unified the roles of actin filament branching and crosslinking during actomyosin contraction. Specifically, we introduce the concept of "network connectivity" and show that the contractions of distinct actin architectures are described by the same master curve when considering their degree of connectivity. This makes it possible to predict the dynamic response of defined actin structures to transient changes in connectivity. We propose that, depending on the connectivity and the architecture, network contraction is dominated by either sarcomeric-like or buckling mechanisms. More generally, this study reveals how actin network contractility depends on its architecture under a defined set of biochemical conditions.

Force to divide: structural and mechanical requirements for actomyosin ring contraction

One of the unresolved questions in the field of cell division is how the actomyosin cytoskeleton remains structurally organized while generating the contractile force to divide one cell into two. In analogy to the actomyosin-based force production mechanism in striated muscle, it was originally proposed that contractile stress in the actomyosin ring is generated via a sliding filament mechanism within an organized sarcomere-like array. However, over the last 30 years, ultrastructural and functional studies have noted important distinctions between cytokinetic structures in dividing cells and muscle sarcomeres. Myosin-II motor activity is not always required, and there is evidence that actin depolymerization contributes to contraction. In this Review, the architecture and contractile dynamics of the actomyosin ring at the cell division plane will be discussed. We will report the interdisciplinary advances in the field as well as their integration into a mechanistic understanding of contraction in cell division and in other biological processes that rely on an actomyosin-based force-generating system.

Mathematical modeling of eukaryotic cell migration: insights beyond experiments

A migrating cell is a molecular machine made of tens of thousands of short-lived and interacting parts. Understanding migration means understanding the self-organization of these parts into a system of functional units. This task is one of tackling complexity: First, the system integrates numerous chemical and mechanical component processes. Second, these processes are connected in feedback interactions and over a large range of spatial and temporal scales. Third, many processes are stochastic, which leads to heterogeneous migration behaviors. Early on in the research of cell migration it became evident that this complexity exceeds human intuition. Thus, the cell migration community has led the charge to build mathematical models that could integrate the diverse experimental observations and measurements in consistent frameworks, first in conceptual and more recently in molecularly explicit models. The main goal of this review is to sift through a series of important conceptual and explicit mathematical models of cell migration and to evaluate their contribution to the field in their ability to integrate critical experimental data.

Sarcomeric pattern formation by actin cluster coalescence

Contractile function of striated muscle cells depends crucially on the almost crystalline order of actin and myosin filaments in myofibrils, but the physical mechanisms that lead to myofibril assembly remains ill-defined. Passive diffusive sorting of actin filaments into sarcomeric order is kinetically impossible, suggesting a pivotal role of active processes in sarcomeric pattern formation. Using a one-dimensional computational model of an initially unstriated actin bundle, we show that actin filament treadmilling in the presence of processive plus-end crosslinking provides a simple and robust mechanism for the polarity sorting of actin filaments as well as for the correct localization of myosin filaments. We propose that the coalescence of crosslinked actin clusters could be key for sarcomeric pattern formation. In our simulations, sarcomere spacing is set by filament length prompting tight length control already at early stages of pattern formation. The proposed mechanism could be generic and apply both to premyofibrils and nascent myofibrils in developing muscle cells as well as possibly to striated stress-fibers in non-muscle cells.

Muscle contraction driving voluntary movements and the beating of the heart relies on the contraction of highly regular bundles of actin and myosin filaments, which share a periodic, sarcomeric pattern. We know little about the mechanisms by which these ‘biological crystals’ are assembled and it is a general question how order on a scale of 100 micrometers can emerge from the interactions of micrometer-sized building blocks, such as actin and myosin filaments. In our paper, we consider a computational model for a bundle of actin filaments and discuss physical mechanisms by which periodic order emerges spontaneously. Mutual crosslinking of actin filaments results in the formation and coalescence of growing actin clusters. Active elongation and shrinkage dynamics of actin filaments generates polymerization forces and causes local actin flow that can act like a conveyor belt to sort myosin filaments in place.