From QCEIMS to QCxMS: A Tool to Routinely Calculate CID Mass Spectra Using Molecular Dynamics

QC-GN2oMS2: a Graph Neural Net for High Resolution Mass Spectra Prediction

Predicting the mass spectrum of a molecular ion is often accomplished via three generalized approaches: rules-based methods for bond breaking, deep learning, or quantum chemical (QC) modeling. Rules-based approaches are often limited by the conditions for different chemical subspaces and perform poorly under chemical regimes with few defined rules. QC modeling is theoretically robust but requires significant amounts of computational time to produce a spectrum for a given target. Among deep learning techniques, graph neural networks (GNNs) have performed better than previous work with fingerprint-based neural networks in mass spectra prediction. To explore this technique further, we investigate the effects of including quantum chemically derived information as edge features in the GNN to increase predictive accuracy. The models we investigated include categorical bond order, bond force constants derived from extended tight-binding (xTB) quantum chemistry, and acyclic bond dissociation energies. We evaluated these models against a control GNN with no edge features in the input graphs. Bond dissociation enthalpies yielded the best improvement with a cosine similarity score of 0.462 relative to the baseline model (0.437). In this work we also apply dynamic graph attention which improves performance on benchmark problems and supports the inclusion of edge features. Between implementations, we investigate the nature of the molecular embedding for spectra prediction and discuss the recognition of fragment topographies in distinct chemistries for further development in tandem mass spectrometry prediction.

Emerging contaminants: A One Health perspective

Environmental pollution is escalating due to rapid global development that often prioritizes human needs over planetary health. Despite global efforts to mitigate legacy pollutants, the continuous introduction of new substances remains a major threat to both people and the planet. In response, global initiatives are focusing on risk assessment and regulation of emerging contaminants, as demonstrated by the ongoing efforts to establish the UN's Intergovernmental Science-Policy Panel on Chemicals, Waste, and Pollution Prevention. This review identifies the sources and impacts of emerging contaminants on planetary health, emphasizing the importance of adopting a One Health approach. Strategies for monitoring and addressing these pollutants are discussed, underscoring the need for robust and socially equitable environmental policies at both regional and international levels. Urgent actions are needed to transition toward sustainable pollution management practices to safeguard our planet for future generations.

Studying the Intrinsic Reactivity of Chromanes by Gas-Phase Infrared Spectroscopy

Tandem mass spectrometry is routinely used for the structural analysis of organic molecules, but many fragmentation reactions are not well understood. Because several potential structures can correspond to a measured mass, the assignment of product ions is ambiguous using mass spectrometry alone. Here, we combine mass spectrometry with high-resolution gas-phase infrared spectroscopy and computational chemistry tools to identify product ion structures and derive collision-induced fragmentation mechanisms of the chromane derivatives Trolox and Methyltrolox. We find that protonated Trolox and Methyltrolox fragment identically via dehydration and decarbonylation, while deprotonated ions display substantially diverging reactivities. For deprotonated Methyltrolox, we observe unusual radical fragmentation reactions and suggest a [1,2]-Wittig rearrangement involving aryl migration in the gas phase. Overall, the combined experimental and theoretical approach presented here revealed complex proton dynamics and intramolecular rearrangement reactions, which expand our understanding on structure-reactivity relationships of isolated molecules in different protonation states.

Identifying isomers in Chinese traditional medicine via density functional theory and ion fragmentation simulation software QCxMS

In the realm of electrospray ionization mass spectrometry (ESI-MS), distinguishing among isomers poses a significant challenge due to the minimal spectral differences that often arise from their subtle structural differences. This makes the accurate identification of these compounds through solely experimental spectra a daunting task. Computational chemistry has emerged as a pivotal tool in bridging the gap between experimental observations and theoretical understanding. This study used the MS fragmentation simulation software, QCxMS, to model the spectra of five groups of isomers, encompassing 11 compounds, found in the traditional Chinese medicine, Zhishi Xiebai Guizhi Decoction. By comparing the spectra predicted through computational methods with those derived from Ultra-high performance liquid chromatography-quadrupole-time of flight-mass spectrometry (UPLC-Q-TOF-MS) experiments, it was observed that, following the optimization of simulation parameters, QCxMS was capable of generating reliable spectra for all examined compounds. Notably, the data calculated under both GFN1-xTB and GFN2-xTB levels exhibited no significant discrepancies. Further analysis enabled the identification of the principal fragments of the 11 compounds from the theoretical data, facilitating the deduction of their fragmentation pathways. The Density Functional Theory (DFT) method was subsequently applied to compute the primary fragmentation energies of these compounds. The findings revealed a congruence between the energy data calculated using both thermodynamic and kinetic approaches and the observed fragment abundance of the isomers. This alignment providing a more precise theoretical framework for understanding the mechanisms underlying the generation of fragment ion differences among isomers.

In-Electrospray source Hydrogen/Deuterium exchange coupled to multistage fragmentation for the investigation of the protonation and fragmentation pathways of gas phase ions

Identification of molecules in complex natural matrices relies on matching the fragmentation spectra of ions under investigation and the spectra acquired for the corresponding analytical standards. Currently, there are many databases of experimentally measured tandem mass spectrometry spectra (such as NIST, MzCloud, and Metlin), and considerable progress has been made in the development of software for predicting tandem mass spectrometry fragments in silico using combinatorial, machine learning, and quantum chemistry approaches (such as MetFrag, CFM-ID, and QCxMS). However, the electrospray ionization molecules can be ionized at different sites (protonated or deprotonated), and the fragmentation spectra of such ions are different. Here, we are using the combination of the in-ESI source hydrogen/deuterium exchange reaction and MSn fragmentation for the investigation of the fragmentation pathways for different protomers of organic molecules. It is shown that the distribution of the deuterium in the fragment ions reflects the presence of different protomers. For several molecules, the distribution of deuterium was traced up to the MS5 level of fragmentation revealing many unusual and unexpected effects. For example, we investigated the loss of HF from the ciprofloxacin and norfloxacin ions and observed that for ions protonated at -COOH group, the eliminating hydrogen always comes from -NH group. When ions are protonated at another site, the elimination of hydrogen with a probability of 30% occurs from the -NH group, and with a probability of 70%, it originates from other sites on the molecule. Such effects were not described previously. Quantum chemical simulation was used for the verification of the protonated structures and simulation of the corresponding fragmentation spectra.

IDSL.GOA: gene ontology analysis for interpreting metabolomic datasets

Biological interpretation of metabolomic datasets often ends at a pathway analysis step to find the over-represented metabolic pathways in the list of statistically significant metabolites. However, definitions of biochemical pathways and metabolite coverage vary among different curated databases, leading to missed interpretations. For the lists of genes, transcripts and proteins, Gene Ontology (GO) terms over-presentation analysis has become a standardized approach for biological interpretation. But, GO analysis has not been achieved for metabolomic datasets. We present a new knowledgebase (KB) and the online tool, Gene Ontology Analysis by the Integrated Data Science Laboratory for Metabolomics and Exposomics (IDSL.GOA) to conduct GO over-representation analysis for a metabolite list. The IDSL.GOA KB covers 2393 metabolic GO terms and associated 3144 genes, 1,492 EC annotations, and 2621 metabolites. IDSL.GOA analysis of a case study of older versus young female brain cortex metabolome highlighted 82 GO terms being significantly overrepresented (FDR < 0.05). We showed how IDSL.GOA identified key and relevant GO metabolic processes that were not yet covered in other pathway databases. Overall, we suggest that interpretation of metabolite lists should not be limited to only pathway maps and can also leverage GO terms as well. IDSL.GOA provides a useful tool for this purpose, allowing for a more comprehensive and accurate analysis of metabolite pathway data. IDSL.GOA tool can be accessed at https://goa.idsl.me/ .

Robustness of Threshold Collision-Induced Dissociation Simulations for Bond Dissociation Energies

The threshold collision-induced dissociation (T-CID) method is the workhorse for gas-phase bond dissociation energy (BDE) measurements. However, T-CID does not measure BDEs directly; instead, BDEs are obtained by fitting simulated data to the experimental data. We previously observed several large discrepancies between the computed and experimental BDEs. To analyze the reliability of the experimental values, we previously reported a study of the dissociation rate models in the simulation. Here, we report a study of the collision simulation part, specifically in the L-CID (ligand CID) program. We show that the BDE values are robust even to intentionally introduced mistakes in the simulations, varying in most cases by less than 3 kcal mol-1. The most significant exception is the collisional energy transfer (CET) simulation, which led to deviations larger than 10 kcal mol-1. However, we found that the BDEs obtained with explicitly simulated CET distributions deviated by only 3 kcal mol-1 from those simulated with the original model. Collectively, our results suggest that the T-CID-derived BDE values are robust and are likely to be accurate.

IDSL.GOA: Gene Ontology Analysis for Interpreting Metabolomic datasets

Biological interpretation of metabolomic datasets often ends at a pathway analysis step to find the over-represented metabolic pathways in the list of statistically significant metabolites. However, definitions of biochemical pathways and metabolite coverage vary among different curated databases, leading to missed interpretations. For the lists of genes, transcripts and proteins, Gene Ontology (GO) terms over-presentation analysis has become a standardized approach for biological interpretation. But, GO analysis has not been achieved for metabolomic datasets. We present a new knowledgebase (KB) and the online tool, Gene Ontology Analysis by the Integrated Data Science Laboratory for Metabolomics and Exposomics (IDSL.GOA) to conduct GO over-representation analysis for a metabolite list. The IDSL.GOA KB covers 2,393 metabolic GO terms and associated 3,144 genes, 1,492 EC annotations, and 2,621 metabolites. IDSL.GOA analysis of a case study of older vs young female brain cortex metabolome highlighted 82 GO terms being significantly overrepresented (FDR <0.05). We showed how IDSL.GOA identified key and relevant GO metabolic processes that were not yet covered in other pathway databases. Overall, we suggest that interpretation of metabolite lists should not be limited to only pathway maps and can also leverage GO terms as well. IDSL.GOA provides a useful tool for this purpose, allowing for a more comprehensive and accurate analysis of metabolite pathway data. IDSL.GOA tool can be accessed at https://goa.idsl.me/.

Evolutionary metabolomics of specialized metabolism diversification in the genus Nicotiana highlights N-acylnornicotine innovations

Specialized metabolite (SM) diversification is a core process to plants' adaptation to diverse ecological niches. Here, we implemented a computational mass spectrometry-based metabolomics approach to exploring SM diversification in tissues of 20 species covering Nicotiana phylogenetics sections. To markedly increase metabolite annotation, we created a large in silico fragmentation database, comprising >1 million structures, and scripts for connecting class prediction to consensus substructures. Together, the approach provides an unprecedented cartography of SM diversity and section-specific innovations in this genus. As a case study and in combination with nuclear magnetic resonance and mass spectrometry imaging, we explored the distribution of N-acylnornicotines, alkaloids predicted to be specific to Repandae allopolyploids, and revealed their prevalence in the genus, albeit at much lower magnitude, as well as a greater structural diversity than previously thought. Together, the data integration approaches provided here should act as a resource for future research in plant SM evolution.

Highly accurate and large-scale collision cross sections prediction with graph neural networks

The collision cross section (CCS) values derived from ion mobility spectrometry can be used to improve the accuracy of compound identification. Here, we have developed the Structure included graph merging with adduct method for CCS prediction (SigmaCCS) based on graph neural networks using 3D conformers as inputs. A model was trained, evaluated, and tested with >5,000 experimental CCS values. It achieved a coefficient of determination of 0.9945 and a median relative error of 1.1751% on the test set. The model-agnostic interpretation method and the visualization of the learned representations were used to investigate the chemical rationality of SigmaCCS. An in-silico database with 282 million CCS values was generated for three different adduct types of 94 million compounds. Its source code is publicly available at https://github.com/zmzhang/SigmaCCS . Altogether, SigmaCCS is an accurate, rational, and off-the-shelf method to directly predict CCS values from molecular structures.

Comparative Evaluation of Electron Ionization Mass Spectral Prediction Methods

During the past decade promising methods for computational prediction of electron ionization mass spectra have been developed. The most prominent ones are based on quantum chemistry (QCEIMS) and machine learning (CFM-EI, NEIMS). Here we provide a threefold comparison of these methods with respect to spectral prediction and compound identification. We found that there is no unambiguous way to determine the best of these three methods. Among other factors, we find that the choice of spectral distance functions play an important role regarding the performance for compound identification.

Collision Cross Section Prediction Based on Machine Learning

Ion mobility-mass spectrometry (IM-MS) is a powerful separation technique providing an additional dimension of separation to support the enhanced separation and characterization of complex components from the tissue metabolome and medicinal herbs. The integration of machine learning (ML) with IM-MS can overcome the barrier to the lack of reference standards, promoting the creation of a large number of proprietary collision cross section (CCS) databases, which help to achieve the rapid, comprehensive, and accurate characterization of the contained chemical components. In this review, advances in CCS prediction using ML in the past 2 decades are summarized. The advantages of ion mobility-mass spectrometers and the commercially available ion mobility technologies with different principles (e.g., time dispersive, confinement and selective release, and space dispersive) are introduced and compared. The general procedures involved in CCS prediction based on ML (acquisition and optimization of the independent and dependent variables, model construction and evaluation, etc.) are highlighted. In addition, quantum chemistry, molecular dynamics, and CCS theoretical calculations are also described. Finally, the applications of CCS prediction in metabolomics, natural products, foods, and the other research fields are reflected.

Ab initio molecular dynamics calculations on electron ionization induced fragmentations of C4F7N and C5F10O for understanding their decompositions under discharge conditions

C₄F₇N and C₅F₁₀O are the most promising SF₆ alternatives as eco-friendly insulating gaseous mediums in electrical engineering. It is necessary to clarify their electrical stability and decomposition mechanisms. In this work, we first introduced our experimental results for decomposition products of C₄F₇N/CO₂ and C₅F₁₀O/synthetic air mixtures under partial discharge and spark discharge conditions. Then, we performed ab initio molecular dynamics (AIMD) simulations on the typical decomposition products. The simulations were performed under standard electron impact mass spectrometry (EI-MS); thus, the statistical results of the mass spectra were compared with those of the experimentally obtained standard mass spectra from the NIST database. The AIMD simulation method in simulating the electron-induced ionization process was verified and found to be reliable. Finally, the calculations were also performed for C₄F₇N and C₅F₁₀O with incident electron energies of 20 eV and 70 eV, respectively. The dominant pathway for both gases is the formation of CF₃⁺ with the fracture of the C-C bond. The AIMD simulation is able to predict the decomposition channels after electron-impact ionization without any preconceived knowledge of fragmentation pathways, which provides a novel insight into understanding the decomposition mechanisms of C₄F₇N and C₅F₁₀O under different discharge conditions with different energies.

Simple In Vitro 18O Labeling for Improved Mass Spectrometry-Based Drug Metabolites Identification: Deep Drug Metabolism Study

The identification of drug metabolites formed with different in vitro systems by HPLC-MS is a standard step in preclinical research. In vitro systems allow modeling of real metabolic pathways of a drug candidate. Despite the emergence of various software and databases, identification of compounds is still a complex task. Measurement of the accurate mass, correlation of chromatographic retention times and fragmentation spectra are often insufficient for identification of compounds especially in the absence of reference materials. Metabolites can "slip under the nose", since it is often not possible to reliably confirm that a signal belongs to a metabolite and not to other compounds in complex systems. Isotope labeling has proved to be a tool that aids in small molecule identification. The introduction of heavy isotopes is done with isotope exchange reactions or with complicated synthetic schemes. Here, we present an approach based on the biocatalytic insertion of oxygen-18 isotope under the action of liver microsomes enzymes in the presence of ¹⁸O₂. Using the local anesthetic bupivacaine as an example, more than 20 previously unknown metabolites were reliably discovered and annotated in the absence of the reference materials. In combination with high-resolution mass spectrometry and modern methods of mass spectrometric metabolism data processing, we demonstrated the ability of the proposed approach to increase the degree of confidence in interpretating metabolism data.

Understanding of protomers/deprotomers by combining mass spectrometry and computation

Multifunctional compounds may form different prototropic isomers under different conditions, which are known as protomers/deprotomers. In biological systems, these protomer/deprotomer isomers affect the interaction modes and conformational landscape between compounds and enzymes and thus present different biological activities. Study on protomers/deprotomers is essentially the study on the acidity/basicity of each intramolecular functional group and its effect on molecular structure. In recent years, the combination of mass spectrometry (MS) and computational chemistry has been proven to be a powerful and effective means to study prototropic isomers. MS-based technologies are developed to discriminate and characterize protomers/deprotomers to provide structural information and monitor transformations, showing great superiority than other experimental methods. Computational chemistry is used to predict the thermodynamic stability of protomers/deprotomers, provide the simulated MS/MS spectra, infrared spectra, and calculate collision cross-section values. By comparing the theoretical data with the corresponding experimental results, the researchers can not only determine the protomer/deprotomer structure, but also investigate the structure-activity relationship in a given system. This review covers various MS methods and theoretical calculations and their devotion to isomer discrimination, structure identification, conformational transformation, and phase transition investigation of protomers/deprotomers.

Calculation of Mass Spectra with the QCxMS Method for Negatively and Multiply Charged Molecules

Analysis and validation of a mass spectrometry (MS) experiment are usually performed by comparison to reference spectra. However, if references are missing, measured spectra cannot be properly matched. To close this gap, the Quantum Chemical Mass Spectrometry (QCxMS) program has been developed. It enables fully automatic calculations of electron ionization (EI) and positive ion collision-induced dissociation (CID) mass spectra of singly charged molecular ions. In this work, the extension to negative and multiple ion charge for the CID run mode is presented. QCxMS is now capable of calculating structures carrying any charge, without the need for pretabulated fragmentation pathways or machine learning of database spectra. Mass spectra of four single negatively charged and two multiple positively charged organic ions with molecular sizes from 12 to 92 atoms were computed and compared to reference spectra. The underlying Born-Oppenheimer molecular dynamics (MD) calculations were conducted using the semiempirical quantum mechanical GFN2-xTB method, while for some small molecules, ab initio DFT-based MD simulations were performed. Detailed insights into the fragmentation pathways were gained, and the effects of the computed charge assignments on the resulting spectrum are discussed. Especially for the negative ion mode, the influence of the deprotonation site to create the anion was found to be substantial. Doubly charged fragments could successfully be calculated fully automatically for the first time, while higher charged structures introduced severe assignment problems. Overall, this extension of the QCxMS program further enhances its applicability and underlines its value as a sophisticated toolkit for CID-based tandem MS structure elucidation.

Computer-aided drug design, quantum-mechanical methods for biological problems

Quantum chemistry enables to study systems with chemical accuracy (<1 kcal/mol from experiment) but is restricted to a handful of atoms due to its computational expense. This has led to ongoing interest to optimize and simplify these methods while retaining accuracy. Implementing quantum mechanical (QM) methods on modern hardware such as multiple-GPUs is one example of how the field is optimizing performance. Multiscale approaches like the so-called QM/molecular mechanical method are gaining popularity in drug discovery because they focus the application of QM methods on the region of choice (e.g., the binding site), while using efficient MM models to represent less relevant areas. The creation of simplified QM methods is another example, including the use of machine learning to create ultra-fast and accurate QM models. Herein, we summarize recent advancements in the development of optimized QM methods that enhance our ability to use these methods in computer aided drug discovery.

Quantum Chemistry-based Molecular Dynamics Simulations as a Tool for the Assignment of ESI-MS/MS Spectra of Drug Molecules

In organic mass spectrometry, fragment ions provide important information on the analyte as a central part of its structure elucidation. With increasing molecular size and possible protonation sites, the potential energy surface (PES) of the analyte can become very complex, which results in a large number of possible fragmentation patterns. Quantum chemical (QC) calculations can help here, enabling the fast calculation of the PES and thus enhancing the mass spectrometry-based structure elucidation processes. In this work, the previously unknown fragmentation pathways of the two drug molecules Nateglinide (45 atoms) and Zopiclone (51 atoms) were investigated using a combination of generic formalisms and calculations conducted with the Quantum Chemical Mass Spectrometry (QCxMS) program. The computations of the de novo fragment spectra were conducted with the semi-empirical GFNn-xTB (n=1, 2) methods and compared against Orbitrap measured electrospray ionization (ESI) spectra in positive ion mode. It was found that the unbiased QC calculations are particularly suitable to predict non-evident fragment ion structures, sometimes contrasting the accepted generic formulation of fragment ion structures from electron migration rules, where the "true" ion fragment structures are approximated. For the first time, all fragment and intermediate structures of these large-sized molecules could be elucidated completely and routinely using this merger of methods, finding new undocumented mechanisms, that are not considered in common rules published so far. Given the importance of ESI for medicinal chemistry, pharmacokinetics, and metabolomics, this approach can significantly enhance the mass spectrometry-based structure elucidation processes and contribute to the understanding of previously unknown fragmentation pathways.

In Silico Collision Cross Section Calculations to Aid Metabolite Annotation

The interpretation of ion mobility coupled to mass spectrometry (IM-MS) data to predict unknown structures is challenging and depends on accurate theoretical estimates of the molecular ion collision cross section (CCS) against a buffer gas in a low or atmospheric pressure drift chamber. The sensitivity and reliability of computational prediction of CCS values depend on accurately modeling the molecular state over accessible conformations. In this work, we developed an efficient CCS computational workflow using a machine learning model in conjunction with standard DFT methods and CCS calculations. Furthermore, we have performed Traveling Wave IM-MS (TWIMS) experiments to validate the extant experimental values and assess uncertainties in experimentally measured CCS values. The developed workflow yielded accurate structural predictions and provides unique insights into the likely preferred conformation analyzed using IM-MS experiments. The complete workflow makes the computation of CCS values tractable for a large number of conformationally flexible metabolites with complex molecular structures.

PyFragMS-A Web Tool for the Investigation of the Collision-Induced Fragmentation Pathways

Dissociation induced by the accumulation of internal energy via collisions of ions with neutral molecules is one of the most important fragmentation techniques in mass spectrometry (MS), and the identification of small singly charged molecules is based mainly on the consideration of the fragmentation spectrum. Many research studies have been dedicated to the creation of databases of experimentally measured tandem mass spectrometry (MS/MS) spectra (such as MzCloud, Metlin, etc.) and developing software for predicting MS/MS fragments in silico from the molecular structure (such as MetFrag, CFM-ID, CSI:FingerID, etc.). However, the fragmentation mechanisms and pathways are still not fully understood. One of the limiting obstacles is that protomers (positive ions protonated at different sites) produce different fragmentation spectra, and these spectra overlap in the case of the presence of different protomers. Here, we are proposing to use a combination of two powerful approaches: computing fragmentation trees that carry information of all consecutive fragmentations and consideration of the MS/MS data of isotopically labeled compounds. We have created PyFragMS-a web tool consisting of a database of annotated MS/MS spectra of isotopically labeled molecules (after H/D and/or ¹⁶O/¹⁸O exchange) and a collection of instruments for computing fragmentation trees for an arbitrary molecule. Using PyFragMS, we investigated how the site of protonation influences the fragmentation pathway for small molecules. Also, PyFragMS offers capabilities for performing database search when MS/MS data of the isotopically labeled compounds are taken into account.

Evaluating the Accuracy of the QCEIMS Approach for Computational Prediction of Electron Ionization Mass Spectra of Purines and Pyrimidines

Mass spectrometry is the most commonly used method for compound annotation in metabolomics. However, most mass spectra in untargeted assays cannot be annotated with specific compound structures because reference mass spectral libraries are far smaller than the complement of known molecules. Theoretically predicted mass spectra might be used as a substitute for experimental spectra especially for compounds that are not commercially available. For example, the Quantum Chemistry Electron Ionization Mass Spectra (QCEIMS) method can predict 70 eV electron ionization mass spectra from any given input molecular structure. In this work, we investigated the accuracy of QCEIMS predictions of electron ionization (EI) mass spectra for 80 purine and pyrimidine derivatives in comparison to experimental data in the NIST 17 database. Similarity scores between every pair of predicted and experimental spectra revealed that 45% of the compounds were found as the correct top hit when QCEIMS predicted spectra were matched against the NIST17 library of >267,000 EI spectra, and 74% of the compounds were found within the top 10 hits. We then investigated the impact of matching, missing, and additional fragment ions in predicted EI mass spectra versus ion abundances in MS similarity scores. We further include detailed studies of fragmentation pathways such as retro Diels–Alder reactions to predict neutral losses of (iso)cyanic acid, hydrogen cyanide, or cyanamide in the mass spectra of purines and pyrimidines. We describe how trends in prediction accuracy correlate with the chemistry of the input compounds to better understand how mechanisms of QCEIMS predictions could be improved in future developments. We conclude that QCEIMS is useful for generating large-scale predicted mass spectral libraries for identification of compounds that are absent from experimental libraries and that are not commercially available.

Fragmentation efficiency of phenethylamines in electrospray ionization source estimated by theoretical chemistry calculation

Small molecules with polar functional groups, including substituted phenethylamines, are commonly analyzed by liquid chromatography-mass spectrometry (LC-MS) with electrospray ionization (ESI). Analyte molecules are mostly detected in protonated and cation-adducted forms through positive-ion electrospray ionization-mass spectrometry (ESI-MS). However, the ESI of substituted phenethylamines commonly provides an intense signal of fragment ions by ESI in-source collision-induced dissociation (IS-CID), which hinders the unambiguous identification of phenethylamines. This phenomenon was approximated as a unimolecular dissociation model, and the dissociation efficiency was evaluated by various quantum chemistry calculations to determine the ESI IS-CID efficiency. The calculated results were consistent with the experimental data, when the dissociation threshold energy of phenethylamines was calculated using the post-Hartree-Fock (post-HF) method, CCSD(t)/cc-pVTZ//MP2(full)/6-311++G(d,p). In contrast to post-HF methods, the utilization of density functional theory calculations with a suitable functional is recognized as an accurate and competitive low-cost approach. In particular, ωB97-XD, M06-2X-D3, and recently developed Minnesota functionals, such as M11, MN12-SX, and MN15, provided reliable results, as in the case of the post-HF method. The results obtained by the recently developed double hybrid functionals, DSD-PEBP86-D3(BJ), PBE0-DH, and PBE-QIDH, were also reliable. The consideration of ESI IS-CID can facilitate the identification of analyte molecules because most phenethylamines, except for N-methylated analogs, provide an intense signal in the ESI mass spectrum.

Nontargeted Screening Using Gas Chromatography-Atmospheric Pressure Ionization Mass Spectrometry: Recent Trends and Emerging Potential

Gas chromatography–high-resolution mass spectrometry (GC–HRMS) is a powerful nontargeted screening technique that promises to accelerate the identification of environmental pollutants. Currently, most GC–HRMS instruments are equipped with electron ionization (EI), but atmospheric pressure ionization (API) ion sources have attracted renewed interest because: (i) collisional cooling at atmospheric pressure minimizes fragmentation, resulting in an increased yield of molecular ions for elemental composition determination and improved detection limits; (ii) a wide range of sophisticated tandem (ion mobility) mass spectrometers can be easily adapted for operation with GC–API; and (iii) the conditions of an atmospheric pressure ion source can promote structure diagnostic ion–molecule reactions that are otherwise difficult to perform using conventional GC–MS instrumentation. This literature review addresses the merits of GC–API for nontargeted screening while summarizing recent applications using various GC–API techniques. One perceived drawback of GC–API is the paucity of spectral libraries that can be used to guide structure elucidation. Herein, novel data acquisition, deconvolution and spectral prediction tools will be reviewed. With continued development, it is anticipated that API may eventually supplant EI as the de facto GC–MS ion source used to identify unknowns.