No Arabic abstract
Cortical actin networks are highly dynamic and play critical roles in shaping the mechanical properties of cells. The actin cytoskeleton undergoes significant reorganization over the course of the cell cycle, when cortical actin transitions between open patched meshworks, homogeneous distributions, and aligned bundles. Several types of myosin motor proteins, characterized by different kinetic parameters, have been involved in this reorganization of actin filaments. Given the limitations in studying the interactions of actin with myosin in vivo, we propose stochastic agent-based model simulations and develop a set of data analysis measures to assess how myosin motor proteins mediate various actin organizations. In particular, we identify individual motor parameters, such as motor binding rate and step size, that generate actin networks with different levels of contractility and different patterns of myosin motor localization. In simulations where two motor populations with distinct kinetic parameters interact with the same actin network, we find that motors may act in a complementary way, by tuning the actin network organization, or in an antagonistic way, where one motor emerges as dominant. This modeling and data analysis framework also uncovers parameter regimes where spatial segregation between motor populations is achieved. By allowing for changes in kinetic rates during the actin-myosin dynamic simulations, our work suggests that certain actin-myosin organizations may require additional regulation beyond mediation by motor proteins in order to reconfigure the cytoskeleton network on experimentally-observed timescales.
The living cell is an open nonequilibrium biochemical system, where ATP hydrolysis serves as the energy source for a wide range of intracellular processes including the assurance for decision-making. In the fission yeast cell cycle, the transition from G2 phase to M phase is triggered by the activation of Cdc13/Cdc2 and Cdc25, and the deactivation of Wee1. Each of these three events involves a phosphorylation-dephosphorylation (PdP) cycle, and together they form a regulatory circuit with feedback loops. Almost all quantitative models for cellular networks in the past have invalid thermodynamics due to the assumption of irreversible enzyme kinetics. We constructed a thermodynamically realistic kinetic model of the G2/M circuit, and show that the phosphorylation energy ($Delta G$), which is determined by the cellular ATP/ADP ratio, critically controls the dynamics and the bistable nature of Cdc2 activation. Using fission yeast nucleoplasmic extract (YNPE), we are able to experimentally verify our model prediction that increased , being synergistic to the accumulation of Cdc13, drives the activation of Cdc2. Furthermore, Cdc2 activation exhibits bistability and hysteresis in response to changes in phosphorylation energy. These findings suggest that adequate maintenance of phosphorylation energy ensures the bistability and robustness of the activation of Cdc2 in the G2/M transition. Free energy might play a widespread role in biological decision-making processes, connecting thermodynamics with information processing in biology.
The interaction between actin filaments and microtubules is crucial for many eukaryotic cellular processes, such as, among others, cell polarization, cell motility and cellular wound healing. The importance of this interaction has long been recognised, yet very little is understood about both the underlying mechanisms and the consequences for the spatial (re)organization of the cellular cytoskeleton. At the same time, understanding the causes and the consequences of the interaction between different biomolecular components are key questions for emph{in vitro} research involving reconstituted biomolecular systems, especially in the light of current interest in creating minimal synthetic cells. In this light, recent emph{in vitro} experiments have shown that the actin-microtubule interaction mediated by the cytolinker TipAct, which binds to actin lattice and microtubule tip, causes the directed transport of actin filaments. We develop an analytical theory of dynamically unstable microtubules, nucleated from the center of a spherical cell, in interaction with actin filaments. We show that, depending on the balance between the diffusion of unbound actin filaments and propensity to bind microtubules, actin is either concentrated in the center of the cell, where the density of microtubules is highest, or becomes localized to the cell cortex.
In cells and in vitro assays the number of motor proteins involved in biological transport processes is far from being unlimited. The cytoskeletal binding sites are in contact with the same finite reservoir of motors (either the cytosol or the flow chamber) and hence compete for recruiting the available motors, potentially depleting the reservoir and affecting cytoskeletal transport. In this work we provide a theoretical framework to study, analytically and numerically, how motor density profiles and crowding along cytoskeletal filaments depend on the competition of motors for their binding sites. We propose two models in which finite processive motor proteins actively advance along cytoskeletal filaments and are continuously exchanged with the motor pool. We first look at homogeneous reservoirs and then examine the effects of free motor diffusion in the surrounding medium. We consider as a reference situation recent in vitro experimental setups of kinesin-8 motors binding and moving along microtubule filaments in a flow chamber. We investigate how the crowding of linear motor proteins moving on a filament can be regulated by the balance between supply (concentration of motor proteins in the flow chamber) and demand (total number of polymerised tubulin heterodimers). We present analytical results for the density profiles of bound motors, the reservoir depletion, and propose novel phase diagrams that present the formation of jams of motor proteins on the filament as a function of two tuneable experimental parameters: the motor protein concentration and the concentration of tubulins polymerized into cytoskeletal filaments. Extensive numerical simulations corroborate the analytical results for parameters in the experimental range and also address the effects of diffusion of motor proteins in the reservoir.
Proteins from the kinesin-8 family promote microtubule (MT) depolymerization, a process thought to be important for the control of microtubule length in living cells. In addition to this MT shortening activity, kinesin 8s are motors that show plus-end directed motility on MTs. Here we describe a simple model that incorporates directional motion and destabilization of the MT plus end by kinesin 8. Our model quantitatively reproduces the key features of length-vs-time traces for stabilized MTs in the presence of purified kinesin 8, including length-dependent depolymerization. Comparison of model predictions with experiments suggests that kinesin 8 depolymerizes processively, i.e., one motor can remove multiple tubulin dimers from a stabilized MT. Fluctuations in MT length as a function of time are related to depolymerization processivity. We have also determined the parameter regime in which the rate of MT depolymerization is length dependent: length-dependent depolymerization occurs only when MTs are sufficiently short; this crossover is sensitive to the bulk motor concentration.
The assumption of linear response of protein molecules to thermal noise or structural perturbations, such as ligand binding or detachment, is broadly used in the studies of protein dynamics. Conformational motions in proteins are traditionally analyzed in terms of normal modes and experimental data on thermal fluctuations in such macromolecules is also usually interpreted in terms of the excitation of normal modes. We have chosen two important protein motors - myosin V and kinesin KIF1A - and performed numerical investigations of their conformational relaxation properties within the coarse-grained elastic network approximation. We have found that the linearity assumption is deficient for ligand-induced conformational motions and can even be violated for characteristic thermal fluctuations. The deficiency is particularly pronounced in KIF1A where the normal mode description fails completely in describing functional mechanochemical motions. These results indicate that important assumptions of the theory of protein dynamics may need to be reconsidered. Neither a single normal mode, nor a superposition of such modes yield an approximation of strongly nonlinear dynamics.