Interpretable neural architecture search and transfer learning for understanding CRISPR/Cas9 off-target enzymatic reactions.

Finely-tuned enzymatic pathways control cellular processes, and their dysregulation can lead to disease. Creating predictive and interpretable models for these pathways is challenging because of the complexity of the pathways and of the cellular and genomic contexts. Here we introduce Elektrum, a deep learning framework which addresses these challenges with data-driven and biophysically interpretable models for determining the kinetics of biochemical systems. First, it uses in vitro kinetic assays to rapidly hypothesize an ensemble of high-quality Kinetically Interpretable Neural Networks (KINNs) that predict reaction rates. It then employs a novel transfer learning step, where the KINNs are inserted as intermediary layers into deeper convolutional neural networks, fine-tuning the predictions for reaction-dependent in vivo outcomes. Elektrum makes effective use of the limited, but clean in vitro data and the complex, yet plentiful in vivo data that captures cellular context. We apply Elektrum to predict CRISPR-Cas9 off-target editing probabilities and demonstrate that Elektrum achieves state-of-the-art performance, regularizes neural network architectures, and maintains physical interpretability.

Interpretable neural architecture search and transfer learning for understanding CRISPR-Cas9 off-target enzymatic reactions.

Finely tuned enzymatic pathways control cellular processes, and their dysregulation can lead to disease. Developing predictive and interpretable models for these pathways is challenging because of the complexity of the pathways and of the cellular and genomic contexts. Here we introduce Elektrum, a deep learning framework that addresses these challenges with data-driven and biophysically interpretable models for determining the kinetics of biochemical systems. First, it uses in vitro kinetic assays to rapidly hypothesize an ensemble of high-quality kinetically interpretable neural networks (KINNs) that predict reaction rates. It then employs a transfer learning step, where the KINNs are inserted as intermediary layers into deeper convolutional neural networks, fine-tuning the predictions for reaction-dependent in vivo outcomes. We apply Elektrum to predict CRISPR-Cas9 off-target editing probabilities and demonstrate that Elektrum achieves improved performance, regularizes neural network architectures and maintains physical interpretability.

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.

Active condensation of filaments under spatial confinement.

Living systems exhibit self-organization, a phenomenon that enables organisms to perform functions essential for life. The interior of living cells is a crowded environment in which the self-assembly of cytoskeletal networks is spatially constrained by membranes and organelles. Cytoskeletal filaments undergo active condensation in the presence of crosslinking motor proteins. In past studies, confinement has been shown to alter the morphology of active condensates. Here, we perform simulations to explore systems of filaments and crosslinking motors in a variety of confining geometries. We simulate spatial confinement imposed by hard spherical, cylindrical, and planar boundaries. These systems exhibit non-equilibrium condensation behavior where crosslinking motors condense a fraction of the overall filament population, leading to coexistence of vapor and condensed states. We find that the confinement lengthscale modifies the dynamics and condensate morphology. With end-pausing crosslinking motors, filaments self-organize into half asters and fully-symmetric asters under spherical confinement, polarity-sorted bilayers and bottle-brush-like states under cylindrical confinement, and flattened asters under planar confinement. The number of crosslinking motors controls the size and shape of condensates, with flattened asters becoming hollow and ring-like for larger motor number. End pausing plays a key role affecting condensate morphology: systems with end-pausing motors evolve into aster-like condensates while those with non-end-pausing crosslinking motor proteins evolve into disordered clusters and polarity-sorted bundles.

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.

Mechanisms of chromosome biorientation and bipolar spindle assembly analyzed by computational modeling.

The essential functions required for mitotic spindle assembly and chromosome biorientation and segregation are not fully understood, despite extensive study. To illuminate the combinations of ingredients most important to align and segregate chromosomes and simultaneously assemble a bipolar spindle, we developed a computational model of fission-yeast mitosis. Robust chromosome biorientation requires progressive restriction of attachment geometry, destabilization of misaligned attachments, and attachment force dependence. Large spindle length fluctuations can occur when the kinetochore-microtubule attachment lifetime is long. The primary spindle force generators are kinesin-5 motors and crosslinkers in early mitosis, while interkinetochore stretch becomes important after biorientation. The same mechanisms that contribute to persistent biorientation lead to segregation of chromosomes to the poles after anaphase onset. This model therefore provides a framework to interrogate key requirements for robust chromosome biorientation, spindle length regulation, and force generation in the spindle.

Before a cell divides, it must make a copy of its genetic material and then promptly split in two. This process, called mitosis, is coordinated by many different molecular machines. The DNA is copied, then the duplicated chromosomes line up at the middle of the cell. Next, an apparatus called the mitotic spindle latches onto the chromosomes before pulling them apart. The mitotic spindle is a bundle of long, thin filaments called microtubules. It attaches to chromosomes at the kinetochore, the point where two copied chromosomes are cinched together in their middle. Proper cell division is vital for the healthy growth of all organisms, big and small, and yet some parts of the process remain poorly understood despite extensive study. Specifically, there is more to learn about how the mitotic spindle self-assembles, and how microtubules and kinetochores work together to correctly orient and segregate chromosomes into two sister cells. These nanoscale processes are happening a hundred times a minute, so computer simulations are a good way to test what we know. Edelmaier et al. developed a computer model to simulate cell division in fission yeast, a species of yeast often used to study fundamental processes in the cell. The model simulates how the mitotic spindle assembles, how its microtubules attach to the kinetochore and the force required to pull two sister chromosomes apart. Building the simulation involved modelling interactions between the mitotic spindle and kinetochore, their movement and forces applied. To test its accuracy, model simulations were compared to recordings of the mitotic spindle ­ including its length, structure and position ­ imaged from dividing yeast cells. Running the simulation, Edelmaier et al. found that several key effects are essential for the proper movement of chromosomes in mitosis. This includes holding chromosomes in the correct orientation as the mitotic spindle assembles and controlling the relative position of microtubules as they attach to the kinetochore. Misaligned attachments must also be readily deconstructed and corrected to prevent any errors. The simulations also showed that kinetochores must begin to exert more force (to separate the chromosomes) once the mitotic spindle is attached correctly. Altogether, these findings improve the current understanding of how the mitotic spindle and its counterparts control cell division. Errors in chromosome segregation are associated with birth defects and cancer in humans, and this new simulation could potentially now be used to help make predictions about how to correct mistakes in the process.

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.

Kinesin-5-independent mitotic spindle assembly requires the antiparallel microtubule crosslinker Ase1 in fission yeast.

Bipolar spindle assembly requires a balance of forces where kinesin-5 produces outward pushing forces to antagonize the inward pulling forces from kinesin-14 or dynein. Accordingly, Kinesin-5 inactivation results in force imbalance leading to monopolar spindle and chromosome segregation failure. In fission yeast, force balance is restored when both kinesin-5 Cut7 and kinesin-14 Pkl1 are deleted, restoring spindle bipolarity. Here we show that the cut7Δpkl1Δ spindle is fully competent for chromosome segregation independently of motor activity, except for kinesin-6 Klp9, which is required for anaphase spindle elongation. We demonstrate that cut7Δpkl1Δ spindle bipolarity requires the microtubule antiparallel bundler PRC1/Ase1 to recruit CLASP/Cls1 to stabilize microtubules. Brownian dynamics-kinetic Monte Carlo simulations show that Ase1 and Cls1 activity are sufficient for initial bipolar spindle formation. We conclude that pushing forces generated by microtubule polymerization are sufficient to promote spindle pole separation and the assembly of bipolar spindle in the absence of molecular motors.

Physical determinants of bipolar mitotic spindle assembly and stability in fission yeast.

Mitotic spindles use an elegant bipolar architecture to segregate duplicated chromosomes with high fidelity. Bipolar spindles form from a monopolar initial condition; this is the most fundamental construction problem that the spindle must solve. Microtubules, motors, and cross-linkers are important for bipolarity, but the mechanisms necessary and sufficient for spindle assembly remain unknown. We describe a physical model that exhibits de novo bipolar spindle formation. We began with physical properties of fission-yeast spindle pole body size and microtubule number, kinesin-5 motors, kinesin-14 motors, and passive cross-linkers. Our model results agree quantitatively with our experiments in fission yeast, thereby establishing a minimal system with which to interrogate collective self-assembly. By varying the features of our model, we identify a set of functions essential for the generation and stability of spindle bipolarity. When kinesin-5 motors are present, their bidirectionality is essential, but spindles can form in the presence of passive cross-linkers alone. We also identify characteristic failed states of spindle assembly-the persistent monopole, X spindle, separated asters, and short spindle, which are avoided by the creation and maintenance of antiparallel microtubule overlaps. Our model can guide the identification of new, multifaceted strategies to induce mitotic catastrophes; these would constitute novel strategies for cancer chemotherapy.