Cytoplasmic stirring by active carpets.

Large cells often rely on cytoplasmic flows for intracellular transport, maintaining homeostasis, and positioning cellular components. Understanding the mechanisms of these flows is essential for gaining insights into cell function, developmental processes, and evolutionary adaptability. Here, we focus on a class of self-organized cytoplasmic stirring mechanisms that result from fluid-structure interactions between cytoskeletal elements at the cell cortex. Drawing inspiration from streaming flows in late-stage fruit fly oocytes, we propose an analytically tractable active carpet theory. This model deciphers the origins and three-dimensional spatiotemporal organization of such flows. Through a combination of simulations and weakly nonlinear theory, we establish the pathway of the streaming flow to its global attractor: a cell-spanning vortical twister. Our study reveals the inherent symmetries of this emergent flow, its low-dimensional structure, and illustrates how complex fluid-structure interaction aligns with classical solutions in Stokes flow. This framework can be easily adapted to elucidate a broad spectrum of self-organized, cortex-driven intracellular flows.

An adaptive h-refinement method for the boundary element fast multipole method for quasi-static electromagnetic modeling.

Objective.In our recent work pertinent to modeling of brain stimulation and neurophysiological recordings, substantial modeling errors in the computed electric field and potential have sometimes been observed for standard multi-compartment head models. The goal of this study is to quantify those errors and, further, eliminate them through an adaptive mesh refinement (AMR) algorithm. The study concentrates on transcranial magnetic stimulation (TMS), transcranial electrical stimulation (TES), and electroencephalography (EEG) forward problems.Approach.We propose, describe, and systematically investigate an AMR method using the boundary element method with fast multipole acceleration (BEM-FMM) as the base numerical solver. The goal is to efficiently allocate additional unknowns to critical areas of the model, where they will best improve solution accuracy. The implemented AMR method's accuracy improvement is measured on head models constructed from 16 Human Connectome Project subjects under problem classes of TES, TMS, and EEG. Errors are computed between three solutions: an initial non-adaptive solution, a solution found after applying AMR with a conservative refinement rate, and a 'silver-standard' solution found by subsequent 4:1 global refinement of the adaptively-refined model.Main results.Excellent agreement is shown between the adaptively-refined and silver-standard solutions for standard head models. AMR is found to be vital for accurate modeling of TES and EEG forward problems for standard models: an increase of less than 25% (on average) in number of mesh elements for these problems, efficiently allocated by AMR, exposes electric field/potential errors exceeding 60% (on average) in the solution for the unrefined models.Significance.This error has especially important implications for TES dosing prediction-where the stimulation strength plays a central role-and for EEG lead fields. Though the specific form of the AMR method described here is implemented for the BEM-FMM, we expect that AMR is applicable and even required for accurate electromagnetic simulations by other numerical modeling packages as well.

An Adaptive H-Refinement Method for the Boundary Element Fast Multipole Method for Quasi-static Electromagnetic Modeling.

Objective: In our recent work pertinent to modeling of brain stimulation and neurophysiological recordings, substantial modeling errors in the computed electric field and potential have sometimes been observed for standard multi-compartment head models. The goal of this study is to quantify those errors and, further, eliminate them through an adaptive mesh refinement (AMR) algorithm. The study concentrates on transcranial magnetic stimulation (TMS), transcranial electrical stimulation (TES), and electroencephalography (EEG) forward problems. Approach: We propose, describe, and systematically investigate an AMR method using the Boundary Element Method with Fast Multipole Acceleration (BEM-FMM) as the base numerical solver. The goal is to efficiently allocate additional unknowns to critical areas of the model, where they will best improve solution accuracy.The implemented AMR method's accuracy improvement is measured on head models constructed from 16 Human Connectome Project subjects under problem classes of TES, TMS, and EEG. Errors are computed between three solutions: an initial non-adaptive solution, a solution found after applying AMR with a conservative refinement rate, and a "silver-standard" solution found by subsequent 4:1 global refinement of the adaptively-refined model. Main Results: Excellent agreement is shown between the adaptively-refined and silver-standard solutions for standard head models. AMR is found to be vital for accurate modeling of TES and EEG forward problems for standard models: an increase of less than 25% (on average) in number of mesh elements for these problems, efficiently allocated by AMR, exposes electric field/potential errors exceeding 60% (on average) in the solution for the unrefined models. Significance: This error has especially important implications for TES dosing prediction - where the stimulation strength plays a central role - and for EEG lead fields. Though the specific form of the AMR method described here is implemented for the BEM-FMM, we expect that AMR is applicable and even required for accurate electromagnetic simulations by other numerical modeling packages as well.

Quantum electrodynamics of chiral and antichiral waveguide arrays.

We consider the quantum electrodynamics of single photons in arrays of one-way waveguides, each containing many atoms. We investigate both chiral and antichiral arrays, in which the group velocities of the waveguides are the same or alternate in sign, respectively. We find that in the continuum limit, the one-photon amplitude obeys a Dirac equation. In the chiral case, the Dirac equation is hyperbolic, while in the antichiral case it is elliptic. This distinction has implications for the nature of photon transport in waveguide arrays. Our results are illustrated by numerical simulations.

Fast interpolation-based t-SNE for improved visualization of single-cell RNA-seq data.

t-distributed stochastic neighbor embedding (t-SNE) is widely used for visualizing single-cell RNA-sequencing (scRNA-seq) data, but it scales poorly to large datasets. We dramatically accelerate t-SNE, obviating the need for data downsampling, and hence allowing visualization of rare cell populations. Furthermore, we implement a heatmap-style visualization for scRNA-seq based on one-dimensional t-SNE for simultaneously visualizing the expression patterns of thousands of genes. Software is available at https://github.com/KlugerLab/FIt-SNE and https://github.com/KlugerLab/t-SNE-Heatmaps .

A coarse-grained computational model of the nuclear pore complex predicts Phe-Gly nucleoporin dynamics.

The phenylalanine-glycine-repeat nucleoporins (FG-Nups), which occupy the lumen of the nuclear pore complex (NPC), are critical for transport between the nucleus and cytosol. Although NPCs differ in composition across species, they are largely conserved in organization and function. Transport through the pore is on the millisecond timescale. Here, to explore the dynamics of nucleoporins on this timescale, we use coarse-grained computational simulations. These simulations generate predictions that can be experimentally tested to distinguish between proposed mechanisms of transport. Our model reflects the conserved structure of the NPC, in which FG-Nup filaments extend into the lumen and anchor along the interior of the channel. The lengths of the filaments in our model are based on the known characteristics of yeast FG-Nups. The FG-repeat sites also bind to each other, and we vary this association over several orders of magnitude and run 100-ms simulations for each value. The autocorrelation functions of the orientation of the simulated FG-Nups are compared with in vivo anisotropy data. We observe that FG-Nups reptate back and forth through the NPC at timescales commensurate with experimental measurements of the speed of cargo transport through the NPC. Our results are consistent with models of transport where FG-Nup filaments are free to move across the central channel of the NPC, possibly informing how cargo might transverse the NPC.