A predictive science approach to aid understanding of electrospray ionisation tandem mass spectrometric fragmentation pathways of small molecules using density functional calculations

Fragmentation pathways of deprotonated ortho-hydroxybenzyl alcohol

The ortho, meta, and para isomers of hydroxybenzyl alcohol can be unequivocally distinguished by the collision-induced dissociation mass spectra of their anions. The presence of a prominent peak at m/z 121 for an elimination of a dihydrogen molecule renders the ortho-isomer spectrum markedly different from those of its meta and para congeners. Investigations carried out with deuterium-labeled isotopologues of the ortho isomer verified that the labile hydrogen atom on the hydroxyl group and one of the benzylic hydrogen atoms are specifically removed in the formation of the m/z 121 ion. The ortho-isomer spectrum also showed a prominent peak at m/z 93. Experimental data indicated that the m/z 93 product ion originates either from a two-step H2and CO elimination mechanism or from a direct loss of a HCHO molecule from the precursor anion. The intensity ratio of the m/z 93 and 94 peaks in the spectrum recorded from the m/z 124 ion generated from a sample of o-hydroxybenzyl alcohol dissolved in D2O supported the notion that the direct HCHO loss is the more dominant pathway for the generation of the phenolate ion under low activation conditions. In contrast, the two-step mechanism becomes the more dominant pathway under high collisional activation conditions. The spectrum also showed a weak peak at m/z 105 for a water loss. Based on computational data, the m/z 105 ion generated in this way appears to be a composite generated from a common ion-neutral complex intermediate in which a hydroxyl anion is positioned equidistantly between one of the benzylic hydrogens and a nearby hydrogen atom of the benzene ring. Upon activation, the complex dissociates to form either a phenide or a quinone methide anion. The reaction forming a carbon dioxide adduct under ion-mobility conditions was used to support the proposed water-loss mechanism.

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.

Relationship between Hammett's parameters and in silico density functional with tandem mass ESI-CID fragmentation: Dihydropyridines as prototypes

Over the years, with the instrumental analysis evolution, the relationships between the carried-out results with the data of theoretical analysis in silico and the Hammett's parameters have been reported. They have been very useful for chemical characterization of small organic molecules. Thus, this work aims at showing the feasibility and limitations for Hammett's and density functional theory applications in electrospray ionization-collision-induced dissociation (ESI-CID) fragmentation provision. For this, 13 dihydropyrimidinones para, meta, and orto monosubstituted were studied using ESI and CID in positive mode. As a result, it was observed that the main fragmentation includes the isocyanate and ethanol loses at low energy. Nevertheless, at higher energies, radical ions formed by McLafferty rearrangement were observed. The Hammett plots were correlated fragmentation profiles, showing good linearity for the [M + H]⁺ , which does not occur to radical ions and carbocation's. These tendencies had demonstrated that the stability of protonate and activation energy of secondary ions changes with the pKa. The density functional theory studies indicated that, both nitrogen atoms in the dihydropyrimidinone's prototypes are capable of being protonated. However, the activation energy of fragmentation products is not changed. Therefore, this work has shown information, which can be useful to understand tandem mass spectrometry in ESI-CID conditions for small organic molecules series. This is the first step for normalization of fragmentation pathway.

Competing fragmentation processes of β-substituted propanoate ions upon collision-induced dissociation

RATIONALE

When subjected to collisional activation, gas-phase carboxylate ions typically undergo decarboxylation. However, alternative fragmentation processes dominate when the carboxylate group is located within certain structural motifs. In this work, the fragmentation processes of β-substituted carboxylate ions are characterized to improve correlations between reactivity and structure.

METHODS

Mass spectra were collected using both ion trap and triple quadrupole mass spectrometers operating in the negative ion mode; collision-induced dissociation (CID) of ions was used to study the relationship between product ions and the structures of their precursor ions. Quantum mechanical computations were performed on a full range of reaction geometries at the MP2/6-311++G(2d,p)//B3LYP/6-31++G(2d,p) level of theory.

RESULTS

For a series of β-substituted carboxylate ions, a product ion corresponding to the anion of the β-substituent was obtained upon CID. Detailed computations indicated that decarboxylative elimination and at least one other fragmentation mechanism had feasible energetics for the formation of substituent anions differing in their gas-phase basicities. Predicted energetics for anti- and synperiplanar alignments in the transition structures for decarboxylative elimination correlated with the positions of crossover points in breakdown curves acquired for conformationally constrained ions.

CONCLUSIONS

The feasibility of more than one mechanism was established for the fragmentation of β-substituted propanoates. The contribution of each mechanistic pathway to the formation of the substituent anion was influenced by structural variations and conformational constraints, but mostly depended on the nature of the substituent. Copyright © 2016 John Wiley & Sons, Ltd.

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

The electron ionization mass spectra of four organic compounds are predicted based on the results of quantum chemical calculations at the DFT/B3LYP/6-311 + G* level of theory. This prediction is performed 'ab initio', i.e. without any prior knowledge of the thermodynamics or kinetics of the reactions under consideration. Using a set of rules determining which routes will be followed, the fragmentation of the molecules' bonds and the complete resulting fragmentation pathways are studied. The most likely fragmentation pathways are identified based on calculated reaction energies ΔE when bond cleavage is considered and on activation energies ΔE^‡ when rearrangements are taken into account; the final intensities of the peaks in the spectrum are estimated from these values. The main features observed in the experimental mass spectra are correctly predicted, as well as a number of minor peaks. In addition, the results of the calculations allow us to propose fragmentation pathways new to empirical mass spectrometry, which have been experimentally verified using tandem mass spectrometry measurements. Copyright © 2016 John Wiley & Sons, Ltd.

Non-targeted screening approaches for contaminants and adulterants in food using liquid chromatography hyphenated to high resolution mass spectrometry

The majority of analytical methods for food safety monitor the presence of a specific compound or defined set of compounds. Non-targeted screening methods are complementary to these approaches by detecting and identifying unexpected compounds present in food matrices that may be harmful to public health. However, the development and implementation of generalized non-targeted screening workflows are particularly challenging, especially for food matrices due to inherent sample complexity and diversity and a large analyte concentration range. One approach that can be implemented is liquid chromatography coupled to high-resolution mass spectrometry, which serves to reduce this complexity and is capable of generating molecular formulae for compounds of interest. Current capabilities, strategies, and challenges will be reviewed for sample preparation, mass spectrometry, chromatography, and data processing workflows. Considerations to increase the accuracy and speed of identifying unknown molecular species will also be addressed, including suggestions for achieving sufficient data quality for non-targeted screening applications.

McLafferty-type rearrangement of protonated N-[nicotinoyl]phenylethyl amines and consequent elimination of styrene

RATIONALE

McLafferty rearrangements occur in radical cations of molecules containing a carbonyl group and a γ hydrogen atom but are not common in the [M+H](+) ions of carbonyl compounds. We propose to investigate the collision-induced dissociation (CID) of the [M+H](+) ions of nicotinoyl and picolinoyl amides of 1- and 2-phenylethylamines to explore the possibility of McLafferty-type rearrangement.

METHODS

The compounds for study were synthesized by the reaction of methyl nicotinate or methyl picolinate with 1- and 2-phenylethylamines. The CID mass spectra of electrospray ionization (ESI)-generated protonated molecules were obtained using a QSTAR XL quadrupole time-of-flight (QTOF) mass spectrometer, and density functional theory (DFT) calculations using the B3LYP method were employed to elucidate the fragmentation mechanisms. The total electronic and thermal energies of intermediate transition states (TSs) and product ions are reported relative to those of the [M+H](+) ions.

RESULTS

CID of the [M+H](+) ions of N-[nicotinoyl]-2-phenylethylamine (1) yielded product ions of m/z 105 (1-phenylethyl cation) and 123 ([M+H-styrene](+) cation). The competitive formation of the ions of m/z 123 and 105 is proposed to involve a McLafferty-type rearrangement. Similarly, the [M+H](+) ions of the isomeric compound 2 and the N-[picolinoyl] phenylethyl amines (3 and 4) dissociate to yield ions of m/z 123 and 105.

CONCLUSIONS

A molecule of styrene was eliminated from the ESI-generated [M+H](+) ions of N-[nicotinoyl]phenylethylamines and the isomeric N-[picolinoyl]phenylethylamines, through a mechanism involving a McLafferty-type 1,5-H shift. The transition state energy for the 1,5-H shift is less for the amides of 1-phenylethylamine than for the amides of 2-phenylethylamine. The process occurs as a charge remote process and the presence of the pyridine ring is essential for the process.

Examination of small molecule losses in 5-methylpyranopelargonidin MS/MS CID spectra by DFT calculations

Pyranoanthocyanins are formed during food treatment and maturation (e.g. wine, juices), and they can be considered a natural alternative to artificial food colorants. Tandem mass spectrometry (MS/MS) is perhaps the most important technique in analysis of anthocyanin dyes. Knowledge of fragmentation pattern is a key aspect of their successful structural characterization. Polyphenolic compounds are known to lose small molecules during collision-induced dissociation (CID) in MS/MS experiments. However, the specific positions where such losses occur preferentially are unknown. The aim of this communication is to investigate the energetically most preferred places for H2 O and CO losses during the fragmentation of 5-methylpyranopelargonidin molecule by the means of computational chemistry (employing density functional theory) combined with CID MS/MS experiments and infrared multiphoton dissociation spectroscopy. Mechanisms responsible for the fragmentations were investigated, and optimal geometries and transition states were obtained. Cleavage of water as well as carbon monoxide occurs preferentially from the C-ring of flavonoid skeleton. In the most stable structure of 5-methylpyranopelargonidin, B-ring was found to be tilted with respect to the rest of the molecule. Planarization effort of the parent molecule contributes both to its decarbonylation and dehydration.

Predicting collision-induced dissociation spectra: semi-empirical calculations as a rapid and effective tool in software-aided mass spectral interpretation

RATIONALE

Fifteen molecules were modelled using quantum chemistry, prior to interpreting their collision-induced dissociation (CID) product ion spectra, in a 'blind trial' to establish if calculated protonation-induced bond elongation could be used to predict which bonds cleaved during CID. Bond elongation has the potential to be used as a descriptor predicting bond cleavage.

METHODS

The 15 molecules were modelled with respect to protonation-induced bond length changes using Density Functional Theory (DFT). Significant bond elongations were highlighted to flag potential bond cleavages. CID product ion spectra, obtained using positive ion electrospray ionisation (Waters Synapt G1), were interpreted to establish if observed bond cleavages correlated with calculated bond elongations. Calculations were also undertaken using AM1 (Austin Model 1) to see if this rapid approach gave similar results to the computationally demanding DFT.

RESULTS

The AM1-calculated bond elongations were found to be similar to those generated by DFT. All the polarised bonds observed to cleave (n = 82) had been calculated to elongate significantly. Protonation, possibly via proton migration, on the most electronegative atom in the bond appeared to initiate cleavage, leading to a 100% success rate in predicting the bonds that broke as a result of protonation on a heteroatom. Cleavage of carbon-carbon bonds was not predicted.

CONCLUSIONS

Cleavage of the polarised bonds appears to result from protonation on the more electronegative atom of the bond, inducing conformational changes leading to bond weakening. AM1-calculated bond length changes act as a descriptor for predicting bond cleavage. However, the impetus for cleavage of the unpolarised bonds may be product ion stability rather than bond weakening.

A systematic investigation of the fragmentation pattern of two furanoheliangolide C-8 stereoisomers using electrospray ionization mass spectrometry

RATIONALE

Budlein A is a sesquiterpene lactone (STL) with some reported biological activities. Pre-clinical studies to identify in vivo metabolites often employ hyphenated techniques such as liquid chromatography/tandem mass spectrometry (LC/MS/MS). It is also possible to use the fragmentation pattern obtained by Collision-Induced Dissociation (CID) and Higher Energy Collision-Induced Dissociation (HCD) to distinguish between the stereoisomers budlein A and centratherin.

METHODS

The experiments were carried out in the positive mode using four different spectrometers with an electrospray ionization (ESI) source: (a) Waters ACQUITY(®) TQD triple quadrupole mass spectrometer (QqQ), (b) AB Sciex API 4000 QTrap(®) (QqQ), (c) Bruker Daltonics micrOTOF™-Q II (time-of-flight, QTOF), and (d) Thermo Scientific LTQ Orbitrap XL hybrid FTMS (Fourier transform mass spectrometer). Computational chemistry studies helped to identify the protonation sites. The B3LYP/6-31G(d) model furnished the equilibrium geometries and energies.

RESULTS

The stereochemistry (α- or β-orientation) of the centratherin and budlein A side-chain esters influences the fragmentation pattern recorded on QqQ, QTOF, and Orbitrap-HCD. On QqQ, centratherin releases the side chain, to generate the m/z 275 fragment ion, whereas budlein A gives the m/z 83 fragment ion. On QTOF and Orbitrap-HCD, only budlein A affords the m/z 293 and 83 fragment ions, respectively.

CONCLUSIONS

The data suggest that proton migration governs the fragmentation pathways under α- or β-orientation. The difference in the QqQ, QTOF, and Orbitrap-HCD spectral profiles of each isomer can help to distinguish between centratherin and budlein A using MS/MS, which becomes an alternative to nuclear magnetic resonance (NMR) analysis.

Rapid and simultaneous in situ assessment of aflatoxins and stilbenes using silica plate imprinting mass spectrometry imaging

A fast and direct combination of techniques for simultaneous mycotoxin and phytoalexin identification in peanut skin and kernel is described. Silica Plate Imprinting Laser Desorption/Ionization Mass Spectrometry Imaging (SPILDI-MSI) is a powerful technique that exhibits great advantages, such as solvent-free and matrix-free characteristics, as well as no sample preparation or separation steps. It also permits accurate identification of mycotoxins and phytoalexins with unique fingerprint profiles in just a few seconds. Results are expressed as chemical images of the 4 identified types of aflatoxins (B1, B2, G1 and G2) and a stilbenoid (resveratrol). Also, SPILDI-MSI allows the comparison between the spatial distribution of aflatoxins and resveratrol found in kernel and skin. This novel application has proven to be useful for instantaneous qualitative assessment of aflatoxins and stilbenoids both in the peanut skin and kernel and offers precise tracking of fungal contamination in nuts and other foodstuffs.