# Strength and stability of active ligand-receptor bonds: a microtubule attached to a wall by molecular motor tethers

## Abstract

We develop a stochastic kinetic model of a pre-formed attachment of a mictrotuble (MT) with a cell cortex, in which the MT is tethered to the cell by a group of active motor proteins. Such an attachment is a particularly unique case of ligand-receptor bonds: The MT ligand changes its length (and thus binding sites) with time by polymerization-depolymerization kinetics, while multiple motor receptors tend to walk actively along the MT length. These processes, combined with force-mediated unbinding of the motors, result in an elaborate behavior of the MT connection to the cell cortex. We present results for the strength and lifetime of the system through the well-established force-clamp and force-ramp protocols when external tension is applied to the MT. The simulation results reveal that the MT-cell attachment behaves as a catch-bond or slip-bond depending on system parameters. We provide analytical approximations of the lifetime and discuss implications of our results on in-vitro experiments.

Correspondence: debch@iitk.ac.in , rwfridd@sandia.gov

## Introduction

Chromosome segregation is the most important process during the mitosis phase of cell cycle pollardbook (). In eukaryotic cells, sister chromatids, that result from chromosome replication, are segregated by a complex multi-component machine called mitotic spindle mcintosh12 (); bouck08 (); helmke13 (). Microtubule (MT) lawson13 (), a stiff tubular filament, whose typical diameter is about 25 nm, forms a major component of the scaffolding of the spindle. A MT is a polar filament in the sense that its two ends are dissimilar; the plus end is more dynamic than the minus end. During the morphogenesis of the spindle petry16 (); kapoor17 (), MTs form transient molecular joints (non-covant bonds) with specific partners. The attachments of the MTs with the chromosomes, mediated by a proteineous complex called kinetochore cheeseman08 (); dewulf09 (), has been under intense investigation in recent years. The forces exerted by these MTs on the kinetochores eventually pull the two sister chromatids apart in the late stages of mitosis thereby driving their journeys to the opposite poles of the spindle asbury17 (); scholey16 (). In contrast, another set of MTs, called astral MT, diverge from the spindle poles and their distal ends (the so-called plus ends) form contacts with the cell cortex. The kinetochore-MT attachment has been under intense investigation in recent years for the obviously important role it plays in timely and accurate segregation of the chromosomes. However, in this paper we focus exclusively on the MT-cortex attachment.

Understanding the physics of the attachment formed by a single MT with the cortex is the first step in ultimately understanding how forces generated by all such attachments collectively determine the position and orientation of the spindle laan12 (); kotak13 (); lu13 (); mcnally13 (); pietro16 (). The studies reported in this paper are important also from the perspective of research on MTs and cytoskeletal motor proteins that use MT as the ‘track’ chowdhury13a (); kolomeisky15 (), particularly those which play crucial force-coupling roles by residing at the plus end tamura12 (); wade09 (); akhmanova15 (). Moreover, in spite of the simplicity of the system of our study where only a single MT is attached to the cortex by active tethers, our analysis reveals the cooperative effects of multiple motors that give rise to the emergent collective properties of the attachment. Such collective phenomena are of general interest in several branches of physics and biology mclaughlin16 ().

The main aim of this paper is to develop a minimal mathematical model of an attachment formed by a single astral MT with the cell cortex by capturing the essential roles of only the key components identified experimentally till now. The plus end of an astral MT is linked with the cortex of a living eukaryotic cell by dynein molecules tuncay16 (). There are equispaced dynein-binding sites on the surface of a MT where the head of a dynein can bind specifically. Dynein is a motor protein whose natural tendency to walk towards the minus end of the MT track is powered by input chemical energy extracted from ATP hydrolysis chowdhury13a (); kolomeisky15 (). If the tail of a dynein is anchored on the cortex, its attempt to step towards the minus end of the MT results in an effective pull on the MT directed towards the cortex. Thus, the dyneins function as active tethers linking the MT with the cortex. The model developed in this paper captures only some of the key features of the real MT-cortex attachments. It is worth pointing out that the MT-wall model developed here differs fundamentally from another MT-wall model reported earlier by Sharma, Shtylla and Chowdhury sharma14 () (onwards referred to as the SSC model). Unlike the active (i,e., energy consuming) linkers, which here mimic dynein motors, the dominant tethers in the SSC model are passive. Moreover, the rigid wall in the SSC model represents a kinetochore cheeseman08 () whereas the wall in the model developed in this paper represents the cell cortex.

An attachment of the type studied in this paper can be viewed as an analog of a ‘ligand-receptor bond’ bongrand99 (); karplus10 () where the MT filament is the analog of a ligand while the specific binding partners that link its plus end with the cortex are the analogs of receptors. However, unlike common ligands, a MT exhibits a unique polymerization-depolymerization kinetics desai97 () and the corresponding receptor proteins are ‘active’ in the sense that these consume chemical fuel for their mechanical function. Drawing analogy with ligand-receptor bonds evans07 (), we analyze the mathematical model to make quantitative predictions on the strength and stability of the attachment.

We consider a pre-formed attachment to investigate its strength and stability using a kinetic model that mimics the protocols of dynamic force spectroscopy bizzarri12 (); arya16 (). Two distinct protocols are routinely used in force spectroscopy for measuring the strength and stability of ligand-receptor bonds bongrand99 (). In the force-clamp protocol, a time-independent load tension is applied against the bond; the time duration after which the bond is just broken is the life time of the bond. Since the underlying physical process is dominated by thermal fluctuations different values of life time are observed upon repetition of the experiment. Therefore, the stability of the bond is characterized by the life time distribution (LTD). On the other hand, in the force-ramp protocol the magnitude of the time-dependent load tension is ramped up at a pre-decided rate till the bond just gets ruptured; the rupture force distribution (RFD) characterizes the strength of the bond. For a slip-bond the mean life time (MLT) decreases monotonically with the increasing magnitude of the tension in the force-clamp experiment. In contrast, a non-monotonic variation (an initial increase followed by decrease) of the MLT with increasing load tension in the force-clamp experiments is a characteristic feature of catch-bonds thomas08a (); thomas08b ().

## Model

Microtubules are cylindrical hollow tubes with, approximately, 25 nm diameter. Hetero-dimers, formed by globular proteins called and tubulins, assemble sequentially to form a protofilament. The length of each - dimer is about 8 nm. Normally 13 such protofilaments, arranged parallel to each other, form a microtubule. However, there is a small offset of about 0.92 nm between the dimers of the neighboring protofilaments. In the single protofilament model anderson13 () each MT is viewed as a single protofilament that grows helically with an effective dimer size nm. Thus, we model the MT as a strictly one-dimensional stiff filament. The plus end of the MT is oriented along the +X-direction of the one-dimensional coordinate system chosen for the model.

In our model the cortex is represented by a rigid wall oriented perpendicular to the X-axis. Since a dynein motor, fueled by ATP, has a natural tendency to walk towards the minus end of the MT, the molecular motors in our model are also assumed to be minus-end directed, while consuming input energy, in the absence of external load tension. We assume that the tail of each dynein motor is permanently anchored on that face of the rigid wall which faces the plus end of the MT. The heads of a motor, however, can either attach or detach from the MT. Each dynein motor has two heads that are, effectively, connected at the midpoint by a hinge. We denote the position of a motor by the location of its midpoint; the midpoint is assumed to be linked to the point of anchoring on the wall by an elastic element that is assumed to be a Hookean spring.

In the one-dimensional coordinate system the origin is fixed on the wall. With respect to this origin denotes the position of the midpoint of th molecular motor at time while denotes the corresponding position the MT tip.

is the total number of dynein motors that can simultaneously attach to the MT whereas denotes the number of motors actually attached at any arbitrary instant of time (i.e., . For any given unbound motor, denotes the rate of binding of its head to the MT. Therefore, the rate at which any unbound motor binds to the MT is muller08 ()

(1) |

Because of being anchored on the fixed rigid wall, the motors bound to the MT cannot walk freely along the MT track. Note that even if a motor remains bound to a fixed site on the MT and steps neither forward nor backward, the spring force acting on it changes in response to the externally applied force . Moreover, forward or backward stepping of a motor on the MT extends or compresses the spring thereby altering the spring force.

In our model represents the motor location (). Applying a tensile force, , to the MT minus-end will move all motors by a distance (Fig. 2). The displacement due to external force and due to walking are contained within , and thus the force experienced by the motor is given by the corresponding force of the elastic linkage at that extension,

(2) |

Following Kramers (or Bell) theory, bell78 (); kramers40 (), we assume that the effective rate of unbinding one of the motors attached to the MT at a given instant of time is approximately muller08 ()

(3) |

where denotes the rate of unbinding of the head of a single motor in the absence of load force. The characteristic ‘detachment force’ can be expressed as where is the distance from the energy minimum to maximum for the motor/MT interaction potential.

Similarly, the effective rate of forward stepping of a motor is given by

(4) |

Based on experimental observations, we assume that a motor can also step towards the positive end of MT mallik04 (); vale06 (); toba06 (); thus we model stepping in the ‘reverse’ direction by

(5) |

where the constant parameter () is the fraction of the path over which work is done and the characteristic force can be expressed as where is the length of a single subunit of MT. The rate of stepping of a motor towards the plus end of the MT in the absence of load force is very small because the natural direction of these motors is the minus end of MT. The form (5) captures the intuitive expectation that would increase with increasing force .

The probability that (the midpoint of) a motor is located at and MT tip is at , while the total number of motors are bound to the MT simultaneously at that instant of time, is given by . Velocity of whole MT body is given by . Note that is a continuous variable whereas can take only non-negative integer values. Let be the conditional probability that, given a MT-bound motor located at site , there is another MT-bound motor at site on the MT . Then is the conditional probability that, given a motor at site , the site is empty. Let be the probability that site is not occupied by any motor, irrespective of the state of occupation of any other site.

Under mean field approximation (MFA), the equations governing the time evolution of is given by

(6) |

where

(7) |

and is the viscous drag coefficient.

The rates of polymerization and de-polymerization of a MT tip are given by and , respectively. The rate of depolymerization of MT is suppressed by externally applied tension franck07 (). We assume that the MT-bound minus-end directed motors at the tip (plus-end) of the MT prevents MT protofilaments from curling outwards, thereby slowing down or speeding up depolymerization rate depending upon the position of the MT tip laan12 ():

(8) |

where is the characteristic load force at which the MT depolymerization rate is an exponentially small fraction of . The Kronecker delta function ensures that the external force affects the depolymerization rate only if the motor is bound to the tip of the MT.

The force balance equation is given by

(9) |

where is a Gaussian white noise and , the length of each subunit of MT is also the spacing between the successive motor-binding sites on the MT.

## Simulation method

Parameter | Values |
---|---|

Spacing between binding sites on MT joglekar02 (); hill85 (); shtylla11 () | nm |

Rate of polymerization of MT joglekar02 (); hill85 (); shtylla11 (); waters96 () | 30 s |

rate of Load-free depolymerization of MT joglekar02 (); hill85 (); shtylla11 (); waters96 () | 350 s |

Rate of Binding of motor to MT | 3 s |

rate of Unbinding of motor from MT | 3 s |

rate of Forward stepping of motor | 6 s |

rate of Backward stepping of motor | 0.1 s |

Characteristic depolymerization force | 0.1 pN |

Characteristic detachment force | 1 pN |

Characterstic spring force of motor | 1 pN |

Rest length of elastic linkage | 6.2 nm |

Linkage spring constant | 1 pN/nm |

Stepping parameter | 0.5 |

Effective drag coefficient joglekar02 (); hill85 (); shtylla11 (); marshall01 () | 6 pNsm |

The simulations based on our theoretical model was carried out using the Gillespie algorithm gillespie07 (). In each time step , eight types of events are possible, namely, binding/ unbinding, forward/backward hopping of any motor, polymerization and depolymerization of MT tip and forward/backward movement of the whole MT body.

Initially all the motors are attached to randomly selected positions on the MT, also the MT tip is placed adjacent to the wall. Rate constants are then determined based on each motor position. A motor can unbind from its occupied site with the unbinding rate given by eq.3. Similarly, a motor can jump forward or backward with the transition rate eq.4 and eq.5 respectively, provided the neighboring site is empty. Finally a depolymerization event, goverened by eq.8, can carry away a motor if it is concurrently located at the MT tip.

At equilibrium, unbound motors are spatially distributed in 1D by the Boltzmann-weighted energy of their tethering springs. (see Fig: 3)

(10) |

with 1 pN/nm the standard deviation of unbound motor fluctuations is only nm. That is, rarely would the spring naturally stretch more than nm. A suitable location for binding is drawn from the cumulative distribution by inverse transform sampling, and checked if the chosen site is empty. In this way, the position for binding is selected.

By the equation of motion (9), the polymerization and depolymerization rates control by changing the length of the MT. But external forces influence movement of the whole MT by the resultant force acting on it, i.e. . To treat movement of the whole MT within our discretized system, we define the corresponding rate constant by . The sign of the expression will decide the movement of MT in forward (positive) or backward (negative) direction.

Since the portion of available MT length is quickly depolymerized, and the instance of zero bound motors will allow the MT drift away from the wall, we neglect rebinding from the 0 state. Thus the time evolution of the motor-MT attachment is monitored until, for the first time, all the motors are detached from the MT. This first passage time is identified as the lifetime of the attachment for both the force clamp and force ramp conditions. We have generated trajectories of up to time steps which were then averaged to arrive at the results of interest. The common parameter values used in the simulation are listed in table 1.

## Results and Discussion

## Force clamp condition

Figure 4(a) shows that the mean lifetime initially increases then decreases with external tension ; such non-monotonic variation of the lifetime represents catch-bond-like behavior of cortical dynein-MT attachment. This can be explained by the two pathways which lead to removing a motor from the MT: ) By breaking the motor/MT bond, or ) by depolymerization of a tip subunit when a motor is bound to it.

We now to turn to finding an analytically tractable approximations of the lifetime. We coarse-grain the relevant microscopic states of the system, , into two macrostates: the completely attached () and completely detached () configurations. The quantity of interest is therefore the effective transition rate (lifetime) representing a one-step process .

For small forces, depolymerization is faster than any other process in the system and is therefore the primary cause of motor detachment, while spontaneous unbinding adds a small contribution to the loss of motors. Neglecting rebinding and stepping dynamics for simplicity, we consider the zero-force rates of detaching all motors by spontaneous unbinding, , or by depolymerization of motors close to the tip, . In either case, decay of the bound state can be viewed as a sequence of independent Poisson processes Erdmann04 (); Williams03 (),

(11) | |||

(12) |

where is the harmonic number. We note that including the polymerization rate adds negligible correction to . The mean lifetime follows as the reciprocal of the sum of the two rates of leaving the completely attached macrostate,

(13) |

For the parameters in Table 1 with equation (13) evaluates to ms, which corresponds to the simulated lifetime at zero force of 27.6 ms (figure 4a).

At large forces, the detachment of motors is no longer governed by depolymerization because the load of a motor at the MT tip will stall depolymerization according to Eq. (8). Assuming the motors are close to one another on the MT, we approximate unbinding by the usual multiple-parallel bond treatment Seifert00 (); Erdmann04 (); Williams03 () , which for the rate of motor removal is dictated by the time to rupture the first bond,

(14) |

At 10 pN and equation (14) gives ms, which coincides with the tail of all three curves at 10 pN in figure 4a.

In between these limits the detachment of motors is governed by one pathway that is enhanced by force () and another that is suppressed by force (). As mentioned above we heuristically treat the intermediate force regime with effective rates of removing all motors from the fully attached state. Again neglecting the slow rebinding and motor stepping we approximate the mean detachment time as,

(15) |

where the parameters of the two terms are effective parameters for the unbinding and depolymerization processes respectively. The best fit to the 0.1 pN data in Fig. 4(a) corresponds to , and , . For fixed force, 1 pN, the first term of equation (15) depends most strongly on ,

(16) |

The inset of Fig:4(a) shows the simulated trend of lifetime with motor number. At large the lifetime becomes independent of A fit of Eq. (16) at 1 pN yields 7.4 s and 0.75 pN.

We also plot the survival probability , the probability that till time the MT-motor attachment survives, in the inset of Fig:4(b) the survival probability of the MT-motor attachment is plotted for different values of external tension . The attachment survives longer at intermediate forces (for example, at pN (blue square)) than at high and low forces. Each of the survival probabilities shown in Fig.4(b), for which has a fixed value, has been calculated using friddle12 (); arya16 (),

## Force ramp condition

Here we present results of our simulation under the common force spectroscopy protocol of measuring the rupture force of the attachment when increasing force with time , and repeating over a range of loading rates, . In Fig: 5(a) and (b) the probability distribution of rupture force and survival probability are shown for the same range of loading rates. It is important to note that at low loading rate the most probable rupture force is not zero, but instead is peaked around 0.2 pN (Fig. 5(a), 10 pNs). This indicates a near-equilibrium regime in which the wall-MT attachment is stabilized by the rebinding rate . As loading rate increases, the rupture transitions from a near-equilbrium to a kinetic regime as a second peak emerges at high rupture force and becomes prominent (Fig. 5(a), 50 pNs).

The fitting of simulation data has been done using a function friddle12 (); arya16 (),

(18) |

We compare the approximation for the lifetime in Eq.15 against the simulated data through calculating the survival probability and corresponding probability distribution as friddle12 (); arya16 (),

(19) |

Equation 19 coincides with the general form of the simulated data at high forces, but diverges at low forces. This can be expected since eqs. 18,19 assume a first-passage process and does not account for binding reversibility at low loading rate and low force. The mean rupture force with loading rate shows the familiar behavior of a force spectrum. At low loading rates the mean force plateaus at a very small equilibrium force, then transitions into a kinetic regime where the force goes approximately as the log of the loading rate (Fig:5(c)). In the irreversible approximation the mean rupture force, , is given by friddle12 (); arya16 (),

(20) |

We substitute from Eq.(15) into Eq.(19) to calculate mean rupture force from (20). We find good agreement at both low and intermediate loading rates. At large loading rates the data becomes non-linear with log-loading rate, growing to larger rupture forces than the prediction of Eq. (20). This is likely explained by considering that faster loading rates allow less time for the dissociation and depolymerization processes to remove bound motors, ultimately leading to rupture of more motor-MT bonds.

## Conclusions

Extending the earlier generalizations karplus10 () of the concept of a ligand, we have treated a microtubule (MT) as a ‘ligand’ that is tethered to a ‘receptor’ wall by a group of minus-end directed molecular motors chowdhury13a (). The tails of the motors are permanently anchored on the wall while their motor heads can bind to- and unbind from the MT. This model of MT-wall attachment captures only a few key ingredients of the MT-cortex attachments in eukaryotic cells, particularly those formed during chromosome segregation. This minimal model incorporates the polymerization and depolymerization kinetics of MT. But, for the sake of simplicity, it does not include the processes of ‘catastrophe’ and ‘rescue’ that are caused by the ’dynamic instability of MT filaments desai97 () although these can be captured in an extended version of this model ghanti16 (). We consider a pre-formed MT-wall attachment and carry out computer simulations to study statistical properties of its rupture under conditions that mimic the protocols of force-clamp and force-ramp experiments in-vitro bizzarri12 (); arya16 (). The simulation results that we report are interpreted in the light of the theory of single-molecule force spectrocopy, popularized by Bell bell78 () and some of its later generalizations arya16 ().

In-vitro experiments have been designed which remarkably resemble the conceptual model depicted in Fig.1, albeit without application of a controlled force. Laan et al. laan12 () used a microfabricated vertical barrier that mimics the cell cortex. Dynein motors were anchored on the barrier and captured the MT that grew from a centrosome which was fixed on a horizontal glass surface. A slightly different experimental set up was used by Hendricks et al. hendricks12 () in which a dynein coated bead was used to mimic the cell cortex. These authors demonstrated the stabilization of the MT within the broader context of the role of MT-cortex interaction in positioning of the mitotic spindle laan12 (); hendricks12 (). However the possibility of controlling the force on the MT by an Atomic Force Microscope or Optical Trap should motivate extending these experimental setups to test the results of our theoretical model in the near future.

## Acknowledgments

DC thanks Gaurav Arya for valuable comments on an earlier preliminary version of the work. This work has been supported by a J.C. Bose National Fellowship (DC) and by the “Prof. S. Sampath Chair” Professorship (DC). Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energyâs National Nuclear Security Administration under contract DE-NA-0003525.

### References

- Pollard, T.D., Earnshaw, W.C., Lippincott-Schwartz, J., and Johnson, G.T. 2016. Cell Biology (3rd ed.). Elsevier.
- McIntosh, J. R., Molodotsov, M. I., and Ataullakhanov, F. I., 2012. Quart. Rev. Biophys. 45: 147-207.
- Bouck, D. C., Joglekar, A. P., and Bloom, K. S., 2008. Annu. Rev. Genet. 42: 13.1-13.25 .
- Helmke, K.J., Heald, R. and Wilbur, J.D., 2013. Int. Rev. Cell Mol. Biol. 306: 83-125.
- Lawson, J. D. and R.E. Carazo Salas. 2013. Biochem. Soc. Trans. 41:1736–1744.
- Petry, S. 2016. Annu. Rev. Biochem. 85: 659-683.
- Kapoor, T.M. 2017. Biology, 6: 8.
- Cheeseman, I. M. and A. Desai. 2008. Nat. Rev. Mol. Cell. Biol. 9:33–46.
- DeWulf, P., Earnshaw, W.C. (Eds.) 2009. The Kinetochore: From Molecular Discoveries to Cancer Therapy. Springer.
- Asbury, C. L. 2017. Biology 6:15.
- Scholey, J.M., Civelekoglu-Scholey, G. and Brust-mascher, I. 2016. Biology 5:15.
- Laan, L., S. Roth, M. Dogterom. 2012. Cell Cycle. 11:3750–3757.
- Kotak, S. and Gönczy, P. 2013. Curr. Opin. Cell Biol. 25: 741-748.
- Lu, M.S. and Johnston, C.A. 2013. Development 140:1843-1856.
- McNally, F. J., 2013. J. Cell Biol., 200: 131-140.
- di Pietro, F., Echard A,. Morin, X. 2016. EMBO Rep., 17: 1106-1130.
- Chowdhury, D. 2013. Phys. Rep. 529:1–198.
- Kolomeisky, A. B. 2015. Motor proteins and Molecular Motors. CRC Press.
- Tamura, N. and Draviam, V.M. 2012. Open Biol. 2: 120132.
- Wade, R. H. 2009. Mol. Biotechnol. 43: 177-191.
- Akhmanova A, Steinmetz, M.O. 2015. Nat. Rev. Mol. Cell Biol., 16: 711-726.
- McLaughlin R. T., Diehl, M.R., Kolomeisky, A.B. 2016. Soft Matter 12:14-21.
- Tuncay, H. and K. Ebnet. 2016. Cell. Mol. Life Sci. 73:1195–1207.
- Sharma, A. K., B. Shtylla and D. Chowdhury. 2014. Phys. Biol. 11:1478–3975.
- Bongrand, P. 1999. Rep. Prog. Phys. 62:921–968.
- Karplus, M., 2010. J. Mol. Recognit. 23:102-104.
- Desai, A. and T.J. Mitchison. 1997. Annu. Rev. Cell Dev. Biol. 13:83–117.
- Evans, E. A., Calderwood, D.A. 2007. Science, 316: 1148-1153.
- Bizzarri, A. R. and S. Cannistraro (eds.). 2012. Dynamic Force Spectroscopy and Biomolecular Recognition. CRC Press.
- Arya, G. 2016. Molecular Simulation. 42:1102–1115.
- Thomas, W. E. 2008. Annu. Rev. Biomed. Eng. 10:39–57.
- Thomas, W. E., V. Vogel and E. Sokurenko. 2008. Annu. Rev. Biophys. 37:399–416.
- Muller, M. S. Klumpp, and R Lipowsky. 2008. PNAS. 105:4609–4614.
- Bell, G. I. 1978. Science. 200:618–627.
- Kramers, H. A. 1940. Physica. 7:284–304.
- Bowne-Anderson, H., M. Zanic, M. Kauer and J. Howard. 2013. Bioessays. 35:452–461.
- Mallik, R., B.C. Carter, S.A. Lex, S.J. King, and S.P. Gross. 2004. Nature. 427:649–652.
- Reck-Peterson, S. L., A. Yildiz, A.P. Carter, A. Gennerich, N. Zhang, and R.D. Vale. 2006. Cell. 126:335–348.
- Toba, S., T.M. Watanabe, L. Yamaguchi-Okimoto, Y.Y. Toyoshima, and H. Higuchi. 2006. Proc. Natl. Acad. Sci. USA. 103:5741–5745.
- Franck, A. D., A. F. Powers, D. R. Gestaut, T. Gonen, T. N. Davis, C. L. Asbury. 2007. Nat Cell Biol. 9:832–837.
- Redner, S. 2001. A Guide to First-Passage Processes. Cambridge University Press. Cambridge.
- Redner, S., R. Metzler and G. Oshanin. (eds) 2014. First-passage Phenomena and their Applications. World Scientific Publishing Co.
- Iyer-Biswas, S. and A. Zilman. 2016. Adv. Chem. Phys. 160:261–356.
- Polizzi, N., M.J. Therien, D.N. Beratan. 2016. Israel J. Chem. 56:816–824.
- Wang, H., C. Peskin and T. Elston. 2003. J. Theo. Biol. 221:491–511.
- Wang, H. and T. C. Elston. 2007. J. Stat. Phys. 128:35–76.
- Gillespie, D. T. 2007. Annu. Rev. Phys. Chem. 58:35–55.
- Friddle, R. W. in: bizzarri12 ().
- Makarov, D. E. 2016. J. Chem. Phys. 144:030901–13.
- Konda, S. S. M., J.N. Brantley, C.W. Bielawski and D.E. Makarov. 2011. J. Chem. Phys. 135:164103–8.
- Ghanti, D. 2017. Ph.D. Thesis, IIT Kanpur (in preparation).
- Hendricks, A. G., J.E. Lazarus, E. Perlson, M.K. Gardner, D.J. Odde, Y.E. Goldman and E.L.F. Holzbaur. 2012. Curr. Biol. 22:632–637.
- Joglekar, A. P. and A.J. Hunt. 2002. Biophys. J. 83:42–58.
- Shtylla, B. and J. P. Keener. 2011. SIAM J. Appl. Math. 71:1821–1848.
- Hill, T. 1985. Proc. Natl. Acad. Sci. U.S.A. 82:4404–4408.
- Waters, J. C., T.J. Mitchison, C.L. Rieder and E. D. Salmon. 1996. Mol. Biol. Cell. 7:1547–1558.
- Marshall, W. F., J. F. Marko, D. A. Agard and J. W. Sedat. 2001. Curr. Biol. 11:569–578.
- Erdmann, T. and U. S. Schwarz. 2004. J. Chem. Phys. 121:8997–9017.
- Williams, P. M. 2003. Ana. Chimi. Acta. 479:107–115.
- Seifert, U. 2000. Phys. Rev. Lett. 84:2750–2753.