NPJ Systems Biology and Applications最新文献

The identification of antibody-specific epitopes on virus proteins is crucial for vaccine development and drug design. Nonetheless, traditional wet-lab approaches for the identification of epitopes are both costly and labor-intensive, underscoring the need for the development of efficient and cost-effective computational tools. Here, EpiScan, an attention-based deep learning framework for predicting antibody-specific epitopes, is presented. EpiScan adopts a multi-input and single-output strategy by designing independent blocks for different parts of antibodies, including variable heavy chain (V_H), variable light chain (V_L), complementary determining regions (CDRs), and framework regions (FRs). The block predictions are weighted and integrated for the prediction of potential epitopes. Using multiple experimental data samples, we show that EpiScan, which only uses antibody sequence information, can accurately map epitopes on specific antigen structures. The antibody-specific epitopes on the receptor binding domain (RBD) of SARS coronavirus 2 (SARS-CoV-2) were located by EpiScan, and the potentially valuable vaccine epitope was identified. EpiScan can expedite the epitope mapping process for high-throughput antibody sequencing data, supporting vaccine design and drug development. Availability: For the convenience of related wet-experimental researchers, the source code and web server of EpiScan are publicly available at https://github.com/gzBiomedical/EpiScan.

CRISPR is a precise and effective genome editing technology; but despite several advancements during the last decade, our ability to computationally design gRNAs remains limited. Most predictive models have relatively low predictive power and utilize only the sequence of the target site as input. Here we suggest a new category of features, which incorporate the target site genomic position and the presence of genes close to it. We calculate four features based on gene expression and codon usage bias indices. We show, on CRISPR datasets taken from 3 different cell types, that such features perform comparably with 425 state-of-the-art predictive features, ranking in the top 2-12% of features. We trained new predictive models, showing that adding expression features to them significantly improves their r² by up to 0.04 (relative increase of 39%), achieving average correlations of up to 0.38 on their validation sets; and that these features are deemed important by different feature importance metrics. We believe that incorporating the target site's position, in addition to its sequence, in features such as we have generated here will improve our ability to predict, design and understand CRISPR experiments going forward.

In several carcinomas, including hepatocellular carcinoma, it has been demonstrated that cancer stem cells (CSCs) have enhanced invasiveness and therapy resistance compared to differentiated cancer cells. Mathematical-computational tools could be valuable for integrating experimental results and understanding the phenotypic plasticity mechanisms for CSCs emergence. Based on the literature review, we constructed a Boolean model that recovers eight stable states (attractors) corresponding to the gene expression profile of hepatocytes and mesenchymal cells in senescent, quiescent, proliferative, and stem-like states. The epigenetic landscape associated with the regulatory network was analyzed. We observed that the loss of p53, p16, RB, or the constitutive activation of β-catenin and YAP1 increases the robustness of the proliferative stem-like phenotypes. Additionally, we found that p53 inactivation facilitates the transition of proliferative hepatocytes into stem-like mesenchymal phenotype. Thus, phenotypic plasticity may be altered, and stem-like phenotypes related to CSCs may be easier to attain following the mutation acquisition.

In this review paper we summarize a recent progress on the problem of describing range of dynamics supported by a network. We show that there is natural connection between network models consisting of collections of multivalued monotone boolean functions and ordinary differential equations models. We show how to construct such collections and use them to answer questions about prevalence of cellular phenotypes that correspond to equilibria of network models.

Single-cell omics technologies can measure millions of cells for up to thousands of biomolecular features, enabling data-driven studies of complex biological networks. However, these high-throughput experimental techniques often cannot track individual cells over time, thus complicating the understanding of dynamics such as time trajectories of cell states. These "dynamical phenotypes" are key to understanding biological phenomena such as differentiation fates. We show by mathematical analysis that, in spite of high dimensionality and lack of individual cell traces, three time-points of single-cell omics data are theoretically necessary and sufficient to uniquely determine the network interaction matrix and associated dynamics. Moreover, we show through numerical simulations that an interaction matrix can be accurately determined with three or more time-points even in the presence of sampling and measurement noise typical of single-cell omics. Our results can guide the design of single-cell omics time-course experiments, and provide a tool for data-driven phase-space analysis.

Ovarian cancer is an aggressive, heterogeneous disease, burdened with late diagnosis and resistance to chemotherapy. Clinical features of ovarian cancer could be explained by investigating its metabolism, and how the regulation of specific pathways links to individual phenotypes. Ovarian cancer is of particular interest for metabolic research due to its heterogeneous nature, with five distinct subtypes having been identified, each of which may display a unique metabolic signature. To elucidate metabolic differences, constraint-based modelling (CBM) represents a powerful technology, inviting the integration of 'omics' data, such as transcriptomics. However, many CBM methods have not prioritised accurate growth rate predictions, and there are very few ovarian cancer genome-scale studies. Here, a novel method for CBM has been developed, employing the genome-scale model Human1 and flux balance analysis, enabling the integration of in vitro growth rates, transcriptomics data and media conditions to predict the metabolic behaviour of cells. Using low- and high-grade ovarian cancer, subtype-specific metabolic differences have been predicted, which have been supported by publicly available CRISPR-Cas9 data from the Cancer Cell Line Encyclopaedia and an extensive literature review. Metabolic drivers of aggressive, invasive phenotypes, as well as pathways responsible for increased chemoresistance in low-grade cell lines have been suggested. Experimental gene dependency data has been used to validate areas of the pentose phosphate pathway as essential for low-grade cellular growth, highlighting potential vulnerabilities for this ovarian cancer subtype.

Unraveling how cellular signaling is remodeled upon perturbation is crucial for understanding disease mechanisms and identifying potential drug targets. In this pursuit, computational tools generating mechanistic hypotheses from multi-omics data have invaluable potential. Here, we present a newly implemented version (2.0) of SignalingProfiler, a multi-step pipeline to draw mechanistic hypotheses on the signaling events impacting cellular phenotypes. SignalingProfiler 2.0 derives context-specific signaling networks by integrating proteogenomic data with the prior knowledge-causal network. This is a freely accessible and flexible tool that incorporates statistical, footprint-based, and graph algorithms to accelerate the integration and interpretation of multi-omics data. Through a benchmarking process on three proof-of-concept studies, we demonstrate the tool's ability to generate hierarchical mechanistic networks recapitulating novel and known perturbed signaling and phenotypic outcomes, in both human and mice contexts. In summary, SignalingProfiler 2.0 addresses the emergent need to derive biologically relevant information from complex multi-omics data by extracting interpretable networks.

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.