Prediction Of Retention Times Research Articles

Generic and accurate prediction of retention times in liquid chromatography by post–projection calibration

Retention time predictions from molecule structures in liquid chromatography (LC) are increasingly used in MS–based targeted and untargeted analyses, providing supplementary evidence for molecule annotation and reducing experimental measurements. Nevertheless, different LC setups (e.g., differences in gradient, column, and/or mobile phase) give rise to many prediction models that can only accurately predict retention times for a specific chromatographic method (CM). Here, a generic and accurate method is present to predict retention times across different CMs, by introducing the concept of post–projection calibration. This concept builds on the direct projections of retention times between different CMs and uses 35 external calibrants to eliminate the impact of LC setups on projection accuracy. Results showed that post–projection calibration consistently achieved a median projection error below 3.2% of the elution time. The ranking results of putative candidates reached similar levels among different CMs. This work opens up broad possibilities for coordinating retention times between different laboratories and developing extensive retention databases.

RT-Transformer: retention time prediction for metabolite annotation to assist in metabolite identification.

Liquid chromatography retention times prediction can assist in metabolite identification, which is a critical task and challenge in nontargeted metabolomics. However, different chromatographic conditions may result in different retention times for the same metabolite. Current retention time prediction methods lack sufficient scalability to transfer from one specific chromatographic method to another. Therefore, we present RT-Transformer, a novel deep neural network model coupled with graph attention network and 1D-Transformer, which can predict retention times under any chromatographic methods. First, we obtain a pre-trained model by training RT-Transformer on the large small molecule retention time dataset containing 80038 molecules, and then transfer the resulting model to different chromatographic methods based on transfer learning. When tested on the small molecule retention time dataset, as other authors did, the average absolute error reached 27.30 after removing not retained molecules. Still, it reached 33.41 when no samples were removed. The pre-trained RT-Transformer was further transferred to 5 datasets corresponding to different chromatographic conditions and fine-tuned. According to the experimental results, RT-Transformer achieves competitive performance compared to state-of-the-art methods. In addition, RT-Transformer was applied to 41 external molecular retention time datasets. Extensive evaluations indicate that RT-Transformer has excellent scalability in predicting retention times for liquid chromatography and improves the accuracy of metabolite identification. The source code for the model is available at https://github.com/01dadada/RT-Transformer. The web server is available at https://huggingface.co/spaces/Xue-Jun/RT-Transformer.

Retention modeling of therapeutic peptides in sub‐/supercritical fluid chromatography

AbstractChromatographic modeling software packages are valuable tools during method optimization steps. These are well established for reversed‐phase applications utilizing retention time (RT) prediction to optimize separations and receive robust methods, which is of high interest for the analysis of pharmaceuticals. In contrast to liquid chromatography, the knowledge of RT prediction in supercritical fluid chromatography is limited to a manageable number of studies.This study uses the software DryLab to predict the RTs of the peptides bacitracin (Bac), colistin, tyrothricin (Tyro), and insulin analogs. Gradient time, column temperature, and the ternary composition (terC) of carbon dioxide, methanol (MeOH), and acetonitrile (ACN) in the gradient elution are varied in a feasibility approach using a neutral (Viridis BEH) and an amino‐derivatized aromatic (Torus 2‐PIC) stationary phase. In the second part, chromatographic optimization is performed in silico through gradient adjustments to optimize the separation of the fingerprint of Bac. The final gradient method utilizes a Viridis BEH column (100 × 3.0 mm, 1.7 μm), carbon dioxide, and a modifier consisting of ACN/MeOH/water/methanesulfonic acid (60:40:2:0.1, v:v:v:v). In addition, changes in the retention order of Tyro compounds with the proportion of the terC in combination with a Torus Diol column are investigated.

Retention time prediction for post-translationally modified peptides: Ser, Thr, Tyr-phosphorylation

The development of a peptide retention prediction model for reversed-phase chromatography applications in proteomics is reported for peptides carrying phosphorylated Ser, Thr and Tyr-residues. The major retention features have been assessed using a collection of over 10,000 phosphorylated/non-phosphorylated peptide pairs identified in a series 1D and 2D LC-MS/MS acquisitions using formic acid as ion pairing modifier. Single modification event on average results in increased peptide retention for phosphorylation of Ser (+ 1.46), Thr (+1.33), Tyr (+0.93% acetonitrile, ACN) on gradient elution scale for Luna C18(2) stationary phase. We established several composition and sequence specific features, which drive deviations from these average values. Thus, single phosphorylation of serine results in retention shifts ranging from -2.4 to 5.5% ACN depending on position of the residue, nature of nearest neighbour residues, peptide length, hydrophobicity and pI value, and its propensity to form amphipathic helical structures. We established that the altered ion-pairing environment upon phosphorylation is detrimental for this variability. Hydrophobicity of ion-pairing modifier directly informs the magnitude of expected shifts: (most hydrophilic) 0.5 % acetic acid (larger positive shift upon phosphorylation) > 0.1 % formic acid (positive) > 0.1 % trifluoroacetic (negative) > 0.1 % heptafluorobutyric acid (larger negative shift). The effect of phosphorylation has been also evaluated for several separation conditions used in the first dimension of 2D LC applications: high pH reversed-phase (RP), hydrophilic interaction liquid chromatography (HILIC), strong cation- and strong anion exchange separations.

Structure Based Machine Learning Prediction of Retention Times for LC Method Development of Pharmaceuticals.

Significant resources are spent on developing robust liquid chromatography (LC) methods with optimum conditions for all project in the pipeline. Although, data-driven computer assisted modelling has been implemented to shorten the method development timelines, these modelling approaches require project-specific screening data to model retention time (RT) as function of method parameters. Sometimes method re-development is required, leading to additional investments and redundant laboratory work. Cheminformatics techniques have been successfully used to predict the RT of metabolites & other component mixtures for similar use cases. Here we will show that these techniques can be used to model structurally diverse molecules and predictions of these models trained on multiple LC conditions can be used for downstream data-driven modelling. The Molecular Operating Environment(MOE) was used to calculate over 800 descriptors using the strucutres of the analytes. These descriptors were used to model the RT of the analytes under four chromatographic conditions. These models were then used to create data-driven models using LC-SIM. A structural-based Random Forest (RF) model outperformed other techniques in cross-validation studies and predicted the RTs of a randomized test set with a median percentage error less than 4% for all LC conditions. RTs predicted by this structure-based model were used to fit a data-driven model that identifies optimum LC conditions without any additional experimental work. These results show that small training sets yield pharmaceutically relevant models when used in a combination of structure-based and data-driven model.

MARS an improved de novo peptide candidate selection method for non-canonical antigen target discovery in cancer

Understanding the nature and extent of non-canonical human leukocyte antigen (HLA) presentation in tumour cells is a priority for target antigen discovery for the development of next generation immunotherapies in cancer. We here employ a de novo mass spectrometric sequencing approach with a refined, MHC-centric analysis strategy to detect non-canonical MHC-associated peptides specific to cancer without any prior knowledge of the target sequence from genomic or RNA sequencing data. Our strategy integrates MHC binding rank, Average local confidence scores, and peptide Retention time prediction for improved de novo candidate Selection; culminating in the machine learning model MARS. We benchmark our model on a large synthetic peptide library dataset and reanalysis of a published dataset of high-quality non-canonical MHC-associated peptide identifications in human cancer. We achieve almost 2-fold improvement for high quality spectral assignments in comparison to de novo sequencing alone with an estimated accuracy of above 85.7% when integrated with a stepwise peptide sequence mapping strategy. Finally, we utilize MARS to detect and validate lncRNA-derived peptides in human cervical tumour resections, demonstrating its suitability to discover novel, immunogenic, non-canonical peptide sequences in primary tumour tissue.

Simulation of gas chromatographic separations and estimation of distribution-centric retention parameters using linear solvation energy relationships

For method development in gas chromatography, suitable computer simulations can be very helpful during the optimization process. For such computer simulations retention parameters are needed, that describe the interaction of the analytes with the stationary phase during the separation process. There are different approaches to describe such an interaction, e.g. thermodynamic models like Blumberg's distribution-centric 3-parameter model (K-centric model) or models using chemical properties like the Linear Solvation Energy Relationships (LSER).In this work LSER models for a Rxi-17Sil MS and a Rxi-5Sil MS GC column are developed for different temperatures. The influences of the temperature to the LSER system coefficients are shown in a range between 40 and 200 °C and can be described with Clark and Glew's ABC model as fit function. A thermodynamic interpretation of the system constants is given and its contribution to enthalpy and entropy is calculated. An estimation method for the retention parameters of the K-centric model via LSER models were presented. The predicted retention parameters for a selection of 172 various compounds, such as FAMEs, PCBs and PAHs are compared to isothermal determined values. 40 measurements of temperature programmed GC separations are compared to computer simulations using the differently determined or estimated K-centric retention parameters. The mean difference (RSME) between the measured and predicted retention time is less than 8 s for both stationary phases using the isothermal retention parameters. With the LSER predicted parameters the difference is 20 s for the Rxi-5Sil MS and 38 s for the Rxi-17Sil MS. Therefore, the presented estimation method can be recommended for first method development in gas chromatography.

An improved workflow for the development of MS-compatible liquid chromatography assay purity and purification methods by using automated LC Screening instrumentation and in silico modeling.

The development of liquid chromatography UV and mass spectrometry (LC-UV-MS) assays in pharmaceutical analysis is pivotal to improve quality control by providing critical information about drug purity, stability, and presence and identity of byproducts and impurities. Analytical method development of these assays is time-consuming, which often causes it to become a bottle neck in drug development and poses a challenge for process chemists to quickly improve the chemistry. In this study, a systematic and efficient workflow was designed to develop purity assay and purification methods for a wide range of compounds including peptides, proteins, and small molecules with MS-compatible mobile phases (MP) by using automated LC screening instrumentation and in silico modeling tools. Initial LC MPs and chromatography column screening experiments enabled quick identification of conditions which provided the best resolution in the vicinity of the target compounds, which is further optimized using computer-assisted modeling (LC Simulator from ACD/Labs). The experimental retention times were in good agreement with the predicted retention times from LC Simulator (ΔtR < 7%). This workflow presents a practical workflow to significantly expedite the time needed to develop optimized LC-UV-MS methods, allowing for a facile, automatic method optimization and reducing the amount of manual work involved in developing new methods during drug development.

Retention time prediction and MRM validation reinforce the biomarker identification of LC-MS based phospholipidomics.

Dysfunctional lipid metabolism plays a crucial role in the development and progression of various diseases. Accurate measurement of lipidomes can help uncover the complex interactions between genes, proteins, and lipids in health and diseases. The prediction of retention time (RT) has become increasingly important in both targeted and untargeted metabolomics. However, the potential impact of RT prediction on targeted LC-MS based lipidomics is still not fully understood. Herein, we propose a simplified workflow for predicting RT in phospholipidomics. Our approach involves utilizing the fatty acyl chain length or carbon-carbon double bond (DB) number in combination with multiple reaction monitoring (MRM) validation. We found that our model's predictive capacity for RT was comparable to that of a publicly accessible program (QSRR Automator). Additionally, MRM validation helped in further mitigating the interference in signal recognition. Using this developed workflow, we conducted phospholipidomics of sorafenib resistant hepatocellular carcinoma (HCC) cell lines, namely MHCC97H and Hep3B. Our findings revealed an abundance of monounsaturated fatty acyl (MUFA) or polyunsaturated fatty acyl (PUFA) phospholipids in these cell lines after developing drug resistance. In both cell lines, a total of 29 lipids were found to be co-upregulated and 5 lipids were co-downregulated. Further validation was conducted on seven of the upregulated lipids using an independent dataset, which demonstrates the potential for translation of the established workflow or the lipid biomarkers.

Deep Learning-Assisted Analysis of Immunopeptidomics Data.

Liquid chromatography-coupled mass spectrometry (LC-MS/MS) is the primary method to obtain direct evidence for the presentation of disease- or patient-specific human leukocyte antigen (HLA). However, compared to the analysis of tryptic peptides in proteomics, the analysis of HLA peptides still poses computational and statistical challenges. Recently, fragment ion intensity-based matching scores assessing the similarity between predicted and observed spectra were shown to substantially increase the number of confidently identified peptides, particularly in use cases where non-tryptic peptides are analyzed. In this chapter, we describe in detail three procedures on how to benefit from state-of-the-art deep learning models to analyze and validate single spectra, single measurements, and multiple measurements in mass spectrometry-based immunopeptidomics. For this, we explain how to use the Universal Spectrum Explorer (USE), online Oktoberfest, and offline Oktoberfest. For intensity-based scoring, Oktoberfest uses fragment ion intensity and retention time predictions from the deep learning framework Prosit, a deep neural network trained on a very large number of synthetic peptides and tandem mass spectra generated within the ProteomeTools project. The examples shown highlight how deep learning-assisted analysis can increase the number of identified HLA peptides, facilitate the discovery of confidently identified neo-epitopes, or provide assistance in the assessment of the presence of cryptic peptides, such as spliced peptides.

Deep Learning Prediction Boosts Phosphoproteomics-Based Discoveries Through Improved Phosphopeptide Identification

Shotgun phosphoproteomics enables high-throughput analysis of phosphopeptides in biological samples. One of the primary challenges associated with this technology is the relatively low rate of phosphopeptide identification during data analysis. This limitation hampers the full realization of the potential offered by shotgun phosphoproteomics. Here we present DeepRescore2, a computational workflow that leverages deep learning-based retention time and fragment ion intensity predictions to improve phosphopeptide identification and phosphosite localization. Using a state-of-the-art computational workflow as a benchmark, DeepRescore2 increases the number of correctly identified peptide-spectrum matches by 17% in a synthetic dataset and identifies 19% to 46% more phosphopeptides in biological datasets. In a liver cancer dataset, 30% of the significantly altered phosphosites between tumor and normal tissues and 60% of the prognosis-associated phosphosites identified from DeepRescore2-processed data could not be identified based on the state-of-the-art workflow. Notably, DeepRescore2-processed data uniquely identifies EGFR hyperactivation as a new target in poor-prognosis liver cancer, which is validated experimentally. Integration of deep learning prediction in DeepRescore2 improves phosphopeptide identification and facilitates biological discoveries.

Machine learning prediction for constructing a universal multidimensional information library of Panax saponins (ginsenosides)

Accurate characterization of Panax herb ginsenosides is challenging because of the isomers and lack of sufficient reference compounds. More structural information could help differentiate ginsenosides and their isomers, enabling more accurate identification. Based on the VionTM ion-mobility high-resolution LC-MS platform, a multidimensional information library for ginsenosides, namely GinMIL, was established by predicting retention time (tR) and collision cross section (CCS) through machine learning. Robustness validation experiments proved tR and CCS were suitable for database construction. Among three machine learning models we attempted, gradient boosting machine (GBM) exhibited the best prediction performance. GinMIL included the multidimensional information (m/z, molecular formula, tR, CCS, and some MS/MS fragments) for 579 known ginsenosides. Accuracy in identifying ginsenosides from diverse ginseng products was greatly improved by a unique LC-MS approach and searching GinMIL, demonstrating a universal Panax saponins library constructed based on hierarchical design. GinMIL could improve the accuracy of isomers identification by approximately 88%.

Retention Time Prediction through Learning from a Small Training Data Set with a Pretrained Graph Neural Network.

Graph neural networks (GNNs) have shown remarkable performance in predicting the retention time (RT) for small molecules. However, the training data set for a particular target chromatographic system tends to exhibit scarcity, which poses a challenge because the experimental process for measuring RT is costly. To address this challenge, transfer learning has been used to leverage an abundant training data set from a related source task. In this study, we present an improved transfer learning method to better predict the RT of molecules for a target chromatographic system by learning from a small training data set with a pretrained GNN. We use a graph isomorphism network as the architecture of the GNN. The GNN is pretrained on the METLIN-SMRT data set and is then fine-tuned on the target training data set for a fixed number of training iterations using the limited-memory Broyden-Fletcher-Goldfarb-Shanno optimizer with a learning rate decay. We demonstrate that the proposed method achieves superior predictive performance on various chromatographic systems compared with that of the existing transfer learning methods, especially when only a small training data set is available for use. A potential avenue for future research is to leverage multiple small training data sets from different chromatographic systems to further enhance the generalization performance.

The Applicability of Chromatographic Retention Modeling on Chiral Stationary Phases in Reverse-Phase Mode: A Case Study for Ezetimibe and Its Impurities.

Mechanistic modeling is useful for predicting and modulating selectivity even in early chromatographic method development. This approach is also in accordance with current analytical quality using design principles and is highly welcomed by the authorities. The aim of this study was to investigate the separation behavior of two different types of chiral stationary phases (CSPs) for the separation of ezetimibe and its related substances using the mechanistic retention modeling approach offered by the Drylab software (version 4.5) package. Based on the obtained results, both CSPs presented with chemoselectivity towards the impurities of ezetimibe. The cyclodextrin-based CSP displayed a higher separation capacity and was able to separate seven related substances from the active pharmaceutical ingredient, while the cellulose-based column enabled the baseline resolution of six impurities from ezetimibe. Generally, the accuracy of predicted retention times was lower for the polysaccharide CSP, which could indicate the presence of additional secondary interactions between the analytes and the CSP. It was also demonstrated that the combination of mechanistic modeling and an experimental design approach can be applied to method development on CSPs in reverse-phase mode. The applicability of the methods was tested on spiked artificial placebo samples, while intraday and long-term (2 years) method repeatability was also challenged through comparing the obtained retention times and resolution values. The results indicated the excellent robustness of the selected setpoints. Overall, our findings indicate that the chiral columns could offer orthogonal selectivity to traditional reverse-phase columns for the separation of structurally similar compounds.

High-throughput multitarget quantitative assay to profile the whole grain-specific phytochemicals alkylresorcinols, benzoxazinoids and avenanthramides in whole grain and grain-based foods.

Currently, there is a growing interest in using whole grain (WG)-specific phytochemicals to perform WG research, including research on dietary assessment, health mechanisms, and quality control. However, the current approaches used for WG-specific phytochemical analysis cannot simultaneously achieve coverage, specificity, and sensitivity. In the present study, a series of WG-specific phytochemicals (alkylresorcinols (ARs), benzoxazinoids (BXs) and avenanthramides (AVAs)) were identified, and their mass spectrometry (MS) fragmentation mechanism was studied by TOF MS. Based on diagnostic fragmentation ions and retention time prediction models, a LC-MS/MS method was developed. Through this method, 56 ARs, 13 BXs, and 19 AVAs in WGs and grain-based foods were quantified for the first time. This method was validated and yielded excellent specificity, high sensitivity and negligible matrix effects. Finally, we established WG-specific phytochemical fingerprints in a variety of WG and grain-based foods. This method can be used for WG quality control and WG precision nutrition research.

Global Analysis of 2-Hydroxy Fatty Acids by Gas Chromatography-Tandem Mass Spectrometry Reveals Species-Specific Enrichment in Echinoderms.

Abnormal levels of 2-hydroxy fatty acids (2-OH FAs) are characterized in multiple diseases, and their quantification in foodstuffs is critical to identify the sources of supplementation for potential treatment. However, due to the structural complexity and limited available standards, the comprehensive profiling of 2-OH FAs remains an ongoing challenge. Herein, an innovative approach based on gas chromatography-tandem mass spectrometry (GC-MS/MS) was developed to determine the full profile of these FA metabolites. MS and MS/MS spectra of the trimethylsilyl (TMS) derivatives of 2-OH fatty acid methyl esters (FAMEs) were collected for peak annotation by their signature fragmentation patterns. The structures were further confirmed by validated structure-dependent retention time (RT) prediction models, taking advantage of the correlation between the RT, carbon chain length, and double bond number from commercial standards and pseudostandards identified in the whole-brain samples from mice. An in-house database containing 50 saturated and monounsaturated 2-OH FAs was established, which is expandible when additional molecular species with different chain lengths and backbone structures are identified in the future. A quantitation method was then developed by scheduled multiple reaction monitoring (MRM) and applied to investigate the profiling of 2-OH FAs in echinoderms. Our results revealed that the levels of total 2-OH FAs in sea cucumber Apostichopus japonicas (8.40 ± 0.28 mg/g dry weight) and starfish Asterias amurensis (7.51 ± 0.18 mg/g dry weight) are much higher than that in sea urchin Mesocentrotus nudus (531 ± 108 μg/g dry weight). Moreover, 2-OH C24:1 is the predominant molecular species accounting for 67.9% of the total 2-OH FA in sea cucumber, while 2-OH C16:0 is the major molecular species in starfish. In conclusion, the current innovative GC-MS approach has successfully characterized distinct molecular species of 2-OH FAs that are highly present in sea cucumbers and starfish. Thus, these findings suggest the possibility of developing future feeding strategies for preventing and treating diseases associated with 2-OH FA deficiency.

Simultaneously quantifying hundreds of acylcarnitines in multiple biological matrices within ten minutes using ultrahigh-performance liquid-chromatography and tandem mass spectrometry

Acylcarnitines are metabolic intermediates of fatty acids and branched-chain amino acids having vital biofunctions and pathophysiological significances. Here, we developed a high-throughput method for quantifying hundreds of acylcarnitines in one run using ultrahigh performance liquid chromatography and tandem mass spectrometry (UPLC-MS/MS). This enabled simultaneous quantification of 1136 acylcarnitines (C0–C26) within 10-min with good sensitivity (limit of detection < 0.7 fmol), linearity (correlation coefficient > 0.992), accuracy (relative error < 20%), precision (coefficient of variation (CV), CV < 15%), stability (CV < 15%), and inter-technician consistency (CV < 20%, n = 6). We also established a quantitative structure-retention relationship (goodness of fit > 0.998) for predicting retention time (tR) of acylcarnitines with no standards and built a database of their multiple reaction monitoring parameters (tR, ion-pairs, and collision energy). Furthermore, we quantified 514 acylcarnitines in human plasma and urine, mouse kidney, liver, heart, lung, and muscle. This provides a rapid method for quantifying acylcarnitines in multiple biological matrices.

Deep graph convolutional network for small-molecule retention time prediction

The retention time (RT) is a crucial source of data for liquid chromatography-mass spectrometry (LCMS). A model that can accurately predict the RT for each molecule would empower filtering candidates with similar spectra but differing RT in LCMS-based molecule identification. Recent research shows that graph neural networks (GNNs) outperform traditional machine learning algorithms in RT prediction. However, all of these models use relatively shallow GNNs. This study for the first time investigates how depth affects GNNs' performance on RT prediction. The results demonstrate that a notable improvement can be achieved by pushing the depth of GNNs to 16 layers by the adoption of residual connection. Additionally, we also find that graph convolutional network (GCN) model benefits from the edge information. The developed deep graph convolutional network, DeepGCN-RT, significantly outperforms the previous state-of-the-art method and achieves the lowest mean absolute percentage error (MAPE) of 3.3% and the lowest mean absolute error (MAE) of 26.55 s on the SMRT test set. We also finetune DeepGCN-RT on seven datasets with various chromatographic conditions. The mean MAE of the seven datasets largely decreases 30% compared to previous state-of-the-art method. On the RIKEN-PlaSMA dataset, we also test the effectiveness of DeepGCN-RT in assisting molecular structure identification. By 30% lessening the number of potential structures, DeepGCN-RT is able to improve top-1 accuracy by about 11%.

Suspect and nontarget screening of per- and polyfluoroalkyl substances based on ion mobility mass spectrometry and machine learning techniques

High-resolution mass spectrometry (HRMS)-based suspect and nontarget screening techniques are powerful tools for the comprehensive identification of per- and polyfluoroalkyl substances (PFASs), but the interference of complex matrices (especially for wastewater) pose an analytical challenge. This study explored the potential of combining ion mobility spectrometry (IMS) with HRMS and machine learning techniques to achieve the rapid and accurate suspect and nontarget screening of PFAS in wastewater. There were fewer interfering peaks and a clearer spectrum in the data acquired by IMS-HRMS than conventional HRMS. The introduction of collision cross section (CCS) in PFAS homologous series search could filter out 63% of false positive results. Retention time and CCS prediction models were helpful in improving the confidence for PFAS qualitative identification and the random forest algorithm combined with RDKit descriptor performed best for CCS prediction. With the inclusion of extra dimensional information, this study also proposed a comprehensive and concise confidence assignment criterion to better convey the certainty of the qualitative identification of PFAS. Finally, a total of 56 potential PFASs were identified in the wastewater sample using the newly developed method and 45 of them were identified outside reference standards, emphasizing the importance of suspect and nontarget screening for PFAS.

Liquid chromatographic retention time prediction models to secure and improve the feature annotation process in high-resolution mass spectrometry

The development of quantitative structure-retention relationship (QSRR) models has, until recently, required an adequate selection of molecular descriptors necessarily obtained based on a known chemical structure. However, these complex descriptors are not always available nor calculable when the high-resolution mass spectrometry (HRMS) annotation process is underway. Depending on the level of annotation, many structures or even various molecular formulas could be candidates. To secure and improve the annotation process and to save time, a QSRR model (using only 0D molecular descriptors) to predict retention times in reverse-phase liquid chromatography (RPLC) based on the molecular formula was developed, and a general QSRR annotation-based methodology was also proposed.