Underlying Mechanisms of Chromatographic H/D, H/F, cis/trans and Isomerism Effects in GC-MS


1. Introduction

The boiling points (Bps) of chemical compounds are constant under fixed environmental conditions, notably temperature and atmospheric pressure. The Bps of many protiated and deuterated solvents such as acetone and benzene are known. The Bp of protiated acetone (CH3COCH3) is 56.0 °C and that of deuterated acetone (CD3COCD3) is 55.2 °C. The Bp of protiated benzene (C6H6) is 80.0 °C and that of deuterated benzene (C6D6) is 79.1 °C. When we define the boiling point isotope effect, bpIE, as the ratio of the Bp of protiated to deuterated solvents (Formula (1)), the bpIE values would be 1.0145 for acetone and 1.0114 for benzene due to the lower Bp of the deuterated solvents (Figure 1).

The Bps of thousands of natural and synthetic analytes are unknown. Their Bp values are not only difficult to determine, but their determination would also be extremely expensive. Consequently, the boiling points, and most likely many other physicochemical parameters of analytes, seem not to be useful physicochemical measures to quantify isotope effects of analytes. On the other hand, data available for many other solvents, notably for purely inorganic water (Bp, 99.98 °C for H2O, Bp, 101.4 °C for D2O), suggest that the utility of the Bp would be of limited value.

In gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) including their tandem variants, analytes and their stable isotope labeled analogs, i.e., the isotopologs, which are used as internal standards (ISs), are separated mass spectrometrically due to their different mass-to-charge ratio (m/z) values. Isotopologs may also differ in their retention times, with the deuterated (D) analytes having, as a rule, shorter retention times, i.e., tR(D), than their protiated (H), i.e., tR(H), analogs (Figure 2) [1]. In order to discriminate this effect from other kinds of isotope effects, such as the kinetic isotope effect, this phenomenon is named in the present work the chromatographic isotope effect, IEC, i.e., hdIEC for the chromatographic isotope effect due to H/D. Its extent can be determined by dividing the experimentally observed tR(H) values by the tR(D) values (see Formula (2)), analogous to the bpIE values mentioned above. The difference in the tR(H) and tR(D) values, i.e., δ(H/D) (Figure 2), may also be useful for the determination of the extent of the chromatographic isotope effect (see Formula (3)).

δ(H/D) = tR(H)tR(D)

Based on the theory of the partition chromatography [2], the smaller retention times of 2H-isotopologs compared to their 1H-isotopologs would result from a weaker interaction strength of the 2H-isotopologs with the stationary phase. The closer the hdIEC value to the unity (1.0000) is, and the smaller the δ(H/D) is, the weaker is the hdIEC.

The underlying mechanisms of chromatographic isotope effects have been little explored. Many factors are likely to contribute to hdIEC, notably including differences in the physical and chemical properties of the isotopologs (e.g., C-H/C-D bond energies, boiling points), as described above for acetone and benzene. The lower Bp values of CD3COCD3 and C6D6 compared to those of CH3COCH3 and C6H6 could be, in analogy, due to somewhat weaker inter-molecular interactions of deuterated analytes compared to their protiated congeners.

In GC-MS and LC-MS, isobaric analytes, i.e., analytes with identical m/z values but different stereochemistry, can be possibly separated by chromatography but not mass spectrometry. This is, for instance, the case of cis– and trans-fatty acids [3]. Empirically, trans-fatty acids have shorter retention times than their congeners cis-fatty acids in GC, thus resembling the hdIEC effect. In analogy, the chromatographic cis/trans-effect is named in the present work the chromatographic cis/trans-effect, ctEC, and can be calculated by Formula (4). In this Formula, tR(cis) and tR(trans) are the retention times of the cis– respectively trans-analytes.

ctEC = tR(cis)/tR(trans)

Again, based on the theory of the partition chromatography, the smaller retention times of trans-analytes compared to their cis-analytes would result from a weaker interaction strength of the trans-analytes with the GC stationary compared with the cis-analytes.

The author of the present article has hypothesized that the underlying mechanisms in hdIEC and ctEC are principally the same. These effects were examined by considering my own experimentally observed data for different classes of analytes and data published by other groups for unlabeled cis– and trans-fatty acids and for their 2H- and 13C-isotopologs [3,4,5,6]. As retention times of analytes were not available from many published articles (see Table 1) in order to determine hdIEC and ctEC values, only a small number of papers have been considered in the present work. All data used in the present study have been generated by validated GC-MS methods that included proper chemical derivatization of the analytes. The focus was on free amino acids and free fatty acids. Thus, amino acids (AAs) were analyzed as methyl ester-pentafluoropropionyl derivatives (AA-Me-PFP), and free fatty acids (FAs) either as methyl ester (FA-Me) or pentafluorobenzyl ester (FA-PFB) derivatives. The results of the present examination support the compiled hypothesis.
GC-MS and LC-MS are widely used in “classical” quantitative analyses and in metabolomics studies (Table 1). Nowadays, quantitative analyses are increasingly named targeted analyses. The expressions untargeted, nontargeted, non-targeted and semi-targeted analyses are widely used in metabolomics studies (Table 1), without a precise definition and often arbitrarily. The present work focused on quantitative GC-MS analyses, in which isotopologs had been used as internal standards. Chromatographic isotope effects and other types of chromatographic effects also occur in LC-MS and can even reach higher values. The majority of the studies considered in the present work are based on “classical” quantitative analyses using isotopologs and quadrupole-based, i.e., low-resolution mass spectrometers. Generally, chromatographic effects are independent of mass spectrometric resolution.

2. Methods

2.1. GC-MS Analyses Performed in the Author’s Group

Analyses were performed by selected-ion monitoring (SIM) or selected-reaction monitoring (SRM) on a quadrupole (Q) type GC-MS and GC-MS/MS apparatus, namely TSQ 45, TSQ 7000 and ISQ (ThermoFisher, Dreieich, Germany). Fused silica capillary GC columns were used as specified below and in the Results section. Plasma, serum, and urine samples of human subjects had been collected and analyzed in previous peer-reviewed studies [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. Helium and methane served as carrier and reagent gases for negative-ion chemical ionization (NICI), respectively. Aliquots (1 µL) of toluene, ethyl acetate or other extracts of analyte derivatives had been injected in the splitless mode. The concentrations of the stable-isotope labelled analogs (isotopologs) of the analytes added to the plasma and urine samples (at the start of the sample workup) were all relevant for the respective biological samples.
Amino acids were analyzed as methyl ester pentafluoropropionyl (Me-PFP) derivatives using Optima 17 (15 m length, 0.25 mm i.d., and 0.25-µm film thickness) [26]. The oven used was kept at 40 °C for 0.5 min and ramped to 210 °C at a rate of 15 °C/min and then to 320 °C at a rate 35 °C/min. Commercially available isotopologs and in situ prepared trideutero-methyl esters of amino acids were used as internal standards.
Prostaglandin E2 (PGE2) was analyzed as a pentafluorobenzyl (PFB) ester methoxime (MO) trimethylsilyl (TMS) derivative (PFB-MO-TMS). Fatty acids and hydroxylated fatty acids were analyzed as PFB and PFB-TMS derivatives, respectively. Dimethylamine (DMA) was analyzed as a pentafluorobenzoyl derivative. Malondialdehyde (MDA) was analyzed as a di-PFB derivative. Information about the GC-MS methods of these and other analytes can be found in the References [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25].

2.2. GC-MS Analyses Performed by Other Groups

The experimental conditions of the work performed by other authors are reported in the Results section. For better understanding, the nomenclature used in the original articles was retained. The original data and the newly calculated hdIEC and ctEC values are presented in the Results section of the present work in form of Tables and Figures as appropriate.

2.3. Calculations and Data Presentation

hdIEC and c/tEC values were calculated using Formulas (2) and (4), respectively. GraphPad Prism Version 7 for Windows (GraphPad Software, San Diego, CA, USA) was used for the statistical analyses and preparation of graphs. Chemical structures of the investigated analytes and their derivatives were drawn using ChemDraw 15.0 Professional (PerkinElmer Informatics, Germany).

3. Results

3.1. GC-MS Analysis of Miscellaneous Analytes—Author’s Group

The IEC values for the PFB derivatives of 14N-nitrite, 15N-nitrite, 14N-nitrate, and 15N-nitrate in plasma and urine were close to 1.0000, with the difference in the retention times being not higher than 0.18 s [1].
Table 2 summarizes the GC-NICI-MS retention times of unlabeled (d0Me) amino acids (AAs) as methyl ester pentafluoropropionyl (AA-Me-PFP) derivatives, the chromatographic isotope effect (hdIEC) values of the AA-Me-PFP derivatives of unlabeled and deuterium labeled (d3Me) amino acid derivatives as measured simultaneously in urine samples of healthy humans. The retention times of all derivatives were measured with high precision (coefficient of variation). Under the same GC-MS conditions, very similar results were obtained from analyses of plasma samples of 38 healthy humans. The hdIEC values of the AA-Me-PFP derivatives ranged between 1.002 and 1.006 (Table 2).
Figure 3 shows that the difference (δ(H/D)) in the retention times of unlabeled (AA-d0Me-PFP) and deuterated amino acids (AA-d3Me-PFP) does not correlate with the chromatographic isotopic H/D effect (hdIEC). The greatest δ(H/D) values were observed for glutamate (Glu), aspartate (Asp), and methionine (Met).
Figure 4 shows the relationship between the molecular weight (MW) of native and derivatized amino acids in human urine and the retention time of the AA-Me-PFP derivatives. The retention time increased with increasing molecular weight of the non-derivatized amino acids (r2 = 0.83) and of their Me-PFP derivatives (r2 = 0.55). The weaker linearity of the Me-(PFP)1 and Me-(PFP)2 derivatives compared to the Me-(PFP)3 derivatives is presumably due to the combined effects of the structures of the side chain and of their derivatized amine and hydroxyl groups of the amino acids. Figure 4D suggests that the extent of the chromatographic isotope effect decreases with the decreasing molecular weight of the AA-Me-PFP derivative.
Similar relationships were observed for the same amino acids analyzed in human plasma samples (Figure 5). The number of F atoms in the AA-Me-PFP derivatives decreases the chromatographic isotope effect: the bigger the number of the F atoms is, the weaker is hdIEC.
Figure 6 shows that δ(H/D) is largely independent of the van der Waals volume, the relative molecular mass and the hydrophobicity of the unlabeled amino acids analyzed by GC-NICI-MS as Me-PFP derivatives. The greatest δ(H/D) values were observed for Glu, Asp, and Met. Glu and Asp are dicarboxylic amino acids and form dimethyl ester. Met is a monocarboxylic sulfur-containing (thioether) amino acid.
Table 3 compiles the calculated hdIEC values observed from GC-NICI-MS and GC-NICI-MS/MS analyses of perfluorinated derivatives of biogenic amines, polyamines, and other physiological substances including creatinine.
The hdIEC values ranged between 0.9993 for the malondialdehyde (MDA) derivative, MDA-(PFB)2, and 1.0125 for the 5-hydroxyeicosanoic acid (HEA) derivative, HEA-PFP-TMS (Table 3). There was a positive correlation after Spearman (r = 0.7075, p = 0.0002) between the hdIEC values and the number of D atoms in the derivatives of the analytes, but not between the hdIEC and the MW values of the derivatives (r = 0.354, p = 0.106). Yet, hdIEC and MW of the native, non-derivatized analytes correlated with each other (r = 0.499, p = 0.018). These results suggest that the PFB, PFBz, and PFP residues, which were introduced by the derivatization of the analytes, decrease the extent of the hdIEC.

3.2. hdIEC in the GC-MS/MS Analysis of Fatty Acid Methyl Esters—Tintrop et al. 2022, 2023 [5,6]

Fatty acids (FAs) occur in biological samples in free form (FFA) and as methyl esters (FAMEs). The simultaneous GC-EI-MS/MS determination of FAs and FAMEs was reported [5,6]. This method includes solvent-free solid-phase microextraction arrow headspace extraction and in situ isotope-labeling of FFAs with deuterated methanol (CD3OD). It utilizes the chromatographic isotope effect and the 3 Da-shift due to d3Me compared to d0Me. Sulfuric acid-catalyzed derivatization of FFAs with CD3OD was performed for 20 min with 4 vol% CD3OD and a pH value of 2.1.

The simultaneous analysis of FAMEs was performed on a GC 2010 with a MS/MS TQ8040 (Shimadzu Deutschland GmbH, Duisburg, Germany) in the MRM (SRM) mode. A Zebron ZB-FAME capillary column (30 m × 0.25 mm × 0.20 µm, Phenomenex, Torrance, USA) was used. The oven temperature program started at 40 °C, which was kept for 5 min, and then raised with a rate of 5 °C/min to 210 °C, where it was held for 5 min. Helium was used as a carrier gas (1.8 mL/min) and argon as a collision gas. Injection of the analytes was performed by splitless thermal desorption.

The investigated fatty acids included saturated (C6:0 to C22:0) and unsaturated C16 and C18 fatty acids, as well as perdeuterated C16:0, i.e., C16:0-d31, and C17:0, i.e., C17:0-d33. The retention times of the FA-d0Me and FA-d3Me derivatives reported by the authors [6] were used to calculate the hdIEC values. The newly calculated data are summarized in Table 4 and illustrated in Figure 7. The retention time of the FA-Me derivatives did not increase linearly with the number of the C-atoms of the native fatty acids when considering all fatty acids (r2 = 0.9813).

The H/D isotope effect is evident. Expectedly, the strongest chromatographic isotope effects are seen for the methyl esters of the perdeuterated fatty acids C16:0-d31 and C17:0-d33. The IEC values decrease with the increasing length of the fatty acid molecule. These observations suggest that the entire H/D effect of the d3Me group decreases with the increasing length of the fatty acid molecules by a factor of 10 (from 0.6% to 0.06% with respect to the retention time). According to the manufacturer’s information, Zebron™ ZB-FAME is based on a high-cyanopropyl (G48) chemistry (O-Si-(C3H6CN)n) and has high polarity. Analogous to the isotope effect seen in the GC-MS analysis of the methyl esters of amino derivatives, one may assume that the deuterium atoms of the d3Me ester group of the fatty acids weaken the interaction of the lipophilic fatty acid molecules with the immobilized liquid phase of the GC column. This effect is maximum in C16:0-d31 (IEC, 1.0171) and C17:0-d33 (hdIEC, 1.0126). The ctEC values for C18:1c and C18:1t are calculated to be each 1.0041 for d0Me and d3Me, i.e., they are higher than the hdIEC values of 1.0011.

Similar effects have been reported by Alexander et al. [28] for 45 unlabeled cis/trans fatty acid methyl esters using a highly polar cyanopropylsiloxane SP 2560 fused-silica capillary column (100 m × 0.25 mm × 0.20 µm). The GC column oven temperature was linearly programmed from 80 °C to 220 °C at 8 °C/min. Initial and final temperatures were held for 2 min and 32 min, respectively, with a total run time of 55 min, using helium as the carrier gas (0.67 mL/min). For example, the retention times were reported to be 35.7 min for C22:0, 37.2 min for C22:1, 44.6 min for C22:4, 46.5 min for C22:5, and 52.0 min for C22:6.
Ecker et al. [29] used a highly polar BPX70 GC column (10 m × 0.10 mm × 0.20 µm; 70% cyanopropyl polysilphenyl-siloxane) for the GC-MS analysis of cis- and trans-fatty acid methyl esters. The initial oven temperature was held for 0.75 min, then programmed to increase with 40 °/min to 155 °C, with 6 °C/min to 210 °C, to finally reach 250 °C with 15 °C/min, which was held for 2 min. FAMEs eluted between 2.54 min (C8:0) and 14.16 min (C28:0). Unsaturated FAMEs had higher retention times than saturated FAMEs. The ctEC values were 1.0154 for C18:1 (c9) and C18:1 (t9), and 1.0364 for C18:2 (c9,c12) and C18:2 (t9,t12) [29]. Within unsaturated FAMEs, the retention time was higher when the –C=C-double bonds were near to the terminal methyl group of the fatty acids.

3.3. hdIEC in the GC-MS Analysis of Fatty Acid-Pentafluorobenzyl Esters—Quehenberger et al. 2011 [4]

Unlabeled and deuterium-labeled FAs (Table 5) were analyzed by GC-NICI-MS as PFB ester derivatives. GC-NICI-MS analysis was carried out on an Agilent 6890N gas chromatograph equipped with an Agilent 7683 autosampler (Santa Clara, CA, USA) using an Agilent 5973 mass selective detector. The FA-PFB esters dissolved in isooctane were injected with a pulsed (25 psi) splitless injection mode onto a Zebron ZB-1 column (15 m × 0.25 mm i.d., coated with 100% dimethylpolysiloxane; Phenomenex, Torrance, CA, USA). Helium (0.9 mL/min) was used as a carrier gas. The GC oven temperature was programmed from 150 °C to 270 °C at 10 °C/min, ramped to 240 °C at 40 °C/min and held at 240 °C for 1 min. The injector and the transfer line were kept at 250 °C and 280 °C, respectively. Methane was used as the reagent gas for NICI at an ion-source temperature of 150 °C. SIM of the [M-PFB] anions of the FA-PFB esters was performed. The [M-PFB] ions correspond to the carboxylate anions.
The retention time increased linearly with the increasing number of C-atoms of the native fatty acids (Figure 8). The highest hdIEC value was observed for the shortest fatty acid (i.e., C12:0), and the lowest hdIEC value for the longest fatty acids (i.e., C24:0 and C26:0). There was no dependency of the hdIEC value upon the number of the deuterium atoms in the fatty acid molecules.

These observations suggest that the H/D isotope effect is evident in the GC-MS analysis of the FA-PFB ester derivatives. However, the isotope effect is weak. A possible explanation could be that the interaction of the PFB ester residue of the fatty acid derivatives with the stationary phase outweighs the interaction of the H/D atoms of the fatty acid skeleton.

3.4. Isotope Effects in the GC-MS Analysis of Cis- and Trans-Fatty Acid Pentafluorobenzyl Esters—Kuiper et al. 2018 [3]

Kuiper and colleagues developed a GC-NICI-MS method for the quantitative analysis of trans-fatty acids in human plasma, serum, and red blood cells (RBCs). As the trans-fatty acids coexist in biological samples at considerably lower concentrations, a “special” GC column known to allow simultaneous analysis of trans– and cis-fatty acids was used, i.e., an Agilent Select FAME.

Samples (100 µL) were combined with 100 μL of an internal standard (IS) solution that contained isotopologs. Subsequently, the samples were hydrolyzed (2 mL of 10% v/v 6 M HCl in acetonitrile followed by 2 mL of 10% v/v 10 M NaOH in methanol), each carried out at 104 °C for 45 min. After neutralization with 6 M HCl, the free FAs were extracted with hexane (three times, 2 mL each). The solvent was removed under vacuum (Genevac, Stone Ridge, NY, USA) and the samples were derivatized at room temperature for 15 min with 100 μL of 7% PFB-Br in acetonitrile and 10 μL triethylamine as the base catalyst. The fatty acid PFB (FA-PFB) esters were extracted with hexane (500 μL) and transferred to autosampler vials for GC-MS analysis. Samples were handled in glass vials to minimize contamination of samples with FAs from plastic supplies.

GC-MS analysis of 27 FAs and 18 ISs was carried out on a 7890/5975C GC–MSD from Agilent Technologies (Santa Clara, CA, USA). The inlet temperature was 240 °C and 1-μL aliquots were injected with a 100:1 split ratio using a Gerstel Multipurpose Sampler MPS (Gerstel, Mülheim an der Ruhr, Germany) equipped with a cool drawer set at 10 °C. The carrier gas was hydrogen at a flow rate of 2 mL/min. An Agilent Select FAME (200 m × 250 μm × 0.25 μm) GC column on a cyanopropyl basis (CP-7421) was used. The column temperature was ramped from 50 °C to 160 °C at 40 °C/min, held at 160 °C for 10 min, increased by 1 °C/min to 175 °C, then by 0.5 °C/min to 210 °C, and finally by 35 °C/min to 260 °C, where it was held for 25 min. NICI was performed using methane as the reagent gas. The ions [M-PFB] were used in the SIM mode (Table 6). The transfer line, source, and quadrupole temperatures were 260 °C, 230 °C, and 150 °C, respectively. The electron multiplier voltage was adjusted throughout the chromatographic run to increase sensitivity for lowabundance FAs and prevent detector saturation for high-abundance FAs.
The fatty acids analyzed and the m/z values of the ions [M-PFB] used in GC-NICI-MS analyses are reported in Table 6. This Table also summarizes the retention times of the FA-PFB derivatives obtained from the GC-NICI-MS analyses of the indicated trans-fatty acids (t), the regular fatty acids including the cis-fatty acids (c), and their 2H- and 13C-isotopologs. The retention times reported by Kuiper et al. [3] were used in the present work to newly calculate the cis-FA-PFB-to trans-FA-PFB ratios (c/t) and the hdIEC values, as applicable.
There were three cis/trans-FA pairs. The FA-PFB derivatives of the cis-fatty acids had longer retention times than the corresponding trans-fatty acids. The c/t values amounted to 1.0233, 1.0171, and 1.0163 (Table 6). The retention times of the FA-PFB derivatives of the trans-fatty acids and their corresponding 13C-isotopes were almost identical, indicating no chromatographic isotope effect due to 13C/12C. The FA-PFB derivatives of di-unsaturated fatty acids had higher retention times than their respective mono-unsaturated fatty acids. Note the considerably longer retention times compared to those in the previous studies (Table 4, Table 5 and Table 6).

There are nine pairs of FAs and their 2H-labeled isotopologs. The FA-PFB derivatives of the unlabeled fatty acids had longer retention times than their 2H-isotopologs. The calculated hdIEC values ranged between 1.0428 for myristic acid (C14:0 and D27-C14:0) and 1.0009 for arachidonic acid (C20:4n-6,9,12,15 and D8-C20:4n-6,9,12,15).

These observations suggest that the FA-PFB ester derivatives of trans-fatty acids behave towards their cis-fatty acids in the same manner as behave non-deuterated towards deuterated fatty acids with respect to gas chromatography in fused-silica capillary columns. Note the long retention time range of the investigated FA-PFB ester derivatives in the study (i.e., 56 min to 116 min).

A further interesting observation is that the retention times of the FA-PFB ester derivatives increased with the increasing grade of unsaturation of fatty acids, e.g., within C18– and C20-fatty acids (Figure 9, Table 6 and Table 7), suggesting enhancement of interaction between the immobilized stationary phase of the GC column and the backbone of the fatty acids. In that study, the retention times of the FA-PFB ester derivatives did not increase linearly with the number of the C-atoms of the native fatty acids (Figure 9 vs. Figure 8).
The median concentrations of trans-fatty acids in plasma, serum, and RBCs of 66 healthy donors ranged between 1.8 µM and 2.1 µM for C16:1n-7 t, 6.2 µM and 8.1 µM for C18:1n-9 t, 8.8 µM and 11.4 µM for C18:1n-7 t, and 0.6 µM and 1.0 µM for C18:2n-6 t,9 t [3]. These concentrations account for 0.01–0.10% of total fatty acids in blood.

4. Discussion

Chromatographic separation of derivatized analytes in GC-MS is based on continuously occurring interactions of the analyte derivatives between a mostly silicone-based stationary phase (of small inner diameter and film-thickness) immobilized inside the GC column and mostly helium as the mobile phase. The just now “free“ analyte derivatives are “carried” by helium through the GC column until they are released into the ion-source of the GC-MS apparatus. The physicochemical properties of the analytes (e.g., boiling point, molecular weight, chemical structure, and possibly changing shape and orientation in the gas phase) contribute to the chromatographic effects, which are finally manifested in the retention times of the analytes. The retention time of analytes is a useful measure to quantify chromatographic effects. In quantitative GC-MS analyses, stable-isotope labeled analytes are used as internal standards.

Kinetic isotope effects, chromatographic isotope effects, and chromatographic cis/trans-effects have been known for several decades, but the underlying mechanisms are incompletely understood. Chromatographic isotope effects (IECs) are assumed to originate from differences in physicochemical properties of unlabeled and stable-isotope labeled analytes that result from the differences introduced into the analytes mainly by the heavy isotopes of H (hdIEC), C, N and O, i.e., 2H (D), 13C and 15N, respectively. Perhaps easier to understand are the chromatographic cis/trans-effects (ctEC), as the orientation of the molecules in the space may differ considerably, for instance in dependence on the residues in olefinic analytes R1-C=C-R2.

The retention time (tR) in gas chromatography (GC) and liquid chromatography (LC) is an experimentally ascertainable integral parameter, which incorporates all factors that are involved in the chromatographic process. In GC, they include boiling points vapor pressure, chemical composition, and three-dimensional structure of the analyte derivatives, the chemical composition of the stationary phase, adsorption/desorption processes, temperature of the injector port, initial and subsequently increasing temperature of the GC column and its dimensions (length, inner diameter, film thickness), as well as the nature and flow rate of the carrier gas. The present work addressed the issue of hdIEC and ctEC effects, and investigated these phenomena on a quantitative basis by re-examining data reported in the literature mainly by four groups including the author’s group. We hypothesized that the chromatographic isotope- and cis/trans-effects can be quantified by using the retention times of analyte derivatives in GC and introduced the parameters hdIEC and ctEC. The hdIEC and ctEC values were calculated by Formulas (2) and (4), respectively. The main working hypothesis was that both effects can be explained by differences in the interaction of the analytes with the stationary phase of fused-silica capillary GC columns. One class of analytes were amino acids that were analyzed as methyl ester pentafluoropropionic (PFP) derivatives. The second class free fatty acids that were analyzed as methyl ester (Me) or pentafluorobenzyl ester (PFB) derivatives.

The utilized GC columns differed in the chemistry of the stationary phase, and in part in length (15 m, 30 m, 200 m), had a comparable internal diameter, and small differences in film thickness(0.20 µm and 0.25 µm) of the immobilized stationary phase. The volume of the GC lumen was calculated to be about 0.7 mL, 1.5 mL and 9.6 mL, respectively. Under consideration of the reported carrier gas, and presumably laminar than turbulent flow, it is calculated that the carrier gas would need about 0.7 min to 5 min for a passage through the GC columns. The shape of the analytes inside the GC columns is unknown. In the case of hydrophobic analyte derivatives such as the FA-PFB derivatives, oblongness could possibly prevail. The interaction of the analyte derivatives and the stationary phase is considered to be due to different forces including van der Waals forces. Other effects, such as the organic solvent extract used in GC-MS analyses, “condensation” of analyte derivatives at lower GC column temperatures, and competition between the derivatives of analytes and their isotopologs, are possible but were not considered in this work.

Schmarr et al. [30] and Thakur et al. [31] reported on normal and inverse chromatographic isotope effects for various analytes using different stationary phases. In the case of normal isotopic effects, deuterated analytes elute later than the protiated analogs. In the case of inverse isotopic effects, deuterated analytes elute earlier than the protiated analogs. Thus, the observations described in the present work can be categorized as inverse isotopic effects. Schmarr et al. [30] and Thakur et al. [31] did not report retention times of the isotopologs pairs so that the hdIEC values could not be reported in this article. It was observed that nonpolar stationary phases often exhibited an inverse isotope effect, whereas polar stationary phases often showed a normal isotope effect. For instance, by using the stationary phase SPB-20 (bonded; poly(20% diphenyl/80% dimethyl siloxane) phase; polarity number, 111), C6D4-(CD3)2 (o-xylene-d10) eluted in front of C6H4-(CH3)2 (o-xylene-d0) with a peak resolution of 2.70 (complete baseline separation, i.e., inverse hdIEC). By using the stationary phase IL111i (non-bonded; 1,5-di(2,3-dimethylimidazolium)pentane bis(trifluoromethanesulfonyl)imide phase; polarity number, 12), C6D4-(CD3)2 (o-xylene-d10) eluted behind C6H4-(CH3)2 (o-xylene-d0) with a peak resolution of 1.00 (i.e., normal hdIEC). All analytes investigated with the stationary phase IL111i were found to be associated with normal hdIEC [31]. The stationary phase IL111i was found to be the only stationary phase that interacts with molecules by π–π, dipole–dipole, and dipole–induced dipole interactions in addition to typical hydrogen bonding, dispersive, and acid-base type interactions [31].

4.1. Chromatographic Isotope Effects

The IEC values observed for PFB-O*NO2 and PFB-*NO2 (an asterisk * indicates 14N/15N; difference of 7% in atom mass) were very close to 1.0000 and indicate lack of a measurable IEC effect and virtually the same interaction extent of the central 14N and 15N atoms of PFB-O*NO2 and PFB-*NO2 due to sterical hindrance through the neighboring C and O atoms.

The IEC values observed for the 13C isotopologs of the fatty acid derivatives are almost identical to those of the 12C isotopologs, while the simultaneously analyzed 2H isotopologs clearly caused IEC. The C atoms of analytes do not interact directly with the stationary phase of the GC column. Furthermore, the 13C atoms of the 13C-labeled analytes differ by only 8% in atom mass from the 12C atoms. These two factors are likely responsible for the missing isotope effect in the 13C-labeled fatty acid derivatives.

The two N-methyl groups of dimethylamine (DMA) and metformin (Metf) in their PFBz, respectively, PFP derivatives caused strong hdIEC effects due to the six H/D atoms of the methyl groups which are not hindered sterically in their interaction with the stationary phase of the GC column. The comparably stronger hdIEC effects are likely to result from the greater differences in the mass of the H and D atoms (100% increase in atom mass).

In agreement with the observations presented in the Results section, Benchekroun et al. [32] reported that analytes with CD3 groups influence the hdIEC of derivatives of caffeine and of its metabolites obtained on a 14% cyanopropylphenyl)methylpolysiloxane fused-silica capillary column in GC-EI-MS (EI, electron ionization). The hdIEC (range, 1.0012 to 1.0045) was also found to increase with the number of the D atoms and their position in the molecules [32].
The CODATA-recommended value for the charge radius of D is 2.5 times bigger than the charge radius of H: 2.13 fm vs. 0.84 fm (https://physics.nist.gov/cuu/Constants/index.html; accessed on 10 December 2024). The Van der Waals radius is 120 pm for H, 135 pm for F, 152 pm for O, and 160 pm for C. The rate of a reaction involving a C–H bond is typically 6–10 times faster than the corresponding C–D bond, meaning that the heavier atom favors a stronger bond [33,34]. These two physicochemical factors are likely to be dominant contributors to the chromatographic isotope effects observed in the present study, with the van der Waals forces being more relevant for fatty acids [35], but not for amino acids.
The bigger charge radius of the D atom lends to the deuterated analytes larger molecular volumes compared to the protiated analytes. At the same time, the stronger C-D bonds tend to interact more weakly with the stationary phase of GC columns such as those made of 50% methylpolysiloxane/50% phenylpolysiloxane or 14% (cyanopropylphenyl)methylpolysiloxane. The synergetic effects of these factors is supported by the highest hdIEC values observed for the perdeuterated fatty acids. Other non-deuterated voluminous structures of the analytes such as in the PFB esters are likely to attenuate the hdIEC [4].
It is reasonable to assume that deuterated analytes may have lower boiling points and higher volatilities analogous to organic solvents such as acetone and benzene, because of stronger intra-molecular compared to inter-molecular interactions. Similar differences in the boiling points also prevail in perfluorinated chemicals including toluene (Bp of CH3-C6H5, 110.5 °C; Bp of CF3-C6F5, 104 °C), which would result in H/F-related effects (bpEH/F) of 1.0625. The bpEH/F values for PFB-Br and pentafluoropropionic anhydride (PFPA) used for the derivatization and GC-MS analysis of fatty acids, respectively, amino acids are calculated by Formula (5) to be 1.1375 and 2.3857. Thus, it can be reasonably expected that the physicochemical properties of the PFB and PFP residues are transferred to the analyte derivatives.
whereas Bp(H) and Bp(F) are the boiling points of the protiated and perfluorinated solvent, respectively.
The hdIEC values observed in the present work may have been influenced by different experimental parameters, mainly including the length and chemistry of the GC columns used, the starting GC oven temperature and the GC oven temperature gradient. Thus, the starting oven temperature was 150 °C in the GC-MS analysis of the FA-PFB esters by Quehenberger et al. [4]), but only 40 °C in the GC-MS analysis of the FA-Me esters by Tintrop et al. [6]. At (very) high starting GC oven temperatures and high GC oven temperature gradients, chromatographic separation may be inadequate for studying chromatographic isotope effects.

4.2. Chromatographic cis/trans-Effects in Fatty Acid Derivatives

Many of the issues discussed above on the chromatographic isotope effects can be translated to the chromatographic cis/trans-effects. Possibly the greatest difference between IEC and c/tEC refers to the “interaction surface” of the analytes due to the cis/trans isomerism as observed in the case of unsaturated fatty acids. The cis/trans configuration “subdivides” the analyte derivative, depending on the position of –C=C-double bond(s), into two or more cis– and trans-oriented structures in the lumen of the GC column. The trans-isomers probably do not to interact with the GC stationary phase to the same degree as the cis-isomers, because the cis-configuration allows for a stronger interaction of the –C=C-double bonds with the stationary phase. This is schematically illustrated for the PFB ester derivatives of the saturated stearic acid (C18:0) and the unsaturated cis-fatty acids oleic acid (C18:1), α-linolenic acid and γ-linolenic acid (both C18:3) (Figure 10).

The shape of the FA-PFB derivatives could be elongated for saturated fatty acids but curvy for unsaturated cis-fatty acids. Curvy fatty acids may possibly collide more often with the GC stationary phase through their –C=C-double bonds, thus increasing the residence time in the lumen of the GC column.

4.3. Chromatographic Effects Due to Other Types of Isomerism

Similar to the chromatographic cis/trans-isomerism in unsaturated fatty acids is the stereo-isomerism of isobaric prostaglandins (PG). The most commonly analyzed isomeric prostaglandins are the F2-isoprostanes (four types, 64 isomers) [36]. The best investigated F2-isoprostanes are prostaglandin F (PGF, 9α,11α,15S-trihydroxy-()-prosta-5Z,13E-dien-1-oic-acid) and 8-iso-prostaglandin F (8-iso-PGF,9α,11α,15S-trihydroxy-()-prosta-5Z,13E-dien-1-oic-acid), and prostaglandin E2 (PGE2,9-oxo-11α,15S-dihydroxy-()-prosta-5Z,13E-dien-1-oic acid) and 8-iso-prostaglandin E2 (8-iso-PGE2,9-oxo-11α,15S-dihydroxy-()-prosta-5Z,13E-dien-1-oic acid) [23,37]. The two epimers have an opposite configuration at only one stereogenic center out of at least two. PGF and 8-iso-PGF differ in the space orientation of the α-chain on the cyclopentane C-8 atom in the molecules. The PFB-TMS derivatives of PGF and 8-iso-PGF behave like a cis/trans isomer with respect to the plane of the cyclopentane ring. The retention time of PGF-PFB-(TMS)3 in GC-MS/MS was determined to be 23.70 min and that of 8-iso-PGF-PFB-(TMS)3 22.97 min [37], resulting in a retention time ratio of 1.0318. This value is even greater than the c/tEC values of cis/trans FA-PFB derivatives. No such effect was observed for the PFB-MO-(TMS)2 derivatives of PGE2 and 8-iso-PGE2, but it occurred in the bigger pentafluorobenzyloxime (PFBO) derivatives, i.e., PFB-PFBO-(TMS)2 derivatives [23]. hdIEC effects of similar extent were observed in both prostaglandins [23,37]. This effect could be named the chromatographic epimer effect epiEC, and the underlying mechanism is assumed to be the same as the ctEc.
Isomerism is not limited to endogenous substances, but it also applies to synthetic compounds such as the polychlorinated terphenyls (PCTs), which include ortho-, meta– and para-homologs. In GC-MS (DB-5 column, (5%-phenyl)-methylpolysiloxane), the elution order of isomers was found to be ortho-PCTs, meta-PCTs and para-PCTs [38], obviously in the order of increasing melting points (58 °C, 86 °C, 212 °C, respectively) and boiling points (337 °C, 379 °C, 389 °C, respectively). Similar effects have been observed by GC-MS (Elite–624 capillary column; 6% cyanopropyl phenyl 94% methyl, 30 m × 0.25 mm id, 1.4 µm film-thickness) for the xylene isomers ortho-xylene, meta-xylene and para-xylene, with the meta-xylene and para-xylene being more closely and earlier eluting than ortho-xylene (tR: 12.2 min, 12.00 min, 12.98 min, respectively; Bp: 139 °C, 137 °C, 144 °C, respectively) [39]. In the case of the o/m/p-isomers, the position of the two methyl groups and the size of the aromatic ring available for interaction with the stationary phase of the GC column are likely to influence more strongly the retention time of the xylene isomers. In the case of the o/m/p-isomers of the polychlorinated terphenyls, the surface of the three aromatic rings of the longish para-PCTs isomers is likely to be more strongly exposed to the stationary phase of the GC column than that of ortho– and meta-PCTs isomers.
In GC-MS, GC separation depends on the relative strength of the analytes–stationary phase interaction. It was calculated that the orientation polarization of meta-xylene is 80 times greater than of para-xylene, and that the permanent dipole moment of meta-xylene is 15 times greater than of para-xylene [22]. Kanai and colleagues assumed that compared with para-xylene, meta-xylene has a greater Debye–Keesom (both belong to the van der Waals forces) interaction tendency with the polar polyethylene glycol (PEG) capillary column [40].

4.4. Sterical Effects on Chromatographic Isotope Effects

The hdIEC of the PFB derivatives of malondialdehyde (MDA) (OHC-CH2-CHO) and [1,3-2H2]malondialdehyde (O2HC-CH2-C2HO), i.e., OHC-C(PFP)2-CHO and O2HC-C(PFB)2-C2HO is very close to 1.0000 despite the involvement of two D atoms. This is likely to be due to the sterical effect of the two voluminous PFB residues, which avoid almost completely the interaction of the D atoms with the stationary phase of the GC column (Figure 11). This observation suggests that the melting and boiling points may have a lower impact on hdIEC.

4.5. Possible Implications of Chromatographic Isotope Effects in Metabolomics

GC-MS and GC-MS/MS are mature technologies and widely used for the quantitative determination of known endogenous substances in biological samples. Nowadays, quantitative GC-MS-based analyses are increasingly called targeted-based metabolomics when several analytes are simultaneously analyzed, as described in the present work for amino acids and fatty acids. Targeted GC-MS analyses require the use of isotopologs as an internal standard for “absolute” quantification and commonly one (e.g., fatty acids), two (e.g., amino acids), and even three different derivatization reactions (e.g., prostaglandins). GC-MS is increasingly used in untargeted metabolomics studies, similar to untargeted LC-MS/MS, aiming to discover novel compounds. Analytical protocols combining targeted and untargeted GC-MS metabolomics have been proposed [41].
In particular, untargeted metabolomics faces formidable analytical challenges, mainly due to the chemical diversity of physiological substances and their expanding concentration range over several orders of magnitude in biological samples (e.g., pM to mM, i.e., 10−12 to 10−3 M). Quantitative GC-MS or targeted GC-MS metabolomics referring to a class of analytes, such as fatty acids or amino acids, are based on the use of analytical methods that yield relatively homogenous extracts. Untargeted GC-MS metabolomics are constrained to make analytical compromises at the cost of analytical quality in order to include as many analytes as possible. One potential drawback is the infamous matrix effect, which is especially strong in LC-MS, notably in the shotgun mode, i.e., without chromatography. In quantitative GC-MS (i.e., targeted GC-MS metabolomics), the use of isotopologs and specific analytical procedures allows for highly reliable quantitative analyses. Moreover, the chromatographic isotope effect, which is easily available, may be suitable for implementation in quality control (QC) systems as shown for amino acids in human plasma, serum, and urine samples from clinical studies [1]. Consideration of the chromatographic isotope effect will increase the quality assurance as it combines two main closely interrelated quality criteria, which are the retention times of the analytes themselves and their isotopologs.
For example, study human plasma samples (n = 297–353) and QC samples (n = 54–64) were analyzed by targeted GC-MS for amino acids within eight runs [42]. hdIEC and δ(H/D) values were close in the study and the QC plasma samples and ranged between 1.002 and 1.006, and 0.84 s and 2.64 s, respectively.

5. Conclusions

The greatest IEC effects in isotopologs are exerted by the deuterium atoms of the lipophilic analytes, because they interact directly with the mostly hydrophobic stationary phase of fused-silica capillary GC columns. This interaction is weaker than that of the H atoms of the analytes, thus resulting in a weaker retention and shorter retention times of 2H-isotopologs. Even a small number of D atoms causes clearly and highly reproducibly quantifiable hdEC effects. Even thoroughly 13C-labeled analytes do not exert measurable IEC effects, because their “inner” C atoms do not interact directly with the stationary phase of the GC column and their possibly indirect action is rather negligible. This is likely to be true for 15N- and 18O-isotopologs even for derivatives such as PFB-*NO2, PFB-O*NO2, and FA-C*O2-PFB.

The hdIEC effects are clearly observable in derivatives of cis– and trans-fatty acids and are accompanied by stronger but equally directed ctEC effects caused by the cis/trans-isomerism thus resembling the hdIEC effects. The comparably shorter retention times of the FA-PFB derivatives are the result of an interaction of a structural part of the elongate trans-isomers, which is of a shorter duration than that of the curvy cis-isomers. Such an explanation is applicable to the epimer stereo-isomerism of prostaglandins. A further explanation for the longer retention times of cis-fatty acids could be that the –C=C-double bond is exposed for a longer time to the stationary phase within GC columns. Yet, it must be noted that the hdIEC effects are appreciably small with respect to the GC-MS retention times, for instance, about 0.8% for d6-metformin (Figure 2).
Protiated and perfluorinated organic solvents and derivatization reagents used in GC-MS analysis such as PFB-Br differ in their boiling points. These differences are transmitted to the derivatized analytes. The H/F effect is stronger than the H/D effect due to the stronger shielding of the analytes by F atoms from the stationary phase, a kind of “Teflon effect”. This is supported by findings by Wang and colleagues who found that benzyl, mono-fluorobenzyl, and di-fluorobenzyl derivatives have longer retention times in GC-ECD than PFB derivatives of 4-hydroxy-acetophenone (EH/F 1.0417, 1.0484, 1.0886, respectively) and an up to seven-fold less strong electron-capture response compared to PFB derivatives due to the higher number of F atoms [43].
Of the so-far-tested stationary GC columns, the ionic-liquid stationary phase IL111i is the only stationary phase that interacts with analytes mainly by π–π, dipole–dipole, and dipole–induced dipole interactions and may be the reason for the elution of deuterated analytes later than their protiated isotopologs [31].
The present work focused on GC-MS, but the chromatographic isotope effects [44] and other types of chromatographic effects such as isomerism [37,45], also apply to LC-MS even to a greater extent [3]. A method has been proposed in HILIC-LC-MS/MS metabolomics for the analysis of fatty acids derivatized with unlabeled and tetradeuterated 2-dimethylaminoethylamine and the authors of that work expected that the proposed approach will improve metabolite annotation in HILIC-MS-based metabolomics analysis [46]. Given the complexity of the chromatography processes and of the biological samples, the plethora of greatly heterogeneous endogenous and exogenous analytes, targeted (i.e., quantitative) GC-MS- and LC-MS-based metabolomics using stable isotope labeled analytes are better prepared for analytical challenges.
The particular analytical challenge for both, GC-MS and LC-MS, namely the analysis of L-amino acids (major fraction), their counterparts D-amino acids (minor fraction), and chiral secondary amino acids, has not been addressed in the present work, not because they are less important from a physiological point of view. Analysis of enantiomeric amino acids and secondary amino acids requires special analytical techniques, notably the use of special chiral columns and derivatization reagents, which are usually performed by targeted, i.e., quantitative GC-MS- and LC-MS-based metabolomics [47,48,49,50]. A potential limitation of the present work is the inclusion of a relatively small number of papers published by authors from other groups because of the lack of retention times of analytes such as polyamines [51].
Deuterium-labeled physiological substances, such as vitamins and polyunsaturated fatty acids, and non-physiological compounds, such as drugs, have been repeatedly used in in vitro and in vivo investigations [51,52,53,54]. For instance, deuterated alpha-tocopherols have been used to study by GC-MS the biokinetics and bioavailability of vitamin E and the role of the liver in their secretion in the plasma [52]. In animal studies, deuterated polyunsaturated fatty acids have been shown to exert protective effects against atherosclerosis by lowering lipid peroxidation and hypercholesterolemia [53], and to reduce hippocampal amyloid β-peptide levels [54]. The underlying mechanisms of the biological activity of deuterium-labeled substances remain elusive. GC-MS and other MS-based methods such as H/D exchange mass spectrometry [55,56] and non-MS techniques, notably NMR [57], should be useful to elucidate such “biological” H/D isotope effects.



Source link

Dimitrios Tsikas www.mdpi.com