Abstract

Mass spectrometer parameters such as Resolving Power, type of fragmentation, and mass calibration mode were optimized in the analysis of 24 pesticide metabolites in human urine using Ultra-High Pressure Liquid Chromatography coupled to Orbitrap High-Resolution Mass Spectrometer (UHPLC-HRMS). The best results were achieved with a Resolving Power of 25,000 FWHM and by applying Collision Induced Dissociation fragmentation mode (40 eV).

1. Introduction

The ever growing number of chemicals being used, such as pesticides, care products, UV filters, parabens, and so on, has an impact on the environment and therefore on humans, especially in vulnerable populations [1]. To protect consumers from such contaminants, an evaluation of the exposure to these chemicals and the subsequent risk assessment is necessary. The National Report on Human Exposure to Environmental Chemicals in the United States [2] and the Human Biomonitoring Report of Environmental Chemicals in Canada [3] are both examples of biomonitoring programs that assess the exposition of a population to environmental chemicals over time. The European Environment and Health Strategy also encourages the adoption of human biomonitoring studies across Europe [4].

Some pesticide metabolites are biomarkers of pesticide exposure. These metabolites are present in urine at concentrations of few ng mL−1 [2], and consequently metabolite determination requires sensitive and selective analytical methods. The usual analytical technique for polar metabolites in human biomonitoring studies is the LC-MS/MS [5]. However, the introduction of the high-resolution (>10,000 FWHM) mass spectrometers, which allow to obtain mass accuracies lower than 5 ppm, such as Q-TOF or Orbitrap, has allowed the implementation of combined quantitative target and postrun target analytical strategies for comprehensive determination of pesticides and other emerging contaminant metabolites [6].

Over the last few years, liquid chromatography coupled to Orbitrap high-resolution mass spectrometry (LC-HRMS) has been applied in human biomonitoring studies [6]. Recently, Cortejade et al. (2016) developed an analytical method for the targeted screening and multiresidue quantification of 38 contaminant metabolites in urine, including 12 pesticides, a pesticide metabolite (tributyl phosphate), and other compounds of different families [7]. Likewise, Roca et al. (2014) developed an analytical method that combined the quantitative target analysis of urinary metabolites of pesticides with a retrospective analysis using liquid chromatography coupled to HRMS. In this study, the main factors governing the ion-source ionization were optimized [6]. In addition, López et al. (2016) developed a retrospective analytical methodology for the analysis of pesticide metabolites in urine by LC-HRMS [8].

In the existing literature on application of HRMS to food and feed contaminants [914] and in the field of human biomonitoring [6, 1521], the emphasis so far has mainly been the detectability of the analytes based on retention time and the exact mass of the most abundant analyte ion and on quantitative determination. However, Resolving Power (R) and mass calibration should be further investigated in order to avoid interferences and improve mass accuracy. Likewise, the selection of the fragmentation mode (High-energy Collision Dissociation, HCD, or Collision Induced Dissociation, CID) and their optimization could improve the sensitivity of the confirming ions (fragments).

How much Resolving Power is necessary to apply to a specific problem should be a frequent analytical question [22]. Orbitrap mass spectrometers allow a wide range of Resolving Powers. Theoretically, a higher Resolving Power provides a better resolution of analytes from isobaric interferences present in complex matrices such as urine. However, in Orbitrap, greater Resolving Power requires longer measurement time. A higher Resolving Power leads to the monitoring of fewer data points per time unit [23]. Therefore, for quantitative analysis, the Resolving Power for each matrix-analyte combination should be selected taking into account the isobaric interferences and the number of data points required for provide good peak shape. The Resolving Power has been studied in the literature for veterinary drugs in food and animal samples [22] as well as for other residues and contaminants in food and feed [24].

In order to obtain a suitable mass accuracy (< 5 ppm) in an Orbitrap spectrometer, a proper mass calibration must be employed. External mass calibration is performed previously to the analysis by direct infusion of a mix of compounds with known masses; the experimental m/z values obtained are then corrected with the theoretical m/z values in order to fit the accuracy of the analyzer. Until now, external mass calibration in Orbitrap has been widely employed in biological samples, food, and feed [6, 7, 11, 12, 18, 20]. However, as Leendet et al. (2015) have pointed out “improper external mass calibration can lead to large systematic errors in mass measurements” and “external mass calibration range must include the mass range of interest” [25]. Another strategy is the use of internal mass calibration, which is achieved by introducing a compound with a known theoretical m/z value (lock mass) during the analysis. The comparison between the theoretical and experimental errors is used to normalize the m/z values of the rest of peaks [26]. Thus, Strano-Rossi et al. (2015) successfully determined stimulants and drugs in food supplements using internal mass calibration in an Orbitrap system [27]. Internal calibration is recommended to avoid drifts of the external calibration over time [25].

Unlike conventional quadrupole (QqQ) instruments, Orbitrap Exactive™ users do not implement a compound-specific fragmentation optimization. In Exactive™, two all ion fragmentation (AIF) modes are allowed: HCD and CID [28]. In HCD, the ions are fragmented in a collision cell using N2 as a collision gas, while CID allows the ions dissociation through interaction with neutral target species. Consequently, some authors have studied the fragmentation, optimizing the HCD cell for different substances in biological samples, food and feed [6, 11, 12, 18, 20] in an Orbitrap detector. Optimization of HRMS (Orbitrap) fragmentation mode has also been employed in areas such as proteomics [28] or for the identification of oligosaccharides [29]. Optimization of the different types of fragmentation modes (HCD and CID) will probably improve the sensitivity of these methods.

In a previous work [6], we studied the influence of the HCD collision energy on the fragmentation of various pesticide metabolites in urine, working at 50,000 FWHM. In the present paper, we want to increase the speed of the analysis using the polarity switching function (ESI+ and ESI- in the same injection) and to study more in depth other factors that can have a decisive influence on the sensitivity, accuracy and speed of the analysis of pesticide metabolites in urine. Consequently, the aims of the present work are to (i) study the influence of Resolving Power (R) on the signal intensity and mass accuracy of the different ions; (ii) compare the two all ions fragmentations (AIF) modes (HCD; CID); (iii) evaluate the two options for mass calibration (internal and external) to improve mass accuracy for this specific application.

2. Materials and Methods

This study has been developed in the framework of the DENAMIC project, which included all the required ethical approvals.

2.1. Reagents and Chemicals

Solvents were specific for pesticide residue analysis and of analytical grade. Acetonitrile and methanol were supplied by Scharlab (Barcelona, Spain). Acetic acid (purity 98-100%), β-glucuronidase arylsulfatase enzyme, and anhydrous sodium acetate were obtained from Merck (KGaA, Darmstadt, Germany). Deionized water was organically and biologically purified by using a Milli-Q Ultrapure System (Millipore, Darmstadt, Germany). QuEChERS EN extraction kits, containing 4 g MgSO4; 1 g NaCl; 1 g Sodium Citrate; 0.5g Sodium Hydrogencitrate Sesquihydrate, were obtained from Agilent Technologies (Madrid, Spain).

Standards of pesticide metabolites (Table 1) were achieved. All commercial standards were of high purity and were obtained from Dr. Ehrenstorfer (Augsburg, Germany), Sigma-Aldrich (Barcelona, Spain), Cerilliant-Certificated Reference Materials (Texas, USA), Cambridge Isotope Laboratories (Massachusetts, USA), Santa Cruz Biotechnology (Heidelberg, Germany), and AccuStandard (New Haven, USA). Stock standard solutions containing 20–500 mg L−1 of the individual compounds were prepared by weighing each compound and dissolving it in acetonitrile. Solutions were stored at −20°C. Multianalyte intermediate standard solutions were prepared by diluting the individual stock solutions in acetonitrile and used for preparing working mixed-standard solutions in acetonitrile: water (10:90, v/v). The concentration of the analytes in working solutions ranged from 1000 to 5000 ng·mL−1 depending on the compound.

2.2. Sample Preparation and UHPLC-HRMS

A previously developed sample preparation was used [6]. Briefly, in order to hydrolyze the possible glucuronide- or sulfate-conjugated metabolites, 5 mL of urine was mixed with 1 mL of 0.2 M acetate buffer 20 μL of β-glucuronidase aryl sulfatase enzyme and the internal standard solution mix in the analysis of real samples. The samples were maintained at 37°C overnight for hydrolysis.

Metabolites were extracted from the urine samples employing the dispersive solid phase extraction QuEChERS kits. In a 50 mL polypropylene tube, the urine was mixed with 10 mL of acetonitrile, a QuEChERS pouch, and 2 ceramic pieces. After centrifugation, the acetonitrile phase was transferred and evaporated to dryness in a water bath at 37°C under a nitrogen stream. Subsequently, 200 μL of methanol: water (10:90, v/v) containing 0.1% of acetic acid was added and the solution was transferred into a Millipore 0.2 μm Eppendorf and ultra-centrifuged. The final extract was transferred into an injection vial and analyzed in the UHPLC-HRMS system.

Chromatographic separation was performed on an ultra-high performance liquid chromatography (UHPLC) system Accela™ equipped with a Hypersil Gold column (100 mm x 2.1 mm, 1.9 μ) from ThermoFisher Scientific (Bremen, Germany). The chromatographic separation was optimized in a previous study [6]. Briefly, the flow rate used was 400 μL min−1 and the injection volume was 10 μL. A binary gradient was used: acetic acid 0.1% (v/v) in water was used as mobile phase A, while acetic acid 0.1% (v/v) in methanol was used as mobile phase B. The analysis started with 95% mobile phase A. After 1 min, this percentage was linearly decreased down to 45% within 5 min. After that, solvent A was decreased quickly to 0% in 0.5 min and maintained for 1.5 min. The total run time was 20 min.

Mass analysis was performed on the Orbitrap mass spectrometer Exactive™ analyzer (Thermo Scientific, Bremen, Germany). The system was equipped with a heat electrospray ionization interface (HESI-II). The ion-source parameters were previously optimized as follows: spray voltage: 3.5 kV (positive mode) and 2.5 kV (negative mode); sheath gas flow rate: 55; auxiliary gas flow rate: 10; skimmer voltage: 23 V; heater temperature: 300°C; capillary temperature: 150°C; capillary voltage: 45 V and tube lens voltage: 120 V. For more details of the HRMS analysis see Roca et al. (2014) [6].

2.3. Compound Identification Criteria

The criteria for target compound identification were established following the SANTE/11813/2017 guideline [30]: (i) mass accuracy of the molecular ion< 5 ppm; (ii) mass accuracy of the fragment ion < 5 ppm; (iii) isotopic pattern similar to the theoretical isotopic pattern (the relative intensity of the A+1 and/or A+2 isotope peaks in the real sample shall correspond to the theoretical relative intensities). For confirmation we used the reference standard solutions of those compounds available in the market. In this case the confirmation criteria were included: (iv) retention time () similar to that of the reference standard ± 0.1 min.

2.4. HRMS Orbitrap Parameters Optimization
2.4.1. Resolving Power Optimization

To optimize the Resolving Power (R), the system operating in full-scan mode (50-800 m/z) was tested at the R of 10,000; 25,000; and 50,000 FWHM. 6 blank matrix urine aliquotes spiked with a mixed-standard solution of 24 target pesticide metabolites (see Table 1) (50 ng·mL−1) were analyzed.

The Resolving Power was evaluated measuring the peak area (signal intensity) and the mass accuracy (Δm) for the diagnostic and fragment ions of each metabolite. A scheme of the Resolving Power optimization study is detailed in Table SI.1. Resolving Power optimization data were acquired with ESI+/- in separated injections, using HCD fragmentation 20 eV and external mass calibration.

2.4.2. Fragmentation Optimization

After the selection of the most suitable R, HCD and CID fragmentations were evaluated. Five spiked urine samples were injected with CID energies of 10, 20, 30, and 40 eV. We previously set the energy for HCD fragmentation to 20 eV. Once the CID energy was optimized, five different methods, in six spiked samples (50 ng·mL−1), were studied using or not HCD and CID fragmentations and using ESI+ and ESI- in the same or in different injections. The response was evaluated measuring the peak area of the fragment ions in (i) ESI+ with and without HCD (HCD= 20 eV); (ii) ESI- with and without HCD (HCD= 20 eV); (iii) ESI+ with and without CID; (iv) ESI- with and without CID; (v) ESI+ and ESI- in the same injection with and without CID. Fragmentation optimization data were acquired using the previously optimized R and external mass calibration.

2.4.3. Internal and External Mass Calibration Study

With respect to mass calibration, both external and internal mass calibrations were evaluated. External calibration was performed using the mixtures ProteoMass™ LTQ/FT-Hybrid ESI Cal Mix in Pos an Neg Mode (Supelco, Bellefonte, PA, USA). Internal mass calibration was achieved introducing caffeine (M+H+ m/z = 195.08765 Da) in the mobile phase as a lock mass for positive ionization (ESI+).

In total, five spiked aliquotes were analyzed using internal and external mass calibration separately. The analytical response was evaluated measuring peak areas and mass accuracies (Δm) for the diagnostic and fragment ions. All mass calibration study data were acquired using the previously optimized R and fragmentation settings.

Data were processed using the TraceFinder™ 3.1 (Thermo Scientific, Bremen, Germany) and Xcalibur™ 2.2 (Thermo Scientific, Bremen, Germany) software.

3. Results and Discussion

3.1. Resolving Power Optimization

In order to select the most appropriate R for the determination of pesticide metabolites in urine, the influence of this parameter on the signal and mass accuracy of the 24 compounds was investigated. Tables 2 and 3 show the signal intensity (peak area) and the mass accuracy (∆m), respectively, of the diagnostic and fragment ions.

Regarding diagnostic ions, some specific metabolites such as DEAMPY, IMPY, PNP, TCPy, MNP, DEP, DETP, or DIMET presented similar areas at the three R checked (Table 2). However, R = 25,000 FWHM provided the largest areas for DMDTP, MMP, OMET, cis DCCA, DBCA, ATZM, ALAM, METM, and 2,4,5-T. For the remaining 7 compounds a R= 25,000 FWHM presented areas close to that provided by the best R for each compound (10,000 FWHM or 50,000 FWHM).

For fragments, 17 out of 20 compounds presented higher intensities with R = 50,000 FWHM (Table 2). In some cases, as for METM, fragment ions areas are very low in comparison with the diagnostic ion response. It is relatively frequent that AIF produces fragment ions with low abundances, partly because the optimal fragmentation energy is not equal for all diagnostic ions. It could affect the analysis of real samples, because some diagnostic ions may not be confirmed because their fragment ions signal is below the noise. To sum up, R = 25,000 FWHM seems more suitable for diagnostic ions and 50,000 FWHM for fragment ions. However, taking into account the total ions (diagnostic plus fragments) and the variability of the response, a R = 25,000 FWHM was selected because it presented more number of ions with high area (see Figure SI.1). We have not found a clear explanation why some ions present higher responses when acquired at a R (e.g., 25,000 FWHM) while other ions present higher responses when acquired at another R (e.g., 50,000 FWHM).

In addition, the R of 25,000 FWHM presented more data points (scans) per peak than R = 50,000 FWHM because in an Orbitrap analyser the scan speed decrease when R increase (i.e., 2 Hz at R = 50,000 FWHM; 4 Hz at R = 25,000 FWHM). The number of data points is important for peak shape and quantification.

Table 3 shows the average and range of mass accuracy (∆m, ppm) of pesticide metabolites diagnostic and fragment ions obtained at the three R tested. Table 4 summarizes the results showing the number of ions in each of the six ranges of mass accuracy considered. As can be observed (Table 4), good mass accuracies (<3 ppm) in most of the ions were obtained applying R of 25,000 and 50,000 FWHM.

Taking into account these results, a R = 25,000 FWHM was chosen because (i) it gave the highest signal (peak area) for more analytes; (ii) it presented good mass accuracy; and (iii) the scan time gave a suitable number of data points. An added advantage of using this intermediate R is that the instrument presents a sufficient speed to be able to use the detection of positive (ESI+) and negative (ESI-) ions in the same injection, increasing the speed of analysis and the throughput of the method. Martínez-Dominguez et al. (2016) also applied R = 25,000 FWHM for the analysis of organic contaminants in food using a UHPLC-Orbitrap detector [12]. QTOF systems also employed similar Resolving Powers (18,000-22,500 FWHM) for the analysis of pesticides, organophosphate flame retardants, and chemical agents in urine and drinking water [15, 16, 31]. However, in Orbitrap analyzer a higher Resolving Power (50,000 FWHM) was used for the analysis of pesticides and drugs in urine, food, and feed [6, 11, 27] and for the analysis of mycotoxins in food (70,000 FWHM) [14]. In contrast, a lower Resolving Power was employed in other studies ranging from 7,000 to 17,500 FWHM in biological and food samples [7, 12, 18]. Kaufmann et al. (2010) compared liquid chromatography selectivity in LC-MS/MS and LC-HRMS, applying different Resolving Powers (10,000-100,000 FWHM). They indicated that a Resolving Power of 50,000 FWHM was the most suitable for the analysis of veterinary drugs in food using LC-Orbitrap detector [22].

3.2. Fragmentation Optimization

Regarding the optimization of the CID energy, the highest fragment areas were obtained using a fragmentation energy of 40 eV. As an example, Figure 1 shows the fragment ion areas obtained using CID at different energies for PBA.

Comparing the results obtained using CID and HCD, Figure 2(a) shows similar areas for fragment ions obtained with CID and HCD modes. However, some compounds such as CMHC and METM were only fragmented applying CID. Consequently, CID energy (40 eV) was selected for all compounds. Similar results were obtained using ESI+ and ESI in the same injection than in two separate injections (see Figure 2(b)). Consequently, a method with ionization mode in ESI+ and ESI in the same injection and CID fragmentation was selected. Strano-Rossi et al. (2015) also selected CID energy (40 eV) in LC-Orbitrap for the determination of stimulants and drugs in food supplements [27]. Tak et al. (2014) employed higher CID energy (150 eV) in an LC-QTOF detector for the determination of chemical warfare agents in drinking water [31]. However, HCD is the most widely fragmentation mode used in HRMS Orbitrap detectors. HCD energies from 20 to 30 eV have been applied to fragment compounds such as pesticides, biotoxins, mycotoxins, veterinary drugs, and other toxins in biological, food, and feed samples [6, 11, 12, 18, 20]. Optimization of the mode of fragmentation in Orbitrap detectors is also carried out in other fields as proteomics. Jedrychowski et al. (2011) have evaluated CID/HCD fragmentation for murine phosphoproteomics [28].

3.3. Internal and External Mass Calibration Study

In general, the use of lock mass (internal calibration) improves mass accuracy. In addition, in the present method some fragments are below 138.06619 Da, which is the low mass in the ESI+ external calibration solution; consequently a continuous correction of the acquired masses (internal calibration) could avoid an excessive mass drift.

In order to check whether this general rule is applicable to this particular application, we studied the influence of the calibration mode on the mass accuracy of eight substances analyzed in ESI+ mode, using caffeine as internal standard. Table 5 shows the mass accuracies (Δm, ppm) for the diagnostic and fragment ions obtained with internal and external mass calibration. As it can be observed, a better Δm was obtained when working with internal mass calibration. Hence, the use of caffeine as lock mass was selected.

In general, external mass calibration is used in the literature for HRMS Orbitrap detector and internal mass calibration for HRMS QTOF detectors. External mass calibration of Orbitrap in urine, plasma, food, and feed has been widely performed for pesticides, parabens, veterinary drugs, biotoxins, mycotoxins, and other substances [6, 8, 11, 12, 18, 20]. However, internal mass calibration has also been used, as in the method described by Strano-Rossi et al. (2015), for the analysis of various stimulants in food supplements. In this case the authors used diisooctyl phthalate ionic species, m/z 391.2843 Da, in order to compensate any possible mass axis drifts, obtaining ∆m < 1 ppm for all ionic species [27]. In contrast, QTOF detectors normally calibrate with negative internal mass calibration (leucine enkephalin) [15, 16, 19] or positive internal mass calibration (purine and hexakis-(1H,1H,3H-tetrafluoro-pentoxy)phosphazene) [31].

4. Conclusions

Mass spectrometry parameters such as Resolving Power, fragmentation mode, and type of mass calibration have been optimized for the analysis of 24 pesticide metabolites in urine. A Resolving Power of 25,000 FWHM, internal calibration, and CID fragmentation were selected as the best options to improve the signal intensity and the mass accuracy of diagnostic and fragment ions. This Resolving Power provides enough resolution to avoid isotopic interferences and allows the use of the polarity switching function (ESI+ and ESI in the same injection), hence reducing the analysis time.

The optimized HRMS parameters allow the determination of pesticide metabolites in urine samples; however a further validation of the method is required to determine the LOD and other performance parameters.

5. Study Limitation

We have not tested CID energies higher than 40 eV.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are very grateful to Marta Roca for her help in the experimental work. This study had the support of FP7-ENV-2011 DENAMIC Project (Code 282957).

Supplementary Materials

Table SI.1: scheme of the Resolving Power optimization study. Figure SI.1: number of ions with high (100-80%) and low (<80%) peak areas at R=10,000, R=25,000, and R=50,000 (FWHM), for total, diagnostic, and fragment ions. (Supplementary Materials)