# Distribution of directional change as a signature of complex dynamics

See allHide authors and affiliations

Contributed by Stuart A. Rice, October 18, 2013 (sent for review July 30, 2013)

## Significance

Since Einstein’s seminal work in 1905, the main means of characterizing stochastic processes has been the mean square displacement (MSD). However, this order parameter fails to capture many features of dynamics at the forefront of science today, ranging from glassy relaxation to active transport in biological cells. Although there have been several studies seeking to go beyond the MSD, such studies have not made full use of the information available in individual trajectories in two (or more) dimensions, as are now commonly obtained in particle tracking experiments. Here, we introduce an approach that quantifies directional properties of complex motions and discover striking correlations in a number of condensed phase systems.

## Abstract

Analyses of random walks traditionally use the mean square displacement (MSD) as an order parameter characterizing dynamics. We show that the distribution of relative angles of motion between successive time intervals of random walks in two or more dimensions provides information about stochastic processes beyond the MSD. We illustrate the behavior of this measure for common models and apply it to experimental particle tracking data. For a colloidal system, the distribution of relative angles reports sensitively on caging as the density varies. For transport mediated by molecular motors on filament networks in vitro and in vivo, we discover self-similar properties that cannot be described by existing models and discuss possible scenarios that can lead to the elucidated statistical features.

Complex dynamics often emerge from ensembles of interacting constituents. Trajectories that are obtained by tracking individual constituents contain information beyond the evolution of ensemble properties, and these data can thus reveal new mechanistic features of the system studied. Examples cut across disciplines and include quantum dots (1), colloidal beads (2), features in cells (3, 4), fish in schools (5), birds in flocks (6), and primates in social groups (7, 8). These data (individual trajectories) demand theoretical frameworks for characterizing and interpreting them.

The standard reporter for different forms of motion is the mean square displacement (MSD)where brackets and overlines denote ensemble and time averages, respectively. In simple Brownian motion (9), the MSD grows linearly with the separation in time between two observation points (the lag time, Δ) and does not depend on the amount of data included in averages (the measurement time, *T*)—i.e., there is ergodicity. Anomalous (i.e., non-Brownian) dynamics can arise from correlations in the walk steps. Correlations in the step sizes [e.g., as in fractional Brownian motion (FBM)] (10) give rise to nonlinear scaling of the MSD with lag time, retaining ergodicity (11). By contrast, a power-law distribution of dwell times [e.g., as in a continuous time random walk (CTRW)] (12) is associated with linear scaling with lag time but nonergodicity (13, 14). The two types of correlations can exist together (4, 15).

Despite the success of the MSD as an order parameter for dynamics, it is essentially a 1D measure. We expect random walks in two and more dimensions to contain information beyond the MSD, and various alternative analyses have been suggested (16) (*Conclusions*). In this paper, we introduce a statistical measure of such information. Specifically, we consider the relative angle, which quantifies the direction of motion over successive time intervals. We show that different models of stochastic processes give rise to different distributions of relative angles and how the intervals can be varied to probe contributing time scales. We apply our order parameter to 2D experimental data obtained for mesoscopic systems. We examine colloidal suspensions at two densities and show that the distribution of relative angles is a sensitive means of detecting and quantifying caging. Two different cytoskeletal systems are considered: insulin-containing vesicles (granules) in a pancreatic cell line (4) and in vitro mixtures of purified myosin motors and actin filaments. These active systems exhibit a common signature of directional motion that cannot be understood in terms of existing models. We propose possible scenarios that could lead to the observed statistics. Together, these examples show that the distribution of relative angles is a straightforward probe of complex dynamics that provides information beyond the MSD.

## Construction of the Relative Angle

In this section, we introduce our order parameter for quantifying directional change. As shown in Fig. 1, each trajectory comprises a set of recorded positions of the particles . Given these data, we form the vectors that connect positions separated by Δ steps, . We define the relative angle as the angle between and In analogy to Eq. **1**, the parameter Δ plays the role of a lag time; it controls the degree of temporal coarse-graining. For a given Δ, we compute across all times in each trajectory and then across all particles. We build a histogram of these values and normalize such that it integrates to one. The resulting probability density function, , serves as our order parameter. We repeat this procedure for different Δ and examine how the order parameter changes. The quantity that we calculate, Eq. **2**, is reminiscent of a velocity autocorrelation function (17), but it differs in that the **V** vectors are not instantaneous velocities and, importantly, we do not ensemble average.

## Dictionary of Model Signatures

The variation in with Δ can serve as a signature of different types of motion. As a simple example, consider an inertial Langevin dynamics (17)Here *m* is the mass of the particle, is a drag, and is a 2D noise vector where and are Gaussian distributed uncorrelated noises , where is the Dirac δ-function, is the Kronecker δ-function, and is a constant noise strength (the brackets indicate ensemble averaging). We integrate this equation of motion to obtain a 2D trajectory. When (Fig. 2, *Upper*), the distribution is peaked at , which is equivalent to owing to the periodicity of the angle. The U-shape in this representation is a signature of inertial motion: at short times the particle tends to continue in the same direction. As Δ increases, the noise randomly shifts the direction of travel, causing the distribution to flatten steadily (Fig. 2, lower three panels). For sufficiently large temporal coarse-graining, successive vectors between sampled points are completely independent and the distribution is that of a simple random walk (Fig. 2, Δ = 300). The distribution remains flat for all further times, a signature of the self-similarity of diffusion (18).

We can similarly understand the signatures of motions that are influenced by various potentials. In Fig. 3*A*, we show for various Δ for diffusion of a particle in a box. At short times, the particle does not feel the boundary and the motion is just a simple random walk. As the temporal coarse-graining increases, a peak grows at . This peak indicates that the successive vectors between sampled points are anticorrelated in direction—the reflecting boundaries force the particle to reverse its motion. Eventually the distribution converges to a stationary profile, again characteristic of a self-similarity. Analogous behavior is observed for other confining potentials, with the shape of the peak at dependent on the details of the specific form.

In Fig. 3*B*, we show the behavior of of a particle under the influence of a harmonic potential with a steadily translating minimum (i.e., a steered dynamics, as in an optical trap experiment). At finite times, the distribution is peaked at due to the trap, just like the situation immediately above. However, as the confining potential moves, the particle dynamics are positively correlated overall, causing the peak to break in two and ultimately give rise to an inertial signature (compare Fig. 3*B*, Δ = 24000, with Fig. 2, Δ = 1).

Complex media can give rise to more complex dynamics with correlated steps, as discussed in the Introduction. For example, transport in a viscoelastic environment is expected to give rise to FBM (19). FBM is a stationary Gaussian process where the spatial position is correlated such that , where *H* is the Hurst exponent. When , the motion is subdiffusive because the steps are negatively correlated, and we find to be akin to that of a confining potential (Fig. 4*A*; and 0.4)—i.e., there is a tendency for reversals. In contrast, when , the motion is superdiffusive because the steps are positively correlated and we find an interial-like signature (Fig. 4*A*; ).

As already discussed, alternative anomalous dynamics arise from a power-law distribution of dwell times for steps (i.e., a CTRW). A pure CTRW signature is shown in Fig. 4*B* . The directional anticorrelation peak at decreases with Δ. When the step sizes are described by FBM and the dwell times are described by CTRW (i.e., when ergodic and nonergodic processes coexist), we say that FBM is subordinated to CTRW (see refs. 4, 15, and 20 for examples). Examples of for such dynamics are shown in the remaining panels in Fig. 4*B*. The case is quite similar to the CTRW (*H* = 0.5), but, in the case (i.e., when there are positively correlated steps, as in the case of active transport), there is a progression from a confined signature to an inertial one with increasing temporal coarse-graining. However, note that, in contrast to Fig. 3*B*, the peak at emerges without the peak at splitting. In other words, the qualitative shape of the profiles at intermediate values of Δ can be used to distinguish different motions.

## Application to Experimental Data

We examine for particle tracking data from three different experimental situations.

### Colloidal Suspensions.

The first system is driven only by thermal energy. It is a colloidal suspension of 1.58 ± 0.04-μm-diameter silica spheres in water between parallel glass plates that are separated by 1.76 ± 0.05 μm (21). The data are those for packing fractions of and from ref. 21. We expect the colloidal particles to move freely at the former *ϕ* and to be close to jammed at the latter *ϕ*. Consistent with these ideas, we find that for is almost flat (Fig. 5*B*, *Upper*), whereas that for exhibits a pronounced peak at (Fig. 5*B*, *Lower*), similar to the profile for a particle in a box (Fig. 3*A*). The peak in Fig. 5*B*, *Lower* decays slightly as the temporal coarse-graining interval is lengthened; we interpret this decay to reflect the fact that particles are not completely confined by their neighbors—they eventually escape their neighbors and the system mixes after many such events. The dependence on Δ provides quantitative information about the time scale of such dynamics. Caging and escape were identified previously in colloidal suspensions (21, 22) and can be directly observed for specific trajectories. Likewise, backscattering of neighboring particles in simulations has long been visualized through the velocity autocorrelation function (23). However, provides a much more sensitive statistical probe: for shows a peak at when the corresponding MSD is already linear (Fig. 5*A*).

### Insulin Granules.

More complex dynamics can arise in systems with active elements, such as molecular motors. Here, we study the transport of insulin-containing vesicles (granules). In addition to active elements, these dynamics are influenced by fluctuations of the crowded cellular environment. Analysis of the MSD recently indicated that these granules combine ergodic and nonergodic random walk processes (4), and we showed that a hybrid model that subordinated a FBM with to a CTRW with dwell time distribution (with ) accounted for a host of statistics of the motion (4). This model provides a simple physical mechanism for obtaining characteristic insulin secretion profiles comprised of a burst followed by sustained release (see ref. 4 for further discussion of the biological implications). In the present study, we obtain similar scalings with different fluorescent constructs [proinsulin-enhanced green fluorescent protein (EGFP) (24) vs. syncollin-EGFP (4)]. Specifically, the dwell-time distribution exponent *β* is obtained from (Fig. 6, *Upper Right*) (13) and then the Hurst exponent *H* is obtained from (Fig. 6, *Upper Left*) and (25).

We plot for various temporal coarse-grainings in Fig. 6. At small Δ, there are two peaks: a larger one at and a smaller one at . The former shrinks and the latter grows with Δ, ultimately becoming equal in size. Given the model in ref. 4, it is worth comparing the profile in Fig. 6 with those in Fig. 4. As discussed above, CTRW gives rise to an effective confinement, leading to a peak at (Fig. 4*B*), which decays to a uniform distribution of *θ* as particles escape their traps, in analogy to the colloidal suspension discussed above. The profile for FBM depends on the Hurst exponent. FBM with (subdiffusive or negatively correlated steps) exhibits a peak at (Fig. 4*A*), whereas FBM with (superdiffusive or positively correlated steps) exhibits a peak at . The subordinated model in ref. 4 combines CTRW with FBM with and thus captures the peak at , including its decay. However, that model fails to reproduce the peak at that is characteristic of active transport. This limitation of the simple model in ref. 4 is not surprising, as a minority of the granules are superdiffusive (4).

We can reproduce the observed behavior (Fig. 7*A*, *Upper*) by expansion of the subordinated model to include Hurst exponents above and below 0.5 for different particles (Fig. 7*A*, *Lower*), as could arise from spatial heterogeneity within the cells. The main limitation of this approach is that is very sensitive to the distribution of *H*; the amplitude of the variation in is also restricted to the range shown. An alternative way to account for the observed shape of is to assume that the motion is anisotropic, again from spatial heterogeneity. We show the behavior for regular diffusion with two different ratios of diffusion coefficients in the Cartesian directions *x* and *y* in Fig. 7*B*. In this case, the profile can be tuned to make the peaks more pronounced. We do not detect global differences in properties of the motion, but it is possible for local differences to exist. A deeper understanding of the structure of for this system requires exploration of microscopic models with explicit molecular features, which is beyond the scope of the present study. Rather, the key point is that the relative angle distribution reveals directional correlations and self-similarity that are not evident in the MSD (or measures of angular correlation that are averages) (16), and these features can be used to further constrain models of transport.

### Reconstituted Filament-Motor System.

The final system that we consider comprises a mixture of filamentous actin (F-actin), myosin thick filaments (motors), and α-actinin (passive) cross-linkers on a model membrane substrate. This system with motor densities μm^{−2} was used recently to show that F-actin buckling breaks the symmetry between extensile and compressive forces to generate contractility in cytoskeletal networks lacking sarcomeric organization (26), as suggested previously (27, 28). Here, we study the motion of the motors when they are at a lower density (0.04 μm^{−2}) such that they do not perturb the F-actin structure.

In contrast to the cellular system discussed above (4), the MSD varies linearly with lag time for this system (Fig. 8, *Upper Left*). The MSD thus suggests that the dynamics of this system are simpler than those of the insulin granules. However, the MSD decays with measurement time with power-law scaling (Fig. 8, *Upper Right*), as in the data for the insulin granules. We interpret this aging as evidence for local trapping. This picture and, more precisely, the MSD statistics are consistent with pure CTRW motion and glassy dynamics.

Nevertheless, for the actomyosin system resembles that for the pancreatic cell line, including the dependence on the temporal coarse-graining interval (Fig. 8, *Lower*). Again, this profile can be fit by a model with anisotropic diffusion, and the amplitudes of the peaks sets the difference in diffusion coefficient (here, ; cf. Fig. 7*B*). This system thus provides an opportunity for investigating the directional correlations in a context in which all of the contributing elements are known. One possible scenario is that the heads of each myosin thick filament bind to different actin filaments and undergo a “tug-of-war” that ultimately gives way to movement in one direction. However, the fact that the peak at does not grow without bound indicates that the directed motion cannot be sustained. The stationarity of the distribution (within experimentally accessible time scales) imposes strong restrictions on the microscopic dynamics and the structure of the filament network. Our goal is not to investigate these dynamics here. Rather, it is to show that the relative angle analysis reveals common features of the filament-motor dynamics (Figs. 6 and 8), in contrast to the MSD, which shows very different behaviors for the two systems (diffusive vs. subdiffusive scaling).

## Conclusions

We have introduced an order parameter for quantifying directional motion in stochastic trajectories, the distribution of relative angles for different temporal coarse-grainings, , and we have shown that it gives unique insights into particle tracking data. There is a long history of methods that seek to go beyond the MSD to characterize dynamics (29, 30), particularly since the advent of single-molecule tracking experiments. Recent innovations include comparing exchange and persistence time distributions to detect glassy behavior (31), *p*-variation (32), the mean maximal excursion method for anomalous diffusion (33), and the diffusivity distribution (34). However, these methods focus on distributions of extents of changes (e.g., distances traveled) and their scaling with time. By contrast, takes advantage of the fact that multidimensional measurements contain additional information about direction of motion. We find this order parameter to be a very sensitive measure of changes in dynamics, and it provides quantitative information about the time scales over which different dynamics contribute. These features are illustrated by our analysis of colloidal suspensions: clearly shows the transition from unconfined to confined dynamics with packing fraction, and it yields insights into the time scales for particles to escape their neighbors and mix. One recent study of lipid droplets in live cells did examine direction of motion but characterized it through the time correlation function for relative angles (16), which is an ensemble average. Although it has been shown that the velocity autocorrelation can yield information about subtle features, such as anharmonicity (30), in general, averaging obscures the heterogeneity in a population and the self-similarity of the motion. Our single-particle analysis here reveals such features in the filament-motor dynamics considered and thus provides constraints for the design of microscopic models of these systems. As 3D data become available, additional (dihedral) angles can be studied. Beyond that, an outstanding challenge is development of a theoretical framework for the order parameter and its behavior.

## Acknowledgments

We thank Monika Kauer for helpful discussions and a critical reading of the manuscript. This work was supported by the WM Keck Foundation, Novo Nordisk, National Science Foundation (NSF) Division of Materials Research (DMR) Materials Research Science and Engineering Center (MRSEC) Grant 0820054, and National Institutes of Health (NIH) R01-DK092616 and P30-DK020595. S.B. wishes to acknowledge the hospitality of the Aspen Center for Physics, which is partially supported by NSF Division of Physics (PHY) Grant 1066293.

## Footnotes

- ↵
^{1}To whom correspondence may be addressed: E-mail: s-rice{at}uchicago.edu or dinner{at}uchicago.edu.

Author contributions: S.B., S.M.A.T., S.A.R., N.F.S., and A.R.D. designed research; S.B., S.M.A.T., T.H., and M.P.M. performed research; L.H.P. and M.L.G. contributed new reagents/analytic tools; and S.B., N.F.S., and A.R.D. wrote the paper.

The authors declare no conflict of interest.

## References

- ↵
- ↵
- Wang B,
- Anthony SM,
- Bae SC,
- Granick S

- ↵
- ↵
- Tabei SM,
- et al.

- ↵
- ↵
- ↵
- ↵
- Song C,
- Qu Z,
- Blumm N,
- Barabási AL

- ↵Einstein A (1905) On the movement of small particles suspended in stationary liquids required by the molecular-kinetic theory of heat.
*Ann Physik*17(8):549–560. - ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Weigel AV,
- Simon B,
- Tamkun MM,
- Krapf D

- ↵
- ↵
- Kubo R,
- Toda N,
- Hashitsume N

- ↵
- Bunde A,
- Havlin S

- Stanely HE

- ↵
- Goychuk I

- ↵
- ↵
- ↵
- Weeks ER,
- Crocker JC,
- Levitt AC,
- Schofield A,
- Weitz DA

- ↵
- ↵
- Rajan S,
- et al.

- ↵
- ↵
- Murrell MP,
- Gardel ML

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Applied Physical Sciences

- Biological Sciences
- Biophysics and Computational Biology