Aswathy K. Raghu, Indumathi Palanikumar, Karthik Raman
{"title":"Designing function-specific minimal microbiomes from large microbial communities","authors":"Aswathy K. Raghu, Indumathi Palanikumar, Karthik Raman","doi":"10.1038/s41540-024-00373-1","DOIUrl":"https://doi.org/10.1038/s41540-024-00373-1","url":null,"abstract":"<p>Microorganisms exist in large communities of diverse species, exhibiting various functionalities. The mammalian gut microbiome, for instance, has the functionality of digesting dietary fibre and producing different short-chain fatty acids. Not all microbes present in a community contribute to a given functionality; it is possible to find a <i>minimal</i> microbiome, which is a subset of the large microbiome, that is capable of performing the functionality while maintaining other community properties such as growth rate and metabolite production. Such a minimal microbiome will also contain keystone species for SCFA production in that community. In this work, we present a systematic constraint-based approach to identify a minimal microbiome from a large community for a user-proposed function. We employ a top-down approach with sequential deletion followed by solving a mixed-integer linear programming problem with the objective of minimising the <i>L</i><sub>1</sub>-norm of the membership vector. Notably, we consider quantitative measures of community growth rate and metabolite production rates. We demonstrate the utility of our algorithm by identifying the minimal microbiomes corresponding to three model communities of the gut, and discuss their validity based on the presence of the keystone species in the community. Our approach is generic, flexible and finds application in studying a variety of microbial communities. The algorithm is available from https://github.com/RamanLab/minMicrobiome.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140829447","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}
Xiulan Lai, Xiaopei Jiao, Haojian Zhang, Jinzhi Lei
{"title":"Computational modeling reveals key factors driving treatment-free remission in chronic myeloid leukemia patients","authors":"Xiulan Lai, Xiaopei Jiao, Haojian Zhang, Jinzhi Lei","doi":"10.1038/s41540-024-00370-4","DOIUrl":"https://doi.org/10.1038/s41540-024-00370-4","url":null,"abstract":"<p>Patients with chronic myeloid leukemia (CML) who receive tyrosine kinase inhibitors (TKIs) have been known to achieve treatment-free remission (TFR) upon discontinuing treatment. However, the underlying mechanisms of this phenomenon remain incompletely understood. This study aims to elucidate the mechanism of TFR in CML patients, focusing on the feedback interaction between leukemia stem cells and the bone marrow microenvironment. We have developed a mathematical model to explore the interplay between leukemia stem cells and the bone marrow microenvironment, allowing for the simulation of CML progression dynamics. Our proposed model reveals a dichotomous response following TKI discontinuation, with two distinct patient groups emerging: one prone to early molecular relapse and the other capable of achieving long-term TFR after treatment cessation. This finding aligns with clinical observations and underscores the essential role of feedback interaction between leukemic cells and the tumor microenvironment in sustaining TFR. Notably, we have shown that the ratio of leukemia cells in peripheral blood (PBLC) and the tumor microenvironment (TME) index can be a valuable predictive tool for identifying patients likely to achieve TFR after discontinuing treatment. This study provides fresh insights into the mechanism of TFR in CML patients and underscores the significance of microenvironmental control in achieving TFR.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140808827","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}
Renee Ti Chou, Amed Ouattara, Matthew Adams, Andrea A. Berry, Shannon Takala-Harrison, Michael P. Cummings
{"title":"Positive-unlabeled learning identifies vaccine candidate antigens in the malaria parasite Plasmodium falciparum","authors":"Renee Ti Chou, Amed Ouattara, Matthew Adams, Andrea A. Berry, Shannon Takala-Harrison, Michael P. Cummings","doi":"10.1038/s41540-024-00365-1","DOIUrl":"https://doi.org/10.1038/s41540-024-00365-1","url":null,"abstract":"<p>Malaria vaccine development is hampered by extensive antigenic variation and complex life stages of <i>Plasmodium</i> species. Vaccine development has focused on a small number of antigens, many of which were identified without utilizing systematic genome-level approaches. In this study, we implement a machine learning-based reverse vaccinology approach to predict potential new malaria vaccine candidate antigens. We assemble and analyze <i>P. falciparum</i> proteomic, structural, functional, immunological, genomic, and transcriptomic data, and use positive-unlabeled learning to predict potential antigens based on the properties of known antigens and remaining proteins. We prioritize candidate antigens based on model performance on reference antigens with different genetic diversity and quantify the protein properties that contribute most to identifying top candidates. Candidate antigens are characterized by gene essentiality, gene ontology, and gene expression in different life stages to inform future vaccine development. This approach provides a framework for identifying and prioritizing candidate vaccine antigens for a broad range of pathogens.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140808821","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 Allee-based distributed algorithm for microbial whole-cell sensors","authors":"Fabricio Cravo, Matthias Függer, Thomas Nowak","doi":"10.1038/s41540-024-00363-3","DOIUrl":"https://doi.org/10.1038/s41540-024-00363-3","url":null,"abstract":"<p>Reliable detection of substances present at potentially low concentrations is a problem common to many biomedical applications. Complementary to well-established enzyme-, antibody-antigen-, and sequencing-based approaches, so-called microbial whole-cell sensors, i.e., synthetically engineered microbial cells that sense and report substances, have been proposed as alternatives. Typically these cells operate independently: a cell reports an analyte upon local detection.</p><p>In this work, we analyze a distributed algorithm for microbial whole-cell sensors, where cells communicate to coordinate if an analyte has been detected. The algorithm, inspired by the Allee effect in biological populations, causes cells to alternate between a logical 0 and 1 state in response to reacting with the particle of interest. When the cells in the logical 1 state exceed a threshold, the algorithm converts the remaining cells to the logical 1 state, representing an easily-detectable output signal. We validate the algorithm through mathematical analysis and simulations, demonstrating that it works correctly even in noisy cellular environments.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140803850","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}
Patrick C. Kinnunen, Brock A. Humphries, Gary D. Luker, Kathryn E. Luker, Jennifer J. Linderman
{"title":"Characterizing heterogeneous single-cell dose responses computationally and experimentally using threshold inhibition surfaces and dose-titration assays","authors":"Patrick C. Kinnunen, Brock A. Humphries, Gary D. Luker, Kathryn E. Luker, Jennifer J. Linderman","doi":"10.1038/s41540-024-00369-x","DOIUrl":"https://doi.org/10.1038/s41540-024-00369-x","url":null,"abstract":"<p>Single cancer cells within a tumor exhibit variable levels of resistance to drugs, ultimately leading to treatment failures. While tumor heterogeneity is recognized as a major obstacle to cancer therapy, standard dose-response measurements for the potency of targeted kinase inhibitors aggregate populations of cells, obscuring intercellular variations in responses. In this work, we develop an analytical and experimental framework to quantify and model dose responses of individual cancer cells to drugs. We first explore the connection between population and single-cell dose responses using a computational model, revealing that multiple heterogeneous populations can yield nearly identical population dose responses. We demonstrate that a single-cell analysis method, which we term a threshold inhibition surface, can differentiate among these populations. To demonstrate the applicability of this method, we develop a dose-titration assay to measure dose responses in single cells. We apply this assay to breast cancer cells responding to phosphatidylinositol-3-kinase inhibition (PI3Ki), using clinically relevant PI3Kis on breast cancer cell lines expressing fluorescent biosensors for kinase activity. We demonstrate that MCF-7 breast cancer cells exhibit heterogeneous dose responses with some cells requiring over ten-fold higher concentrations than the population average to achieve inhibition. Our work reimagines dose-response relationships for cancer drugs in an emerging paradigm of single-cell tumor heterogeneity.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140616976","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}
Emil Andersson, Ellen V. Rothenberg, Carsten Peterson, Victor Olariu
{"title":"T-cell commitment inheritance—an agent-based multi-scale model","authors":"Emil Andersson, Ellen V. Rothenberg, Carsten Peterson, Victor Olariu","doi":"10.1038/s41540-024-00368-y","DOIUrl":"https://doi.org/10.1038/s41540-024-00368-y","url":null,"abstract":"<p>T-cell development provides an excellent model system for studying lineage commitment from a multipotent progenitor. The intrathymic development process has been thoroughly studied. The molecular circuitry controlling it has been dissected and the necessary steps like programmed shut off of progenitor genes and T-cell genes upregulation have been revealed. However, the exact timing between decision-making and commitment stage remains unexplored. To this end, we implemented an agent-based multi-scale model to investigate inheritance in early T-cell development. Treating each cell as an agent provides a powerful tool as it tracks each individual cell of a simulated T-cell colony, enabling the construction of lineage trees. Based on the lineage trees, we introduce the concept of the last common ancestors (LCA) of committed cells and analyse their relations, both at single-cell level and population level. In addition to simulating wild-type development, we also conduct knockdown analysis. Our simulations predicted that the commitment is a three-step process that occurs on average over several cell generations once a cell is first prepared by a transcriptional switch. This is followed by the loss of the Bcl11b-opposing function approximately two to three generations later. This is when our LCA analysis indicates that the decision to commit is taken even though in general another one to two generations elapse before the cell actually becomes committed by transitioning to the DN2b state. Our results showed that there is decision inheritance in the commitment mechanism.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140616854","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}
Julia Corell-Sierra, Joan Marquez-Molins, María-Carmen Marqués, Andrea Gabriela Hernandez-Azurdia, Roser Montagud-Martínez, María Cebriá-Mendoza, José M. Cuevas, Eliseo Albert, David Navarro, Guillermo Rodrigo, Gustavo Gómez
{"title":"SARS-CoV-2 remodels the landscape of small non-coding RNAs with infection time and symptom severity","authors":"Julia Corell-Sierra, Joan Marquez-Molins, María-Carmen Marqués, Andrea Gabriela Hernandez-Azurdia, Roser Montagud-Martínez, María Cebriá-Mendoza, José M. Cuevas, Eliseo Albert, David Navarro, Guillermo Rodrigo, Gustavo Gómez","doi":"10.1038/s41540-024-00367-z","DOIUrl":"https://doi.org/10.1038/s41540-024-00367-z","url":null,"abstract":"<p>The COVID-19 pandemic caused by the coronavirus SARS-CoV-2 has significantly impacted global health, stressing the necessity of basic understanding of the host response to this viral infection. In this study, we investigated how SARS-CoV-2 remodels the landscape of small non-coding RNAs (sncRNA) from a large collection of nasopharyngeal swab samples taken at various time points from patients with distinct symptom severity. High-throughput RNA sequencing analysis revealed a global alteration of the sncRNA landscape, with abundance peaks related to species of 21-23 and 32-33 nucleotides. Host-derived sncRNAs, including microRNAs (miRNAs), transfer RNA-derived small RNAs (tsRNAs), and small nucleolar RNA-derived small RNAs (sdRNAs) exhibited significant differential expression in infected patients compared to controls. Importantly, miRNA expression was predominantly down-regulated in response to SARS-CoV-2 infection, especially in patients with severe symptoms. Furthermore, we identified specific tsRNAs derived from Glu- and Gly-tRNAs as major altered elements upon infection, with 5’ tRNA halves being the most abundant species and suggesting their potential as biomarkers for viral presence and disease severity prediction. Additionally, down-regulation of C/D-box sdRNAs and altered expression of tinyRNAs (tyRNAs) were observed in infected patients. These findings provide valuable insights into the host sncRNA response to SARS-CoV-2 infection and may contribute to the development of further diagnostic and therapeutic strategies in the clinic.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140616899","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":"Radiopharmaceutical transport in solid tumors via a 3-dimensional image-based spatiotemporal model","authors":"Anahita Piranfar, Farshad Moradi Kashkooli, Wenbo Zhan, Ajay Bhandari, Babak Saboury, Arman Rahmim, M. Soltani","doi":"10.1038/s41540-024-00362-4","DOIUrl":"https://doi.org/10.1038/s41540-024-00362-4","url":null,"abstract":"<p>Lutetium-177 prostate-specific membrane antigen (<sup>177</sup>Lu-PSMA)-targeted radiopharmaceutical therapy is a clinically approved treatment for patients with metastatic castration-resistant prostate cancer (mCRPC). Even though common practice reluctantly follows “one size fits all” approach, medical community believes there is significant room for deeper understanding and personalization of radiopharmaceutical therapies. To pursue this aim, we present a 3-dimensional spatiotemporal radiopharmaceutical delivery model based on clinical imaging data to simulate pharmacokinetic of <sup>177</sup>Lu-PSMA within the prostate tumors. The model includes interstitial flow, radiopharmaceutical transport in tissues, receptor cycles, association/dissociation with ligands, synthesis of PSMA receptors, receptor recycling, internalization of radiopharmaceuticals, and degradation of receptors and drugs. The model was studied for a range of values for injection amount (100–1000 nmol), receptor density (10–500 nmol•l<sup>–1</sup>), and recycling rate of receptors (10<sup>–4</sup> to 10<sup>–1</sup> min<sup>–1</sup>). Furthermore, injection type, different convection-diffusion-reaction mechanisms, characteristic time scales, and length scales are discussed. The study found that increasing receptor density, ligand amount, and labeled ligands improved radiopharmaceutical uptake in the tumor. A high receptor recycling rate (0.1 min<sup>–1</sup>) increased radiopharmaceutical concentration by promoting repeated binding to tumor cell receptors. Continuous infusion results in higher radiopharmaceutical concentrations within tumors compared to bolus administration. These insights are crucial for advancing targeted therapy for prostate cancer by understanding the mechanism of radiopharmaceutical distribution in tumors. Furthermore, measures of characteristic length and advection time scale were computed. The presented spatiotemporal tumor transport model can analyze different physiological parameters affecting <sup>177</sup>Lu-PSMA delivery.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140594336","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}
Ataur Katebi, Xiaowen Chen, Daniel Ramirez, Sheng Li, Mingyang Lu
{"title":"Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML","authors":"Ataur Katebi, Xiaowen Chen, Daniel Ramirez, Sheng Li, Mingyang Lu","doi":"10.1038/s41540-024-00366-0","DOIUrl":"https://doi.org/10.1038/s41540-024-00366-0","url":null,"abstract":"<p>Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140594482","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}
Elyas Mohammadi, Samira Dashti, Neda Shafizade, Han Jin, Cheng Zhang, Simon Lam, Mojtaba Tahmoorespur, Adil Mardinoglu, Mohammad Hadi Sekhavati
{"title":"Drug repositioning for immunotherapy in breast cancer using single-cell analysis","authors":"Elyas Mohammadi, Samira Dashti, Neda Shafizade, Han Jin, Cheng Zhang, Simon Lam, Mojtaba Tahmoorespur, Adil Mardinoglu, Mohammad Hadi Sekhavati","doi":"10.1038/s41540-024-00359-z","DOIUrl":"https://doi.org/10.1038/s41540-024-00359-z","url":null,"abstract":"<p>Immunomodulatory peptides, while exhibiting potential antimicrobial, antifungal, and/or antiviral properties, can play a role in stimulating or suppressing the immune system, especially in pathological conditions like breast cancer (BC). Thus, deregulation of these peptides may serve as an immunotherapeutic strategy to enhance the immune response. In this meta-analysis, we utilized single-cell RNA sequencing data and known therapeutic peptides to investigate the deregulation of these peptides in malignant versus normal human breast epithelial cells. We corroborated our findings at the chromatin level using ATAC-seq. Additionally, we assessed the protein levels in various BC cell lines. Moreover, our in-house drug repositioning approach was employed to identify potential drugs that could positively impact the relapse-free survival of BC patients. Considering significantly deregulated therapeutic peptides and their role in BC pathology, our approach aims to downregulate B2M and SLPI, while upregulating PIGR, DEFB1, LTF, CLU, S100A7, and SCGB2A1 in BC epithelial cells through our drug repositioning pipeline. Leveraging the LINCS L1000 database, we propose BRD-A06641369 for B2M downregulation and ST-4070043 and BRD-K97926541 for SLPI downregulation without negatively affecting the MHC complex as a significantly correlated pathway with these two genes. Furthermore, we have compiled a comprehensive list of drugs for the upregulation of other selected immunomodulatory peptides. Employing an immunotherapeutic approach by integrating our drug repositioning pipeline with single-cell analysis, we proposed potential drugs and drug targets to fortify the immune system against BC.</p>","PeriodicalId":19345,"journal":{"name":"NPJ Systems Biology and Applications","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140594483","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}