Flexible modeling of regulatory networks improves transcription factor activity estimation

Transcriptional regulation plays a crucial role in determining cell fate and disease, yet inferring the key regulators from gene expression data remains a significant challenge. Existing methods for estimating transcription factor (TF) activity often rely on static TF-gene interaction databases and cannot adapt to changes in regulatory mechanisms across different cell types and disease conditions. Here, we present a new algorithm - Transcriptional Inference using Gene Expression and Regulatory data (TIGER) - that overcomes these limitations by flexibly modeling activation and inhibition events, up-weighting essential edges, shrinking irrelevant edges towards zero through a sparse Bayesian prior, and simultaneously estimating both TF activity levels and changes in the underlying regulatory network. When applied to yeast and cancer TF knock-out datasets, TIGER outperforms comparable methods in terms of prediction accuracy. Moreover, our application of TIGER to tissue- and cell-type-specific RNA-seq data demonstrates its ability to uncover differences in regulatory mechanisms. Collectively, our findings highlight the utility of modeling context-specific regulation when inferring transcription factor activities.

Identification of Molecular Mechanisms of Ameloblastoma and Drug Repositioning by Integration of Bioinformatics Analysis and Molecular Docking Simulation

Background

Ameloblastoma (AM) is a benign tumor locally originated from odontogenic epithelium that is commonly found in the jaw. This tumor makes aggressive invasions and has a high recurrence rate. This study aimed to investigate the differentially expressed genes (DEGs), biological function alterations, disease targets, and existing drugs for AM using bioinformatics analysis.

Methods

The data set of AM was retrieved from the GEO database (GSE132474) and identified the DEGs using bioinformatics analysis. The biological alteration analysis was applied to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Protein-protein interaction (PPI) network analysis and hub gene identification were screened through NetworkAnalyst. The transcription factor-protein network was constructed via OmicsNet. We also identified candidate compounds from L1000CDS2 database. The target of AM and candidate compounds were verified using docking simulation.

Results

Totally, 611 DEGs were identified. The biological function enrichment analysis revealed glycosaminoglycan and GABA (γ-aminobutyric acid) signaling were most significantly up-regulated and down-regulated in AM, respectively. Subsequently, hub genes and transcription factors were screened via the network and showed FOS protein was found in both networks. Furthermore, we evaluated FOS protein to be a therapeutic target in AMs. Candidate compounds were screened and verified using docking simulation. Tanespimycin showed the greatest affinity binding value to bind FOS protein.

Conclusions

This study presented the underlying molecular mechanisms of disease pathogenesis, biological alteration, and important pathways of AMs and provided a candidate compound, Tanespimycin, targeting FOS protein for the treatment of AMs.

HANSynergy: Heterogeneous Graph Attention Network for Drug Synergy Prediction

Drug synergy therapy is a promising strategy for cancer treatment. However, the extensive variety of available drugs and the time-intensive process of determining effective drug combinations through clinical trials pose significant challenges. It requires a reliable method for the rapid and precise selection of drug synergies. In response, various computational strategies have been developed for predicting drug synergies, yet the exploitation of heterogeneous biological network features remains underexplored. In this study, we construct a heterogeneous graph that encompasses diverse biological entities and interactions, utilizing rich data sets from sources, such as DrugCombDB, PubChem, UniProt, and cancer cell line encyclopedia (CCLE). We initialize node feature representations and introduce a novel virtual node to enhance drug representation. Our proposed method, the heterogeneous graph attention network for drug-drug synergy prediction (HANSynergy), has been experimentally validated to demonstrate that the heterogeneous graph attention network can extract key node features, efficiently harness the diversity of information, and further enhance network functionality through the incorporation of a multihead attention mechanism. In the comparative experiment, the highest accuracy (Acc) and area under the curve (AUC) are 0.877 and 0.947, respectively, in DrugCombDB_early data set, demonstrating the superiority of HANSynergy over the competing methods. Moreover, protein-protein interactions are important in understanding the mechanism of action of drugs. The heterogeneous attention mechanism facilitates protein-protein interaction analysis. By analyzing the changes of attention weight before and after heterogeneous network training, we investigated proteins that may be associated with drug combinations. Additionally, case studies align our findings with existing research, underscoring the potential of HANSynergy in drug synergy prediction. This advancement not only contributes to the burgeoning field of drug synergy prediction but also holds the potential to provide valuable insights and uncover new drug synergies for combating cancer.

Prediction of inhibitor development in previously untreated and minimally treated children with severe and moderately severe hemophilia A using a machine-learning network

BACKGROUND

Prediction of inhibitor development in patients with hemophilia A (HA) remains a challenge.

OBJECTIVES

To construct a predictive model for inhibitor development in HA using a network of clinical variables and biomarkers based on the individual similarity network.

METHODS

Previously untreated and minimally treated children with severe/moderately severe HA, participants of the HEMFIL Cohort Study, were followed up until reaching 75 exposure days (EDs) without inhibitor (INH-) or upon inhibitor development (INH+). Clinical data and biological samples were collected before the start of factor (F)VIII replacement (T0). A predictive model (HemfilNET) was built to compare the networks and potential global topological differences between INH- and INH+ at T0, considering the network robustness. For validation, the "leave-one-out" cross-validation technique was employed. Accuracy, precision, recall, and F1-score were used as evaluation metrics for the machine-learning model.

RESULTS

We included 95 children with HA (CHA), of whom 31 (33%) developed inhibitors. The algorithm, featuring 37 variables, identified distinct patterns of networks at T0 for INH+ and INH-. The accuracy of the model was 74.2% for CHA INH+ and 98.4% for INH-. By focusing the analysis on CHA with high-risk F8 mutations for inhibitor development, the accuracy in identifying CHA INH+ increased to 82.1%.

CONCLUSION

Our machine-learning algorithm demonstrated an overall accuracy of 90.5% for predicting inhibitor development in CHA, which further improved when restricting the analysis to CHA with a high-risk F8 genotype. However, our model requires validation in other cohorts. Yet, missing data for some variables hindered more precise predictions.

Extending PROXIMAL to predict degradation pathways of phenolic compounds in the human gut microbiota

Despite significant advances in reconstructing genome-scale metabolic networks, the understanding of cellular metabolism remains incomplete for many organisms. A promising approach for elucidating cellular metabolism is analysing the full scope of enzyme promiscuity, which exploits the capacity of enzymes to bind to non-annotated substrates and generate novel reactions. To guide time-consuming costly experimentation, different computational methods have been proposed for exploring enzyme promiscuity. One relevant algorithm is PROXIMAL, which strongly relies on KEGG to define generic reaction rules and link specific molecular substructures with associated chemical transformations. Here, we present a completely new pipeline, PROXIMAL2, which overcomes the dependency on KEGG data. In addition, PROXIMAL2 introduces two relevant improvements with respect to the former version: i) correct treatment of multi-step reactions and ii) tracking of electric charges in the transformations. We compare PROXIMAL and PROXIMAL2 in recovering annotated products from substrates in KEGG reactions, finding a highly significant improvement in the level of accuracy. We then applied PROXIMAL2 to predict degradation reactions of phenolic compounds in the human gut microbiota. The results were compared to RetroPath RL, a different and relevant enzyme promiscuity method. We found a significant overlap between these two methods but also complementary results, which open new research directions into this relevant question in nutrition.

Relevance of the Adjuvant Effect between Cellular Homeostasis and Resistance to Antibiotics in Gram-Negative Bacteria with Pathogenic Capacity: A Study of Klebsiella pneumoniae

Antibiotic resistance has become a global issue. The most significant risk is the acquisition of these mechanisms by pathogenic bacteria, which can have a severe clinical impact and pose a public health risk. This problem assumes that bacterial fitness is a constant phenomenon and should be approached from an evolutionary perspective to develop the most appropriate and effective strategies to contain the emergence of strains with pathogenic potential. Resistance mechanisms can be understood as adaptive processes to stressful conditions. This review examines the relevance of homeostatic regulatory mechanisms in antimicrobial resistance mechanisms. We focus on the interactions in the cellular physiology of pathogenic bacteria, particularly Gram-negative bacteria, and specifically Klebsiella pneumoniae. From a clinical research perspective, understanding these interactions is crucial for comprehensively understanding the phenomenon of resistance and developing more effective drugs and treatments to limit or attenuate bacterial sepsis, since the most conserved adjuvant phenomena in bacterial physiology has turned out to be more optimized and, therefore, more susceptible to alterations due to pharmacological action.

Insulin receptor at the blood-brain barrier: Transport and signaling

The blood-brain barrier (BBB) is a unique system of the brain microvasculature that limits the exchange between the blood and the brain. Brain microvascular endothelial cells form the BBB as part of the neurovascular unit and express insulin receptors. The insulin receptor at the BBB has been studied in two different functional aspects. These functions include (1) the supplying of blood insulin to the brain and (2) the modulation of BBB function via insulin signaling. The first function involves drug delivery to the brain, while the second function is related to the association between central nervous system diseases and type 2 diabetes through insulin resistance. This chapter summarizes recent progress in research on the function of insulin receptors at the BBB.

Personalized metabolic whole-body models for newborns and infants predict growth and biomarkers of inherited metabolic diseases

Comprehensive whole-body models (WBMs) accounting for organ-specific dynamics have been developed to simulate adult metabolism, but such models do not exist for infants. Here, we present a resource of 360 organ-resolved, sex-specific models of newborn and infant metabolism (infant-WBMs) spanning the first 180 days of life. These infant-WBMs were parameterized to represent the distinct metabolic characteristics of newborns and infants, including nutrition, energy requirements, and thermoregulation. We demonstrate that the predicted infant growth was consistent with the recommendation by the World Health Organization. We assessed the infant-WBMs' reliability and capabilities for personalization by simulating 10,000 newborns based on their blood metabolome and birth weight. Furthermore, the infant-WBMs accurately predicted changes in known biomarkers over time and metabolic responses to treatment strategies for inherited metabolic diseases. The infant-WBM resource holds promise for personalized medicine, as the infant-WBMs could be a first step to digital metabolic twins for newborn and infant metabolism.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

The role of elasticity on adhesion and clustering of neurons on soft surfaces

The question of whether material stiffness enhances cell adhesion and clustering is still open to debate. Results from the literature are seemingly contradictory, with some reports illustrating that adhesion increases with surface stiffness and others suggesting that the performance of a system of cells is curbed by high values of elasticity. To address the role of elasticity as a regulator in neuronal cell adhesion and clustering, we investigated the topological characteristics of networks of neurons on polydimethylsiloxane (PDMS) surfaces - with values of elasticity (E) varying in the 0.55-2.65 MPa range. Results illustrate that, as elasticity increases, the number of neurons adhering on the surface decreases. Notably, the small-world coefficient - a topological measure of networks - also decreases. Numerical simulations and functional multi-calcium imaging experiments further indicated that the activity of neuronal cells on soft surfaces improves for decreasing E. Experimental findings are supported by a mathematical model, that explains adhesion and clustering of cells on soft materials as a function of few parameters - including the Young's modulus and roughness of the material. Overall, results indicate that - in the considered elasticity interval - increasing the compliance of a material improves adhesion, improves clustering, and enhances communication of neurons.

Unlocking Glioblastoma Vulnerabilities with CRISPR-Based Genetic Screening

Glioblastoma (GBM) is the most common malignant brain tumor in adults. Despite advancements in treatment, the prognosis for patients with GBM remains poor due to its aggressive nature and resistance to therapy. CRISPR-based genetic screening has emerged as a powerful tool for identifying genes crucial for tumor progression and treatment resistance, offering promising targets for tumor therapy. In this review, we provide an overview of the recent advancements in CRISPR-based genetic screening approaches and their applications in GBM. We highlight how these approaches have been used to uncover the genetic determinants of GBM progression and responsiveness to various therapies. Furthermore, we discuss the ongoing challenges and future directions of CRISPR-based screening methods in advancing GBM research.

Information Density Enhancement Using Lossy Compression in DNA Data Storage

This study develops two deoxyribonucleic acid (DNA) lossy compression models, Models A and B, to encode grayscale images into DNA sequences, enhance information density, and enable high-fidelity image recovery. These models, distinguished by their handling of pixel domains and interpolation methods, offer a novel approach to data storage for DNA. Model A processes pixels in overlapped domains using linear interpolation (LI), whereas Model B uses non-overlapped domains with nearest-neighbor interpolation (NNI). Through a comparative analysis with Joint Photographic Experts Group (JPEG) compression, the DNA lossy compression models demonstrate competitive advantages in terms of information density and image quality restoration. The application of these models to the Modified National Institute of Standards and Technology (MNIST) dataset reveals their efficiency and the recognizability of decompressed images, which is validated by convolutional neural network (CNN) performance. In particular, Model B2, a version of Model B, emerges as an effective method for balancing high information density (surpassing over 20 times the typical densities of two bits per nucleotide) with reasonably good image quality. These findings highlight the potential of DNA-based data storage systems for high-density and efficient compression, indicating a promising future for biological data storage solutions.

Distinct Network Patterns Emerge from Cartesian and XOR Epistasis Models: A Comparative Network Science Analysis

Background

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, can significantly contribute to the observed phenotypic variance of complex traits. To date, it has been generally assumed that genetic interactions can be detected using a Cartesian, or multiplicative, interaction model commonly utilized in standard regression approaches. However, a recent study investigating epistasis in obesity-related traits in rats and mice has identified potential limitations of the Cartesian model, revealing that it only detects some of the genetic interactions occurring in these systems. By applying an alternative approach, the exclusive-or (XOR) model, the researchers detected a greater number of epistatic interactions and identified more biologically relevant ontological terms associated with the interacting loci. This suggests that the XOR model may provide a more comprehensive understanding of epistasis in these species and phenotypes. To further explore these findings and determine if different interaction models also make up distinct epistatic networks, we leverage network science to provide a more comprehensive view into the genetic interactions underlying BMI in this system.

Results

Our comparative analysis of networks derived from Cartesian and XOR interaction models in rats (Rattus norvegicus) uncovers distinct topological characteristics for each model-derived network. Notably, we discover that networks based on the XOR model exhibit an enhanced sensitivity to epistatic interactions. This sensitivity enables the identification of network communities, revealing novel trait-related biological functions through enrichment analysis. Furthermore, we identify triangle network motifs in the XOR epistatic network, suggestive of higher-order epistasis, based on the topology of lower-order epistasis.

Conclusions

These findings highlight the XOR model's ability to uncover meaningful biological associations as well as higher-order epistasis from lower-order epistatic networks. Additionally, our results demonstrate that network approaches not only enhance epistasis detection capabilities but also provide more nuanced understandings of genetic architectures underlying complex traits. The identification of community structures and motifs within these distinct networks, especially in XOR, points to the potential for network science to aid in the discovery of novel genetic pathways and regulatory networks. Such insights are important for advancing our understanding of phenotype-genotype relationships.

Gut microbial metabolite targets HDAC3-FOXK1-interferon axis in fibroblast-like synoviocytes to ameliorate rheumatoid arthritis

Rheumatoid arthritis (RA) is an autoimmune disease. Early studies hold an opinion that gut microbiota is environmentally acquired and associated with RA susceptibility. However, accumulating evidence demonstrates that genetics also shape the gut microbiota. It is known that some strains of inbred laboratory mice are highly susceptible to collagen-induced arthritis (CIA), while the others are resistant to CIA. Here, we show that transplantation of fecal microbiota of CIA-resistant C57BL/6J mice to CIA-susceptible DBA/1J mice confer CIA resistance in DBA/1J mice. C57BL/6J mice and healthy human individuals have enriched B. fragilis than DBA/1J mice and RA patients. Transplantation of B. fragilis prevents CIA in DBA/1J mice. We identify that B. fragilis mainly produces propionate and C57BL/6J mice and healthy human individuals have higher level of propionate. Fibroblast-like synoviocytes (FLSs) in RA are activated to undergo tumor-like transformation. Propionate disrupts HDAC3-FOXK1 interaction to increase acetylation of FOXK1, resulting in reduced FOXK1 stability, blocked interferon signaling and deactivation of RA-FLSs. We treat CIA mice with propionate and show that propionate attenuates CIA. Moreover, a combination of propionate with anti-TNF etanercept synergistically relieves CIA. These results suggest that B. fragilis or propionate could be an alternative or complementary approach to the current therapies.

Lipidomic Analysis Reveals Systemic Alterations in Servicemen Exposed to Repeated Occupational Low-Level Blast Waves

INTRODUCTION

Occupational exposure to blast is a prevalent risk experienced by military personnel. While low-level exposure may not manifest immediate signs of illness, prolonged and repetitive exposure may result in neurophysiological dysfunction. Such repeated exposure to occupational blasts has been linked to structural and functional modifications in the brain, adversely affecting the performance of servicemen in the field. These neurological changes can give rise to symptoms resembling concussion and contribute to the development of post-traumatic stress disorder.

MATERIALS AND METHODS

To understand long-term effects of blast exposure, the study was conducted to assess memory function, serum circulatory protein and lipid biomarkers, and associated concussive symptomology in servicemen. Concussion-like symptoms were assessed using the Rivermead Post-Concussion Symptoms Questionnaire (RPSQ) along with memory function using PGI memory scale. The serum protein biomarkers were quantified using a sandwich ELISA assay, and the serum lipid profile was measured using liquid chromatography-mass spectrometer.

RESULTS

The findings revealed that repeated low-level blast exposure resulted in impaired memory function, accompanied by elevated levels of serum neurofilament light chain (neuroaxonal injury) and C-reactive protein. Furthermore, alterations in the lipid profile were observed, with an increase in lipid species associated with immune activation. These changes collectively point to systemic inflammation, neuronal injury, and memory dysfunction as pathological characteristics of repeated low-level blast exposure.

CONCLUSION

The results of our preliminary investigation offer valuable insights for further large-scale study and provide a guiding principle that necessitates a suitable mitigation approach to safeguard the health of personnel against blast overpressure.

Integrating gene expression and splicing dynamics across dose-response oxidative modulators

Toxicological risk assessment increasingly utilizes transcriptomics to derive point of departure (POD) and modes of action (MOA) for chemicals. One essential biological process that allows a single gene to generate several different RNA isoforms is called alternative splicing. To comprehensively assess the role of splicing dysregulation in toxicological evaluation and elucidate its potential as a complementary endpoint, we performed RNA-seq on A549 cells treated with five oxidative stress modulators across a wide dose range. Differential gene expression (DGE) showed limited pathway enrichment except at high concentrations. However, alternative splicing analysis revealed variable intron retention events affecting diverse pathways for all chemicals in the absence of significant expression changes. For instance, diazinon elicited negligible gene expression changes but progressive increase in the number of intron retention events, suggesting splicing alterations precede expression responses. Benchmark dose modeling of intron retention data highlighted relevant pathways overlooked by expression analysis. Systematic integration of splicing datasets should be a useful addition to the toxicogenomic toolkit. Combining both modalities paint a more complete picture of transcriptomic dose-responses. Overall, evaluating intron retention dynamics afforded by toxicogenomics may provide biomarkers that can enhance chemical risk assessment and regulatory decision making. This work highlights splicing-aware toxicogenomics as a possible additional tool for examining cellular responses.

NetActivity enhances transcriptional signals by combining gene expression into robust gene set activity scores through interpretable autoencoders

Grouping gene expression into gene set activity scores (GSAS) provides better biological insights than studying individual genes. However, existing gene set projection methods cannot return representative, robust, and interpretable GSAS. We developed NetActivity, a machine learning framework that generates GSAS based on a sparsely-connected autoencoder, where each neuron in the inner layer represents a gene set. We proposed a three-tier training that yielded representative, robust, and interpretable GSAS. NetActivity model was trained with 1518 GO biological processes terms and KEGG pathways and all GTEx samples. NetActivity generates GSAS robust to the initialization parameters and representative of the original transcriptome, and assigned higher importance to more biologically relevant genes. Moreover, NetActivity returns GSAS with a more consistent definition and higher interpretability than GSVA and hipathia, state-of-the-art gene set projection methods. Finally, NetActivity enables combining bulk RNA-seq and microarray datasets in a meta-analysis of prostate cancer progression, highlighting gene sets related to cell division, key for disease progression. When applied to metastatic prostate cancer, gene sets associated with cancer progression were also altered due to drug resistance, while a classical enrichment analysis identified gene sets irrelevant to the phenotype. NetActivity is publicly available in Bioconductor and GitHub.

A Causal Regulation Modeling Algorithm for Temporal Events with Application to Escherichia coli's Aerobic to Anaerobic Transition

Time-series experiments are crucial for understanding the transient and dynamic nature of biological phenomena. These experiments, leveraging advanced classification and clustering algorithms, allow for a deep dive into the cellular processes. However, while these approaches effectively identify patterns and trends within data, they often need to improve in elucidating the causal mechanisms behind these changes. Building on this foundation, our study introduces a novel algorithm for temporal causal signaling modeling, integrating established knowledge networks with sequential gene expression data to elucidate signal transduction pathways over time. Focusing on Escherichia coli's (E. coli) aerobic to anaerobic transition (AAT), this research marks a significant leap in understanding the organism's metabolic shifts. By applying our algorithm to a comprehensive E. coli regulatory network and a time-series microarray dataset, we constructed the cross-time point core signaling and regulatory processes of E. coli's AAT. Through gene expression analysis, we validated the primary regulatory interactions governing this process. We identified a novel regulatory scheme wherein environmentally responsive genes, soxR and oxyR, activate fur, modulating the nitrogen metabolism regulators fnr and nac. This regulatory cascade controls the stress regulators ompR and lrhA, ultimately affecting the cell motility gene flhD, unveiling a novel regulatory axis that elucidates the complex regulatory dynamics during the AAT process. Our approach, merging empirical data with prior knowledge, represents a significant advance in modeling cellular signaling processes, offering a deeper understanding of microbial physiology and its applications in biotechnology.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte- like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

The biomedical knowledge graph of symptom phenotype in coronary artery plaque: machine learning-based analysis of real-world clinical data

A knowledge graph can effectively showcase the essential characteristics of data and is increasingly emerging as a significant means of integrating information in the field of artificial intelligence. Coronary artery plaque represents a significant etiology of cardiovascular events, posing a diagnostic challenge for clinicians who are confronted with a multitude of nonspecific symptoms. To visualize the hierarchical relationship network graph of the molecular mechanisms underlying plaque properties and symptom phenotypes, patient symptomatology was extracted from electronic health record data from real-world clinical settings. Phenotypic networks were constructed utilizing clinical data and protein‒protein interaction networks. Machine learning techniques, including convolutional neural networks, Dijkstra's algorithm, and gene ontology semantic similarity, were employed to quantify clinical and biological features within the network. The resulting features were then utilized to train a K-nearest neighbor model, yielding 23 symptoms, 41 association rules, and 61 hub genes across the three types of plaques studied, achieving an area under the curve of 92.5%. Weighted correlation network analysis and pathway enrichment were subsequently utilized to identify lipid status-related genes and inflammation-associated pathways that could help explain the differences in plaque properties. To confirm the validity of the network graph model, we conducted coexpression analysis of the hub genes to evaluate their potential diagnostic value. Additionally, we investigated immune cell infiltration, examined the correlations between hub genes and immune cells, and validated the reliability of the identified biological pathways. By integrating clinical data and molecular network information, this biomedical knowledge graph model effectively elucidated the potential molecular mechanisms that collude symptoms, diseases, and molecules.

Therapeutic Effects of Plant Anthocyanin against Alzheimer's Disease and Modulate Gut Health, Short-Chain Fatty Acids

Alzheimer's disease (AD) is the most common form of dementia and neurogenerative disease (NDD), and it is also one of the leading causes of death worldwide. The number of AD patients is over 55 million according to 2020 Alzheimer's Disease International (ADI), and the number is increasing drastically without any effective cure. In this review, we discuss and analyze the potential role of anthocyanins (ACNs) against AD while understanding the molecular mechanisms. ACNs have been reported as having neuroprotective effects by mitigating cognitive impairments, apoptotic markers, neuroinflammation, aberrant amyloidogenesis, and tauopathy. Taken together, ACNs could be an important therapeutic agent for combating or delaying the onset of AD.

A gray box framework that optimizes a white box logical model using a black box optimizer for simulating cellular responses to perturbations

Predicting cellular responses to perturbations requires interpretable insights into molecular regulatory dynamics to perform reliable cell fate control, despite the confounding non-linearity of the underlying interactions. There is a growing interest in developing machine learning-based perturbation response prediction models to handle the non-linearity of perturbation data, but their interpretation in terms of molecular regulatory dynamics remains a challenge. Alternatively, for meaningful biological interpretation, logical network models such as Boolean networks are widely used in systems biology to represent intracellular molecular regulation. However, determining the appropriate regulatory logic of large-scale networks remains an obstacle due to the high-dimensional and discontinuous search space. To tackle these challenges, we present a scalable derivative-free optimizer trained by meta-reinforcement learning for Boolean network models. The logical network model optimized by the trained optimizer successfully predicts anti-cancer drug responses of cancer cell lines, while simultaneously providing insight into their underlying molecular regulatory mechanisms.

HMMF: a hybrid multi-modal fusion framework for predicting drug side effect frequencies

BACKGROUND

The identification of drug side effects plays a critical role in drug repositioning and drug screening. While clinical experiments yield accurate and reliable information about drug-related side effects, they are costly and time-consuming. Computational models have emerged as a promising alternative to predict the frequency of drug-side effects. However, earlier research has primarily centered on extracting and utilizing representations of drugs, like molecular structure or interaction graphs, often neglecting the inherent biomedical semantics of drugs and side effects.

RESULTS

To address the previously mentioned issue, we introduce a hybrid multi-modal fusion framework (HMMF) for predicting drug side effect frequencies. Considering the wealth of biological and chemical semantic information related to drugs and side effects, incorporating multi-modal information offers additional, complementary semantics. HMMF utilizes various encoders to understand molecular structures, biomedical textual representations, and attribute similarities of both drugs and side effects. It then models drug-side effect interactions using both coarse and fine-grained fusion strategies, effectively integrating these multi-modal features.

CONCLUSIONS

HMMF exhibits the ability to successfully detect previously unrecognized potential side effects, demonstrating superior performance over existing state-of-the-art methods across various evaluation metrics, including root mean squared error and area under receiver operating characteristic curve, and shows remarkable performance in cold-start scenarios.

Comprehensive molecular interaction map of TGFβ induced epithelial to mesenchymal transition in breast cancer

Breast cancer is one of the prevailing cancers globally, with a high mortality rate. Metastatic breast cancer (MBC) is an advanced stage of cancer, characterised by a highly nonlinear, heterogeneous process involving numerous singling pathways and regulatory interactions. Epithelial-mesenchymal transition (EMT) emerges as a key mechanism exploited by cancer cells. Transforming Growth Factor-β (TGFβ)-dependent signalling is attributed to promote EMT in advanced stages of breast cancer. A comprehensive regulatory map of TGFβ induced EMT was developed through an extensive literature survey. The network assembled comprises of 312 distinct species (proteins, genes, RNAs, complexes), and 426 reactions (state transitions, nuclear translocations, complex associations, and dissociations). The map was developed by following Systems Biology Graphical Notation (SBGN) using Cell Designer and made publicly available using MINERVA ( http://35.174.227.105:8080/minerva/?id=Metastatic_Breast_Cancer_1 ). While the complete molecular mechanism of MBC is still not known, the map captures the elaborate signalling interplay of TGFβ induced EMT-promoting MBC. Subsequently, the disease map assembled was translated into a Boolean model utilising CaSQ and analysed using Cell Collective. Simulations of these have captured the known experimental outcomes of TGFβ induced EMT in MBC. Hub regulators of the assembled map were identified, and their transcriptome-based analysis confirmed their role in cancer metastasis. Elaborate analysis of this map may help in gaining additional insights into the development and progression of metastatic breast cancer.

Emergence of Orchestrated and Dynamic Metabolism of Saccharomyces cerevisiae

Microbial metabolism is a fundamental cellular process that involves many biochemical events and is distinguished by its emergent properties. While the molecular details of individual reactions have been increasingly elucidated, it is not well understood how these reactions are quantitatively orchestrated to produce collective cellular behaviors. Here we developed a coarse-grained, systems, and dynamic mathematical framework, which integrates metabolic reactions with signal transduction and gene regulation to dissect the emergent metabolic traits of Saccharomyces cerevisiae. Our framework mechanistically captures a set of characteristic cellular behaviors, including the Crabtree effect, diauxic shift, diauxic lag time, and differential growth under nutrient-altered environments. It also allows modular expansion for zooming in on specific pathways for detailed metabolic profiles. This study provides a systems mathematical framework for yeast metabolic behaviors, providing insights into yeast physiology and metabolic engineering.

Nitrogen starvation causes lipid remodeling in Rhodotorula toruloides

BACKGROUND

The oleaginous yeast Rhodotorula toruloides is a promising chassis organism for the biomanufacturing of value-added bioproducts. It can accumulate lipids at a high fraction of biomass. However, metabolic engineering efforts in this organism have progressed at a slower pace than those in more extensively studied yeasts. Few studies have investigated the lipid accumulation phenotype exhibited by R. toruloides under nitrogen limitation conditions. Consequently, there have been only a few studies exploiting the lipid metabolism for higher product titers.

RESULTS

We performed a multi-omic investigation of the lipid accumulation phenotype under nitrogen limitation. Specifically, we performed comparative transcriptomic and lipidomic analysis of the oleaginous yeast under nitrogen-sufficient and nitrogen deficient conditions. Clustering analysis of transcriptomic data was used to identify the growth phase where nitrogen-deficient cultures diverged from the baseline conditions. Independently, lipidomic data was used to identify that lipid fractions shifted from mostly phospholipids to mostly storage lipids under the nitrogen-deficient phenotype. Through an integrative lens of transcriptomic and lipidomic analysis, we discovered that R. toruloides undergoes lipid remodeling during nitrogen limitation, wherein the pool of phospholipids gets remodeled to mostly storage lipids. We identify specific mRNAs and pathways that are strongly correlated with an increase in lipid levels, thus identifying putative targets for engineering greater lipid accumulation in R. toruloides. One surprising pathway identified was related to inositol phosphate metabolism, suggesting further inquiry into its role in lipid accumulation.

CONCLUSIONS

Integrative analysis identified the specific biosynthetic pathways that are differentially regulated during lipid remodeling. This insight into the mechanisms of lipid accumulation can lead to the success of future metabolic engineering strategies for overproduction of oleochemicals.

Enhancement of vitamin B6 production driven by omics analysis combined with fermentation optimization

BACKGROUND

Microbial engineering aims to enhance the ability of bacteria to produce valuable products, including vitamin B6 for various applications. Numerous microorganisms naturally produce vitamin B6, yet the metabolic pathways involved are rigorously controlled. This regulation by the accumulation of vitamin B6 poses a challenge in constructing an efficient cell factory.

RESULTS

In this study, we conducted transcriptome and metabolome analyses to investigate the effects of the accumulation of pyridoxine, which is the major commercial form of vitamin B6, on cellular processes in Escherichia coli. Our omics analysis revealed associations between pyridoxine and amino acids, as well as the tricarboxylic acid (TCA) cycle. Based on these findings, we identified potential targets for fermentation optimization, including succinate, amino acids, and the carbon-to-nitrogen (C/N) ratio. Through targeted modifications, we achieved pyridoxine titers of approximately 514 mg/L in shake flasks and 1.95 g/L in fed-batch fermentation.

CONCLUSION

Our results provide insights into pyridoxine biosynthesis within the cellular metabolic network for the first time. Our comprehensive analysis revealed that the fermentation process resulted in a remarkable final yield of 1.95 g/L pyridoxine, the highest reported yield to date. This work lays a foundation for the green industrial production of vitamin B6 in the future.

SARS-CoV-2 NSP5 antagonizes MHC II expression by subverting histone deacetylase 2

SARS-CoV-2 interferes with antigen presentation by downregulating major histocompatibility complex (MHC) II on antigen-presenting cells, but the mechanism mediating this process is unelucidated. Herein, analysis of protein and gene expression in human antigen-presenting cells reveals that MHC II is downregulated by the SARS-CoV-2 main protease, NSP5. This suppression of MHC II expression occurs via decreased expression of the MHC II regulatory protein CIITA. CIITA downregulation is independent of the proteolytic activity of NSP5, and rather, NSP5 delivers HDAC2 to the transcription factor IRF3 at an IRF-binding site within the CIITA promoter. Here, HDAC2 deacetylates and inactivates the CIITA promoter. This loss of CIITA expression prevents further expression of MHC II, with this suppression alleviated by ectopic expression of CIITA or knockdown of HDAC2. These results identify a mechanism by which SARS-CoV-2 limits MHC II expression, thereby delaying or weakening the subsequent adaptive immune response.

GRouNdGAN: GRN-guided simulation of single-cell RNA-seq data using causal generative adversarial networks

We introduce GRouNdGAN, a gene regulatory network (GRN)-guided reference-based causal implicit generative model for simulating single-cell RNA-seq data, in silico perturbation experiments, and benchmarking GRN inference methods. Through the imposition of a user-defined GRN in its architecture, GRouNdGAN simulates steady-state and transient-state single-cell datasets where genes are causally expressed under the control of their regulating transcription factors (TFs). Training on six experimental reference datasets, we show that our model captures non-linear TF-gene dependencies and preserves gene identities, cell trajectories, pseudo-time ordering, and technical and biological noise, with no user manipulation and only implicit parameterization. GRouNdGAN can synthesize cells under new conditions to perform in silico TF knockout experiments. Benchmarking various GRN inference algorithms reveals that GRouNdGAN effectively bridges the existing gap between simulated and biological data benchmarks of GRN inference algorithms, providing gold standard ground truth GRNs and realistic cells corresponding to the biological system of interest.

Negative feedback systems for modelling NF-κB transcription factor oscillatory activity

Low-dimensional negative feedback systems (NFSs) were developed within a signal flow model to describe the oscillatory activities of NF-κB caused by interactions with its inhibitor IκBα. The NFSs were established as 3rd- and 4th-order linear systems containing unperturbed and perturbed negative feedback (NF) loops with constant or time-varying NF strengths and a feed-forward loop. NF-related analytical solutions to the NFSs representing the time courses of NF-κB and IκBα were determined and their exact mathematical relationship was found. The NFS's parameters were determined to fit the experimental time courses of NF-κB in TNF-α-stimulated embryonic fibroblasts, rela-/- embryonic fibroblasts reconstituted with RelA, C9L cells, GFP-p65 knock-in embryonic fibroblasts and embryogenic fibroblasts lacking Iκβ and IκBε, LPS-stimulated IC-21 macrophages treated or not with DCPA, and anti-IgM-stimulated DT40 B-lymphocytes. The unperturbed and perturbed NFSs describing the above biosystems generated isochronous and non-isochronous solutions, depending on a constant or time-varying NF strength, respectively. The oscillation period of the NF-coupled solutions, the phase difference between them and the time delays in the appearance of cytoplasmic IκBα after stimulation of NF-κB were determined. A significant divergence between the IκBα solutions to the NFSs and the IκBα experimental courses led to a rejection of the NF coupling between NF-κB and IκBα in the above biosystems. It was shown that neither the linearity nor the low dimensionality of the NFSs altered the NF relationship and the divergence between the IκBα solutions to the NFS and IκBα experimental time courses. Although the NF relationship between IκBα and NF-κB was not confirmed in all the experimental data analyzed, delayed negative feedback was found in some cases.

Data Valuation with Gradient Similarity

High-quality data is crucial for accurate machine learning and actionable analytics, however, mislabeled or noisy data is a common problem in many domains. Distinguishing low- from high-quality data can be challenging, often requiring expert knowledge and considerable manual intervention. Data Valuation algorithms are a class of methods that seek to quantify the value of each sample in a dataset based on its contribution or importance to a given predictive task. These data values have shown an impressive ability to identify mislabeled observations, and filtering low-value data can boost machine learning performance. In this work, we present a simple alternative to existing methods, termed Data Valuation with Gradient Similarity (DVGS). This approach can be easily applied to any gradient descent learning algorithm, scales well to large datasets, and performs comparably or better than baseline valuation methods for tasks such as corrupted label discovery and noise quantification. We evaluate the DVGS method on tabular, image and RNA expression datasets to show the effectiveness of the method across domains. Our approach has the ability to rapidly and accurately identify low-quality data, which can reduce the need for expert knowledge and manual intervention in data cleaning tasks.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte-like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

Beyond VEGF: Angiopoietin-Tie Signaling Pathway in Diabetic Retinopathy

Complications from diabetic retinopathy such as diabetic macular edema (DME) and proliferative diabetic retinopathy (PDR) constitute leading causes of preventable vision loss in working-age patients. Since vascular endothelial growth factor (VEGF) plays a major role in the pathogenesis of these complications, VEGF inhibitors have been the cornerstone of their treatment. Anti-VEGF monotherapy is an effective but burdensome treatment for DME. However, due to the intensive and burdensome treatment, most patients in routine clinical practice are undertreated, and therefore, their outcomes are compromised. Even in adequately treated patients, persistent DME is reported anywhere from 30% to 60% depending on the drug used. PDR is currently treated by anti-VEGF, panretinal photocoagulation (PRP) or a combination of both. Similarly, a number of eyes, despite these treatments, continue to progress to tractional retinal detachment and vitreous hemorrhage. Clearly there are other molecular pathways other than VEGF involved in the pathogenesis of DME and PDR. One of these pathways is the angiopoietin-Tie signaling pathway. Angiopoietin 1 (Ang1) plays a major role in maintaining vascular quiescence and stability. It acts as a molecular brake against vascular destabilization and inflammation that is usually promoted by angiopoietin 2 (Ang2). Several pathological conditions including chronic hyperglycemia lead to Ang2 upregulation. Recent regulatory approval of the bi-specific antibody, faricimab, may improve long term outcomes in DME. It targets both the Ang/Tie and VEGF pathways. The YOSEMITE and RHINE were multicenter, double-masked, randomized non-inferiority phase 3 clinical trials that compared faricimab to aflibercept in eyes with center-involved DME. At 12 months of follow-up, faricimab demonstrated non-inferior vision gains, improved anatomic outcomes and a potential for extended dosing when compared to aflibercept. The 2-year results of the YOSEMITE and RHINE trials demonstrated that the anatomic and functional results obtained at the 1 year follow-up were maintained. Short term outcomes of previously treated and treatment-naive eyes with DME that were treated with faricimab during routine clinical practice suggest a beneficial effect of faricimab over other agents. Targeting of Ang2 has been reported by several other means including VE-PTP inhibitors, integrin binding peptide and surrobodies.

An integrative multi-omic analysis defines gut microbiota, mycobiota, and metabolic fingerprints in ulcerative colitis patients

Background

Ulcerative colitis (UC) is a multifactorial chronic inflammatory bowel disease (IBD) that affects the large intestine with superficial mucosal inflammation. A dysbiotic gut microbial profile has been associated with UC. Our study aimed to characterize the UC gut bacterial, fungal, and metabolic fingerprints by omic approaches.

Methods

The 16S rRNA- and ITS2-based metataxonomics and gas chromatography-mass spectrometry/solid phase microextraction (GC-MS/SPME) metabolomic analysis were performed on stool samples of 53 UC patients and 37 healthy subjects (CTRL). Univariate and multivariate approaches were applied to separated and integrated omic data, to define microbiota, mycobiota, and metabolic signatures in UC. The interaction between gut bacteria and fungi was investigated by network analysis.

Results

In the UC cohort, we reported the increase of Streptococcus, Bifidobacterium, Enterobacteriaceae, TM7-3, Granulicatella, Peptostreptococcus, Lactobacillus, Veillonella, Enterococcus, Peptoniphilus, Gemellaceae, and phenylethyl alcohol; and we also reported the decrease of Akkermansia; Ruminococcaceae; Ruminococcus; Gemmiger; Methanobrevibacter; Oscillospira; Coprococus; Christensenellaceae; Clavispora; Vishniacozyma; Quambalaria; hexadecane; cyclopentadecane; 5-hepten-2-ol, 6 methyl; 3-carene; caryophyllene; p-Cresol; 2-butenal; indole, 3-methyl-; 6-methyl-3,5-heptadiene-2-one; 5-octadecene; and 5-hepten-2-one, 6 methyl. The integration of the multi-omic data confirmed the presence of a distinctive bacterial, fungal, and metabolic fingerprint in UC gut microbiota. Moreover, the network analysis highlighted bacterial and fungal synergistic and/or divergent interkingdom interactions.

Conclusion

In this study, we identified intestinal bacterial, fungal, and metabolic UC-associated biomarkers. Furthermore, evidence on the relationships between bacterial and fungal ecosystems provides a comprehensive perspective on intestinal dysbiosis and ecological interactions between microorganisms in the framework of UC.

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

Reverse engineering morphogenesis through Bayesian optimization of physics-based models

Morphogenetic programs coordinate cell signaling and mechanical interactions to shape organs. In systems and synthetic biology, a key challenge is determining optimal cellular interactions for predicting organ shape, size, and function. Physics-based models defining the subcellular force distribution facilitate this, but it is challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the experimentally observed organ shapes. This integrative framework employs Gaussian Process Regression, a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that maintain the final organ shape. We calibrated and tested the method on Drosophila wing imaginal discs to study mechanisms that regulate epithelial processes ranging from development to cancer. The parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with imaging data of wing discs perturbed with collagenase. The computational pipeline identifies distinct parameter sets mimicking wild-type shapes. It enables a global sensitivity analysis to support the regulation of actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with experimental imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This workflow is extensible toward reverse-engineering morphogenesis across organ systems and for real-time control of complex multicellular systems.

Chronic Stress-Induced Neuroinflammation: Relevance of Rodent Models to Human Disease

The brain is the central organ of adaptation to stress because it perceives and determines threats that induce behavioral, physiological, and molecular responses. In humans, chronic stress manifests as an enduring consistent feeling of pressure and being overwhelmed for an extended duration. This can result in a persistent proinflammatory response in the peripheral and central nervous system (CNS), resulting in cellular, physiological, and behavioral effects. Compounding stressors may increase the risk of chronic-stress-induced inflammation, which can yield serious health consequences, including mental health disorders. This review summarizes the current knowledge surrounding the neuroinflammatory response in rodent models of chronic stress-a relationship that is continually being defined. Many studies investigating the effects of chronic stress on neuroinflammation in rodent models have identified significant changes in inflammatory modulators, including nuclear factor-κB (NF-κB) and toll-like receptors (TLRs), and cytokines, including tumor necrosis factor-alpha (TNF-α), interleukin (IL)-1β, and IL-6. This suggests that these are key inflammatory factors in the chronic stress response, which may contribute to the establishment of anxiety and depression-like symptoms. The behavioral and neurological effects of modulating inflammatory factors through gene knockdown (KD) and knockout (KO), and conventional and alternative medicine approaches, are discussed.

Risk of Dementia and Alzheimer's Disease Associated With Antidiabetics: A Bayesian Network Meta-Analysis

INTRODUCTION

Dementia risk is substantially elevated in patients with diabetes. However, evidence on dementia risk associated with various antidiabetic regimens is still limited. This study aims to comprehensively investigate the risk of dementia and Alzheimer's disease (AD) associated with various antidiabetic classes.

METHODS

Cochrane Central Register of Controlled Trials, Embase, MEDLINE (PubMed), and Scopus were searched from inception to March 2024 (PROSPERO CRD 42022365927). Observational studies investigating dementia and AD incidences after antidiabetic initiation were identified. Bayesian network meta-analysis was performed to determine dementia and AD risks associated with antidiabetics. Preferred Reporting Items for Systematic Reviews-Network Meta-Analyses (PRISMA-NMA) guidelines were followed. Statistical analysis was performed and updated in November 2023 and March 2024, respectively.

RESULTS

A total of 1,565,245 patients from 16 studies were included. Dementia and AD risks were significantly lower with metformin and sodium glucose co-transporter-2 inhibitors (SGLT2i). Metformin displayed the lowest risk of dementia across diverse antidiabetics, whereas α-glucosidase inhibitors demonstrated the highest risk. SGLT2i exhibited the lowest dementia risk across second-line antidiabetics. Dementia risk was significantly higher with dipeptidyl peptidase-4 inhibitor (DPP4i), metformin, sulfonylureas, and thiazolidinediones (TZD) compared to SGLT2i in the elderly (≥75 years). Dementia risk associated with metformin was substantially lower, regardless of diabetic complication status or baseline A1C.

DISCUSSION

Metformin and SGLT2i demonstrated lower dementia risk than other antidiabetic classes. Patient-specific factors may affect this relationship and cautious interpretation is warranted as metformin is typically initiated at an earlier stage with fewer complications. Hence, further large-scaled clinical trials are required.

De Novo Multi-Omics Pathway Analysis Designed for Prior Data Independent Inference of Cell Signaling Pathways

New tools for cell signaling pathway inference from multi-omics data that are independent of previous knowledge are needed. Here, we propose a new de novo method, the de novo multi-omics pathway analysis (DMPA), to model and combine omics data into network modules and pathways. DMPA was validated with published omics data and was found accurate in discovering reported molecular associations in transcriptome, interactome, phosphoproteome, methylome, and metabolomics data, and signaling pathways in multi-omics data. DMPA was benchmarked against module discovery and multi-omics integration methods and outperformed previous methods in module and pathway discovery especially when applied to datasets of relatively low sample sizes. Transcription factor, kinase, subcellular location, and function prediction algorithms were devised for transcriptome, phosphoproteome, and interactome modules and pathways, respectively. To apply DMPA in a biologically relevant context, interactome, phosphoproteome, transcriptome, and proteome data were collected from analyses carried out using melanoma cells to address gamma-secretase cleavage-dependent signaling characteristics of the receptor tyrosine kinase TYRO3. The pathways modeled with DMPA reflected the predicted function and its direction in validation experiments.

Biomarkers of Alzheimer's Disease Associated with Programmed Cell Death Reveal Four Repurposed Drugs

Alzheimer's disease (AD) is a neurodegenerative disorder and the most common cause of dementia. Programmed cell death (PCD) is mainly characterized by unique morphological features and energy-dependent biochemical processes. The predominant pathway leading to cell death in AD has not been thoroughly analyzed, although there is evidence of neuron loss in AD and numerous pathways of PCD have been associated with this process. A better understanding of the systems biology underlying the relationship between AD and PCD could lead to the development of new therapeutic approaches. To this end, publicly available transcriptome data were examined using bioinformatic methods such as differential gene expression and weighted gene coexpression network analysis (WGCNA) to find PCD-related AD biomarkers. The diagnostic significance of these biomarkers was evaluated using a logistic regression-based predictive model. Using these biomarkers, a multifactorial regulatory network was developed. Last, a drug repositioning study was conducted to propose new drugs for the treatment of AD targeting PCD. The development of 3PM (predictive, preventive, and personalized) drugs for the treatment of AD would be enabled by additional research on the effects of these drugs on this disease.

Scalable integration of multiomic single-cell data using generative adversarial networks

MOTIVATION

Single-cell profiling has become a common practice to investigate the complexity of tissues, organs, and organisms. Recent technological advances are expanding our capabilities to profile various molecular layers beyond the transcriptome such as, but not limited to, the genome, the epigenome, and the proteome. Depending on the experimental procedure, these data can be obtained from separate assays or the very same cells. Yet, integration of more than two assays is currently not supported by the majority of the computational frameworks avaiable.

RESULTS

We here propose a Multi-Omic data integration framework based on Wasserstein Generative Adversarial Networks suitable for the analysis of paired or unpaired data with a high number of modalities (>2). At the core of our strategy is a single network trained on all modalities together, limiting the computational burden when many molecular layers are evaluated.

AVAILABILITY AND IMPLEMENTATION

Source code of our framework is available at https://github.com/vgiansanti/MOWGAN.

TransGEM: a molecule generation model based on Transformer with gene expression data

MOTIVATION

It is difficult to generate new molecules with desirable bioactivity through ligand-based de novo drug design, and receptor-based de novo drug design is constrained by disease target information availability. The combination of artificial intelligence and phenotype-based de novo drug design can generate new bioactive molecules, independent from disease target information. Gene expression profiles can be used to characterize biological phenotypes. The Transformer model can be utilized to capture the associations between gene expression profiles and molecular structures due to its remarkable ability in processing contextual information.

RESULTS

We propose TransGEM (Transformer-based model from gene expression to molecules), which is a phenotype-based de novo drug design model. A specialized gene expression encoder is used to embed gene expression difference values between diseased cell lines and their corresponding normal tissue cells into TransGEM model. The results demonstrate that the TransGEM model can generate molecules with desirable evaluation metrics and property distributions. Case studies illustrate that TransGEM model can generate structurally novel molecules with good binding affinity to disease target proteins. The majority of genes with high attention scores obtained from TransGEM model are associated with the onset of the disease, indicating the potential of these genes as disease targets. Therefore, this study provides a new paradigm for de novo drug design, and it will promote phenotype-based drug discovery.

AVAILABILITY AND IMPLEMENTATION

The code is available at https://github.com/hzauzqy/TransGEM.

Use of multiple pharmacodynamic measures to deconstruct the Nix-TB regimen in a short-course murine model of tuberculosis

A major challenge for tuberculosis (TB) drug development is to prioritize promising combination regimens from a large and growing number of possibilities. This includes demonstrating individual drug contributions to the activity of higher-order combinations. A BALB/c mouse TB infection model was used to evaluate the contributions of each drug and pairwise combination in the clinically relevant Nix-TB regimen [bedaquiline-pretomanid-linezolid (BPaL)] during the first 3 weeks of treatment at human equivalent doses. The rRNA synthesis (RS) ratio, an exploratory pharmacodynamic (PD) marker of ongoing Mycobacterium tuberculosis rRNA synthesis, together with solid culture CFU counts and liquid culture time to positivity (TTP) were used as PD markers of treatment response in lung tissue; and their time-course profiles were mathematically modeled using rate equations with pharmacologically interpretable parameters. Antimicrobial interactions were quantified using Bliss independence and Isserlis formulas. Subadditive (or antagonistic) and additive effects on bacillary load, assessed by CFU and TTP, were found for bedaquiline-pretomanid and linezolid-containing pairs, respectively. In contrast, subadditive and additive effects on rRNA synthesis were found for pretomanid-linezolid and bedaquiline-containing pairs, respectively. Additionally, accurate predictions of the response to BPaL for all three PD markers were made using only the single-drug and pairwise effects together with an assumption of negligible three-way drug interactions. The results represent an experimental and PD modeling approach aimed at reducing combinatorial complexity and improving the cost-effectiveness of in vivo systems for preclinical TB regimen development.

Topological benchmarking of algorithms to infer gene regulatory networks from single-cell RNA-seq data

MOTIVATION

In recent years, many algorithms for inferring gene regulatory networks from single-cell transcriptomic data have been published. Several studies have evaluated their accuracy in estimating the presence of an interaction between pairs of genes. However, these benchmarking analyses do not quantify the algorithms' ability to capture structural properties of networks, which are fundamental, e.g., for studying the robustness of a gene network to external perturbations. Here, we devise a three-step benchmarking pipeline called STREAMLINE that quantifies the ability of algorithms to capture topological properties of networks and identify hubs.

RESULTS

To this aim, we use data simulated from different types of networks as well as experimental data from three different organisms. We apply our benchmarking pipeline to four inference algorithms and provide guidance on which algorithm should be used depending on the global network property of interest.

AVAILABILITY AND IMPLEMENTATION

STREAMLINE is available at https://github.com/ScialdoneLab/STREAMLINE. The data generated in this study are available at https://doi.org/10.5281/zenodo.10710444.

BN-BacArena: Bayesian network extension of BacArena for the dynamic simulation of microbial communities

MOTIVATION

Simulating gut microbial dynamics is extremely challenging. Several computational tools, notably the widely used BacArena, enable modeling of dynamic changes in the microbial environment. These methods, however, do not comprehensively account for microbe-microbe stimulant or inhibitory effects or for nutrient-microbe inhibitory effects, typically observed in different compounds present in the daily diet.

RESULTS

Here, we present BN-BacArena, an extension of BacArena consisting on the incorporation within the native computational framework of a Bayesian network model that accounts for microbe-microbe and nutrient-microbe interactions. Using in vitro experiments, 16S rRNA gene sequencing data and nutritional composition of 55 foods, the output Bayesian network showed 23 significant nutrient-bacteria interactions, suggesting the importance of compounds such as polyols, ascorbic acid, polyphenols and other phytochemicals, and 40 bacteria-bacteria significant relationships. With test data, BN-BacArena demonstrates a statistically significant improvement over BacArena to predict the time-dependent relative abundance of bacterial species involved in the gut microbiota upon different nutritional interventions. As a result, BN-BacArena opens new avenues for the dynamic modeling and simulation of the human gut microbiota metabolism.

AVAILABILITY AND IMPLEMENTATION

MATLAB and R code are available in https://github.com/PlanesLab/BN-BacArena.

SLIDE: Significant Latent Factor Interaction Discovery and Exploration across biological domains

Modern multiomic technologies can generate deep multiscale profiles. However, differences in data modalities, multicollinearity of the data, and large numbers of irrelevant features make analyses and integration of high-dimensional omic datasets challenging. Here we present Significant Latent Factor Interaction Discovery and Exploration (SLIDE), a first-in-class interpretable machine learning technique for identifying significant interacting latent factors underlying outcomes of interest from high-dimensional omic datasets. SLIDE makes no assumptions regarding data-generating mechanisms, comes with theoretical guarantees regarding identifiability of the latent factors/corresponding inference, and has rigorous false discovery rate control. Using SLIDE on single-cell and spatial omic datasets, we uncovered significant interacting latent factors underlying a range of molecular, cellular and organismal phenotypes. SLIDE outperforms/performs at least as well as a wide range of state-of-the-art approaches, including other latent factor approaches. More importantly, it provides biological inference beyond prediction that other methods do not afford. Thus, SLIDE is a versatile engine for biological discovery from modern multiomic datasets.

Realizing the promise of Project Optimus: Challenges and emerging opportunities for dose optimization in oncology drug development

Project Optimus is a US Food and Drug Administration Oncology Center of Excellence initiative aimed at reforming the dose selection and optimization paradigm in oncology drug development. This project seeks to bring together pharmaceutical companies, international regulatory agencies, academic institutions, patient advocates, and other stakeholders. Although there is much promise in this initiative, there are several challenges that need to be addressed, including multidimensionality of the dose optimization problem in oncology, the heterogeneity of cancer and patients, importance of evaluating long-term tolerability beyond dose-limiting toxicities, and the lack of reliable biomarkers for long-term efficacy. Through the lens of Totality of Evidence and with the mindset of model-informed drug development, we offer insights into dose optimization by building a quantitative knowledge base integrating diverse sources of data and leveraging quantitative modeling tools to build evidence for drug dosage considering exposure, disease biology, efficacy, toxicity, and patient factors. We believe that rational dose optimization can be achieved in oncology drug development, improving patient outcomes by maximizing therapeutic benefit while minimizing toxicity.

Pathomics Signature for Prognosis and Chemotherapy Benefits in Stage III Colon Cancer

Importance

The current TNM staging system may not provide adequate information for prognostic purposes and to assess the potential benefits of chemotherapy for patients with stage III colon cancer.

Objective

To develop and validate a pathomics signature to estimate prognosis and benefit from chemotherapy using hematoxylin-eosin (H-E)-stained slides.

Design, Setting, and Participants

This retrospective prognostic study used data from consecutive patients with histologically confirmed stage III colon cancer at 2 medical centers between January 2012 and December 2015. A total of 114 pathomics features were extracted from digital H-E-stained images from Nanfang Hospital of Southern Medical University, Guangzhou, China, and a pathomics signature was constructed using a least absolute shrinkage and selection operator Cox regression model in the training cohort. The associations of the pathomics signature with disease-free survival (DFS) and overall survival (OS) were evaluated. Patients at the Sixth Affiliated Hospital, Sun Yat-sen University, Guangzhou, China, formed the validation cohort. Data analysis was conducted from September 2022 to March 2023.

Main Outcomes and Measures

The prognostic accuracy of the pathomics signature as well as its association with chemotherapy response were evaluated.

Results

This study included 785 patients (mean [SD] age, 62.7 [11.1] years; 437 [55.7%] male). A pathomics signature was constructed based on 4 features. Multivariable analysis revealed that the pathomics signature was an independent factor associated with DFS (hazard ratio [HR], 2.46 [95% CI, 2.89-4.13]; P < .001) and OS (HR, 2.78 [95% CI, 2.34-3.31]; P < .001) in the training cohort. Incorporating the pathomics signature into pathomics nomograms resulted in better performance for the estimation of prognosis than the traditional model in a concordance index comparison in the training cohort (DFS: HR, 0.88 [95% CI, 0.86-0.89] vs HR, 0.73 [95% CI, 0.71-0.75]; P < .001; OS: HR, 0.85 [95% CI, 0.84-0.86] vs HR, 0.74 [95% CI, 0.72-0.76]; P < .001) and validation cohort (DFS: HR, 0.83 [95% CI, 0.82-0.85] vs HR, 0.70 [95% CI, 0.67-0.72]; P < .001; OS: HR, 0.80 [95% CI, 0.78-0.82] vs HR, 0.69 [0.67-0.72]; P < .001). Further analysis revealed that patients with a low pathomics signature were more likely to benefit from chemotherapy (eg, combined cohort: DFS: HR, 0.44 [95% CI, 0.28-0.69]; P = .001; OS: HR, 0.43 [95% CI, 0.29-0.64]; P < .001).

Conclusions and Relevance

These findings suggest that a pathomics signature could help identify patients most likely to benefit from chemotherapy in stage III colon cancer.

Organoid-guided precision hepatology for metabolic liver disease

Metabolic dysfunction-associated steatotic liver disease affects millions of people worldwide. Progress towards a definitive cure has been incremental and treatment is currently limited to lifestyle modification. Hepatocyte-specific lipid accumulation is the main trigger of lipotoxic events, driving inflammation and fibrosis. The underlying pathology is extraordinarily heterogenous, and the manifestations of steatohepatitis are markedly influenced by metabolic communications across non-hepatic organs. Synthetic human tissue models have emerged as powerful platforms to better capture the mechanistic diversity in disease progression, while preserving person-specific genetic traits. In this review, we will outline current research efforts focused on integrating multiple synthetic tissue models of key metabolic organs, with an emphasis on organoid-based systems. By combining functional genomics and population-scale en masse profiling methodologies, human tissues derived from patients can provide insights into personalised genetic, transcriptional, biochemical, and metabolic states. These collective efforts will advance our understanding of steatohepatitis and guide the development of rational solutions for mechanism-directed diagnostic and therapeutic investigation.

Elucidating uptake and metabolic fate of dipeptides in CHO cell cultures using 13C labeling experiments and kinetic modeling

The rapidly growing market of biologics including monoclonal antibodies has stimulated the need to improve biomanufacturing processes including mammalian host systems such as Chinese Hamster Ovary (CHO) cells. Cell culture media formulations continue to be enhanced to enable intensified cell culture processes and optimize cell culture performance. Amino acids, major components of cell culture media, are consumed in large amounts by CHO cells. Due to their low solubility and poor stability, certain amino acids including tyrosine, leucine, and phenylalanine can pose major challenges leading to suboptimal bioprocess performance. Dipeptides have the potential to replace amino acids in culture media. However, very little is known about the cleavage, uptake, and utilization kinetics of dipeptides in CHO cell cultures. In this study, replacing amino acids, including leucine and tyrosine by their respective dipeptides including but not limited to Ala-Leu and Gly-Tyr, supported similar cell growth, antibody production, and lactate profiles. Using 13C labeling techniques and spent media studies, dipeptides were shown to undergo both intracellular and extracellular cleavage in cultures. Extracellular cleavage increased with the culture duration, indicating cleavage by host cell proteins that are likely secreted and accumulate in cell culture over time. A kinetic model was built and for the first time, integrated with 13C labeling experiments to estimate dipeptide utilization rates, in CHO cell cultures. Dipeptides with alanine at the N-terminus had a higher utilization rate than dipeptides with alanine at the C-terminus and dipeptides with glycine instead of alanine at N-terminus. Simultaneous supplementation of more than one dipeptide in culture led to reduction in individual dipeptide utilization rates indicating that dipeptides compete for the same cleavage enzymes, transporters, or both. Dipeptide utilization rates in culture and cleavage rates in cell-free experiments appeared to follow Michaelis-Menten kinetics, reaching a maximum at higher dipeptide concentrations. Dipeptide utilization behavior was found to be similar in cell-free and cell culture environments, paving the way for future testing approaches for dipeptides in cell-free environments prior to use in large-scale bioreactors. Thus, this study provides a deeper understanding of the fate of dipeptides in CHO cell cultures through an integration of cell culture, 13C labeling, and kinetic modeling approaches providing insights in how to best use dipeptides in media formulations for robust and optimal mammalian cell culture performance.