Publications

From the Golgi apparatus to endosomes, organelles in the endomembrane system exhibit complex and varied morphologies that are often related to their function. Such membrane-bound organelles operate far from equilibrium due to directed fluxes of smaller trafficking vesicles; the physical principles governing the emergence and maintenance of these structures have thus remained elusive. By understanding individual fission and fusion events in terms of active mechano-chemical cycles, we show how such trafficking manifests at the hydrodynamic scale, resulting not only in fluxes of material -- such as membrane area and encapsulated volume -- but also in active stresses that drive momentum transfer between an organelle and its cytosolic environment. Due to the fluid and deformable nature of the bounding membrane, this gives rise to novel physics, coupling nonequilibrium forces to organelle composition, morphology and hydrodynamic flows. We demonstrate how both stable compartment drift and ramified sac-like morphologies, each reminiscent of Golgi-cisternae, emerge naturally from the same underlying nonequilibrium dynamics of fission and fusion.

Machine learning based surrogate models offer researchers powerful tools for accelerating simulation-based workflows. However, as standard datasets in this space often cover small classes of physical behavior, it can be difficult to evaluate the efficacy of new approaches. To address this gap, we introduce the Well: a large-scale collection of datasets containing numerical simulations of a wide variety of spatiotemporal physical systems. The Well draws from domain experts and numerical software developers to provide 15TB of data across 16 datasets covering diverse domains such as biological systems, fluid dynamics, acoustic scattering, as well as magneto-hydrodynamic simulations of extra-galactic fluids or supernova explosions. These datasets can be used individually or as part of a broader benchmark suite. To facilitate usage of the Well, we provide a unified PyTorch interface for training and evaluating models. We demonstrate the function of this library by introducing example baselines that highlight the new challenges posed by the complex dynamics of the Well. The code and data is available at https://github.com/PolymathicAI/the_well.

During eukaryotic cell division, a microtubule-based structure called the spindle exerts forces on chromosomes. The best-studied spindle forces, including those responsible for the separation of sister chromatids, are directed parallel to the spindle’s long axis. By contrast, little is known about forces perpendicular to the spindle axis, which determine the metaphase plate configuration and thus the location of chromosomes in the subsequent nucleus. Using live-cell microscopy, we find that metaphase chromosomes are spatially anti-correlated in mouse oocyte spindles, evidence of previously unknown long-range forces acting perpendicular to the spindle axis. We explain this observation by showing that the spindle’s microtubule network behaves as a nematic liquid crystal and that deformation of the nematic field around embedded chromosomes causes long-range repulsion between them.

Viral infections can induce changes in innate immunity that persist after virus clearance. Here, we used blood samples from a human influenza H3N2 challenge study to perform comprehensive multi-omic analyses. We detected remodeling of immune programs in innate immune cells after resolution of the infection that was proportional in magnitude to the level of prior viral load. We found changes associated with suppressed inflammation including decreased cytokine and AP-1 gene expression as well as decreased accessibility at AP-1 targets and interleukin-related gene promoter regions. We also found decreased histone deacetylase gene expression, increased MAP kinase gene expression, and increased accessibility at interferon-related gene promoter regions. Genes involved in inflammation and epigenetic-remodeling showed modulation of gene-chromatin site regulatory circuit activity. These results reveal a coordinated rewiring of the epigenetic landscape in innate immune cells induced by mild influenza virus infection.

The forces that orient the spindle in human cells remain poorly understood due to a lack of direct mechanical measurements in mammalian systems. We use magnetic tweezers to measure the force on human mitotic spindles. Combining the spindle’s measured resistance to rotation, the speed at which it rotates after laser ablating astral microtubules, and estimates of the number of ablated microtubules reveals that each microtubule contacting the cell cortex is subject to ∼5 pN of pulling force, suggesting that each is pulled on by an individual dynein motor. We find that the concentration of dynein at the cell cortex and extent of dynein clustering are key determinants of the spindle’s resistance to rotation, with little contribution from cytoplasmic viscosity, which we explain using a biophysically based mathematical model. This work reveals how pulling forces on astral microtubules determine the mechanics of spindle orientation and demonstrates the central role of cortical dynein clustering.

Phosphotriesterases (PTEs) represent a class of enzymes capable of efficient neutralization of organophosphates (OPs), a dangerous class of neurotoxic chemicals. PTEs suffer from low catalytic activity, particularly at higher temperatures, due to low thermostability and low solubility. Supercharging, a protein engineering approach via selective mutation of surface residues to charged residues, has been successfully employed to generate proteins with increased solubility and thermostability by promoting charge–charge repulsion between proteins. We set out to overcome the challenges in improving PTE activity against OPs by employing a computational protein supercharging algorithm in Rosetta. Here, we discover two supercharged PTE variants, one negatively supercharged (with −14 net charge) and one positively supercharged (with +12 net charge) and characterize them for their thermodynamic stability and catalytic activity. We find that positively supercharged PTE possesses slight but significant losses in thermostability, which correlates to losses in catalytic efficiency at all temperatures, whereas negatively supercharged PTE possesses increased catalytic activity across 25°C–55°C while offering similar thermostability characteristic to the parent PTE. The impact of supercharging on catalytic efficiency will inform the design of shelf-stable PTE and criteria for enzyme engineering.

This work probes the role of cell geometry in orienting self-organized fluid flows in the late stage Drosophila oocyte. Recent theoretical work has shown that a model, which relies only on hydrodynamic interactions of flexible, cortically anchored microtubules (MTs) and the mechanical loads from molecular motors moving upon them, is sufficient to generate observed flows. While the emergence of flows has been studied in spheres, oocytes change shape during streaming and it was unclear how robust these flows are to the geometry of the cell. Here we use biophysical theory and computational analysis to investigate the role of geometry and find that the axis of rotation is set by the shape of the domain and that the flow is robust to biologically relevant perturbations of the domain shape. Using live imaging and 3D flow reconstruction, we test the predictions of the theory/simulation, finding consistency between the model and live experiments, further demonstrating a geometric dependence on flow direction in late-stage Drosophila oocytes.

We augment the `QUAsi-symmetric Stellarator Repository' (QUASR) to include vacuum field stellarators with quasihelical symmetry using a globalized optimization workflow. The database now has almost 370,000 quasisaxisymmetry and quasihelically symmetric devices along with coil sets, optimized for a variety of aspect ratios, rotational transforms, and discrete rotational symmetries. This paper outlines a couple of ways to explore and characterize the data set. We plot devices on a near-axis quasisymmetry landscape, revealing close correspondence to this predicted landscape. We also use principal component analysis to reduce the dimensionality of the data so that it can easily be visualized in two or three dimensions. Principal component analysis also gives a mechanism to compare the new devices here to previously published ones in the literature. We are able to characterize the structure of the data, observe clusters, and visualize the progression of devices in these clusters. These techniques reveal that the data has structure, and that typically one, two or three principal components are sufficient to characterize it. QUASR is archived at this https URL and can be explored online at this http URL.

Signaling pathways induce stereotyped transcriptional changes as stem cells progress into mature cell types during embryogenesis. Signaling perturbations are necessary to discover which genes are responsive or insensitive to pathway activity. However, gene regulation is additionally dependent on cell state-specific factors like chromatin modifications or transcription factor binding. Thus, transcriptional profiles need to be assayed in single cells to identify potentially multiple, distinct perturbation responses among heterogeneous cell states in an embryo. In perturbation studies, comparing heterogeneous transcriptional states among experimental conditions often requires samples to be collected over multiple independent experiments. Datasets produced in such complex experimental designs can be confounded by batch effects. We present Design-Aware Integration of Single Cell ExpEriments (DAISEE), a new algorithm that models perturbation responses in single-cell datasets with a complex experimental design. We demonstrate that DAISEE improves upon a previously available integrative non-negative matrix factorization framework, more efficiently separating perturbation responses from confounding variation. We use DAISEE to integrate newly collected single-cell RNA-sequencing datasets from 5-hour old zebrafish embryos expressing optimized photoswitchable MEK (psMEK), which globally activates the extracellular signal-regulated kinase (ERK), a signaling molecule involved in many cell specification events. psMEK drives some cells that are normally not exposed to ERK signals towards other wild type states and induces novel states expressing a mixture of transcriptional programs, including precociously activated endothelial genes. ERK signaling is therefore capable of introducing profoundly new gene expression states in developing embryos.

We introduce OpenRAND, a C++17 library aimed at facilitating reproducible scientific research by generating statistically robust yet replicable random numbers in as little as two lines of code, overcoming some of the unnecessary complexities of existing RNG libraries. OpenRAND accommodates single and multi-threaded applications on CPUs and GPUs and offers a simplified, user-friendly API that complies with the C++ standard’s random number engine interface. It is lightweight; provided as a portable, header-only library. It is statistically robust: a suite of built-in tests ensures no pattern exists within single or multiple streams. Despite its simplicity and portability, it remains performant—matching and sometimes outperforming native libraries. Our tests, including a Brownian walk simulation, affirm its reproducibility and ease-of-use while highlight its computational efficiency, outperforming CUDA’s cuRAND by up to 1.8 times.