Functional Dissection and Hierarchy of Tubulin-folding Cofactor Homologues in Fission Yeast

We describe the isolation of fission yeast homologues of tubulin-folding cofactors B (Alp11) and E (Alp21), which are essential for cell viability and the maintenance of microtubules. Alp11^B contains the glycine-rich motif (the CLIP-170 domain) involved in microtubular functions, whereas, unlike mammalian cofactor E, Alp21^E does not. Both mammalian and yeast cofactor E, however, do contain leucine-rich repeats. Immunoprecipitation analysis shows that Alp11^B interacts with both α-tubulin and Alp21^E, but not with the cofactor D homologue Alp1, whereas Alp21^E also interacts with Alp1^D. The cellular amount of α-tubulin is decreased in both alp1 and alp11 mutants. Overproduction of Alp11^B results in cell lethality and the disappearance of microtubules, which is rescued by co-overproduction of α-tubulin. Both full-length Alp11^B and the C-terminal third containing the CLIP-170 domain localize in the cytoplasm, and this domain is required for efficient binding to α-tubulin. Deletion of alp11 is suppressed by multicopy plasmids containing either alp21⁺ or alp1⁺, whereas alp21 deletion is rescued by overexpression of alp1⁺ but not alp11⁺. Finally, the alp1 mutant is not complemented by either alp11⁺ or alp21⁺. The results suggest that cofactors operate in a linear pathway (Alp11^B-Alp21^E-Alp1^D), each with distinct roles.     

Tubulin Subunits Exist in an Activated Conformational State Generated and Maintained by Protein Cofactors

The production of native α/β tubulin heterodimer in vitro depends on the action of cytosolic chaperonin and several protein cofactors. We previously showed that four such cofactors (termed A, C, D, and E) together with native tubulin act on β-tubulin folding intermediates generated by the chaperonin to produce polymerizable tubulin heterodimers. However, this set of cofactors generates native heterodimers only very inefficiently from α-tubulin folding intermediates produced by the same chaperonin. Here we describe the isolation, characterization, and genetic analysis of a novel tubulin folding cofactor (cofactor B) that greatly enhances the efficiency of α-tubulin folding in vitro. This enabled an integrated study of α- and β-tubulin folding: we find that the pathways leading to the formation of native α- and β-tubulin converge in that the folding of the α subunit requires the participation of cofactor complexes containing the β subunit and vice versa. We also show that sequestration of native α-or β-tubulins by complex formation with cofactors results in the destabilization and decay of the remaining free subunit. These data demonstrate that tubulin folding cofactors function by placing and/or maintaining α-and β-tubulin polypeptides in an activated conformational state required for the formation of native α/β heterodimers, and imply that each subunit provides information necessary for the proper folding of the other.     

Mzt1/Tam4, a fission yeast MOZART1 homologue,is an essential component of the γ-tubulin complex and directly interacts with GCP3Alp6

In humans, MOZART1 plays an essential role in mitotic spindle formation as a component of the γ-tubulin ring complex. We report that the fission yeast homologue of MOZART1, Mzt1/Tam4, is located at microtubule-organizing centers (MTOCs) and coimmunoprecipitates with γ-tubulin Gtb1 from cell extracts. We show that mzt1/tam4 is an essential gene in fission yeast, encoding a 64–amino acid peptide, depletion of which leads to aberrant microtubule structure, including malformed mitotic spindles and impaired interphase microtubule array. Mzt1/Tam4 depletion also causes cytokinesis defects, suggesting a role of the γ-tubulin complex in the regulation of cytokinesis. Yeast two-hybrid analysis shows that Mzt1/Tam4 forms a complex with Alp6, a fission yeast homologue of γ-tubulin complex protein 3 (GCP3). Biophysical methods demonstrate that there is a direct interaction between recombinant Mzt1/Tam4 and the N-terminal region of GCP3^Alp6. Together our results suggest that Mzt1/Tam4 contributes to the MTOC function through regulation of GCP3^Alp6.     

Dynamic Subunit Exchange and the Regulation of Microtubule Assembly by the Stress Response Protein Human αB Crystallin

Background

The small heat shock protein (sHSP), human αB crystallin, forms large, polydisperse complexes that modulate the tubulin-microtubule equilibrium using a dynamic mechanism that is poorly understood. The interactive sequences in αB crystallin for tubulin are surface exposed, and correspond to interactive sites for the formation of αB crystallin complexes.

Methodology/Principal Findings

There is sequence homology between tubulin and the interactive domains in the β8-strand of the core domain and the C-terminal extension of αB crystallin. This study investigated the hypothesis that the formation of tubulin and αB crystallin quaternary structures was regulated through common interactive domains that alter the dynamics of their assembly. Size exclusion chromatography (SEC), SDS-PAGE, microtubule assembly assays, aggregation assays, multiple sequence alignment, and molecular modeling characterized the dynamic response of tubulin assembly to increasing concentrations of αB crystallin. Low molar ratios of αB crystallin∶tubulin were favorable for microtubule assembly and high molar ratios of αB crystallin∶tubulin were unfavorable for microtubule assembly. Interactions between αB crystallin and unassembled tubulin were observed using SEC and SDS-PAGE.

Conclusions/Significance

Subunits of αB crystallin that exchange dynamically with the αB crystallin complex can interact with tubulin subunits to regulate the equilibrium between tubulin and microtubules.     

The Free Energy Profile of Tubulin Straight-Bent Conformational Changes,with Implications for Microtubule Assembly and Drug Discovery

αβ-tubulin dimers need to convert between a ‘bent’ conformation observed for free dimers in solution and a ‘straight’ conformation required for incorporation into the microtubule lattice. Here, we investigate the free energy landscape of αβ-tubulin using molecular dynamics simulations, emphasizing implications for models of assembly, and modulation of the conformational landscape by colchicine, a tubulin-binding drug that inhibits microtubule polymerization. Specifically, we performed molecular dynamics, potential-of-mean force simulations to obtain the free energy profile for unpolymerized GDP-bound tubulin as a function of the ∼12° intradimer rotation differentiating the straight and bent conformers. Our results predict that the unassembled GDP-tubulin heterodimer exists in a continuum of conformations ranging between straight and bent, but, in agreement with existing structural data, suggests that an intermediate bent state has a lower free energy (by ∼1 kcal/mol) and thus dominates in solution. In agreement with predictions of the lattice model of microtubule assembly, lateral binding of two αβ-tubulins strongly shifts the conformational equilibrium towards the straight state, which is then ∼1 kcal/mol lower in free energy than the bent state. Finally, calculations of colchicine binding to a single αβ-tubulin dimer strongly shifts the equilibrium toward the bent states, and disfavors the straight state to the extent that it is no longer thermodynamically populated.     

Detailed Per-residue Energetic Analysis Explains the Driving Force for Microtubule Disassembly

Microtubules are long filamentous hollow cylinders whose surfaces form lattice structures of αβ-tubulin heterodimers. They perform multiple physiological roles in eukaryotic cells and are targets for therapeutic interventions. In our study, we carried out all-atom molecular dynamics simulations for arbitrarily long microtubules that have either GDP or GTP molecules in the E-site of β-tubulin. A detailed energy balance of the MM/GBSA inter-dimer interaction energy per residue contributing to the overall lateral and longitudinal structural stability was performed. The obtained results identified the key residues and tubulin domains according to their energetic contributions. They also identified the molecular forces that drive microtubule disassembly. At the tip of the plus end of the microtubule, the uneven distribution of longitudinal interaction energies within a protofilament generates a torque that bends tubulin outwardly with respect to the cylinder''s axis causing disassembly. In the presence of GTP, this torque is opposed by lateral interactions that prevent outward curling, thus stabilizing the whole microtubule. Once GTP hydrolysis reaches the tip of the microtubule (lateral cap), lateral interactions become much weaker, allowing tubulin dimers to bend outwards, causing disassembly. The role of magnesium in the process of outward curling has also been demonstrated. This study also showed that the microtubule seam is the most energetically labile inter-dimer interface and could serve as a trigger point for disassembly. Based on a detailed balance of the energetic contributions per amino acid residue in the microtubule, numerous other analyses could be performed to give additional insights into the properties of microtubule dynamic instability.     

Clostridium difficile Toxin A Decreases Acetylation of Tubulin,Leading to Microtubule Depolymerization through Activation of Histone Deacetylase 6, and This Mediates Acute Inflammation

Clostridium difficile toxin A is known to cause actin disaggregation through the enzymatic inactivation of intracellular Rho proteins. Based on the rapid and severe cell rounding of toxin A-exposed cells, we speculated that toxin A may be involved in post-translational modification of tubulin, leading to microtubule instability. In the current study, we observed that toxin A strongly reduced α-tubulin acetylation in human colonocytes and mouse intestine. Fractionation analysis demonstrated that toxin A-induced α-tubulin deacetylation yielded monomeric tubulin, indicating the presence of microtubule depolymerization. Inhibition of the glucosyltransferase activity against Rho proteins of toxin A by UDP-2′,3′-dialdehyde significantly abrogated toxin A-induced α-tubulin deacetylation. In colonocytes treated with trichostatin A (TSA), an inhibitor of the HDAC6 tubulin deacetylase, toxin A-induced α-tubulin deacetylation and loss of tight junction were completely blocked. Administration of TSA also attenuated proinflammatory cytokine production, mucosal damage, and epithelial cell apoptosis in mouse intestine exposed to toxin A. These results suggest that toxin A causes microtubule depolymerization by activation of HDAC6-mediated tubulin deacetylation. Indeed, blockage of HDAC6 by TSA markedly attenuates α-tubulin deacetylation, proinflammatory cytokine production, and mucosal damage in a toxin A-induced mouse enteritis model. Tubulin deacetylation is an important component of the intestinal inflammatory cascade following toxin A-mediated Rho inactivation in vitro and in vivo.     

TBCD Links Centriologenesis,Spindle Microtubule Dynamics,and Midbody Abscission in Human Cells

Microtubule-organizing centers recruit α- and β-tubulin polypeptides for microtubule nucleation. Tubulin synthesis is complex, requiring five specific cofactors, designated tubulin cofactors (TBCs) A–E, which contribute to various aspects of microtubule dynamics in vivo. Here, we show that tubulin cofactor D (TBCD) is concentrated at the centrosome and midbody, where it participates in centriologenesis, spindle organization, and cell abscission. TBCD exhibits a cell-cycle-specific pattern, localizing on the daughter centriole at G1 and on procentrioles by S, and disappearing from older centrioles at telophase as the protein is recruited to the midbody. Our data show that TBCD overexpression results in microtubule release from the centrosome and G1 arrest, whereas its depletion produces mitotic aberrations and incomplete microtubule retraction at the midbody during cytokinesis. TBCD is recruited to the centriole replication site at the onset of the centrosome duplication cycle. A role in centriologenesis is further supported in differentiating ciliated cells, where TBCD is organized into “centriolar rosettes”. These data suggest that TBCD participates in both canonical and de novo centriolar assembly pathways.     

Tubulin Acetyltransferase αTAT1 Destabilizes Microtubules Independently of Its Acetylation Activity

Acetylation of α-tubulin at lysine 40 (K40) is a well-conserved posttranslational modification that marks long-lived microtubules but has poorly understood functional significance. Recently, αTAT1, a member of the Gcn5-related N-acetyltransferase superfamily, has been identified as an α-tubulin acetyltransferase in ciliated organisms. Here, we explored the function of αTAT1 with the aim of understanding the consequences of αTAT1-mediated microtubule acetylation. We demonstrate that α-tubulin is the major target of αTAT1 but that αTAT1 also acetylates itself in a regulatory mechanism that is required for effective modification of tubulin. We further show that in mammalian cells, αTAT1 promotes microtubule destabilization and accelerates microtubule dynamics. Intriguingly, this effect persists in an αTAT1 mutant with no acetyltransferase activity, suggesting that interaction of αTAT1 with microtubules, rather than acetylation per se, is the critical factor regulating microtubule stability. Our data demonstrate that αTAT1 has cellular functions that extend beyond its classical enzymatic activity as an α-tubulin acetyltransferase.     

C-Terminal Tail Polyglycylation and Polyglutamylation Alter Microtubule Mechanical Properties

Microtubules are biopolymers that perform diverse cellular functions. Microtubule behavior regulation occurs in part through post-translational modification of both the α- and β-subunits of tubulin. One class of modifications is the heterogeneous addition of glycine and/or glutamate residues to the disordered C-terminal tails (CTTs) of tubulin. Because of their prevalence in stable, high-stress cellular structures such as cilia, we sought to determine if these modifications alter microtubules’ intrinsic stiffness. Here, we describe the purification and characterization of differentially modified pools of tubulin from Tetrahymena thermophila. We found that post-translational modifications do affect microtubule stiffness but do not affect the number of protofilaments incorporated into microtubules. We measured the spin dynamics of nuclei in the CTT backbone by NMR spectroscopy to explore the mechanism of this change. Our results show that the α-tubulin CTT does not protrude out from the microtubule surface, as is commonly depicted in models, but instead interacts with the dimer’s surface. This suggests that the interactions of the α-tubulin CTT with the tubulin body contributes to the stiffness of the assembled microtubule, thus providing insight into the mechanism by which polyglycylation and polyglutamylation can alter microtubule mechanical properties.     

Centrosomes Isolated from Spisula solidissima Oocytes Contain Rings and an Unusual Stoichiometric Ratio of α/β Tubulin

Centrosome-dependent microtubule nucleation involves the interaction of tubulin subunits with pericentriolar material. To study the biochemical and structural basis of centrosome-dependent microtubule nucleation, centrosomes capable of organizing microtubules into astral arrays were isolated from parthenogenetically activated Spisula solidissima oocytes. Intermediate voltage electron microscopy tomography revealed that each centrosome was composed of a single centriole surrounded by pericentriolar material that was studded with ring-shaped structures ~25 nm in diameter and <25 nm in length. A number of proteins copurified with centrosomes including: (a) proteins that contained M-phase–specific phosphoepitopes (MPM-2), (b) α-, β-, and γ-tubulins, (c) actin, and (d) three low molecular weight proteins of <20 kD. γ-Tubulin was not an MPM-2 phosphoprotein and was the most abundant form of tubulin in centrosomes. Relatively little α- or β-tubulin copurified with centrosomes, and the ratio of α- to β-tubulin in centrosomes was not 1:1 as expected, but rather 1:4.6, suggesting that centrosomes contain β-tubulin that is not dimerized with α-tubulin.     

Conformational changes in tubulin upon binding cryptophycin-52 reveal its mechanism of action

Cryptophycin-52 (Cp-52) is potentially the most potent anticancer drug known, with IC₅₀ values in the low picomolar range, but its binding site on tubulin and mechanism of action are unknown. Here, we have determined the binding site of Cp-52, and its parent compound, cryptophycin-1, on HeLa tubulin, to a resolution of 3.3 Å and 3.4 Å, respectively, by cryo-EM and characterized this binding further by molecular dynamics simulations. The binding site was determined to be located at the tubulin interdimer interface and partially overlap that of maytansine, another cytotoxic tubulin inhibitor. Binding induces curvature both within and between tubulin dimers that is incompatible with the microtubule lattice. Conformational changes occur in both α-tubulin and β-tubulin, particularly in helices H8 and H10, with distinct differences between α and β monomers and between Cp-52-bound and cryptophycin-1-bound tubulin. From these results, we have determined: (i) the mechanism of action of inhibition of both microtubule polymerization and depolymerization, (ii) how the affinity of Cp-52 for tubulin may be enhanced, and (iii) where linkers for targeted delivery can be optimally attached to this molecule.     

Fission yeast Alp14 is a dose-dependent plus end-tracking microtubule polymerase

XMAP215/Dis1 proteins are conserved tubulin-binding TOG-domain proteins that regulate microtubule (MT) plus-end dynamics. Here we show that Alp14, a XMAP215 orthologue in fission yeast, Schizosaccharomyces pombe, has properties of a MT polymerase. In vivo, Alp14 localizes to growing MT plus ends in a manner independent of Mal3 (EB1). alp14-null mutants display short interphase MTs with twofold slower assembly rate and frequent pauses. Alp14 is a homodimer that binds a single tubulin dimer. In vitro, purified Alp14 molecules track growing MT plus ends and accelerate MT assembly threefold. TOG-domain mutants demonstrate that tubulin binding is critical for function and plus end localization. Overexpression of Alp14 or only its TOG domains causes complete MT loss in vivo, and high Alp14 concentration inhibits MT assembly in vitro. These inhibitory effects may arise from Alp14 sequestration of tubulin and effects on the MT. Our studies suggest that Alp14 regulates the polymerization state of tubulin by cycling between a tubulin dimer-bound cytoplasmic state and a MT polymerase state that promotes rapid MT assembly.     

Effects of tubulin acetylation and tubulin acetyltransferase binding on microtubule structure

Tubulin undergoes posttranslational modifications proposed to specify microtubule subpopulations for particular functions. Most of these modifications occur on the C-termini of tubulin and may directly affect the binding of microtubule-associated proteins (MAPs) or motors. Acetylation of Lys-40 on α-tubulin is unique in that it is located on the luminal surface of microtubules, away from the interaction sites of most MAPs and motors. We investigate whether acetylation alters the architecture of microtubules or the conformation of tubulin, using cryo–electron microscopy (cryo-EM). No significant changes are observed based on protofilament distributions or microtubule helical lattice parameters. Furthermore, no clear differences in tubulin structure are detected between cryo-EM reconstructions of maximally deacetylated or acetylated microtubules. Our results indicate that the effect of acetylation must be highly localized and affect interaction with proteins that bind directly to the lumen of the microtubule. We also investigate the interaction of the tubulin acetyltransferase, αTAT1, with microtubules and find that αTAT1 is able to interact with the outside of the microtubule, at least partly through the tubulin C-termini. Binding to the outside surface of the microtubule could facilitate access of αTAT1 to its luminal site of action if microtubules undergo lateral opening between protofilaments.     

Activation of Microtubule Dynamics Increases Neuronal Growth via the Nerve Growth Factor (NGF)- and Gαs-mediated Signaling Pathways

Signals that activate the G protein Gα_s and promote neuronal differentiation evoke Gα_s internalization in rat pheochromocytoma (PC12) cells. These agents also significantly increase Gα_s association with microtubules, resulting in an increase in microtubule dynamics because of the activation of tubulin GTPase by Gα_s. To determine the function of Gα_s/microtubule association in neuronal development, we used real-time trafficking of a GFP-Gα_s fusion protein. GFP-Gα_s concentrates at the distal end of the neurites in differentiated living PC12 cells as well as in cultured hippocampal neurons. Gα_s translocates to specialized membrane compartments at tips of growing neurites. A dominant-negative Gα chimera that interferes with Gα_s binding to tubulin and activation of tubulin GTPase attenuates neurite elongation and neurite number both in PC12 cells and primary hippocampal neurons. This effect is greatest on differentiation induced by activated Gα_s. Together, these data suggest that activated Gα_s translocates from the plasma membrane and, through interaction with tubulin/microtubules in the cytosol, is important for neurite formation, development, and outgrowth. Characterization of neuronal G protein dynamics and their contribution to microtubule dynamics is important for understanding the molecular mechanisms by which G protein-coupled receptor signaling orchestrates neuronal growth and differentiation.     

Orientation and structure of the Ndc80 complex on the microtubule lattice

The four-subunit Ndc80 complex, comprised of Ndc80/Nuf2 and Spc24/Spc25 dimers, directly connects kinetochores to spindle microtubules. The complex is anchored to the kinetochore at the Spc24/25 end, and the Ndc80/Nuf2 dimer projects outward to bind to microtubules. Here, we use cryoelectron microscopy and helical image analysis to visualize the interaction of the Ndc80/Nuf2 dimer with microtubules. Our results, when combined with crystallography data, suggest that the globular domain of the Ndc80 subunit binds strongly at the interface between tubulin dimers and weakly at the adjacent intradimer interface along the protofilament axis. Such a binding mode, in which the Ndc80 complex interacts with sequential α/β-tubulin heterodimers, may be important for stabilizing kinetochore-bound microtubules. Additionally, we define the binding of the Ndc80 complex relative to microtubule polarity, which reveals that the microtubule interaction surface is at a considerable distance from the opposite kinetochore-anchored end; this binding geometry may facilitate polymerization and depolymerization at kinetochore-attached microtubule ends.     

Structural basis of tubulin tyrosination by tubulin tyrosine ligase

Tubulin tyrosine ligase (TTL) catalyzes the post-translational retyrosination of detyrosinated α-tubulin. Despite the indispensable role of TTL in cell and organism development, its molecular mechanism of action is poorly understood. By solving crystal structures of TTL in complex with tubulin, we here demonstrate that TTL binds to the α and β subunits of tubulin and recognizes the curved conformation of the dimer. Biochemical and cellular assays revealed that specific tubulin dimer recognition controls the activity of the enzyme, and as a consequence, neuronal development. The TTL–tubulin structure further illustrates how the enzyme binds the functionally crucial C-terminal tail sequence of α-tubulin and how this interaction catalyzes the tyrosination reaction. It also reveals how TTL discriminates between α- and β-tubulin, and between different post-translationally modified forms of α-tubulin. Together, our data suggest that TTL has specifically evolved to recognize and modify tubulin, thus highlighting a fundamental role of the evolutionary conserved tubulin tyrosination cycle in regulating the microtubule cytoskeleton.     

Structural Basis for the Association of MAP6 Protein with Microtubules and Its Regulation by Calmodulin

Microtubules are highly dynamic αβ-tubulin polymers. In vitro and in living cells, microtubules are most often cold- and nocodazole-sensitive. When present, the MAP6/STOP family of proteins protects microtubules from cold- and nocodazole-induced depolymerization but the molecular and structure determinants by which these proteins stabilize microtubules remain under debate. We show here that a short protein fragment from MAP6-N, which encompasses its Mn1 and Mn2 modules (MAP6(90–177)), recapitulates the function of the full-length MAP6-N protein toward microtubules, i.e. its ability to stabilize microtubules in vitro and in cultured cells in ice-cold conditions or in the presence of nocodazole. We further show for the first time, using biochemical assays and NMR spectroscopy, that these effects result from the binding of MAP6(90–177) to microtubules with a 1:1 MAP6(90–177):tubulin heterodimer stoichiometry. NMR data demonstrate that the binding of MAP6(90–177) to microtubules involve its two Mn modules but that a single one is also able to interact with microtubules in a closely similar manner. This suggests that the Mn modules represent each a full microtubule binding domain and that MAP6 proteins may stabilize microtubules by bridging tubulin heterodimers from adjacent protofilaments or within a protofilament. Finally, we demonstrate that Ca²⁺-calmodulin competes with microtubules for MAP6(90–177) binding and that the binding mode of MAP6(90–177) to microtubules and Ca²⁺-calmodulin involves a common stretch of amino acid residues on the MAP6(90–177) side. This result accounts for the regulation of microtubule stability in cold condition by Ca²⁺-calmodulin.     

The fission yeast transforming acidic coiled coil-related protein Mia1p/Alp7p is required for formation and maintenance of persistent microtubule-organizing centers at the nuclear envelope

Microtubule-organizing centers (MTOCs) concentrate microtubule nucleation, attachment and bundling factors and thus restrict formation of microtubule arrays in spatial and temporal manner. How MTOCs occur remains an exciting question in cell biology. Here, we show that the transforming acidic coiled coil–related protein Mia1p/Alp7p functions in emergence of large MTOCs in interphase fission yeast cells. We found that Mia1p was a microtubule-binding protein that preferentially localized to the minus ends of microtubules and was associated with the sites of microtubule attachment to the nuclear envelope. Cells lacking Mia1p exhibited less microtubule bundles. Microtubules could be nucleated and bundled but were frequently released from the nucleation sites in mia1Δ cells. Mia1p was required for stability of microtubule bundles and persistent use of nucleation sites both in interphase and postanaphase array dynamics. The γ-tubulin–rich material was not organized in large perinuclear or microtubule-associated structures in mia1Δ cells. Interestingly, absence of microtubules in dividing wild-type cells prevented appearance of large γ-tubulin–rich MTOC structures in daughters. When microtubule polymerization was allowed, MTOCs were efficiently assembled de novo. We propose a model where MTOC emergence is a self-organizing process requiring the continuous association of microtubules with nucleation sites.     

Analysis of the Strength of Interfacial Hydrogen Bonds between Tubulin Dimers Using Quantum Theory of Atoms in Molecules

Microtubules are key structural elements that, among numerous biological functions, maintain the cytoskeleton of the cell and have a major role in cell division, which makes them important cancer chemotherapy targets. Understanding the energy balance that brings tubulin dimers, the building blocks of microtubules, together to form a microtubule is especially important for revealing the mechanism of their dynamic instability. Several studies have been conducted to estimate various contributions to the free energy of microtubule formation. However, the hydrogen-bond contribution was not studied before as a separate component. In this work, we use concepts such as the quantum theory of atoms in molecules to estimate the per-residue strength of hydrogen bonds contributing to the overall stability that brings subunits together in pair of tubulin heterodimers, across both the longitudinal and lateral interfaces. Our study shows that hydrogen bonding plays a major role in the stability of tubulin systems. Several residues that are crucial to the binding of vinca alkaloids are shown to be strongly involved in longitudinal microtubule stabilization. This indicates a direct relation between the binding of these agents and the effect on the interfacial hydrogen-bonding network, and explains the mechanism of their action. Lateral contacts showed much higher stability than longitudinal ones (–462 ± 70 vs. –392 ± 59 kJ/mol), which suggests a dramatic lateral stabilization effect of the GTP cap in the β-subunit. The role of the M-loop in lateral stability in absence of taxol was shown to be minor. The B-lattice lateral hydrogen bonds are shown to be comparable in strength to the A-lattice ones (−462 ± 70 vs. –472 ± 46 kJ/mol). These findings establish the importance of hydrogen bonds to the stability of tubulin systems.