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 dinners, 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/nnol), which suggests a dramatic lateral stabilization effect of the GTP cap in the beta-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.