Microtubule assembly governed by tubulin allosteric gain in flexibility and lattice induced fit
Figures

Allosteric and lattice models of MT assembly represented by a three-dimensional thermodynamic cycle.
GTP- (orange) and GDP-dimers (blue) can be free (no box), MT-cap-bound (dotted box), or integrated into the MT body (solid box). The dimmed states denote energetically unfavored states.

PCA on the set of experimental structures.
(A) Projection of the set onto the plane constituted by the first and second largest-amplitude conformational modes. Each structure is represented by a point on this plane. Known experimental structures involved in this study (PDB ID: 3JAT, 3JAS, 5JQG, 5FNV, 4ZOL, 4ZHQ) as well as other structures widely referenced in literature (PDB ID: 1TUB, 1JFF, 1SA0) are highlighted. (B) Same data as in (A) but with straight (TUB) and kinked (TUB) subpopulations highlighted in green and pink, respectively. (C) Curvature differences between the subpopulations in (B) displayed in terms of structure. Shown are the extreme conformations along conformational mode 1 (x-axis). -tubulin is shown in gray and -tubulin in green or pink, depending on the subpopulation. The anchor point that does not move during the straight-to-kinked transition is marked with a back dashed circle/line.
-
Figure 2—source data 1
Coordinates of the PDB structures projected onto conformational modes 1 and 2.
- https://doi.org/10.7554/eLife.34353.005

A simplified 2D example of the principal component analysis (PCA).
Given a set of observations, each being a point on a plane described by only two coordinates (left panel), the method calculates ‘collective’ coordinates on this plane that account for as much of the variation of the data as possible. These coincide with the directions of the largest variance in the data set and are always orthogonal to each other. The new coordinates are termed principal components or conformational modes (right panel). In a general case, individual observations in the data set may be described by an arbitrary number of coordinates (e.g., atoms each described by three coordinates in case of atomic structures), and only a few conformational modes are usually enough to describe the variation in the entire data set.
Merged set of the analyzed PDB structures played in consecutive order.
https://doi.org/10.7554/eLife.34353.006Structural motion along conformational mode 1.
https://doi.org/10.7554/eLife.34353.007Structural motion along conformational mode 2.
https://doi.org/10.7554/eLife.34353.008
Simulated tubulin ensembles in two nucleotide states (left vs.right panel) plotted as in Figure 2A,B.
Experimental structures are shown as black squares, whereas the simulated ensembles are shown as green (started from straight conformation) or pink (started from kinked conformation) point clouds. Experimental structures that were used for simulations are indicated with larger colored squares and vertical gray lines.
-
Figure 3—source data 1
Coordinates of the merged GTP-tubulin ensemble started from a straight structure (PDB ID: 3JAT, two independent simulations) projected onto conformational modes 1 and 2 (left panel, green).
- https://doi.org/10.7554/eLife.34353.011
-
Figure 3—source data 2
Coordinates of the merged GTP-tubulin ensemble started from kinked structures (PDB IDs: 5JQG, 5FNV) projected onto conformational modes 1 and 2 (left panel, pink).
- https://doi.org/10.7554/eLife.34353.012
-
Figure 3—source data 3
Coordinates of the merged GDP-tubulin ensemble started from a straight structure (PDB ID: 3JAS, two independent simulations) projected onto conformational modes 1 and 2 (right panel, green).
- https://doi.org/10.7554/eLife.34353.013
-
Figure 3—source data 4
Coordinates of the merged GDP-tubulin ensemble started from kinked structures (PDB IDs: 4ZOL, 4ZHQ) projected onto conformational modes 1 and 2 (right panel, pink).
- https://doi.org/10.7554/eLife.34353.014

Convergence of the simulated ensembles assessed by overlaying point clouds originating from different independent simulations (two independent simulations per bound nucleotide per curvature state).
The data are plotted as in Figure 2A,B and Figure 3. Dynamics of the simulated dimers stared from the straight MT conformations (PDB ID: 3JAT and 3JAS) during the -ns equilibration phase are shown as zoomed-in regions (top; red points).

Free energy profiles derived from the projections of the simulated ensembles onto a plane, where the x-axis denotes the position along the ensemble separation RC and the y-axis denotes an orthogonal coordinate describing the largest-amplitude bending motion in the combined TUB + TUB ensemble.
Experimental structures used for the free simulations are highlighted with gray dashed lines and green (straight structures) and pink squares (kinked structures).
-
Figure 4—source data 1
Coordinates of the PDB structures projected onto the ensemble separation and common bending coordinates (GTP).
- https://doi.org/10.7554/eLife.34353.018
-
Figure 4—source data 2
Coordinates of the merged GTP-tubulin ensemble started from a straight structure (PDB ID: 3JAT, two independent simulations) projected onto the ensemble separation and common bending coordinates (GTP).
- https://doi.org/10.7554/eLife.34353.019
-
Figure 4—source data 3
Coordinates of the merged GTP-tubulin ensemble started from kinked structures (PDB ID: 5JQG, 5FNV) projected onto the ensemble separation and common bending coordinates (GTP).
- https://doi.org/10.7554/eLife.34353.020
-
Figure 4—source data 4
Coordinates of the PDB structures projected onto the ensemble separation and common bending coordinates (GDP).
- https://doi.org/10.7554/eLife.34353.021
-
Figure 4—source data 5
Coordinates of the merged GDP-tubulin ensemble started from a straight structure (PDB ID: 3JAS, two independent simulations) projected onto the ensemble separation and common bending coordinates (GDP).
- https://doi.org/10.7554/eLife.34353.022
-
Figure 4—source data 6
Coordinates of the merged GDP-tubulin ensemble started from kinked structures (PDB ID: 4ZOL, 4ZHQ) projected onto the ensemble separation and common bending coordinates (GDP).
- https://doi.org/10.7554/eLife.34353.023

A 2D example of the ensemble separation search.
Given a system with two energy minima in the 2D conformational space (state 1 and 2), the method aims at rotating the coordinate system such that the x-axis coincides with the axis along which the system would encounter the highest barrier while transitioning from state 1 to state 2 (black arrow). A poorly chosen coordinate (left column) may result in a misleading free energy profile with an underestimated energy barrier and, hence, wrong conclusions. In a general case, the conformational space in which the two minima (state 1 and 2) exist has a dimension much higher than two, and the search is performed in all of these dimensions.

Dot products of the three-dimensional ensemble separation vector with vectors of higher dimensions for the tubulin ensembles in two different nucleotide states.
With the increasing number of PCA eigenvectors in the search space, the ensemble separation vector becomes independent of the search space dimension. This yields the minimal search space in which the ensembles can be efficiently separated.
Structural motion along the ensemble separation RC.
https://doi.org/10.7554/eLife.34353.024
Nucleotide-dependent motions of tubulin bending.
(A) Molecular representation of the twist-bend and splay-bend tubulin bending motions derived from the TUB and TUB simulated ensembles. The anchor point around the H8 helix is indicated with dashed lines/circles. See also Figure 5—video 1 for a side-by-side comparison. (B) Same data as in Figure 4 (left) with the derived bending motions shown as solid lines overlayed onto the free energy distributions. Experimental structures used for the free simulations are highlighted with gray dashed lines and green (straight structures) and pink squares (kinked structures).

Functional mode analysis of the TUB and TUB ensembles.
(A) Straight (reference) and kinked tubulin structures aligned by -subunit. The RMSD to the reference increases as the tubulin conformation deviates from the straight conformation, irrespective of the bending direction. (B) Pearson correlation coefficients between model building (, red) and cross-validation (, cyan) sets as functions of the FMA basis dimension. The optimal number of FMA components is indicated with an arrow in both cases. (C) Overlays of the RMSD (black) and the FMA fit with the model building and cross-correlation parts shown in red and cyan, respectively.
Side-by-side comparison of the twist-bend (green) and splay-bend (pink) conformational motions.
https://doi.org/10.7554/eLife.34353.027
Free energy profiles along the bending RCs defined for the TB and SB bending motions calculated for GTP- and GDP-tubulin.
Smaller (larger) values on the x-axis correspond to straighter (more kinked) dimer conformations. Squares denote the positions of the starting PDB structures in the GTP- (black) and GDP-state (gray) used for the free MD simulations. Note that the scaling for both RCs does not necessarily coincide with that of the common bending RC in Figure 4, as the vectors defining the TB and SB motions and the common bending motion were derived differently (see Materials and methods).
-
Figure 6—source data 1
Free energy values (second and third columns, mean +/- SD) along the TB bending coordinate (first column) for GTP-tubulin.
- https://doi.org/10.7554/eLife.34353.030
-
Figure 6—source data 2
Free energy values (second and third columns, mean +/- SD) along the TB bending coordinate (first column) for GDP-tubulin.
- https://doi.org/10.7554/eLife.34353.031
-
Figure 6—source data 3
Free energy values (second and third columns, mean +/- SD) along the SB bending coordinate (first column) for GTP-tubulin.
- https://doi.org/10.7554/eLife.34353.032
-
Figure 6—source data 4
Free energy values (second and third columns, mean +/- SD) along the SB bending coordinate (first column) for GDP-tubulin.
- https://doi.org/10.7554/eLife.34353.033

Convergence properties and robustness of the error evaluation with respect to the method used.
(A) Unnormalized individual umbrella histograms, each derived from a restrained simulation of at least 120 ns or longer. (B) Error evaluation using a block average method (Zhu and Hummer, 2012). The errors are consistent with those obtained through Bayesian bootstrapping in Figure 6. Gray circles depict the points on which the restraining harmonic potentials were centered. Here, the reference points are located at the leftmost edge of the respective free energy profile (black circles). The same color-coding is used as in Figure 5 and Figure 6.

Dynamics and energetics of GTP-tubulin along the ensemble separation RC.
(A) Collective mode of motion along the ensemble separation RC represented as a linear interpolation between the TUB and TUB simulated ensembles. The largest-amplitude intrinsic rearrangements in -tubulin are indicated with dashed lines and involve the H7 and H8 helices as well as a loop connecting them. (B) Free energy landscape as a function of the ensemble separation RC. Orientation of this coordinate with respect to the TB and SB free energy basins is schematically shown as an insert. Experimental structures used for the free simulations are highlighted with green (straight structures) and pink squares (kinked structures).
-
Figure 7—source data 1
Free energy values (second and third columns, mean +/- SD) along the ensemble separation coordinate (first column).
- https://doi.org/10.7554/eLife.34353.036

Convergence properties and robustness of the error evaluation with respect to the method used.
(A) Unnormalized individual umbrella histograms, each derived from a restrained simulation of at least 120 ns or longer. The region around the barrier at nm was sampled more extensively. (B) Error evaluation using a block average method (Zhu and Hummer, 2012). The errors are consistent with those obtained through Bayesian bootstrapping in Figure 7B (, , and ). Gray circles depict the points on which the restraining harmonic potentials were centered. Here, the reference point is located at the rightmost edge of the free energy profile (black circle).

The proposed, combined model of MT assembly represented by a thermodynamic cycle.
Here, TB and SB denote not only the different modes of tubulin bending but also the belonging to the respective free energy basin, i.e. anchor point state.
Additional files
-
Transparent reporting form
- https://doi.org/10.7554/eLife.34353.038