Force Dependence of Velocity and Run Length of Kinesin-1, Kinesin-2 and Kinesin-5 Family Molecular Motors

Identification of a Novel Gene Signature Based on Kinesin Family Members to Predict Prognosis in Glioma

Background and Objectives: Extensive research indicates that the kinesin superfamily (KIFs) regulates tumor progression. Nonetheless, the potential prognostic and therapeutic role of KIFs in glioma has been limited. Materials and Methods: Four independent cohorts from The Cancer Genome Atlas (TCGA) database and the Chinese Glioma Genome Atlas (CGGA) database were generated into a large combination cohort for identification of the prognostic signature. Following that, systematic analyses of multi-omics data were performed to determine the differences between the two groups. In addition, IDH1 was selected for the differential expression analysis. Results: The signature consists of five KIFs (KIF4A, KIF26A, KIF1A, KIF13A, and KIF13B) that were successfully identified. Receiver operating characteristic (ROC) curves indicated the signature had a suitable performance in prognosis prediction with the promising predictive area under the ROC curve (AUC) values. We then explored the genomic features differences, including immune features and tumor mutation status between high- and low-risk groups, from which we found that patients in the high-risk group had a higher level of immune checkpoint modules, and IDH1 was identified mutated more frequently in the low-risk group. Results of gene set enrichment analysis (GSEA) analysis showed that the E2F target, mitotic spindle, EMT, G2M checkpoint, and TNFa signaling were significantly activated in high-risk patients, partially explaining the differential prognosis between the two groups. Moreover, we also verified the five signature genes in the Human Protein Atlas (HPA) database. Conclusion: According to this study, we were able to classify glioma patients based on KIFs in a novel way. More importantly, the discovered KIFs-based signature and related characteristics may serve as a candidate for stratification indicators in the future for gliomas.

A model for the catalytic activity of microtubule polymerases

A XMAP215/Stu2/Alp14 polymerase can catalyze processively the tubulin addition to the microtubule (MT) plus end. In this work, a model is proposed for the underlying molecular mechanism of the polymerase activity, where the polymerase can not only catalyze processively the tubulin addition to but also promote the tubulin removal from the MT plus end. Based on the model the dynamics of both the wild-type and mutant polymerases is studied theoretically, explaining consistently and well various available experimental data. To further test the model, predicted results are provided.

Energetic constraints on filament-mediated cell polarization

Cell polarization underlies many cellular processes, such as differentiation, migration, and budding. Many living cells, such as budding yeast and fission yeast, use cytoskeletal structures to actively transport proteins to one location on the membrane and create a high-density spot of membrane-bound proteins. Yet, the thermodynamic constraints on filament-based cell polarization remain unknown. We show by mathematical modeling that cell polarization requires detailed balance to be broken, and we quantify the free-energy cost of maintaining a polarized state of the cell. Our study reveals that detailed balance cannot only be broken via the active transport of proteins along filaments but also via a chemical modification cycle, allowing detailed balance to be broken by the shuttling of proteins between the filament, membrane, and cytosol. Our model thus shows that cell polarization can be established via two distinct driving mechanisms, one based on active transport and one based on nonequilibrium binding. Furthermore, the model predicts that the driven binding process dissipates orders of magnitude less free energy than the transport-based process to create the same membrane spot. Active transport along filaments may be sufficient to create a polarized distribution of membrane-bound proteins, but an additional chemical modification cycle of the proteins themselves is more efficient and less sensitive to the physical exclusion of proteins on the transporting filaments, providing insight in the design principles of the Pom1/Tea1/Tea4 system in fission yeast and the Cdc42 system in budding yeast.

Dynamics of kinesin motor proteins under longitudinal and sideways loads

The available single-molecule data showed that different species of N-terminal kinesin molecular motors have very different features on dependences of run length and dissociation rate upon longitudinal load acting on stalks of the motors. The prior single-molecule data for Loligo pealei kinesin-1 indicated that the sideways load has only a weak effect on the velocity, but even a small sideways load can cause a large reduction in the run length. However, these puzzling experimental data remain to be explained and the underlying physical mechanisms are unclear. Here, based on our proposed model we study analytically the dynamics of the N-terminal kinesin motors such as Loligo pealei kinesin-1, Drosophila kinesin-1, truncated kinesin-5/Eg5, truncated kinesin-12/Kif15, kinesin-2/Kif17 and kinesin-2/Kif3AB dimers under both longitudinal and sideways loads. The theoretical results explain quantitatively the available experimental data and provide predictions. The physical mechanism of different kinesin species showing very different features on the load-dependent dynamics and the physical mechanism of the effect of the sideways load on the dynamics are revealed.

Effect of Kinesin-5 Tail Domain on Motor Dynamics for Antiparallel Microtubule Sliding

Kinesin-5 motor consists of two pairs of heads and tail domains, which are situated at the opposite ends of a common stalk. The two pairs of heads can bind to two antiparallel microtubules (MTs) and move on the two MTs independently towards the plus ends, sliding apart the two MTs, which is responsible for chromosome segregation during mitosis. Prior experimental data showed that the tails of kinesin-5 Eg5 can modulate the dynamics of single motors and are critical for multiple motors to generate high steady forces to slide apart two antiparallel MTs. To understand the molecular mechanism of the tails modulating the ability of Eg5 motors, based on our proposed model the dynamics of the single Eg5 with the tails and that without the tails moving on single MTs is studied analytically and compared. Furthermore, the dynamics of antiparallel MT sliding by multiple Eg5 motors with the tails and that without the tails is studied numerically and compared. Both the analytical results for single motors and the numerical results for multiple motors are consistent with the available experimental data.

Studies of Conformational Changes of Tubulin Induced by Interaction with Kinesin Using Atomistic Molecular Dynamics Simulations

The transition between strong and weak interactions of the kinesin head with the microtubule, which is regulated by the change of the nucleotide state of the head, is indispensable for the processive motion of the kinesin molecular motor on the microtubule. Here, using all-atom molecular dynamics simulations, the interactions between the kinesin head and tubulin are studied on the basis of the available high-resolution structural data. We found that the strong interaction can induce rapid large conformational changes of the tubulin, whereas the weak interaction cannot. Furthermore, we found that the large conformational changes of the tubulin have a significant effect on the interaction of the tubulin with the head in the weak-microtubule-binding ADP state. The calculated binding energy of the ADP-bound head to the tubulin with the large conformational changes is only about half that of the tubulin without the conformational changes.

A model of processive walking and slipping of kinesin-8 molecular motors

Kinesin-8 molecular motor can move with superprocessivity on microtubules towards the plus end by hydrolyzing ATP molecules, depolymerizing microtubules. The available single molecule data for yeast kinesin-8 (Kip3) motor showed that its superprocessive movement is frequently interrupted by brief stick–slip motion. Here, a model is presented for the chemomechanical coupling of the kinesin-8 motor. On the basis of the model, the dynamics of Kip3 motor is studied analytically. The analytical results reproduce quantitatively the available single molecule data on velocity without including the slip and that with including the slip versus external load at saturating ATP as well as slipping velocity versus external load at saturating ADP and no ATP. Predicted results on load dependence of stepping ratio at saturating ATP and load dependence of velocity at non-saturating ATP are provided. Similarities and differences between dynamics of kinesin-8 and that of kinesin-1 are discussed.

A common ATP-dependent stepping model for kinesin-5 and kinesin-1: Mechanism of bi-directionality of kinesin-5

Kinesin-5 and kinesin-1 proteins are two families of kinesin superfamily molecular motors that can move processively on microtubules powered by ATP hydrolysis. Kinesin-1 is a unidirectional motor. By contrast, some yeast kinesin-5 motors are bidirectional and the directionality can be switched by changing the experimental conditions. Here, on the basis of a common chemomechanical coupling model, the dynamics of kinesin-1 and in particular the dynamics of kinesin-5 is studied theoretically, explaining the available experimental data. For example, the experimental data about different movement directions under different experimental conditions for kinesin-5 are explained well. The origin of why kinesin-1 can only make unidirectional movement and kinesin-5 can make bidirectional movements is revealed. The origin of mutations or deletions of several structural elements affecting the directionality of kinesin-5 is revealed. Moreover, some predicted results for kinesin-5 are provided.

A common chemomechanical coupling model for orphan and conventional kinesin molecular motors

Orphan and conventional kinesin dimers represent two families of the kinesin superfamily molecular motors. Conventional kinesin, having a 14-residue neck linker (NL) in each head, can step processively on microtubule (MT), with an ATP hydrolysis being coupled with a mechanical stepping under no load. Orphan kinesin phragmoplast-associated kinesin-related protein 2 (PAKRP2) dimer, despite having a NL of 32 residues in each head, can also step processively on MT and exhibits tight chemomechanical coupling under no load. However, the dynamic properties of the wild type PAKRP2 and the mutant one with each NL truncated to 14 residues are very different from those of the wild type conventional kinesin and the mutant one with each NL being replaced by the 32-residue NL from PAKRP2. Here, based on a common chemomechanical coupling model we study computationally the dynamics of the two families of the kinesin dimers, with the simulated results explaining quantitatively the available experimental data. The large differences in the dynamics between the two families of kinesin dimers arise mainly from different rate constants of NL docking and ATPase activity and different weak affinities of the head in ADP state for MT. The studies indicate that both the orphan kinesin PAKRP2 and conventional kinesin use the same mechanism for processive motility.

Theoretical Analysis of Dynamics of Kinesin Molecular Motors

Kinesin is a typical molecular motor that can step processively on microtubules powered by hydrolysis of adenosine triphosphate (ATP) molecules, playing a critical role in intracellular transports. Its dynamical properties such as its velocity, stepping ratio, run length, dissociation rate, etc. as well as the load dependencies of these quantities have been well documented through single-molecule experimental methods. In particular, the run length shows a dramatic asymmetry with respect to the direction of the load, and the dissociation rate exhibits a slip-catch-slip bond behavior under the backward load. Here, an analytic theory was provided for the dynamics of kinesin motors under both forward and backward loads, explaining consistently and quantitatively the diverse available experimental results.

Non-tight and tight chemomechanical couplings of biomolecular motors under hindering loads

Biomolecular motors make use of free energy released from chemical reaction (typically ATP hydrolysis) to perform mechanical motion or work. An important issue is whether a molecular motor exhibits tight or non-tight chemomechanical (CM) coupling. The tight CM coupling refers to that each ATPase activity is coupled with a mechanical step, while the non-tight CM coupling refers to that an ATPase activity is not necessarily coupled with a mechanical step. Here, we take kinesin, monomeric DNA helicase, ring-shaped hexameric DNA helicase and ribosome as examples to study this issue. Our studies indicate that some motors such as kinesin, monomeric helicase and ribosome exhibit non-tight CM coupling under hindering forces, while others such as the ring-shaped hexameric helicase exhibit tight or nearly tight CM coupling under any force. For the former, the reduction of the velocity caused by the hindering force arises mainly from the reduction of the CM coupling efficiency, while the ATPase rate is independent or nearly independent of the force. For the latter, the reduction of the velocity caused by the hindering force arises mainly from the reduction of the ATPase rate, while the CM coupling efficiency is independent or nearly independent of the force.

Run length distribution of dimerized kinesin-3 molecular motors: comparison with dimeric kinesin-1

Kinesin-3 and kinesin-1 molecular motors are two families of the kinesin superfamily. It has been experimentally revealed that in monomeric state kinesin-3 is inactive in motility and cargo-mediated dimerization results in superprocessive motion, with an average run length being more than 10-fold longer than that of kinesin-1. In contrast to kinesin-1 showing normally single-exponential distribution of run lengths, dimerized kinesin-3 shows puzzlingly Gaussian distribution of run lengths. Here, based on our proposed model, we studied computationally the dynamics of kinesin-3 and compared with that of kinesin-1, explaining quantitatively the available experimental data and revealing the origin of superprocessivity and Gaussian run length distribution of kinesin-3. Moreover, predicted results are provided on ATP-concentration dependence of run length distribution and force dependence of mean run length and dissociation rate of kinesin-3.

All-atom molecular dynamics simulations reveal how kinesin transits from one-head-bound to two-heads-bound state

Kinesin dimer walks processively along a microtubule (MT) protofilament in a hand-over-hand manner, transiting alternately between one-head-bound (1HB) and two-heads-bound (2HB) states. In 1HB state, one head bound by adenosine diphosphate (ADP) is detached from MT and the other head is bound to MT. Here, using all-atom molecular dynamics simulations we determined the position and orientation of the detached ADP-head relative to the MT-bound head in 1HB state. We showed that in 1HB state when the MT-bound head is in ADP or nucleotide-free state, with its neck linker being undocked, the detached ADP-head and the MT-bound head have the high binding energy, and after adenosine triphosphate (ATP) binds to the MT-bound head, with its neck linker being docked, the binding energy between the two heads is reduced greatly. These results reveal how the kinesin dimer retains 1HB state before ATP binding and how the dimer transits from 1HB to 2HB state after ATP binding. Key residues involved in the head-head interaction in 1HB state were identified.

Epistatic Analysis of the Contribution of Rabs and Kifs to CATCHR Family Dependent Golgi Organization

Multisubunit members of the CATCHR family: COG and NRZ complexes, mediate intra-Golgi and Golgi to ER vesicle tethering, respectively. We systematically addressed the genetic and functional interrelationships between Rabs, Kifs, and the retrograde CATCHR family proteins: COG3 and ZW10, which are necessary to maintain the organization of the Golgi complex. We scored the ability of siRNAs targeting 19 Golgi-associated Rab proteins and all 44 human Kifs, microtubule-dependent motor proteins, to suppress CATCHR-dependent Golgi fragmentation in an epistatic fluorescent microscopy-based assay. We found that co-depletion of Rab6A, Rab6A’, Rab27A, Rab39A and two minus-end Kifs, namely KIFC3 and KIF25, suppressed both COG3- and ZW10-depletion-induced Golgi fragmentation. ZW10-dependent Golgi fragmentation was suppressed selectively by a separate set of Rabs: Rab11A, Rab33B and the little characterized Rab29. 10 Kifs were identified as hits in ZW10-depletion-induced Golgi fragmentation, and, in contrast to the double suppressive Kifs, these were predominantly plus-end motors. No Rabs or Kifs selectively suppressed COG3-depletion-induced Golgi fragmentation. Protein-protein interaction network analysis indicated putative direct and indirect links between suppressive Rabs and tether function. Validation of the suppressive hits by EM confirmed a restored organization of the Golgi cisternal stack. Based on these outcomes, we propose a three-way competitive model of Golgi organization in which Rabs, Kifs and tethers modulate sequentially the balance between Golgi-derived vesicle formation, consumption, and off-Golgi transport.

Force dependence of unbinding rate of kinesin motor during its processive movement on microtubule

Kinesin is a biological molecular motor that can move continuously on microtubule until it unbinds. Here, we studied computationally the force dependence of the unbinding rate of the motor. Our results showed that while the unbinding rate under the forward load has the expected characteristic of "slip bond", with the unbinding rate increasing monotonically with the increase of the forward load, the unbinding rate under the backward load shows counterintuitive characteristic of "slip-catch-slip bond": as the backward load increases, the unbinding rate increases exponentially firstly, then drops rapidly and then increases again. Our calculated data are in agreement with the available single-molecule data from different research groups. The mechanism of the slip-catch-slip bond was revealed.