Microtubule dynamics: Caps, catastrophes, and coupled hydrolysis

EB3-informed dynamics of the microtubule stabilizing cap during stalled growth

Microtubule stability is known to be governed by a stabilizing GTP/GDP-Pi cap, but the exact relation between growth velocity, GTP hydrolysis, and catastrophes remains unclear. We investigate the dynamics of the stabilizing cap through in vitro reconstitution of microtubule dynamics in contact with microfabricated barriers, using the plus-end binding protein GFP-EB3 as a marker for the nucleotide state of the tip. The interaction of growing microtubules with steric objects is known to slow down microtubule growth and accelerate catastrophes. We show that the lifetime distributions of stalled microtubules, as well as the corresponding lifetime distributions of freely growing microtubules, can be fully described with a simple phenomenological 1D model based on noisy microtubule growth and a single EB3-dependent hydrolysis rate. This same model is furthermore capable of explaining both the previously reported mild catastrophe dependence on microtubule growth rates and the catastrophe statistics during tubulin washout experiments.

Minimal Mechanisms of Microtubule Length Regulation in Living Cells

The microtubule cytoskeleton is responsible for sustained, long-range intracellular transport of mRNAs, proteins, and organelles in neurons. Neuronal microtubules must be stable enough to ensure reliable transport, but they also undergo dynamic instability, as their plus and minus ends continuously switch between growth and shrinking. This process allows for continuous rebuilding of the cytoskeleton and for flexibility in injury settings. Motivated by in vivo experimental data on microtubule behavior in Drosophila neurons, we propose a mathematical model of dendritic microtubule dynamics, with a focus on understanding microtubule length, velocity, and state-duration distributions. We find that limitations on microtubule growth phases are needed for realistic dynamics, but the type of limiting mechanism leads to qualitatively different responses to plausible experimental perturbations. We therefore propose and investigate two minimally-complex length-limiting factors: limitation due to resource (tubulin) constraints and limitation due to catastrophe of large-length microtubules. We combine simulations of a detailed stochastic model with steady-state analysis of a mean-field ordinary differential equations model to map out qualitatively distinct parameter regimes. This provides a basis for predicting changes in microtubule dynamics, tubulin allocation, and the turnover rate of tubulin within microtubules in different experimental environments.

Modelling microtubule dynamic instability: Microtubule growth, shortening and pause

Microtubules (MTs) are protein polymers found in all eukaryotic cells. They are crucial for normal cell development, providing structural support for the cell and aiding in the transportation of proteins and organelles. In order to perform these functions, MTs go through periods of relatively slow polymerization (growth) and very fast depolymerization (shortening), where the switch from growth to shortening is called a catastrophe and the switch from shortening to growth is called a rescue. Although MT dynamic instability has traditionally been described solely in terms of growth and shortening, MTs have been shown to pause for extended periods of time, however the reason for pausing is not well understood. Here, we present a new mathematical model to describe MT dynamics in terms of growth, shortening, and pausing. Typically, MT dynamics are defined by four key parameters which include the MT growth rate, shortening rate, frequency of catastrophe, and the frequency of rescue. We derive a mathematical expression for the catastrophe frequency in the presence of pausing, as well as expressions to describe the total time that MTs spend in a state of growth and pause. In addition to exploring MT dynamics in a control-like setting, we explore the implicit effect of stabilizing MT associated proteins (MAPs) and stabilizing and destabilizing chemotherapeutic drugs that target MTs on MT dynamics through variations in model parameters.

Coupling of mitochondrial population evolution to microtubule dynamics in fission yeast cells: a kinetic Monte Carlo study

Mitochondrial populations in cells are maintained by cycles of fission and fusion events. Perturbation of this balance has been observed in several diseases such as cancer and neurodegeneration. In fission yeast cells, the association of mitochondria with microtubules inhibits mitochondrial fission [Mehta et al., J. Biol. Chem., 2019, 294, 3385], illustrating the intricate coupling between mitochondria and the dynamic population of microtubules within the cell. In order to understand this coupling, we carried out kinetic Monte Carlo (KMC) simulations to predict the evolution of mitochondrial size distributions for different cases; wild-type cells, cells with short and long microtubules, and cells without microtubules. Comparisons are made with mitochondrial distributions reported in experiments with fission yeast cells. Using experimentally determined mitochondrial fission and fusion frequencies, simulations implemented without the coupling of microtubule dynamics predicted an increase in the mean number of mitochondria, equilibrating within 50 s. The mitochondrial length distribution in these models also showed a higher occurrence of shorter mitochondria, implying a greater tendency for fission, similar to the scenario observed in the absence of microtubules and cells with short microtubules. Interestingly, this resulted in overestimating the mean number of mitochondria and underestimating mitochondrial lengths in cells with wild-type and long microtubules. However, coupling mitochondria's fission and fusion events to the microtubule dynamics effectively captured the mitochondrial number and size distributions in wild-type and cells with long microtubules. Thus, the model provides greater physical insight into the temporal evolution of mitochondrial populations in different microtubule environments, allowing one to study both the short-time evolution as observed in the experiments (<5 minutes) as well as their transition towards a steady-state (>15 minutes). Our study illustrates the critical role of microtubules in mitochondrial dynamics and coupling microtubule growth and shrinkage dynamics is critical to predicting the evolution of mitochondrial populations within the cell.

Universal catastrophe time distributions of dynamically unstable polymers

Dynamic instability-the growth, catastrophe, and shrinkage of quasi-one-dimensional filaments-has been observed in multiple biopolymers. Scientists have long understood the catastrophic cessation of growth and subsequent depolymerization as arising from the interplay of hydrolysis and polymerization at the tip of the polymer. Here we show that for a broad class of catastrophe models, the expected catastrophe time distribution is exponential. We show that the distribution shape is insensitive to noise, but that depletion of monomers from a finite pool can dramatically change the distribution shape by reducing the polymerization rate. We derive a form for this finite-pool catastrophe time distribution and show that finite-pool effects can be important even when the depletion of monomers does not greatly alter the polymerization rate.

Polymerization kinetics of tubulin from mung seedlings modeled as a competition between nucleation and GTP-hydrolysis rates

Microtubules (MTs) form physiologically important cytoskeletal structures that assemble by tubulin polymerization in a nucleation- and GTP-dependent manner. GTP-hydrolysis competes with the addition of monomers, to determine the GTP-cap size, and onset of shrinkage, which alternates with growth. Multiple theoretical models of MT polymerization dynamics have been reconciled to the kinetics of animal brain tubulins, but more recently rapid kinetics seen in Arabidopsis tubulin polymerization suggest the need to sample a wider diversity in tubulin polymerization kinetics and reconcile it to theory. Here, we have isolated tubulin from seedlings of Vigna sp. (mung bean), compared polymerization kinetics to animal brain tubulin and used a computational model to understand the di_erences. We _nd that activity isolated mung tubulin polymerizes in a nucleation-dependent manner, based on turbidimetry, qualitatively similar to brain tubulin, but with a ten-fold smaller critical critical concentration. GTP-dependent polymerization kinetics also appear to be transient, indicative of high rates of GTP-hydrolysis. Computational modeling of tubulin nucleation and vectorial GTP-hydrolysis to examine the e_ect of high nucleation and GTP-hydrolysis rates predicts a dominance of the latter in determining MT lengths and numbers. Microscopy of mung tubulin _laments stabilized by GMPCPP or taxol result in few and short MTs, compared to the many long MTs arising from goat tubulin, qualitatively matching the model predictions. We _nd GTP-hydrolysis outcompetes nucleation rates in determining MT lengths and numbers. This article is protected by copyright. All rights reserved.

Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule

Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.

Microtubule catastrophe under force: mathematical and computational results from a Brownian ratchet model

In the intracellular environment, the intrinsic dynamics of microtubule filaments is often hindered by the presence of barriers of various kind, such as kinetochore complexes and cell cortex, which impact their polymerisation force and dynamical properties such as catastrophe frequency. We present a theoretical study of the effect of a forced barrier, also subjected to thermal noise, on the statistics of catastrophe events in a single microtubule as well as a 'bundle' of two parallel microtubules. For microtubule dynamics, which includes growth, detachment, hydrolysis and the consequent dynamic instability, we employ a one-dimensional discrete stochastic model. The dynamics of the barrier is captured by over-damped Langevin equation, while its interaction with a growing filament is assumed to be hard-core repulsion. A unified treatment of the continuum dynamics of the barrier and the discrete dynamics of the filament is realized using a hybrid Fokker-Planck equation. An explicit mathematical formula for the force-dependent catastrophe frequency of a single microtubule is obtained by solving the above equation, under some assumptions. The prediction agrees well with results of numerical simulations in the appropriate parameter regime. More general situations are studied via numerical simulations. To investigate the extent of 'load-sharing' in a microtubule bundle, and its impact on the frequency of catastrophes, the dynamics of a two-filament bundle is also studied. Here, two parallel, non-interacting microtubules interact with a common, forced barrier. The equations for the two-filament model, when solved using a mean-field assumption, predicts equal sharing of load between the filaments. However, numerical results indicate the existence of a wide spectrum of load-sharing behaviour, which is characterized using a dimensionless parameter.

Regulation of microtubule disassembly by spatially heterogeneous patterns of acetylation

Microtubules (MTs) are bio-polymers, composed of tubulin proteins, involved in several functions such as cell division, transport of cargoes within cells, maintaining cellular structures etc. Their kinetics are often affected by chemical modifications on the filament known as Post Translational Modifications (PTMs). Acetylation is a PTM which occurs on the luminal surface of the MT lattice and has been observed to reduce the lateral interaction between tubulins on adjacent protofilaments. Depending on the properties of the acetylase enzyme αTAT1 and the structural features of MTs, the patterns of acetylation formed on MTs are observed to be quite diverse. In this study, we present a multi-protofilament model with spatially heterogeneous patterns of acetylation, and investigate how the local kinetic differences arising from heterogeneity affect the global kinetics of MT filaments. From the computational study we conclude that a filament with spatially uniform acetylation is least stable against disassembly, while ones with more clustered acetylation patterns may provide better resistance against disassembly. The increase in disassembly times for clustered pattern as compared to uniform pattern can be up to fifty percent for identical amounts of acetylation. Given that acetylated MTs affect several cellular functions as well as diseases such as cancer, our study indicates that spatial patterns of acetylation need to be focused on, apart from the overall amount of acetylation.

Behaviors of individual microtubules and microtubule populations relative to critical concentrations: dynamic instability occurs when critical concentrations are driven apart by nucleotide hydrolysis

The concept of critical concentration (CC) is central to understanding the behavior of microtubules (MTs) and other cytoskeletal polymers. Traditionally, these polymers are understood to have one CC, measured in multiple ways and assumed to be the subunit concentration necessary for polymer assembly. However, this framework does not incorporate dynamic instability (DI), and there is work indicating that MTs have two CCs. We use our previously established simulations to confirm that MTs have (at least) two experimentally relevant CCs and to clarify the behavior of individuals and populations relative to the CCs. At free subunit concentrations above the lower CC (CC_Elongation), growth phases of individual filaments can occur transiently; above the higher CC (CC_NetAssembly), the population’s polymer mass will increase persistently. Our results demonstrate that most experimental CC measurements correspond to CC_NetAssembly, meaning that “typical” DI occurs below the concentration traditionally considered necessary for polymer assembly. We report that [free tubulin] at steady state does not equal CC_NetAssembly, but instead approaches CC_NetAssembly asymptotically as [total tubulin] increases, and depends on the number of stable MT nucleation sites. We show that the degree of separation between CC_Elongation and CC_NetAssembly depends on the rate of nucleotide hydrolysis. This clarified framework helps explain and unify many experimental observations.

Bacterial cell division: modeling FtsZ assembly and force generation from single filament experimental data

The bacterial cytoskeletal protein FtsZ binds and hydrolyzes GTP, self-aggregates into dynamic filaments and guides the assembly of the septal ring on the inner side of the membrane at midcell. This ring constricts the cell during division and is present in most bacteria. Despite exhaustive studies undertaken in the last 25 years after its discovery, we do not yet know the mechanism by which this GTP-dependent self-aggregating protein exerts force on the underlying membrane. This paper reviews recent experiments and theoretical models proposed to explain FtsZ filament dynamic assembly and force generation. It highlights how recent observations of single filaments on reconstituted model systems and computational modeling are contributing to develop new multiscale models that stress the importance of previously overlooked elements as monomer internal flexibility, filament twist and flexible anchoring to the cell membrane. These elements contribute to understand the rich behavior of these GTP consuming dynamic filaments on surfaces. The aim of this review is 2-fold: (1) to summarize recent multiscale models and their implications to understand the molecular mechanism of FtsZ assembly and force generation and (2) to update theoreticians with recent experimental results.

A Growth-Fragmentation Approach for Modeling Microtubule Dynamic Instability

Microtubules (MTs) are protein filaments found in all eukaryotic cells which are crucial for many cellular processes including cell movement, cell differentiation, and cell division. Due to their role in cell division, they are often used as targets for chemotherapy drugs used in cancer treatment. Experimental studies of MT dynamics have played an important role in the development and administration of many novel cancer drugs; however, a complete description of MT dynamics is lacking. Here, we propose a new mathematical model for MT dynamics, that can be used to study the effects of chemotherapy drugs on MT dynamics. Our model consists of a growth-fragmentation equation describing the dynamics of a length distribution of MTs, coupled with two ODEs that describe the dynamics of free GTP- and GDP-tubulin concentrations (the individual dimers that comprise of MTs). Here, we prove the well-posedness of our system and perform a numerical exploration of the influence of certain model parameters on the systems dynamics. In particular, we focus on a qualitative description for how a certain class of destabilizing drugs, the vinca alkaloids, alter MT dynamics. Through variation of certain model parameters which we know are altered by these drugs, we make comparisons between simulation results and what is observed in in vitro studies.

SPR2 protects minus ends to promote severing and reorientation of plant cortical microtubule arrays

The cortical microtubule arrays of higher plants are organized without centrosomes and feature treadmilling polymers that are dynamic at both ends. The control of polymer end stability is fundamental for the assembly and organization of cytoskeletal arrays, yet relatively little is understood about how microtubule minus ends are controlled in acentrosomal microtubule arrays, and no factors have been identified that act at the treadmilling minus ends in higher plants. Here, we identify Arabidopsis thaliana SPIRAL2 (SPR2) as a protein that tracks minus ends and protects them against subunit loss. SPR2 function is required to facilitate the rapid reorientation of plant cortical arrays as stimulated by light perception, a process that is driven by microtubule severing to create a new population of microtubules. Quantitative live-cell imaging and computer simulations reveal that minus protection by SPR2 acts by an unexpected mechanism to promote the lifetime of potential SPR2 severing sites, increasing the likelihood of severing and thus the rapid amplification of the new microtubule array.

Physical principles of filamentous protein self-assembly kinetics

The polymerization of proteins and peptides into filamentous supramolecular structures is an elementary form of self-organization of key importance to the functioning biological systems, as in the case of actin biofilaments that compose the cellular cytoskeleton. Aberrant filamentous protein self-assembly, however, is associated with undesired effects and severe clinical disorders, such as Alzheimer's and Parkinson's diseases, which, at the molecular level, are associated with the formation of certain forms of filamentous protein aggregates known as amyloids. Moreover, due to their unique physicochemical properties, protein filaments are finding extensive applications as biomaterials for nanotechnology. With all these different factors at play, the field of filamentous protein self-assembly has experienced tremendous activity in recent years. A key question in this area has been to elucidate the microscopic mechanisms through which filamentous aggregates emerge from dispersed proteins with the goal of uncovering the underlying physical principles. With the latest developments in the mathematical modeling of protein aggregation kinetics as well as the improvement of the available experimental techniques it is now possible to tackle many of these complex systems and carry out detailed analyses of the underlying microscopic steps involved in protein filament formation. In this paper, we review some classical and modern kinetic theories of protein filament formation, highlighting their use as a general strategy for quantifying the molecular-level mechanisms and transition states involved in these processes.

Signatures of a macroscopic switching transition for a dynamic microtubule

Characterising complex kinetics of non-equilibrium self-assembly of bio-filaments is of general interest. Dynamic instability in microtubules, consisting of successive catastrophes and rescues, is observed to occur as a result of the non-equilibrium conversion of GTP-tubulin to GDP-tubulin. We study this phenomenon using a model for microtubule kinetics with GTP/GDP state-dependent polymerisation, depolymerisation and hydrolysis of subunits. Our results reveal a sharp switch-like transition in the mean velocity of the filaments, from a growth phase to a shrinkage phase, with an associated co-existence of the two phases. This transition is reminiscent of the discontinuous phase transition across the liquid-gas boundary. We probe the extent of discontinuity in the transition quantitatively using characteristic signatures such as bimodality in velocity distribution, variance and Binder cumulant, and also hysteresis behaviour of the system. We further investigate ageing behaviour in catastrophes of the filament, and find that the multi-step nature of catastrophes is intensified in the vicinity of the switching transition. This assumes importance in the context of Microtubule Associated Proteins which have the potential of altering kinetic parameter values.

Steady-state EB cap size fluctuations are determined by stochastic microtubule growth and maturation

Growing microtubules are protected from depolymerization by the presence of a GTP or GDP/Pi cap. End-binding proteins of the EB1 family bind to the stabilizing cap, allowing monitoring of its size in real time. The cap size has been shown to correlate with instantaneous microtubule stability. Here we have quantitatively characterized the properties of cap size fluctuations during steady-state growth and have developed a theory predicting their timescale and amplitude from the kinetics of microtubule growth and cap maturation. In contrast to growth speed fluctuations, cap size fluctuations show a characteristic timescale, which is defined by the lifetime of the cap sites. Growth fluctuations affect the amplitude of cap size fluctuations; however, cap size does not affect growth speed, indicating that microtubules are far from instability during most of their time of growth. Our theory provides the basis for a quantitative understanding of microtubule stability fluctuations during steady-state growth.

Regulation of Microtubule Growth and Catastrophe: Unifying Theory and Experiment

Recent studies have found that microtubule-associated proteins (MAPs) can regulate the dynamical properties of microtubules in unexpected ways. For most MAPs, there is an inverse relationship between their effects on the speed of growth and the frequency of catastrophe, the conversion of a growing microtubule to a shrinking one. Such a negative correlation is predicted by the standard GTP-cap model, which posits that catastrophe is due to loss of a stabilizing cap of GTP-tubulin at the end of a growing microtubule. However, many other MAPs, notably Kinesin-4 and combinations of EB1 with XMAP215, contradict this general rule. In this review, we show that a more nuanced, but still simple, GTP-cap model, can account for the diverse regulatory activities of MAPs.

Microtubule aging probed by microfluidics-assisted tubulin washout

Microtubule aging—the decrease of stability with age—is an interesting but mechanistically not understood property of microtubules. It constrains possible mechanisms of catastrophe induction and is believed to be crucial for length regulation. New in vitro experiments and model fits provide insight into the origin of microtubule aging.

Microtubules switch stochastically between phases of growth and shrinkage. The molecular mechanism responsible for the end of a growth phase, an event called catastrophe, is still not understood. The probability for a catastrophe to occur increases with microtubule age, putting constraints on the possible molecular mechanism of catastrophe induction. Here we used microfluidics-assisted fast tubulin washout experiments to induce microtubule depolymerization in a controlled manner at different times after the start of growth. We found that aging can also be observed in this assay, providing valuable new constraints against which theoretical models of catastrophe induction can be tested. We found that the data can be quantitatively well explained by a simple kinetic threshold model that assumes an age-dependent broadening of the protective cap at the microtubule end as a result of an evolving tapered end structure; this leads to a decrease of the cap density and its stability. This analysis suggests an intuitive picture of the role of morphological changes of the protective cap for the age dependence of microtubule stability.

Dynamics of protein aggregation and oligomer formation governed by secondary nucleation

The formation of aggregates in many protein systems can be significantly accelerated by secondary nucleation, a process where existing assemblies catalyse the nucleation of new species. In particular, secondary nucleation has emerged as a central process controlling the proliferation of many filamentous protein structures, including molecular species related to diseases such as sickle cell anemia and a range of neurodegenerative conditions. Increasing evidence suggests that the physical size of protein filaments plays a key role in determining their potential for deleterious interactions with living cells, with smaller aggregates of misfolded proteins, oligomers, being particularly toxic. It is thus crucial to progress towards an understanding of the factors that control the sizes of protein aggregates. However, the influence of secondary nucleation on the time evolution of aggregate size distributions has been challenging to quantify. This difficulty originates in large part from the fact that secondary nucleation couples the dynamics of species distant in size space. Here, we approach this problem by presenting an analytical treatment of the master equation describing the growth kinetics of linear protein structures proliferating through secondary nucleation and provide closed-form expressions for the temporal evolution of the resulting aggregate size distribution. We show how the availability of analytical solutions for the full filament distribution allows us to identify the key physical parameters that control the sizes of growing protein filaments. Furthermore, we use these results to probe the dynamics of the populations of small oligomeric species as they are formed through secondary nucleation and discuss the implications of our work for understanding the factors that promote or curtail the production of these species with a potentially high deleterious biological activity.

Feedback mechanism for microtubule length regulation by stathmin gradients

We formulate and analyze a theoretical model for the regulation of microtubule (MT) polymerization dynamics by the signaling proteins Rac1 and stathmin. In cells, the MT growth rate is inhibited by cytosolic stathmin, which, in turn, is inactivated by Rac1. Growing MTs activate Rac1 at the cell edge, which closes a positive feedback loop. We investigate both tubulin sequestering and catastrophe promotion as mechanisms for MT growth inhibition by stathmin. For a homogeneous stathmin concentration in the absence of Rac1, we find a switchlike regulation of the MT mean length by stathmin. For constitutively active Rac1 at the cell edge, stathmin is deactivated locally, which establishes a spatial gradient of active stathmin. In this gradient, we find a stationary bimodal MT-length distribution for both mechanisms of MT growth inhibition by stathmin. One subpopulation of the bimodal length distribution can be identified with fast-growing and long pioneering MTs in the region near the cell edge, which have been observed experimentally. The feedback loop is closed through Rac1 activation by MTs. For tubulin sequestering by stathmin, this establishes a bistable switch with two stable states: one stable state corresponds to upregulated MT mean length and bimodal MT length distributions, i.e., pioneering MTs; the other stable state corresponds to an interrupted feedback with short MTs. Stochastic effects as well as external perturbations can trigger switching events. For catastrophe-promoting stathmin, we do not find bistability.

The size of the EB cap determines instantaneous microtubule stability

The function of microtubules relies on their ability to switch between phases of growth and shrinkage. A nucleotide-dependent stabilising cap at microtubule ends is thought to be lost before this switch can occur; however, the nature and size of this protective cap are unknown. Using a microfluidics-assisted multi-colour TIRF microscopy assay with close-to-nm and sub-second precision, we measured the sizes of the stabilizing cap of individual microtubules. We find that the protective caps are formed by the extended binding regions of EB proteins. Cap lengths vary considerably and longer caps are more stable. Nevertheless, the trigger of instability lies in a short region at the end of the cap, as a quantitative model of cap stability demonstrates. Our study establishes the spatial and kinetic characteristics of the protective cap and provides an insight into the molecular mechanism by which its loss leads to the switch from microtubule growth to shrinkage.

DOI:http://dx.doi.org/10.7554/eLife.13470.001

Much like the skeleton supports the human body, a structure called the cytoskeleton provides support and structure to cells. Part of this cytoskeleton is made up of small tubes called microtubules that – unlike bones – can shrink and grow very quickly. This allows the cell to change shape, move and split into two new cells.

Exactly how the microtubules switch between growing and shrinking was not clear. One suggestion is that a protective cap at the end of microtubule allows it to keep growing and prevents it from shrinking. However, the nature and size of this cap have been debated.

Now, Duellberg et al. have measured the caps of microtubules with high precision by combining the techniques of microfluidics, TIRF microscopy and recently developed image analysis tools. This revealed that the cap sizes change, with longer caps being more stable. In addition, proteins called end-binding proteins can destabilize the cap by binding to it. This allows microtubules to switch from a growing to a shrinking state more often.

Future work could now investigate how changes in cap length cause the microtubules to switch from growing to shrinking. It also remains to be seen whether other proteins also influence the cap to control this switching behaviour.

DOI:http://dx.doi.org/10.7554/eLife.13470.002

When cell biology meets theory

Cell biologists now have tools and knowledge to generate useful quantitative data. But how can we make sense of these data, and are we measuring the correct parameters? Moreover, how can we test hypotheses quantitatively? To answer these questions, the theory of physics is required and is essential to the future of quantitative cell biology.

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.

Thermodynamics and kinetics of DNA nanotube polymerization from single-filament measurements

Single-filament measurement of the thermodynamic and kinetic parameters of DNA nanotube assembly supports a polymerization/depolymerization model sharing common features with cytoskeletal polymer models.

DNA nanotubes provide a programmable architecture for molecular self-assembly and can serve as model systems for one-dimensional biomolecular assemblies. While a variety of DNA nanotubes have been synthesized and employed as models for natural biopolymers, an extensive investigation of DNA nanotube kinetics and thermodynamics has been lacking. Using total internal reflection microscopy, DNA nanotube polymerization was monitored in real time at the single filament level over a wide range of free monomer concentrations and temperatures. The measured polymerization rates were subjected to a global nonlinear fit based on polymerization theory in order to simultaneously extract kinetic and thermodynamic parameters. For the DNA nanotubes used in this study, the association rate constant is (5.99 ± 0.15) × 10⁵ M^–1 s^–1, the enthalpy is 87.9 ± 2.0 kcal mol^–1, and the entropy is 0.252 ± 0.006 kcal mol^–1 K^–1. The qualitative and quantitative similarities between the kinetics of DNA nanotubes, actin filaments, and microtubules polymerization highlight the prospect of building complex dynamic systems from DNA molecules inspired by biological architecture.

Force-induced dynamical properties of multiple cytoskeletal filaments are distinct from that of single filaments

How cytoskeletal filaments collectively undergo growth and shrinkage is an intriguing question. Collective properties of multiple bio-filaments (actin or microtubules) undergoing hydrolysis have not been studied extensively earlier within simple theoretical frameworks. In this paper, we study the collective dynamical properties of multiple filaments under force, and demonstrate the distinct properties of a multi-filament system in comparison to a single filament. Comparing stochastic simulation results with recent experimental data, we show that multi-filament collective catastrophes are slower than catastrophes of single filaments. Our study also shows further distinctions as follows: (i) force-dependence of the cap-size distribution of multiple filaments are quantitatively different from that of single filaments, (ii) the diffusion constant associated with the system length fluctuations is distinct for multiple filaments, and (iii) switching dynamics of multiple filaments between capped and uncapped states and the fluctuations therein are also distinct. We build a unified picture by establishing interconnections among all these collective phenomena. Additionally, we show that the collapse times during catastrophes can be sharp indicators of collective stall forces exceeding the additive contributions of single filaments.

Asymptotic solutions of the Oosawa model for the length distribution of biofilaments

Nucleated polymerisation phenomena are general linear growth processes that underlie the formation of a range of biofilaments in nature, including actin and tubulin that are key components of the cellular cytoskeleton. The conventional theoretical framework for describing this process is the Oosawa model that takes into account homogeneous nucleation coupled to linear growth. In his original work, Oosawa provided an analytical solution to the total mass concentration of filaments; the time evolution of the full length distribution has, however, been challenging to access, in large part due to the nonlinear nature of the rate equations inherent in the description of such phenomena and to date analytical solutions for the filament distribution are known only in certain special cases. Here, by exploiting a technique based on the method of matched asymptotics, we present an analytical treatment of the Oosawa model that describes the shape of the length distribution of biofilaments reversibly growing through primary nucleation and filament elongation. Our work highlights the power of matched asymptotics for obtaining closed-form analytical solutions to nonlinear master equations in biophysics and allows us to identify the key time scales that characterize biological polymerization processes.

Brownian dynamics of subunit addition-loss kinetics and thermodynamics in linear polymer self-assembly

The structure and free energy of multistranded linear polymer ends evolves as individual subunits are added and lost. Thus, the energetic state of the polymer end is not constant, as assembly theory has assumed. Here we utilize a Brownian dynamics approach to simulate the addition and loss of individual subunits at the polymer tip. Using the microtubule as a primary example, we examined how the structure of the polymer tip dictates the rate at which units are added to and lost from individual protofilaments. We find that freely diffusing subunits arrive less frequently to lagging protofilaments but bind more efficiently, such that there is no kinetic difference between leading and lagging protofilaments within a tapered tip. However, local structure at the nanoscale has up to an order-of-magnitude effect on the rate of addition. Thus, the kinetic on-rate constant, integrated across the microtubule tip (kon,MT), is an ensemble average of the varying individual protofilament on-rate constants (kon,PF). Our findings have implications for both catastrophe and rescue of the dynamic microtubule end, and provide a subnanoscale framework for understanding the mechanism of action of microtubule-associated proteins and microtubule-directed drugs. Although we utilize the specific example of the microtubule here, the findings are applicable to multistranded polymers generally.

Effects of microtubule mechanics on hydrolysis and catastrophes

We introduce a model for microtubule (MT) mechanics containing lateral bonds between dimers in neighboring protofilaments, bending rigidity of dimers, and repulsive interactions between protofilaments modeling steric constraints to investigate the influence of mechanical forces on hydrolysis and catastrophes. We use the allosteric dimer model, where tubulin dimers are characterized by an equilibrium bending angle, which changes from 0° to 22° by hydrolysis of a dimer. This also affects the lateral interaction and bending energies and, thus, the mechanical equilibrium state of the MT. As hydrolysis gives rise to conformational changes in dimers, mechanical forces also influence the hydrolysis rates by mechanical energy changes modulating the hydrolysis rate. The interaction via the MT mechanics then gives rise to correlation effects in the hydrolysis dynamics, which have not been taken into account before. Assuming a dominant influence of mechanical energies on hydrolysis rates, we investigate the most probable hydrolysis pathways both for vectorial and random hydrolysis. Investigating the stability with respect to lateral bond rupture, we identify initiation configurations for catastrophes along the hydrolysis pathways and values for a lateral bond rupture force. If we allow for rupturing of lateral bonds between dimers in neighboring protofilaments above this threshold force, our model exhibits avalanche-like catastrophe events.

Microtubule catastrophe from protofilament dynamics

The disappearance of the guanosine triphosphate- (GTP) tubulin cap is widely believed to be the forerunner event for the growth-shrinkage transition ("catastrophe") in microtubule filaments in eukaryotic cells. We study a discrete version of a stochastic model of the GTP cap dynamics, originally proposed by Flyvbjerg, Holy, and Leibler [Phys. Rev. Lett. 73, 2372 (1994)]. Our model includes both spontaneous and vectorial hydrolysis, as well as dissociation of a nonhydrolyzed dimer from the filament after incorporation. In the first part of the paper, we apply this model to a single protofilament of a microtubule. A catastrophe transition is defined for each protofilament, similarly to the earlier one-dimensional models, the frequency of occurrence of which is then calculated under various conditions but without explicit assumption of steady-state conditions. Using a perturbative approach, we show that the leading asymptotic behavior of the protofilament catastrophe in the limit of large growth velocities is remarkably similar across different models. In the second part of the paper, we extend our analysis to the entire filament by making a conjecture that a minimum number of such transitions are required to occur for the onset of microtubule catastrophe. The frequency of microtubule catastrophe is then determined using numerical simulations and compared with analytical and semianalytical estimates made under steady-state and quasi-steady-state assumptions, respectively, for the protofilament dynamics. A few relevant experimental results are analyzed in detail and compared with predictions from the model. Our results indicate that loss of GTP cap in two to three protofilaments is necessary to trigger catastrophe in a microtubule.

Microtubule dynamic instability: a new model with coupled GTP hydrolysis and multistep catastrophe

A key question in understanding microtubule dynamics is how GTP hydrolysis leads to catastrophe, the switch from slow growth to rapid shrinkage. We first provide a review of the experimental and modeling literature, and then present a new model of microtubule dynamics. We demonstrate that vectorial, random, and coupled hydrolysis mechanisms are not consistent with the dependence of catastrophe on tubulin concentration and show that, although single-protofilament models can explain many features of dynamics, they do not describe catastrophe as a multistep process. Finally, we present a new combined (coupled plus random hydrolysis) multiple-protofilament model that is a simple, analytically solvable generalization of a single-protofilament model. This model accounts for the observed lifetimes of growing microtubules, the delay to catastrophe following dilution and describes catastrophe as a multistep process.

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.

Cooperative dynamics of microtubule ensembles: polymerization forces and rescue-induced oscillations

We investigate the cooperative dynamics of an ensemble of N microtubules growing against an elastic barrier. Microtubules undergo so-called catastrophes, which are abrupt stochastic transitions from a growing to a shrinking state, and rescues, which are transitions back to the growing state. Microtubules can exert pushing or polymerization forces on an obstacle, such as an elastic barrier, if the growing end is in contact with the obstacle. We use dynamical mean-field theory and stochastic simulations to analyze a model where each microtubule undergoes catastrophes and rescues and where microtubules interact by force sharing. For zero rescue rate, cooperative growth terminates in a collective catastrophe. The maximal polymerization force before catastrophes grows linearly with N for small N or a stiff elastic barrier, in agreement with available experimental results, whereas it crosses over to a logarithmic dependence for larger N or a soft elastic barrier. For a nonzero rescue rate and a soft elastic barrier, the dynamics becomes oscillatory with both collective catastrophe and rescue events, which are part of a robust limit cycle. Both the average and maximal polymerization forces then grow linearly with N, and we investigate their dependence on tubulin on-rates and rescue rates, which can be involved in cellular regulation mechanisms. We further investigate the robustness of the collective catastrophe and rescue oscillations with respect to different catastrophe models.

Dynamics and length distribution of microtubules under force and confinement

We investigate the microtubule polymerization dynamics with catastrophe and rescue events for three different confinement scenarios, which mimic typical cellular environments: (i) The microtubule is confined by rigid and fixed walls, (ii) it grows under constant force, and (iii) it grows against an elastic obstacle with a linearly increasing force. We use realistic catastrophe models and analyze the microtubule dynamics, the resulting microtubule length distributions, and force generation by stochastic and mean field calculations; in addition, we perform stochastic simulations. Freely growing microtubules exhibit a phase of bounded growth with finite microtubule length and a phase of unbounded growth. The main results for the three confinement scenarios are as follows: (i) In confinement by fixed rigid walls, we find exponentially decreasing or increasing stationary microtubule length distributions instead of bounded or unbounded phases, respectively. We introduce a realistic model for wall-induced catastrophes and investigate the behavior of the average length as a function of microtubule growth parameters. (ii) Under a constant force, the boundary between bounded and unbounded growth is shifted to higher tubulin concentrations and rescue rates. The critical force f(c) for the transition from unbounded to bounded growth increases logarithmically with tubulin concentration and the rescue rate, and it is smaller than the stall force. (iii) For microtubule growth against an elastic obstacle, the microtubule length and polymerization force can be regulated by microtubule growth parameters. For zero rescue rate, we find that the average polymerization force depends logarithmically on the tubulin concentration and is always smaller than the stall force in the absence of catastrophes and rescues. For a nonzero rescue rate, we find a sharply peaked steady-state length distribution, which is tightly controlled by microtubule growth parameters. The corresponding average microtubule length self-organizes such that the average polymerization force equals the critical force f(c) for the transition from unbounded to bounded growth. We also investigate the force dynamics if growth parameters are perturbed in dilution experiments. Finally, we show the robustness of our results against changes of catastrophe models and load distribution factors.

On the nature and shape of tubulin trails: implications on microtubule self-organization

Microtubules, major elements of the cell skeleton are, most of the time, well organized in vivo, but they can also show self-organizing behaviors in time and/or space in purified solutions in vitro. Theoretical studies and models based on the concepts of collective dynamics in complex systems, reaction-diffusion processes and emergent phenomena were proposed to explain some of these behaviors. In the particular case of microtubule spatial self-organization, it has been advanced that microtubules could behave like ants, self-organizing by 'talking to each other' by way of hypothetic (because never observed) concentrated chemical trails of tubulin that are expected to be released by their disassembling ends. Deterministic models based on this idea yielded indeed like-looking spatio-temporal self-organizing behaviors. Nevertheless the question remains of whether microscopic tubulin trails produced by individual or bundles of several microtubules are intense enough to allow microtubule self-organization at a macroscopic level. In the present work, by simulating the diffusion of tubulin in microtubule solutions at the microscopic scale, we measure the shape and intensity of tubulin trails and discuss about the assumption of microtubule self-organization due to the production of chemical trails by disassembling microtubules. We show that the tubulin trails produced by individual microtubules or small microtubule arrays are very weak and not elongated even at very high reactive rates. Although the variations of concentration due to such trails are not significant compared to natural fluctuations of the concentration of tubuline in the chemical environment, the study shows that heterogeneities of biochemical composition can form due to microtubule disassembly. They could become significant when produced by numerous microtubule ends located in the same place. Their possible formation could play a role in certain conditions of reaction. In particular, it gives a mesoscopic basis to explain the collective dynamics observed in excitable microtubule solutions showing the propagation of concentration waves of microtubules at the millimeter scale, although we doubt that individual microtubules or bundles can behave like molecular ants.

Random hydrolysis controls the dynamic instability of microtubules

Uncovering mechanisms that control the dynamics of microtubules is fundamental for our understanding of multiple cellular processes such as chromosome separation and cell motility. Building on previous theoretical work on the dynamic instability of microtubules, we propose here a stochastic model that includes all relevant biochemical processes that affect the dynamics of microtubule plus-end, namely, the binding of GTP-bound monomers, unbinding of GTP- and GDP-bound monomers, and hydrolysis of GTP monomers. The inclusion of dissociation processes, present in our approach but absent from many previous studies, is essential to guarantee the thermodynamic consistency of the model. Our theoretical method allows us to compute all dynamic properties of microtubules explicitly. Using experimentally determined rates, it is found that the cap size is ∼3.6 layers, an estimate that is compatible with several experimental observations. In the end, our model provides a comprehensive description of the dynamic instability of microtubules that includes not only the statistics of catastrophes but also the statistics of rescues.

Modeling the effects of drug binding on the dynamic instability of microtubules

We propose a stochastic model that accounts for the growth, catastrophe and rescue processes of steady-state microtubules assembled from MAP-free tubulin in the possible presence of a microtubule-associated drug. As an example of the latter, we both experimentally and theoretically study the perturbation of microtubule dynamic instability by S-methyl-D-DM1, a synthetic derivative of the microtubule-targeted agent maytansine and a potential anticancer agent. Our model predicts that among the drugs that act locally at the microtubule tip, primary inhibition of the loss of GDP tubulin results in stronger damping of microtubule dynamics than inhibition of GTP tubulin addition. On the other hand, drugs whose action occurs in the interior of the microtubule need to be present in much higher concentrations to have visible effects.

Macroscopic simulations of microtubule dynamics predict two steady-state processes governing array morphology

Microtubule polymers typically function through their collective organization into a patterned array. The formation of the pattern, whether it is a relatively simple astral array or a highly complex mitotic spindle, relies on controlled microtubule nucleation and the basal dynamics parameters governing polymer growth and shortening. We have investigated the interaction between the microtubule nucleation and dynamics parameters, using macroscopic Monte Carlo simulations, to determine how these parameters contribute to the underlying microtubule array morphology (i.e. polymer density and length distribution). In addition to the well-characterized steady state achieved between free tubulin subunits and microtubule polymer, we propose that microtubule nucleation and extinction constitute a second, interdependent steady state process. Our simulation studies show that the magnitude of both nucleation and extinction additively impacts the final steady state free subunit concentration. We systematically varied individual microtubule dynamics parameters to survey the effects on array morphology and find specific sensitivity to perturbations of catastrophe frequency. Altering the cellular context for the microtubule array, we find that nucleation template number plays a defining role in shaping the microtubule length distribution and polymer density.

Prolonging assembly through dissociation: a self-assembly paradigm in microtubules

We study a one-dimensional model of microtubule assembly and disassembly in which GTP bound to tubulins within the microtubule undergoes stochastic hydrolysis. In contrast to models that consider only a cap of GTP-bound tubulin, stochastic hydrolysis allows GTP-bound tubulin remnants to exist within the microtubule. We find that these buried GTP remnants enable an alternative mechanism of recovery from shrinkage and enhances fluctuations of filament lengths. Under conditions for which this alternative mechanism dominates, an increasing depolymerization rate leads to a decrease in dissociation rate and thus a net increase in assembly.

A first-passage-time theory for search and capture of chromosomes by microtubules in mitosis

The mitotic spindle is an important intermediate structure in eukaryotic cell division, in which each of a pair of duplicated chromosomes is attached through microtubules to centrosomal bodies located close to the two poles of the dividing cell. Several mechanisms are at work toward the formation of the spindle, one of which is the 'capture' of chromosome pairs, held together by kinetochores, by randomly searching microtubules. Although the entire cell cycle can be up to 24 hours long, the mitotic phase typically takes only less than an hour. How does the cell keep the duration of mitosis within this limit? Previous theoretical studies have suggested that the chromosome search and capture is optimized by tuning the microtubule dynamic parameters to minimize the search time. In this paper, we examine this conjecture. We compute the mean search time for a single target by microtubules from a single nucleating site, using a systematic and rigorous theoretical approach, for arbitrary kinetic parameters. The result is extended to multiple targets and nucleating sites by physical arguments. Estimates of mitotic time scales are then obtained for different cells using experimental data. In yeast and mammalian cells, the observed changes in microtubule kinetics between interphase and mitosis are beneficial in reducing the search time. In Xenopus extracts, by contrast, the opposite effect is observed, in agreement with the current understanding that large cells use additional mechanisms to regulate the duration of the mitotic phase.

Role of ATP-hydrolysis in the dynamics of a single actin filament

We study the stochastic dynamics of growth and shrinkage of single actin filaments taking into account insertion, removal, and ATP hydrolysis of subunits either according to the vectorial mechanism or to the random mechanism. In a previous work, we developed a model for a single actin or microtubule filament where hydrolysis occurred according to the vectorial mechanism: the filament could grow only from one end, and was in contact with a reservoir of monomers. Here we extend this approach in two ways--by including the dynamics of both ends and by comparing two possible mechanisms of ATP hydrolysis. Our emphasis is mainly on two possible limiting models for the mechanism of hydrolysis within a single filament, namely the vectorial or the random model. We propose a set of experiments to test the nature of the precise mechanism of hydrolysis within actin filaments.

Suprastructures and dynamic properties of Mycobacterium tuberculosis FtsZ

Tuberculosis causes the most death in humans by any bacterium. Drug targeting of bacterial cytoskeletal proteins requires detailed knowledge of the various filamentous suprastructures and dynamic properties. Here, we have investigated by high resolution electron microscopy the assembly of cell division protein and microtubule homolog FtsZ from Mycobacterium tuberculosis (MtbFtsZ) in vitro in the presence of various monovalent salts, crowding agents and polycations. Supramolecular structures, including two-dimensional rings, three-dimensional toroids, and multistranded helices formed in the presence of molecular crowding, were similar to those observed by fluorescence microscopy in bacteria in vivo. Dynamic properties of MtbFtsZ filaments were visualized by light scattering and real time total internal reflection fluorescence microscopy. Interestingly, MtbFtsZ revealed a form of dynamic instability at steady state. Cation-induced condensation phenomena of bacterial cytomotive polymers have not been investigated in any detail, although it is known that many bacteria can contain high amounts of polycations, which may modulate the prokaryotic cytoskeleton. We find that above a threshold concentration of polycations which varied with the valence of the cation, ionic strength, and pH, MtbFtsZ mainly formed sheets. The general features of these cation-induced condensation phenomena could be explained in the framework of the Manning condensation theory. Chirality and packing defects limited the dimensions of sheets and toroids at steady state as predicted by theoretical models. In first approximation simple physical principles seem to govern the formation of MtbFtsZ suprastructures.

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.

FtsZ filament dynamics at steady state: subunit exchange with and without nucleotide hydrolysis

We have measured three aspects of FtsZ filament dynamics at steady state: rates of GTP hydrolysis, subunit exchange between protofilaments, and disassembly induced by dilution or excess GDP. All three reactions were slowed with an increase in the potassium concentration from 100 to 500 mM, via replacement of potassium with rubidium, or with an increase in the magnesium concentration from 5 to 20 mM. Electron microscopy showed that the polymers assembled under the conditions of fastest assembly were predominantly short, one-stranded protofilaments, whereas under conditions of slower dynamics, the protofilaments tended to associate into long, thin bundles. We suggest that exchange of subunits between protofilaments at steady state involves two separate mechanisms: (1) fragmentation or dissociation of subunits from protofilament ends following GTP hydrolysis and (2) reversible association and dissociation of subunits from protofilament ends independent of hydrolysis. Exchange of nucleotides on these recycling subunits could give the appearance of exchange directly into the polymer. Several of our observations suggest that exchange of nucleotide can take place on these recycling subunits, but not directly into the FtsZ polymer. Annealing of protofilaments was demonstrated for the L68W mutant in EDTA buffer but not in Mg buffer, where rapid cycling of subunits may obscure the effect of annealing. We also reinvestigated the nucleotide composition of FtsZ polymers at steady state. We found that the GDP:GTP ratio was 50:50 for concentrations of GTP >100 microM, significantly higher than the 20:80 ratio previously reported at 20 microM GTP.

Nonequilibrium self-assembly of a filament coupled to ATP/GTP hydrolysis

We study the stochastic dynamics of growth and shrinkage of single actin filaments or microtubules taking into account insertion, removal, and ATP/GTP hydrolysis of subunits. The resulting phase diagram contains three different phases: two phases of unbounded growth: a rapidly growing phase and an intermediate phase, and one bounded growth phase. We analyze all these phases, with an emphasis on the bounded growth phase. We also discuss how hydrolysis affects force-velocity curves. The bounded growth phase shows features of dynamic instability, which we characterize in terms of the time needed for the ATP/GTP cap to disappear as well as the time needed for the filament to reach a length of zero (i.e., to collapse) for the first time. We obtain exact expressions for all these quantities, which we test using Monte Carlo simulations.