Skip to main content
Advertisement
  • Loading metrics

Microtubule instability driven by longitudinal and lateral strain propagation

Abstract

Tubulin dimers associate longitudinally and laterally to form metastable microtubules (MTs). MT disassembly is preceded by subtle structural changes in tubulin fueled by GTP hydrolysis. These changes render the MT lattice unstable, but it is unclear exactly how they affect lattice energetics and strain. We performed long-time atomistic simulations to interrogate the impacts of GTP hydrolysis on tubulin lattice conformation, lateral inter-dimer interactions, and (non-)local lateral coordination of dimer motions. The simulations suggest that most of the hydrolysis energy is stored in the lattice in the form of longitudinal strain. While not significantly affecting lateral bond stability, the stored elastic energy results in more strongly confined and correlated dynamics of GDP-tubulins, thereby entropically destabilizing the MT lattice.

Author summary

The dynamic nature of microtubules, long and hollow tubes formed by αβ-tubulin proteins, is crucial for their function is cells, and its precise characterization has been a long-standing problem for cell scientists. Microtubules are essential for cargo transport and provide mechanical forces in chromosome segregation when they disassemble. The disassembly proceeds via changes in the shapes of tubulins upon consumption of a chemical fuel called GTP that binds to every tubulin molecule. This leads to the accumulation of mechanical tension inside the microtubule and ultimately drives it beyond the stability threshold. However, it is still elusive how and where these shape changes contribute to the rapid release of the stored elastic energy. Here, we investigate the behavior of tubulin dimers in a microtubule-like environment using extensive atomistic simulations and show that tubulins locked in the microtubule operate as both ‘loadable springs’ and ‘conformational switches’, tightly controlled by their surrounding neighbours. We further show how these shape changes potentially control the overall stability of the microtubule, providing quantitative estimates of the system’s energetics.

Introduction

Microtubules (MTs) are one of the major components of the eukaryotic cytoskeleton and essential for intracellular transport, cell motility, and chromosome separation during mitosis. These are filamentous assemblies of αβ-tubulin dimers stacked head-to-tail in polar protofilaments (PFs) and folded into hollow tubes via lateral interactions [1, 2] (Fig 1A). Each dimer binds two GTP molecules of which only the one bound to β-tubulin is hydrolyzed in the MT lattice over time [3, 4]. This hydrolysis reaction is fundamental to MT dynamic instability [5], i.e. random switching between phases of growth and shrinkage (Fig 1A). Remarkably, both slow assembly and rapid disassembly of MTs—the latter termed catastrophe—are able to perform mechanical work because each tubulin dimer is a storage of chemical energy [68].

thumbnail
Fig 1. Tubulin life cycle and lattice compaction upon GTP hydrolysis.

(A), Cartoon representation of structural intermediates in MT assembly and disassembly. Individual dimers are composed of α-tubulins (gray circles) and β-tubulins (orange circles when GTP-bound or cyan circles when GDP-bound). Lattice cross-sections (bottom) indicate the location of the seam interface. (B), Local conformational changes proposed to accompany GTP hydrolysis are shown schematically (viewed from within the lumen). Each monomer is illustrated as two domains: intermediate or I and nucleotide-binding or N (C-terminal domains are not shown for simplicity). Rearrangements in α-tubulin around the nucleotide-binding pocket at the inter-dimer interface result in a ∼0.2-nm lattice compaction. The PFs are aligned with respect to monomer βi (marked with a circle). Other more subtle changes (e.g., PF twisting) or intermediate nucleotide states (e.g., GDP-Pi) are not shown for simplicity.

https://doi.org/10.1371/journal.pcbi.1008132.g001

The switch from a relaxed ‘curved’ conformation of tubulin favored in solution to a higher-energy ‘straight’ one is inherent to MT assembly [915]. It allows growing MTs to recruit and temporarily stabilize GTP-tubulin in the straight form, most likely due to the greater bending flexibility of GTP-PFs at intra- and inter-dimer interfaces [13, 1618]. It is therefore conceivable that collapsing MTs would follow a reverse pathway during disassembly; namely, they would release the conformational tension stored in GDP-tubulins that lateral bonds can no longer counteract. However, due to the system complexity and together with the inability of modern structural methods to directly visualize all sequential steps in the GTPase cycle in the straight MT body at high resolutions, it is still unknown exactly how and where the hydrolysis energy is converted to mechanical strain in the lattice.

Recent high- and low-resolution structural studies have revealed, in line with the early finding [19], that the use of a non-hydrolyzable GTP analog, GMPCPP, for MT assembly results in a more expanded MT lattice compared to a fully hydrolyzed GDP-lattice [2024], which is commonly interpreted as the lattice response to GTP hydrolysis (Fig 1B). Because by itself this global rearrangement does not fully indicate how and whether at all it is linked to GTP hydrolysis and strain accumulation at the single-dimer level, several competing models of MT cap maturation and MT disassembly have been proposed. According to the seam-centric or strain model [2022], the gradual build-up of longitudinal tension along the lattice upon GTP hydrolysis is the primary source of MT instability, where the lateral interfaces play only a passive role. In this model, lattice rupture is initiated at the seam because of the greater distance, and presumably weaker interactions, between PFs at this interface observed in unsymmetrized cryo-EM reconstructions. The role of the MT seam as the weakest interface is supported by recent computational evidence [25], but has been challenged experimentally [26]. In contrast, the holistic or bond model [23, 27] assumes that MT catastrophe can be explained by a sequential weakening of lateral inter-dimer contacts accompanied by a simultaneous strengthening of longitudinal contacts. Also here, recent atomistic simulations of full MT lattices suggest that lateral interactions in GDP-MTs might be weaker than those in GTP-MTs [25]. Finally, the most recent ‘no expansion’ model [24] provides an alternative view of the cap maturation process in which both pure GTP- and pure GDP-MTs have equally compacted lattices, while the higher-energy expanded lattice induced by GTP hydrolysis (mimicked by GMPCPP) corresponds to an intermediate, phosphate-releasing state. This model is partially supported by the observation that the extent of lattice compaction differs across different eukaryotic species [28].

The coexistence of the three models originates from the fact that the interplay between tubulin intrinsic strain and lateral binding inside the straight MT body is largely unclear. Indeed, the subtle changes in lattice compaction and dimer-dimer contacts would be best studied within straight PF assemblies in the presence or absence of lateral neighbors and conditioned on a fixed nucleotide state, which has not yet been achieved. This has prompted us to assess the mechanochemistry of both lattice compaction and lateral inter-dimer coupling using extensive molecular dynamics (MD) simulations of (i) isolated PFs, (ii) standard (homotypic) double-PF systems, as well as (iii) three-PF lattice patches. In all cases the PFs are locked in the straight conformation due to the use of periodic boundaries along the MT axis. Essentially and by construction, this setup allows to probe both mechanics and lateral cooperativity of individual dimers embedded in straight lattice regions distant from the dynamic tip. Because the three models of MT disassembly assume different properties of assembled tubulins and their lateral interactions, it is hence possible to test all models directly. By focusing on small, controllable MT-like subsystems, our simulations provide new insights into the lattice mechanics and energetics that drive MT disassembly.

Results

GTP hydrolysis in β-tubulin stiffens individual PFs

If MTs accumulate longitudinal elastic strain upon GTP hydrolysis, one would expect them to change the mechanical properties of individuals PFs also in the absence of lateral interactions. We therefore asked how the nucleotide state affects both equilibrium conformation and elasticity of isolated PFs and how much mechanical energy can be potentially stored in a single dimer upon GTP hydrolysis. The recent cryo-EM reconstructions of MTs in non-hydrolyzable GMPCPP- (mimicking GTP) and GDP-state [21, 22] enabled us to construct atomistic models of isolated PFs (Fig 2A) using correlation-driven molecular dynamics [31] and to assess their equilibrium and elastic properties by atomistic MD simulations (see Methods for details on system preparation, cryo-EM model refinement, simulation protocol).

thumbnail
Fig 2. Elastic properties of isolated ‘infinite’ PFs.

(A), Simulation setup for the single-PF system. α-tubulin (gray) and β-tubulin (cyan) are shown in surface representation. Potassium and chloride ions are shown as orange and cyan spheres, respectively. Water molecules are hidden for clarity. Periodic box with the axial dimension Lz is marked by a black rectangle. (B), Equilibrium probability distributions of the dimer rise in the PFs obtained from stress-free simulations of the system in (A). Shaded areas show statistical uncertanties of the distributions estimated with umbrella sampling. Dashed lines indicate the dimer rise values observed in the cryo-EM densities of GMPCPP- and GDP-MTs. (C), Stress-strain curves calculated for the system in (A) in both GTP- (orange) and GDP-state (cyan). Strain is computed relative to the equilibrium dimer length of GDP-PF, and negative (positive) stresses correspond to PF compression (extension). Only separate fits to the positive and negative stress ranges are shown. (D), Bending stiffness parameters of GTP- and GDP-MTs calculated using the elastic moduli in (B) (all stress values) and for varying PF numbers (orange and cyan dots, respectively). Experimental values (dashed lines with shaded areas) represent inverse-variance weighted means and standard deviations that combine multiple independent thermal fluctuation measurements summarized in [29] and recently updated in [30].

https://doi.org/10.1371/journal.pcbi.1008132.g002

To assess the equilibrium properties of isolated PFs at room temperature, we first performed multiple simulations of GTP- and GDP-PFs totaling ∼23 μs and ∼13 μs, respectively, following a previously published protocol [32]. We monitored the dynamics of both tubulin dimer shape and axial periodic box size Lz, i.e. the lattice spacing (Fig 2A). There are three possibilities how the conformation of the PF could contribute to an increase or decrease of Lz: (a) changes at the intra-dimer interface between α- and β-subunits belonging to the same dimer (referred to as ‘dimer spacing’), (b) changes at the inter-dimer interface between α- and β-subunits belonging to neighboring dimers along the PF axis (referred to as ‘PF spacing’), and (c) changes in the shapes of α- and β-subunits due to elastic deformations. We then employed the Functional Mode Analysis [33, 34] to train a regression model on the dynamics of Lz and to derive a reaction coordinate that best describes the compaction/expansion dynamics of the PF in equilibrium (see Methods). This reaction coordinate, that we termed ‘dimer rise’ for consistency with cryo-EM experiments, was used in all subsequent free energy calculations.

Fig 2B shows the equilibrium probability distributions of the dimer rise as a function of the nucleotide state computed with additional umbrella sampling simulations (∼114 μs of cumulative simulation time; see Methods). Both GTP- and GDP-PFs slightly increased the lattice period during the simulations relative to their initial cryo-EM conformations. This slight elongation might be caused by thermal expansion of the 70–80 K cryo-EM structures after re-equilibration at room temperature. This possibility is supported by the observation that the relative elongation is largely independent of system size and sampling, as will be seen further. However, GTP-PFs maintained a significantly longer dimer rise compared with GDP-PFs (+0.25 ± 0.07 nm), consistent with the experimentally observed difference of ∼0.2 nm. In the following, we will refer to these two states as expanded and compacted. Hence, it is likely that the difference between the global states of MT structure seen by cryo-EM reflects a local response of tubulin to GTP hydrolysis; otherwise, it would vanish in the absence of lateral contacts.

Further, GTP-PFs sampled a wider range of dimer rise values as indicated by the distribution widths in Fig 2B, which suggests that GTP-PFs are mechanically more flexible. We previously showed that, when placed in solution, GTP-tubulin exhibits higher bending flexibility than GDP-tubulin [17]. It was therefore surprising that, when tubulin was locked in the straight MT-like conformation, also the longitudinal elasticity of the dimer was affected by the nucleotide state.

To quantify the mechanical elasticity of the system in Fig 2A, we performed a set of steady-state compression/extension simulations at constant values of the axial component Pzz of the pressure tensor (along the PF). Fig 2C shows the obtained strain-stress curves, where the strain was computed relative to the equilibrium conformation of GDP-PFs. Irrespective of the nucleotide state, the stress-strain data clearly falls into two elastic regimes: a rather soft response for positive stresses (extension) and a much stiffer response for negative stresses (compression). This previously observed behavior of GDP-tubulin [32], which we here confirmed using higher quality structures of both nucleotide states and wider strain ranges, emerges likely because different parts of the heterodimer are involved in the mechanical response upon compression or extension. Whereas extension mainly stretches the inter-dimer and, to a lesser extent, intra-dimer interfaces, compression forces individual monomers to change their shapes, causing much more resistance. We therefore analyzed the positive and negative strain ranges separately.

A linear fit to the negative stress data of the GTP- and GDP-PF simulations yielded elastic moduli of 0.89 ± 0.07 GPa and 1.77 ± 0.13 GPa, respectively. Fitting to the positive stress data yielded systematically smaller moduli of 0.37 ± 0.07 GPa (GTP-PF) and 0.53 ± 0.03 GPa (GDP-PF). A fit to the entire stress range analyzed in our simulations resulted in values of 0.82 ± 0.06 GPa (GTP-PF) and 1.55 ± 0.15 GPa (GDP-PF) that agreed better with the moduli obtained by fitting the negative stress data only. Whereas the single-PF system tolerated high compression stresses up to Pzz = +200 bar without undergoing plastic deformations and irrespective of the nucleotide state, this was not the case for extension stresses. GTP-PFs withstood stretching up to Pzz = −65 bar without rupturing at the inter-dimer interface in the course of our simulations (∼1 μs each). In contrast, GDP-PFs ruptured already at stress values below Pzz = −40 bar, implying that a lower force is likely sufficient to break the longitudinal bond. Although more sampling would be required to investigate PF rupture pathways, our stress-strain data (Fig 2C) together with the equilibrium free energy calculations (Fig 2B) support the interpretation that straight GDP-PFs are stiffer than GTP-PFs while they might possess more fragile inter-dimer longitudinal bonds.

The factor of two difference in the elastic moduli of GTP- and GDP-PFs is remarkable given the high structural similarity of the two conformational states. Careful review of existing measurements of MT bending mechanics reveals that, despite variations in the experimental protocols and theoretical models used to analyze such data, MTs are intrinsically softer when polymerized in the presence of GMPCPP and/or Taxol [29, 30]. This is reflected, e.g., in a significantly distinct bending stiffness E × I, where E is the elastic modulus and I is the second moment of the cross-sectional area of the MT. We therefore asked if the nucleotide-dependent elasticity of PFs (Fig 2B and 2C) might explain the experimentally observed differences in coarse-grained elastic properties of MTs.

Bending stiffness of MTs is typically obtained by monitoring and quantifying their equilibrium fluctuations or by directly applying a force to bend MTs and then measure their resistance (e.g., by using optical tweezers) [29]. It is known that in thermal fluctuation experiments, MTs behave on average stiffer than in force-probing experiments [29, 30]. It has been proposed that this discrepancy can be reconciled by taking into account that large deformations caused by external forces acting on MTs could surpass the elastic limit, which would lead to non-elastic deformations of tubulin dimers and/or breakage of inter-dimer contacts [30]. Because such events are, by construction, unlikely to happen in our simulation setup, we compared our results only with equilibrium fluctuation experiments, where tubulin dimers are mostly subject to small-strain elastic deformations.

Fig 2D compares bending stiffnesses of hypothetical MTs with varying PF numbers calculated using our data in Fig 2C (all stress values) and consensus values calculated as precision-weighted averages from a pool of independent experimental measurements (reviewed in [29] and recently updated in [30]; see Methods). The comparison revealed that a good agreement between the experimental data and our calculations can only be achieved if the PF number is approximately 14 for both GTP- and GDP-MTs. It is known that MT mechanics is highly sensitive to changes in the PF number (see Discussion in [35]). Most MTs polymerized in vitro without co-factors and MT-binding drugs possess 14 PFs [36], with ratios of 13-PF to 14-PF MTs reaching approximately 1:9 for GDP-MTs and 1:3 for GMPCPP-MTs [21, 22]. Assuming that tubulin axial elasticity does not depend on the PF number, the two-fold higher bending stiffness of GDP-MTs can be accounted for almost entirely by a two-fold higher elastic modulus of GDP-PFs, at least for small-strain deformations.

Finally, the good agreement of our elasticity calculations with experimental knowledge allowed us to estimate that a free energy of ΔGel ≈ 11.6 kBT, where kB is the Boltzmann constant, would be stored in a GTP-PF per dimer when mechanically compressed to the state of a GDP-PF (see S1 Material). Remarkably, this energy is very close to both the energy harvested by MTs upon GTP hydrolysis [19, 37] and the maximal excess energy that can be stored in an MT lattice to maintain one of the most favorable configurations (∼11 kBT per dimer for MTs with 13 or 14 PFs [3840]). Together with the consistency of our calculated elastic moduli with the observed softening of GMPCPP-MTs or Taxol-stabilized MTs vs. GDP-MTs, this strongly suggests that almost the entire energy available from GTP hydrolysis is stored in the MT lattice in the form of longitudinal elastic strain. We note that, during model preparation, the GMPCPP molecules were manually converted to GTP, and the starting structures were allowed to adapt to this change in the subsequent production simulations (see Methods). However, it cannot be ruled out that some effects of GMPCPP on the dimer conformation may still remain. Nevertheless, in the absence of alternative high-resolution models of the putative GTP state, we consider our tubulin model a sufficiently good approximation of the true GTP-tubulin structure in MTs.

Lateral coupling and GTP hydrolysis reduce conformational freedom of tubulin in PFs

One of the key unanswered questions is how the MT lattice would accommodate laterally coupled dimers in conflicting conformational states (expanded vs. compacted), a situation that is very likely to arise downstream from the growing MT tip. It was previously speculated that such a structural conflict would either weaken the lateral interactions between incompatible dimers or increase the rate of GTPase activity [41, 42]. In the latter case, the hydrolysis-triggered compaction of an expanded dimer located next to a compacted dimer would be more favorable. However, testing these hypotheses experimentally is currently challenging.

To get insight into how the presence and conformation of a lateral neighbor affects the compaction-expansion dynamics of tubulin in PFs, we constructed atomistic models of double-PF systems in both nucleotide states (Fig 3A; see Methods for model refinement and simulation protocol). We then computed free energy surfaces of the double-PF systems as a function of dimer rise and nucleotide state using the umbrella sampling approach with ∼80 μs of cumulative simulation time (Fig 3B and 3C; see Figure A in S1 Material for statistical uncertainties). Like the isolated PFs (Fig 2B), the double-PF systems adopted, on average, slightly more expanded conformations relative to their starting cryo-EM structures, most likely, due to thermal expansion. Also, the constant shift between the two distributions by 0.19 ± 0.05 nm was preserved, which was close to the experimentally observed difference of ∼0.2 nm and, within statistical error, consistent with the difference of 0.25 ± 0.07 nm calculated for the isolated PFs (Fig 2B).

thumbnail
Fig 3. Lateral coupling and nucleotide state affect PF dynamics.

(A), Simulation setup for the double-PF system mimicking a standard (homotypic) lateral interface. Color coding as in Fig 2A. Water molecules are hidden for clarity. Periodic box is marked by a black rectangle. Individual PFs are labeled as (1) and (2). (B) and (C), Free energy surfaces of the system in (A) as a function of dimer rise and nucleotide state obtained by umbrella sampling. The surfaces are color-coded by free energy values with an increment of 1 kBT (dark red to gray). Black solid lines additionally show isoenergetic contours. Orange and cyan circles indicate the dimer rise values observed in the cryo-EM densities of GMPCPP- and GDP-MTs, respectively. Cartooned dimers in (B) schematically show the extreme conformations of the double-PF system in which both are similarly expanded or compacted (along the diagonal) or in conflicting conformations (along the anti-diagonal). The relative shift of 0.19 nm between the minima of the free energy sufraces in (B) and (C) is additionaly indicated.

https://doi.org/10.1371/journal.pcbi.1008132.g003

As described above, one would expect each PF in the double-PF system to behave differently depending on both conformational state of the neighbor and own nucleotide state. In particular, their motion should be statistically correlated due to lateral coupling. In addition, the substantial difference in mechanical flexibility of isolated GTP- and GDP-PFs (Fig 2B and 2C) should be reflected in the dynamics of coupled PFs as well. To test these expectations, we introduced two metrics to quantify the changes in the double-PF free energy profiles upon nucleotide exchange. First, we used the normalized mutual information (NMI) which is a mere statistical measure of both linear and nonlinear correlation between two stochastic variables (see Methods for the rigorous definition). For the particular system in Fig 3A, NMI would be zero if the PFs moved fully independently or unity if their dimer rise fluctuations were fully synchronized. Second, we used the confinement entropy Hconf that quantifies the conformational space ‘volume’ available to both PFs, irrespective of how much the PF motions are correlated (see Methods for the rigorous definition). For the particular system in Fig 3A, Hconf would be close to zero if the dimer rise fluctuations were strongly localized around fixed values or maximal if all dimer rise values were equally likely. Both stronger inter-PF correlation (higher NMI) and higher PF stiffness (lower Hconf) would naturally result into the PFs having a more restrictive influence on each other, moving the double-PF system up in free energy due to the associated loss of conformational entropy. Vice versa, the joint conformational space increases when the PFs become more flexible and their fluctuations become less correlated, hence moving the double-PF system down in free energy.

As visible from the free energy surfaces having elliptic shapes extended along the diagonal (Fig 3B and 3C) and supported by the calculated NMI values, the conformations of the double-PF system in which the PFs were similarly expanded or compacted were lower in free energy than those in which the PFs adopted conflicting conformations. Thus, the PFs have a mutually restrictive influence on each other, penalizing configurations in which the PF conformations are too different. As a result, the double-PF system exhibits more correlated dynamics than would be the case if the PFs were isolated. Furthermore, the correlation and confinement effect was stronger for the system in GDP-state (ΔNMI = +0.06 and ΔHconf = −1.0 bits compared to GTP-state). Together with the PF stiffening upon GTP hydrolysis (Fig 2B and 2C), this suggests that GTP hydrolysis further reduces the conformational space available to tubulin dimers in the double-PF system, making it thermodynamically less favorable than that in GTP-state.

Nearest-neighbor interactions between PFs modulate GTPase response of tubulin

Our observation that the double-PF system favors conformations in which the PFs are similarly expanded/compacted suggests that the system is less stable when there is a conformational mismatch between the PFs, likely because the lateral bond would be under excessive shear tension. To quantify the extent of lateral bond destabilization by the conformational mismatch between the PFs, we considered a thermodynamic cycle shown in Fig 4A, following a previous scheme [43]. We assume that the equilibrium conformation of the double-PF system can be changed into the one with a conformational mismatch between the PFs (vertical transitions in Fig 4A). The free energy cost associated with this transformations (ΔGeq→mis) was calculated using our previous umbrella sampling results for both single- and double-PF system (Figs 2B and 3B and 3C; see S1 Material). Because the sum over all transition paths in the cycle must vanish, the difference between these values, , equals the bond stability of the mismatched double-PF system relative to the equilibrium case, (horizontal transitions in Fig 4A). Hence, a positive ΔΔGassoc equally implies that (a) the PF association is less favorable when the PFs are in conflicting compaction states, or that (b) the GTPase response of an expanded dimer is stimulated if its nearest neighbor is in the compacted state.

thumbnail
Fig 4. Relative thermodynamic stability of the lateral bond in the double-PF system.

(A), Thermodynamic cycle demonstrating the idea behind estimating the effect of unequal PF conformations on the association free energy between the PFs. While simulating the horizontal transitions (PF association) is computationally more expensive, the free energy changes linked to the vertical transitions (PF compaction) have already been obtained (Figs 2 and 3). (B) and (C), Distributions of the relative stability of the double-PF systems with respect to their equilibrium conformations marked with orange and cyan circles for GTP- and GDP-state, respectively, as a function of dimer rise and nucleotide state. White circles denote conformations with the strongest observed dimer rise mismatch. Free energy color coding is adjusted such that red (blue) areas correspond to conformations of the double-PF system in which the lateral bond is destabilized (stabilized) relative to equilibrium. White areas correspond to no change in the lateral bond stability.

https://doi.org/10.1371/journal.pcbi.1008132.g004

Fig 4B and 4C shows ΔΔGassoc relative to the lowest-energy system configuration (free energy minima in Fig 3B and 3C, respectively) as a function of dimer rise and nucleotide state. The calculations suggest that a conformational mismatch between the PFs would have a statistically significant effect on the thermodynamic stability of the double-PF system, corresponding to a change of ΔΔGassoc = +4.0 ± 1.6 kBT (equilibrium constant fold-change by ∼55). In contrast, simultaneous compaction/expansion of the two PFs has no statistically significant effect on the stability of the double-PF system with a relative change of ΔΔGassoc = −1.0 ± 1.5 kBT (equilibrium constant fold-change by ∼0.37), implying that the lateral bond is stabilized once the conformational mismatch is resolved. Our results, therefore, provide quantitative evidence for the previous ideas that a structural conflict at the lateral interface due to unequal nucleotide states would either weaken it or locally increase the rate of GTPase activity [41, 42], i.e. locally facilitate the compaction transition. However, not only do we propose that both ideas would be equivalent, but we also estimate the magnitude of lateral bond destabilization and predict that it would be a transient and reversible effect. Our estimate for the bond destabilization energy in the absence of any lateral mismatch, ΔΔGassoc = −1.0 ± 1.5 kBT, also agrees well with a recent computational study by the Odde lab [44], where only a weak nucleotide dependence of the lateral bond stability was found, though using a finite PF setup.

We note that the association free energy differences in Fig 4A (representing lateral transitions) refer to a situation in which one long and straight PF fully associates with another. The considered scheme differs from how PFs most likely associate/dissociate at the dynamic MT tip, namely, dimer by dimer while bending away from the MT lumen. Hence, and describe per-dimer contributions to the lateral thermodynamic stability of MT lattices in regions distant from the dynamic MT tip.

Nearest-neighbor interactions between PFs cause long-range correlations in the lattice

The finding that lateral coupling leads to more confined and correlated dynamics of tubulin in the double-PF system is explained by the nearest-neighbor interaction that prevents dimers in the adjacent PFs from adopting conflicting conformations by energetically penalizing local mismatches. It is therefore clear that also the motions of dimers situated in distant PFs should be correlated as a consequence of the elementary short-range interactions shown in Fig 3B and 3C. However, it is unclear to what extent the nucleotide state would affect such long-range correlations. To quantify their magnitude and the dependence on the bound nucleotide in similar minimalist but computationally feasible settings, we constructed a larger PF system comprising 3 × 1 dimers per periodic box dimension Lz (Fig 5A and 5B), which allowed us to quantify the statistical correlation between a pair of non-adjacent PFs. From the equilibrium dynamics of this three-PF system, similarly as above, we estimated its free energy landscape as a function of dimer rise and nucleotide state and subsequently disentangled nearest-neighbor interactions and long-range correlations.

thumbnail
Fig 5. Lateral coupling induces long-range correlations between distant PFs.

(A) and (B), Side and top views of the simulation setup for the three-PF system mimicking a larger segment of the MT lattice. Color coding as in Figs 2A and 3A. Water molecules are hidden for clarity. Individual PFs are labeled as (1), (2) and (3). (C) and (D), Free energy energy landscapes of the system in (A) as a function of dimer rise and nucleotide state. The 3D landscapes were pairwise projected onto planes corresponding to 2D free energy landscapes of adjacent (α(1)β(1)α(2)β(2) and α(2)β(2)α(3)β(3), left and right, respectively) and non-adjacent PFs (α(1)β(1)α(3)β(3), center). Orange and cyan circles indicate the dimer rise values observed in the cryo-EM densities of GMPCPP- and GDP-MTs, respectively. Note the shift between the GTP and GDP distributions by ∼0.2 nm along both reaction coordinates, consistent with the other simulations in Figs 2 and 3.

https://doi.org/10.1371/journal.pcbi.1008132.g005

As it was unfeasible to perform sufficiently accurate free energy calculations for such a large system, we instead resorted to a Bayesian inference approach that integrates prior knowledge about the energetics of the smaller subsystems (Figs 2A and 3A) to infer the joint free energy distribution of the three-PF system from unbiased MD simulations (see S1 Material). To this end, six independent, 600-ns long equilibrium simulations of the three-PF system in each nucleotide state were performed, yielding a total of ∼7.2 μs of sampling time. The inferred three-dimensional (3D) joint free energy distributions were then pairwise projected onto planes corresponding to two-dimensional (2D) free energy landscapes of adjacent and non-adjacent double-PF subsystems (Fig 5C and 5D). Consistent with the single-PF and double-PF systems analyzed above (Figs 2 and 3), the conformation of the three-PF system in our simulations was more expanded than the underlying cryo-EM structures, while the nucleotide-dependent difference in lattice compaction, again, was preserved.

The NMI values for the non-adjacent free energy landscapes were calculated (NMI13) and compared with those for the adjacent landscapes in the same system (NMI12 and NMI23). If the non-adjacent PFs did not interfere, NMI13 would be negligible relative to both NMI12 and NMI23, yielding almost circular free energy landscapes in Fig 5C and 5D (center). However, we found that NMI13 is only by a factor ∼0.5 and ∼0.85 smaller than the values for the directly interacting PFs in GTP- and GDP-state, respectively. This suggests that the correlations between non-adjacent PFs induced by the nearest-neighbor PF interactions are enhanced upon GTP hydrolysis.

In fact, several recent findings provide intriguing evidence that weaker intra-lattice correlations might stabilize the MT. First, some MT-stabilizing drugs such as Taxol have been recently shown to increase the lattice heterogeneity of GDP-MTs as compared to drug-free GDP-MTs [45], which resonates with the ability of Taxol to restore the bending flexibility of GDP-MTs [46, 47]. Second, a very similar effect on MT stability and mechanical resilience has been reported for acetylated vs. wild-type MTs [48, 49], likely due to a small but additive allosteric effect of α-tubulin acetylation at residue K40 [50]. In light of our drug- and acetylation-free simulation results, we propose that GTP hydrolysis reduces tubulin axial flexibility and enhances short- and long-range correlations between PFs, thereby leading to a loss of conformational entropy by the MT lattice.

Discussion

Tubulin dimers locked in the MT lattice operate as ‘loadable springs’ and ‘conformational switches’ whose load and conformation critically depend on both nucleotide state and lattice surrounding. The result is a metastable behavior of MTs because they have to reconcile the favorable dimer-to-lattice binding and the internal strain build-up that is fueled by GTP hydrolysis and has a destabilizing effect. Not surprisingly, it is hard to reach a consensus on precisely where this chemical energy is converted to mechanical strain and spread over the lattice because the system complexity allows various interpretations.

To clarify this issue, we have aimed at a quantitative understanding of the interplay between tubulin intrinsic strain and lateral binding inside straight MT-like compartments. Our results support the following conclusions: (i) there is a two-fold increase in longitudinal lattice tension upon nucleotide hydrolysis, which we attribute to the increased stiffness of GDP-PFs; (ii) lateral coupling between PFs reduces the conformational flexibility of tubulin by entropically penalizing PF conformations that are too different; (iii) restrictive interactions between neighboring PFs induce long-range correlated motions of non-adjacent PFs; (iv) both short- and long-range cooperativity of PF motions is stronger for GDP-PFs suggesting a loss of conformational entropy by MT lattices upon GTP hydrolysis.

Our computational findings provide quantitative insights into tubulin mechanochemistry and, therefore, the structural and energetic basis of MT dynamic instability. The results presented here enable us to test the three major models of MT cap maturation and MT catastrophe: (i) the strain model proposed by Nogales, Alushin, Zhang and colleagues [2022], (ii) the bond model proposed by the Moores lab [23, 27], and (iii) the most recent ‘no expansion’ model by Estévez-Gallego et al. [24]. Our results do not support the bond model because we do not find a statistically significant effect of the nucleotide state on the lateral bond stability, which is key to the bond model. This has also been confirmed by an earlier computational study from the Odde lab using a finite PF setup [44]. Rather, our results are consistent with a weakening of longitudinal bonds upon GTP hydrolysis, based on the stronger response of single GDP-PFs to mechanical stretching (Fig 2C), which is indicative of more fragile longitudinal bonds between GDP dimers. Further, our results do not support the most recent model by Estévez-Gallego et al. [24] postulating that both pre- and post-hydrolysis MT segments are equally compact, while the lattice undergoes a transient, hydrolysis-energy-consuming expansion (approximated by the GMPCPP lattice) to release the γ-phosphate. Our combined evidence (Figs 25) rather suggests: (a) whatever the GMPCPP state corresponds to in the GTPase cycle of tubulin, it most likely precedes the GDP state; and (b) a pure GDP lattice is higher in free energy than a pure GMPCPP one because the transition to the GDP lattice can only be achieved by investing a per-dimer energy on the order of 11 kBT. Thus, the conformational cycle proposed by Estévez-Gallego et al. would involve around twice the energy available from GTP hydrolysis to first expand and then compact the lattice, which is energetically implausible according to our estimates.

Overall, our results are currently most consistent with the strain model by Nogales, Alushin, Zhang et al [2022]. However, until the status of the GMPCPP lattice is entirely clear, the strain model remains incomplete. This model does not describe the behavior of mixed nucleotide MT lattices, which we have now predicted using long-time MD simulations and free energy calculations. At present, the strain model cannot preclude the possibility of tubulin adopting unknown pre-hydrolysis conformations prior to that mimicked by GMPCPP. Although our simulations do not show large rearrangements of the GMPCPP-tubulin structure upon replacement of GMPCPP with GTP, we still have to assume that the GMPCPP-MT lattice is sufficiently similar to the unknown pre-hydrolysis GTP-MT lattice. Whatever the precise conformational cycle, our results agree best with the view of GMPCPP-tubulin being one of the cap-stabilizing, expanded, and flexible conformations that are unlikely to be preceded by stiffer and compacted ones. Perhaps a more promising approach toward ultimately resolving this issue in future structure determination efforts would be to use knowledge-based point mutations that selectively uncouple the tubulin conformational and GTPase cycles [51, 52].

Taken together, a new picture emerges in which the MT lattice stability is not exclusively determined by the nucleotide-dependent dynamics of individual dimers, but more generally, by their non-additive collective behavior. In this work, we provide a thermodynamic explanation for the intrinsic destabilization (distant from the dynamic tip) which precedes MT breakdown and which relies on the idea that MT lattices gradually accumulate mechanical strain and lose conformational entropy as GTP hydrolysis proceeds. In other words, the MT becomes thermodynamically less and less stable already during the growing phase, which predisposes it to explosive strain release. Exascale atomistic simulations (≫106 atoms) and coarse-grained kinetic models can now be used to extrapolate how the results of our study will apply to the time evolution of the MT plus-end tip at much larger spatiotemporal scales.

By construction, our lattice simulations focus on straight ‘infinite’ PF systems with only one dimer layer per periodic box length. As a result, each dimer interacts with itself along the MT axis. Such setups are established and well-tested in the MD field, and the resulting artifacts are well-characterised [32, 5355]. For the periodic tubulin systems at hand, two types of artifacts warrant attention. First, fluctuations with wavelengths larger than the box size are suppressed and, therefore, fluctuations of the dimer rise may be smaller than in a simulation with a much larger box size or in reality. Second, ‘diagonal’ correlations are not present in the double- and three-PF systems, i.e. a dimer is unable to influence the conformations of other dimers in the neighboring PFs located in the layer above or below that dimer. Another possible issue is the fact that our simulations did not include a closed segment of the MT body with 14 PFs so that possible ‘edge’ effects cannot be fully discounted. Contrary to the periodic boundaries, this simplification of the MT geometry might lead to more relaxed fluctuations of the dimer rise for those PFs having only a single neighbor. It is conceivable that, to some extent, the periodic boundaries mitigate the absence a closed MT lattice due to error cancellation.

It is therefore justified to ask what the simulated PF models actually represent in a real MT system and how strongly the chosen simulation setup affects the conclusions of our study. Two observations are important here. First, our PF elasticity calculations provide estimates that are largely consistent with previous experimental knowledge, indicating that longitudinal PF mechanics does not significantly depend on the choice of simulation protocol. Second, we primarily focus on the effect of the nucleotide state and always compare the results of GTP and GDP simulations. It is likely that, by considering relative changes, some of these artifacts cancel out and their effect on the conclusions is smaller than it would be for absolute values. Therefore, the periodic and finite size effects described above are unlikely to significantly affect the conclusions. Overall, we assume that our PF models sufficiently accurately describe the dynamics of straight MT lattice regions distant from the dynamic MT tip, with the important simplification that intra-lattice correlations are restricted to the lateral dimension.

Methods

Force-field parameters and protonation states

The CHARMM22* force field [56] and the CHARMM-modified TIP3P water model [57] were used in all simulations. GTP/GDP parameters were adapted from those for ATP/ADP implemented in the CHARMM22/CMAP force field [57, 58]. Titration curves of histidines were calculated using the GMCT package [59] and assigned as described previously [17].

Simulation system preparation and cryo-EM refinement

Initial models for the tubulin dimers were obtained from PDB IDs 3JAT (GMPCPP) and 3JAS (GDP) [21] by extracting the central dimer from the 3 × 2 lattice patches (chains A and H in the original PDBs). GMPCPP was converted into GTP by replacing the carbon atom between α- and β-phosphate with an oxygen atom. The missing loop in the α-subunit (residues 38-46) was modelled in for structure consistency using MODELLER version 9.17 [60] but excluded from further refinement. Unlike in our previous study [17], we did not include the missing C-termini (α:437–451 and β:426–445) in our simulations to reduce the system size and reach the best possible sampling. Unless differently specified, all structure and map manipulations were performed using UCSF Chimera [61] or VMD [62].

In all refinement simulations, the following data sets were used: EMD-6352 and EMD-6353 for symmetrized cryo-EM reconstructions of 14-PF GMPCPP- and GDP-MTs decorated with kinesin [21]. To create ‘infinite’ single-, double-, and three-PF systems, where the actual simulated part comprises exactly one layer of dimers and is coupled to copies of itself through axial periodic boundaries, we first constructed finite PF systems comprising two layers of dimers in the axial direction. To this end, subsections of the cryo-EM maps with the desired PF topology were extracted using an orthorhombic box, and the single dimer models were rigid-body fitted into the PF maps. The constructed PF systems were solvated in a triclinic water box of size 8.0 × 8.0 × 22.0 nm3 (single-PF), 12.7 × 12.7 × 22.0 nm3 (double-PF), or 19.0 × 19.0 × 22.0 nm3 (three-PF). The systems were then neutralized with 150mM KCl.

Refinement was done with correlation-driven molecular dynamics implemented as a custom module in the GROMACS 5.0.7 package [63], following our previously published protocols [31]. Briefly, we used the cold-fitting protocol with the longest refinement time (i.e. T = 100 K and total run time of 50 ns) followed by 15 ns of simulated annealing. The starting values for the biasing strength and the simulated map resolution were set to 1 × 105 kJ mol−1 and 0.6 nm and linearly ramped to 5 × 105 kJ mol−1 and 0.2 nm, respectively. The quality of the resulting models and the goodness of fit were ensured by calculating common stereochemical and correlation metrics (Table A in S1 Material).

MD simulations

The finite PF models were converted into ‘infinite’ PF models by removing the extra tubulin monomers and nucleotides. Water and ion atoms were then trimmed to conform to the experimental value of the axial periodic dimension Lz, namely, 8.31 nm for GMPCPP-MTs and 8.15 nm for GDP-MTs [21]. The number of ions in the trimmed water shell was fixed such as to keep the systems neutral and to maintain the ionic strength of 150mM KCl. All subsequent MD simulations were carried out with GROMACS 5.0.7 [63]. Lennard-Jones and short-range electrostatic interactions were calculated with a 0.95-nm cutoff, while long-range electrostatic interactions were treated using particle-mesh Ewald summation [64] with a 0.12-nm grid spacing. The bond lengths were constrained using the LINCS algorithm [65] (hydrogen bonds during equilibration and all bonds in the production runs). Velocity rescaling [66] with a heat bath coupling constant of 0.5 ps was used to control the temperature for solute and solvent separately. Applying virtual site constraints [67] allowed us to increase the integration step size to 4 fs in the production runs. Center-of-mass correction was applied to solute and solvent separately every 100 steps.

With the above parameters fixed, the equilibration protocol consisted of the following steps: (i) energy minimization using steepest descent; (ii) short NVT equilibration for 1ns at T = 100 K with position restraints on heavy atoms and using a 1-fs integration time step; (iii) gradually heating up the system to 300 K within 10 ns in the NPT ensemble (Berendsen barostat [68] with a 5-ps coupling constant) using a 2-fs integration time step; (iv) equilibration in the NPT ensemble for 30 ns using isotropic Parrinello-Rahman barostat [69] with a 5-ps coupling constant and using a 2-fs integration time step; (v) equilibration in the NPT ensemble for 100 ns using semi-isotropic Parrinello-Rahman barostat with a 5-ps coupling constant and using a 2-fs time step. The last frame of step (v) was used to spawn stress-free production runs, stress-strain calculations, and umbrella sampling simulations.

Derivation of the reaction coordinate

We carried out 20 independent, 1-μs long equilibrium simulations of the single-PF system in GTP-state, where the starting structure for each simulation was drawn every 150 ns from a ‘seeding’ simulation trajectory of 3 μs. For the single-PF system in GDP-state, we carried out 10 independent simulations (1 μs each) with the starting configurations drawn every 300 ns from a 3-μs ‘seeding’ trajectory. We then extracted backbone atoms (N, Cα, C and O) and excluded flexible protein regions (α: 38-46, α: 278-284 and β: 276-284) from further analysis.

Partial least-squares (PLS) functional mode analysis [33, 34] was then applied to the combined simulation set (both GTP- and GDP-state) to derive the collective mode of motion that correlated best with the fluctuations of the axial periodic dimension Lz and had the largest variance in terms of molecular motion. The linear regression model was trained on the first half of the GTP data set (∼13 μs) and the second half of the GDP data set (∼7 μs), and the remaining halves were used for cross-validation. The cross-validation revealed that the ensemble-weighted collective mode (corresponds to the solution with one PLS component by construction) had correlation coefficients of 0.9 (training set) and 0.85 (validation set), hence yielding a robust representation of the conformational transition between the expanded GTP- and compacted GDP-state (Fig 2B). A visualization of this transition is shown in a supplementary movie (see S1 Material).

Normalized mutual information and confinement entropy

In theory, the mutual information (MI) between two stochastic quantities χ1 and χ2 is (1) where H(χi) = − ∫pi(χi)logpii is the entropy and pi(χi) is the probability density of χi (i = 1, 2). The joint entropy H(χ1, χ2) is defined similarly and requires knowledge of the joint probability density p12(χ1, χ2).

In practice, calculation of the MI is very sensitive to how the underlying probability densities are discretized. Too coarse-grained discretization leads to an underestimation and too detailed discretization leads to an overestimation of the MI. We therefore used the Jack Knifed estimate that is known to be a low bias estimate of the MI and robust to discretization bin size [70]. It is defined by substituting the entropy in Eq 1 with the following estimate , where is the entropy calculated by a straightforward discretization and is the same as but when leaving out bin value j, and N is the total number of bins. The confinement entropy used to estimate the conformational space ‘volume’ is then , whereas the normalized mutual information is defined as: (2)

Stress-strain simulations of isolated PFs

To measure the response of the single-PF systems to external axial strain, we let the prepared systems equilibrate under anisotropic pressure conditions Pxx = PyyPzz until convergence of Lz, where Pzz ranged from − 65 bar to + 200 bar. All equilibration simulations were run for at least 1 μs, and the last 200 ns were used for further analysis. Due to the pressure difference maintained by the barostat, the simulated system (both solute and solvent) was subjected to an axial force fz such that the net stress on the PF along the z-axis, σzz, is: (3) where P = (Pxx + Pyy)/2, Lx and Ly are the lateral dimensions of the simulation box, and Az is the PF cross-section area (see next section). The axial strain was computed as: (4) where Lz,eq is the mean axial periodic dimension. We also note that Pzz, P, and Lx,y, z are, generally speaking, stochastic quantities. Therefore, block averaging with five blocks per trajectory and basic error propagation rules were used to estimate the mean and standard deviation of εzz and σzz.

Calculation of per-dimer elastic strain energy and flexural rigidity

According to linear elasticity theory, the per-dimer energy stored in a GTP-PF subjected to an axial elastic deformation by work required to compress it to the equilibrium state of a GDP-PF is ΔGel = ΔgelV, where Δgel is the elastic energy density and V is the effective dimer volume. Using the generalized Hooke’s law, we estimated the elastic energy density as in which EGTP ≈ 0.89 GPa (see Fig 2 in the main text) and is the axial strain tensor component reflecting the difference in the equilibrium dimer lengths of GTP- and GDP-PFs (derived from stress-free simulations). We then calculated the effective dimer volume by requiring that the PF cross-section area Az matches the mass per PF unit length m, i.e. m = ρAzLz,eq, where m ≈ 100 kDa and ρ ≈ 1.41 g/cm3 is the mass density of globular proteins with molecular weights M > 30 kDa [71]. This yielded Az ≈ 14.2 nm2, which allowed us to directly compute the sought elastic strain energy ΔGel ≈ 28.9 kJ/mol ≈ 11.6 kBT at T = 300 K.

Following previous work [26], the flexural rigidity (or bending stiffness) of a long hollow cylindrical filament is a product of its axial elasticity modulus E and the second moment of the cross-sectional area , where Rin and Rout are the inner and outer radii of the cylinder, respectively. Using the estimate Rin ≈ 10.19 nm from [72] (for 14_3 type MTs) and requiring that the MT cross-section conforms with the mass per MT unit length, i.e. , we obtained an estimate for the outer radius Rout ≈ 12.95 nm and, hence, for the second moment I ≈ 1.31 × 10−32 m4. This value allowed us to directly compute the flexural rigidities of GTP- and GDP-MTs (see Fig 2 in the main text).

Estimating MT bending stiffness from previous experimental data

The experimental values for MT bending stiffnesses and the respective uncertainties shown in Fig 2D were calculated using inverse-variance weighting [73]. Given a set of independent measurements yi with variances , the consensus inverse-variance mean and standard deviation are given by and , where the weights . For GDP-MTs and GDP-MTs stabilized with Taxol, we used the E × I values estimated by quantifying thermal fluctuations of MTs, as summarized in [29] and [30]. As there were only few measurements of GMPCPP-MTs in the cited publications, we extended the set by considering further thermal fluctuation studies [46, 47, 74].

Code and supplementary data

All refined starting structures are provided as Supplementary Data Sets. Unless explicitly specified, all numerical calculations were carried out using Python 2.7 [75] and Cython [76].

Supporting information

S1 Material. Supplementary text that includes all supplementary figures and tables as well as detailed information on stress-strain calculations, estimation of the per-dimer elastic strain energy, umbrella sampling simulations, estimation of the relative lateral bond stability, and Bayesian inference of the joint free energy distribution for the three-PF system.

https://doi.org/10.1371/journal.pcbi.1008132.s001

(PDF)

S1 Movie. Animation showing the compaction transition derived from equilibrium simulations of the single-PF system in both GTP- and GDP-state (see Fig 2).

https://doi.org/10.1371/journal.pcbi.1008132.s002

(MOV)

S1 Data Set. Archive containing refined single-, double- and three-PF structures in both GTP- and GDP-state.

S1: Refined single-PF structure in the GTP state. S2: Refined single-PF structure in the GDP state. S3: Refined homotypic double-PF structure in the GTP state. S4: Refined homotypic double-PF structure in the GDP state. S5: Refined three-PF structure in the GTP state. S6: Refined three-PF structure in the GDP state.

https://doi.org/10.1371/journal.pcbi.1008132.s003

(ZIP)

Acknowledgments

We thank Rui Zhang (Washington University in St. Louis, USA) and Eva Nogales (UC Berkeley, USA) for insightful discussions and for kindly providing the microtubule cryo-EM reconstructions; Gregory Bubnis (UC San Francisco, USA) and Thomas Ullmann (MPI-BPC, Göttingen, Germany) for suggestions about free energy calculations and error estimation. Computational resources were provided by the North-German Supercomputing Alliance (Berlin/Göttingen, Germany) as well as by the Max Planck Computing and Data Facility and the Leibniz Supercomputing Centre (Garching, Germany).

References

  1. 1. Sosa H. and Milligan R. A. Three-dimensional Structure of ncd-decorated Microtubules Obtained by a Back-projection Method. J. Mol. Biol., 260(5):743–755, 1996. pmid:8709152
  2. 2. Nogales E., Whittaker M., Milligan R. A., and Downing K. H. High-Resolution Model of the Microtubule. Cell, 96(1):79–88, 1999. pmid:9989499
  3. 3. Nogales E., Wolf S. G., and Downing K. H. Structure of the αβ Tubulin Dimer by Electron Crystallography. Nature, 391(6663):199–203, 1998. pmid:9428769
  4. 4. Carlier M. F., Didry D., and Pantaloni D. Microtubule Elongation and Guanosine 5’-Triphosphate Hydrolysis. Role of Guanine Nucleotides in Microtubule Dynamics. Biochemistry, 26(14):4428–4437, 1987. pmid:3663597
  5. 5. Mitchison T. and Kirschner M. Dynamic Instability of Microtubule Growth. Nature, 312(5991):237–42, 1984. pmid:6504138
  6. 6. Dogterom M. and Yurke B. Measurement of the Force-Velocity Relation for Growing Microtubules. Science, 278(5339):856–60, 1997. pmid:9346483
  7. 7. Grishchuk E. L., Molodtsov M. I., Ataullakhanov F. I., and McIntosh J. R. Force Production by Disassembling Microtubules. Nature, 438(7066):384–388, 2005. pmid:16292315
  8. 8. Driver J. W., Geyer E. A., Bailey M. E., Rice L. M., and Asbury C. L. Direct Measurement of Conformational Strain Energy in Protofilaments Curling Outward from Disassembling Microtubule Tips. Elife, 6:e28433, 2017. pmid:28628007
  9. 9. Simon J. R. and Salmon E. D. The Structure of Microtubule Ends During the Elongation and Shortening Phases of Dynamic Instability Examined by Negative-Stain Electron Microscopy. J. Cell Sci., 96 (Pt 4):571–82, 1990. pmid:2283357
  10. 10. Mandelkow E. M., Mandelkow E., and Milligan R. A. Microtubule Dynamics and Microtubule Caps: A Time-resolved Cryo- Electron Microscopy Study. J. Cell Biol., 114(5):977–991, 1991. pmid:1874792
  11. 11. Melki R., Carlier M. F., Pantaloni D., and Timasheff S. N. Cold Depolymerization of Microtubules to Double Rings: Geometric Stabilization of Assemblies. Biochemistry, 28(23):9143–9152, 1989. pmid:2605248
  12. 12. Chrétien D., Fuller S. D., and Karsenti E. Structure of Growing Microtubule Ends: Two-dimensional sheets Close into Tubes at Variable Rates. J. Cell Biol., 129(5):1311–1328, 1995. pmid:7775577
  13. 13. Müller-Reichert T., Chretien D., Severin F., and Hyman A. A. Structural Changes at Microtubule Ends Accompanying GTP Hydrolysis: Information from a Slowly Hydrolyzable Analogue of GTP, Guanylyl (α,β)methylenediphosphonate. Proc. Natl. Acad. Sci., 95(7):3661–3666, 1998. pmid:9520422
  14. 14. Atherton J., Stouffer M., Francis F., and Moores C. A. Microtubule Architecture in vitro and in Cells Revealed by Cryo-Electron Tomography. Acta Crystallogr. Sect. D Struct. Biol., 74(6):572–584, 2018.
  15. 15. McIntosh J. R., O’Toole E., Morgan G., Austin J., Ulyanov E., Ataullakhanov F., and Gudimchuk N. Microtubules Grow by the Addition of Bent Guanosine Triphosphate Tubulin to the Tips of Curved Protofilaments. J. Cell Biol., 217(8):2691–2708, 2018. pmid:29794031
  16. 16. Wang H. W. and Nogales E. Nucleotide-Dependent Bending Flexibility of Tubulin Regulates Microtubule Assembly. Nature, 435(7044):911–915, 2005. pmid:15959508
  17. 17. Igaev M. and Grubmüller H. Microtubule Assembly Governed by Tubulin Allosteric Gain in Flexibility and Lattice Induced Fit. Elife, 7:e34353, 2018. pmid:29652248
  18. 18. Fedorov V. A., Orekhov P. S., Kholina E. G., Zhmurov A. A., Ataullakhanov F. I., Kovalenko I. B., and Gudimchuk N. B. Mechanical Properties of Tubulin Intra- and Inter-dimer Interfaces and Their Implications for Microtubule Dynamic Instability. PLOS Comp. Biol., 15(8):e1007327, 2019.
  19. 19. Hyman A. A., Chrétien D., Arnal I., and Wade R. H. Structural Changes Accompanying GTP Hydrolysis in Microtubules: Information from a Slowly Hydrolyzable Analogue Guanylyl-(alpha,beta)-methylene-diphosphonate. J. Cell Biol., 128(1-2):117–25, 1995. pmid:7822409
  20. 20. Alushin G. M., Lander G.C., Kellogg E. H., Zhang R., Baker D., and Nogales E. High-Resolution Microtubule Structures Reveal the Structural Transitions in αβ-Tubulin upon GTP Hydrolysis. Cell, 157(5):1117–1129, 2014. pmid:24855948
  21. 21. Zhang R., Alushin G. M., Brown A., and Nogales E. Mechanistic Origin of Microtubule Dynamic Instability and Its Modulation by EB Proteins. Cell, 162(4):849–859, 2015. pmid:26234155
  22. 22. Zhang R., LaFrance B., and Nogales E. Separating the Effects of Nucleotide and EB Binding on Microtubule Structure. Proc. Natl. Acad. Sci., 115(27):E6191–E6200, 2018. pmid:29915050
  23. 23. Manka S. W. and Moores C. A. The Role of Tubulin-Tubulin Lattice Contacts in the Mechanism of Microtubule Dynamic Instability. Nat. Struct. Mol. Biol., 25(7):607–615, 2018. pmid:29967541
  24. 24. Estévez-Gallego J., Josa-Prado F., Ku S., Buey R. M., Balaguer F. A., Prota A. E., Lucena-Agell D., Kamma-Lorger C., Yagi T., Iwamoto H., Duchesne L., Barasoain I., Steinmetz M. O., Chrétien D., Kamimura S., Díaz J. F., and Oliva M. A. Structural Model for Differential Cap Maturation at Growing Microtubule Ends. Elife, 9:e50155, 2020. pmid:32151315
  25. 25. Tong D. and Voth G. A. Microtubule Simulations Provide Insight into the Molecular Mechanism Underlying Dynamic Instability. Biophys. J., 118(12):2938–2951, 2020. pmid:32413312
  26. 26. Harris B. J., Ross J. L., and Hawkins T. L. Microtubule Seams Are Not Mechanically Weak Defects. Phys. Rev. E, 97(6):1–7, 2018.
  27. 27. Manka S. W. and Moores C. A. Microtubule Structure by Cryo-EM: Snapshots of Dynamic Instability. Essays Biochem., 62(6):737–751, 2018. pmid:30315096
  28. 28. Howes S. C, Geyer E. A., LaFrance B., Zhang R., Kellogg E. H., Westermann S., Rice L. M., and Nogales E. Structural Differences between Yeast and Mammalian Microtubules Revealed by Cryo-EM. J. Cell Biol., 216(9): 2669–2677, 2017. pmid:28652389
  29. 29. Hawkins T., Mirigian M., Yasar M. S., and Ross J. L. Mechanics of Microtubules. J. Biomech., 43(1):23–30, 2010. pmid:19815217
  30. 30. Memet E., Hilitski F., Morris M. A., Schwenger W. J., Dogic Z., and Mahadevan L. Microtubules Soften due to Cross-Sectional Flattening. Elife, 7:e34695, 2018. pmid:29856317
  31. 31. Igaev M., Kutzner C., Bock L. V., Vaiana A. C., and Grubmüller H. Automated Cryo-EM Structure Refinement Using Correlation-Driven Molecular Dynamics. Elife, 8:e43542, 2019. pmid:30829573
  32. 32. Wells D. B. and Aksimentiev A. Mechanical Properties of a Complete Microtubule Revealed through Molecular Dynamics Simulation. Biophys. J., 99(2):629–637, 2010. pmid:20643083
  33. 33. Hub J. S., De Groot B. L., and van der Spoel D. g_wham—a free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. J. Chem. Theory Comput., 6(12):3713–3720, 2010.
  34. 34. Krivobokova T., Briones R., Hub J. S., Munk A., and De Groot B. L. Partial Least-Squares Functional Mode Analysis: Application to the Membrane Proteins AQP1, Aqy1, and CLC-ec1. Biophys. J., 103(4):786–796, 2012. pmid:22947940
  35. 35. Kikumoto M., Kurachi M., Tosa V., and Tashiro H. Flexural Rigidity of Individual Microtubules Measured by a Buckling Force with Optical Traps. Biophys. J., 90(5):1687–1696, 2006. pmid:16339879
  36. 36. Ray S., Meyhöfer E., Milligan R. A., and Howard J. Kinesin Follows the Microtubule’s Protofilament Axis. J. Cell Biol., 121(5):1083–1093, 1993. pmid:8099076
  37. 37. Caplow M., Ruhlen R. L., and Shanks J. The Free Energy for Hydrolysis of a Microtubule-Bound Nucleotide Triphosphate Is Near Zero: All of the Free Energy for Hydrolysis Is Stored in the Microtubule Lattice. J. Cell Biol., 127(3):779–788, 1994. pmid:7962059
  38. 38. Chrétien D., Metoz F., Verde F., Karsenti E., and Wade R. H. Lattice Defects in Microtubules: Protofilament Numbers Vary within Individual Microtubules. J. Cell Biol., 117(5):1031–40, 1992. pmid:1577866
  39. 39. Chrétien D. and Fuller S.D. Microtubules Switch Occasionally into Unfavorable Configurations during Elongation. J. Mol. Biol., 298(4):663–676, 2000. pmid:10788328
  40. 40. Hunyadi V., Chrétien D., and Jánosi I. M. Mechanical Stress Induced Mechanism of Microtubule Catastrophes. J. Mol. Biol., 348(4):927–938, 2005. pmid:15843023
  41. 41. Brouhard G. J. and Rice L. M. Microtubule Dynamics: An Interplay of Biochemistry and Mechanics. Nat. Rev. Mol. Cell Biol., 19(7):451–463, 2018. pmid:29674711
  42. 42. Kim T. and Rice L. M. Long-Range, Through-Lattice Coupling Improves Predictions of Microtubule Catastrophe. Mol. Biol. Cell, 30(12):1451–1462, 2019. pmid:30943103
  43. 43. Yee A. W., Aldeghi M., Blakeley M. P., Ostermann A., Mas P. J., Moulin M., de Sanctis D., Bowler M. W., Mueller-Dieckmann C., Mitchell E. P., Haertlein M., de Groot B. L., Boeri Erba E., and Forsyth V. T. A Molecular Mechanism for Transthyretin Amyloidogenesis. Nat. Commun., 10(1):925, 2019. pmid:30804345
  44. 44. Hemmat M., Castle B. T., Sachs J. N., and Odde D. J. Multiscale Computational Modeling of Tubulin-Tubulin Lateral Interaction. Biophys. J., 117(7):1234–1249, 2019. pmid:31493861
  45. 45. Kellogg E. H., Hejab N. M.A., Howes S., Northcote P., Miller J. H., Díaz J. F., Downing K. H., and Nogales E. Insights into the Distinct Mechanisms of Action of Taxane and Non-Taxane Microtubule Stabilizers from Cryo-EM Structures. J. Mol. Biol., 429(5):633–646, 2017. pmid:28104363
  46. 46. Mickey B. and Howard J. Rigidity of Microtubules Is Increased by Stabilizing Agents. J. Cell Biol., 130(4):909–917, 1995. pmid:7642706
  47. 47. Hawkins T. L., Sept D., Mogessie B., Straube A., and Ross J. L. Mechanical Properties of Doubly Stabilized Microtubule Filaments. Biophys. J., 104(7):1517–1528, 2013. pmid:23561528
  48. 48. Portran D., Schaedel L., Xu Z., Théry M., and Nachury M. V. Tubulin Acetylation Protects Long-Lived Microtubules Against Mechanical Ageing. Nat. Cell Biol., 19(4):391–398, 2017. pmid:28250419
  49. 49. Xu Z., Schaedel L., Portran D., Aguilar A., Gaillard J., Marinkovich M. P., Théry M., and Nachury M. V. Microtubules Acquire Resistance from Mechanical Breakage through Intralumenal Acetylation. Science, 356(6335):328–332, 2017. pmid:28428427
  50. 50. Eshun-Wilson L., Zhang R., Portran D., Nachury M. V., Toso D. B., Löhr T., Vendruscolo M., Bonomi M., Fraser J. S., and Nogales E. Effects of α-Tubulin Acetylation on Microtubule Structure and Stability. Proc. Natl. Acad. Sci., 116(21):10366–10371, 2019. pmid:31072936
  51. 51. Geyer E. A., Burns A., Lalonde B. A., Ye X., Piedra F.-A., Huffaker T. C., and Rice L. M. A Mutation Uncouples the Tubulin Conformational and GTPase Cycles, Revealing Allosteric Control of Microtubule Dynamics. Elife, 4:e10113, 2015. pmid:26439009
  52. 52. Roostalu J., Thomas C., Cade N. I., Kunzelmann S., Taylor I. A., and Surrey T. The Speed of GTP Hydrolysis Determines GTP Cap Size and Controls Microtubule Stability. Elife, 9:e51992, 2020. pmid:32053491
  53. 53. Chu J.-W. and Voth G. A. Allostery of Actin Filaments: Molecular Dynamics Simulations and Coarse-Grained Analysis. Proc. Natl. Acad. Sci., 102(37):13111–13116, 2005. pmid:16135566
  54. 54. Luan B. and Aksimentiev A. Strain Softening in Stretched DNA. Phys. Rev. Lett., 101(11):118101, 2008. pmid:18851334
  55. 55. Bochicchio D., Kwangmettatam S., Kudernac T., and Pavan G. M. How Defects Control the Out-of-Equilibrium Dissipative Evolution of a Supramolecular Tubule. ACS Nano, 13(4):4322–4334, 2019. pmid:30875196
  56. 56. Piana S., Lindorff-Larsen K., and Shaw D. E. How Robust Are Protein Folding Simulations with Respect to Force Field Parameterization? Biophys. J., 100(9):L47–L49, 2011. pmid:21539772
  57. 57. MacKerell A. D., Bashford D., Bellott M., Dunbrack R. L., Evanseck J. D., Field M. J., Fischer S., Gao J., Guo H., Ha S., Joseph-McCarthy D., Kuchnir L., Kuczera K., Lau F. T. K., Mattos C., Michnick S., Ngo T., Nguyen D. T., Prodhom B., Reiher W. E., Roux B., Schlenkrich M., Smith J. C., Stote R., Straub J., Watanabe M., Wiórkiewicz-Kuczera J., Yin D., and Karplus M. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B, 102(18):3586–3616, 1998. pmid:24889800
  58. 58. Mackerell A. D., Feig M., and Brooks C. L. Extending the Treatment of Backbone Energetics in Protein Force Fields: Limitations of Gas-Phase Quantum Mechanics in Reproducing Protein Conformational Distributions in Molecular Dynamics Simulations. J. Comput. Chem., 25(11):1400–1415, 2004. pmid:15185334
  59. 59. Ullmann R. T. and Ullmann G. M. GMCT: A Monte Carlo Simulation Package for Macromolecular Receptors. J. Comput. Chem., 33(8):887–900, 2012. pmid:22278916
  60. 60. Fiser A., Do R. K., and Sali A. Modeling of Loops in Protein Structures. Protein Sci., 9(9):1753–73, 2000. pmid:11045621
  61. 61. Pettersen E. F., Goddard T. D., Huang C. C., Couch G. S., Greenblatt D. M., Meng E. C., and Ferrin T. E. UCSF Chimera—A Visualization System for Exploratory Research and Analysis. J. Comput. Chem., 25(13):1605–1612, 2004. pmid:15264254
  62. 62. Humphrey W., Dalke A., and Schulten K. VMD: Visual Molecular Dynamics. J. Mol. Graph., 14(1):33–38, 1996. pmid:8744570
  63. 63. Abraham M. J., Murtola T., Schulz R., Páll S., Smith J. C., Hess B., and Lindahl E. GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX, 1-2:19–25, 2015.
  64. 64. Essmann U., Perera L., Berkowitz M. L., Darden T., Lee H., and Pedersen L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys., 103(19):8577–8593, 1995.
  65. 65. Hess B., Kutzner C., van der Spoel D., and Lindahl E. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput., 4(3):435–47, 2008. pmid:26620784
  66. 66. Bussi G., Donadio D., and Parrinello M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys., 126(1):014101, 2007. pmid:17212484
  67. 67. Feenstra K. A., Hess B., and Berendsen H. J. C. Improving Efficiency of Large Time-Scale Molecular Dynamics Simulations of Hydrogen-Rich Systems. J. Comput. Chem., 20(8):786–798, 1999.
  68. 68. Berendsen H. J. C., Postma J. P. M., van Gunsteren W. F., DiNola A., and Haak J. R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys., 81(8):3684–3690, 1984.
  69. 69. Parrinello M. and Rahman A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys., 52(12):7182–7190, 1981.
  70. 70. Paninski L. Estimation of Entropy and Mutual Information. Neural Comput., 15(6):1191–1253, 2003.
  71. 71. Fischer H., Polikarpov I., and Craievich A. F. Average Protein Density Is a Molecular-Weight-Dependent Function. Protein Sci., 13(10):2825–2828, 2009.
  72. 72. Donhauser Z. J., Jobs W. B., and Binka E. C. Mechanics of Microtubules: Effects of Protofilament Orientation. Biophys. J., 99(5):1668–1675, 2010. pmid:20816081
  73. 73. Hartung J., Knapp G., and Sinha B. K. Statistical Meta-Analysis with Applications. Wiley Series in Probability and Statistics. John Wiley & Sons, Inc., 2008.
  74. 74. Kawaguchi K. and Yamaguchi A. Temperature Dependence Rigidity of Non-taxol Stabilized Single Microtubules. Biochem. Biophys. Res. Commun., 402(1):66–69, 2010. pmid:20920471
  75. 75. Python Software Foundation. Python Language Reference, version 2.7. February 6, 2017, 2017.
  76. 76. Behnel S., Bradshaw R., Citro C., Dalcin L., Seljebotn D. S., and Smith K. Cython: The Best of Both Worlds. Comput. Sci. Eng., 13(2):31–39, 2011.