Quantum chemical mass spectrometry: ab initio prediction of electron ionization mass spectra and identification of new fragmentation pathways

Computational methods for processing and interpreting mass spectrometry-based metabolomics

Metabolomics has emerged as an indispensable tool for exploring complex biological questions, providing the ability to investigate a substantial portion of the metabolome. However, the vast complexity and structural diversity intrinsic to metabolites imposes a great challenge for data analysis and interpretation. Liquid chromatography mass spectrometry (LC-MS) stands out as a versatile technique offering extensive metabolite coverage. In this mini-review, we address some of the hurdles posed by the complex nature of LC-MS data, providing a brief overview of computational tools designed to help tackling these challenges. Our focus centers on two major steps that are essential to most metabolomics investigations: the translation of raw data into quantifiable features, and the extraction of structural insights from mass spectra to facilitate metabolite identification. By exploring current computational solutions, we aim at providing a critical overview of the capabilities and constraints of mass spectrometry-based metabolomics, while introduce some of the most recent trends in data processing and analysis within the field.

Computed Mass-Fragmentation Energy Profiles of Some Acetalized Monosaccharides for Identification in Mass Spectrometry

Our study found that quantum calculations can differentiate fragmentation energies into isomeric structures with asymmetric carbon atoms, such as those of acetalized monosaccharides. It was justified by the good results that have been published in recent years on the discrimination of structural isomers and diastereomers by correlating the calculated mass energy fragmentation profiles with their mass spectra. Based on the quantitative structure–fragmentation relationship (QSFR), this technique compares the intensities of primary ions from the experimental spectrum using the mass energy profiles calculated for the candidate structures. Maximum fit is obtained for the true structure. For a preliminary assessment of the accuracy of the identification of some di-O-isopropylidene monosaccharide diastereomers, we used fragmentation enthalpies (ΔfH) and Gibbs energies (ΔfG) as the energetic descriptors of fragmentation. Four quantum chemical methods were used: RM1, PM7, DFT ΔfH and DFT ΔfG. The mass energy database shows that the differences between the profiles of the isomeric candidate structures could be large enough to be distinguished from each other. This database allows the optimization of energy descriptors and quantum computing methods that can ensure the correct identification of these isomers.

Simulation of the electron ionization mass spectra of the Novichok nerve agent

Novichok is one of the most feared and controversial nerve agents, which existence was confirmed only after the Salisbury attack in 2018. A new attack on August 2020, in Russia, was confirmed. After the 2018 attack, the agent was included in the list of the most dangerous chemicals of the Chemical Weapons Convention (CWC). However, information related to its electron ionization mass spectrometry (EI/MS), essential for unambiguous identification, is scarce. Therefore, investigations about Novichok EI/MS are urgent. In this work, we employed Born-Oppenheimer molecular dynamics through the Quantum Chemistry Electron Ionization Mass Spectrometry (QCEIMS) method to simulate and rationalize the EI/MS spectra and fragmentation pathways of 32 Novichok molecules recently incorporated into the CWC. The comparison of additional simulations with the measured EI spectrum of another Novichok analog is very favorable. A general scheme of the fragmentation pathways derived from simulation results was presented. The present results will be useful for elucidation and prediction of the EI spectra and fragmentation pathways of the dangerous Novichok nerve agent.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

A map of mass spectrometry-based in silico fragmentation prediction and compound identification in metabolomics

Metabolomics, the comprehensive study of the metabolome, and lipidomics-the large-scale study of pathways and networks of cellular lipids-are major driving forces in enabling personalized medicine. Complicated and error-prone data analysis still remains a bottleneck, however, especially for identifying novel metabolites. Comparing experimental mass spectra to curated databases containing reference spectra has been the gold standard for identification of compounds, but constructing such databases is a costly and time-demanding task. Many software applications try to circumvent this process by utilizing cutting-edge advances in computational methods-including quantum chemistry and machine learning-and simulate mass spectra by performing theoretical, so called in silico fragmentations of compounds. Other solutions concentrate directly on experimental spectra and try to identify structural properties by investigating reoccurring patterns and the relationships between them. The considerable progress made in the field allows recent approaches to provide valuable clues to expedite annotation of experimental mass spectra. This review sheds light on individual strengths and weaknesses of these tools, and attempts to evaluate them-especially in view of lipidomics, when considering complex mixtures found in biological samples as well as mass spectrometer inter-instrument variability.

QCMS2 as a new method for providing insight into peptide fragmentation: The influence of the side-chain and inter-side-chain interactions

The identification of peptides and proteins from tandem mass spectra is a difficult task and multiple tools have been developed to aid this identification. We present a new method called quantum chemical mass spectrometry for materials science (QCMS² ), which is based on quantum chemical calculations of bond orders, reaction, and transition-state energies at the DFT/B3LYP/6-311+G* level of theory. The method was used to describe the fragmentation pathways of five X-His-Ser tripeptides with X = Asn, Asp, Glu, Ser, and Trp, thereby focusing on the influence of the side chain and inter-side-chain interactions on the fragmentation. The main features in the mass spectra of the five tripeptides were correctly reproduced, and a number of fragments were assigned to fragmentations involving the side chain and the influence of inter-side-chain interactions. Product ion spectra were recorded to evaluate the capabilities and limitations of QCMS² and a number of conventional tools.

Quantum chemical mass spectrometry: Ab initio study of b2 -ion formation mechanisms for the singly protonated Gln-His-Ser tripeptide

RATIONALE

Both amide bond protonation triggering peptide fragmentations and the controversial b₂ -ion structures have been subjects of intense research. The involvement of histidine (H), with its imidazole side chain that induces specific dissociation patterns involving inter-side-chain (ISC) interactions, in b₂ -ion formation was investigated, focusing on the QHS model tripeptide.

METHODS

To identify the effect of histidine on fragmentations issued from ISC interactions, QHS was selected for a comprehensive analysis of the pathways leading to the three possible b₂ -ion structures, using quantum chemical calculations performed at the DFT/B3LYP/6-311+G* level of theory. Electrospray ionization ion trap mass spectrometry allowed the recording of MS² and MS³ tandem mass spectra, whereas the Quantum Chemical Mass Spectrometry for Materials Science (QCMS² ) method was used to predict fragmentation patterns.

RESULTS

Whereas it is very difficult to differentiate among protonated oxazolone, diketopiperazine, or lactam b₂ -ions using MS² and MS³ mass spectra, the calculations indicated that the QH b₂ -ion (detected at m/z 266) is probably a mixture of the lactam and oxazolone structures formed after amide nitrogen protonation, making the formation of diketopiperazine less likely as it requires an additional step for its formation.

CONCLUSIONS

In contrast to glycine-histidine-containing b₂ -ions, known to be issued from the backbone-imidazole cyclization, we found that interactions between the side chains were not obvious to perceive, neither from a thermodynamics nor from a fragmentation perspective, emphasizing the importance of the whole sequence on the dissociation behavior usually demonstrated from simple glycine-containing tripeptides.

Molecular dynamics simulation of the electron ionization mass spectrum of tabun

Tabun (ethyl N,N-dimethylphosphoramidocyanidate), or GA, is a chemical warfare nerve agent produced during the World War II. The synthesis of its analogs is rather simple; thus, it is a significant threat. Furthermore, experiments with tabun and other nerve agents are greatly limited by the involved life risks and the severe restrictions imposed by the Chemical Weapons Convention. For these reasons, accurate theoretical assignment of fragmentation pathways can be especially important. In this work, we employ the Quantum Chemistry Electron Ionization Mass Spectra method, which combines molecular dynamics, quantum chemistry methods, and stochastic approaches, to accurately investigate the electron ionization/mass spectrometry (EI/MS) fragmentation spectrum and pathways of the tabun molecule. We found that different rearrangement reactions occur including a McLafferty involving the nitrile group. An essential and characteristic pathway for identification of tabun and analogs, a two-step fragmentation producing the m/z 70 ion, was confirmed. The present results will be also useful to predict EI/MS spectrum and fragmentation pathways of other members of the tabun family, namely, the O-alkyl/cycloalkyl N,N-dialkyl (methyl, ethyl, isopropyl, or propyl) phosphoramidocyanidates.

An overview of tools, software, and methods for natural product fragment and mass spectral analysis

One major challenge in natural product (NP) discovery is the determination of the chemical structure of unknown metabolites using automated software tools from either GC–mass spectrometry (MS) or liquid chromatography–MS/MS data only. This chapter reviews the existing spectral libraries and predictive computational tools used in MS-based untargeted metabolomics, which is currently a hot topic in NP structure elucidation. We begin by focusing on spectral databases and the general workflow of MS annotation. We then describe software and tools used in MS, particularly those used to predict fragmentation patterns, mass spectral classifiers, and tools for fragmentation trees analysis. We then round up the chapter by looking at more advanced approaches implemented in tools for competitive fragmentation modeling and quantum chemical approaches.

The use of fast molecular descriptors and artificial neural networks approach in organochlorine compounds electron ionization mass spectra classification

Developing of theoretical tools can be very helpful for supporting new pollutant detection. Nowadays, a combination of mass spectrometry and chromatographic techniques are the most basic environmental monitoring methods. In this paper, two organochlorine compound mass spectra classification systems were proposed. The classification models were developed within the framework of artificial neural networks (ANNs) and fast 1D and 2D molecular descriptor calculations. Based on the intensities of two characteristic MS peaks, namely, [M] and [M-35], two classification criterions were proposed. According to criterion I, class 1 comprises [M] signals with the intensity higher than 800 NIST units, while class 2 consists of signals with the intensity lower or equal than 800. According to criterion II, class 1 consists of [M-35] signals with the intensity higher than 100, while signals with the intensity lower or equal than 100 belong to class 2. As a result of ANNs learning stage, five models for both classification criterions were generated. The external model validation showed that all ANNs are characterized by high predicting power; however, criterion I-based ANNs are much more accurate and therefore are more suitable for analytical purposes. In order to obtain another confirmation, selected ANNs were tested against additional dataset comprising popular sunscreen agents disinfection by-products reported in previous works.

ChemFrag: Chemically meaningful annotation of fragment ion mass spectra

Identification and structural determination of small molecules by mass spectrometry is an important step in chemistry and biochemistry. However, the chemically realistic annotation of a fragment ion spectrum can be a difficult challenge. We developed ChemFrag, for the detection of fragmentation pathways and the annotation of fragment ions with chemically reasonable structures. ChemFrag combines a quantum chemical with a rule-based approach. For different doping substances as test instances, ChemFrag correctly annotates fragment ions. In most cases, the predicted fragments are chemically more realistic than those from purely combinatorial approaches, or approaches based on machine learning. The annotation generated by ChemFrag often coincides with spectra that have been manually annotated by experts. This is a major advance in peak annotation and allows a more precise automatic interpretation of mass spectra.

Identification of small molecules using accurate mass MS/MS search

Tandem mass spectral library search (MS/MS) is the fastest way to correctly annotate MS/MS spectra from screening small molecules in fields such as environmental analysis, drug screening, lipid analysis, and metabolomics. The confidence in MS/MS-based annotation of chemical structures is impacted by instrumental settings and requirements, data acquisition modes including data-dependent and data-independent methods, library scoring algorithms, as well as post-curation steps. We critically discuss parameters that influence search results, such as mass accuracy, precursor ion isolation width, intensity thresholds, centroiding algorithms, and acquisition speed. A range of publicly and commercially available MS/MS databases such as NIST, MassBank, MoNA, LipidBlast, Wiley MSforID, and METLIN are surveyed. In addition, software tools including NIST MS Search, MS-DIAL, Mass Frontier, SmileMS, Mass++, and XCMS2 to perform fast MS/MS search are discussed. MS/MS scoring algorithms and challenges during compound annotation are reviewed. Advanced methods such as the in silico generation of tandem mass spectra using quantum chemistry and machine learning methods are covered. Community efforts for curation and sharing of tandem mass spectra that will allow for faster distribution of scientific discoveries are discussed.

Software Tools and Approaches for Compound Identification of LC-MS/MS Data in Metabolomics

The annotation of small molecules remains a major challenge in untargeted mass spectrometry-based metabolomics. We here critically discuss structured elucidation approaches and software that are designed to help during the annotation of unknown compounds. Only by elucidating unknown metabolites first is it possible to biologically interpret complex systems, to map compounds to pathways and to create reliable predictive metabolic models for translational and clinical research. These strategies include the construction and quality of tandem mass spectral databases such as the coalition of MassBank repositories and investigations of MS/MS matching confidence. We present in silico fragmentation tools such as MS-FINDER, CFM-ID, MetFrag, ChemDistiller and CSI:FingerID that can annotate compounds from existing structure databases and that have been used in the CASMI (critical assessment of small molecule identification) contests. Furthermore, the use of retention time models from liquid chromatography and the utility of collision cross-section modelling from ion mobility experiments are covered. Workflows and published examples of successfully annotated unknown compounds are included.

Quantum Chemical Mass Spectrometry: Verification and Extension of the Mobile Proton Model for Histidine

The quantum chemical mass spectrometry for materials science (QCMS²) method is used to verify the proposed mechanism for proton transfer - the Mobile Proton Model (MPM) - by histidine for ten XHS tripeptides, based on quantum chemical calculations at the DFT/B3LYP/6-311+G* level of theory. The fragmentations of the different intermediate structures in the MPM mechanism are studied within the QCMS² framework, and the energetics of the proposed mechanism itself and those of the fragmentations of the intermediate structures are compared, leading to the computational confirmation of the MPM. In addition, the calculations suggest that the mechanism should be extended from considering only the formation of five-membered ring intermediates to include larger-ring intermediates. Graphical Abstract ᅟ.