CyberSec Research

Background: In vitro endothelial cell culture is widely used to study angiogenesis. Histomicrographic images of cell networks are often analyzed manually, a process that is time-consuming and subjective. Automated tools like ImageJ (NIH) can assist, but are often slow and inaccurate. Additionally, as endothelial networks grow more complex, traditional architectural metrics may not fully reflect network maturity. To address these limitations, we developed tubuleTracker, a software tool that quantifies endothelial network architecture and maturity rapidly and objectively. Methods: Human umbilical vein endothelial cells were cultured in an extracellular matrix, and 54 images were acquired using phase contrast microscopy. Each image was analyzed manually by three independent reviewers, and by both ImageJ and tubuleTracker. Key metrics included tubule count, total length, node count, tubule area, and vessel circularity. In parallel, trained scientists rated each image for angiogenesis maturity on a 1-5 scale (1 = most mature). Results: Analysis time per image differed significantly: manual (8 min), ImageJ (58+/-4 s), and tubuleTracker (6+/-2 s) (p<0.0001). Significant differences were also found in tubule count (manual 168+/-SD, tubuleTracker 92+/-SD, ImageJ 433+/-SD), length, and node count (all p<0.0001). tubuleTracker's metrics varied significantly across angiogenesis maturity scores, including tubule count, length, node count, area, and circularity (all p<0.0001). Conclusions: tubuleTracker was faster and more consistent than both manual and ImageJ-based analysis. Vessel circularity proved especially effective in capturing angiogenesis maturity. tubuleTracker is available as free shareware for the biomedical research community.

We study how simple eukaryotic organisms make decisions in response to competing stimuli in the context of phototaxis by the unicellular alga $Chlamydomonas~reinhardtii$. While negatively phototactic cells swim directly away from a collimated light beam, when presented with two beams of adjustable intersection angle and intensities, we find that cells swim in a direction given by an intensity-weighted average of the two light propagation vectors. This geometrical law is a fixed point of an adaptive model of phototaxis and minimizes the average light intensity falling on the anterior pole of the cell. At large angular separations, subpopulations of cells swim away from one source or the other, or along the direction of the geometrical law, with some cells stochastically switching between the three directions. This behavior is shown to arise from a population-level distribution of photoreceptor locations that breaks front-back symmetry of photoreception.

Biological image analysis has traditionally focused on measuring specific visual properties of interest for cells or other entities. A complementary paradigm gaining increasing traction is image-based profiling - quantifying many distinct visual features to form comprehensive profiles which may reveal hidden patterns in cellular states, drug responses, and disease mechanisms. While current tools like CellProfiler can generate these feature sets, they pose significant barriers to automated and reproducible analyses, hindering machine learning workflows. Here we introduce cp_measure, a Python library that extracts CellProfiler's core measurement capabilities into a modular, API-first tool designed for programmatic feature extraction. We demonstrate that cp_measure features retain high fidelity with CellProfiler features while enabling seamless integration with the scientific Python ecosystem. Through applications to 3D astrocyte imaging and spatial transcriptomics, we showcase how cp_measure enables reproducible, automated image-based profiling pipelines that scale effectively for machine learning applications in computational biology.

Recent measurements of Norway spruce have revealed stress-state-dependent normalized creep behavior, highlighting a gap in our fundamental understanding. This study examines whether the anisotropic response originates from the micro-structural, cellular nature of composite cell walls with varying tracheid types. Cell wall creep parameters are identified via surrogate-based inverse parameter identification, applied to hierarchical micro-mechanical and FEM models of increasing topological complexity up to the growth ring scale. Despite microstructural disorder, simulated creep curves converge toward a universal set of proportionality factors. The results indicate that directional creep behavior cannot be attributed solely to tissue-scale topology, and that realistic predictions require the inclusion of non-linear material responses at stress concentration sites.

This study focuses on the synthesis and characterization of advanced polymeric composite electrospun nanofibers (NFs) containing magnetic oxide nanoparticles (NPs). By leveraging the method of electrospinning, the research aims to investigate polymer composites with enhanced interfacial properties, improved double-layer capacitance, and adequate biocompatibility. Electrospun polyacrylonitrile (PAN) NFs embedded with Fe2O3 and MnZn ferrite NPs were comprehensively characterized using advanced techniques, i.e., Fourier transform infrared spectroscopy (FTIR), X-ray photoelectron spectroscopy (XPS), high-resolution scanning electron microscopy (HR-SEM), X-ray diffraction (XRD), and alternating gradient field magnetometry (AGFM). The incorporation of metal oxide NPs led to significant changes in the thermal, spectroscopic, and morphological properties of the NFs. XPS analysis confirmed increased oxidation, graphitic carbon content, and the formation of new nitrogen functionalities after heat treatment. Furthermore, interactions between nitrile groups and metal ions were observed, indicating the influence of nanoparticles on surface chemistry. Magnetic characterization demonstrated the potential of these composite NFs to generate magnetic fields for biomedical manipulation. Cytocompatibility studies revealed no significant impact on the viability or morphology of human mesenchymal stromal cells, highlighting their biocompatibility. These findings suggest the promising use of PAN-magnetic NFs in applications including targeted drug administration, magnetic resonance imaging (MRI), and magnetic hyperthermia for cancer treatment.

The objective of this paper is to investigate the structural stability, dynamic properties, and potential interactions among Amyloid Precursor Protein (APP), Tau, and Alpha-synuclein through a series of molecular dynamics simulations that integrate publicly available structural data, detailed force-field parameters, and comprehensive analytical protocols. By focusing on these three proteins, which are each implicated in various neurodegenerative disorders, the study aims to elucidate how their conformational changes and interprotein contact sites may influence larger biological processes. Through rigorous evaluation of their folding behaviors, energetic interactions, and residue-specific functions, this work contributes to the broader understanding of protein aggregation mechanisms and offers insights that may ultimately guide therapeutic intervention strategies.

The advent of single-cell multi-omics technologies has enabled the simultaneous profiling of diverse omics layers within individual cells. Integrating such multimodal data provides unprecedented insights into cellular identity, regulatory processes, and disease mechanisms. However, it remains challenging, as current methods often rely on selecting highly variable genes or peaks during preprocessing, which may inadvertently discard crucial biological information. Here, we present scMamba, a foundation model designed to integrate single-cell multi-omics data without the need for prior feature selection while preserving genomic positional information. scMamba introduces a patch-based cell tokenization strategy that treats genomics regions as words (tokens) and cells as sentences. Building upon the concept of state space duality, scMamba distills rich biological insights from high-dimensional, sparse single-cell multi-omics data. Additionally, our novel contrastive learning approach, enhanced with cosine similarity regularization, enables superior alignment across omics layers compared to traditional methods. Systematic benchmarking across multiple datasets demonstrates that scMamba significantly outperforms state-of-the-art methods in preserving biological variation, aligning omics layers, and enhancing key downstream tasks such as clustering, cell type annotation, and trajectory inference. Our findings position scMamba as a powerful tool for large-scale single-cell multi-omics integration, capable of handling large-scale atlases and advancing biological discovery.

Tip growth in filamentous cells, such as root hairs, moss protonemata, and fungal hyphae, depends on coordinated cell wall extension driven by turgor pressure, wall mechanics, and exocytosis. We introduce a dual-configuration model that incorporates both turgid and unturgid states to describe cell wall growth as the combined effect of elastic deformation and irreversible extension. This framework infers exocytosis profiles directly from cell morphology and elastic stretches, formulated as an initial value problem based on the self-similarity condition. Applying the model to Medicago truncatula root hairs, moss Physcomitrium patens protonemata, and hyphoid-like shapes, we find that exocytosis peaks at the tip in tapered cells but shifts to an annular region away from the apex in flatter-tip cells beyond a threshold. The model generalizes previous fluid models and provides a mechanistic link between exocytosis distribution and cell shape, explaining observed variations in tip-growing cells across species.

Mutual information is a theoretically grounded metric for quantifying cellular signaling pathways. However, its measurement demands characterization of both input and output distributions, limiting practical applications. Here, we present alternative method that alleviates this requirement using dual reporter systems. By extending extrinsic-intrinsic noise analysis, we derive a mutual information estimator that eliminates the need to measure input distribution. We demonstrate our method by analyzing the bacterial chemotactic pathway, regarding multiple flagellar motors as natural dual reporters. We show the biological relevance of the measured information flow by comparing it with theoretical bounds on sensory information. This framework opens new possibilities for quantifying information flow in cellular signaling pathways.

Generating high-fidelity and biologically plausible synthetic single-cell RNA sequencing (scRNA-seq) data, especially with conditional control, is challenging due to its high dimensionality, sparsity, and complex biological variations. Existing generative models often struggle to capture these unique characteristics and ensure robustness to structural noise in cellular networks. We introduce LapDDPM, a novel conditional Graph Diffusion Probabilistic Model for robust and high-fidelity scRNA-seq generation. LapDDPM uniquely integrates graph-based representations with a score-based diffusion model, enhanced by a novel spectral adversarial perturbation mechanism on graph edge weights. Our contributions are threefold: we leverage Laplacian Positional Encodings (LPEs) to enrich the latent space with crucial cellular relationship information; we develop a conditional score-based diffusion model for effective learning and generation from complex scRNA-seq distributions; and we employ a unique spectral adversarial training scheme on graph edge weights, boosting robustness against structural variations. Extensive experiments on diverse scRNA-seq datasets demonstrate LapDDPM's superior performance, achieving high fidelity and generating biologically-plausible, cell-type-specific samples. LapDDPM sets a new benchmark for conditional scRNA-seq data generation, offering a robust tool for various downstream biological applications.

We present a coarse-grained stochastic model for axonal extension on periodic arrays of parallel micropatterns that integrates three key biophysical mechanisms: (i) the molecular clutch that couples actin retrograde flow to substrate adhesions, (ii) an active biopolymer-based mechanism generating traction-force fluctuations, and (iii) the mechanical interaction of the growth cone with the micropatterned substrate. Using the Shannon-Jaynes maximum entropy principle with constraints derived from experimental observations, we derive a unique probability distribution for the colored acceleration noise that enters the Langevin equation. The resulting stationary process exhibits power-law temporal correlations with negative exponent, which accounts for the observed superdiffusive dynamics of axons. For biologically relevant parameters the model predicts this exponent to be -1/2, in close quantitative agreement with measurements of cortical neurons cultured on patterned substrates.

Single-cell transcriptomics has become a great source for data-driven insights into biology, enabling the use of advanced deep learning methods to understand cellular heterogeneity and transcriptional regulation at the single-cell level. With the advent of spatial transcriptomics data we have the promise of learning about cells within a tissue context as it provides both spatial coordinates and transcriptomic readouts. However, existing models either ignore spatial resolution or the gene regulatory information. Gene regulation in cells can change depending on microenvironmental cues from neighboring cells, but existing models neglect gene regulatory patterns with hierarchical dependencies across levels of abstraction. In order to create contextualized representations of cells and genes from spatial transcriptomics data, we introduce HEIST, a hierarchical graph transformer-based foundation model for spatial transcriptomics and proteomics data. HEIST models tissue as spatial cellular neighborhood graphs, and each cell is, in turn, modeled as a gene regulatory network graph. The framework includes a hierarchical graph transformer that performs cross-level message passing and message passing within levels. HEIST is pre-trained on 22.3M cells from 124 tissues across 15 organs using spatially-aware contrastive learning and masked auto-encoding objectives. Unsupervised analysis of HEIST representations of cells, shows that it effectively encodes the microenvironmental influences in cell embeddings, enabling the discovery of spatially-informed subpopulations that prior models fail to differentiate. Further, HEIST achieves state-of-the-art results on four downstream task such as clinical outcome prediction, cell type annotation, gene imputation, and spatially-informed cell clustering across multiple technologies, highlighting the importance of hierarchical modeling and GRN-based representations.

Intervertebral discs are avascular and maintain immune privilege. However, during intervertebral disc degeneration (IDD), this barrier is disrupted, leading to extensive immune cell infiltration and localized inflammation. In degenerated discs, macrophages, T lymphocytes, neutrophils, and granulocytic myeloid-derived suppressor cells (G-MDSCs) are key players, exhibiting functional heterogeneity. Dysregulated activation of inflammatory pathways, including nuclear factor kappa-B (NF-kappaB), interleukin-17 (IL-17), and nucleotide-binding oligomerization domain-like receptor protein 3 (NLRP3) inflammasome activation, drives local pro-inflammatory responses, leading to cell apoptosis and extracellular matrix (ECM) degradation. Innovative immunotherapies, including exosome-based treatments, CRISPR/Cas9-mediated gene editing, and chemokine-loaded hydrogel systems, have shown promise in reshaping the immunological niche of intervertebral discs. These strategies can modulate dysregulated immune responses and create a supportive environment for tissue regeneration. However, current studies have not fully elucidated the mechanisms of inflammatory memory and the immunometabolic axis, and they face challenges in balancing tissue regeneration with immune homeostasis. Future studies should employ interdisciplinary approaches such as single-cell and spatial transcriptomics to map a comprehensive immune atlas of IDD, elucidate intercellular crosstalk and signaling networks, and develop integrated therapies combining targeted immunomodulation with regenerative engineering, thereby facilitating the clinical translation of effective IDD treatments.

Understanding the interactions between cells and the extracellular matrix (ECM) during collective cell invasion is crucial for advancements in tissue engineering, cancer therapies, and regenerative medicine. This study focuses on the roles of contact guidance and ECM remodelling in directing cell behaviour, with a particular emphasis on exploring how differences in cell phenotype impact collective cell invasion. We present a computationally tractable two-dimensional hybrid model of collective cell migration within the ECM, where cells are modelled as individual entities and collagen fibres as a continuous tensorial field. Our model incorporates random motility, contact guidance, cell-cell adhesion, volume filling, and the dynamic remodelling of collagen fibres through cellular secretion and degradation. Through a comprehensive parameter sweep, we provide valuable insights into how differences in the cell phenotype, in terms of the ability of the cell to migrate, secrete, degrade, and respond to contact guidance cues from the ECM, impacts the characteristics of collective cell invasion.

We study single cell E.coli chemotaxis in a spatio-temporally varying attractant environment. Modeling the attractant concentration in the form of a traveling sine wave, we measure in our simulations, the chemotactic drift velocity of the cell for different propagation speed of the attractant wave. We find a highly non-trivial dependence where the chemotactic drift velocity changes sign, and also shows multiple peaks. For slowly moving attractant wave, drift velocity is negative, i.e. the drift motion is directed opposite to wave propagation. As the wave speed increases, drift velocity shows a negative peak, then changes sign, reaches a positive peak and finally becomes zero when the wave moves too fast for the cell to respond. We explain this rich behavior from the difference in attractant gradient perceived by the cell during its run along the propagation direction and opposite to it. In particular, when the cell moves in the same direction as the wave, the relative velocity of the cell with respect to the wave becomes zero when the wave speed matches the run speed. In this limit, the cell is able to ride the wave and experiences no concentration gradient during these runs. On the contrary, for runs in the opposite direction, no such effect is present and the effective gradient increases monotonically with the wave speed. We show, using detailed quantitative measurements, how this difference gives rise to the counter-intuitive behavior of chemotactic drift velocity described above.

Multicellular organisms develop a wide variety of highly-specialized cell types. The consistency and robustness of developmental cell fate trajectories suggests that complex gene regulatory networks effectively act as low-dimensional cell fate landscapes. A complementary set of works draws on the theory of dynamical systems to argue that cell fate transitions can be categorized into universal decision-making classes. However, the theory connecting geometric landscapes and decision-making classes to high-dimensional gene expression space is still in its infancy. Here, we introduce a phenomenological model that allows us to identify gene expression signatures of decision-making classes from single-cell RNA-sequencing time-series data. Our model combines low-dimensional gradient-like dynamical systems and high-dimensional Hopfield networks to capture the interplay between cell fate, gene expression, and signaling pathways. We apply our model to the maturation of alveolar cells in mouse lungs to show that the transient appearance of a mixed alveolar type 1/type 2 state suggests the triple cusp decision-making class. We also analyze lineage-tracing data on hematopoetic differentiation and show that bipotent neutrophil-monocyte progenitors likely undergo a heteroclinic flip bifurcation. Our results suggest it is possible to identify universal decision-making classes for cell fate transitions directly from data.

Many biological processes can be thought of as the result of an underlying dynamics in which the system repeatedly undergoes distinct and abortive trajectories with the dynamical process only ending when some specific process, purpose, structure or function is achieved. A classic example is the way in which microtubules attach to kinetochores as a prerequisite for chromosome segregation and cell division. In this example, the dynamics is characterized by apparently futile time histories in which microtubules repeatedly grow and shrink without chromosomal attachment. We hypothesize that for biological processes for which it is not the initial conditions that matter, but rather the final state, this kind of exploratory dynamics is biology's unique and necessary solution to achieving these functions with high fidelity. This kind of cause and effect relationship can be contrasted to examples from physics and chemistry where the initial conditions determine the outcome. In this paper, we examine the similarities of many biological processes that depend upon random trajectories starting from the initial state and the selection of subsets of these trajectories to achieve some desired functional final state. We begin by reviewing the long history of the principles of dynamics, first in the context of physics, and then in the context of the study of life. These ideas are then stacked against the broad categories of biological phenomenology that exhibit exploratory dynamics. We then build on earlier work by making a quantitative examination of a succession of increasingly sophisticated models for exploratory dynamics, all of which share the common feature of being a series of repeated trials that ultimately end in a "winning" trajectory. We also explore the ways in which microscopic parameters can be tuned to alter exploratory dynamics as well as the energetic burden of performing such processes.

The emerging field of immunometabolism has underscored the central role of metabolic pathways in orchestrating immune cell function. Far from being passive background processes, metabolic activities actively regulate key immune responses. Fundamental pathways such as glycolysis, the tricarboxylic acid (TCA) cycle, and oxidative phosphorylation critically shape the behavior of immune cells, influencing macrophage polarization, T cell activation, and dendritic cell function. In this review, we synthesize recent advances in immunometabolism, with a focus on the metabolic mechanisms that govern the responses of both innate and adaptive immune cells to bacterial, viral, and fungal pathogens. Drawing on experimental, computational, and integrative methodologies, we highlight how metabolic reprogramming contributes to host defense in response to infection. These findings reveal new opportunities for therapeutic intervention, suggesting that modulation of metabolic pathways could enhance immune function and improve pathogen clearance.