Abstract
Microtubules (MTs) are dynamic cytoskeletal filaments composed of α- and β-tubulin protein dimers. They are crucial for maintaining cell structure, facilitating intracellular transport, and ensuring proper chromosome segregation among other things. These biological functions are influenced by the dynamic instability of the MT plus-end tip. Recent simulations have discovered formation of protofilament (PF) clusters at the MT plus-end tip, but reliable extrapolation of PF cluster dynamics and detailed microscopic mechanism are still needed to understand their behavior thoroughly. In this work, we have constructed, from "bottom up," a relatively high-resolution coarse-grained (CG) molecular dynamics (MD) model for tubulins with 20 CG sites per tubulin monomer, performed extensive CG MD simulations on MT lattices with 8 and 40 layers of heterodimers, and conducted comprehensive atomistic-level analysis. Our findings demonstrate that, in both GTP and GDP states, PF clusters are stable up to tens of microseconds of CG MD simulation time during spontaneous outward bending relaxation. PF clustering is initiated by longitudinal relaxation, stabilized by residual lateral interaction in the PF clusters. This process is thermodynamically driven by intrinsic lattice instability. In longer microtubules, this instability accumulates and further facilitates PF bending and clustering at the plus-end tip, but it can also be released via lattice curvature and supertwist. GDP-MTs form more PF clusters than GTP-MT on average and undergo more lateral cleavage and faster bending relaxation due to weaker lateral interactions, which facilitates MT catastrophe. GTP-MT forms flatter and more rigid PF clusters that favor nucleotide addition. Our findings highlight the critical role of lattice instability in microtubule dynamics and offer new insights on the conformational variability of MT plus-end tips.