MHCpLogics: an interactive machine learning-based tool for unsupervised data visualization and cluster analysis of immunopeptidomes.

The major histocompatibility complex (MHC) encodes a range of immune response genes, including the human leukocyte antigens (HLAs) in humans. These molecules bind peptide antigens and present them on the cell surface for T cell recognition. The repertoires of peptides presented by HLA molecules are termed immunopeptidomes. The highly polymorphic nature of the genres that encode the HLA molecules confers allotype-specific differences in the sequences of bound ligands. Allotype-specific ligand preferences are often defined by peptide-binding motifs. Individuals express up to six classical class I HLA allotypes, which likely present peptides displaying different binding motifs. Such complex datasets make the deconvolution of immunopeptidomic data into allotype-specific contributions and further dissection of binding-specificities challenging. Herein, we developed MHCpLogics as an interactive machine learning-based tool for mining peptide-binding sequence motifs and visualization of immunopeptidome data across complex datasets. We showcase the functionalities of MHCpLogics by analyzing both in-house and published mono- and multi-allelic immunopeptidomics data. The visualization modalities of MHCpLogics allow users to inspect clustered sequences down to individual peptide components and to examine broader sequence patterns within multiple immunopeptidome datasets. MHCpLogics can deconvolute large immunopeptidome datasets enabling the interrogation of clusters for the segregation of allotype-specific peptide sequence motifs, identification of sub-peptidome motifs, and the exportation of clustered peptide sequence lists. The tool facilitates rapid inspection of immunopeptidomes as a resource for the immunology and vaccine communities. MHCpLogics is a standalone application available via an executable installation at: https://github.com/PurcellLab/MHCpLogics.

SAPrIm, a semi-automated protocol for mid-throughput immunopeptidomics.

Human leukocyte antigen (HLA) molecules play a crucial role in directing adaptive immune responses based on the nature of their peptide ligands, collectively coined the immunopeptidome. As such, the study of HLA molecules has been of major interest in the development of cancer immunotherapies such as vaccines and T-cell therapies. Hence, a comprehensive understanding and profiling of the immunopeptidome is required to foster the growth of these personalised solutions. We herein describe SAPrIm, an Immunopeptidomics tool for the Mid-Throughput era. This is a semi-automated workflow involving the KingFisher platform to isolate immunopeptidomes using anti-HLA antibodies coupled to a hyper-porous magnetic protein A microbead, a variable window data independent acquisition (DIA) method and the ability to run up to 12 samples in parallel. Using this workflow, we were able to concordantly identify and quantify ~400 - 13000 unique peptides from 5e5 - 5e7 cells, respectively. Overall, we propose that the application of this workflow will be crucial for the future of immunopeptidome profiling, especially for mid-size cohorts and comparative immunopeptidomics studies.

Immunolyser: A web-based computational pipeline for analysing and mining immunopeptidomic data.

Immunopeptidomics has made tremendous contributions to our understanding of antigen processing and presentation, by identifying and quantifying antigenic peptides presented on the cell surface by Major Histocompatibility Complex (MHC) molecules. Large and complex immunopeptidomics datasets can now be routinely generated using Liquid Chromatography-Mass Spectrometry techniques. The analysis of this data - often consisting of multiple replicates/conditions - rarely follows a standard data processing pipeline, hindering the reproducibility and depth of analysis of immunopeptidomic data. Here, we present Immunolyser, an automated pipeline designed to facilitate computational analysis of immunopeptidomic data with a minimal initial setup. Immunolyser brings together routine analyses, including peptide length distribution, peptide motif analysis, sequence clustering, peptide-MHC binding affinity prediction, and source protein analysis. Immunolyser provides a user-friendly and interactive interface via its webserver and is freely available for academic purposes at https://immunolyser.erc.monash.edu/. The open-access source code can be downloaded at our GitHub repository: https://github.com/prmunday/Immunolyser. We anticipate that Immunolyser will serve as a prominent computational pipeline to facilitate effortless and reproducible analysis of immunopeptidomic data.

Benchmarking Bioinformatics Pipelines in Data-Independent Acquisition Mass Spectrometry for Immunopeptidomics.

Immunopeptidomes are the peptide repertoires bound by the molecules encoded by the major histocompatibility complex [human leukocyte antigen (HLA) in humans]. These HLA-peptide complexes are presented on the cell surface for immune T-cell recognition. Immunopeptidomics denotes the utilization of tandem mass spectrometry to identify and quantify peptides bound to HLA molecules. Data-independent acquisition (DIA) has emerged as a powerful strategy for quantitative proteomics and deep proteome-wide identification; however, DIA application to immunopeptidomics analyses has so far seen limited use. Further, of the many DIA data processing tools currently available, there is no consensus in the immunopeptidomics community on the most appropriate pipeline(s) for in-depth and accurate HLA peptide identification. Herein, we benchmarked four commonly used spectral library-based DIA pipelines developed for proteomics applications (Skyline, Spectronaut, DIA-NN, and PEAKS) for their ability to perform immunopeptidome quantification. We validated and assessed the capability of each tool to identify and quantify HLA-bound peptides. Generally, DIA-NN and PEAKS provided higher immunopeptidome coverage with more reproducible results. Skyline and Spectronaut conferred more accurate peptide identification with lower experimental false-positive rates. All tools demonstrated reasonable correlations in quantifying precursors of HLA-bound peptides. Our benchmarking study suggests a combined strategy of applying at least two complementary DIA software tools to achieve the greatest degree of confidence and in-depth coverage of immunopeptidome data.

FTICR mass spectrometry-based multivariate analysis to explore distinctive metabolites and metabolic pathways: A comprehensive bioanalytical strategy toward time-course metabolic profiling of Thymus vulgaris plants responding to drought stress.

In this research, metabolic profiling/pathways of Thymus vulgaris (thyme) plant were assessed during a water deficit stress using an FTICR mass spectrometry-based metabolomics strategy incorporating multivariate data analysis and bioinformatics techniques. Herein, differences of MS signals in specific time courses after water deficit stress and control cases without any timing period were distinguished significantly by common pattern recognition techniques, i.e., PCA, HCA-Heatmap, and PLS-DA. Subsequently, the results were compared with supervised Kohonen neural network (SKN) ones as a non-linear data visualization and capable mapping tool. The classification models showed excellent performance to predict the level of drought stress. By assessing variances contribution on the PCA-loadings of the MS data, the discriminant variables related to the most critical metabolites were identified and then confirmed by ANOVA. Indeed, FTICR MS-based multivariate analysis strategy could explore distinctive metabolites and metabolic pathways/profiles, grouped into three metabolism categories including amino acids, carbohydrates (i.e., galactose, glucose, fructose, sucrose, and mannose), and other metabolites (rosmarinic acid and citrate), to indicate biological mechanisms in response to drought stress for thyme. It was achieved and approved through the MS signals, genomics databases, and transcriptomics factors to interpret and predict the plant metabolic behavior. Eventually, a comprehensive pathway analysis was used to provide a pathway enrichment analysis and explore topological pathway characteristics dealing with the remarkable metabolites to demonstrate that galactose metabolism is the most significant pathway in the biological system of thyme.

Multivariate spectrochemical analysis of interactions of three common Isatin derivatives to calf thymus DNA in vitro.

Interactions of Isatin and its derivatives, Isatin-3-isonicotinylhydrazone (IINH) and Isatin-ß-thiosemicarbazone (IBT), with calf thymus DNA (ctDNA) have been investigated to delineate pharmaceutical-physicochemical properties using UV-Vis/fluorescence/circular dichroism (CD) spectroscopy, viscosity measurements, and multivariate chemometrics. IINH and IBT molecules intercalate between base pairs of DNA, hypochromism in UV absorptions, increase in the CD positive band, sharp increase in specific viscosity, and the displacement of the methylene blue and Neutral Red dye in complexes with ctDNA, by the IINH and IBT molecules, respectively. The observed intrinsic binding constants (Kb[IBT-ctDNA] = 1.03 × 105 and Kb[IINH-ctDNA] = 1.09 × 105 L mol-1) were roughly comparable to other intercalators. In contrast, Isatin binds with ctDNA via groove mode (Kb[Isatin-ctDNA] = 7.32 × 104 L mol-1) without any significant enhancement in ctDNA viscosity. The fluorescence quenching of Isatin by ctDNA was observed as static. CD spectra indicated that Isatin effectively absorbs into grooves of ctDNA, leading to transition from B to C form. Thermodynamic parameters like enthalpy changes (∆H < 0) and entropy changes (∆S > 0) were calculated according to Van't Hoff's equation, indicating the spontaneous interactions. The common soft/hard chemometric methods were used not only to resolve pure concentration and spectral profiles of components using the acquired spectra but also to calculate Stern-Volmer quenching constants, binding stoichiometry, apparent binding constants (Ka), binding constants (Kb), and thermodynamic parameters. The Kb values for Isatin, IINH, and IBT were calculated as 9.18 × 103, 1.53 × 105, and 2.45 × 104 L mol-1, respectively. The results obtained from experimental-spectroscopic analyses showed acceptable agreement with chemometric outlines.

Oblique rotation of factors: a novel pattern recognition strategy to classify fluorescence excitation-emission matrices of human blood plasma for early diagnosis of colorectal cancer.

Colorectal cancer (CRC) ranks high in both men and women, accounting for about 13% of all cancers. In this study, a novel pattern recognition strategy is proposed to improve early diagnosis of CRC through visualizing the relationship between different spectral patterns in a case-control research. Partial least squares-discriminant analysis (PLS-DA) and supervised Kohonen network (SKN) were used to classify the fluorescence excitation-emission matrices (EEMs) from 289 human blood plasma samples containing CRC patients, adenomas tumor, other non-malignant findings and healthy individuals. To obtain optimal factors, oblique rotation (OR) and genetic algorithm (GA) were used to rotate the factors by optimizing transformation matrix elements. Transformed factors were introduced to SKN to build a classification model and the model performance was examined via comparison with a common classifier; PLS-DA. Classification models were built for CRC-healthy and adenomas-healthy samples and the best results were obtained through applying GA-OR on PLS factors and introducing them to the classifiers. Non-error rates for SKN and PLS-DA models assisted with GA (for selecting more informative PLS factors) and OR were equal to 0.97 and 0.95 in cross validation and 0.93 and 0.90 for prediction of the external test set, respectively. Moreover, according to the acceptable results for adenomas-healthy cases using optimal factors, CRC can be diagnosed in early stages. Combining classifiers and optimal factors proved to be efficient for distinguishing healthy and malignant samples, and OR can significantly improve performance of the classification model.

QSAR analysis for nano-sized layered manganese-calcium oxide in water oxidation: An application of chemometric methods in artificial photosynthesis.

Water oxidation is among the most important reactions in artificial photosynthesis, and nano-sized layered manganese-calcium oxides are efficient catalysts toward this reaction. Herein, a quantitative structure-activity relationship (QSAR) model was constructed to predict the catalytic activities of twenty manganese-calcium oxides toward water oxidation using multiple linear regression (MLR) and genetic algorithm (GA) for multivariate calibration and feature selection, respectively. Although there are eight controlled parameters during synthesizing of the desired catalysts including ripening time, temperature, manganese content, calcium content, potassium content, the ratio of calcium:manganese, the average manganese oxidation state and the surface of catalyst, by using GA only three of them (potassium content, the ratio of calcium:manganese and the average manganese oxidation state) were selected as the most effective parameters on catalytic activities of these compounds. The model's accuracy criteria such as R(2)test and Q(2)test in order to predict catalytic rate for external test set experiments; were equal to 0.941 and 0.906, respectively. Therefore, model reveals acceptable capability to anticipate the catalytic activity.