Continuum descriptions of cytoskeletal dynamics
© Karsten; licensee BioMed Central Ltd. 2013
Published: 10 December 2013
This tutorial presents an introduction into continuum descriptions of cytoskeletal dynamics. In contrast to discrete models in which each molecule keeps its identity, such descriptions are given in terms of averaged quantities per unit volume like the number density of a certain molecule. Starting with a discrete description for the assembly dynamics of cytoskeletal filaments, we derive the continuity equation, which serves as the basis of many continuum theories. We illustrate the use of this approach with an investigation of spontaneous cytoskeletal polymerization waves. Such waves have by now been observed in various cell types and might help to orchestrate cytoskeletal dynamics during cell spreading and locomotion. Our analysis shows how processes at the scale of single molecules, namely, the nucleation of new filaments and filament treadmilling, can lead to the spontaneous appearance of coherent traveling waves on scales spanning many filament lengths. For readers less familiar with calculus, we include an informal introduction to the Taylor expansion.
With the advancements in microscopy techniques it has become increasingly clear that a true understanding of many cellular phenomena requires to take spatial aspects into account. This is clearly the case for the cytoskeleton as illustrated by the changes in the microtubule network during cell division or the reorganization of the actin meshwork during cell locomotion. To reach a quantitative understanding of the mechanisms underlying these processes, concepts and methods from physics can be extremely valuable. These methods include notably the theoretical analysis of cellular systems. Most biological and medical curricula today lack, unfortunately, a thorough introduction into mathematical and physical tools, which often leads to discomfort on the part of biologists and physicians, when confronted with the results of a theoretical study. This holds in particular for continuum descriptions. This tutorial is meant to familiarize life scientists with the basic ideas underlying this approach.
Continuum theories have had a great success in describing static and dynamic phenomena for large classes of matter. Well-known examples range from simple fluids  and elastic materials  to fluid membranes and vesicles . Less conventional examples include the cytoskeleton  and even flocks of birds . Especially for flocks of birds, one might at first sight be inclined to use a discrete approach instead, in which the position and behavior of each individual bird is considered. Such models indeed exist and have yielded valuable insights, see for example . Similarly, simulation tools like cytosim allow the user to study cytoskeletal dynamics, while keeping track of each individual cytoskeletal filament, of each motor molecule, and of any other cytoskeletal protein possibly present, for example, the passive cross-linker α-actinin . There are, however, two major problems associated with this approach. One is most easily revealed by considering water. In a discrete description, we would need to track the position of each individual water molecule, its momentum, and its interaction with all surrounding water molecules. Given the astronomic number of water molecules even in small volumes, this is a daunting task. Furthermore, one needs to provide details of the interactions between the individual entities and it is a priori not clear, which of these details will be important for the large scale dynamics of a system.
In contrast, in a continuum theory, the ultimately particulate nature of a material or system is neglected. Intuitively, this approximation is valid for structures on length scales that are very much larger than the discrete constituents of the material under investigation. Maybe somewhat surprisingly, though, they often describe a material's behavior correctly even down to length scales that approach the constituents' sizes. The Navier-Stokes equation, for example, describes the flow of water in micro- and to some extent even in nanofluidic devices . From a more technical point of view, continuum theories often have the advantage over more microscopic descriptions to provide a comprehensive account of possible material behavior. Furthermore, many of the details of the interactions between the constituents are lumped together into a comparatively small set of parameters. In this tutorial, we show how continuum theories can be developed and used in a cellular context to describe cytoskeletal dynamics.
Our aim is not to enable you to formulate and solve your own continuum theories. Still, you should get a feeling for the method and for what you can and cannot expect from this approach. This tutorial necessarily contains a certain amount of mathematics. To make the tutorial accessible to an audience not using calculus on an everyday basis, we have tried to make all steps in the calculations explicit. You should have a pencil and a sheet of paper ready to follow them in detail. In addition, some calculations are left for the reader. To perform them you do not need much beyond a sincere interest in these techniques and some basic knowledge of calculus, notably to take derivatives. To keep the mathematics tractable, we have also simplified certain biological processes. The emphasis of this tutorial is really on mathematical and physical concepts rather than on a detailed analysis of cellular processes.
The tutorial is organized as follows. We start with a discrete description of the length dynamics of cytoskeletal filaments from which we derive a continuum theory. This approach will lead us to the continuity equation. We will use this equation to develop a continuum description of treadmilling filaments interacting with nucleation promoting factors. We will use this formalism to study a possible mechanism underlying spontaneous polymerization waves. Two excursions aim at recalling the Taylor expansion and give a general account of the continuity equation.
Mathematical prelude: The Taylor expansion
A mathematical tool that we will make frequent use of when going to the continuum limit is the so-called Taylor expansion of a (real-valued) function. The aim here is not to give a mathematically rigorous account of the Taylor expansion, but rather to introduce the basic idea. Consequently, we assume that the functions we will deal with are all sufficiently well-behaved such that all operations we perform on them are allowed.
Here, we have suppressed the dependence on x0 for the ease of notation.
we see that the derivative of the approximation for f equals the approximation for f'. The degree of the polynomial used to approximate a function f is called the order of the approximation.
Your turn: Calculate the Taylor expansion of the exponential function exp(x) with respect to x0 = 0 up to second order.
Obviously, this scheme can be iterated further to get better and better approximations of f in the vicinity of x0. However, we will content ourselves in the following with approximations up to second order at most.
Length-dynamics of active polar filaments - part I
As a first example for a continuum description of cytoskeletal dynamics, we will consider the assembly and disassembly of cytoskeletal filaments. Actin filaments and microtubules are linear aggregates of non-covalently bound subunits, namely G-actin and tubulin dimers, respectively. G-actin is a molecule of about 42 kDa with linear extensions of about 6.7 nm and 4 nm, whereas a tubulin dimer has a molecular weight of about 110 kDa and is about 8 nm long and 4 nm in diameter. The length distributions of F-actin and of microtubules in cells is currently largely unknown, but typical values have been estimated to be around 500 nm up to a few microns for F-actin and up to several tens of microns for microtubules. These subunits are structurally polar resulting in polar filaments with different subunit attachment and removal rates at the two ends. For further reference, we will call the fast growing end the plus-end and the other end the minus-end. In addition, G-actin and tubulin can bind nucleotides and their rates of attachment to and removal from a filament depend on the state of the nucleotide bound, that is, whether it is a nucleoside triphosphate (NTP) or a nucleoside diphosphate (NDP). Through hydrolysis of an NTP energy can be introduced into the system on a molecular level. This is the defining property of an active material.
The structural polarity together with the activity makes the assembly and disassembly dynamics of cytoskeletal filaments more interesting than that of most polymers that are usually studied. In particular, it can lead to a phenomenon called treadmilling, where the plus-end grows on average, while the minus-end shrinks on average [9, 10]. This is the situation we will consider in the following without further questioning the origin of this phenomenon. Also, we will consider generic aspects of treadmilling filaments rather than specifically actin filaments or microtubules, because here we are mostly interested in introducing the concept of continuum theories.
where k a and k d are the rates of subunit attachment at the plus- and of subunit removal at the minus-end, respectively. These rates depend strongly on the environmental conditions and on the state of the nucleoside bound, but typical values are about 10s-1μ M-1 for the attachment rate of ATP-G-actin at the plus-end of an actin filament and 0.25s-1 for the removal rate of ADP-G-actin at the minus-end . Note, however, that the effective rates of subunit addition at the plus- and subunit removal at the minus-end result from complex dynamic processes inside the filament, such that they can strongly deviate from these values. For the time being, we will consider them to be constant, which amounts to assuming a constant reservoir of filament subunits. We will discuss a possible experimental realization below.
While there is in principle nothing wrong with this discrete description, it turns out that a continuum description is easier to analyze. The idea is to neglect the discrete nature of the subunits and consider a continuous length distribution. Indeed, if you are interested in the behavior of the length distribution on scales that are larger than many subunit sizes, many details are irrelevant.
Here, v a = k a δ and v d = k d δ are the assembly and disassembly velocities. Furthermore, in this expression the derivative with respect to time, d/dt has been replaced by a partial derivative ∂/∂t. The reason is that the density c depends on two continuous variables, ℓ and t, and taking the partial derivative indicates that we want to calculate the derivative with respect to t for fixed ℓ. Analogously, ∂/∂ℓ indicates that we take the derivative with respect to ℓ for fixed t.
Interlude: The continuity equation
where j(ℓ, t) = (v a - v d )c(ℓ, t). Such an equation is known as a continuity equation. It generally describes the evolution of a conserved quantity with j being the associated current. We will now see, how such an equation can be obtained on general grounds.
where q(x) = limΔx→ 0Q(x)/ Δx is the (line) density associated with . This is the one-dimensional continuity equation and has the same form as Eq. (9).
Note, that here q is now a surface density.
Your turn: Derive Equation (14). How does it look in three spatial dimensions?
which is the most general form of the continuity equation.
Example: A reaction-diffusion system
The diffusion constants D a and D b are positive and the minus-signs indicate that the flux is directed from high to low concentrations. Since we assume that species A is only generated and destroyed by conversion into B and vice versa, the source terms in Eqs. (16) and (17) have to have opposite signs. Their functional form depends on the molecular details of the conversion processes. If, for example, conversions occur spontaneously at a rate ω a from A to B and at rate ω b from B to A, then one would write S = -ω a a + ω b b.
Here, Δ is the so-called Laplace operator. In absence of the source terms, each of the two equations is the common diffusion equation. In presence of the source terms, the two equations define a simple reaction-diffusion system as introduced by A. Turing .
To complete the description of the evolution of in a finite domain one has to specify so-called boundary conditions. They determine the behavior of the current and/or the density at the domain boundaries. For example, if you consider particles in a vessel with impenetrable walls, then the current through these walls must vanish, j⊥ = 0, where j⊥ denotes the component of the current perpendicular to the wall.
Length-dynamics of active polar filaments - part II
Let us return to the length distribution of treadmilling filaments and focus on the case that v a > v d such that filaments grow on average. If v a and v d are constants, then filament growth would be unbounded. In a real system, though, the pool of subunits is exhausted at some point, which leads to a stop of filament growth. We will consider here a formally simpler alternative, namely that filaments vanish at a constant rate ν d independent of their length. To our knowledge, there is no experimental evidence for such disassembly of neither actin filaments nor microtubules. However, a more realistic description would at this point only distract from our main goal which is to introduce the idea of continuum theories. One could, however, realize such a dynamics in vitro by removing in regular intervals let's say half of the solution containing the filaments and subsequently adding new buffer containing filament subunits. We neglect possible breaking of filaments during this process as well as filament annealing.
Since partitioning and distributing filaments results in diluting them, we also need to provide a source of new filaments. At physiological monomer concentrations, the rate at which new actin filaments form spontaneously is negligible. Instead cells rely on proteins that assist the formation of new filaments, so-called nucleation promoting factors (NPFs). Common examples of such proteins are the Arp2/3 complex, which when itself bound to an existing filaments generates a new filament that has its pointed end connected to the complex, or members of the formin family, which stay attached to the barbed end of a newly generated filament and assist further elongation. In the following, we will neglect these aspects of NPFs and just consider them to generate new filaments.
where λ = (v a - v d )/ν d is a characteristic length.
In the present example, we have derived the continuum theory from an underlying discrete microscopic picture. In general, different microscopic pictures will lead to the same continuum theory. An alternative possibility is to start directly with the continuity equation. In that case, the definition of the current and the source terms reflect essential features of the physics underlying the molecular processes of interest. We will give now an example of this approach and treat the case of spatially heterogeneous distributions of treadmilling filaments interacting with NPFs.
Theoretical description of a network of filaments
where and are the respective filament fluxes in the x-direction. The source and sink terms S+ and S - account for the generation of new filaments and their disappearance.
where v is the treadmilling velocity. As in the previous section, we also have to introduce boundary conditions. We will use in the following periodic boundary conditions, that is, for a system of size L, the densities and currents are L-periodic functions, for example, j+(x + L) = j+(x). We use these boundary conditions mostly for computational convenience, but they have also some biological relevance, for example, in the case of contractile actin rings. Let us recall, however, that in this tutorial, we do not want to give a detailed account of biological processes, but rather focus on formal aspects. For a more realistic application of a continuum description of cytoskeletal waves, see for example Refs. [15–18].
Here, ν is the nucleation rate and ε a numerical factor that accounts for the weaker nucleation activity of cytosolic NPFs. Finally, similarly to the previous section, ν d is the rate of filament disassembly, see Figure 4b.
Also the dynamics of the NPFs is governed by continuity equations. The transport of NPFs is purely diffusive. This holds also for bound NPFs as they stay attached to a specific subunit of a filament and as treadmilling does not move filament subunits. Note, however, that the system behavior does not change qualitatively if NPFs are transported with the filaments, as could be the case in presence of molecular motors. For bound NPFs we assume a constant unbinding rate ω d , see Figure 4e. The attachment dynamics of NPFs to filaments is more involved. First of all there is a base rate ω a per filament at which cytosolic NPFs bind. In addition, we assume a cooperative effect among bound NPFs, namely, that bound NPFs favor the binding of additional NPFs, see Figure 4d. We refrain from suggesting a possible molecular picture underlying this effect and simply note that some cooperative effect is necessary to spontaneously generate waves. Altogether, we arrive at the following dynamic equations:
In this expression, Dc is the diffusion constant for cytosolic NPFs and ω1 is a measure of the cooperative effect on NPF attachment. We have assumed that NPFs are neither generated nor do they degrade.
Your turn: Show that the total number of NPFs is conserved.
This completes our description of treadmilling filaments. There are now several ways to analyze the above equations. An analytic solution of such coupled non-linear partial differential equations is usually not possible and one thus relies on numerical solutions. Before describing one method to numerically solve the equations, let us briefly comment on a powerful method to get an overview about the parameter values for which we might expect interesting behavior.
Linear stability analysis
This method starts by realizing that often, there is one rather simple solution to a continuum theory. Indeed, as is the case for the above equations, there is a spatially homogenous steady state, which is obtained by setting all derivatives equal to zero. This leaves us with an algebraic set of equations, which can be solved numerically. The corresponding homogenous densities are denoted by , , , , and nc,0
Your turn: Derive the equations for the homogenous steady state of Eqs. (30), (31), and (36)-(38).
A real system is permanently subject to fluctuations. What happens to the spatially homogenous state if it is slightly perturbed? To study this question, we can write the distributions in the form and analogously for c - , n+, n - , and nc. Inserting these expressions into the dynamic equations and keeping only terms in linear order in the perturbation, we arrive at a set of linear differential equations.
Your turn: Derive the differential equations up to linear order in the perturbations δc+, δc - , δn+, δn - , and δni. What happens to the zeroth order terms?
The linearized equations are easier to solve than the nonlinear dynamic equations, sometimes even analytically. We do not want to do this explicitly here, as it requires some advanced mathematical techniques. Suffice it here to state that the linear stability analysis not only allows to identify interesting regions in parameter space, but often even yields some properties of the new state. In particular, it can indicate, whether the states that the system assumes under conditions of an unstable homogenous state is stationary or oscillatory. However, only a nonlinear analysis can give the real answer. To this end one usually solves the dynamic equations numerically.
There are many sophisticated methods to numerically solve non-linear partial differential equations. We will briefly introduce one very simple method that makes direct contact to the considerations in the mathematical interlude. While it may be slow, it works in a remarkably large number of cases.
In this expression j i denotes the current from site i - 1 to site i. Since c i Δx is number of filaments in the bin i, this equation is the formal analog of Figure 2a.
Your turn: What is the condition on Δx and Δt resulting from this current?
In this tutorial, we have given an introduction into continuum theories as a tool for analyzing cytoskeletal processes. We gave one example of a description that was derived from a discrete (microscopic) picture, and one example that was developed entirely within the frame of the continuum approach and was based on the continuity equation. Let us mention, that not all relevant quantities obey the continuity equation. If you consider, for example, a cytoskeletal network in two or three spatial dimensions, then one might need to take into account that, on large scales, the network can present orientational order. Since the filament orientation is not a conserved quantity, its dynamics is not given by the continuity equation. To derive dynamic equations for such fields one either relies on microscopic theories or on symmetry considerations.
In the examples we discussed, we assumed a deterministic evolution of the system. In reality, the molecular processes underlying cellular processes are stochastic. In general, it is not an easy task to obtain the correct order of magnitude of stochastic effects. As long as one is only interested in qualitative behavior, it is often sufficient to capture molecular noise by effective diffusion terms.
Probably a more serious point is that we neglected in our examples mechanical forces altogether. However, in cells and often also in reconstituted systems forces play an important role. These can be incorporated into continuum theories by accounting also for momentum conservation in addition to matter conservation. Sources and sinks are in this case given by external forces, while the momentum flux density is given by the stress tensor. In addition to viscous or visco-elastic properties, the stress tensor also accounts for the stress generated by active processes . In contrast, energy conservation is usually not an issue as cellular systems are embedded in a heat reservoir which maintains the system at constant temperature.
In addition to the formal aspects of the descriptions developed above, let us comment on the biological significance of the processes considered. There is at least indirect experimental evidence for a feedback between the activity of nucleation promoting factors and actin filaments in neutrophils , although the molecular players involved are probably not completely known and mechanisms remain to be identified. As mentioned already above, the assembly dynamics of cytoskeletal filaments was also greatly simplified to not distract the reader from the focus of this tutorial. Indeed, neither the assembly nor the disassembly rates of filaments are expected to remain constant in time as typically they depend effectively on the filament length . The assumption of constant assembly and disassembly rates prompted us to also assume a rate at which filaments spontaneously disintegrate. Although this is again unlikely to be the case in living cells, the qualitative behavior of the artificial system is comparable to specific situations in more complicated scenarios . A thorough discussion of realistic applications of continuum theories to cytoskeletal or other-cell-biological systems goes beyond the aim of the present tutorial.
In summary, continuum descriptions provide powerful tools to analyze cellular processes. In a number of cases, they have provided valuable information about biological and, in particular, cellular structures and processes. It will be exciting to see, whether they can also play a role in designing new therapeutic approaches.
I thank A. Dreher for carefully reading the manuscript and him as well as K. Doubrovinski, C. Erlenkämper, and D. Johann for countless discussions on the length dynamics of filaments and on the emergence of spontaneous polymerization waves.
This article has been published as part of Journal of Nanobiotechnology Volume 11 Supplement 1, 2013: Nanophysics for Health. The full contents of the supplement are available online at http://www.jnanobiotechnology.com/supplements/11/S1. Publication charges for this tutorial were funded by the CNRS School "Nanophysics for Health", 5 - 9 November 2012, Mittelwhir, France
- Landau LD, Lifshitz EM: Fluid Mechanics. Course of Theoretical Physics, Butterworth Heinemann, 6: 1987.Google Scholar
- Landau LD, Lifshitz EM: Theory of Elasticity. Course of Theoretical Physics, Butterworth Heinemann, 7: 1987.Google Scholar
- Seifert U: Configurations of fluid membranes and vesicles. Advances in Physics. 1997, 46: 13-137. 10.1080/00018739700101488.View ArticleGoogle Scholar
- Jülicher F, Kruse K, Prost J, Joanny JF: Active behavior of the cytoskeleton. Physics Reports. 2007, 449: 3-28. 10.1016/j.physrep.2007.02.018.View ArticleGoogle Scholar
- Toner J, Tu Y, Ramaswamy S: Hydrodynamics and phases of flocks. Annals of Physics. 2005, 318: 170-244. 10.1016/j.aop.2005.04.011.View ArticleGoogle Scholar
- Vicsek T, Czirók A, Ben-Jacob E, Cohen I, Shochet O: Novel type of phase transition in a system of self-driven particles. Physical Review Letters. 1995, 75: 1226-1229. 10.1103/PhysRevLett.75.1226.View ArticleGoogle Scholar
- Nedelec F, Foethke D: Collective Langevin dynamics of flexible cytoskeletal fibers. New Journal of Physics. 2007, 9: 427-427. 10.1088/1367-2630/9/11/427.View ArticleGoogle Scholar
- Bonn D, Eggers J, Indekeu J, Meunier J, Rolley E: Wetting and spreading. Reviews of Modern Physics. 2009, 81: 739-805. 10.1103/RevModPhys.81.739.View ArticleGoogle Scholar
- Wegner A: Head to tail polymerization of actin. Journal of Molecular Biology. 1976, 108: 139-150. 10.1016/S0022-2836(76)80100-3.View ArticleGoogle Scholar
- Rodionov VI, Borisy GG: Microtubule treadmilling in vivo. Science. 1997, 275: 215-218. 10.1126/science.275.5297.215.View ArticleGoogle Scholar
- Jégou A, Niedermayer T, Orbán J, Didry D, Lipowsky R, Carlier MF, Romet-Lemonne G: Individual Actin Filaments in a Microfluidic Flow Reveal the Mechanism of ATP Hydrolysis and Give Insight Into the Properties of Profilin. PLoS Biology. 2011, 9: e1001161-10.1371/journal.pbio.1001161.View ArticleGoogle Scholar
- Turing AM: The Chemical Basis of Morphogenesis. Philosophical Transactions of the Royal Society B-Biological Sciences. 1952, 237: 37-72. 10.1098/rstb.1952.0012.View ArticleGoogle Scholar
- Carlsson AE: Actin Dynamics: From Nanoscale to Microscale. Annual Review Of Biophysics. 2010, 39: 91-110. 10.1146/annurev.biophys.093008.131207.View ArticleGoogle Scholar
- Weiner OD, Marganski WA, Wu LF, Altschuler SJ, Kirschner MW: An actin-based wave generator organizes cell motility. PLoS Biology. 2007, 5: e221-10.1371/journal.pbio.0050221.View ArticleGoogle Scholar
- Doubrovinski K, Kruse K: Cell motility resulting from spontaneous polymerization waves. Physical Review Letters. 2011, 107: 258103-View ArticleGoogle Scholar
- Doubrovinski K, Kruse K: Cytoskeletal waves in the absence of molecular motors. EPL. 2008, 83: 18003-10.1209/0295-5075/83/18003.View ArticleGoogle Scholar
- Whitelam S, Bretschneider T, Burroughs NJ: Transformation from Spots to Waves in a Model of Actin Pattern Formation. Physical Review Letters. 2009, 102: 198103-View ArticleGoogle Scholar
- Carlsson AE: Dendritic Actin Filament Nucleation Causes Traveling Waves and Patches. Physical Review Letters. 2010, 104: 228102-View ArticleGoogle Scholar
- Erlenkämper C, Kruse K: Uncorrelated changes of subunit stability can generate length-dependent disassembly of treadmilling filaments. Physical Biology. 2009, 6: 046016-10.1088/1478-3975/6/4/046016.View ArticleGoogle Scholar
- Bretschneider T, Anderson K, Ecke M, Müller-Taubenberger A, Schroth-Diez B, Ishikawa-Ankerhold HC, Gerisch G: The three-dimensional dynamics of actin waves, a model of cytoskeletal self-organization. Biophysical Journal. 2009, 96: 2888-2900. 10.1016/j.bpj.2008.12.3942.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.