Microtubule assembly and disassembly dynamics model: Exploring dynamic instability and identifying features of Microtubules’ Growth, Catastrophe, Shortening, and Rescue

Microtubules (MTs), a cellular structure element, exhibit dynamic instability and can switch stochastically from growth to shortening; but the factors that trigger these processes at the molecular level are not understood. We developed a 3D Microtubule Assembly and Disassembly DYnamics (MADDY) model, based upon a bead-per-monomer representation of the αβ-tubulin dimers forming an MT lattice, stabilized by the lateral and longitudinal interactions between tubulin subunits. The model was parameterized against the experimental rates of MT growth and shortening, and pushing forces on the Dam1 protein complex due to protofilaments splaying out. Using the MADDY model, we carried out GPU-accelerated Langevin simulations to access dynamic instability behavior. By applying Machine Learning techniques, we identified the MT characteristics that distinguish simultaneously all four kinetic states: growth, catastrophe, shortening, and rescue. At the cellular 25 μM tubulin concentration, the most important quantities are the MT length [Formula: see text] , average longitudinal curvature [Formula: see text] , MT tip width [Formula: see text] , total energy of longitudinal interactions in MT lattice [Formula: see text] , and the energies of longitudinal and lateral interactions required to complete MT to full cylinder [Formula: see text] and [Formula: see text]. At high 250 μM tubulin concentration, the most important characteristics are [Formula: see text] , [Formula: see text] , number of hydrolyzed αβ-tubulin dimers [Formula: see text] and number of lateral interactions per helical pitch [Formula: see text] in MT lattice, energy of lateral interactions in MT lattice [Formula: see text] , and energy of longitudinal interactions in MT tip [Formula: see text]. These results allow greater insights into what brings about kinetic state stability and the transitions between states involved in MT dynamic instability behavior.