$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
To determine successful preparation of active composites (Figure 1), and to characterize their dynamics and structure, a laser scanning fluorescence microscope with at least two fluorescence channels is used to visualize the actin filaments and microtubules simultaneously (Figure 2 and Figure 6). All actin filaments and microtubules in the composites are sparsely labeled, rather than doping in tracer bright filaments, as is often done in in vitro studies. This method ensures that the measured dynamics and structure are representative of the composite itself rather than the tracers which are formed under different conditions than the composites. For this reason, individual actin filaments and microtubules cannot typically be resolved, rather images portray mesoscale network structure (Figure 2 and Figure 6).
This labeling approach was optimized for spatial image autocorrelation (SIA) and differential dynamic microscopy (DDM) analyses that examine the dynamics and structure in reciprocal Fourier space (Figure 4, Figure 5, and Figure 8)52,53,54,55. Particle-image velocimetry (PIV) can also be used to depict and characterize dynamics and flow fields (Figure 3 and Figure 7), but it requires pixel-binning (lower spatial resolution) and larger lag-time increments (lower temporal resolution) than SIA and DDM to eliminate erroneous vectors that arise from noise in the dense, low-signal images. Nevertheless, PIV is recommended for qualitative examination of flow fields and corroboration of DDM results (Figure 4 and Figure 8)26,50.
Sample characterization of the described networks using these analyses (i.e., DDM, SIA, PIV) is provided to aid researchers in adopting similar analyses to benchmark and characterize their samples. However, detailed descriptions of these techniques are outside the scope of this work. For detailed descriptions of how to perform DDM on these and other similar systems, including user-friendly Python code, refer to previous works17,26,49,50 and the references there within. For details regarding how to perform SIA and PIV on the systems described here, the reader is directed to previous works17,50.
Several controls, described below, should be done to ensure that the composites are functioning as expected. A composite without myosin or kinesin should appear essentially static with minimal thermal fluctuations or drift. Actin filaments and microtubules should appear co-entangled and homogeneously distributed, with minimal bundling, aggregation, or phase separation of actin and microtubules throughout a field of view of ~200 µm x 200 µm (Figure 2, far left)17. One should expect a similar result for composites that contain myosin but are not exposed to 488 nm light (to deactivate the blebbistatin).
Upon incorporation of myosin and exposure to 488 nm light, the composites undergo contraction that is largely isotropic and similar for actin and microtubules, as seen in microscope images taken before and after myosin activity (Figure 2), as well as corresponding PIV flow fields for varying times during activity (Figure 3). To determine if the motion is ballistic, diffusive, subdiffusive, etc., the characteristic decorrelation time τ(q) determined from DDM is evaluated as a function of wave vector (i.e., reciprocal space). See as described in detail previously17,26,49. Figure 4 also demonstrates how to use DDM to characterize these composites. Power-law scaling τ(q)~1/vqβ, with β = 1, indicates ballistic motion with speed v. For reference, β = 2 represents diffusive dynamics with v being the diffusion coefficient. All active composites exhibit ballistic scaling (Figure 4A) with speeds that are tuned by the concentrations of actin and myosin (Figure 4B), and that can vary in time during activity, either accelerating or decelerating (Figure 4C,D).
Network restructuring and clustering, visible in Figure 2 and more evident for higher actin and myosin concentrations, can be characterized using SIA, as depicted in Figure 5, and described previously17,48,50. Briefly, a correlation length ξ, which is a measure of the characteristic size of features in an image, can be determined by fitting each spatial intensity autocorrelation curve g(r) to an exponential function of distance r between pixels. Larger g(r) peaks that persist for longer distances indicate larger structural features (i.e., bundling, clustering of the individual filaments). As shown in Figure 5, for higher actin fractions and myosin concentrations, significant restructuring and aggregation is reflected in the increase in ξ over time.
The viscoelastic properties and nonlinear mechanical response of the active composites can also be measured using optical tweezers microrheology (OTM). However, protocols and representative results for these experiments are outside the scope of this work. Interested readers are referred to previous works48,56 which thoroughly describe how to perform OTM measurements and the expected results.
Using the same program of experimental and analysis tools described above, the following section describes how the dynamics and structure change when kinesin motors and biotin-NA crosslinkers are incorporated into the composites (Figure 6, Figure 7, and Figure 8). Figure 6 shows representative confocal images of composites driven by either kinesin-only (K) or kinesin and myosin (K+M), with and without passive crosslinking (XL) of actin filaments or microtubules.
Incorporating kinesin into composites initially results in similar dynamics and restructuring as myosin-driven composites as seen in the top row of Figure 7 (Class 1). However, the dynamics typically transition to large-scale anisotropic flow (Figure 7 middle row, Class 2), acceleration, and deceleration (Figure 7 bottom row, Class 3). These characteristics couple with mesoscale clustering and aggregation after 5-30 min (Figure 6 and Figure 8B). PIV-generated flow fields and temporal color maps shown in Figure 7 depict examples of isotropic restructuring (Class 1, top panel), directed flow (Class 2, middle panels), and bi-directional acceleration (Class 3, bottom panels).
Speeds of actin and microtubules at varying time points during activity, determined via fits to τ(q) curves, illustrate acceleration followed by deceleration (Figure 8), which depends on crosslinking. As also shown in Figure 8, when both motor proteins are incorporated, the dynamics are actually slower than kinesin-only composites, and there is delayed onset of mesoscale flow. Myosin also supports more homogeneous interpenetration of actin and microtubule networks throughout the duration of activity, as well as less aggregation and restructuring. These effects can be seen in the images in Figure 6 and are quantified by the time-varying correlation lengths computed via SIA, which are generally smaller in the presence of myosin (Figure 8B).

Figure 1. Design and characterization of active actin-microtubule composites with multiple force-generating motors and passive crosslinkers. (A) Actin monomers and tubulin dimers are co-polymerized at molar concentrations cA and cT of 0.73-11.6 μM and molar fractions of actin ΦA = cA / (cA + cT) = 0, 0.25, 0.5, 0.75, and 1, to form co-entangled networks of actin filaments (green) and microtubules (red). Passive crosslinking is achieved using NA to link biotinylated actin filaments (Actin XL) or microtubules (MT XL) at crosslinker:protein molar ratios of RA = 0.01-0.08 and RMT = 0.001-0.01 for actin and microtubules, respectively. Myosin-II mini-filaments (purple) and kinesin clusters (orange), at concentrations of cM = 0.12 - 0.48 μM and cK = 0.2 - 0.7 μM, push and pull on the filaments to drive the composites out of steady-state. (B) Schematic of formulation space. Myosin II mini-filaments (M), kinesin clusters (K), or both motors (K+M) are incorporated into composites with no passive crosslinkers (No XL), actin-actin crosslinks (Actin XL), and microtubule-microtubule crosslinks (MT XL). All cartoons are not drawn to scale. Please click here to view a larger version of this figure.

Figure 2. Two-color confocal imaging of myosin-driven cytoskeleton composites with varying myosin concentrations cM and molar actin fractions ΦA. (A) 256 x 128 square-pixel (212 x 106 μm2) two-color confocal microscopy images show how composites of actin filaments (green) and microtubules (red) are rearranged via myosin motor activity. No kinesin motors or passive crosslinkers are present. In each panel, images taken at the beginning (left, before) and end (right, after) of the 45 min myosin activation (via illumination with 488 nm light to deactivate blebbistatin) are shown. Panels are ordered by increasing molar concentration of myosin (cM), going from left to right, and increasing molar fraction of actin (ΦA), going from top to bottom. The colors outlining each panel match the color coding used in Figure 4 and Figure 5. Scale bars are 50 μM. To best capture dynamics and structure for analysis, we use frame rates of 1-5 fps, ROIs with 50-250 µm sides, and time-series durations of 5-45 min, depending on the rate of contraction and rearrangement. Panels in which the before and after images look similar indicate minimal restructuring, as seen in the pink, magenta, and cyan panels. Small-scale clustering, evidenced by increased heterogeneity and the presence of bright punctate features, can be seen in the orange, green, and red panels. Large-scale contraction, seen as a uniformly shrinking network, is evident in the blue and purple panels. This figure has been modified from reference17. Please click here to view a larger version of this figure.

Figure 3. Particle image velocimetry (PIV) shows that actomyosin activity triggers coordinated contractile dynamics of actin and microtubules in co-entangled composites. PIV flow fields for actin (top row) and microtubules (bottom row) in a myosin-driven composite with (ΦA, cM) = (0.5, 0.24) at increasing times during a 6 min time-series. Flow fields were generated using the Fiji/ImageJ PIV plugin with a lag-time of 20 s and 2 pixel x 2 pixel binning. Both actin and microtubules show consistent motion directed towards the center region of the field-of-view throughout the movie duration. Scale bars in all images are 50 μm. Different arrow colors correspond to different speeds as indicated in the color scale to the right of vector fields. This figure has been modified from reference26. Please click here to view a larger version of this figure.

Figure 4. Time-resolved differential dynamic microscopy (DDM) measures the rate and type of motion of actin and microtubules in active composites. (A) DDM is performed on microtubule (top, open symbols) and actin (bottom, filled symbols) channels of time-series to determine characteristic decay times τ vs wavenumber q for both actin (filled symbols) and microtubules (open symbols) as described previously17,26. All curves follow τ ~ q-1 scaling, indicating ballistic motion, with speeds v that are determined via fits to τ(q) = (vq)-1. Faster speeds correspond to smaller τ(q) values for any given q. Symbol colors and shapes correspond to (ΦA, cM) combinations shown in B. (B) Contraction speeds v are determined via fits to τ(q) curves shown in A, which are averaged over all lag-times for the duration of each 45 min time-series. (C) Time-resolved DDM (trDDM) quantifies how the dynamics vary over time by evaluating τ(q) for actin (filled symbols, left) and microtubules (open symbols, right) for consecutive 6 min intervals (denoted by different shades of the same color) during the 45 min activation time. trDDM is performed for each (ΦA, cM) combination (denoted by different symbols and colors) as described in the legend in lower right. τ(q) curves shown in C follow similar scaling and trends as those in A but also show time-dependence for certain (ΦA, cM) compositions, most notably for ΦA = 0.75. (D) Contraction speeds for actin filaments (closed symbols) and microtubules (open symbols) are determined from fits to corresponding τ(q) curves. Error bars in all plots represent the standard error of values across three to five replicates. This figure has been modified from reference17. Please click here to view a larger version of this figure.

Figure 5. Spatial image autocorrelation (SIA) analysis quantifies the motor-driven restructuring of active cytoskeletal composites. (A) Autocorrelation g(r) for the microtubules at the beginning (left, t = 0 min, dark shades) and end (right, t = 42 min, light shades) of the experiment for (ΦA, cM) formulations listed in the legend. Inset: example fits of data to
at the initial and final times for (ΦA, cM) = (0.75, 0.12). (B) Average correlation lengths ξ for actin (closed symbols) and microtubules (open symbols) for each (ΦA, cM) determined via exponential fits of each g(r) curve, as shown in the inset in A. Data is divided into those that exhibit minimal (left) versus substantial (right) restructuring. Error bars in A and B represent the standard error across three to five replicates. This figure has been modified from reference17. Please click here to view a larger version of this figure.

Figure 6. Incorporating kinesin motors and passive crosslinkers into active composites to increase programmability and expand the phase space of dynamics and structure. (A) Two-color confocal images of actin (green) and microtubules (red) in active composites show complex formulation-dependent restructuring over time (listed in min). The five images in each row correspond to five frames of a 2000 frame time-series acquired for a composite driven by kinesin (K, rows 1, 3, 5) or kinesin and myosin (K+M, rows 2, 4, 6), and including either no passive crosslinkers (No XL, rows 1, 2), actin-actin crosslinks (Actin XL, rows 3, 4), or microtubule-microtubule crosslinks (MT XL, rows 5, 6). Scale bars are all 50 µm. Outline colors match color scheme in Figure 8. (B) Separate actin and microtubule fluorescence channels for the kinesin-only composites show varied structures with both actin-MT co-localization and micro-phase separation. The images shown are for composites with cA = 2.32 μM, cT = 3.48 μM, cK = 0.35 μM, cM = 0.47 μM (rows 2, 4, 6), RA = 0.02 (rows 3, 4), and RMT = 0.005 (rows 5, 6). All composites begin with uniformly distributed interpenetrating networks of actin and microtubules (column 1). Kinesin-driven composites without crosslinkers (row 1) form loosely connected amorphous clusters that are MT-rich. Actin co-localizes in the centers of these aggregates initially but then is squeezed out of the MT-rich regions that continue to contract and disconnect from one another. Actin-actin crosslinking (row 3) hinders this microscale actin-MT separation, and instead MT-rich aggregates are connected via long strands of actin. Actin crosslinking also enables slow uptake of actin into the MT-rich regions, such that the composite becomes a connected network of co-localized actin and MT clusters. Microtubule crosslinking (row 5) leads to amorphous clustering of MTs that coalesce over time, resulting in larger scale phase separation of actin and MTs. Adding myosin (rows 2, 4, 6) reduces kinesin-driven de-mixing and restructuring. Without crosslinkers (row 2), composites show little rearrangement over the course of hours. Crosslinking increases restructuring and co-localization of actin and microtubules (rows 4, 6). Specifically, when microtubules are crosslinked (row 6), there is significant interpenetration and reorganization into web-like networks of fibers. This figure has been modified from reference50. Please click here to view a larger version of this figure.

Figure 7. PIV shows that active composites exhibit three classes of spatiotemporally distinct flow fields. (A) PIV flow fields for the first (ti) and last (tf) frames of three representative time-series, showing the different dynamical classes that composites shown in Figure 6 exhibit. PIV flow fields for microtubules (top) and actin (bottom) for class 1 (top, purple), class 2 (middle, orange), and class 3 (bottom, magenta) example videos, with arrow colors corresponding to the universal speed scale at the bottom, and the grayscale colormap showing the spatial speed distribution, normalized separately for each flow field according to the scale shown at bottom. Scale bars are all 50 μM. (B) Angular distributions of velocity vectors from A (in units of radians) with listed initial and final standard deviations σi and σf. (C) Temporal color maps for the videos analyzed in A and B show the frame-to-frame position of each pixel relative to its starting point. Class 1 maps show small-scale random motion; class 2 maps depict fast unidirectional motion with minimal spatial or temporal variation; class 3 maps exhibit features of both class 1 and 2. This figure has been modified from reference50. Please click here to view a larger version of this figure.

Figure 8. DDM and SIA measure the time-varying dynamics and structure of two-motor actin-microtubule composites. (A) Speeds for composites described in Figure 6 and Figure 7, measured via DDM, show acceleration and deceleration of composites, programmed by crosslinking and myosin activity. Speeds of microtubules (MT, closed circles) and actin (A, open circles) are plotted as a function of activity time in composites with no crosslinking (top, blue), actin crosslinking (middle, green), microtubule crosslinking (bottom, red), without myosin (K, darker shades), and with myosin (K+M, lighter shades). For class 3 cases, which have two speeds, the slower speed is indicated by a star. Data points enclosed by dashed black circles correspond to maximum speed vmax for each formulation. Error bars (most too small to see) are the standard error over the power-law fits of the corresponding τ(q). (B) Structural correlation lengths ξ, determined via SIA, versus activity time, for the same set of time-series evaluated in A. Each data point is an average of the correlation lengths determined for the first and last frame of the corresponding time-series. In general, ξ increases in time for both actin and microtubules in all composites systems, and composites driven solely by kinesin have greater correlation lengths than those in which myosin is also present. Data points in A and B that correspond to the three time-series analyzed in Figure 7 are circled in the corresponding class color (1 = purple, 2 = orange, 3 = magenta). This figure has been modified from reference50. Please click here to view a larger version of this figure.