NPJ Systems Biology and Applications最新文献

Integrating multi-omics data into predictive models has the potential to enhance accuracy, which is essential for precision medicine. In this study, we developed interpretable predictive models for multi-omics data by employing neural networks informed by prior biological knowledge, referred to as visible networks. These neural networks offer insights into the decision-making process and can unveil novel perspectives on the underlying biological mechanisms associated with traits and complex diseases. We tested the performance, interpretability and generalizability for inferring smoking status, subject age and LDL levels using genome-wide RNA expression and CpG methylation data from the blood of the BIOS consortium (four population cohorts, N_total = 2940). In a cohort-wise cross-validation setting, the consistency of the diagnostic performance and interpretation was assessed. Performance was consistently high for predicting smoking status with an overall mean AUC of 0.95 (95% CI: 0.90-1.00) and interpretation revealed the involvement of well-replicated genes such as AHRR, GPR15 and LRRN3. LDL-level predictions were only generalized in a single cohort with an R² of 0.07 (95% CI: 0.05-0.08). Age was inferred with a mean error of 5.16 (95% CI: 3.97-6.35) years with the genes COL11A2, AFAP1, OTUD7A, PTPRN2, ADARB2 and CD34 consistently predictive. For both regression tasks, we found that using multi-omics networks improved performance, stability and generalizability compared to interpretable single omic networks. We believe that visible neural networks have great potential for multi-omics analysis; they combine multi-omic data elegantly, are interpretable, and generalize well to data from different cohorts.

A complex interplay between various processes underlies the neuropathology of Alzheimer's disease (AD) and its progressive course. Several lines of evidence point to the coupling between Aβ aggregation and neuroinflammation and its role in maintaining brain homeostasis during the long prodromal phase of AD. Little is however known about how this protective mechanism fails and as a result, an irreversible and progressive transition to clinical AD occurs. Here, we introduce a minimal model of a coupled system of Aβ aggregation and inflammation, numerically simulate its dynamical behavior, and analyze its bifurcation properties. The introduced model represents the following events: generation of Aβ monomers, aggregation of Aβ monomers into oligomers and fibrils, induction of inflammation by Aβ aggregates, and clearance of various Aβ species. Crucially, the rates of Aβ generation and clearance are modulated by inflammation level following a Hill-type response function. Despite its relative simplicity, the model exhibits enormously rich dynamics ranging from overdamped kinetics to sustained oscillations. We then specify the region of inflammation- and coupling-related parameters space where a transition to oscillatory dynamics occurs and demonstrate how changes in Aβ aggregation parameters could shift this oscillatory region in parameter space. Our results reveal the propensity of coupled Aβ aggregation-inflammation systems to oscillatory dynamics and propose prolonged sustained oscillations and their consequent immune system exhaustion as a potential mechanism underlying the transition to a more progressive phase of amyloid pathology in AD. The implications of our results in regard to early diagnosis of AD and anti-AD drug development are discussed.

Atrial fibrillation (AF) is the most common form of cardiac arrhythmia, often evolving from paroxysmal episodes to persistent stages over an extended timeframe. While various factors contribute to this progression, the precise biophysical mechanisms driving it remain unclear. Here we explore how rapid firing of cardiomyocytes at the outlet of the pulmonary vein of the left atria can create a substrate for a persistent re-entry wave. This is grounded in a recently formulated mathematical model of the regulation of calcium ion channel density by intracellular calcium concentration. According to the model, the number of calcium channels is controlled by the intracellular calcium concentration. In particular, if the concentration increases above a certain target level, the calcium current is weakened to restore the target level of calcium. During rapid pacing, the intracellular calcium concentration of the cardiomyocytes increases leading to a substantial reduction of the calcium current across the membrane of the myocytes, which again reduces the action potential duration. In a spatially resolved cell-based model of the outlet of the pulmonary vein of the left atria, we show that the reduced action potential duration can lead to re-entry. Initiated by rapid pacing, often stemming from paroxysmal AF episodes lasting several days, the reduction in calcium current is a critical factor. Our findings illustrate how such episodes can foster a conducive environment for persistent AF through electrical remodeling, characterized by diminished calcium currents. This underscores the importance of promptly addressing early AF episodes to prevent their progression to chronic stages.

Biological signal transduction networks are central to information processing and regulation of gene expression across all domains of life. Dysregulation is known to cause a wide array of diseases, including cancers. Here I introduce self-consistent signal transduction analysis, which utilizes genome-scale -omics data (specifically transcriptomics and/or proteomics) in order to predict the flow of information through these networks in an individualized manner. I apply the method to the study of endocrine therapy in breast cancer patients, and show that drugs that inhibit estrogen receptor α elicit a wide array of antitumoral effects, and that their most clinically-impactful ones are through the modulation of proliferative signals that control the genes GREB1, HK1, AKT1, MAPK1, AKT2, and NQO1. This method offers researchers a valuable tool in understanding how and why dysregulation occurs, and how perturbations to the network (such as targeted therapies) effect the network itself, and ultimately patient outcomes.

Energy metabolism is crucial for all living cells, especially during fast growth or stress scenarios. Many cancer and activated immune cells (Warburg effect) or yeasts (Crabtree effect) mostly rely on aerobic glucose fermentation leading to lactate or ethanol, respectively, to generate ATP. In recent years, several mathematical models have been proposed to explain the Warburg effect on theoretical grounds. Besides glucose, glutamine is a very important substrate for eukaryotic cells-not only for biosynthesis, but also for energy metabolism. Here, we present a minimal constraint-based stoichiometric model for explaining both the classical Warburg effect and the experimentally observed respirofermentation of glutamine (WarburQ effect). We consider glucose and glutamine respiration as well as the respective fermentation pathways. Our resource allocation model calculates the ATP production rate, taking into account enzyme masses and, therefore, pathway costs. While our calculation predicts glucose fermentation to be a superior energy-generating pathway in human cells, different enzyme characteristics in yeasts reduce this advantage, in some cases to such an extent that glucose respiration is preferred. The latter is observed for the fungal pathogen Candida albicans, which is a known Crabtree-negative yeast. Further, optimization results show that glutamine is a valuable energy source and important substrate under glucose limitation, in addition to its role as a carbon and nitrogen source of biomass in eukaryotic cells. In conclusion, our model provides insights that glutamine is an underestimated fuel for eukaryotic cells during fast growth and infection scenarios and explains well the observed parallel respirofermentation of glucose and glutamine in several cell types.

Predicting olfactory perceptions from odorant molecules is challenging due to the complex and potentially discontinuous nature of the perceptual space for smells. In this study, we introduce a deep learning model, Mol-PECO (Molecular Representation by Positional Encoding of Coulomb Matrix), designed to predict olfactory perceptions based on molecular structures and electrostatics. Mol-PECO learns the efficient embedding of molecules by utilizing the Coulomb matrix, which encodes atomic coordinates and charges, as an alternative of the adjacency matrix and its Laplacian eigenfunctions as positional encoding of atoms. With a comprehensive dataset of odor molecules and descriptors, Mol-PECO outperforms traditional machine learning methods using molecular fingerprints and graph neural networks based on adjacency matrices. The learned embeddings by Mol-PECO effectively capture the odor space, enabling global clustering of descriptors and local retrieval of similar odorants. This work contributes to a deeper understanding of the olfactory sense and its mechanisms.

Mathematical models of biochemical reaction networks are an important and emerging tool for the study of cell signaling networks involved in disease processes. One promising potential application of such mathematical models is the study of how disease-causing mutations promote the signaling phenotype that contributes to the disease. It is commonly assumed that one must have a thorough characterization of the network readily available for mathematical modeling to be useful, but we hypothesized that mathematical modeling could be useful when there is incomplete knowledge and that it could be a tool for discovery that opens new areas for further exploration. In the present study, we first develop a mechanistic mathematical model of a G-protein coupled receptor signaling network that is mutated in almost all cases of uveal melanoma and use model-driven explorations to uncover and explore multiple new areas for investigating this disease. Modeling the two major, mutually-exclusive, oncogenic mutations (Gα_q/11 and CysLT₂R) revealed the potential for previously unknown qualitative differences between seemingly interchangeable disease-promoting mutations, and our experiments confirmed oncogenic CysLT₂R was impaired at activating the FAK/YAP/TAZ pathway relative to Gα_q/11. This led us to hypothesize that CYSLTR2 mutations in UM must co-occur with other mutations to activate FAK/YAP/TAZ signaling, and our bioinformatic analysis uncovers a role for co-occurring mutations involving the plexin/semaphorin pathway, which has been shown capable of activating this pathway. Overall, this work highlights the power of mechanism-based computational systems biology as a discovery tool that can leverage available information to open new research areas.

There is increasing evidence that biological systems are modular in both structure and function. Complex biological signaling networks such as gene regulatory networks (GRNs) are proving to be composed of subcategories that are interconnected and hierarchically ranked. These networks contain highly dynamic processes that ultimately dictate cellular function over time, as well as influence phenotypic fate transitions. In this work, we use a stochastic multicellular signaling network of pancreatic cancer (PC) to show that the variance in topological rankings of the most phenotypically influential modules implies a strong relationship between structure and function. We further show that induction of mutations alters the modular structure, which analogously influences the aggression and controllability of the disease in silico. We finally present evidence that the impact and location of mutations with respect to PC modular structure directly corresponds to the efficacy of single agent treatments in silico, because topologically deep mutations require deep targets for control.

Immunoglobulins (Ig), which exist either as B-cell receptors (BCR) on the surface of B cells or as antibodies when secreted, play a key role in the recognition and response to antigenic threats. The capability to jointly characterize the BCR and antibody repertoire is crucial for understanding human adaptive immunity. From peripheral blood, bulk BCR sequencing (bulkBCR-seq) currently provides the highest sampling depth, single-cell BCR sequencing (scBCR-seq) allows for paired chain characterization, and antibody peptide sequencing by tandem mass spectrometry (Ab-seq) provides information on the composition of secreted antibodies in the serum. Yet, it has not been benchmarked to what extent the datasets generated by these three technologies overlap and complement each other. To address this question, we isolated peripheral blood B cells from healthy human donors and sequenced BCRs at bulk and single-cell levels, in addition to utilizing publicly available sequencing data. Integrated analysis was performed on these datasets, resolved by replicates and across individuals. Simultaneously, serum antibodies were isolated, digested with multiple proteases, and analyzed with Ab-seq. Systems immunology analysis showed high concordance in repertoire features between bulk and scBCR-seq within individuals, especially when replicates were utilized. In addition, Ab-seq identified clonotype-specific peptides using both bulk and scBCR-seq library references, demonstrating the feasibility of combining scBCR-seq and Ab-seq for reconstructing paired-chain Ig sequences from the serum antibody repertoire. Collectively, our work serves as a proof-of-principle for combining bulk sequencing, single-cell sequencing, and mass spectrometry as complementary methods towards capturing humoral immunity in its entirety.

Cooking, a quintessential creative pursuit, holds profound significance for individuals, communities, and civilizations. Food and cooking transcend mere sensory pleasure to influence nutrition and public health outcomes. Inextricably linked to culinary and cultural heritage, food systems play a pivotal role in sustainability and the survival of life on our planet. Computational Gastronomy is a novel approach for investigating food through a data-driven paradigm. It offers a systematic, rule-based understanding of culinary arts by scrutinizing recipes for taste, nutritional value, health implications, and environmental sustainability. Probing the art of cooking through the lens of computation will open up a new realm of possibilities for culinary creativity. Amidst the ongoing quest for imitating creativity through artificial intelligence, an interesting question would be, 'Can a machine think like a Chef?' Capturing the experience and creativity of a chef in an AI algorithm presents an exciting opportunity for generating a galaxy of hitherto unseen recipes with desirable culinary, flavor, nutrition, health, and carbon footprint profiles.