Basic Principles Underlying Cellular Processes
Daniel M. Zuckerman

Molecular Motors

$ \newcommand{\avg}[1]{\langle #1 \rangle} \newcommand{\cc}[1]{[\mathrm{#1}]^{\mathrm{cell}}} \newcommand{\cgdp}{\mathrm{C \! \cdot \! GDP}} \newcommand{\cgtp}{\mathrm{C \! \cdot \! GTP}} \newcommand{\comb}[1]{{#1}^{\mathrm{comb}}} \newcommand{\conc}[1]{[\mathrm{#1}]} \newcommand{\conceq}[1]{[\mathrm{#1}]^{\mathrm{eq}}} \newcommand{\concss}[1]{[\mathrm{#1}]^{\mathrm{ss}}} \newcommand{\conctot}[1]{[\mathrm{#1}]_{\mathrm{tot}}} \newcommand{\cu}{\conc{U}} \newcommand{\dee}{\partial} \newcommand{\dgbind}{\Delta G_0^{\mathrm{bind}}} \newcommand{\dgdp}{\mathrm{D \! \cdot \! GDP}} \newcommand{\dgtp}{\mathrm{D \! \cdot \! GTP}} \newcommand{\dmu}{\Delta \mu} \newcommand{\dphi}{\Delta \Phi} \newcommand{\dplus}[1]{\mbox{#1}^{++}} \newcommand{\eq}[1]{{#1}^{\mathrm{eq}}} \newcommand{\fidl}{F^{\mathrm{idl}}} \newcommand{\idl}[1]{{#1}^{\mathrm{idl}}} \newcommand{\inn}[1]{{#1}_{\mathrm{in}}} \newcommand{\ka}{k_a} \newcommand{\kcat}{k_{\mathrm{cat}}} \newcommand{\kf}{k_f} \newcommand{\kfc}{k_{fc}} \newcommand{\kftot}{k_f^{\mathrm{tot}}} \newcommand{\kd}{K_{\mathrm{d}}} \newcommand{\kdt}{k_{\mathrm{dt}}} \newcommand{\kdtsol}{k^{\mathrm{sol}}_{\mathrm{dt}}} \newcommand{\kgtp}{K_{\mathrm{GTP}}} \newcommand{\kij}{k_{ij}} \newcommand{\kji}{k_{ji}} \newcommand{\kkeq}{K^{\mathrm{eq}}} \newcommand{\kmmon}{\kon^{\mathrm{ES}}} \newcommand{\kmmoff}{\koff^{\mathrm{ES}}} \newcommand{\kconf}{k_{\mathrm{conf}}} \newcommand{\konf}{k^{\mathrm{on}}_{\mathrm{F}}} \newcommand{\koff}{k_{\mathrm{off}}} \newcommand{\kofff}{k^{\mathrm{off}}_{\mathrm{F}}} \newcommand{\konu}{k^{\mathrm{on}}_{\mathrm{U}}} \newcommand{\koffu}{k^{\mathrm{off}}_{\mathrm{U}}} \newcommand{\kon}{k_{\mathrm{on}}} \newcommand{\kr}{k_r} \newcommand{\ks}{k_s} \newcommand{\ku}{k_u} \newcommand{\kuc}{k_{uc}} \newcommand{\kutot}{k_u^{\mathrm{tot}}} \newcommand{\ktd}{k_{\mathrm{td}}} \newcommand{\ktdsol}{k^{\mathrm{sol}}_{\mathrm{td}}} \newcommand{\minus}[1]{\mbox{#1}^{-}} \newcommand{\na}{N_A} \newcommand{\nai}{N_A^i} \newcommand{\nao}{N_A^o} \newcommand{\nb}{N_B} \newcommand{\nbi}{N_B^i} \newcommand{\nbo}{N_B^o} \newcommand{\nc}{N_{C}} \newcommand{\nl}{N_L} \newcommand{\nltot}{N_L^{\mathrm{tot}}} \newcommand{\nr}{N_R} \newcommand{\nrl}{N_{RL}} \newcommand{\nrtot}{N_R^{\mathrm{tot}}} \newcommand{\out}[1]{{#1}_{\mathrm{out}}} \newcommand{\plus}[1]{\mbox{#1}^{+}} \newcommand{\rall}{\mathbf{r}^N} \newcommand{\rn}[1]{\mathrm{r}^N_{#1}} \newcommand{\rdotc}{R \!\! \cdot \! C} \newcommand{\rstarc}{R^* \! \! \cdot \! C} \newcommand{\rstard}{R^* \! \! \cdot \! D} \newcommand{\rstarx}{R^* \! \! \cdot \! X} \newcommand{\ss}{\mathrm{SS}} \newcommand{\totsub}[1]{{#1}_{\mathrm{tot}}} \newcommand{\totsup}[1]{{#1}^{\mathrm{tot}}} \newcommand{\ztot}{Z^{\mathrm{tot}}} % Rate notation: o = 1; w = two; r = three; f = four \newcommand{\aow}{\alpha_{f}} \newcommand{\awo}{\alpha_{u}} \newcommand{\kow}{\kf} % {\kf(12)} \newcommand{\kwo}{\ku} % {\ku(21)} \newcommand{\kor}{\conc{C} \, \konu} % \konu(13)} \newcommand{\kwf}{\conc{C} \, \konf} % \konf(24)} \newcommand{\kro}{\koffu} % {\koffu(31)} \newcommand{\kfw}{\kofff} % {\kofff(42)} \newcommand{\krf}{\kfc} % {\kfc(34)} \newcommand{\kfr}{\kuc} % {\kuc(43)} \newcommand{\denom}{ \krf \, \kfw + \kro \, \kfw + \kro \, \kfr } $

Kinesin: A Prototypical Molecular Motor

Molecular motors generate motion by coupling conformational changes to ATP hydrolysis. In doing so, motors use the free energy stored in ATP (and ATP itself stores free energy because it is maintained out of equilibrium as an activated carrier).

Kinesin is a two-headed motor, like many others (dynein, and most myosins). Each motor head is capable of hydrolyzing ATP, which results in a conformational change. The heads are also coupled to one another via "neck linker" domains. A simplified version of kinesin's cycle is shown below.

Kinesin Cycle

In the figure, the green triangle represents a fixed location along the microtubule, which is shown as a simple horizontal line. The colored ovals are the two heads, and the colored lines are the neck linkers. The chemical state is indicated in the oval, with a subscript noting whether a head is forward (f) or in the rear (r). The top-most transition, which is simply unbinding of the rear head, is somewhat confusing in the four-state diagram because the identity of the forward/rear heads has switched: note the green triangle and see figure below.

Kinesin moves along microtubulue filaments in the cell and plays a key role in mitosis. The steps are 8nm in length and discrete based on the polymeric structure of the microtubule. Hence we can label the positions as (..., i-2, i-1, i, i+1, i+2, ...). The cycle of chemical and mechanical steps is probably best visualized in the context of a microtubule. The pink line below shows the order of steps in purely forward motion.

Kinesin Microbubule

A simple quantitative model for kinesin motion

Kinesin Cycle Numbered

To turn our original cartoon into a concrete model, we simply need to assign rate constants. We will assign rates derived by Liepelt and Lipowski ("LL" - see References) on the basis of experiments. Notice that the state numbers (1, 2, 5, 6) are not consecutive: this is simply for consistency with LL; our states are a subset of theirs, although this does not lead to significant deviations from the original model's behavior for the cases examined here.

The first thing to do with a kinetic model is write out the basic mass-action equations. For brevity, we'll note the first two here,


where $\conc{X}$ is the concnetration of state X and $k_{ij}$ is the rate constant from state $i$ to $j$. Remember, you can always consult the Notation page.

Equilibrium analysis of the cycle

It is necessary to perform an equilibrium analysis of the cycle for two reasons: (i) to ensure the rates permit detailed balance to hold in equilibrium, and (ii) to ensure that the rates are consistent with the thermodynamic balance holding between ATP and its hyrolysis product, ADP.

We apply the condition of detailed balance (equal and opposite flows) to all legs of the cycle. For the 1-2 and 2-5 steps, we have


To be clear, we are not assuming our system is always in equilibrium. Rather, we are establishing that our system should be capable of reaching equilibrium consisent with physical principles.

Although each of the (3) - (6) indidvidually represents a detailed-balance condition for a particular process, together they yield a single constraint on the rate constants. To see this, solve for $\conceq{2}$ in (4) and substitute the result in (3); then continue this process and successively eliminate all the state concentrations. You will be left with a relation between the rate constants and the concentrations of ATP, ADP, and Pi:


where we have noted that the ratio of concentrations on the right is nothing other than the equilibrium constant for ATP, its effective "dissociation" constant. The relation (7) implies that once 7 of the rate constants are known, the 8th is determined by thermodynnamic considerations.

Following LL, we set $K_d = 4.9 \times 10^5$ M, which is a large value reflecting that, in solution, ATP hydrolysis is strongly favored over synthesis - ADP will predominate over ATP unless ATP is synthsesized. In a cell, of course, ATP is constantly being synthesized to maintain a non-equilibrium condition in which ATP stores free energy available for use by many processes.

A Subtlety. As shown in the diagram with the "cycle" superimposed on a microtubule, the states really are specific to the position (i, i+1, ...). Nevertheless, the equilibrium analysis we just performed ignoring the position is correct. To see this, consider an ensemble of single kinesin molecules, each walking on a separate circular microtubule. In such a case, it is clear that the average population of each chemical/configurational state is independent of the positional state i - e.g., the population of state 6 is the same at every position. Hence, the equilibrium properties can be ascertained by assuming we have a simple cycle: the correct populations are obtained, aside from a trivial factor accounting for the number of positions, which is the same for each state.

The equilibrium constant and directionality.

When the concentrations of ATP, ADP and Pi are at equilibrium values, which satisfy (7), the system will not be driven. That is, on average, the number of forward and reverse steps per unit time will be the same. In typical cellular conditions, however, the ratio of concentrations will greatly favor ATP - and hence forward steps will predominate. That is, we expect


See the discussion of simulation data, below, for some examples. Although reverse stepping is possible in principle, and in the present model, note that realistic concentrations in fact make reversals highly unlikely - see exercises.

Kinesin as Catalyst

An alternative point of view is that kinesin may be seen merely as a catalyst for ATP hydrolysis (and synthesis, under suitable conditions), as sketched above. In this picture, the stepping of kinesin can be considered an unnecessary consequence of catalyzed event: see also the "thought experiment" proposed for ATP in its role as an activated carrier. The basic view of ATP is correct although to be fully accurate, one would have to account for work done against viscous forces in typical steps - see below.

The issue of work. In the present model, for simplicity, there is no "load" or force on our motor - and so mechanical forward/backward steps perform no work. This is unphysical because a viscous/frictional force will always oppose the direction of motion. Nevertheless, the zero-load case is of interest for building understanding of motor systems, and compared to experimental situations, the viscous force may be considerably smaller than forces which are applied on such motors.

Model parameters

Specific rate constants are given in the table below, which are derived from the LL paper.

ATP binding$k_{12}$$2 \times 10^6$ / (M s)
ATP unbinding$k_{21}$$100$ / s
Forward step$k_{25}$$3 \times 10^5$ / s
Reverse step$k_{52}$$0.24$ / s
ADP unbinding$k_{56}$$100$ / s
ADP binding$k_{65}$$2 \times 10^4$ / (M s)
Hydrolysis and Pi unbinding$k_{61}$$0.74$ / s
Pi binding and synthesis$k_{16}$$2 \times 10^4$ / (M s)

Note that $k_{61}$ was calculated from other rate constants and $K_d$ using (7).

Simulating the model

We can simulate the simple kinetic model using BioNetGen, a rule-based platform for kinetic modeling. If we assign each mechanical step (the 2-to-5 transition) the experimentally determined value of 8 nm, we can graph the simulation results for different concentrations of ATP, ADP and Pi.

Kinesin Forward

Not surprisingly, higher ATP concentration leads to faster motion. But there are limits on how fast the motor can move - see the exercises for hints. In the graph above, the ATP concentration increased 100-fold, but the speed of the motor clearly did not.

Although it is safe to say that a given kinesin motor (there are many kinds) evolved to move in a single direction, it is of physical interest to push our model to extremes. In particular, given that all molecular processes are reversible, can we get our motor to run in reverse? You should recognize that forward steps are expected whenever the nucleotide concentrations are shifted toward larger ATP - i.e., when the actual concentrations differ from their equilibrium values and excess ATP is present as in (8). By the same argument, reverse steps are expected when there is an excess of ADP and/or Pi - when the opposite of (8) holds.

Kinesin Backwards

Note the timescale and nucleotide concentrations compared to previous data with forward stepping. Such concentrations are unphysical as can be seen from the exercises.

These simulations were performed using BioNetGen, a rule-based platform for kinetic modeling. The source code for the model (a .bngl file) can be downloaded by right-clicking here.

Making a more complete model

Our discussion has avoided two issues important to detailed consideration of motor proteins: the effect of a significant opposing force which may be present if a cargo is attached to the motor, and a more complete view of the potential state space. Qualitatively, it is clear that an opposing force or "load" will decrease the speed of the motor at given nucleotide concentrations; for further details, the reader is referred to the LL paper and references therein.

Kinesin Cycle Big Space

Using a conceptual framework similar to what was presented in our advanced discussion of cycle logic, we can also consider a more complete state space for the kinesin motor. In the figure above, D, T, and $\phi$ represent possible states of an indivdual motor head: ADP-bound, ATP-bound, and empty, respectively. A pair of symbols represents the rear/forward combination - e.g., "DT" means the rear head is ADP-bound and the forward head is ATP-bound - using the convention that right is forward.

The space shown here may seem large and complicated, but really we have omitted a number of possibilities: that both heads could become detached, the position on the filament, additional configurational states, and states with both ADP and Pi bound.

  • General references
    • B. Alberts et al., "Molecular Biology of the Cell," Garland Science (many editions available).
    • T.L. Hill, "Free Energy Transduction and Biochemical Cycle Kinetics," (Dover, 2005).
    • J. Howard, "Mechanics of Motor Proteins and the Cytoskeleton," (Sinauer, 2001).
  • The specific kinesin model used is based on the following paper:
  • Experimental studies of kinesin's chemomechanical cycle - only a tiny selection of a vast literature:
    • N. Carter and R. Cross, "Mechanics of the kinesin step," Nature 435, 308-312 (2005).
    • N. Guydosh and S. Block, "Backsteps induced by nucleotide analogs suggest the front head of kinesin is gated by strain," PNAS 103 (21) 8054-8059 (2006).
    • D. Hackney, "The tethered motor domain of a kinesin-microtubule complex catalyzes reversible synthesis of bound ATP," PNAS 102:18338-18343 (2005).
    • W. Schief, R. Clark, A. Crevenna, and J. Howard, "Inhibition of kinesin motility by ADP and phosphate supports a hand-over-hand mechanism," PNAS 101:1183-1188 (2004).
  1. Derive (7).
  2. In a mass-action model, the overall rate of ATP binding and hydrolysis will increase in proportion to the the ATP concentration without limit. In reality, a maximum saturating value is expected for the binding step: Why? [Hint - review the mass-action page.]
  3. Even if the world truly obeyed mass-action principles, there would still be a limit on the speed with which the motor could move given the present model. Discuss this issue by examining rate constants and limiting cases.
  4. Use a geometric argument based on the size of an ADP molecule and the molecular spacing at $\conc{ADP} = 10^3$ M to show that concentration is unphysical. Explain why the mass-action picture breaks down.