Prediction of peptide cleavage sites using protein language models and graph neural networks

Prediction of peptide cleavage sites using protein language models and graph neural networks

October 30, 2025

Paula Cifuentes, Ramon Adàlia, Ismael Zamora

Abstract

The growing interest in using peptide molecules as therapeutic agents, driven by their high selectivity and efficacy, has become a significant trend in the pharmaceutical industry. However, their oral administration remains challenging due to their low bioavailability and vulnerability to proteases, which produce the cleavage of peptide bonds. To optimize peptide drug development, in silico tools based on machine learning algorithms have been developed for site of cleavage prediction. These tools, which rely on manual feature extraction, have limitations in capturing complex peptide structures, especially those involving non-natural amino acids or cyclic peptides. This study presents two novel in silico approaches for cleavage site prediction. The first approach uses protein language models, specifically ESM-2, which has been fine- tuned to leverage its learned peptide structure embeddings for accurate cleavage site prediction, eliminating the need for manual feature engineering. The second approach employs graph neural networks, representing peptides via hierarchical graphs at the atom and amino acid levels, effectively handling cyclic peptide structures, including those containing non-natural amino acids. The applicability of this second approach is shown through a case study on a set of four cyclic peptides containing non-natural amino acids, comparing in silico predictions with experimental data.

Scalable Peptide MRM Transition Prediction for High-Throughput Proteomics via Hashing-Based Sequence Encoding

Scalable Peptide MRM Transition Prediction for High-Throughput Proteomics via Hashing-Based Sequence Encoding

Peptide analysis via Multiple Reaction Monitoring (MRM) is indispensable for quantification and/or biomarker validation and drug development, yet its reliance on experimental transition optimization limits scalability. Current computational models for small molecules fail to address peptide-specific complexities, such as sequence-dependent fragmentation and charge-state variability. We introduce a novel framework that combines hashing-based peptide fragment encoding with gradient-boosted decision trees to predict MRM transitions efficiently. This method eliminates bottlenecks in experimental workflows, enabling rapid, resource-efficient transition identification without compromising accuracy—a critical advancement for high-throughput proteomics pipelines.

Molecular Structure and Mass Spectral Data Quality Driven Processing of High-Resolution Mass Spectrometry Data for Pharmacokinetics Studies

Molecular Structure and Mass Spectral Data Quality Driven Processing of High-Resolution Mass Spectrometry Data for Pharmacokinetics Studies

Our inability to comprehensively process high resolution mass spectrometry data for quantitative analysis has long been an impediment to the broader adoption of this powerful technique. We have developed an approach that agnostically and automatically identifies all ions related to the compound in both the MS and MSMS data. The algorithm uses the structure of the molecule to automatically select the optimal compound related MS and MSMS signals, and parameters (extraction window, S/N) to provide the best overall method to meet the assay acceptance criteria defined by the user. Results using this structure and data driven approach are presented for pharmacokinetic data that were collected using the same set of samples analyzed on both QQQ and HRMS instruments.

Plasma lipidomics analysis reveals altered profile of triglycerides and phospholipids in children with Medium-Chain Acyl-CoA dehydrogenase deficiency

Plasma lipidomics analysis reveals altered profile of triglycerides and phospholipids in children with Medium-Chain Acyl-CoA dehydrogenase deficiency

July 2024

Inês M S GuerraHelena B FerreiraTatiana MaurícioMarisa PinhoLuísa DiogoSónia MoreiraLaura GoracciStefano BonciarelliTânia MeloPedro Domingues, M Rosário DominguesAna S P Moreira

Abstract

Medium-chain acyl-CoA dehydrogenase deficiency (MCADD) is the most prevalent mitochondrial fatty acid β-oxidation disorder. In this study, we assessed the variability of the lipid profile in MCADD by analysing plasma samples obtained from 25 children with metabolically controlled MCADD (following a normal diet with frequent feeding and under l-carnitine supplementation) and 21 paediatric control subjects (CT). Gas chromatography-mass spectrometry was employed for the analysis of esterified fatty acids, while high-resolution C18-liquid chromatography-mass spectrometry was used to analyse lipid species. We identified a total of 251 lipid species belonging to 15 distinct lipid classes. Principal component analysis revealed a clear distinction between the MCADD and CT groups. Univariate analysis demonstrated that 126 lipid species exhibited significant differences between the two groups. The lipid species that displayed the most pronounced variations included triacylglycerols and phosphatidylcholines containing saturated and monounsaturated fatty acids, specifically C14:0 and C16:0, which were found to be more abundant in MCADD. The observed changes in the plasma lipidome of children with non-decompensated MCADD suggest an underlying alteration in lipid metabolism. Therefore, longitudinal monitoring and further in-depth investigations are warranted to better understand whether such alterations are specific to MCADD children and their potential long-term impacts.

Keywords: Lipid profile; Lipidomics; Mass spectrometry; Medium‐chain acyl‐CoA dehydrogenase deficiency (MCADD); Phospholipids (PL); Plasma analysis; Triacylglycerols (TG).

Unmasking the lipid landscape: carbamazepine induces alterations in Leydig cell lipidome

Unmasking the lipid landscape: carbamazepine induces alterations in Leydig cell lipidome

February 2025

Inês Nobre, Inês M. S. Guerra, Marisa Pinho, Ana D. Martins, Laura Goracci, Stefano Bonciarelli, Tânia Melo, Pedro Domingues, Artur Paiva, Pedro F. Oliveira, M. Rosário Domingues

Abstract

Leydig cells rely on lipids and fatty acids (FA) for essential functions like maintaining structural integrity, energy metabolism, and steroid hormone synthesis, including testosterone production. Carbamazepine (CBZ), a common anticonvulsant medication, can influence lipid metabolism and profiles, potentially impacting Leydig cell function and testosterone levels. Understanding this interplay is crucial to optimize treatment strategies for individuals requiring CBZ therapy while mitigating any adverse effects on male reproductive health. This study focuses on evaluating the effects of selected CBZ concentrations on the lipid homeostasis of BLTK-1 murine Leydig cells. By employing liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS), we aimed to uncover the specific changes in lipid profiles induced by CBZ exposure (25 and 200 μM). FA analysis demonstrated a significant decrease in FA 22:6 n-3 with increasing CBZ concentration and an increase in the n-6/n-3 ratio. Furthermore, changes in the lipidome, particularly in lipid species belonging to phosphatidylethanolamine (PE), phosphatidylcholine (PC), phosphatidylglycerol (PG), and sphingomyelin (SM) classes were observed. PE and PC lipid species were significantly elevated in Leydig cells exposed to 200 μM CBZ, whereas PG and SM species were downregulated. CBZ treatment significantly altered the Leydig cell phospholipidome, suggesting specific phospholipids such as PG 40:4, PG 34:1, PC O-32:1, PC 32:2, and PE P-38:6, which exhibited the lowest p-values, as potential biomarkers for clinical assessment of CBZ’s impact on Leydig cells. These findings underscore the intricate relationship between CBZ exposure and alterations in lipid profiles, offering potential insights for monitoring and mitigating the drug’s effects on male reproductive health.