Yujing Liu, Adriana Zanca, Michael P H Stumpf, Lucy Ham
{"title":"Molecular noise modulates transitions in the cell-fate differentiation landscape.","authors":"Yujing Liu, Adriana Zanca, Michael P H Stumpf, Lucy Ham","doi":"10.1038/s41540-026-00762-8","DOIUrl":"10.1038/s41540-026-00762-8","url":null,"abstract":"<p><p>Waddington's epigenetic landscape has become one of biology's cornerstone metaphors, widely used both conceptually and computationally. Cell types are often associated with the qualitatively stable points or valleys of this landscape. In previous work, we showed that the molecular noise dominating sub-cellular dynamics can distort and profoundly reshape this landscape. In non-equilibrium systems, an equally profound question arises: to what extent does noise alter the transition paths between valleys in such dynamic landscapes? We tackle this question using a set of illustrative exemplars and show that noise gives rise to paths that differ substantially from the canonical least-action paths calculated under deterministic dynamics. We dissect the dynamics of these exemplars and determine the reactive density and transition currents, which show us, respectively, where and how transitions occur for different realizations of stochastic dynamics. Our analysis unambiguously demonstrates that reaction paths for stochastic dynamics diverge non-trivially from their deterministic least-action paths or simple barrier crossing models.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13646251/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148295856","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sidrah Maryam, Martin Manolov, Rafael Kramann, Sikander Hayat
{"title":"Pathway-centric visualization of cell-cell communication in single-cell transcriptomics data.","authors":"Sidrah Maryam, Martin Manolov, Rafael Kramann, Sikander Hayat","doi":"10.1038/s41540-026-00768-2","DOIUrl":"10.1038/s41540-026-00768-2","url":null,"abstract":"<p><p>Single-nuclei transcriptomics enables investigating ligand-receptor mediated cell-cell communication between different cell-types. However, current tools do not allow for non-programmatic means to access this analysis. Additionally, most methods can not account for molecular pathways involved in downstream receptor signaling in cell-cell communication. We developed scVizComm, a ShinyApp based web-portal that can be used to interactively visualize ligand-receptor networks across cell-types and different conditions. scVizComm can be used to host pathway centric pre-calculated cell-cell communication results. Using human kidney and heart single-nuclei transcriptomics data, we showcase the utility of scVizComm in understanding ligand-receptor interactions involved in fibrosis related pathways. Additionally using mouse anterior brain data, we demonstrate the major ligand-receptor expression exploration in space. scVizComm allows interactive visualization and pathway centric prioritization of cell-cell communication analyses.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"12 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13283212/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148295784","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Growth rate-dependent mutations and phenotypic switching facilitate rapid adaptation.","authors":"Anton Grishechkin, Omer Karin","doi":"10.1038/s41540-026-00758-4","DOIUrl":"10.1038/s41540-026-00758-4","url":null,"abstract":"<p><p>The mutation rates of many species are regulated in response to internal and external conditions. Studies in several model organisms have shown that mutation rates can be modulated by the organism's growth rate, particularly through stress-induced mutagenesis. Here, we investigate the functional consequences of this coupling. Using a model of asexual evolution on a high-dimensional, single-peaked fitness landscape, we show that regulating mutation rate through growth rate provides an efficient strategy for rapid adaptation. This mechanism maximizes growth by balancing mutation and selection pressures. We approximate the optimal coupling as a biphasic response function that prescribes a low, constant mutation rate at high growth rates, which gradually increases once growth falls below a critical threshold. This adaptive strategy depends on two key parameters: mutational gain and critical growth rate, both determined by environmental conditions. This raises the question of how organisms establish these parameters in fluctuating environments. We demonstrate that weak selection in well-adapted populations, combined with long-term phenotypic fluctuations, can generate variability in these parameters. Such standing variation prepares populations for rapid evolution following environmental change. Our findings connect this adaptive strategy to specific molecular regulatory circuits and provide predictions for sequential antibiotic treatment regimens.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13635105/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148265322","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Intelligent tool orchestration for rapid mechanistic model prototyping: MCP servers as AI-biology interfaces.","authors":"Marco Ruscone, Miguel Vazquez, Alfonso Valencia","doi":"10.1038/s41540-026-00767-3","DOIUrl":"10.1038/s41540-026-00767-3","url":null,"abstract":"<p><p>Constructing multicellular mechanistic models traditionally requires extensive time and computational expertise. We introduce intelligent tool orchestration via Model Context Protocol (MCP) servers, enabling Large Language Model (LLM) agents to act as AI laboratory assistants for rapid model prototyping. We demonstrate this approach by constructing a multiscale model of cancer cell fate in response to TNF using an AI agent connected to MCP servers interfacing with three complementary tools: NeKo for gene regulatory networks construction, MaBoSS for Boolean models simulation, and PhysiCell for setting up multicellular agent-based models. This workflow was executed entirely through natural language interactions, without manual coding, direct parameter editing, or manual modification of generated model files. Through this use case, we identified key principles for biological AI-tool integration, specifically regarding tool granularity, session management, and flexible orchestration. Testing across multiple LLMs demonstrated our framework's portability, though model-dependent variations emphasize the need for rigorous validation. Ultimately, this work establishes a foundation for AI-assisted rapid prototyping, enabling researchers to explore computational hypotheses more rapidly through natural language interaction.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13631280/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148258963","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Seol Ah Park, Davide Gurnari, Paweł Dłotko, Yongsoo Kim, Jooyoung Hahn
{"title":"Spatial analysis of malignant-immune cell interactions in the tumor microenvironment using topological data analysis.","authors":"Seol Ah Park, Davide Gurnari, Paweł Dłotko, Yongsoo Kim, Jooyoung Hahn","doi":"10.1038/s41540-026-00763-7","DOIUrl":"10.1038/s41540-026-00763-7","url":null,"abstract":"<p><p>The spatial interactions between malignant and immune cells in the tumor microenvironment are important for tumor immunobiology and patient outcomes. However, analytical tools that can extract rigorous yet interpretable spatial features and link them to patient outcomes remain limited. We propose a framework integrating TDA with statistical approaches to extract interpretable spatial features characterizing malignant-immune interactions. We introduce Topological Malignant Region (TopMR), which uses topological persistence to automatically define regions of malignant cells, providing an objective reference for spatial analysis even when tumor boundaries are ambiguous. Global-scale infiltration is quantified using signed distances from immune cells to the TopMR boundary and local-scale interactions are captured via malignant cell density around individual immune cells. These global and local features are integrated into a unified signed distance-density ( <math><mrow><mi>s</mi> <mi>DD</mi></mrow> </math> ) space, enabling comprehensive characterization of spatial patterns. We apply this framework to high-resolution multiplex immunofluorescence images of diffuse large B-cell lymphoma, analyzing both malignant-enriched and tumor border regions. Two-stage hierarchical clustering stratifies patients based on spatial interaction patterns, revealing associations with survival outcomes. This framework provides an end-to-end pipeline from spatial feature extraction to clinical interpretation, suggesting how region-aware spatial analysis can capture biologically meaningful patterns linked to patient survival.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13620163/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148240045","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Determining parameter impact in systems biology models via sensitivity analysis: a comparative approach.","authors":"Kelsey I Gasior","doi":"10.1038/s41540-026-00757-5","DOIUrl":"10.1038/s41540-026-00757-5","url":null,"abstract":"<p><p>Mathematical modeling is a powerful tool to understand biological phenomena, predict behavior, and guide experiments. Model parameters represent experimentally derived rates, and it can be challenging to understand how parameters affect the system's output(s). Global sensitivity analysis (GSA) can help determine how uncertainty in model outcomes can be attributed to parameters. Methods like Sobol' indices provide detailed analysis but at a high computational cost. Thus, it may be necessary to precede it by another method to reduce the parameter space. Morris Method Screening (MMS) is a qualitative analysis that can provide information about sensitivity and parameter interactions, making it a good partner for Sobol'. Partial Rank Correlation Coefficient (PRCC) may also be a viable option as it has a similar cost to MMS. This work tests whether a ranking (PRCC) or screening (MMS) method is better when used in combination with Sobol' analysis. Using a model of the epithelial-mesenchymal transition as a test case, this work shows that PRCC-Sobol' and MMS-Sobol' can produce diverging results, making it difficult to analyze how biological events are informed by specific rates and quantities. However, the differences between these methods mean that these techniques may be most effective when all three are used together.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13623827/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148252536","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A powerful representation learning method for enhanced analysis of incomplete multi-omics data.","authors":"Jenna L Ballard, Zongyu Dai, Li Shen, Qi Long","doi":"10.1038/s41540-026-00759-3","DOIUrl":"10.1038/s41540-026-00759-3","url":null,"abstract":"<p><p>Integrative analysis of multi-omics data provides a more comprehensive and nuanced view of a subject's biological state. However, high-dimensionality and ubiquitous modality missingness present significant analytical challenges. Existing methods for incomplete multi-omics data are scarce, do not fully leverage both modality-specific and shared information, and produce task-biased representations. We propose JASMINE, a self-supervised representation learning method for incomplete multi-omics data that preserves both modality-specific and joint information and enhances sample similarity structure. JASMINE produces embeddings that achieve superior performance across multiple tasks for two different incomplete multi-omics datasets while requiring only a single round of training per dataset.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13612587/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148227010","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Gemma van der Voort, Maike Effern, Michelle C R Yong, Lukas Kiwitz, Roberta Turiello, Sonia Leonardelli, Susanna Ng, Dillon Corvino, Tobias Bald, Nicole Glodde, Kevin Thurley, Jan Hasenauer, Michael Hölzel
{"title":"Longitudinal spatial neutrophil profiling during ACT in murine melanoma reveals distinct lymph node infiltration patterns.","authors":"Gemma van der Voort, Maike Effern, Michelle C R Yong, Lukas Kiwitz, Roberta Turiello, Sonia Leonardelli, Susanna Ng, Dillon Corvino, Tobias Bald, Nicole Glodde, Kevin Thurley, Jan Hasenauer, Michael Hölzel","doi":"10.1038/s41540-026-00765-5","DOIUrl":"10.1038/s41540-026-00765-5","url":null,"abstract":"<p><p>Reactive neutrophil infiltration can restrain CD8<sup>+</sup> T cell expansion in lymph nodes during adoptive T cell therapy (ACT), yet its spatiotemporal regulation remains incompletely understood. Levaraging flow cytometry and multiplex immunofluorescence data, we performed a time-resolved quantitative assessment of immune cell dynamics in tumor-draining lymph node (tdLN) and non-tumor-draining lymph node (non-tdLN) in a melanoma mouse model receiving ACT. Transferred tumor-reactive CD8<sup>+</sup> T cells accumulated and expanded early after treatment initiation, showing the highest frequency of a favorable central memory 13 CD8<sup>+</sup> T cell phenotype in the tdLN. Enhancing innate immune signaling in melanomas increased neutrophil influx into lymph nodes, particularly the non-tdLN; however, within the tdLN, neutrophils were enriched in the T cell zone, which also contained the largest absolute reservoir of transferred CD8<sup>+</sup> T cells. Together, these findings indicate that tdLN and non-tdLN differ in early neutrophil dynamics and compartmentalization during ACT, influenced by the strength of innate immune signaling in the tumor.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"12 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13261029/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148227016","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Reshmi J S Patel, Ayesha Das, Aliya Anil, Gauri Patel, Ryan T Woodall, Chase Christenson, Tarini Thiagarajan, Ernesto A B F Lima, Guillermo Lorenzo, David A Hormuth, Russell C Rockne, Vikram Adhikarla, Chengyue Wu, Thomas E Yankeelov
{"title":"Integrating imaging and mathematical modeling to predict and optimize patient outcomes in oncology.","authors":"Reshmi J S Patel, Ayesha Das, Aliya Anil, Gauri Patel, Ryan T Woodall, Chase Christenson, Tarini Thiagarajan, Ernesto A B F Lima, Guillermo Lorenzo, David A Hormuth, Russell C Rockne, Vikram Adhikarla, Chengyue Wu, Thomas E Yankeelov","doi":"10.1038/s41540-026-00764-6","DOIUrl":"https://doi.org/10.1038/s41540-026-00764-6","url":null,"abstract":"<p><p>Biology-based mathematical models calibrated to patient-specific imaging data can accurately predict and optimize therapeutic response for individual cancer patients. We present the framework, motivation, and examples of using imaging data to calibrate mathematical models and predict patient outcomes. We then discuss patient-specific digital twins for optimizing interventions and outcomes. We aim to demonstrate the importance of imaging-based mathematical models in continuing the paradigm shift from population-based to individual-based cancer care.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":3.5,"publicationDate":"2026-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148218673","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Hyperpathway: visualizing organization of pathway-molecule enriched interactions in omics studies via hyperbolic bipartite network embedding.","authors":"Ilyes Abdelhamid, Ziheng Liao, Yuchi Liu, Armel Lefebvre, Aldo Acevedo, Carlo Vittorio Cannistraci","doi":"10.1038/s41540-026-00752-w","DOIUrl":"10.1038/s41540-026-00752-w","url":null,"abstract":"<p><p>Pathway enrichment analysis (PEA) of omics data identifies significant pathway-molecule associations, yet delivers results as tabular lists in which complex systems-biology insights remain inaccessible. Hyperpathway is an open-access network-based webtool that addresses this limitation through three original innovations: (1) conversion of a PEA results table into a pathway-molecule bipartite network; (2) a minimal artificial linking strategy to resolve structural disconnections; (3) a leaf removal and post-hoc reinsertion pipeline that accelerates coalescent embedding without any loss of geometric fidelity. The resulting network is visualized in a two-dimensional hyperbolic disk with flexible coloring schemes encoding hierarchical relevance, connectivity similarity, statistical significance, or user-defined annotations; revealing latent functional modules that are invisible in conventional tabular outputs. Validated on genomic, metabolomic, and lipidomic datasets, Hyperpathway enables a deeper, systems-level understanding of the interplay between pathways and their molecular components, providing insights that go beyond p-value-based significance testing. Beyond PEA, Hyperpathway can be used as a general-purpose open webtool for fast hyperbolic embedding and interactive visualization of any bipartite network.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":" ","pages":""},"PeriodicalIF":4.4,"publicationDate":"2026-06-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13586325/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148205570","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}