Optimization of metabolomic data processing using NOREVA

An integrated in silico approach for the identification of novel potential drug target and chimeric vaccine against Neisseria meningitides strain 331401 serogroup X by subtractive genomics and reverse vaccinology

Neisseria meningitidis, commonly known as the meningococcus, leads to substantial illness and death among children and young adults globally, revealing as either epidemic or sporadic meningitis and/or septicemia. In this study, we have designed a novel peptide-based chimeric vaccine candidate against the N. meningitidis strain 331,401 serogroup X. Through rigorous analysis of subtractive genomics, two essential cytoplasmic proteins, namely UPI000012E8E0(UDP-3-O-acyl-GlcNAc deacetylase) and UPI0000ECF4A9(UDP-N-acetylglucosamine acyltransferase) emerged as potential drug targets. Additionally, using reverse vaccinology, the outer membrane protein UPI0001F4D537 (Membrane fusion protein MtrC) identified by subcellular localization and recognized for its known indispensable role in bacterial survival was identified as a novel chimeric vaccine target. Following a careful comparison of MHC-I, MHC-II, T-cell, and B-cell epitopes, three epitopes derived from UPI0001F4D537 were linked with three types of linkers-GGGS, EAAAK, and the essential PADRE-for vaccine construction. This resulted in eight distinct vaccine models (V1-V8). Among them V1 model was selected as the final vaccine construct. It exhibits exceptional immunogenicity, safety, and enhanced antigenicity, with 97.7 % of its residues in the Ramachandran plot's most favored region. Subsequently, the vaccine structure was docked with the TLR4/MD2 complex and six different HLA allele receptors using the HADDOCK server. The docking resulted in the lowest HADDOCK score of 39.3 ± 9.0 for TLR/MD2. Immune stimulation showed a strong immune response, including antibodies creation and the activation of B-cells, T Cytotoxic cells, T Helper cells, Natural Killer cells, and interleukins. Furthermore, the vaccine construct was successfully expressed in the Escherichia coli system by reverse transcription, optimization, and ligation in the pET-28a (+) vector for the expression study. The current study proposes V1 construct has the potential to elicit both cellular and humoral responses, crucial for the developing an epitope-based vaccine against N. meningitidis strain 331,401 serogroup X.

Legacy and alternative per- and polyfluoroalkyl substances (PFAS) alter the lipid profile of HepaRG cells

Per- and polyfluoroalkyl substances (PFAS) are synthetic chemicals used in various industrial and consumer products. They have gained attention due to their ubiquitous occurrence in the environment and potential for adverse effects on human health, often linked to immune suppression, hepatotoxicity, and altered cholesterol metabolism. This study aimed to explore the impact of ten individual PFAS, 3H-perfluoro-3-[(3-methoxypropoxy) propanoic acid] (PMPP/Adona), ammonium perfluoro-(2-methyl-3-oxahexanoate) (HFPO-DA/GenX), perfluorobutanoic acid (PFBA), perfluorobutanesulfonic acid (PFBS), perfluorodecanoic acid (PFDA), perfluorohexanoic acid (PFHxA), perfluorohexanesulfonate (PFHxS), perfluorononanoic acid (PFNA), perfluorooctanoic acid (PFOA), and perfluorooctanesulfonic acid (PFOS) on the lipid metabolism in human hepatocyte-like cells (HepaRG). These cells were exposed to different concentrations of PFAS ranging from 10µM to 5000µM. Lipids were extracted and analyzed using liquid chromatography coupled with mass spectrometry (LC- MS-QTOF). PFOS at 10µM and PFOA at 25µM increased the levels of ceramide (Cer), diacylglycerol (DAG), N-acylethanolamine (NAE), phosphatidylcholine (PC), and triacylglycerol (TAG) lipids, while PMPP/Adona, HFPO-DA/GenX, PFBA, PFBS, PFHxA, and PFHxS decreased the levels of these lipids. Furthermore, PFOA and PFOS markedly reduced the levels of palmitic acid (FA 16.0). The present study shows distinct concentration-dependent effects of PFAS on various lipid species, shedding light on the implications of PFAS for essential cellular functions. Our study revealed that the investigated legacy PFAS (PFOS, PFOA, PFBA, PFDA, PFHxA, PFHxS, and PFNA) and alternative PFAS (PMPP/Adona, HFPO-DA/GenX and PFBS) can potentially disrupt lipid homeostasis and metabolism in hepatic cells. This research offers a comprehensive insight into the impacts of legacy and alternative PFAS on lipid composition in HepaRG cells.

Precision healthcare: A deep dive into machine learning algorithms and feature selection strategies for accurate heart disease prediction

This paper presents a comprehensive exploration of machine learning algorithms (MLAs) and feature selection techniques for accurate heart disease prediction (HDP) in modern healthcare. By focusing on diverse datasets encompassing various challenges, the research sheds light on optimal strategies for early detection. MLAs such as Decision Trees (DT), Random Forests (RF), Support Vector Machines (SVM), Gaussian Naive Bayes (NB), and others were studied, with precision and recall metrics emphasized for robust predictions. Our study addresses challenges in real-world data through data cleaning and one-hot encoding, enhancing the integrity of our predictive models. Feature extraction techniques-Recursive Feature Extraction (RFE), Principal Component Analysis (PCA), and univariate feature selection-play a crucial role in identifying relevant features and reducing data dimensionality. Our findings showcase the impact of these techniques on improving prediction accuracy. Optimized models for each dataset have been achieved through grid search hyperparameter tuning, with configurations meticulously outlined. Notably, a remarkable 99.12 % accuracy was achieved on the first Kaggle dataset, showcasing the potential for accurate HDP. Model robustness across diverse datasets was highlighted, with caution against overfitting. The study emphasizes the need for validation of unseen data and encourages ongoing research for generalizability. Serving as a practical guide, this research aids researchers and practitioners in HDP model development, influencing clinical decisions and healthcare resource allocation. By providing insights into effective algorithms and techniques, the paper contributes to reducing heart disease-related morbidity and mortality, supporting the healthcare community's ongoing efforts.

Discovery of the potential biomarkers for early diagnosis of endometrial cancer via integrating metabolomics and transcriptomics

Endometrial cancer (EC) is one of the most common malignant tumors in women, and the increasing incidence and mortality pose a serious threat to the public health. Early diagnosis of EC could prolong the survival period and optimize the survivorship, greatly alleviating patients' suffering and social medical pressure. In this study, we collected urine and serum samples from the recruited patients, analyzed the samples using LC-MS approach, and identified the differential metabolites through metabolomic analysis. Then, the differentially expressed genes were identified through the systematic transcriptomic analysis of EC-related dataset from Gene Expression Omnibus (GEO), followed by network profiling of metabolic-reaction-enzyme-gene. In this experiment, a total of 83 differential metabolites and 19 hub genes were discovered, of which 10 different metabolites and 3 hub genes were further evaluated as more potential biomarkers based on network analysis. According to the KEGG enrichment analysis, the potential biomarkers and gene-encoded proteins were found to be involved in the arginine and proline metabolism, histidine metabolism, and pyrimidine metabolism, which was of significance for the early diagnosis of EC. In particular, the combination of metabolites (histamine, 1-methylhistamine, and methylimidazole acetaldehyde) as well as the combination of RRM2, TYMS and TK1 exerted more accurate discrimination abilities between EC and healthy groups, providing more criteria for the early diagnosis of EC.

MOINER: A Novel Multiomics Early Integration Framework for Biomedical Classification and Biomarker Discovery

In the context of precision medicine, multiomics data integration provides a comprehensive understanding of underlying biological processes and is critical for disease diagnosis and biomarker discovery. One commonly used integration method is early integration through concatenation of multiple dimensionally reduced omics matrices due to its simplicity and ease of implementation. However, this approach is seriously limited by information loss and lack of latent feature interaction. Herein, a novel multiomics early integration framework (MOINER) based on information enhancement and image representation learning is thus presented to address the challenges. MOINER employs the self-attention mechanism to capture the intrinsic correlations of omics-features, which make it significantly outperform the existing state-of-the-art methods for multiomics data integration. Moreover, visualizing the attention embedding and identifying potential biomarkers offer interpretable insights into the prediction results. All source codes and model for MOINER are freely available https://github.com/idrblab/MOINER.

Demonstrating the reliability of in vivo metabolomics based chemical grouping: towards best practice

While grouping/read-across is widely used to fill data gaps, chemical registration dossiers are often rejected due to weak category justifications based on structural similarity only. Metabolomics provides a route to robust chemical categories via evidence of shared molecular effects across source and target substances. To gain international acceptance, this approach must demonstrate high reliability, and best-practice guidance is required. The MetAbolomics ring Trial for CHemical groupING (MATCHING), comprising six industrial, government and academic ring-trial partners, evaluated inter-laboratory reproducibility and worked towards best-practice. An independent team selected eight substances (WY-14643, 4-chloro-3-nitroaniline, 17α-methyl-testosterone, trenbolone, aniline, dichlorprop-p, 2-chloroaniline, fenofibrate); ring-trial partners were blinded to their identities and modes-of-action. Plasma samples were derived from 28-day rat tests (two doses per substance), aliquoted, and distributed to partners. Each partner applied their preferred liquid chromatography-mass spectrometry (LC-MS) metabolomics workflows to acquire, process, quality assess, statistically analyze and report their grouping results to the European Chemicals Agency, to ensure the blinding conditions of the ring trial. Five of six partners, whose metabolomics datasets passed quality control, correctly identified the grouping of eight test substances into three categories, for both male and female rats. Strikingly, this was achieved even though a range of metabolomics approaches were used. Through assessing intrastudy quality-control samples, the sixth partner observed high technical variation and was unable to group the substances. By comparing workflows, we conclude that some heterogeneity in metabolomics methods is not detrimental to consistent grouping, and that assessing data quality prior to grouping is essential. We recommend development of international guidance for quality-control acceptance criteria. This study demonstrates the reliability of metabolomics for chemical grouping and works towards best-practice.

FERREG: ferroptosis-based regulation of disease occurrence, progression and therapeutic response

Ferroptosis is a non-apoptotic, iron-dependent regulatory form of cell death characterized by the accumulation of intracellular reactive oxygen species. In recent years, a large and growing body of literature has investigated ferroptosis. Since ferroptosis is associated with various physiological activities and regulated by a variety of cellular metabolism and mitochondrial activity, ferroptosis has been closely related to the occurrence and development of many diseases, including cancer, aging, neurodegenerative diseases, ischemia-reperfusion injury and other pathological cell death. The regulation of ferroptosis mainly focuses on three pathways: system Xc-/GPX4 axis, lipid peroxidation and iron metabolism. The genes involved in these processes were divided into driver, suppressor and marker. Importantly, small molecules or drugs that mediate the expression of these genes are often good treatments in the clinic. Herein, a newly developed database, named 'FERREG', is documented to (i) providing the data of ferroptosis-related regulation of diseases occurrence, progression and drug response; (ii) explicitly describing the molecular mechanisms underlying each regulation; and (iii) fully referencing the collected data by cross-linking them to available databases. Collectively, FERREG contains 51 targets, 718 regulators, 445 ferroptosis-related drugs and 158 ferroptosis-related disease responses. FERREG can be accessed at https://idrblab.org/ferreg/.

Strategy for Identifying a Robust Metabolomic Signature Reveals the Altered Lipid Metabolism in Pituitary Adenoma

Despite the well-established connection between systematic metabolic abnormalities and the pathophysiology of pituitary adenoma (PA), current metabolomic studies have reported an extremely limited number of metabolites associated with PA. Moreover, there was very little consistency in the identified metabolite signatures, resulting in a lack of robust metabolic biomarkers for the diagnosis and treatment of PA. Herein, we performed a global untargeted plasma metabolomic profiling on PA and identified a highly robust metabolomic signature based on a strategy. Specifically, this strategy is unique in (1) integrating repeated random sampling and a consensus evaluation-based feature selection algorithm and (2) evaluating the consistency of metabolomic signatures among different sample groups. This strategy demonstrated superior robustness and stronger discriminative ability compared with that of other feature selection methods including Student's t-test, partial least-squares-discriminant analysis, support vector machine recursive feature elimination, and random forest recursive feature elimination. More importantly, a highly robust metabolomic signature comprising 45 PA-specific differential metabolites was identified. Moreover, metabolite set enrichment analysis of these potential metabolic biomarkers revealed altered lipid metabolism in PA. In conclusion, our findings contribute to a better understanding of the metabolic changes in PA and may have implications for the development of diagnostic and therapeutic approaches targeting lipid metabolism in PA. We believe that the proposed strategy serves as a valuable tool for screening robust, discriminating metabolic features in the field of metabolomics.

MultiClassMetabo: A Superior Classification Model Constructed Using Metabolic Markers in Multiclass Metabolomics

Multiclass metabolomics has become a popular technique for revealing the mechanisms underlying certain physiological processes, different tumor types, or different therapeutic responses. In multiclass metabolomics, it is highly important to uncover the underlying biological information on biosamples by identifying the metabolic markers with the most associations and classifying the different sample classes. The classification problem of multiclass metabolomics is more difficult than that of the binary problem. To date, various methods exist for constructing classification models and identifying metabolic markers consisting of well-established techniques and newly emerging machine learning algorithms. However, how to construct a superior classification model using these methods remains unclear for a given multiclass metabolomic data set. Herein, MultiClassMetabo has been developed for constructing a superior classification model using metabolic markers identified in multiclass metabolomics. MultiClassMetabo can enable online services, including (a) identifying metabolic markers by marker identification methods, (b) constructing classification models by classification methods, and (c) performing a comprehensive assessment from multiple perspectives to construct a superior classification model for multiclass metabolomics. In summary, MultiClassMetabo is distinguished for its capability to construct a superior classification model using the most appropriate method through a comprehensive assessment, which makes it an important complement to other available tools in multiclass metabolomics. MultiClassMetabo can be accessed at http://idrblab.cn/multiclassmetabo/.

VARIDT 3.0: the phenotypic and regulatory variability of drug transporter

The phenotypic and regulatory variability of drug transporter (DT) are vital for the understanding of drug responses, drug-drug interactions, multidrug resistances, and so on. The ADME property of a drug is collectively determined by multiple types of variability, such as: microbiota influence (MBI), transcriptional regulation (TSR), epigenetics regulation (EGR), exogenous modulation (EGM) and post-translational modification (PTM). However, no database has yet been available to comprehensively describe these valuable variabilities of DTs. In this study, a major update of VARIDT was therefore conducted, which gave 2072 MBIs, 10 610 TSRs, 46 748 EGRs, 12 209 EGMs and 10 255 PTMs. These variability data were closely related to the transportation of 585 approved and 301 clinical trial drugs for treating 572 diseases. Moreover, the majority of the DTs in this database were found with multiple variabilities, which allowed a collective consideration in determining the ADME properties of a drug. All in all, VARIDT 3.0 is expected to be a popular data repository that could become an essential complement to existing pharmaceutical databases, and is freely accessible without any login requirement at: https://idrblab.org/varidt/.

TTD: Therapeutic Target Database describing target druggability information

Target discovery is one of the essential steps in modern drug development, and the identification of promising targets is fundamental for developing first-in-class drug. A variety of methods have emerged for target assessment based on druggability analysis, which refers to the likelihood of a target being effectively modulated by drug-like agents. In the therapeutic target database (TTD), nine categories of established druggability characteristics were thus collected for 426 successful, 1014 clinical trial, 212 preclinical/patented, and 1479 literature-reported targets via systematic review. These characteristic categories were classified into three distinct perspectives: molecular interaction/regulation, human system profile and cell-based expression variation. With the rapid progression of technology and concerted effort in drug discovery, TTD and other databases were highly expected to facilitate the explorations of druggability characteristics for the discovery and validation of innovative drug target. TTD is now freely accessible at: https://idrblab.org/ttd/.

TheMarker: a comprehensive database of therapeutic biomarkers

Distinct from the traditional diagnostic/prognostic biomarker (adopted as the indicator of disease state/process), the therapeutic biomarker (ThMAR) has emerged to be very crucial in the clinical development and clinical practice of all therapies. There are five types of ThMAR that have been found to play indispensable roles in various stages of drug discovery, such as: Pharmacodynamic Biomarker essential for guaranteeing the pharmacological effects of a therapy, Safety Biomarker critical for assessing the extent or likelihood of therapy-induced toxicity, Monitoring Biomarker indispensable for guiding clinical management by serially measuring patients' status, Predictive Biomarker crucial for maximizing the clinical outcome of a therapy for specific individuals, and Surrogate Endpoint fundamental for accelerating the approval of a therapy. However, these data of ThMARs has not been comprehensively described by any of the existing databases. Herein, a database, named 'TheMarker', was therefore constructed to (a) systematically offer all five types of ThMAR used at different stages of drug development, (b) comprehensively describe ThMAR information for the largest number of drugs among available databases, (c) extensively cover the widest disease classes by not just focusing on anticancer therapies. These data in TheMarker are expected to have great implication and significant impact on drug discovery and clinical practice, and it is freely accessible without any login requirement at: https://idrblab.org/themarker.

Evaluating Machine Learning Methods of Analyzing Multiclass Metabolomics

Multiclass metabolomic studies have become popular for revealing the differences in multiple stages of complex diseases, various lifestyles, or the effects of specific treatments. In multiclass metabolomics, there are multiple data manipulation steps for analyzing raw data, which consist of data filtering, the imputation of missing values, data normalization, marker identification, sample separation, classification, and so on. In each step, several to dozens of machine learning methods can be chosen for the given data set, with potentially hundreds or thousands of method combinations in the whole data processing chain. Therefore, a clear understanding of these machine learning methods is helpful for selecting an appropriate method combination for obtaining stable and reliable analytical results of specific data. However, there has rarely been an overall introduction or evaluation of these methods based on multiclass metabolomic data. Herein, detailed descriptions of these machine learning methods in multiple data manipulation steps are reviewed. Moreover, an assessment of these methods was performed using a benchmark data set for multiclass metabolomics. First, 12 imputation methods for imputing missing values were evaluated based on the PSS (Procrustes statistical shape analysis) and NRMSE (normalized root-mean-square error) values. Second, 17 normalization methods for processing multiclass metabolomic data were evaluated by applying the PMAD (pooled median absolute deviation) value. Third, different methods of identifying markers of multiclass metabolomics were evaluated based on the CWrel (relative weighted consistency) value. Fourth, nine classification methods for constructing multiclass models were assessed using the AUC (area under the curve) value. Performance evaluations of machine learning methods are highly recommended to select the most appropriate method combination before performing the final analysis of the given data. Overall, detailed descriptions and evaluation of various machine learning methods are expected to improve analyses of multiclass metabolomic data.

Integrative deep learning framework predicts lipidomics-based investigation of preservatives on meat nutritional biomarkers and metabolic pathways

Preservatives are added as antimicrobial agents to extend the shelf life of meat. Adding preservatives to meat products can affect their flavor and nutrition. This review clarifies the effects of preservatives on metabolic pathways and network molecular transformations in meat products based on lipidomics, metabolomics and proteomics analyses. Preservatives change the nutrient content of meat products via altering ionic strength and pH to influence enzyme activity. Ionic strength in salt triggers muscle triglyceride hydrolysis by causing phosphorylation and lipid droplet splitting in adipose tissue hormone-sensitive lipase and triglyceride lipase. DisoLipPred exploiting deep recurrent networks and transfer learning can predict the lipid binding trend of each amino acid in the disordered region of input protein sequences, which could provide omics analyses of biomarkers metabolic pathways in meat products. While conventional meat quality assessment tools are unable to elucidate the intrinsic mechanisms and pathways of variables in the influences of preservatives on the quality of meat products, the promising application of omics techniques in food analysis and discovery through multimodal learning prediction algorithms of neural networks (e.g., deep neural network, convolutional neural network, artificial neural network) will drive the meat industry to develop new strategies for food spoilage prevention and control.

Integrated proteomic and metabolomic analysis of plasma reveals regulatory pathways and key elements in thyroid cancer

Thyroid cancer (TC) is the most common endocrine malignancy with increasing incidence in recent years. Fine-needle aspiration biopsy (FNAB), as a gold standard for the initial evaluation of thyroid nodules, fails to cover all the cytopathologic conditions resulting in overdiagnosis. There is an urgent need for a better classification of thyroid cancer from benign thyroid nodules (BTNs). Here, data independent acquisition (DIA)-based proteomics and untargeted metabolomics in plasma samples of 10 patients with TC and 15 patients with BTNs were performed. Key proteins and metabolites were identified specific to TC, and an independent cohort was used to validate the potential biomarkers using enzyme-linked immunosorbent assay (ELISA). In total, 1429 proteins and 1172 metabolites were identified. Principal component analysis showed a strong overlap at the proteomic level and a significant discrimination at the metabolomic level between the two groups, indicating a more drastic disturbance in the metabolome of thyroid cancer. Integrated analysis of proteomics and metabolomics shows glycerophospholipid metabolism and arachidonic acid metabolism as key regulatory pathways. Furthermore, a multi-omics biomarker panel was developed consisting of LCAT, GPX3 and leukotriene B4. Based on the AUC value for the discovery set, the classification performance was 0.960. The AUC value of the external validation set was 0.930. Altogether, our results will contribute to the clinical application of potential biomarkers in the diagnosis of thyroid cancer.

Opportunities and challenges for sample preparation and enrichment in mass spectrometry for single-cell metabolomics

Single-cell heterogeneity in metabolism, drug resistance and disease type poses the need for analytical techniques for single-cell analysis. As the metabolome provides the closest view of the status quo in the cell, studying the metabolome at single-cell resolution may unravel said heterogeneity. A challenge in single-cell metabolome analysis is that metabolites cannot be amplified, so one needs to deal with picolitre volumes and a wide range of analyte concentrations. Due to high sensitivity and resolution, MS is preferred in single-cell metabolomics. Large numbers of cells need to be analysed for proper statistics; this requires high-throughput analysis, and hence automation of the analytical workflow. Significant advances in (micro)sampling methods, CE and ion mobility spectrometry have been made, some of which have been applied in high-throughput analyses. Microfluidics has enabled an automation of cell picking and metabolite extraction; image recognition has enabled automated cell identification. Many techniques have been used for data analysis, varying from conventional techniques to novel combinations of advanced chemometric approaches. Steps have been set in making data more findable, accessible, interoperable and reusable, but significant opportunities for improvement remain. Herein, advances in single-cell analysis workflows and data analysis are discussed, and recommendations are made based on the experimental goal.

A critical review of machine-learning for "multi-omics" marine metabolite datasets

During the last decade, genomic, transcriptomic, proteomic, metabolomic, and other omics datasets have been generated for a wide range of marine organisms, and even more are still on the way. Marine organisms possess unique and diverse biosynthetic pathways contributing to the synthesis of novel secondary metabolites with significant bioactivities. As marine organisms have a greater tendency to adapt to stressed environmental conditions, the chance to identify novel bioactive metabolites with potential biotechnological application is very high. This review presents a comprehensive overview of the available "-omics" and "multi-omics" approaches employed for characterizing marine metabolites along with novel data integration tools. The need for the development of machine-learning algorithms for "multi-omics" approaches is briefly discussed. In addition, the challenges involved in the analysis of "multi-omics" data and recommendations for conducting "multi-omics" study were discussed.

Rm-LR: A long-range-based deep learning model for predicting multiple types of RNA modifications

Recent research has highlighted the pivotal role of RNA post-transcriptional modifications in the regulation of RNA expression and function. Accurate identification of RNA modification sites is important for understanding RNA function. In this study, we propose a novel RNA modification prediction method, namely Rm-LR, which leverages a long-range-based deep learning approach to accurately predict multiple types of RNA modifications using RNA sequences only. Rm-LR incorporates two large-scale RNA language pre-trained models to capture discriminative sequential information and learn local important features, which are subsequently integrated through a bilinear attention network. Rm-LR supports a total of ten RNA modification types (m6A, m1A, m5C, m5U, m6Am, Ψ, Am, Cm, Gm, and Um) and significantly outperforms the state-of-the-art methods in terms of predictive capability on benchmark datasets. Experimental results show the effectiveness and superiority of Rm-LR in prediction of various RNA modifications, demonstrating the strong adaptability and robustness of our proposed model. We demonstrate that RNA language pretrained models enable to learn dense biological sequential representations from large-scale long-range RNA corpus, and meanwhile enhance the interpretability of the models. This work contributes to the development of accurate and reliable computational models for RNA modification prediction, providing insights into the complex landscape of RNA modifications.

A diagnostic model for COVID-19 based on proteomics analysis

BACKGROUND AND OBJECTIVE

Early diagnosis of Coronavirus Disease 2019 (COVID-19) can help save patients' lives before the disease turns severe. This can be achieved through an effective and correct treatment protocol. In this paper, a prediction model is proposed to detect infected cases and determine the severity level of the disease.

METHODS

The proposed model is based on utilizing proteins and metabolites as features for each patient, which are then analyzed using feature selection methods such as Principal Component Analysis (PCA), Information Gain (IG), and analysis of Variance (ANOVA) to select the most significant features. The model employs three classifiers, namely K-Nearest Neighbor (KNN), Support Vector Machine (SVM), and Random Forest (RF), to predict and classify the severity level of the COVID-19 infection. The proposed model is evaluated using four performance measures: accuracy, sensitivity, specificity, and precision.

RESULTS

The experiment results show that the proposed model accuracy can reach 80% using RF classifier with PCA. The PCA selects 22 proteins and 10 metabolites. While ANOVA selects 9 proteins and 5 metabolites. The accuracy reaches 92% after applying RF classifier with the ANOVA. Finally, the accuracy reaches 93% using the RF classifier with only ten features. The selected features are 7 proteins and 3 metabolites. Moreover, it shows that the selected features have a relation to the immune system and respiratory systems.

CONCLUSION

The proposed model uses three classifiers and shows promising results by selecting the important features and maximizing the prediction accuracy.

GBA Regulates EMT/MET and Chemoresistance in Squamous Cell Carcinoma Cells by Modulating the Cellular Glycosphingolipid Profile

Glycosphingolipids (GSL) are plasma membrane components that influence molecular processes involved in cancer initiation, progression, and therapeutic responses. They also modulate receptor tyrosine kinases involved in EMT. Therefore, understanding the mechanisms that regulate GSLs in cancer has important therapeutic potential. One critical regulator of GSLs is the lysosomal glucosylceramidase β1 (GBA) that catalyzes the last step in GSL degradation. We show that, in cancer, GBA copy number amplifications and increased expression are widespread. We show that depleting GBA in squamous cell carcinoma cell lines results in a mesenchymal-to-epithelial shift, decreased invasion and migration, increased chemotherapeutic sensitivity, and decreased activation of receptor tyrosine kinases that are involved in regulating EMT. Untargeted lipidomics shows that GBA depletion had significant effects on sphingolipids and GSLs, suggesting that increased GBA activity in cancer sustains EMT and chemoresistance by modulating receptor tyrosine kinase activity and signaling via effects on the cellular lipid profile.

Robust microarray data feature selection using a correntropy based distance metric learning approach

Classification of high-dimensional microarray data is a challenge in bioinformatics and genetic data processing. One of the challenging issues of feature selection is the presence of outliers. The Euclidean distance metric is sensitive to outliers. In this study, a distance metric learning based feature selection approach that uses the correntropy function as the discrimination metric is proposed. For this purpose, the metric learning problem is formulated as an optimization problem and solved using the Lagrange method. The output of the approach signifies the most important and robust features. After feature selection, different classification methods such as SVM, decision trees, and NN classifiers are used to investigate the classification accuracy of the proposed method as well as precision, recall, and F-measure. Experiments are carried out on 13 high-dimensional datasets and show that the proposed method outperforms the previous models in terms of accuracy and robustness.

Advanced network pharmacology study reveals multi-pathway and multi-gene regulatory molecular mechanism of Bacopa monnieri in liver cancer based on data mining, molecular modeling, and microarray data analysis

Liver cancer is a malignant tumor that grows on the surface or inside the liver. The leading cause is a viral infection with hepatitis B or C virus. Natural products and their structural analogues have historically made a major contribution to pharmacotherapy, especially for cancer. A list of studies evidences the therapeutic efficacy of Bacopa monnieri against liver cancer, but the precise molecular mechanism is yet to be discovered. This study combines data mining, network pharmacology, and molecular docking analysis to potentially revolutionize liver cancer treatment by identifying effective phytochemicals. Initially, the information on active constituents of B. monnieri and target genes of both liver cancer and B. monnieri were retrieved from literature as well as from publicly available databases. Based on the matching results between B. monnieri potential targets and liver cancer targets, the protein-protein interaction (PPI) network was constructed using the STRING database and imported into Cytoscape for screening of hub genes based on their degree of connectivity. Later, the interactions network between compounds and overlapping genes was constructed using Cytoscape software to analyze the network pharmacological prospective effects of B. monnieri on liver cancer. Gene Ontology (GO) and KEGG pathway analysis of hub genes revealed that these genes are involved in the cancer-related pathway. Lastly, the expression level of core targets was analyzed using microarray data (GSE39791, GSE76427, GSE22058, GSE87630, and GSE112790). Further, the GEPIA server and PyRx software were used for survival and molecular docking analysis, respectively. In summary, we proposed that quercetin, luteolin, apigenin, catechin, epicatechin, stigmasterol, beta-sitosterol, celastrol, and betulic acid inhibit tumor growth by affecting tumor protein 53 (TP53), interleukin 6 (IL6), RAC-alpha serine/threonine protein kinases 1 (AKT1), caspase-3 (CASP3), tumor necrosis factor (TNF), jun proto-oncogene (JUN), heat shot protein 90 AA1 (HSP90AA1), vascular endothelial growth factor A (VEGFA), epidermal growth factor receptor (EGFR), and SRC proto-oncogene (SRC). Through, microarray data analysis, the expression level of JUN and IL6 were found to be upregulated while the expression level of HSP90AA1 was found to be downregulated. Kaplan-Meier survival analysis indicated that HSP90AA1 and JUN are promising candidate genes that can serve as diagnostic and prognostic biomarkers for liver cancer. Moreover, the molecular docking and molecular dynamic simulation of 60ns well complemented the binding affinity of the compound and revealed strong stability of predicted compounds at the docked site. Calculation of binding free energies using MMPBSA and MMGBSA validated the strong binding affinity between the compound and binding pockets of HSP90AA1 and JUN. Despite that, in vivo and in vitro studies are mandatory to unveil pharmacokinetics and biosafety profiles to completely track the candidature status of B. monnieri in liver cancer.

Multi-omics integration analysis of GPCRs in pan-cancer to uncover inter-omics relationships and potential driver genes

G protein-coupled receptors (GPCRs) are the largest drug target family. Unfortunately, applications of GPCRs in cancer therapy are scarce due to very limited knowledge regarding their correlations with cancers. Multi-omics data enables systematic investigations of GPCRs, yet their effective integration remains a challenge due to the complexity of the data. Here, we adopt two types of integration strategies, multi-staged and meta-dimensional approaches, to fully characterize somatic mutations, somatic copy number alterations (SCNAs), DNA methylations, and mRNA expressions of GPCRs in 33 cancers. Results from the multi-staged integration reveal that GPCR mutations cannot well predict expression dysregulation. The correlations between expressions and SCNAs are primarily positive, while correlations of the methylations with expressions and SCNAs are bimodal with negative correlations predominating. Based on these correlations, 32 and 144 potential cancer-related GPCRs driven by aberrant SCNA and methylation are identified, respectively. In addition, the meta-dimensional integration analysis is carried out by using deep learning models, which predict more than one hundred GPCRs as potential oncogenes. When comparing results between the two integration strategies, 165 cancer-related GPCRs are common in both, suggesting that they should be prioritized in future studies. However, 172 GPCRs emerge in only one, indicating that the two integration strategies should be considered concurrently to complement the information missed by the other such that obtain a more comprehensive understanding. Finally, correlation analysis further reveals that GPCRs, in particular for the class A and adhesion receptors, are generally immune-related. In a whole, the work is for the first time to reveal the associations between different omics layers and highlight the necessity of combing the two strategies in identifying cancer-related GPCRs.

Integrative computational modeling to unravel novel potential biomarkers in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a major health problem around the world. The management of this disease is complicated by the lack of noninvasive diagnostic tools and the few treatment options available. Better clinical outcomes can be achieved if HCC is detected early, but unfortunately, clinical signs appear when the disease is in its late stages. We aim to identify novel genes that can be targeted for the diagnosis and therapy of HCC. We performed a meta-analysis of transcriptomics data to identify differentially expressed genes and applied network analysis to identify hub genes. Fatty acid metabolism, complement and coagulation cascade, chemical carcinogenesis and retinol metabolism were identified as key pathways in HCC. Furthermore, we integrated transcriptomics data into a reference human genome-scale metabolic model to identify key reactions and subsystems relevant in HCC. We conclude that fatty acid activation, purine metabolism, vitamin D, and E metabolism are key processes in the development of HCC and therefore need to be further explored for the development of new therapies. We provide the first evidence that GABRP, HBG1 and DAK (TKFC) genes are important in HCC in humans and warrant further studies.

Clustering ensemble in scRNA-seq data analysis: Methods, applications and challenges

With the rapid development of single-cell RNA-sequencing techniques, various computational methods and tools were proposed to analyze these high-throughput data, which led to an accelerated reveal of potential biological information. As one of the core steps of single-cell transcriptome data analysis, clustering plays a crucial role in identifying cell types and interpreting cellular heterogeneity. However, the results generated by different clustering methods showed distinguishing, and those unstable partitions can affect the accuracy of the analysis to a certain extent. To overcome this challenge and obtain more accurate results, currently clustering ensemble is frequently applied to cluster analysis of single-cell transcriptome datasets, and the results generated by all clustering ensembles are nearly more reliable than those from most of the single clustering partitions. In this review, we summarize applications and challenges of the clustering ensemble method in single-cell transcriptome data analysis, and provide constructive thoughts and references for researchers in this field.

Methodology to identify a gene expression signature by merging microarray datasets

A vast number of microarray datasets have been produced as a way to identify differentially expressed genes and gene expression signatures. A better understanding of these biological processes can help in the diagnosis and prognosis of diseases, as well as in the therapeutic response to drugs. However, most of the available datasets are composed of a reduced number of samples, leading to low statistical, predictive and generalization power. One way to overcome this problem is by merging several microarray datasets into a single dataset, which is typically a challenging task. Statistical methods or supervised machine learning algorithms are usually used to determine gene expression signatures. Nevertheless, statistical methods require an arbitrary threshold to be defined, and supervised machine learning methods can be ineffective when applied to high-dimensional datasets like microarrays. We propose a methodology to identify gene expression signatures by merging microarray datasets. This methodology uses statistical methods to obtain several sets of differentially expressed genes and uses supervised machine learning algorithms to select the gene expression signature. This methodology was validated using two distinct research applications: one using heart failure and the other using autism spectrum disorder microarray datasets. For the first, we obtained a gene expression signature composed of 117 genes, with a classification accuracy of approximately 98%. For the second use case, we obtained a gene expression signature composed of 79 genes, with a classification accuracy of approximately 82%. This methodology was implemented in R language and is available, under the MIT licence, at https://github.com/bioinformatics-ua/MicroGES.

iMSEA: A Novel Metabolite Set Enrichment Analysis Strategy to Decipher Drug Interactions

Drug combinations are commonly used to treat various diseases to achieve synergistic therapeutic effects or to alleviate drug resistance. Nevertheless, some drug combinations might lead to adverse effects, and thus, it is crucial to explore the mechanisms of drug interactions before clinical treatment. Generally, drug interactions have been studied using nonclinical pharmacokinetics, toxicology, and pharmacology. Here, we propose a complementary strategy based on metabolomics, which we call interaction metabolite set enrichment analysis, or iMSEA, to decipher drug interactions. First, a digraph-based heterogeneous network model was constructed to model the biological metabolic network based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Second, treatment-specific influences on all detected metabolites were calculated and propagated across the whole network model. Third, pathway activity was defined and enriched to quantify the influence of each treatment on the predefined functional metabolite sets, i.e., metabolic pathways. Finally, drug interactions were identified by comparing the pathway activity enriched by the drug combination treatments and the single drug treatments. A data set consisting of hepatocellular carcinoma (HCC) cells that were treated with oxaliplatin (OXA) and/or vitamin C (VC) was used to illustrate the effectiveness of the iMSEA strategy for evaluation of drug interactions. Performance evaluation using synthetic noise data was also performed to evaluate sensitivities and parameter settings for the iMSEA strategy. The iMSEA strategy highlighted synergistic effects of combined OXA and VC treatments including the alterations in the glycerophospholipid metabolism pathway and glycine, serine, and threonine metabolism pathway. This work provides an alternative method to reveal the mechanisms of drug combinations from the viewpoint of metabolomics.

Critical Assessment of the Biomarker Discovery and Classification Methods for Multiclass Metabolomics

Multiclass metabolomics has been widely applied in clinical practice to understand pathophysiological processes involved in disease progression and diagnostic biomarkers of various disorders. In contrast to the binary problem, the multiclass classification problem is more difficult in terms of obtaining reliable and stable results due to the increase in the complexity of determining exact class decision boundaries. In particular, methods of biomarker discovery and classification have a significant effect on the multiclass model because different methods with significantly varied theories produce conflicting results even for the same dataset. However, a systematic assessment for selecting the most appropriate methods of biomarker discovery and classification for multiclass metabolomics is still lacking. Therefore, a comprehensive assessment is essential to measure the suitability of methods in multiclass classification models from multiple perspectives. In this study, five biomarker discovery methods and nine classification methods were assessed based on four benchmark datasets of multiclass metabolomics. The performance assessment of the biomarker discovery and classification methods was performed using three evaluation criteria: assessment a (cluster analysis of sample grouping), assessment b (biomarker consistency in multiple subgroups), and assessment c (accuracy in the classification model). As a result, 13 combining strategies with superior performance were selected under multiple criteria based on these benchmark datasets. In conclusion, superior strategies that performed consistently well are suggested for the discovery of biomarkers and the construction of a classification model for multiclass metabolomics.

Metabolomics meets systems immunology

Metabolic processes play a critical role in immune regulation. Metabolomics is the systematic analysis of small molecules (metabolites) in organisms or biological samples, providing an opportunity to comprehensively study interactions between metabolism and immunity in physiology and disease. Integrating metabolomics into systems immunology allows the exploration of the interactions of multilayered features in the biological system and the molecular regulatory mechanism of these features. Here, we provide an overview on recent technological developments of metabolomic applications in immunological research. To begin, two widely used metabolomics approaches are compared: targeted and untargeted metabolomics. Then, we provide a comprehensive overview of the analysis workflow and the computational tools available, including sample preparation, raw spectra data preprocessing, data processing, statistical analysis, and interpretation. Third, we describe how to integrate metabolomics with other omics approaches in immunological studies using available tools. Finally, we discuss new developments in metabolomics and its prospects for immunology research. This review provides guidance to researchers using metabolomics and multiomics in immunity research, thus facilitating the application of systems immunology to disease research.

Identification of key genes in colorectal cancer diagnosis by co-expression analysis weighted gene co-expression network analysis

BACKGROUND

The purpose of this study was using bioinformatics tools to identify biomarkers and molecular factors involved in the diagnosis of colorectal cancer, which are effective for the diagnosis and treatment of the disease.

METHODS

We determined differentially expressed genes (DEGs) related to colorectal cancer (CRC) using the data series retrieved from GEO database. Then the weighted gene co-expression network analysis (WGCNA) was conducted to explore co-expression modules related to CRC diagnosis. Next, the relationship between the integrated modules with clinical features such as the stage of CRC was evaluated. Other downstream analyses were performed on selected module genes.

RESULTS

In this study, after performing the WGCNA method, a module named blue module which was more significantly associated with the CRC stage was selected for further evaluation. Afterward, the Protein-protein interaction network through sting software for 154 genes of the blue module was constructed and eight hub genes were identified through the evaluation of constructed network with Cytoscape. Among these eight hub genes, upregulation of MMP9, SERPINH1, COL1A2, COL5A2, COL1A1, SPARC, and COL5A1 in CRC was validated in other microarray and TCGA data. Based on the results of the mRNA-miRNA interaction network, SERPINH1 was found as a target gene of miR-940. Finally, results of the DGIDB database indicated that Andecaliximab, Carboxylated glucosamine, Marimastat, Tozuleristide, S-3304, Incyclinide, Curcumin, Prinomastat, Demethylwedelolactone, and Bevacizumab, could be used as a therapeutic agent for targeting the MMP9. Furthermore, Ocriplasmin and Collagenase clostridium histolyticum could target COL1A1, COL1A2, COL5A1, and COL5A2.

CONCLUSION

Taken together, the results of the current study indicated that seven hub genes including COL1A2, COL5A1, COL5A2, SERPINH1, MMP9, SPARC, and COL1A1 which were upregulated in CRC could be used as a diagnostic and progression biomarker of CRC. On the other hand, miR-940 which targets SERPINH1 could be used as a potential biomarker of CRC. More ever, Andecaliximab, Carboxylated glucosamine, Marimastat, Tozuleristide, S-3304, Incyclinide, Curcumin, Prinomastat, Demethylwedelolactone, Bevacizumab, Ocriplasmin , and Collagenase clostridium histolyticum were introduced as therapeutic agents for CRC which their therapeutic potential should be evaluated experimentally.

Deep learning automates bidimensional and volumetric tumor burden measurement from MRI in pre- and post-operative glioblastoma patients

Tumor burden assessment by magnetic resonance imaging (MRI) is central to the evaluation of treatment response for glioblastoma. This assessment is, however, complex to perform and associated with high variability due to the high heterogeneity and complexity of the disease. In this work, we tackle this issue and propose a deep learning pipeline for the fully automated end-to-end analysis of glioblastoma patients. Our approach simultaneously identifies tumor sub-regions, including the enhancing tumor, peritumoral edema and surgical cavity in the first step, and then calculates the volumetric and bidimensional measurements that follow the current Response Assessment in Neuro-Oncology (RANO) criteria. Also, we introduce a rigorous manual annotation process which was followed to delineate the tumor sub-regions by the human experts, and to capture their segmentation confidences that are later used while training deep learning models. The results of our extensive experimental study performed over 760 pre-operative and 504 post-operative adult patients with glioma obtained from the public database (acquired at 19 sites in years 2021-2020) and from a clinical treatment trial (47 and 69 sites for pre-/post-operative patients, 2009-2011) and backed up with thorough quantitative, qualitative and statistical analysis revealed that our pipeline performs accurate segmentation of pre- and post-operative MRIs in a fraction of the manual delineation time (up to 20 times faster than humans). Volumetric measurements were in strong agreement with experts with the Intraclass Correlation Coefficient (ICC): 0.959, 0.703, 0.960 for ET, ED, and cavity. Similarly, automated RANO compared favorably with experienced readers (ICC: 0.681 and 0.866) producing consistent and accurate results. Additionally, we showed that RANO measurements are not always sufficient to quantify tumor burden. The high performance of the automated tumor burden measurement highlights the potential of the tool for considerably improving and simplifying radiological evaluation of glioblastoma in clinical trials and clinical practice.

scMAGS: Marker gene selection from scRNA-seq data for spatial transcriptomics studies

Single-Cell RNA sequencing (scRNA-seq) has provided unprecedented opportunities for exploring gene expression and thus uncovering regulatory relationships between genes at the single-cell level. However, scRNA-seq relies on isolating cells from tissues. Therefore, the spatial context of the regulatory processes is lost. A recent technological innovation, spatial transcriptomics, allows for the measurement of gene expression while preserving spatial information. An initial step in the spatial transcriptomic analysis is to identify the cell type, which requires a careful selection of cell-specific marker genes. For this purpose, currently, scRNA-seq data is used to select a limited number of marker genes from among all genes that distinguish cell types from each other. This study proposes scMAGS (single-cell MArker Gene Selection), a novel method for marker gene selection from scRNA-seq data for spatial transcriptomics studies. scMAGS uses a filtering step in which the candidate genes are identified before the marker gene selection step. For the selection of marker genes, cluster validity indices, the Silhouette index, or the Calinski-Harabasz index (for large datasets) are utilized. Experimental results showed that, in comparison to the existing methods, scMAGS is scalable, fast, and accurate. Even for large datasets with millions of cells, scMAGS could find the required number of marker genes in a reasonable amount of time with fewer memory requirements. scMAGS is made freely available at https://github.com/doganlab/scmags and can be downloaded from the Python Package Directory (PyPI) software repository with the command pip install scmags.

GENECI: A novel evolutionary machine learning consensus-based approach for the inference of gene regulatory networks

Gene regulatory networks define the interactions between DNA products and other substances in cells. Increasing knowledge of these networks improves the level of detail with which the processes that trigger different diseases are described and fosters the development of new therapeutic targets. These networks are usually represented by graphs, and the primary sources for their correct construction are usually time series from differential expression data. The inference of networks from this data type has been approached differently in the literature. Mostly, computational learning techniques have been implemented, which have finally shown some specialization in specific datasets. For this reason, the need arises to create new and more robust strategies for reaching a consensus based on previous results to gain a particular capacity for generalization. This paper presents GENECI (GEne NEtwork Consensus Inference), an evolutionary machine learning approach that acts as an organizer for constructing ensembles to process the results of the main inference techniques reported in the literature and to optimize the consensus network derived from them, according to their confidence levels and topological characteristics. After its design, the proposal was confronted with datasets collected from academic benchmarks (DREAM challenges and IRMA network) to quantify its accuracy. Subsequently, it was applied to a real-world biological network of melanoma patients whose results could be contrasted with medical research collected in the literature. Finally, it has been proved that its ability to optimize the consensus of several networks leads to outstanding robustness and accuracy, gaining a certain generalization capacity after facing the inference of multiple datasets. The source code is hosted in a public repository at GitHub under MIT license: https://github.com/AdrianSeguraOrtiz/GENECI. Moreover, to facilitate its installation and use, the software associated with this implementation has been encapsulated in a python package available at PyPI: https://pypi.org/project/geneci/.

Simulated annealing aided genetic algorithm for gene selection from microarray data

In recent times, microarray gene expression datasets have gained significant popularity due to their usefulness to identify different types of cancer directly through bio-markers. These datasets possess a high gene-to-sample ratio and high dimensionality, with only a few genes functioning as bio-markers. Consequently, a significant amount of data is redundant, and it is essential to filter out important genes carefully. In this paper, we propose the Simulated Annealing aided Genetic Algorithm (SAGA), a meta-heuristic approach to identify informative genes from high-dimensional datasets. SAGA utilizes a two-way mutation-based Simulated Annealing (SA) as well as Genetic Algorithm (GA) to ensure a good trade-off between exploitation and exploration of the search space, respectively. The naive version of GA often gets stuck in a local optimum and depends on the initial population, leading to premature convergence. To address this, we have blended a clustering-based population generation with SA to distribute the initial population of GA over the entire feature space. To further enhance the performance, we reduce the initial search space by a score-based filter approach called the Mutually Informed Correlation Coefficient (MICC). The proposed method is evaluated on 6 microarray and 6 omics datasets. Comparison of SAGA with contemporary algorithms has shown that SAGA performs much better than its peers. Our code is available at https://github.com/shyammarjit/SAGA.

Mechanistic insights of Cucumis melo L. seeds for gastrointestinal muscle spasms through calcium signaling pathway-related gene regulation networks in WGCNA and in vitro, in vivo studies

BACKGROUND

In addition to the nutritional benefits of Cucumis melo L., herbalists in Pakistan and India employ seeds to treat various ailments. This study aimed to determine the regulatory role of C. melo seeds in calcium-mediated smooth muscle contraction.

METHODS

We identified and quantified the phytochemicals of C. melo with LC ESI-MS/MS and HPLC, then conducted in vitro and in vivo tests to confirm the involvement in smooth muscle relaxation. Then, diarrhea-predominant irritable bowel syndrome gene datasets from NCBI GEO were acquired, DEGs and WGCNA followed by functional enrichment analysis. Next, molecular docking of key genes was performed.

RESULTS

The quantification of C. melo seeds revealed concentrations of rutin, kaempferol, and quercetin were 702.38 μg/g, 686.29 μg/g, and 658.41 μg/g, respectively. In vitro experiments revealed that C. melo seeds had a dose-dependent relaxant effect for potassium chloride (80 mM)-induced spastic contraction and exhibited calcium antagonistic response in calcium dose-response curves. In in vivo studies, Cm.EtOH exhibited antidiarrheal, antiperistaltic, and antisecretory effects. The functional enrichment of WGCNA and DEGs IBS-associated pathogenic genes, including those involved in calcium-mediated signaling, MAPK cascade, and inflammatory responses. MAPK1 and PIK3CG were identified as key genes with greater binding affinity with rutin, quercitrin, and kaempferol in molecular docking.

CONCLUSIONS

The bronchodilator and antidiarrheal effects of C. melo were produced by altering the regulatory genes of calcium-mediated smooth contraction.

Integrated analysis of multi-omics data for the discovery of biomarkers and therapeutic targets for colorectal cancer

The considerable burden of colorectal cancer and the rising trend in young adults emphasize the necessity of understanding its underlying mechanisms, providing new diagnostic and prognostic markers, and improving therapeutic approaches. Precision medicine is a new trend all over the world and identification of novel biomarkers and therapeutic targets is a step forward towards this trend. In this context, multi-omics data and integrated analysis are being investigated to develop personalized medicine in the management of colorectal cancer. Given the large amount of data from multi-omics approach, data integration and analysis is a great challenge. In this Review, we summarize how statistical and machine learning techniques are applied to analyze multi-omics data and how it contributes to the discovery of useful diagnostic and prognostic biomarkers and therapeutic targets. Moreover, we discuss the importance of these biomarkers and therapeutic targets in the clinical management of colorectal cancer in the future. Taken together, integrated analysis of multi-omics data has great potential for finding novel diagnostic and prognostic biomarkers and therapeutic targets, however, there are still challenges to overcome in future studies.

Computational drug repurposing by exploiting large-scale gene expression data: Strategy, methods and applications

De novo drug development is an extremely complex, time-consuming and costly task. Urgent needs for therapies of various diseases have greatly accelerated searches for more effective drug development methods. Luckily, drug repurposing provides a new and effective perspective on disease treatment. Rapidly increased large-scale transcriptome data paints a detailed prospect of gene expression during disease onset and thus has received wide attention in the field of computational drug repurposing. However, how to efficiently mine transcriptome data and identify new indications for old drugs remains a critical challenge. This review discussed the irreplaceable role of transcriptome data in computational drug repurposing and summarized some representative databases, tools and strategies. More importantly, it proposed a practical guideline through establishing the correspondence between three gene expression data types and five strategies, which would facilitate researchers to adopt appropriate strategies to deeply mine large-scale transcriptome data and discover more effective therapies.

A systematic review of state-of-the-art strategies for machine learning-based protein function prediction

New drug discovery is inseparable from the discovery of drug targets, and the vast majority of the known targets are proteins. At the same time, proteins are essential structural and functional elements of living cells necessary for the maintenance of all forms of life. Therefore, protein functions have become the focus of many pharmacological and biological studies. Traditional experimental techniques are no longer adequate for rapidly growing annotation of protein sequences, and approaches to protein function prediction using computational methods have emerged and flourished. A significant trend has been to use machine learning to achieve this goal. In this review, approaches to protein function prediction based on the sequence, structure, protein-protein interaction (PPI) networks, and fusion of multi-information sources are discussed. The current status of research on protein function prediction using machine learning is considered, and existing challenges and prominent breakthroughs are discussed to provide ideas and methods for future studies.

Investigating cardiotoxicity related with hERG channel blockers using molecular fingerprints and graph attention mechanism

Human ether-a-go-go-related gene (hERG) channel blockade by small molecules is a big concern during drug development in the pharmaceutical industry. Failure or inhibition of hERG channel activity caused by drug molecules can lead to prolonging QT interval, which will result in serious cardiotoxicity. Thus, evaluating the hERG blocking activity of all these small molecular compounds is technically challenging, and the relevant procedures are expensive and time-consuming. In this study, we develop a novel deep learning predictive model named DMFGAM for predicting hERG blockers. In order to characterize the molecule more comprehensively, we first consider the fusion of multiple molecular fingerprint features to characterize its final molecular fingerprint features. Then, we use the multi-head attention mechanism to extract the molecular graph features. Both molecular fingerprint features and molecular graph features are fused as the final features of the compounds to make the feature expression of compounds more comprehensive. Finally, the molecules are classified into hERG blockers or hERG non-blockers through the fully connected neural network. We conduct 5-fold cross-validation experiment to evaluate the performance of DMFGAM, and verify the robustness of DMFGAM on external validation datasets. We believe DMFGAM can serve as a powerful tool to predict hERG channel blockers in the early stages of drug discovery and development.

Innovation in identifying metabolites from complex metabolome-Highlights of recent analytical platforms and protocols

Metabolites are closely intertwined genotypes that can provide clear information about the final phenotype. The high-throughput analysis platform used to identify candidate metabolites and describe their contributions can help to quickly detect metabolic characteristics from large spectral data, which may lead to peak data preprocessing, statistical analysis and functional interpretation. Developing a comprehensive strategy for discovering and verifying bioactive metabolites can provide a large number of new functional biomarkers, and then more closely reveal their functional changes, which has relevant biological significance for disease diagnosis and prognosis treatment.

Label-free proteome quantification and evaluation

The label-free quantification (LFQ) has emerged as an exceptional technique in proteomics owing to its broad proteome coverage, great dynamic ranges and enhanced analytical reproducibility. Due to the extreme difficulty lying in an in-depth quantification, the LFQ chains incorporating a variety of transformation, pretreatment and imputation methods are required and constructed. However, it remains challenging to determine the well-performing chain, owing to its strong dependence on the studied data and the diverse possibility of integrated chains. In this study, an R package EVALFQ was therefore constructed to enable a performance evaluation on >3000 LFQ chains. This package is unique in (a) automatically evaluating the performance using multiple criteria, (b) exploring the quantification accuracy based on spiking proteins and (c) discovering the well-performing chains by comprehensive assessment. All in all, because of its superiority in assessing from multiple perspectives and scanning among over 3000 chains, this package is expected to attract broad interests from the fields of proteomic quantification. The package is available at https://github.com/idrblab/EVALFQ.

Automatic diagnosis and grading of Prostate Cancer with weakly supervised learning on whole slide images

BACKGROUND

The workflow of prostate cancer diagnosis and grading is cumbersome and the results suffer from substantial inter-observer variability. Recent trials have shown potential in using machine learning to develop automated systems to address this challenge. Most automated deep learning systems for prostate cancer Gleason grading focused on supervised learning requiring demanding fine-grained pixel-level annotations.

METHODS

A weakly-supervised deep learning model with slide-level labels is presented in this study for the diagnosis and grading of prostate cancer with whole slide image (WSI). WSIs are first cropped into small patches and then processed with a deep learning model to extract patch-level features. A graph convolution network (GCN) is used to aggregate the features for classifications. Throughout the training process, the noisy labels are progressively filtered out to reduce inter-observer variations in clinical reports. Finally, multi-center independent test cohorts with 6,174 slides are collected to evaluate the prostate cancer diagnosis and grading performance of our model.

RESULTS

The cancer diagnosis (2-level classification) results on two external test sets (n= 4,675, n= 844) show an area under the receiver operating characteristic curve (AUC) of 0.985 and 0.986. The Gleason grading (6-level classification) results reach 0.931 quadratic weighted kappa on the internal test set (n= 531). It generalizes well on the external test dataset (n= 844) with 0.801 quadratic weighted kappa with the reference standard set independently. The model enables pathological meaningful interpretability by visualizing the most attended lesions which are highly consistent with expert annotations.

CONCLUSION

The proposed model incorporates a graph network in weakly supervised learning with only slide-level reports. A robust learning strategy is also employed to correct the label noise. It is highly accurate (>0.985 AUC for diagnosis) and also interpretable with intuitive heatmap visualization. It can be unified with a digital pathology pipeline to deliver prostate cancer metrics for a pathology report.

Network pharmacology and LC-MS approachs to explore the active compounds and mechanisms of Yuanjiang decoction for treating bradyarrhythmia

BACKGROUND

Yuanjiang decoction (YJD), a traditional Chinese medicinal prescription, has been found to have a significant heart rate-increasing effect and is effective in the treatment of symptomatic bradyarrhythmia in previous studies. However, its specific components and potential mechanisms remain unclear.

METHODS

In this study, we detected and identified the main compounds of YJD using liquid chromatography-mass spectrometry (LC-MS). Through the approach of network pharmacology, we predicted the core targets of the active components, bradyarrhythmia targets, and obtained potential anti-bradyarrhythmia targets of YJD. We further performed protein to protein interaction (PPI), gene ontology (GO) enrichment analyses and kyoto encyclopedia of genes and genomes (KEGG) signaling pathway analyses for core targets, and constructed network of key active ingredients-core targets of YJD. Finally, molecular docking and molecular dynamics simulation were performed for key active ingredients and core targets.

RESULTS

The YJD contains a total of 35 main chemical components. The key active ingredients-core targets network contains 36 nodes and 90 edges, including 20 key active ingredients and 16 core targets. The core targets in the PPI network were TP53, TNF, HRAS, PPARG, IL1B, KCNH2, SCN5A, IDH1, LMNA, ACHE, F2, DRD2, CALM1, KCNQ1, TNNI3, IDH2 and TNNT2. KEGG pathway analysis showed that YJD treatment of bradyarrhythmia mainly involves neuroactive ligand-receptor interaction, adrenergic signaling in cardiomyocytes, cAMP signaling pathway, calcium signaling pathway, cholinergic synaptic and serotonergic synapse signaling pathway. The biological processes mainly include regulation of hormone levels, regulation of cardiac contraction, chemical synaptic transmission, circadian rhythm, positive regulation of heart rate, smooth muscle contraction, response to metal ion, oxidation-reduction process, neurotransmitter transport and import across plasma membrane. Molecular docking and molecular dynamics simulation results showed that hesperidin and tetrahydropalmatine had higher affinity with DRD2 and KCNQ1, respectively.

CONCLUSION

This study reveals the pharmacodynamic material basis of YJD and its potential multicomponent-multitarget-multipathway pharmacological effects, predicted its potential anti-bradyarrhythmia mechanism may be related to the regulation of myocardial autonomic nervous function and related ion channels. Our work demonstrates that YJD has great potential for treating bradyarrhythmias as a complementary medicine, and the results can provide a theoretical basis for the development and clinical application of YJD.

A systematic review on the state-of-the-art strategies for protein representation

The study of drug-target protein interaction is a key step in drug research. In recent years, machine learning techniques have become attractive for research, including drug research, due to their automated nature, predictive power, and expected efficiency. Protein representation is a key step in the study of drug-target protein interaction by machine learning, which plays a fundamental role in the ultimate accomplishment of accurate research. With the progress of machine learning, protein representation methods have gradually attracted attention and have consequently developed rapidly. Therefore, in this review, we systematically classify current protein representation methods, comprehensively review them, and discuss the latest advances of interest. According to the information extraction methods and information sources, these representation methods are generally divided into structure and sequence-based representation methods. Each primary class can be further divided into specific subcategories. As for the particular representation methods involve both traditional and the latest approaches. This review contains a comprehensive assessment of the various methods which researchers can use as a reference for their specific protein-related research requirements, including drug research.

Identification of crucial genes for predicting the risk of atherosclerosis with system lupus erythematosus based on comprehensive bioinformatics analysis and machine learning

BACKGROUND

Systemic lupus erythematosus (SLE) has become a major public health problem over the years, and atherosclerosis (AS) is one of the main complications of SLE associated with serious cardiovascular consequences in this patient population. The present study aimed to identify potential biomarkers for SLE patients with AS.

METHODS

Five microarray datasets (GSE50772, GSE81622, GSE100927, GSE28829, GSE37356) were downloaded from the NCBI Gene Expression Omnibus database. The Limma package was used to identify differentially expressed genes (DEGs) in AS. Weighted gene coexpression network analysis (WGCNA) was used to identify significant module genes associated with SLE. Functional enrichment analysis, protein-protein interaction (PPI) network construction, and machine learning algorithms (least absolute shrinkage and selection operator (Lasso, Support Vector Machine-Recursive Feature Elimination (SVM-RFE), and random forest) were applied to identify hub genes. Subsequently, we generated a nomogram and receiver operating characteristic curve (ROC) for predicting the risk of AS in SLE patients. Finally, immune cell infiltrations were analyzed, and Consensus Cluster Analysis was conducted based on Single Sample Gene Set Enrichment Analysis (ssGSEA) scores.

RESULTS

Five hub genes (SPI1, MMP9, C1QA, CX3CR1, and MNDA) were identified and used to establish a nomogram that yielded a high predictive performance (area under the curve 0.900-0.981). Dysregulated immune cell infiltrations were found in AS, with positive correlations with the five hub genes. Consensus clustering showed that the optimal number of subtypes was 3. Compared to subtypes A and B, subtype C presented higher expression of the five hub genes, immune cell infiltration levels and immune checkpoint expression.

CONCLUSION

Our study systematically identified five candidate hub genes (SPI1, MMP9, C1QA, CX3CR1, MNDA) and established a nomogram that could predict the risk of AS with SLE using various bioinformatic analyses and machine learning algorithms. Our findings provide the foothold for future studies on potential crucial genes for AS in SLE patients. Additionally, the dysregulated immune cell proportions and immune checkpoint expressions in AS with SLE were identified.

Therapeutic effects of Crataegus monogyna inhibitors against breast cancer

Breast cancer is a silent killer disorder among women and a serious economic burden in healthcare management. Every 19 s, a woman is diagnosed with breast cancer, and every 74 s, a woman worldwide passes away from the disease. Despite the increase in progressive research, advanced treatment approaches, and preventive measures, breast cancer rates continue to increase. This study provides a combination of data mining, network pharmacology, and docking analysis that surely could revolutionize cancer treatment by exploiting prestigious phytochemicals. Crataegus monogyna is a small, rounded deciduous tree with glossy, deeply lobed leaves and flat sprays of cream flowers, followed by dark red berries in autumn. Various studies demonstrated that C. monogyna is therapeutically effective against breast cancer. However, the particular molecular mechanism is still unknown. This study is credited for locating bioactive substances, metabolic pathways, and target genes for breast cancer treatment. According to the current investigation, which examined compound-target genes-pathway networks, it was found that the bioactive compounds of C. monogyna may operate as a viable solution against breast cancer by altering the target genes implicated in the disease pathogenesis. The expression level of target genes was analyzed using GSE36295 microarray data. Docking analysis and molecular dynamic simulation studies further strengthened the current findings by validating the effective activity of the bioactive compounds against putative target genes. In summary, we propose that six key compounds, luteolin, apigenin, quercetin, kaempferol, ursolic acid, and oleanolic acid, contributed to the development of breast cancer by affecting the MMP9 and PPARG proteins. Integration of network pharmacology and bioinformatics revealed C. monogyna's multitarget pharmacological mechanisms against breast cancer. This study provides convincing evidence that C. monogyna might partially alleviate breast cancer and ultimately lays a foundation for further experimental research on the anti-breast cancer activity of C. monogyna.

Application of Machine Learning in Spatial Proteomics

Spatial proteomics is an interdisciplinary field that investigates the localization and dynamics of proteins, and it has gained extensive attention in recent years, especially the subcellular proteomics. Numerous evidence indicate that the subcellular localization of proteins is associated with various cellular processes and disease progression. Mass spectrometry (MS)-based and imaging-based experimental approaches have been developed to acquire large-scale spatial proteomic data. To allow the reliable analysis of increasingly complex spatial proteomics data, machine learning (ML) methods have been widely used in both MS-based and imaging-based spatial proteomic data analysis pipelines. Here, we comprehensively survey the applications of ML in spatial proteomics from following aspects: (1) data resources for spatial proteome are comprehensively introduced; (2) the roles of different ML algorithms in data analysis pipelines are elaborated; (3) successful applications of spatial proteomics and several analytical tools integrating ML methods are presented; (4) challenges existing in modern ML-based spatial proteomics studies are discussed. This review provides guidelines for researchers seeking to apply ML methods to analyze spatial proteomic data and can facilitate insightful understanding of cell biology as well as the future research in medical and drug discovery communities.

The mechanism and biomarker function of Cavin-2 in lung ischemia-reperfusion injury

BACKGROUND

Lung Ischemia Reperfusion injury(LIRI) is one of the most predominant complications of ischemic lung disease. Cavin-2 emerged as a regulator of a variety of cellular processes, including endocytosis, lipid homeostasis, signal transduction and tumorigenesis, but the function of Cavin-2 in LIRI is unknown. The purpose of this study was to determine the predictive potential of Cavin-2 in protecting lung ischemia-reperfusion injury and its corresponding mechanisms.

METHODS

We found the strong relationship between Cavin-2 and multiple immune-related genes by deep learning method. To reveal the mechanism of Cavin-2 in LIRI, the LIRI SD rat model was constructed to detect the expression of Cavin-2 in the lung tissue of SD rats after LIRI, and the expression of Cavin-2 in lung cell lines was also detected. The expression of IL-6, IL-10 and MDA in cells after Cavin-2 over-expression or knockdown was examined under hypoxic conditions. The expression levels of p-AKT, p-STAT3 and p-ERK1/2 were measured in over-expressing Cavin-2 cells under hypoxic-ischemia conditions, and then the corresponding blockers of AKT, STAT3 and ERK1/2 were given to verify, whether they play a protective role in LIRI.

RESULTS

After hypoxia, the expression of Cavin-2 in rat lung tissues was significantly increased, and the cellular activity and IL-10 in Cavin-2 over-expressing cells were significantly higher than that of the control group, while IL-6 and MDA were significantly lower than that of the control group, while the above results were reversed in Cavin-2 knockdown cells; Meanwhile, the phosphorylation levels of AKT, STAT3, and ERK1/2 were significantly increased in Cavin-2 over-expression cells after hypoxia. When AKT, STAT3, and ERK1/2 specific blockers were given, they lost their protective effect against LIRI.

CONCLUSIONS

Cavin-2 shows biomarker potential in protecting lung from ischemia-reperfusion injury through the survivor activating factor enhancement (SAFE) and reperfusion injury salvage kinase (RISK) pathway.