Nathaniel Linden-Santangeli, Jin Zhang, Boris Kramer, Padmini Rangamani
{"title":"Systems modeling and uncertainty quantification of AMP-activated protein kinase signaling.","authors":"Nathaniel Linden-Santangeli, Jin Zhang, Boris Kramer, Padmini Rangamani","doi":"10.1038/s41540-025-00588-w","DOIUrl":"https://doi.org/10.1038/s41540-025-00588-w","url":null,"abstract":"<p><p>AMP-activated protein kinase (AMPK) plays a key role in restoring cellular metabolic homeostasis after energy stress. Importantly, AMPK acts as a hub of metabolic signaling, integrating multiple inputs and acting on numerous downstream targets to activate catabolic processes and inhibit anabolic ones. Despite the importance of AMPK signaling, unlike other well-studied pathways, such as MAPK/ERK or NF-κB, only a handful of mechanistic models of AMPK signaling have been developed. Epistemic uncertainty in the biochemical mechanism of AMPK activation, combined with the complexity of the AMPK pathway, makes model development particularly challenging. Here, we leveraged uncertainty quantification (UQ) methods and recently developed AMPK biosensors to construct a new, data-informed model of AMPK signaling. Specifically, we applied Bayesian parameter estimation and model selection to ensure that model predictions and assumptions are well-constrained to measurements of AMPK activity using recently developed AMPK biosensors. As an application of the new model, we predicted AMPK activity in response to exercise-like stimuli. We found that AMPK acts as a time- and exercise-dependent integrator of its input. Our results highlight how UQ can facilitate model development and address epistemic uncertainty in a complex signaling pathway, such as AMPK. This work shows the potential for future applications of UQ in systems biology to drive new biological insights by incorporating state-of-the-art experimental data at all stages of model development.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"113"},"PeriodicalIF":3.5,"publicationDate":"2025-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145293025","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}
María Moscardó García, Atte Aalto, Arthur N Montanari, Alexander Skupin, Jorge Gonçalves
{"title":"Multi-omic network inference from time-series data.","authors":"María Moscardó García, Atte Aalto, Arthur N Montanari, Alexander Skupin, Jorge Gonçalves","doi":"10.1038/s41540-025-00591-1","DOIUrl":"https://doi.org/10.1038/s41540-025-00591-1","url":null,"abstract":"<p><p>Biological phenotypes emerge from complex interactions across molecular layers. Yet, data-driven approaches to infer these regulatory networks have primarily focused on single-omic studies, overlooking inter-layer regulatory relationships. To address these limitations, we developed MINIE, a computational method that integrates multi-omic data from bulk metabolomics and single-cell transcriptomics through a Bayesian regression approach that explicitly models the timescale separation between molecular layers. We validate the method on both simulated datasets and experimental Parkinson's disease data. MINIE exhibits accurate and robust predictive performance across and within omic layers, including curated multi-omic networks and the lac operon. Benchmarking demonstrated significant improvements over state-of-the-art methods while ranking among the top performers in comprehensive single-cell network inference analysis. The integration of regulatory dynamics across molecular layers and temporal scales provides a powerful tool for comprehensive multi-omic network inference.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"114"},"PeriodicalIF":3.5,"publicationDate":"2025-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145293088","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":"An integrative molecular systems approach unravels mechanisms underlying biphasic nitrate uptake by plant nitrate transporter NRT1.1.","authors":"Seemadri Subhadarshini, Sarthak Sahoo, Mohit Kumar Jolly, Mubasher Rashid","doi":"10.1038/s41540-025-00587-x","DOIUrl":"10.1038/s41540-025-00587-x","url":null,"abstract":"<p><p>Elucidating the mechanisms of transport kinetics in plants is crucial to develop crops that can use nutrients efficiently. The plant nitrate transporter NRT1.1 rapidly switches between high- and low-affinity transport modes to maintain an optimal uptake amidst fluctuations in nitrate levels. This functional switch is regulated by NRT1.1 phosphorylation, but the precise mechanisms remain poorly understood. Here, using an integrated molecular and systems-level modeling, we identify mechanisms underlying biphasic behaviour of NRT1.1. Phosphorylation of NRT1.1 and its binding to nitrate impacts its overall flexibility and synergistically modulates its global conformation, impacting the nitrate transport rate. Integrating these observations with a regulatory network involving kinases CIPK8/CIPK23 and calcium binding proteins CBL1/9, reveals that in high nitrate conditions, CIPK8-mediated sequestration of CBL1 disrupts the CIPK23-CBL complex required for NRT1.1 phosphorylation, switching NRT1.1 to a low-affinity mode. Together, our findings untangle the molecular complexity enabling NRT1.1 phosphorylation switch with broader implications in nitrate sensing and molecular-level adaption to fluctuating external nutrient levels.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"112"},"PeriodicalIF":3.5,"publicationDate":"2025-10-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12518564/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145286680","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}
Xavier Benedicto, Åsmund Flobak, Miguel Ponce-de-Leon, Alfonso Valencia
{"title":"Constraint based modeling of drug induced metabolic changes in a cancer cell line.","authors":"Xavier Benedicto, Åsmund Flobak, Miguel Ponce-de-Leon, Alfonso Valencia","doi":"10.1038/s41540-025-00586-y","DOIUrl":"10.1038/s41540-025-00586-y","url":null,"abstract":"<p><p>Cancer cells frequently reprogramme their metabolism to support growth and survival, making metabolic pathways attractive targets for therapy. In this study, we investigated the metabolic effects of three kinase inhibitors and their synergistic combinations in the gastric cancer cell line AGS using genome-scale metabolic models and transcriptomic profiling. We applied the tasks inferred from the differential expression (TIDE) algorithm to infer pathway activity changes in the different conditions. We also explored a variant of TIDE that uses task-essential genes to infer metabolic task changes, providing a complementary perspective to the original algorithm. Our results revealed widespread down-regulation of biosynthetic pathways, particularly in amino acid and nucleotide metabolism. Combinatorial treatments induced condition-specific metabolic alterations, including strong synergistic effects in the PI3Ki-MEKi condition affecting ornithine and polyamine biosynthesis. These metabolic shifts provide insight into drug synergy mechanisms and highlight potential therapeutic vulnerabilities. To support reproducibility, we developed an open-source Python package, MTEApy, implementing both TIDE frameworks.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"111"},"PeriodicalIF":3.5,"publicationDate":"2025-10-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12501097/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145239183","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}
Hossein Akbarialiabad, Amirmohammad Pasdar, Dédée F Murrell, Mehrnaz Mostafavi, Farhan Shakil, Ehsan Safaee, Sancy A Leachman, Alireza Haghighi, Michelle Tarbox, Chritopher G Bunick, Ayman Grada
{"title":"Enhancing randomized clinical trials with digital twins.","authors":"Hossein Akbarialiabad, Amirmohammad Pasdar, Dédée F Murrell, Mehrnaz Mostafavi, Farhan Shakil, Ehsan Safaee, Sancy A Leachman, Alireza Haghighi, Michelle Tarbox, Chritopher G Bunick, Ayman Grada","doi":"10.1038/s41540-025-00592-0","DOIUrl":"10.1038/s41540-025-00592-0","url":null,"abstract":"<p><p>Digital twins (DTs) can transform randomized clinical trials by improving ethical standards, including safety, informed consent, equity, and data privacy. They also enhance trial efficiency by enabling early detection of adverse events and streamlined design. This paper explores the role of DTs in personalized medicine, from pre-clinical research to post-marketing, while addressing technological, legal, and ethical challenges in implementation.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"110"},"PeriodicalIF":3.5,"publicationDate":"2025-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12494699/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145225536","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}
Nabia Shahreen, Syed Ahsan Shahid, Mahfuze Subhani, Adil Al-Siyabi, Rajib Saha
{"title":"Minimal gene signatures enable high-accuracy prediction of antibiotic resistance in Pseudomonas aeruginosa.","authors":"Nabia Shahreen, Syed Ahsan Shahid, Mahfuze Subhani, Adil Al-Siyabi, Rajib Saha","doi":"10.1038/s41540-025-00584-0","DOIUrl":"10.1038/s41540-025-00584-0","url":null,"abstract":"<p><p>Antimicrobial resistance (AMR) in Pseudomonas aeruginosa poses a critical global health challenge, with current diagnostics relying on slow, culture-based methods. Here, we present a ML framework leveraging transcriptomic data to predict antibiotic resistance with high accuracy. We applied a genetic algorithm to 414 clinical isolates to identify minimal, highly predictive gene sets (~35-40 genes) distinguishing resistant from susceptible strains for meropenem, ciprofloxacin, tobramycin, and ceftazidime. Automated ML classifiers trained on these sets achieved accuracies of 96-99% on test data (F1 scores: 0.93-0.99), surpassing clinical deployment thresholds. Multiple distinct, non-overlapping gene subsets exhibited comparable performance, suggesting that resistance acquisition is associated with changes in the expression of diverse regulatory and metabolic genes. Comparison with known resistance markers from CARD and operon annotations revealed a substantial number of previously unannotated clusters, highlighting significant knowledge gaps in current AMR understanding. Mapping these genes onto independently modulated gene sets (iModulons) revealed transcriptional adaptations across diverse genetic regions. Overall, this study presents a streamlined machine-learning workflow for transcriptomic data and offers a pathway toward rapid diagnostics and personalized treatment strategies against AMR.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"108"},"PeriodicalIF":3.5,"publicationDate":"2025-10-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12491527/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145213323","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}
Simranjit Grewal, Uwa Iyamu, Daniel Ferrer Vinals, Catherine J Mitran, Nidhi Hegde, Stephanie K Yanow
{"title":"Machine learning framework to extract physicochemical features of B-cell epitopes recognized by a cross-reactive antibody.","authors":"Simranjit Grewal, Uwa Iyamu, Daniel Ferrer Vinals, Catherine J Mitran, Nidhi Hegde, Stephanie K Yanow","doi":"10.1038/s41540-025-00583-1","DOIUrl":"10.1038/s41540-025-00583-1","url":null,"abstract":"<p><p>During infection with Plasmodium falciparum in pregnancy, parasites express a unique virulence factor, VAR2CSA, that mediates binding of infected red blood cells to the placenta. A major goal in designing vaccines to protect pregnant women from malaria is to elicit antibodies to VAR2CSA. The challenge is that VAR2CSA is highly polymorphic and identifying conserved epitopes is essential to elicit strain-transcending immunity. Unexpectedly, a mouse monoclonal antibody, 3D10, raised against region II of the unrelated Duffy binding protein from P. vivax (DBPII) cross-reacts with diverse alleles of VAR2CSA in vitro, suggesting that epitopes may be shared across this family of 'Duffy binding-like' (DBL) proteins. Peptide arrays spanning four DBL proteins from two Plasmodium spp, including two alleles of VAR2CSA, DBPII, and PvEBP2 (as a negative control), were screened with 3D10 but the data were too complex to manually identify common epitope sequences. As such, we designed a machine learning framework to analyse the array data. We applied decision trees to extract features correlated to 3D10 binding and evaluated the model on an independent dataset for a rodent Plasmodium DBL protein (PcDBP). Next, we analysed patterns of the features predicted by the model to be strongly associated with 3D10 binding and designed mutant peptides to test complex sequence motifs. Features associated with 3D10 reactivity were mapped onto predicted 3D structures of Plasmodium proteins and validated based on 3D10 reactivity to the recombinant antigens. While the array data identified certain linear epitopes, the framework predicted other epitopes to be conformational. This was demonstrated with PcDBP; as predicted by the model, no linear peptides reacted strongly with 3D10, yet the folded protein was recognized by the antibody in a conformation-dependent manner. With this approach, peptide array data can be mined to extract physicochemical properties of epitopes recognized by cross-reactive antibodies.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"109"},"PeriodicalIF":3.5,"publicationDate":"2025-10-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12491407/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145213274","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}
Aarón Castillo-Jiménez, Adriana Garay-Arroyo, María de La Paz Sánchez, Juan Carlos Martínez-García, Elena R Álvarez-Buylla
{"title":"Cellular patterns in Arabidopsis root epidermis emerge from gene regulatory network and diffusion dynamical feedback.","authors":"Aarón Castillo-Jiménez, Adriana Garay-Arroyo, María de La Paz Sánchez, Juan Carlos Martínez-García, Elena R Álvarez-Buylla","doi":"10.1038/s41540-025-00551-9","DOIUrl":"10.1038/s41540-025-00551-9","url":null,"abstract":"<p><p>We propose a system biology approach to understand how GRNs' dynamical feedback with diffusion of some molecular components underlie the emergence of spatial cellular patterns. We use experimental data on the GRN underlying cell differentiation and spatial arrangement in the root epidermis of WT and mutant Arabidopsis phenotypes to validate our proposal. We test a generalized model of reaction-diffusion, which includes cell-to-cell interaction through lateral inhibition dynamics. The GRN corresponds to the reactive part, and diffusion involves two of its components. The Arabidopsis thaliana root epidermis has a distinct interspersed spatial pattern of hair and non-hair cells. Central to this process is the diffusion of CPC and GL3/EGL3 proteins, which drive lateral inhibition to coordinate cell identity. Existing models have shown a limited predictive power due to incomplete GRN topologies and the lack of explicit diffusion dynamics. Here, we introduce a diffusion-coupled meta-GRN model that integrates positive and negative feedback loops to simulate root epidermal pattern formation in wild-type and mutant lines under varying diffusion levels. By explicitly simulating CPC and GL3/EGL3 protein diffusion, in addition to recovering 28 single and multiple loss-of-function mutant phenotypes, as well as capturing trichoblast and atrichoblast spatial distributions relative to cortex cells, this study presents a 2-D morphospace or phenotypic landscape for epidermis patterning depending on diffusion levels. The findings highlight the critical role of protein diffusion and its dynamic feedback loops with complex GRN in shaping cellular spatial configurations and offer new insights into an extended reaction-diffusion dynamic patterning mechanism that is surely at play in other biological systems.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"107"},"PeriodicalIF":3.5,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12488856/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145207046","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":"On the mechanisms of epidermal stemness and differentiation.","authors":"Saumya Shukla, Raghvendra Singh","doi":"10.1038/s41540-025-00581-3","DOIUrl":"10.1038/s41540-025-00581-3","url":null,"abstract":"<p><p>High Wnt and low Notch activities characterize epidermal stem cells (SCs), while low Wnt and high Notch activities characterize the terminally differentiated epidermal cells (TDCs). However, the mechanism by which transit amplifying cells (TACs) are induced to become terminally differentiated remains unclear. Our analysis suggests that oscillations in Wnt, Notch, and YAP/TAZ activities lead to the production of TDCs from TACs. Furthermore, the role of stem cell markers in epidermal differentiation, regeneration, and the functional aspects of the epidermis remains unclear. Here, based on the ability of the epidermal SCs to induce the differentiation of TACs, we characterize the SCs as having the expression of Notch ligand, Delta, higher than a critical value. Further, we have functionally defined the critical value of the Delta expression by SCs. Our paper may have general implications for the stemness and differentiation of other tissues.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"106"},"PeriodicalIF":3.5,"publicationDate":"2025-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12485183/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145200510","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}
Stéphanie Chevalier, Julia Becker, Yujuan Gui, Vincent Noël, Cui Su, Sascha Jung, Laurence Calzone, Andrei Zinovyev, Antonio Del Sol, Jun Pang, Lasse Sinkkonen, Thomas Sauter, Loïc Paulevé
{"title":"Data-driven inference of Boolean networks from transcriptomes to predict cellular differentiation and reprogramming.","authors":"Stéphanie Chevalier, Julia Becker, Yujuan Gui, Vincent Noël, Cui Su, Sascha Jung, Laurence Calzone, Andrei Zinovyev, Antonio Del Sol, Jun Pang, Lasse Sinkkonen, Thomas Sauter, Loïc Paulevé","doi":"10.1038/s41540-025-00569-z","DOIUrl":"10.1038/s41540-025-00569-z","url":null,"abstract":"<p><p>Boolean networks provide robust, explainable, and predictive models of cellular dynamics, especially for cellular differentiation and fate decision processes. Yet, the construction of such models is extremely challenging, as it requires integrating prior knowledge with experimental observation of the transcriptome, potentially relating thousands of genes. We present a general methodology for integrating transcriptome data and prior knowledge on the underlying gene regulatory network in order to generate automatically ensembles of Boolean networks able to reproduce the modeled qualitative behavior. Our methodology builds on the software BoNesis, which implements the automatic construction of Boolean networks from a specification of their expected structural and dynamical properties. We show how to transform transcriptome data into such a qualitative specification, and then how to exploit the generated ensembles of Boolean networks for identifying families of candidate models, and for predicting robust cellular reprogramming targets. We illustrate the scalability and versatility of our overall approach with two applications: the modeling of hematopoiesis from single-cell RNA-Seq data, and modeling the differentiation of bone marrow stromal cells into adipocytes and osteoblasts from bulk RNA-seq time series data. For this latter case, we took advantage of ensemble modeling to predict combinations of reprogramming factors for trans-differentiation that are robust to model uncertainties due to variations in experimental replicates and choice of binarization method. Moreover, we performed an in silico assessment of the fidelity and efficiency of the reprogramming and conducted preliminary experimental validation.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":"11 1","pages":"105"},"PeriodicalIF":3.5,"publicationDate":"2025-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12475257/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145176837","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}