Distinct regions of the kinesin-5 C-terminal tail are essential for mitotic spindle midzone localization and sliding force.

Kinesin-5 motor proteins play essential roles during mitosis in most organisms. Their tetrameric structure and plus-end-directed motility allow them to bind to and move along antiparallel microtubules, thereby pushing spindle poles apart to assemble a bipolar spindle. Recent work has shown that the C-terminal tail is particularly important to kinesin-5 function: The tail affects motor domain structure, ATP hydrolysis, motility, clustering, and sliding force measured for purified motors, as well as motility, clustering, and spindle assembly in cells. Because previous work has focused on presence or absence of the entire tail, the functionally important regions of the tail remain to be identified. We have therefore characterized a series of kinesin-5/Cut7 tail truncation alleles in fission yeast. Partial truncation causes mitotic defects and temperature-sensitive growth, while further truncation that removes the conserved BimC motif is lethal. We compared the sliding force generated by cut7 mutants using a kinesin-14 mutant background in which some microtubules detach from the spindle poles and are pushed into the nuclear envelope. These Cut7-driven protrusions decreased as more of the tail was truncated, and the most severe truncations produced no observable protrusions. Our observations suggest that the C-terminal tail of Cut7p contributes to both sliding force and midzone localization. In the context of sequential tail truncation, the BimC motif and adjacent C-terminal amino acids are particularly important for sliding force. In addition, moderate tail truncation increases midzone localization, but further truncation of residues N-terminal to the BimC motif decreases midzone localization.

The kinesin-5 protein Cut7 moves bidirectionally on fission yeast spindles with activity that increases in anaphase.

Kinesin-5 motors are essential to separate mitotic spindle poles and assemble a bipolar spindle in many organisms. These motors crosslink and slide apart antiparallel microtubules via microtubule plus-end-directed motility. However, kinesin-5 localization is enhanced away from antiparallel overlaps. Increasing evidence suggests this localization occurs due to bidirectional motility or trafficking. The purified fission-yeast kinesin-5 protein Cut7 moves bidirectionally, but bidirectionality has not been shown in cells, and the function of the minus-end-directed movement is unknown. Here, we characterized the motility of Cut7 on bipolar and monopolar spindles and observed movement toward both plus- and minus-ends of microtubules. Notably, the activity of the motor increased at anaphase B onset. Perturbations to microtubule dynamics only modestly changed Cut7 movement, whereas Cut7 mutation reduced movement. These results suggest that the directed motility of Cut7 contributes to the movement of the motor. Comparison of the Cut7 mutant and human Eg5 (also known as KIF11) localization suggest a new hypothesis for the function of minus-end-directed motility and spindle-pole localization of kinesin-5s.

Motor guidance by long-range communication on the microtubule highway.

Coupling of motor proteins within arrays drives muscle contraction, flagellar beating, chromosome segregation, and other biological processes. Current models of motor coupling invoke either direct mechanical linkage or protein crowding, which rely on short-range motor-motor interactions. In contrast, coupling mechanisms that act at longer length scales remain largely unexplored. Here we report that microtubules can physically couple motor movement in the absence of detectable short-range interactions. The human kinesin-4 Kif4A changes the run length and velocity of other motors on the same microtubule in the dilute binding limit, when approximately 10-nm-sized motors are much farther apart than the motor size. This effect does not depend on specific motor-motor interactions because similar changes in Kif4A motility are induced by kinesin-1 motors. A micrometer-scale attractive interaction potential between motors is sufficient to recreate the experimental results in a biophysical model. Unexpectedly, our theory suggests that long-range microtubule-mediated coupling affects not only binding kinetics but also motor mechanochemistry. Therefore, the model predicts that motors can sense and respond to motors bound several micrometers away on a microtubule. Our results are consistent with a paradigm in which long-range motor interactions along the microtubule enable additional forms of collective motor behavior, possibly due to changes in the microtubule lattice.

Chiral self-sorting of active semiflexible filaments with intrinsic curvature.

Many-body interactions in systems of active matter can cause particles to move collectively and self-organize into dynamic structures with long-range order. In cells, the self-assembly of cytoskeletal filaments is critical for cellular motility, structure, intracellular transport, and division. Semiflexible cytoskeletal filaments driven by polymerization or motor-protein interactions on a two-dimensional substrate, such as the cell cortex, can induce filament bending and curvature leading to interesting collective behavior. For example, the bacterial cell-division filament FtsZ is known to have intrinsic curvature that causes it to self-organize into rings and vortices, and recent experiments reconstituting the collective motion of microtubules driven by motor proteins on a surface have observed chiral symmetry breaking of the collective behavior due to motor-induced curvature of the filaments. Previous work on the self-organization of driven filament systems have not studied the effects of curvature and filament structure on collective behavior. In this work, we present Brownian dynamics simulation results of driven semiflexible filaments with intrinsic curvature and investigate how the interplay between filament rigidity and radius of curvature can tune the self-organization behavior in homochiral systems and heterochiral mixtures. We find a curvature-induced reorganization from polar flocks to self-sorted chiral clusters, which is modified by filament flexibility. This transition changes filament transport from ballistic to diffusive at long timescales.

Modeling spatiotemporally varying protein-protein interactions in CyLaKS, the Cytoskeleton Lattice-based Kinetic Simulator.

Interaction of cytoskeletal filaments, motor proteins, and crosslinking proteins drives important cellular processes such as cell division and cell movement. Cytoskeletal networks also exhibit nonequilibrium self-assembly in reconstituted systems. An emerging problem in cytoskeletal modeling and simulation is spatiotemporal alteration of the dynamics of filaments, motors, and associated proteins. This can occur due to motor crowding, obstacles along the filament, motor interactions and direction switching, and changes, defects, or heterogeneity in the filament binding lattice. How such spatiotemporally varying cytoskeletal filaments and motor interactions affect their collective properties is not fully understood. We developed the Cytoskeleton Lattice-based Kinetic Simulator (CyLaKS) to investigate such problems. The simulation model builds on previous work by incorporating motor mechanochemistry into a simulation with many interacting motors and/or associated proteins on a discretized lattice. CyLaKS also includes detailed balance in binding kinetics, movement, and lattice heterogeneity. The simulation framework is flexible and extensible for future modeling work and is available on GitHub for others to freely use or build upon. Here we illustrate the use of CyLaKS to study long-range motor interactions, microtubule lattice heterogeneity, motion of a heterodimeric motor, and how changing crosslinker number affects filament separation.

Comparison of explicit and mean-field models of cytoskeletal filaments with crosslinking motors.

In cells, cytoskeletal filament networks are responsible for cell movement, growth, and division. Filaments in the cytoskeleton are driven and organized by crosslinking molecular motors. In reconstituted cytoskeletal systems, motor activity is responsible for far-from-equilibrium phenomena such as active stress, self-organized flow, and spontaneous nematic defect generation. How microscopic interactions between motors and filaments lead to larger-scale dynamics remains incompletely understood. To build from motor-filament interactions to predict bulk behavior of cytoskeletal systems, more computationally efficient techniques for modeling motor-filament interactions are needed. Here, we derive a coarse-graining hierarchy of explicit and continuum models for crosslinking motors that bind to and walk on filament pairs. We compare the steady-state motor distribution and motor-induced filament motion for the different models and analyze their computational cost. All three models agree well in the limit of fast motor binding kinetics. Evolving a truncated moment expansion of motor density speeds the computation by [Formula: see text]-[Formula: see text] compared to the explicit or continuous-density simulations, suggesting an approach for more efficient simulation of large networks. These tools facilitate further study of motor-filament networks on micrometer to millimeter length scales.

Bound-State Diffusion due to Binding to Flexible Polymers in a Selective Biofilter.

Selective biofilters are used by cells to control the transport of proteins, nucleic acids, and other macromolecules. Biological filters demonstrate both high specificity and rapid motion or high flux of proteins. In contrast, high flux comes at the expense of selectivity in many synthetic filters. Binding can lead to selective transport in systems in which the bound particle can diffuse, but the mechanisms that lead to bound diffusion remain unclear. Previous theory has proposed a molecular mechanism of bound-state mobility based only on transient binding to flexible polymers. However, this mechanism has not been directly tested in experiments. We demonstrate that bound mobility via tethered diffusion can be engineered into a synthetic gel using protein fragments derived from the nuclear pore complex. The resulting bound-state diffusion is quantitatively consistent with theory. Our results suggest that synthetic biological filters can be designed to take advantage of tethered diffusion to give rapid, selective transport.

Collective motion of driven semiflexible filaments tuned by soft repulsion and stiffness.

In active matter systems, self-propelled particles can self-organize to undergo collective motion, leading to persistent dynamical behavior out of equilibrium. In cells, cytoskeletal filaments and motor proteins form complex structures important for cell mechanics, motility, and division. Collective dynamics of cytoskeletal systems can be reconstituted using filament gliding experiments, in which cytoskeletal filaments are propelled by surface-bound motor proteins. These experiments have observed diverse dynamical states, including flocks, polar streams, swirling vortices, and single-filament spirals. Recent experiments with microtubules and kinesin motor proteins found that the collective behavior of gliding filaments can be tuned by altering the concentration of the crowding macromolecule methylcellulose in solution. Increasing the methylcellulose concentration reduced filament crossing, promoted alignment, and led to a transition from active, isotropically oriented filaments to locally aligned polar streams. This emergence of collective motion is typically explained as an increase in alignment interactions by Vicsek-type models of active polar particles. However, it is not yet understood how steric interactions and bending stiffness modify the collective behavior of active semiflexible filaments. Here we use simulations of driven filaments with tunable soft repulsion and rigidity in order to better understand how the interplay between filament flexibility and steric effects can lead to different active dynamic states. We find that increasing filament stiffness decreases the probability of filament alignment, yet increases collective motion and long-range order, in contrast to the assumptions of a Vicsek-type model. We identify swirling flocks, polar streams, buckling bands, and spirals, and describe the physics that govern transitions between these states. In addition to repulsion and driving, tuning filament stiffness can promote collective behavior, and controls the transition between active isotropic filaments, locally aligned flocks, and polar streams.

Theory of Cytoskeletal Reorganization during Cross-Linker-Mediated Mitotic Spindle Assembly.

Cells grow, move, and respond to outside stimuli by large-scale cytoskeletal reorganization. A prototypical example of cytoskeletal remodeling is mitotic spindle assembly, during which microtubules nucleate, undergo dynamic instability, bundle, and organize into a bipolar spindle. Key mechanisms of this process include regulated filament polymerization, cross-linking, and motor-protein activity. Remarkably, using passive cross-linkers, fission yeast can assemble a bipolar spindle in the absence of motor proteins. We develop a torque-balance model that describes this reorganization because of dynamic microtubule bundles, spindle-pole bodies, the nuclear envelope, and passive cross-linkers to predict spindle-assembly dynamics. We compare these results to those obtained with kinetic Monte Carlo-Brownian dynamics simulations, which include cross-linker-binding kinetics and other stochastic effects. Our results show that rapid cross-linker reorganization to microtubule overlaps facilitates cross-linker-driven spindle assembly, a testable prediction for future experiments. Combining these two modeling techniques, we illustrate a general method for studying cytoskeletal network reorganization.

Contributions of Microtubule Dynamic Instability and Rotational Diffusion to Kinetochore Capture.

Microtubule dynamic instability allows search and capture of kinetochores during spindle formation, an important process for accurate chromosome segregation during cell division. Recent work has found that microtubule rotational diffusion about minus-end attachment points contributes to kinetochore capture in fission yeast, but the relative contributions of dynamic instability and rotational diffusion are not well understood. We have developed a biophysical model of kinetochore capture in small fission-yeast nuclei using hybrid Brownian dynamics/kinetic Monte Carlo simulation techniques. With this model, we have studied the importance of dynamic instability and microtubule rotational diffusion for kinetochore capture, both to the lateral surface of a microtubule and at or near its end. Over a range of biologically relevant parameters, microtubule rotational diffusion decreased capture time, but made a relatively small contribution compared to dynamic instability. At most, rotational diffusion reduced capture time by 25%. Our results suggest that while microtubule rotational diffusion can speed up kinetochore capture, it is unlikely to be the dominant physical mechanism for typical conditions in fission yeast. In addition, we found that when microtubules undergo dynamic instability, lateral captures predominate even in the absence of rotational diffusion. Counterintuitively, adding rotational diffusion to a dynamic microtubule increases the probability of end-on capture.

Effects of soft interactions and bound mobility on diffusion in crowded environments: a model of sticky and slippery obstacles.

Crowded environments modify the diffusion of macromolecules, generally slowing their movement and inducing transient anomalous subdiffusion. The presence of obstacles also modifies the kinetics and equilibrium behavior of tracers. While previous theoretical studies of particle diffusion have typically assumed either impenetrable obstacles or binding interactions that immobilize the particle, in many cellular contexts bound particles remain mobile. Examples include membrane proteins or lipids with some entry and diffusion within lipid domains and proteins that can enter into membraneless organelles or compartments such as the nucleolus. Using a lattice model, we studied the diffusive movement of tracer particles which bind to soft obstacles, allowing tracers and obstacles to occupy the same lattice site. For sticky obstacles, bound tracer particles are immobile, while for slippery obstacles, bound tracers can hop without penalty to adjacent obstacles. In both models, binding significantly alters tracer motion. The type and degree of motion while bound is a key determinant of the tracer mobility: slippery obstacles can allow nearly unhindered diffusion, even at high obstacle filling fraction. To mimic compartmentalization in a cell, we examined how obstacle size and a range of bound diffusion coefficients affect tracer dynamics. The behavior of the model is similar in two and three spatial dimensions. Our work has implications for protein movement and interactions within cells.

Motor Protein Accumulation on Antiparallel Microtubule Overlaps.

Biopolymers serve as one-dimensional tracks on which motor proteins move to perform their biological roles. Motor protein phenomena have inspired theoretical models of one-dimensional transport, crowding, and jamming. Experiments studying the motion of Xklp1 motors on reconstituted antiparallel microtubule overlaps demonstrated that motors recruited to the overlap walk toward the plus end of individual microtubules and frequently switch between filaments. We study a model of this system that couples the totally asymmetric simple exclusion process for motor motion with switches between antiparallel filaments and binding kinetics. We determine steady-state motor density profiles for fixed-length overlaps using exact and approximate solutions of the continuum differential equations and compare to kinetic Monte Carlo simulations. Overlap motor density profiles and motor trajectories resemble experimental measurements. The phase diagram of the model is similar to the single-filament case for low switching rate, while for high switching rate we find a new (to our knowledge) low density-high density-low density-high density phase. The overlap center region, far from the overlap ends, has a constant motor density as one would naïvely expect. However, rather than following a simple binding equilibrium, the center motor density depends on total overlap length, motor speed, and motor switching rate. The size of the crowded boundary layer near the overlap ends is also dependent on the overlap length and switching rate in addition to the motor speed and bulk concentration. The antiparallel microtubule overlap geometry may offer a previously unrecognized mechanism for biological regulation of protein concentration and consequent activity.

Quantifying yeast microtubules and spindles using the Toolkit for Automated Microtubule Tracking (TAMiT).

Fluorescently labeled proteins absorb and emit light, appearing as Gaussian spots in fluorescence imaging. When fluorescent tags are added to cytoskeletal polymers such as microtubules, a line of fluorescence and even non-linear structures results. While much progress has been made in techniques for imaging and microscopy, image analysis is less well developed. Current analysis of fluorescent microtubules uses either manual tools, such as kymographs, or automated software. As a result, our ability to quantify microtubule dynamics and organization from light microscopy remains limited. Despite development of automated microtubule analysis tools for in vitro studies, analysis of images from cells often depends heavily on manual analysis. One of the main reasons for this disparity is the low signal-to-noise ratio in cells, where background fluorescence is typically higher than in reconstituted systems. Here, we present the Toolkit for Automated Microtubule Tracking (TAMiT), which automatically detects, optimizes and tracks fluorescent microtubules in living yeast cells with sub-pixel accuracy. Using basic information about microtubule organization, TAMiT detects linear and curved polymers using a geometrical scanning technique. Images are fit via an optimization problem for the microtubule image parameters that is solved using non-linear least squares in Matlab. We benchmark our software using simulated images and show that it reliably detects microtubules, even at low signal-to-noise ratios. Then, we use TAMiT to measure monopolar spindle microtubule bundle number, length, and lifetime in a large dataset that includes several S. pombe mutants that affect microtubule dynamics and bundling. The results from the automated analysis are consistent with previous work, and suggest a direct role for CLASP/Cls1 in bundling spindle microtubules. We also illustrate automated tracking of single curved astral microtubules in S. cerevisiae , with measurement of dynamic instability parameters. The results obtained with our fully-automated software are similar to results using hand-tracked measurements. Therefore, TAMiT can facilitate automated analysis of spindle and microtubule dynamics in yeast cells.

Quantifying Yeast Microtubules and Spindles Using the Toolkit for Automated Microtubule Tracking (TAMiT).

Fluorescently labeled proteins absorb and emit light, appearing as Gaussian spots in fluorescence imaging. When fluorescent tags are added to cytoskeletal polymers such as microtubules, a line of fluorescence and even non-linear structures results. While much progress has been made in techniques for imaging and microscopy, image analysis is less well-developed. Current analysis of fluorescent microtubules uses either manual tools, such as kymographs, or automated software. As a result, our ability to quantify microtubule dynamics and organization from light microscopy remains limited. Despite the development of automated microtubule analysis tools for in vitro studies, analysis of images from cells often depends heavily on manual analysis. One of the main reasons for this disparity is the low signal-to-noise ratio in cells, where background fluorescence is typically higher than in reconstituted systems. Here, we present the Toolkit for Automated Microtubule Tracking (TAMiT), which automatically detects, optimizes, and tracks fluorescent microtubules in living yeast cells with sub-pixel accuracy. Using basic information about microtubule organization, TAMiT detects linear and curved polymers using a geometrical scanning technique. Images are fit via an optimization problem for the microtubule image parameters that are solved using non-linear least squares in Matlab. We benchmark our software using simulated images and show that it reliably detects microtubules, even at low signal-to-noise ratios. Then, we use TAMiT to measure monopolar spindle microtubule bundle number, length, and lifetime in a large dataset that includes several S. pombe mutants that affect microtubule dynamics and bundling. The results from the automated analysis are consistent with previous work and suggest a direct role for CLASP/Cls1 in bundling spindle microtubules. We also illustrate automated tracking of single curved astral microtubules in S. cerevisiae, with measurement of dynamic instability parameters. The results obtained with our fully-automated software are similar to results using hand-tracked measurements. Therefore, TAMiT can facilitate automated analysis of spindle and microtubule dynamics in yeast cells.

Agitated ants: regulation and self-organization of incipient nest excavation via collisional cues.

Ants are millimetres in scale yet collectively create metre-scale nests in diverse substrates. To discover principles by which ant collectives self-organize to excavate crowded, narrow tunnels, we studied incipient excavation in small groups of fire ants in quasi-two-dimensional arenas. Excavation rates displayed three stages: initially excavation occurred at a constant rate, followed by a rapid decay, and finally a slower decay scaling in time as t-1/2. We used a cellular automata model to understand such scaling and motivate how rate modulation emerges without global control. In the model, ants estimated their collision frequency with other ants, but otherwise did not communicate. To capture early excavation rates, we introduced the concept of 'agitation'-a tendency of individuals to avoid rest if collisions are frequent. The model reproduced the observed multi-stage excavation dynamics; analysis revealed how parameters affected features of multi-stage progression. Moreover, a scaling argument without ant-ant interactions captures tunnel growth power-law at long times. Our study demonstrates how individual ants may use local collisional cues to achieve functional global self-organization. Such contact-based decisions could be leveraged by other living and non-living collectives to perform tasks in confined and crowded environments.

Bacterial efflux pump modulators prevent bacterial growth in macrophages and under broth conditions that mimic the host environment.

New approaches for combating microbial infections are needed. One strategy for disrupting pathogenesis involves developing compounds that interfere with bacterial virulence. A critical molecular determinant of virulence for Gram-negative bacteria are efflux pumps of the resistance-nodulation-division family, which includes AcrAB-TolC. We previously identified small molecules that bind AcrB, inhibit AcrAB-TolC, and do not appear to damage membranes. These efflux pump modulators (EPMs) were discovered in an in-cell screening platform called SAFIRE (Screen for Anti-infectives using Fluorescence microscopy of IntracellulaR Enterobacteriaceae). SAFIRE identifies compounds that disrupt the growth of a Gram-negative human pathogen, Salmonella enterica serotype Typhimurium (S. Typhimurium), in macrophages. We used medicinal chemistry to iteratively design ~200 EPM35 analogs and test them for activity in SAFIRE, generating compounds with nanomolar potency. Analogs were demonstrated to bind AcrB in a substrate binding pocket by cryo-electron microscopy. Despite having amphipathic structures, the EPM analogs do not disrupt membrane voltage, as monitored by FtsZ localization to the cell septum. The EPM analogs had little effect on bacterial growth in standard Mueller Hinton Broth. However, under broth conditions that mimic the micro-environment of the macrophage phagosome, acrAB is required for growth, the EPM analogs are bacteriostatic, and the EPM analogs increase the potency of antibiotics. These data suggest that under macrophage-like conditions, the EPM analogs prevent the export of a toxic bacterial metabolite(s) through AcrAB-TolC. Thus, compounds that bind AcrB could disrupt infection by specifically interfering with the export of bacterial toxic metabolites, host defense factors, and/or antibiotics.IMPORTANCEBacterial efflux pumps are critical for resistance to antibiotics and for virulence. We previously identified small molecules that inhibit efflux pumps (efflux pump modulators, EPMs) and prevent pathogen replication in host cells. Here, we used medicinal chemistry to increase the activity of the EPMs against pathogens in cells into the nanomolar range. We show by cryo-electron microscopy that these EPMs bind an efflux pump subunit. In broth culture, the EPMs increase the potency (activity), but not the efficacy (maximum effect), of antibiotics. We also found that bacterial exposure to the EPMs appear to enable the accumulation of a toxic metabolite that would otherwise be exported by efflux pumps. Thus, inhibitors of bacterial efflux pumps could interfere with infection not only by potentiating antibiotics, but also by allowing toxic waste products to accumulate within bacteria, providing an explanation for why efflux pumps are needed for virulence in the absence of antibiotics.

Bacterial Efflux Pump Modulators Prevent Bacterial Growth in Macrophages and Under Broth Conditions that Mimic the Host Environment.

New approaches for combatting microbial infections are needed. One strategy for disrupting pathogenesis involves developing compounds that interfere with bacterial virulence. A critical molecular determinant of virulence for Gram-negative bacteria are efflux pumps of the resistance-nodulation-division (RND) family, which includes AcrAB-TolC. We previously identified small molecules that bind AcrB, inhibit AcrAB-TolC, and do not appear to damage membranes. These efflux pump modulators (EPMs) were discovered in an in-cell screening platform called SAFIRE (Screen for Anti-infectives using Fluorescence microscopy of IntracellulaR Enterobacteriaceae). SAFIRE identifies compounds that disrupt the growth of a Gram-negative human pathogen, Salmonella enterica serotype Typhimurium (S. Typhimurium) in macrophages. We used medicinal chemistry to iteratively design ~200 EPM35 analogs and test them for activity in SAFIRE, generating compounds with nanomolar potency. Analogs were demonstrated to bind AcrB in a substrate binding pocket by cryo-electron microscopy (cryo-EM). Despite having amphipathic structures, the EPM analogs do not disrupt membrane voltage, as monitored by FtsZ localization to the cell septum. The EPM analogs had little effect on bacterial growth in standard Mueller Hinton Broth. However, under broth conditions that mimic the micro-environment of the macrophage phagosome, acrAB is required for growth, the EPM analogs are bacteriostatic, and increase the potency of antibiotics. These data suggest that under macrophage-like conditions the EPM analogs prevent the export of a toxic bacterial metabolite(s) through AcrAB-TolC. Thus, compounds that bind AcrB could disrupt infection by specifically interfering with the export of bacterial toxic metabolites, host defense factors, and/or antibiotics.

A new view of how cytoplasmic viscosity affects microtubule dynamics.

How the physical nature of the cytoplasm shapes cellular processes has been increasingly studied in recent years. Work in this issue demonstrates that cytoplasmic viscosity affects the dynamic instability of microtubules. Surprisingly, both polymerization and depolymerization slow with increasing viscosity. This raises interesting questions for future work.

Toward the cellular-scale simulation of motor-driven cytoskeletal assemblies.

The cytoskeleton - a collection of polymeric filaments, molecular motors, and crosslinkers - is a foundational example of active matter, and in the cell assembles into organelles that guide basic biological functions. Simulation of cytoskeletal assemblies is an important tool for modeling cellular processes and understanding their surprising material properties. Here, we present aLENS (a Living Ensemble Simulator), a novel computational framework designed to surmount the limits of conventional simulation methods. We model molecular motors with crosslinking kinetics that adhere to a thermodynamic energy landscape, and integrate the system dynamics while efficiently and stably enforcing hard-body repulsion between filaments. Molecular potentials are entirely avoided in imposing steric constraints. Utilizing parallel computing, we simulate tens to hundreds of thousands of cytoskeletal filaments and crosslinking motors, recapitulating emergent phenomena such as bundle formation and buckling. This simulation framework can help elucidate how motor type, thermal fluctuations, internal stresses, and confinement determine the evolution of cytoskeletal active matter.