This Article Abstract Full Text (PDF) All Versions of this Article: biophysj.106.091793v1 92/9/3092    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Hunyadi, V. Articles by Jánosi, I. M. Search for Related Content PubMed PubMed Citation Articles by Hunyadi, V. Articles by Jánosi, I. M.

Metastability of Microtubules Induced by Competing Internal Forces

Department of Physics of Complex Systems, Eötvös University, Pázmány Péter s. 1/A, H-1117 Budapest, Hungary

Correspondence: Address reprint requests to Imre M. Jánosi, E-mail: janosi{at}lecso.elte.hu.

	ABSTRACT

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Recent modeling efforts to estimate energies of tubulin-tubulin bonds shed light on a delicate balance between competing mechanical forces maintaining microtubule walls. Here we formulate two important refinements to the explanation of bond energetics. First, energy surface calculations in the elastic filament approximation reveal a finite stabilizing barrier assumed a simple Lennard-Jones-like potential for protein bonds. The presence of a guanosine triphosphate (GTP) cap represented by straight segments is necessary, as it is predicted for a long time. In the lack of such a cap, the protofilaments are either in an absolutely stable or absolutely unstable state. Second, our calculations show that this barrier appears only if the mechanical energy associated with the conformational change after GTP hydrolysis (curling energy) is larger than the strength of lateral bonds. The overall energy balance we propose supports continuous assembly of GTP dimers, a metastable state in the presence of a finite GTP cap and energetically driven disassembly of guanosine diphosphate protofilaments.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Tubulin -ß dimers in a microtubule (MT) are most frequently arranged in 13 protofilaments aligned parallel to the central axis (1). When MTs depolymerize, these protofilaments curve outward as a consequence of a conformational change, which is triggered by the hydrolysis of guanosine triphosphate (GTP) bound to ß-tubulin. Direct evidence of this long assumed conformational change has been recently presented by Wang and Nogales (2). They found that in guanosine diphosphate (GDP)-bound dimers, the contacts between the - and ß-tubulins within and between dimers are both affected, resulting in a curved protofilament that cannot form lateral contacts.

According to the "traditional" thermodynamic-kinetic view, GTP-bound tubulin subunits have a high affinity for MT ends, promoting persistent MT growth, whereas GDP-bound tubulin dimers have a low affinity and dissociate quickly. A kinetic lag between polymerization and hydrolysis could generate a "GTP cap" formed by almost straight GTP dimers, which is presumed to stabilize growing ends (3).

A crucial question in understanding dynamic instability (the random switch between continuous growth and rapid shortening) is how the chemical energy from GTP hydrolysis is exploited to power both growth and shrinkage of MTs. Recent modeling studies, especially by VanBuren et al. (4,5) and Molodtsov et al. (6,7), seem to reach a conclusion on quantitative aspects of bond energetics of the MT wall.

The stability of an MT is determined by the bonds between tubulin subunits forming the wall lattice. It is plausible to characterize the strength of bonds by energy parameters associated with the depth of (not precisely known) potentials of protein-protein interaction. For a particular subunit, the energetic balance (total energy), E_tot, contains at least three terms: net longitudinal bond energy, E_long; effective lateral bond energy, E_lat; and curling energy, E_c. (Note that any effect of possible external deformations is omitted in the following considerations.)

The individual terms in the energy balance for tubulin subunits have been estimated by various methods. First of all, stochastic modeling of the chemical kinetics (4) yields an estimate of per dimer for the net standard free-energy promoting spontaneous assembly of GTP tubulin (k_B is the Boltzmann constant, and T is the absolute temperature in the physiological range). This estimate is further supported by a different study on energetically unfavorable MT configurations, concluding that the stabilizing free-energy is around E_tot –10.5 k_BT per dimer for 13 protofilament GDP MTs (8). Note that the difference between and E_tot can be attributed to the contribution of internal curling triggered by the hydrolysis step (see below).

Second, depending on the assembly rate constant, longitudinal and lateral bond energies are estimated by the stochastic assembly model in VanBuren et al. (4) (Table 1). The free-energy shift associated with the conformational change after GTP hydrolysis is estimated at 2.1–2.5 k_BT/dimer by pure kinetic modeling (4) and 3.7–4.2 k_BT/dimer by considering flexural rigidity and measured shortening velocities (5). Similar energy values are difficult to extract from the molecular-mechanical model of Molodtsov et al. (6) and Molodtsov et al. (7) because the parameterization of interactions has no direct link to measured absolute values (e.g., longitudinal bonds are not represented by any potential; they are not extensible and do not break).

The helical dimer arrangement in the lateral direction (1) has a negligible effect on the overall MT stability and mechanical properties, as demonstrated in Molodtsov et al. (6). We will exploit this result in our calculations by considering a single protofilament in a laterally symmetric effective environment.

In this work we refine this picture by two main points. First, we explicitly show that a finite energy barrier can naturally arise from the superposition of a simple quadratic bending- and a Lennard-Jones-like bonding potential. A single GTP tubulin ring at the growing end is necessary and sufficient to maintain a metastable equilibrium, in agreement with many earlier predictions. Second, we illustrate that the driving force of rapid disassembly is easily available when the conformational change destabilizes the lattice, i.e., intrinsic bending forces can break lateral bonds.

Energy contour calculations
We performed computations by means of the elastic filament approximation described in Jánosi et al. (13). In this model, a protofilament is considered as a discretized one-dimensional string bound laterally to neighboring filaments. The tubular symmetry allows a simplified description of lateral bonds by an effective radial potential. Such a filament obeys also a longitudinal bending constraint represented by a nonzero intrinsic curvature for GDP segments and a quadratic potential associated with it. Longitudinal and lateral bonds can be characterized by Lennard-Jones potentials (Fig. 1):

(1)

where V₀ is the depth of the potential valley located at r₀. We also tested the Morse potential, which has a shape very similar to V_LJ with an adjustable width parameter (see Fig. 1):

(2)

View larger version (10K): [in this window] [in a new window]   FIGURE 1  Potential functions implemented in the energy map calculations; see Eqs. 1–3. Parameters are given in Figs. 2–4.

View larger version (13K): [in this window] [in a new window]   FIGURE 2  Illustration of the energy contour map calculation. (Left) Relaxed filament configuration with a unit cap defines the equilibrium position of the free terminal. (Right) The very end is forced to a given position, and the filament is allowed to relax into a minimal energy configuration. The final energy is used to draw contour lines. (The horizontal axes are lengthened; units are in nm.) Parameters: Morse potential, Etot = –10.5 kBT/dimer, Ec = +3.0 kBT/dimer, Elat = –2.0 kBT/dimer, = 1.12 nm–1.

View larger version (10K): [in this window] [in a new window]   FIGURE 4  Reconstruction of the energy barrier (along similar ridges indicated in Fig. 3 middle by a heavy line) for Morse potential maps with Elong = –11.5 kBT/dimer, Elat = –2.0 kBT/dimer, and Ec = +3.0 kBT/dimer, = 1.12 nm–1. Changing parameter is the cap size: no cap (crosses), unit cap (open squares), two-unit cap (solid circles).

View larger version (28K): [in this window] [in a new window]   FIGURE 3  Energy contours (in units of kBT) for a filament of 40 units as a function of the spatial separation from the spontaneous equilibrium terminal position (req, zeq) denoted by a heavy dot. (Top) Lennard-Jones bond potentials (1), with parameters from VanBuren et al. (4,5): Elong = –9.4 kBT/dimer, Elat = –3.2 kBT/dimer, and Ec = +2.1 kBT/dimer, without GTP cap. (Middle) Lennard-Jones bond potentials, Elong = –11.5 kBT/dimer, Elat = –2.0 kBT/dimer, and Ec = +3.0 kBT/dimer, with one unit of GTP cap. (The heavy line indicates the ridge where the resulting potential is reconstructed in Fig. 4.) (Bottom) Morse potentials, the same bond parameters as in the middle, = 1.85 nm–1, with one unit of GTP cap.

(3)

where r₀ defines the location of potential energy minimum again, and is a parameter characterizing the distance and relative height of the built-in energy barrier (see Fig. 1). Such potential entails an "automatic" bond breaking for large enough (r – r₀) spatial separations. Note that an additional term to Eq. 3 easily produces a negative minimum together with an asymptotic zero value (for details, see Molodtsov et al. (6,7)); however this gives only an additive shift in the final energy balance (6). The principal difference between Eq. 3 and the Lennard-Jones-like potentials in Eqs. 1 and 2 is that the former has a long-range repulsive regime, which is not really supported by computed (14–17) or measured (18) results.

We note that VanBuren et al. implemented simple quadratic potentials for each bonds in their mechanochemical model (5); bond breaking was realized "by hand" when external mechanical forces exceeded critical values.

We have computed energy contour maps for elastic filaments represented by a longitudinally joined sequence of straight rods with one clamped and one free terminal. For the sake of simplicity, only radial deflections were considered, as they manifest the lowest energy deformation modes. Longitudinal and lateral bonds were characterized by Lennard-Jones or Morse potentials, and bending constraints were implemented by a prescribed equilibrium angle between segments (22°) and a quadratic potential. The very end of the free terminal was forced to have a fixed spatial position, and then the filament was allowed to be relaxed by global conjugate gradient minimization (see, e.g., Jánosi et al. (19)). The resulting minimum energy was utilized to construct contour maps reflecting energetic stability (Figs. 2 and 3). Note that this energy is not normalized by the number of segments; it measures the difference between the relaxed global minimum value and the energy of a given forced configuration for a filament of arbitrary length.

The parameters of the potentials were tuned to have a minimum at the effective tube radius r₀ = 10.72 nm and at the equilibrium dimer length of 8.08 nm (see, e.g., Chrétien and Fuller (1)). The numerical value for total energy, E_tot = –10.5 k_BT, per segment was set for a GDP unit buried deep in the MT wall, following VanBuren et al. (4) and Hunyadi et al. (8).

The effects of a GTP cap were modeled by setting the equilibrium angle to 0° (perfectly straight limit) or 5° (see Wang and Nogales (2)) with the same bond energy parameters for a few top segments of a filament. Note that such a finite cap does not necessarily adopt a straight equilibrium configuration (see Fig. 2 a) because some curling deformation propagates upward from the GDP parts below (6,13).

	RESULTS AND DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Representative energy maps are shown in Fig. 3. We have scanned Lennard-Jones and Morse potentials with various (E_long, E_lat, and E_c) combinations by keeping the value of E_tot fixed and with straight cap segments of 0–4 units. The spatial resolution of the mapping was r = z = 0.1 nm in both directions. The curling energy, E_c, and the lateral bond strength, E_lat, were changed in the range 1–10 k_BT/dimer (E_long automatically arises from the balance). Combinations giving finite, positive resulting energy barriers are listed in Table 2 for Lennard-Jones potentials with unit and two-unit caps.

View this table: [in this window] [in a new window]   TABLE 2  The height Eb and location (r,z) of the energy barrier for a filament with unit and two-unit GTP caps (under the line)

View larger version (19K): [in this window] [in a new window]   FIGURE 5  Barrier energy Eb (in units of kBT) as a function of lateral bond energy, Elat (kBT/dimer) and curling energy, Ec (kBT/dimer), Lennard-Jones potential, Etot = –10.5 kBT/dimer. Data are from Table 2. Gray shading indicates the metastable region where the finite energy barrier appears (see Fig. 4).

1) In the case of pure GDP filaments, the necessary condition for a stable tube configuration is |E_lat| > E_c. This plausible requirement was also demonstrated in Molodtsov et al. (6). We found that such energy balance always yields to an absolutely stable equilibrium with Lennard-Jones or Morse potentials, as shown in Fig. 3, top. When the relationship is reversed, |E_lat| < E_c, the filament is absolutely unstable: it always curls up spontaneously.

2) Our calculations explicitly show that a metastable equilibrium is possible with |E_lat| < E_c in the presence of a GTP cap (at least one straight segment). In this case, an energy barrier appears as shown in Fig. 3, middle and bottom. The height and location of the barrier depends on the type and parameters of the potentials and the length of the cap segment (see Table 2 and Figs. 3 and 4). Largest barriers arise with long caps or narrow potential valleys at small differences between curling and lateral bond energies.

3) Metastability breaks down when the curling energy, E_c, overly exceeds |E_lat|. This limit is indicated by a steeply decreasing barrier height (illustrated in Fig. 5). The calculated critical lines are depicted in Fig. 6 both for Lennard-Jones and Morse potentials with various cap sizes.

View larger version (12K): [in this window] [in a new window]   FIGURE 6  Stability diagram for (a) Lennard-Jones and (b) Morse potential bond representations with different cap lengths (indicated by numbers at the border of metastable regions). Etot = –10.5 kBT/dimer, and = 1.12 nm–1 for the Morse potential.

An energy barrier implicitly arises with "realistic" bond potentials (without built-in or "handmade" bond breaking) in the presence of a single straight or moderately curved segment. The parameter sets we found support both continuous MT growth with a GTP cap (E_tot < 0) and energetically driven spontaneous disassembly of GDP protofilaments. The essential point is the possibility of |E_lat| < E_c in this case bending forces alone can break lateral bonds in the GDP state. Such bond breaking is possible only at the tube terminals; GDP units deeply inside the wall are fixed by longitudinal bonds (kinetic barrier). This picture suggests that catastrophes can occur also in the case of strong enough thermal fluctuations; the complete loss of the GTP cap is not a precondition.

We emphasize that the values of prescribed equilibrium angles are not important for the existence of metastable equilibrium, whenever their difference is large enough. When the GTP segments are assumed to not be perfectly straight, e.g., the intrinsic angle is 5° (2), the domain of metastability shrinks by 10% in the maps shown in Fig. 6.

Similarly, the parameter values for longitudinal bonds do not affect the general behavior. This is because the model here intends to capture the initial phase of disassembly characterized by the famous "ram's horn" configuration: protofilaments are curling off, but longitudinal bonds are not broken yet. The parameter values E_long < 0 do not permit fast disassembly of single protofilaments. Such delay between the lateral and longitudinal bond-breaking events is supported by the observation of many curled short GDP tubulin segments in a solution after MT catastrophes (21).

Our results are not in full agreement with some conclusions of VanBuren et al. (4,5). Only the parameter set (iii) listed in Table 1 obeys the criterion |E_lat| < E_c, which we found to be crucial for fast disassembly of an energetic driving force. With parameter sets (i), (ii), and (iv) in Table 1, breaking the top ring of lateral bonds (assumed to be the first step of a catastrophe) results in higher energy than the initial state of an intact tube. This is because the gain in curling energy is less than the loss of lateral bond contributions. Then the question naturally arises: what drives fast and continuous disassembly in the absence of permanent external forces? Thermal fluctuations obviously cannot propel depolymerization with a speed 10–20 times greater than MT growth.

Another conclusion was that sheet-like MT tip configurations were much more likely to undergo catastrophes than blunt ends (5), simply because unclosed sheets have many dimers with one-sided lateral bonds. The gain from global structural relaxation (see Jánosi et al. (19), Müller-Reichert et al. (20), and Chrétien et al. (21)) for sheet-like MT tips was not considered in VanBuren et al. (5); however this is important in the energy balance: weaker bond bending frustrations yield to higher energy barriers similarly to the case of larger GTP caps (see Fig. 4). A quantitative analysis of sheets is not possible in the framework of a single filament approximation, but the results of Molodtsov et al. (6), Jánosi et al. (19), and Chrétien et al. (21) strongly support that structural relaxation increases the stability at the MT terminal.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We thank Collegium Budapest Institute for Advanced Study for access to their computing facilities.

This work was supported by the Hungarian Science Foundation (OTKA) under grant No. TS044839. I.M.J. is grateful for a János Bolyai research scholarship of the Hungarian Academy of Sciences.

Submitted on June 20, 2006; accepted for publication January 4, 2007.

	REFERENCES

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
1. Chrétien, D., and S. D. Fuller. 2000. Microtubules switch occasionally into unfavorable configurations during elongation. J. Mol. Biol. 298:663–676.[CrossRef][Medline]

2. Wang, H.-W., and E. Nogales. 2005. Nucleotide-dependent bending flexibility of tubulin regulates microtubule assembly. Nature. 435:911–915.[CrossRef][Medline]

3. Mitchison, T., and M. Kirschner. 1984. Dynamic instability of microtubule growth. Nature. 312:237–242.[CrossRef][Medline]

4. VanBuren, V., D. J. Odde, and L. Cassimeris. 2002. Estimates of lateral and longitudinal bond energies within the microtubule lattice. Proc. Natl. Acad. Sci. USA. 99:6035–6040.[Abstract/Free Full Text]

5. VanBuren, V., L. Cassimeris, and D. J. Odde. 2005. Mechanochemical model of microtubule structure and self-assembly kinetics. Biophys. J. 89:2911–2926.[Abstract/Free Full Text]

6. Molodtsov, M. I., E. A. Ermakova, E. E. Shnol, E. L. Grishchuk, J. R. McIntosh, and F. I. Ataullakhanov. 2005. A molecular-mechanical model of the microtubule. Biophys. J. 88:3167–3179.[Abstract/Free Full Text]

7. Molodtsov, M. I., E. L. Grishchuk, A. K. Efremov, J. R. McIntosh, and F. I. Ataullakhanov. 2005. Force production by depolymerizing microtubules: a theoretical study. Proc. Natl. Acad. Sci. USA. 102:4353–4358.[Abstract/Free Full Text]

8. Hunyadi, V., D. Chrétien, and I. M. Jánosi. 2005. Mechanical stress induced mechanism of microtubule catastrophes. J. Mol. Biol. 348:927–938.[CrossRef][Medline]

9. Erickson, H. P., and E. T. O'Brien. 1992. Microtubule dynamic instability and GTP hydrolysis. Annu. Rev. Biophys. Biomol. Struct. 21:145–166.[Medline]

10. Hyman, A. A., S. Salser, D. N. Drechsel, N. Unwin, and T. J. Mitchison. 1992. Role of GTP hydrolysis in microtubule dynamics: information from a slowly hydrolyzable analogue, GMPCPP. Mol. Biol. Cell. 3:1155–1167.[Abstract]

11. Caplow, M., and J. Shanks. 1996. Evidence that a single monolayer tubulin-GTP cap is both necessary and sufficient to stabilize microtubules. Mol. Biol. Cell. 7:663–675.[Abstract]

12. Desai, A., and T. J. Mitchison. 1997. Microtubule polymerization dynamics. Annu. Rev. Cell Dev. Biol. 13:83–117.[CrossRef][Medline]

13. Jánosi, I. M., D. Chrétien, and H. Flyvbjerg. 2002. Structural microtubule cap: stability, catastrophe, rescue, and third state. Biophys. J. 83:1317–1330.[Abstract/Free Full Text]

14. Mitchell, J. B. O., R. A. Laskowski, A. Alex, and J. M. Thornton. 1999. BLEEP—potential of mean force describing protein-ligand interactions: I. Generating potential. J. Comp. Chem. 20:1165–1176.[CrossRef]

15. Jiang, L., Y. Gao, F. Mao, Z. Liu, and L. Lai. 2002. Potential of mean force for protein-protein interaction studies. Proteins. 46:190–196.[CrossRef][Medline]

16. Sept, D., N. A. Baker, and J. A. McCammon. 2003. The physical basis of microtubule structure and stability. Protein Sci. 12:2257–2261.[Abstract/Free Full Text]

17. Drabik, P., S. Gusarov, and A. Kovalenko. 2007. Microtubule stability studied by three-dimensional molecular theory of solvation. Biophys. J. 92:394–403.[Abstract/Free Full Text]

18. Moon, Y. U., R. A. Curtis, C. O. Anderson, H. W. Blanch, and J. M. Prausnitz. 2000. Protein-protein interactions in aqueous ammonium sulfate solutions. Lysozyme and bovine serum albumin (BSA). J. Solution Chem. 29:699–718.[CrossRef]

19. Jánosi, I. M., D. Chrétien, and H. Flyvbjerg. 1998. Modeling elastic properties of microtubule tips and walls. Eur. Biophys. J. 27:501–513.[CrossRef][Medline]

20. Müller-Reichert, Th., D. Chrétien, F. Severin, and A. A. Hyman. 1998. Structural changes at microtubule ends accompanying GTP hydrolysis: information from a slowly hydrolyzable analogue of GTP, guanylyl (,ß)methylenediphosphonate. Proc. Natl. Acad. Sci. USA. 95:3661–3666.[Abstract/Free Full Text]

21. Chrétien, D., I. M. Jánosi, J. C. Taveau, and H. Flyvbjerg. 1999. Microtubule's conformational cap. Cell Struct. Funct. 24:299–303.[CrossRef][Medline]

This Article Abstract Full Text (PDF) All Versions of this Article: biophysj.106.091793v1 92/9/3092    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Hunyadi, V. Articles by Jánosi, I. M. Search for Related Content PubMed PubMed Citation Articles by Hunyadi, V. Articles by Jánosi, I. M.