Collective oscillations in microtubule growth

Oscillations, travelling fronts and patterns in a supramolecular system

Supramolecular polymers, such as microtubules, operate under non-equilibrium conditions to drive crucial functions in cells, such as motility, division and organelle transport¹. In vivo and in vitro size oscillations of individual microtubules^2,3 (dynamic instabilities) and collective oscillations⁴ have been observed. In addition, dynamic spatial structures, like waves and polygons, can form in non-stirred systems⁵. Here we describe an artificial supramolecular polymer made of a perylene diimide derivative that displays oscillations, travelling fronts and centimetre-scale self-organized patterns when pushed far from equilibrium by chemical fuels. Oscillations arise from a positive feedback due to nucleation-elongation-fragmentation, and a negative feedback due to size-dependent depolymerization. Travelling fronts and patterns form due to self-assembly induced density differences that cause system-wide convection. In our system, the species responsible for the nonlinear dynamics and those that self-assemble are one and the same. In contrast, other reported oscillating assemblies formed by vesicles⁶, micelles⁷ or particles⁸ rely on the combination of a known chemical oscillator and a stimuli-responsive system, either by communication through the solvent (for example, by changing pH^7-9), or by anchoring one of the species covalently (for example, a Belousov-Zhabotinsky catalyst^6,10). The design of self-oscillating supramolecular polymers and large-scale dissipative structures brings us closer to the creation of more life-like materials¹¹ that respond to external stimuli similarly to living cells, or to creating artificial autonomous chemical robots¹².

Effects of aging in catastrophe on the steady state and dynamics of a microtubule population

Several independent observations have suggested that the catastrophe transition in microtubules is not a first-order process, as is usually assumed. Recent in vitro observations by Gardner et al. [M. K. Gardner et al., Cell 147, 1092 (2011)] showed that microtubule catastrophe takes place via multiple steps and the frequency increases with the age of the filament. Here we investigate, via numerical simulations and mathematical calculations, some of the consequences of the age dependence of catastrophe on the dynamics of microtubules as a function of the aging rate, for two different models of aging: exponential growth, but saturating asymptotically, and purely linear growth. The boundary demarcating the steady-state and non-steady-state regimes in the dynamics is derived analytically in both cases. Numerical simulations, supported by analytical calculations in the linear model, show that aging leads to nonexponential length distributions in steady state. More importantly, oscillations ensue in microtubule length and velocity. The regularity of oscillations, as characterized by the negative dip in the autocorrelation function, is reduced by increasing the frequency of rescue events. Our study shows that the age dependence of catastrophe could function as an intrinsic mechanism to generate oscillatory dynamics in a microtubule population, distinct from hitherto identified ones.

Models of the collective behavior of proteins in cells: tubulin, actin and motor proteins

One of the most important issues of molecular biophysics is the complex and multifunctional behavior of the cell's cytoskeleton. Interiors of living cells are structurally organized by the cytoskeleton networks of filamentous protein polymers: microtubules, actin and intermediate filaments with motor proteins providing force and directionality needed for transport processes. Microtubules (MT's) take active part in material transport within the cell, constitute the most rigid elements of the cell and hence found many uses in cell motility (e.g. flagella andcilia). At present there is, however, no quantitatively predictable explanation of how these important phenomena are orchestrated at a molecular level. Moreover, microtubules have been demonstrated to self-organize leading to pattern formation. We discuss here several models which attempt to shed light on the assembly of microtubules and their interactions with motor proteins. Subsequently, an overview of actin filaments and their properties isgiven with particular emphasis on actin assembly processes. The lengths of actin filaments have been reported that were formed by spontaneous polymerization of highly purified actin monomers after labeling with rhodamine-phalloidin. The length distributions are exponential with a mean of about 7 μm. This length is independent of the initial concentration of actin monomer, an observation inconsistent with a simple nucleation-elongation mechanism. However, with the addition of physically reasonable rates of filament annealing and fragmenting, a nucleation-elongation mechanism can reproduce the observed average length of filaments in two types of experiments: (1) filaments formed from a wide range of highly purified actin monomer concentrations, and (2) filaments formed from 24 mM actin over a range of CapZ concentrations. In the final part of the paper we briefly review the stochastic models used to describe the motion of motor proteins on protein filaments. The vast majority of these models are based on ratchet potentials with the presence of thermal noise and forcing due to ATP binding and a subsequent hydrolysis. Many outstanding questions remain to be quantitatively addressed on a molecular level in order to explain the structure-to-function relationship for the key elements of the cytoskeleton discussed in this review.

Models of assembly and disassembly of individual microtubules: stochastic and averaged equations

In this paper we present solutions of the master equations for the microtubule length and show that the local probability for rescues or catastrophes can lead to bell-shaped length histograms. Conversely, as already known, non-local probabilities for these events result in exponential length histograms. We also derive master equations for a stabilizing cap and obtain a new boundary condition which provides an explanation of the results obtained in dilution and cutting experiments.

Continuous model for microtubule dynamics with catastrophe, rescue, and nucleation processes

Microtubules are a major component of the cytoskeleton distinguished by highly dynamic behavior both in vitro and in vivo referred to as dynamic instability. We propose a general mathematical model that accounts for the growth, catastrophe, rescue, and nucleation processes in the polymerization of microtubules from tubulin dimers. Our model is an extension of various mathematical models developed earlier formulated in order to capture and unify the various aspects of tubulin polymerization. While attempting to use a minimal number of adjustable parameters, the proposed model covers a broad range of behaviors and has predictive features discussed in the paper. We have analyzed the range of resultant dynamical behavior of the microtubules by changing each of the parameter values at a time and observing the emergence of various dynamical regimes that agree well with the previously reported experimental data and behavior.

Liénard systems and potential-Hamiltonian decomposition: applications in biology

In separated notes, we described the mathematical aspects of the potential-Hamiltonian (PH) decomposition, in particular, for n-switches and Liénard systems [J. Demongeot, N. Glade, L. Forest, Liénard systems and potential-Hamiltonian decomposition - I. Methodology, II. Algorithm and III. Applications, C. R. Acad. Sci., Paris, Ser. I, in press]. In the present note, we give some examples of biological regulatory systems susceptible to be decomposed. We show that they can be modelled in terms of 2D ordinary differential equations belonging to n-switches and Liénard system families [O. Cinquin, J. Demongeot, High-dimensional switches and the modeling of cellular differentiation, J. Theor. Biol. 233 (2005) 391-411]. Although simplified, these models can be decomposed into a set of equations combining a potential and a Hamiltonian part. We discuss about the advantage of such a PH-decomposition for understanding the mechanisms involved in their regulatory abilities. We suggest a generalized algorithm to deal with differential systems having a second part of rational-fraction type (frequently used in metabolic systems). Finally, we comment what can be interpreted as a precise signification in biological systems from the dynamical behaviours of both the potential and Hamiltonian parts.

Dynamic instability of microtubules: effect of catastrophe-suppressing drugs

Microtubules are stiff filamentary proteins that constitute an important component of the cytoskeleton of cells. These are known to exhibit a dynamic instability. A steadily growing microtubule can suddenly start depolymerizing very rapidly; this phenomenon is known as a "catastrophe." However, often a shrinking microtubule is "rescued" and starts polymerizing again. Here we develop a model for the polymerization-depolymerization dynamics of microtubules in the presence of catastrophe-suppressing drugs. Solving the dynamical equations in the steady state, we derive exact analytical expressions for the length distributions of the microtubules tipped with drug-bound tubulin subunits as well as those of the microtubules, in the growing and shrinking phases, tipped with drug-free pure tubulin subunits. We also examine the stability of the steady-state solutions.

Effect of tubulin diffusion on polymerization of microtubules

The dynamics of microtubules (MT's) growing from a nucleation center is simulated with a kinetic Monte Carlo model that includes tubulin diffusion. In the limit of fast diffusion (homogeneous tubulin concentration), MT growth is synchronous and bounded. The microtubules form an aster with a monotonously decreasing long-time distribution of lengths. Slow tubulin diffusion leads to rapid dephasing in the growth dynamics, unbounded growth of some MT's, spatial inhomogeneities, and morphological change toward a morphology with bounded short MT's located in the nucleation center and unbounded long MT's with narrowly distributed lengths. The transition from unbounded to bounded growth is driven by the competition between the reaction rate of the tubulin assembly and the tubulin's diffusion rate. While the present study reports the effect of the tubulin diffusion coefficient on the transition, the results of the simulations are qualitatively comparable to the morphological and dynamical changes of centrosome-nucleated MT's from interphase to mitosis in cellular systems where the transition is regulated by the reaction rates.

Modeling oscillatory microtubule polymerization

Polymerization of microtubules is ubiquitous in biological cells and under certain conditions it becomes oscillatory in time. Here, simple reaction models are analyzed that capture such oscillations as well as the length distribution of microtubules. We assume reaction conditions that are stationary over many oscillation periods, and it is a Hopf bifurcation that leads to a persistent oscillatory microtubule polymerization in these models. Analytical expressions are derived for the threshold of the bifurcation and the oscillation frequency in terms of reaction rates, and typical trends of their parameter dependence are presented. Both, a catastrophe rate that depends on the density of guanosine triphosphate liganded tubulin dimers and a delay reaction, such as the depolymerization of shrinking microtubules or the decay of oligomers, support oscillations. For a tubulin dimer concentration below the threshold, oscillatory microtubule polymerization occurs transiently on the route to a stationary state, as shown by numerical solutions of the model equations. Close to threshold, a so-called amplitude equation is derived and it is shown that the bifurcation to microtubule oscillations is supercritical.

A Landau-Ginzburg Model of the Co-existence of Free Tubulin and Assembled Microtubules in Nucleation and Oscillations Phenomena

A link is shown between reaction-diffusion kinetics for microtubuleassembly and time-dependent Landau-Ginzburg phenomenology. In the latter,microtubule assembly is treated as a first-order phase transition using apostulated Landau-Ginzburg free energy expansion. The results establish aconnection between the oscillations observed in experiment and the phasediagram for microtubule assembly. The model also predicts a specific heatbehavior which could be verified experimentally.

Model for spatial microtubule oscillations

Under particular in vitro conditions, oscillating spatial and temporal waves of assembled microtubules can be observed. A reaction-diffusion model is presented to reproduce these results. This model is based on a set of chemical reaction equations and extended to include spatial dependence and diffusion. The basic properties of the model are presented and the results are demonstrated to connect the observable waves with turbidimetric measurements. The results of the model are consistent with experimental findings.

Microtubules: strange polymers inside the cell

This paper provides a consistent approach (within a one-dimensional approximation) to the description of the evolution of the microtubule length at both low- and high-density concentrations. We derive general master-type equations which are based on the key chemical reactions involved in the assembly and disassembly of microtubules. The processes included are: polymerization and depolymerization of a single protein dimer, catastrophic disassembly affecting an a piori arbitrary number of dimers, and a rescue event. Solutions of the derived equations are compared with the existing experimental data. Important conclusions linking the emergence of bell-shaped histograms with the nature of catastrophe and rescue phenomena are drawn. Finally, we briefly discuss the emergence of coherent phenomena in microtubule polymerization, i.e., a transition to collective oscillations in the assembly and disassembly effects.