No abstract

A method for simulating the growth of branched actin networks against obstacles has been developed. The method is based on simple stochastic events, including addition or removal of monomers at filament ends, capping of filament ends, nucleation of branches from existing filaments, and detachment of branches; the network structure for several different models of the branching process has also been studied. The models differ with regard to their inclusion of effects such as preferred branch orientations, filament uncapping at the obstacle, and preferential branching at filament ends. The actin ultrastructure near the membrane in lamellipodia is reasonably well produced if preferential branching in the direction of the obstacle or barbed-end uncapping effects are included. Uncapping effects cause the structures to have a few very long filaments that are similar to those seen in pathogen-induced "actin tails." The dependence of the growth velocity, branch spacing, and network density on the rate parameters for the various processes is quite different among the branching models. An analytic theory of the growth velocity and branch spacing of the network is described. Experiments are suggested that could distinguish among some of the branching models.

The growth of an actin network against an obstacle that stimulates branching locally is studied using several variants of a kinetic rate model based on the orientation-dependent number density of filaments. The model emphasizes the effects of branching and capping on the density of free filament ends. The variants differ in their treatment of side versus end branching and dimensionality, and assume that new branches are generated by existing branches (autocatalytic behavior) or independently of existing branches (nucleation behavior). In autocatalytic models, the network growth velocity is rigorously independent of the opposing force exerted by the obstacle, and the network density is proportional to the force. The dependence of the growth velocity on the branching and capping rates is evaluated by a numerical solution of the rate equations. In side-branching models, the growth velocity drops gradually to zero with decreasing branching rate, while in end-branching models the drop is abrupt. As the capping rate goes to zero, it is found that the behavior of the velocity is sensitive to the thickness of the branching region. Experiments are proposed for using these results to shed light on the nature of the branching process.

Patterns of waves, patches, and peaks of actin are observed experimentally in many living cells. Models of this phenomenon have been based on the interplay between filamentous actin (F-actin) and its nucleation promoting factors (NPF’s) such as Arp2/3 complex. Here we present an alternative biologically-motivated model for F-actin-NPF interaction based on properties of GTPases acting as NPFs. The model is a natural extension of a previous mathematical mini-model of small GTPases that generates a static cell polarization. GTPases (such as Cdc42, Rac) are known to promote actin nucleation, and to have active membrane-bound and inactive cytosolic forms. Like others, we assume that F-actin negative feedback shapes the observed patterns by suppressing the trailing edge of NPF-generated wave-fronts (hence localizing the activity spatially). We find that our NPF-actin model generates a rich set of behaviours, spanning a transition from static polarization to single pulses, reflecting waves, wave trains, and oscillations localized at the cell edge. The model is developed with simplicity in mind to investigate the interaction between NPF kinetics and negative feedback (weak vs strong for example), explain distinct types of pattern initiation mechanisms, patterns formed, and parameter regimes corresponding to such behaviours. It is shown that weak actin feedback yields static patterning, moderate feedback yields dynamical behaviour such as travelling waves, and strong feedback can lead to wave trains or total suppression of patterning. We use a novel nonlinear bifurcation analysis to explore the parameter space of this model and predict its behaviour with simulations validating those results.

Membrane deformation during endocytosis in yeast is driven by local, templated assembly of a sequence of proteins including polymerized actin and curvature-generating coat proteins such as clathrin. Actin polymerization is required for successful endocytosis, but it is not known by what mechanisms actin polymerization generates the required pulling forces. To address this issue, we develop a simulation method in which the actin network at the protein patch is modeled as an active gel. The deformation of the gel is treated using a finite-element approach. We explore the effects and interplay of three different types of force driving invagination: 1), forces perpendicular to the membrane, generated by differences between actin polymerization rates at the edge of the patch and those at the center; 2), the inherent curvature of the coat-protein layer; and 3), forces parallel to the membrane that buckle the coat protein layer, generated by an actomyosin contractile ring. We find that with optimistic estimates for the stall stress of actin gel growth and the shear modulus of the actin gel, actin polymerization can generate almost enough force to overcome the turgor pressure. In combination with the other mechanisms, actin polymerization can the force over the critical value.