NPJ Systems Biology and Applications最新文献

Bronchiolitis is the leading cause of infant hospitalization. However, the molecular networks driving bronchiolitis pathobiology remain unknown. Integrative molecular networks, including the transcriptome and metabolome, can identify functional and regulatory pathways contributing to disease severity. Here, we integrated nasopharyngeal transcriptome and metabolome data of 397 infants hospitalized with bronchiolitis in a 17-center prospective cohort study. Using an explainable deep network model, we identified an omics-cluster comprising 401 transcripts and 38 metabolites that distinguishes bronchiolitis severity (test-set AUC, 0.828). This omics-cluster derived a molecular network, where innate immunity-related metabolites (e.g., ceramides) centralized and were characterized by toll-like receptor (TLR) and NF-κB signaling pathways (both FDR < 0.001). The network analyses identified eight modules and 50 existing drug candidates for repurposing, including prostaglandin I₂ analogs (e.g., iloprost), which promote anti-inflammatory effects through TLR signaling. Our approach facilitates not only the identification of molecular networks underlying infant bronchiolitis but the development of pioneering treatment strategies.

Ordinary differential equation (ODE) models are powerful tools for studying the dynamics of metabolic pathways. However, key challenges lie in constructing ODE models for metabolic pathways, specifically in our limited knowledge about which metabolite levels control which reaction rates. Identification of these regulatory networks is further complicated by the limited availability of relevant data. Here, we assess the conditions under which it is feasible to accurately identify regulatory networks in metabolic pathways by computationally fitting candidate network models with biochemical systems theory (BST) kinetics to data of varying quality. We use network motifs commonly found in metabolic pathways as a simplified testbed. Key features correlated with the level of difficulty in identifying the correct regulatory network were identified, highlighting the impact of sampling rate, data noise, and data incompleteness on structural uncertainty. We found that for a simple branched network motif with an equal number of metabolites and fluxes, identification of the correct regulatory network can be largely achieved and is robust to missing one of the metabolite profiles. However, with a bi-substrate bi-product reaction or more fluxes than metabolites in the network motif, the identification becomes more challenging. Stronger regulatory interactions and higher metabolite concentrations were found to be correlated with less structural uncertainty. These results could aid efforts to predict whether the true metabolic regulatory network can be computationally identified for a given stoichiometric network topology and dataset quality, thus helping to identify optimal measures to mitigate such identifiability issues in kinetic model development.

Complex signaling pathways are believed to be responsible for drug resistance. Drug combinations perturbing multiple signaling targets have the potential to reduce drug resistance. The large-scale multi-omic datasets and experimental drug combination synergistic score data are valuable resources to study mechanisms of synergy (MoS) to guide the development of precision drug combinations. However, signaling patterns of MoS are complex and remain unclear, and thus it is challenging to identify synergistic drug combinations in clinical. Herein, we proposed a novel integrative and interpretable graph AI model, DeepSignalingFlow, to uncover the MoS by integrating and mining multi-omic data. The major innovation is that we uncover MoS by modeling the signaling flow from multi-omic features of essential disease proteins to the drug targets, which has not been introduced by the existing models. The model performance was assessed utilizing four distinct drug combination synergy evaluation datasets, i.e., NCI ALMANAC, O'Neil, DrugComb, and DrugCombDB. The comparison results showed that the proposed model outperformed existing graph AI models in terms of synergy score prediction, and can interpret MoS using the core signaling flows. The code is publicly accessible via Github: https://github.com/FuhaiLiAiLab/DeepSignalingFlow.

Glioblastoma is the most common and deadliest brain tumour in adults, with a median survival of 15 months under the current standard of care. Immunotherapies like immune checkpoint inhibitors and oncolytic viruses have been extensively studied to improve this endpoint. However, most thus far have failed. To improve the efficacy of immunotherapies to treat glioblastoma, new single-cell imaging modalities like imaging mass cytometry can be leveraged and integrated with computational models. This enables a better understanding of the tumour microenvironment and its role in treatment success or failure in this hard-to-treat tumour. Here, we implemented an agent-based model that allows for spatial predictions of combination chemotherapy, oncolytic virus, and immune checkpoint inhibitors against glioblastoma. We initialised our model with patient imaging mass cytometry data to predict patient-specific responses and found that oncolytic viruses drive combination treatment responses determined by intratumoral cell density. We found that tumours with higher tumour cell density responded better to treatment. When fixing the number of cancer cells, treatment efficacy was shown to be a function of CD4 + T cell and, to a lesser extent, of macrophage counts. Critically, our simulations show that care must be put into the integration of spatial data and agent-based models to effectively capture intratumoral dynamics. Together, this study emphasizes the use of predictive spatial modelling to better understand cancer immunotherapy treatment dynamics, while highlighting key factors to consider during model design and implementation.

YAP/TAZ signaling pathway is regulated by a multiplicity of feedback loops, crosstalk with other pathways, and both mechanical and biochemical stimuli. Computational modeling serves as a powerful tool to unravel how these different factors can regulate YAP/TAZ, emphasizing biophysical modeling as an indispensable tool for deciphering mechanotransduction and its regulation of cell fate. We provide a critical review of the current state-of-the-art of computational models focused on YAP/TAZ signaling.

Mechanistic mathematical models (MMs) are a powerful tool to help us understand and predict the dynamics of tumour growth under various conditions. In this work, we use 5 MMs with an increasing number of parameters to explore how certain (often overlooked) decisions in estimating parameters from data of experimental tumour growth affect the outcome of the analysis. In particular, we propose a framework for including tumour volume measurements that fall outside the upper and lower limits of detection, which are normally discarded. We demonstrate how excluding censored data results in an overestimation of the initial tumour volume and the MM-predicted tumour volumes prior to the first measurements, and an underestimation of the carrying capacity and the MM-predicted tumour volumes beyond the latest measurable time points. We show in which way the choice of prior for the MM parameters can impact the posterior distributions, and illustrate that reporting the most likely parameters and their 95% credible interval can lead to confusing or misleading interpretations. We hope this work will encourage others to carefully consider choices made in parameter estimation and to adopt the approaches we put forward herein.

We present a study where predictive mechanistic modeling is combined with deep learning methods to predict individual patient survival probabilities under immune checkpoint inhibitor (ICI) immunotherapy. This hybrid approach enables prediction based on both measures that are calculable from mechanistic models of key mechanisms underlying ICI therapy that may not be directly measurable in the clinic and easily measurable quantities or patient characteristics that are not always readily incorporated into predictive mechanistic models. A deep learning time-to-event predictive model trained on a hybrid mechanistic + clinical data set from 93 patients achieved higher per-patient predictive accuracy based on event-time concordance, Brier score, and negative binomial log-likelihood-based criteria than when trained on only mechanistic model-derived values or only clinical data. Feature importance analysis revealed that both clinical and model-derived parameters play prominent roles in increasing prediction accuracy, further supporting the advantage of our hybrid approach.

Network controllability is unifying the traditional control theory with the structural network information rooted in many large-scale biological systems of interest, from intracellular networks in molecular biology to brain neuronal networks. In controllability approaches, the set of minimum driver nodes is not unique, and critical nodes are the most important control elements because they appear in all possible solution sets. On the other hand, a common but largely unexplored feature in network control approaches is the probabilistic failure of edges or the uncertainty in the determination of interactions between molecules. This is particularly true when directed probabilistic interactions are considered. Until now, no efficient algorithm existed to determine critical nodes in probabilistic directed networks. Here we present a probabilistic control model based on a minimum dominating set framework that integrates the probabilistic nature of directed edges between molecules and determines the critical control nodes that drive the entire network functionality. The proposed algorithm, combined with the developed mathematical tools, offers practical efficiency in determining critical control nodes in large probabilistic networks. The method is then applied to the human intracellular signal transduction network revealing that critical control nodes are associated with important biological features and perturbed sets of genes in human diseases, including SARS-CoV-2 target proteins and rare disorders. We believe that the proposed methodology can be useful to investigate multiple biological systems in which directed edges are probabilistic in nature, both in natural systems or when determined with large uncertainties in-silico.

Ligand-receptor systems, covalent modification cycles, and transcriptional networks are the fundamental components of cell signaling and gene expression systems. While their behavior in reaching a steady-state regime under step-like stimulation is well understood, their response under repetitive stimulation, particularly at early time stages is poorly characterized. Yet, early-stage responses to external inputs are arguably as informative as late-stage ones. In simple systems, a periodic stimulation elicits an initial transient response, followed by periodic behavior. Transient responses are relevant when the stimulation has a limited time span, or when the stimulated component's timescale is slow as compared to the timescales of the downstream processes, in which case the latter processes may be capturing only those transients. In this study, we analyze the frequency response of simple motifs at different time stages. We use dose-conserved pulsatile input signals and consider different metrics versus frequency curves. We show that in ligand-receptor systems, there is a frequency preference response in some specific metrics during the transient stages, which is not present in the periodic regime. We suggest this is a general system-level mechanism that cells may use to filter input signals that have consequences for higher order circuits. In addition, we evaluate how the described behavior in isolated motifs is reflected in similar types of responses in cascades and pathways of which they are a part. Our studies suggest that transient frequency preferences are important dynamic features of cell signaling and gene expression systems, which have been overlooked.

The strength of molecular interactions is characterized by their dissociation constants (K_D). Only high-affinity interactions (K_D ≤ 10^-8 M) are extensively investigated and support binary on/off switches. However, such analyses have discounted the presence of low-affinity binders (K_D > 10^-5 M) in the cellular environment. We assess the potential influence of low-affinity binders on high-affinity interactions. By employing Gillespie stochastic simulations and continuous methods, we demonstrate that the presence of low-affinity binders can alter the kinetics and the steady state of high-affinity interactions. We refer to this effect as 'herd regulation' and have evaluated its possible impact in two different contexts including sex determination in Drosophila melanogaster and in signalling systems that employ molecular thresholds. We have also suggested experiments to validate herd regulation in vitro. We speculate that low-affinity binders are prevalent in biological contexts where the outcomes depend on molecular thresholds impacting homoeostatic regulation.