Next Article in Journal
Characterization of Copper(II) and Zinc(II) Complexes of Peptides Mimicking the CuZnSOD Enzyme
Previous Article in Journal
Highly Sensitive Measurement of Horseradish Peroxidase Using Surface-Enhanced Raman Scattering of 2,3-Diaminophenazine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing Analysis Methods: Rapid and Accurate Determination of Cuaminosulfate Residues with LC-MS/MS Based on Box–Behnken Design Study

1
State Key Laboratory for Biology of Plant Disease and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China
2
Guangxi SPR Technology Co., Ltd., Nanning 530000, China
3
State Key Laboratory of Elemento-Organic Chemistry, National Engineering Research Center of Pesticide, College of Chemistry, Nankai University, Tianjin 300071, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2024, 29(4), 794; https://doi.org/10.3390/molecules29040794
Submission received: 6 January 2024 / Revised: 2 February 2024 / Accepted: 6 February 2024 / Published: 8 February 2024
(This article belongs to the Section Analytical Chemistry)

Abstract

:
In view of the defects in the previous detection of cuaminosulfate, which only focused on the analysis of copper ions, there is currently no analysis method available to determine the actual state of cuaminosulfate as chelated or bound. In order to investigate the dissipation and terminal residues in soil and watermelon of cuaminosulfate for food safety and environmental risk, a highly effective technique was developed to detect cuaminosulfate residues in watermelon and soil, and field experiments were conducted in China. After single-factor experiments, residual cuaminosulfate in samples was extracted by pure water, purified using a liquid–liquid approach combined with a dispersive solid-phase extraction, and detected by liquid chromatography tandem mass spectrometry (LC-MS/MS). The Box–Behnken design (BBD) study was used to find the optimal solutions for the time of liquid–liquid purification, the amount of extraction solvent, and the amounts of cleanup sorbents for the analytical method. The average recovery of the method was in the range of 80.0% to 101.1%, the average relative standard deviation (RSD) was 5.3–9.9%, and the detection limit was lower than 0.05 mg/kg. The BBD study not only improved the extraction rate of the method, but also saved time and was operated easily. The final residues of cuaminosulfate in watermelon at different sampling intervals were all lower than 0.05 mg/kg under field conditions. The cuaminosulfate in soils dissipated following exponential kinetics, with half-life values in the range of 9.39 to 12.58 days, which varied by different locations. Based on the validated method, food safety residues and soil residues can be determined rapidly and accurately.

Graphical Abstract

1. Introduction

Watermelon is one of the most important fruit crops consumed throughout the world because it provides high nutritional value and multiple health benefits to humans [1]. For example, watermelon is rich in vitamins, which balance blood pressure and regulate heart function [2]. Approximately 3 million hectares of watermelon are grown on farmland worldwide, with an estimated global production of 101 million tons per year [3]. China is the world’s largest watermelon producer, with an annual output of 60 million tons [3]. However, watermelon is susceptible to diseases such as fusarium wilt during growth, which seriously affects fruit quality [4]. Consequently, pesticides are widely applied at different growth stages of watermelon to prevent a decline in yield [4]. Nowadays, the presence of pesticide residues in different fruits and vegetables raises serious health concerns [5]. Therefore, pesticide residues in fruits, including watermelon, should be continuously monitored to ensure consumers’ food safety and minimize the health risks [6].
Copper fungicide is the oldest class of pesticides [7], and is widely used in agricultural production due to its remarkable effect on prevention and control of fungal and bacterial diseases. However, current research showed that fungicides containing copper can cause adverse damage to human health [8,9,10]. Therefore, it is essential to monitor the residues of copper fungicides in crops to assess their potential threat to human health. Cuaminosulfate (Figure 1) is a broad-spectrum protective cupric fungicide used to control watermelon wilt. It can form a tight protective film on the surface of the crop, and then slowly release copper ions on the surface of the plant. The conversion of copper ions with K+, H+, and other cations on the surface of pathogen cell membrane leads to the coagulation of pathogen cell membrane proteins. At the same time, copper ions can penetrate through the cell membrane and combine with related enzymes, inhibiting the germination of pathogenic bacteria and mycelium development and affecting the activity of pathogenic bacteria [11,12,13,14].
In the past, the monitored residues of cuaminosulfate have been defined as copper ions, which are not consistent with the actual state. For example, the method reported for detecting the dissipation dynamics of cuaminosulfate in tobacco and soil focused only on copper ions, rather than chelated or bound copper compounds, and this method is easily affected by the environment and biological copper background [15]. In recent years, the guidelines of pesticide registration and management have been continually improved. Analysis of the residue of copper ions has mostly used flame atomic absorption spectrometry [15]. Residue monitoring and the toxicity assessment of copper pesticides are no longer focused on copper ions. The detection and evaluation of copper compounds in the chelated or bound state are needed to advance the detection of residues of the pesticide in crops and to establish its maximum residue limits, which are important for human health. As the most widely used analytical tool, liquid chromatography tandem mass spectrometry (LC-MS/MS) has been reported for the analysis of pesticides in the environment and food matrix [16,17,18,19].
Cuaminosulfate is insoluble in most organic solvents and easy to decompose under acidic conditions. Due to its special physical and chemical properties, pretreatments including extraction, enrichment, and clean-up of cuaminosulfate are challenging. Moreover, cuaminosulfate is widely used in fruits and vegetables, which may cause a high risk of residues after application. Therefore, there is an urgent need to establish a high sensitivity and accurate analysis method to determine the residues of cuaminosulfate complexion in various fruits such as watermelon.
Therefore, our objective in this study was to establish a rapid and effective pretreatment procedure with an analytical method to determine the presence of residues of cuaminosulfate in a chelated or bound state in watermelon and soil, while considering the efficiency and cost optimization of sulfamic acid residue analysis. Based on the QuEChERS pretreatment method, the Box–Behnken design (BBD) of the response surface method (RSM) was used to optimize the pretreatment conditions and parameters, so as to obtain better extraction and purification effects and a higher recovery rate. Moreover, the study of residue levels of cuaminosulfate in watermelon and soil under field conditions was also undertaken. Additionally, this study provides data support for the safety evaluation of pesticide residues in soil and watermelon samples, in addition to establishing a reliable residue analysis method and theoretical base.

2. Results and Discussion

2.1. Optimization of Mass Spectrometry Conditions

In the analysis of real samples at trace levels for identifying and quantifying analytes, mass spectrometry (MS) functions are used to provide structural information for identifying non-target contaminants, increase specificity for identifying target pesticides, and achieve highly sensitive trace detection [20,21,22]. The 5 mg/L individual standard solution was directly injected into the mass spectrum, and Q1 parent ion full scanning was performed for sulfamic acid under ESI positive ion mode. To obtain its protonated [M + H]+ as a precursor ion, different fragment voltages and different collision energies were further optimized and explored to facilitate the capture of the most abundant product ions. All of the above screening methods were used to find the most sensitive transition ions in the MRM mode as quantitative ions. Table 1 summarizes the optimal values of cuaminosulfate detection under various conditions.

2.2. Optimization of Chromatography

The mobile phase is critical in LC-MS because it affects retention time, chromatographic selectivity, ionization, detection and quantitation limits, and linear range [21,22]. The effect of pH of the mobile phase on the separation of cuaminosulfate was investigated. Certain amounts of formic acid, ammonium acetate, and ammonia water were added to the aqueous phase to adjust the pH value of the mobile phase. The chromatograms of cuaminosulfate in different mobile phases are shown in Figure 2. When 10 mmol/L ammonium acetate-acetonitrile was selected as the mobile phase, the peak shape was good, but the mass spectrum response was low (Figure 2B). It is possible that the addition of ammonium acetate inhibited the ionization of cuaminosulfate, leading to a significant decrease in MS response and a quick flowing out of the peak. When 0.1% ammonia-acetonitrile (v/v) was selected as the mobile phase, the MS response was enhanced significantly due to the obvious variations in pH. However, the chromatographic peaked shape of cuaminosulfate was not satisfactory, and was not conducive to the separation and detection of compounds (Figure 2C). When 0.1% formic acid aqueous-acetonitrile (v/v) was used, the separation and peak shape of cuaminosulfate were improved, and the ionic strength also increased significantly (Figure 2A). According to the above results, 0.1% formic acid aqueous-acetonitrile (v/v) was selected as the mobile phase for separation, which can obtain a better peak shape and stronger mass spectrum intensity.

2.3. Optimization of Column

The columns selected have different separation effects. Optimizing the chromatographic method primarily involved selecting a suitable chromatographic column that could provide better retention behavior and separate target compounds from impurities [23,24]. In this study, the effects of different chromatographic columns, including Boltimate HILIC (150 mm × 21 mm, 2.7 μm), ZORBAX Eclipse Plus C18 (50 mm × 21 mm 1.8 μm), Boltimate XB-NH2 (100 mm × 30 mm, 3 μm), and Acclaim Trinity P1 (150 mm × 2.1 mm, 3 μm), on chromatographic peak shape and mass spectrum intensity were investigated. The results showed that the type of chromatographic column has a great influence on the analysis of cuaminosulfate. Boltimate XB-NH2 and Boltimate HILIC amino columns showed poor retention of cuaminosulfate and it was difficult to fully separate the impurities (Figure 3B,C). Unfortunately, the strong adsorption of the Acclaim Trinity P1 ion exchange column on cuaminosulfate led to the failure of elution from the chromatographic column, and no target chromatographic peak was observed (Figure 3D). Finally, the ZORBAX Eclipse Plus C18 (50 mm × 21 mm, 1.8 μm) chromatographic column was selected for sample separation, and a better separation effect and higher sensitivity were obtained (Figure 3A).

2.4. Optimization of Extraction and Cleanup Procedure

Watermelon samples are higher in moisture and sugar content than soil samples, and their matrix is more complex, making it more challenging to optimize analytical methods. Therefore, the optimization of analytical methods was focused on watermelon.
Since QuEChERS’s high efficiency, cost effectiveness, safety, and low cost have made it an ideal pretreatment technique for analyzing residues of pesticides in a variety of matrices, it has been widely used to analyze pesticide residues in various matrices [25,26,27,28]. To ensure efficient extraction, it was essential to select a solvent based on an improved QuEChERS pretreatment. In the present experiments, pure water, methanol, 0.1% formic acid water (v/v), and 0.1% ammonia water (v/v) were studied regarding their extraction efficiencies. To evaluate the performance of the proposed method, a recovery assay was conducted, and the average recoveries were 67.9%, 81.2%, 42.7%, 45.7%, and 25.7%, respectively, for pure water, pure water (twice), 0.1% formic acid water (v/v), 0.1% ammonia water (v/v), and methanol (Figure 4A). The results showed only pure water (twice) met the requirement of recovery (70–120%) and RSDs (<20%).
The process of liquid–liquid purification can reduce the co-extraction of insoluble impurities and reduce the difficulty of subsequent purification. Therefore, we further explored the organic solvents, ethyl acetate, n-hexane, acetonitrile, and dichloromethane for liquid–liquid extraction. According to the results shown in Figure 4B, dichloromethane met the requirements of both RSD and recovery. The recoveries of other solvents were below 60%.
Purification agents are used to eliminate the organic impurities in co-extraction and achieve one-step purification. Fruits contain various sugars and carotenoids, which could dissolve in water. Multiple sorbents were selected to optimize the purification process due to their different characteristics. For example, PSA is suitable for the removal of various polar organic acids, polar pigments, and certain sugars. C18 has a strong adsorption effect on non-polar compounds [27,29,30]. Florisil is a weakly polar adsorbent, retaining less polarity as well as fewer polar compounds [31]. Four sorbents, namely, 50 mg PSA + 150 mg MgSO4, 50 mg C18 + 150 mg MgSO4, 50 mg Florish + 150 mg MgSO4, and 50 mg GCB + 150 mg MgSO4, were compared in this study. According to the results shown in Figure 4C, recoveries were 62.3%, 80.5%, 47.8%, and 33.3% for 50 mg PSA + 150 mg MgSO4, 50 mg C18 + 150 mg MgSO4, 50 mg Florish + 150 mg MgSO4, and 50 mg GCB + 150 mg MgSO4. Therefore, 50 mg C18 + 150 mg MgSO4 was used as the sorbent for the following procedure.

2.5. BBD Studies

As a result of single-factor experiments, different extraction solvents (pHs), liquid–liquid purifications, and cleanup procedures influenced recovery. Nevertheless, the efficiency of the pretreatment was the result of a combination of multiple factors, and it is not comprehensive to consider a single factor separately [32]. In order to obtain the optimal solution of QuEChERS conditions, we developed a three-variable, three-level BBD experiment on the volume of extraction solvent, the time of liquid–liquid purification, and the amount of cleanup sorbents (Table S2). The experiments were designed with 15 tests and 3 central points (Table S3) utilizing the software Design Expert 13.0. The results of the data analysis of BBD indicated that a quadratic regression model fits well with the data (R2 = 0.9766). The analysis of variance for the chosen model is shown in Table S4, which includes coefficients, F values, and p values. The F value for the “lack of fit” is not significant at the 95% confidence level, indicating that the model was able to accurately predict the response variables (Table S5).
Figure 5 shows a 3D response surface graph, which represents the regression equation diagram.
R = −178.48333 + 4.83417A + 42.97500B + 3.50700C + 0.115000A × B + 0.003800A × C + 0.082000B × C − 0.110667A2 − 7.61667B2 − 0.023387C2
where R is recovery, and A, B, and C are the volume in milliliters of extraction solvent, time of liquid–liquid purification, and milligrams of cleanup sorbents, respectively.
From 2D response surface graphs, the interaction between two factors can be found by the shape of the contour map. The shape in Figure 5D is a circle, which indicates the interaction between the volume of extraction solvent and the amount of cleanup sorbents was not significant. Alternatively, the shapes in Figure 5B,F are ovals, demonstrating that the interactions between liquid–liquid purification and amount of cleanup sorbents, volume of extraction solvent, and time of liquid–liquid purification were significant.
Table S6 displays the constraints on independent and dependent variables. In order to obtain the highest recovery, the recovery value was set to 100. Cleanup sorbents can reduce matrix effects when a sufficient quantity of sorbents is used, and the goal of factor C was set to “maximize”. Since it also impacted the response value the most (lowest p value), the relevance of factor C was set at 5. Meanwhile, for sake of conserving material and solvent, factors A (volume of extraction solvent) and B (time of liquid–liquid purification) were set to “minimize”. The optimum solution of BBD is displayed in Table S7, and shows that the ideal QuEChERS method for cuaminosulfate in watermelon samples was found to be 20 mL of pure water as the extraction solvent, liquid–liquid purification repeated three times, and 80 mg C18 as the cleanup sorbent with a 10 g sample.

2.6. Method Validation

In this study, the linearity, sensitivity, precision, and accuracy of the analytical method were validated, and the results are shown in Table 2 and Table 3. Good linearity was obtained in the range of 0.005~0.5 mg/L for cuaminosulfate, with R2 > 0.999 (Table 2).
In this experiment, ME was calculated based on the slope of 0.1% ammonium hydroxide and the matrix-matched calibration curves. The results showed that ME of watermelon and soil was −93.8% and −56.6%. ME means that other components in the matrix will interfere with the quantification of target analyte when they are co-ionized [32]. When |ME| ≥ 50%, ME was strong [33,34]. Therefore, the influence of the substrate on the results cannot be ignored, and the calibration should be performed for cuaminosulfate according to the external matrix-matched standards to eliminate the ME and obtain more accurate quantification in all samples.
The precision and accuracy of cuaminosulfate in watermelon and soil are shown in Table 3. The intra-day average recovery of cuaminosulfate in watermelon samples was 79.6~95.1%, and inter-day average recovery was 88.3~101.1%. The RSDr ranged from 5.3 to 9.9%, and RSDR ranged from 3.5 to 9.4%, indicating that the method had high accuracy and precision. For soil samples, the intra-day and inter-day recoveries were 80.8~93.1% and 81.3~93.4%, and the RSDs were 3.4~7.8% and 4.1~7.3%, respectively, which all satisfied the requirements. The recovery results could verify the reliability of the optimal solution from the BBD study, which was a practical operation for optimizing the parameters.
After assessment of AGREEprep, a somewhat high score of 0.58 was obtained (Figure S1A). Although the preparation procedure is a green technique, some particular procedures had a few drawbacks, namely, ex situ mode and sample throughput, more waste, and the use of HPLC-MS/MS. Moreover, the blue applicability grade index (BAGI) was used as a metric tool to evaluate the practicability of method; the score was 65.0, which is higher than 60 (Figure S1B). This demonstrated that this method was considered to be practical.

2.7. Dissipation and Terminal Residue of Cuaminosulfate under Field Conditions

The residues of cuaminosulfate in whole watermelon at different sampling points were determined after the last application under field conditions. The average residual values of the original concentration of cuaminosulfate in watermelon were all below 0.05 mg/kg, so the dissipation kinetic fitting could not be performed. The results of final residue tests at 10 sites in 2019 showed that the median residual value (STMR) and maximum residual value (HR) of cuaminosulfate in watermelon were 0.05 mg/kg at 20, 30, and 40 days after the last application. At different harvesting intervals, the residual concentrations of cuaminosulfate in watermelon were all <0.05 mg/kg. It was also proved that there was a low risk of cuaminosulfate residues when watermelon was treated with root irrigation at the recommended dose.
To better understand the residues of cuaminosulfate in soil, the dissipation studies were undertaken in the provinces of Anhui, Guangxi, Shanxi, and Fujian in China. The initial residues of cuaminosulfate in soils ranged from 4.0 to 10.2 mg/kg, as shown in Table 4. The dissipation kinetics of cuaminosulfate in soils fitted well with the first-order kinetic equation in Fujian, Guangxi, and Anhui (R2 = 0.7615~0.8856). In Fujian, cuaminosulfate degraded from 9.1 to 7.4 mg/kg with a degradation rate of 77.3% in 7 days, and finally decreased to 91.6% at the end of the observation points at 28 days. In Guangxi and Guangdong, the cuaminosulfate degraded to approximately 50% at 7 days after the last application, and finally degraded to 85% and 73.7%. According to the first-order kinetic equation, the DT50 value of cuaminosulfate was 9.4, 12.3, and 12.6 days in Fujian, Guangxi, and Guangdong. It was demonstrated that the degradation rate was different among different locations. In Shanxi, the residues of cuaminosulfate were degraded to approximately 50% of the initial point, and then residues did not show significant differences among the following sampling points. The degradation of pesticides in the matrix may be closely related to different climatic conditions and properties of soil [35,36]. The cuaminosulfate degraded faster in Fujian compared with in other soils, which may relate to pH.

3. Materials and Methods

3.1. Reagents and Chemicals

Standard cuaminosulfate (purity 85.1%) and 15% cuaminosulfate water were obtained from Shandong United Pesticide Industry Co., Ltd. (Jinan, China) and Tianjin Huayu Pesticide Co., Ltd. (Tianjin, China) Acetonitrile (LC-MS grade) was bought from Thermo Fisher SCIENTIFIC (Shanghai, China). LC-grade formic acid was supplied by Tianjin Kermel Chemical Reagent Co., Ltd. (Tianjin, China). Methanol and n-hexane were purchased from Autogenic Bioengineering (Shanghai) Co., Ltd. (Shanghai, China) Ethyl acetate was obtained from Beijing Jinruilin Technology Development Co., Ltd. (Beijing, China) Analytical ammonium hydroxide (NH3·H2O) and dichloromethane were purchased from Xilong Scientific Co., Ltd. (Guangzhou, China) and Tianjin Fuyu Fine Chemical Co., Ltd. (Tianjin, China), respectively. Ultra-pure grade water was obtained from Hangzhou Wahaha Group Co., Ltd. (Hangzhou, China). Primary secondary amine (PSA) and octadecylsilane (C18) were obtained from Welch Materials, Inc. (Shanghai, China).

3.2. Standard Solution Preparation

The stock solution of cuaminosulfate (1000 mg/L) was prepared in 0.1% (v/v) ammonium hydroxide solution. The stock solution was stored in a refrigerator at 4 °C in the dark and remained stable for up to 3 months. Standard working solutions at 0.01, 0.05, 0.1, 0.2, 0.4, and 0.5 mg/L were prepared by taking appropriate aliquots of the stock solution and diluting them serially to volume in 0.1% (v/v) ammonium hydroxide solution. Correspondingly, blank sample extracts were used as solvents to prepare matrix-matched standard solutions with concentrations of 0.01, 0.05, 0.1, 0.2, 0.4, and 0.5 mg/L, respectively.

3.3. Chromatographic and Mass Spectrometric Conditions

Chromatographic analysis was conducted using a Waters Acquity UPLC system (Waters, Milford, MA, USA) coupled to an API 4000 triple quadrupole mass spectrometer (Sciex, Foster City, CA, USA). The separation was achieved using a ZORBAX, Eclipse Plus C18 column (2.1 mm × 50 mm, with particle size of 1.8 μm) (Agilent Technologies, Shanghai, China) at a flow rate of 0.2 mL/min. A 2 μL volume of target solution was injected into the LC–MS/MS system, and the total analysis time was 6.0 min. The mobile phase consisted of 0.1% (v/v) formic acid in water (solvent A) and acetonitrile (solvent B). Elution conditions were as follows: 0 to 1.0 min, 10% solvent B; 1.0 to 2.0 min, a linear gradient change from 10% to 90% B; 2.0 to 4.0 min, 90% B; 4.0 to 4.1 min, a linear gradient change from 90% to 10% B; and 4.1 to 6.0 min, the mobile phase returned to initial conditions. The column temperature was maintained at 35 °C.
MS detection was operated in positive ion mode. Multiple Reaction Monitoring (MRM) mode was performed on the protonated molecule for cuaminosulfate using the following optimized parameters: ion spray voltage (v): 5500; source temperature (°C): 500; curtain gas pressure (psi), 30; nebulizer gas pressure (psi), 55; auxiliary gas pressure (psi): 55; desolvation temperature (°C): 500; sheath gas pressure (psi): 9. MRM transitions and other acquisition parameters for cuaminosulfate are given in Table 1.

3.4. Sample Preparation

Extraction process: A quantity of 10.0 g of watermelon or soil was weighed after being ground and homogenized, then transferred to a 50 mL plastic centrifuge tube with 20 mL ultrapure water. Using a vortex mixer, the tube was vigorously shaken for 2 min. Afterward, the mixture was centrifuged at 5000 rpm for 3 min; then, the supernatant was placed in a graduated tube. The above extraction process was repeated once and the supernatant was mixed. Finally, the supernatant was diluted with ultrapure water to a final volume of 50 mL.
Liquid–liquid purification: A quantity of 4 mL of the extraction sample was transferred to a centrifuge tube, and 4 mL of dichloromethane was added. The mixture was vortexed for 2 min and centrifuged for 5 min at 4000 rpm. The above process was repeated twice and the inorganic layer was mixed.
Purified process: A quantity of 2 mL aqueous phase from the above process was transferred to a centrifuge tube containing 100 mg C18 and 20 mg PSA. The mixture was vortexed for 2 min and centrifuged for 3 min at 6000 rpm. The solution was filtered through a 0.22 μm water membrane prior to analysis.

3.5. Field Trial Design and Sample Collection

A field trial study was conducted at 10 sites representative of China in 2019 according to Pesticide Residue Testing Guidelines (NY/T788–2018) issued by Ministry of Agriculture and Rural Affairs of the People’s Republic of China. The 10 sites were located in Shandong, Hunan, Guangxi, Hebei, Guizhou, Shanxi, Henan, Zhejiang, Anhui, and Fujian, as shown in Figure 6. One treatment plot (50 m2) and one control plot (50 m2) were designed for each place. A conservation area of 20 m2 was set to separate each plot in the same field. There had never been a field trial with cuaminosulfate conducted previously. The background information is shown in Table S1.
Each plot was applied with 0.25 g active indigent/plant three times through root irrigation. Samples of watermelon and soil were collected randomly at 0, 3, 10, 20, 30, and 40 days after the last application. More than 12 disease-free and mature watermelon fruits with a total weight of at least 2 kg and soil around the collected plants were collected at each treated plot. The samples were transported to the laboratory immediately after collection. After being chopped and thoroughly mixed, samples were kept in a deep freeze (−20 °C) before analysis.

3.6. Method Validation

The analytical method was validated by evaluating the linearity, limit of quantification (LOQ), accuracy, precision, selectivity, and matrix effect. The specificity of the method was assessed by extracting and analyzing the blank samples of watermelon and soil.
The linearity was assessed by the correlation coefficient (R2) of working standard solutions with six concentrations in the range of 0.01, 0.05, 0.1, 0.2, 0.4, and 0.5 mg/L, which were obtained from standard stock solutions diluted serially by 0.1% ammonium hydroxide solution (v/v) and different matrix solutions from blank samples. The selectivity of the method was evaluated by LOQ, which was defined as the minimum spiked concentration in matrix.
Recovery tests were conducted to verify the accuracy and precision of the method. Five replications were prepared at four different spiked concentrations. The precision of the method was evaluated by intra-day (repeatability) and inter-day (reproducibility) analysis. The intra-day precision, which was used in the analysis of the spiked samples in five replicates in one day, and was calculated by the percentage of RSD (relative standard deviation), is represented by RSDr. The inter-day (reproducibility) analysis was conducted with the same spiked samples (five replicates) on three different days by following all the identical experimental and instrumental conditions, and the result is denoted RSDR.
Meanwhile, the blue applicability grade index (BAGI) and greenness of analysis were also estimated in this study to meet sustainability requirements. The analytical greenness metric for sample preparation (AGREEprep) [37] and the blue applicability grade index (BAGI) metric tool [38] were selected as the tools for the evaluation of greenness and blueness analysis.

3.7. Statistical Analysis

3.7.1. Matrix Effect and Recovery

The matrix effect was calculated as below [39]:
ME % = (slopem − slopes)/slopes × 100%
where slopem represents the slope of the matrix-matched standard curve, and slopes is the slope of the solvent standard curve.
The recovery was calculated as the ratio between the peak area of the target compound in sample extracts and the peak area of a corresponding matrix-matched calibration solution. The equation is as below [40]:
Recovery (%) = Amount of pesticide obtained/Amount of pesticide added × 100%

3.7.2. Dissipation Kinetics

The dissipation kinetics of cuaminosulfate in soil samples were determined using first-order kinetics according to GB/T 31270.1-2014 [41]. The equations are as follows:
Ct = C0 × e−kt
T1/2 = ln2/k
where Ct is the residue (mg/kg) at time t (days), C0 is the initial residue (mg/kg), k represents a first-order kinetic degradation rate constant, and T1/2 is the half-life [42].

4. Conclusions

In this study, a QuEChERS method based on HPLC-MS/MS for the determination of cuaminosulfate residues in watermelon was established for the first time. The influence of a single factor, i.e., extraction solvent, liquid–liquid extraction, and cleanup sorbent, on the pretreatment process was assessed, and then a BBD study of the response surface methodology determined the efficiency of combining multiple factors to obtain an optimal solution. With the BBD study, we are able to obtain a comprehensive evaluation of the optimal parameters for the multivariable pretreatment conditions. The developed method was successfully applied to determination of cuaminosulfate in watermelon and soil samples. Good linearity and satisfactory recoveries and reproducibility were obtained. The average recoveries ranged from 80.0% to 101.1%, the RSD was 5.3–9.9%, and the detection limit was lower than 0.05 mg/kg. Linearity for target compounds with satisfactory correlation coefficients (R2 > 0.999) was obtained in all of the matrices. Meanwhile, after evaluating the greenness and blueness of the analytical method, the score showed the method can be recommended and is practical. The results of dissipation and final residue studies showed that the residue of cuaminosulfate in watermelon was <0.05 mg/kg. The cuaminosulfate degraded fast in watermelon. In soil, cuaminosulfate dissipated following exponential kinetics with a half-life value of 9.39 to 12.58 days, which varied by location. This method was fast and reliable and should be used for checking the presence of residues in watermelon and soils, and will also help in establishing maximum residue limits (MRLs) of cuaminosulfate in watermelon in China. In future, we will investigate the effects of soil properties on the degradation rate, such as pH and soil texture.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29040794/s1, Table S1: The field trial background information; Table S2: Different parameter and variable factor Settings for BBD; Table S3: Observation and prediction of recovery of cuaminosulfate in watermelon by BBD in 15 cycle designs; Table S4: Analysis of variance (ANOVA) of the established model; Table S5: ANOVA for quadratic model of recovery of cuaminosulfate in watermelon; Table S6: The constraints of independent variables and dependent variable; Table S7: The solution of BBD; Table S8: The supplement information of chemical substances; Figure S1: The results of AGREEprep assessment (A) and BAGI index pictograms (B).

Author Contributions

Conceptualization, M.H. and X.L.; Methodology, M.H., X.L. and C.W.; Validation, L.Z. (Lan Zhang), L.M., L.Z. (Lizhen Zhu) and Y.Z.; Formal analysis, L.M.; Investigation, Y.W.; Resources, M.H.; Data curation, Y.W. and L.Z. (Lizhen Zhu); Writing—original draft, Y.W. and L.Z. (Lizhen Zhu); Writing—review & editing, L.Z. (Lan Zhang), Y.Z. and X.L.; Supervision, C.W.; Project administration, M.H., Y.W. and C.W.; Funding acquisition, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by The National Natural Science Foundation of China (32202367).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

Author Mingyuan He was employed by the company Guangxi SPR Technology Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Sorokina, M.; McCaffrey, K.S.; Deaton, E.E.; Ma, G.Y.; Ordovás, J.M.; Perkins-Veazie, P.M.; Steinbeck, C.; Levi, A.; Parnell, L.D. A Catalog of Natural Products Occurring in Watermelon—Citrullus lanatus. Front. Nutr. 2021, 8, 729822. [Google Scholar] [CrossRef] [PubMed]
  2. Manivannan, A.; Lee, E.S.; Han, K.; Lee, H.E.; Kim, D. Versatile Nutraceutical Potentials of Watermelon—A Modest Fruit Loaded with Pharmaceutically Valuable Phytochemicals. Molecules 2020, 25, 5258. [Google Scholar] [CrossRef] [PubMed]
  3. FAO. World Food and Agriculture—Statistical Yearbook 2022; FAO: Rome, Italy, 2022. [Google Scholar] [CrossRef]
  4. Al-Shamary, N.M.; Al-Ghouti, M.A.; Al-Shaikh, I.; Al-Meer, S.H.; Ahmad, T.A. Evaluation of pesticide residues of organochlorine in vegetables and fruits in Qatar: Statistical analysis. Environ. Monit. Assess. 2016, 188, 198. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, M.; Li, X.; Han, L.; Wang, Q.; Kong, X.; Xu, M.; Wang, K.; Xu, H.; Shen, Y.; Gao, G.; et al. Determination and risk assessment of 31 pesticide residues in apples from China’s major production regions. J. Food Compos. Anal. 2023, 118, 105188. [Google Scholar] [CrossRef]
  6. Mahdavi, V.; Eslami, Z.; Molaee-Aghaee, E.; Peivasteh-Roudsari, L.; Sadighara, P.; Van Nam, T.; Fakhri, Y.; Ravanlou, A.A. Evaluation of pesticide residues and risk assessment in apple and grape from western Azerbaijan Province of Iran. Environ. Res. 2022, 203, 111882. [Google Scholar] [CrossRef] [PubMed]
  7. Lamichhane, J.R.; Osdaghi, E.; Behlau, F.; Köhl, J.; Jones, J.B.; Aubertot, J.N. Thirteen decades of antimicrobial copper compounds applied in agriculture. A review. Agron. Sustain. Dev. 2018, 38, 28. [Google Scholar] [CrossRef]
  8. McNeil, E.M. Antifouling: Regulation of biocides in the UK before and after Brexit. Mar. Policy 2018, 92, 58–60. [Google Scholar] [CrossRef]
  9. Capinera, J.L.; Dickens, K. Some effects of copper-based fungicides on plant-feeding terrestrial molluscs: A role for repellents in mollusc management. Crop Prot. 2016, 83, 76–82. [Google Scholar] [CrossRef]
  10. Otter, J.; Cantrell, F.L. Death by Copper Ammonium Complex-Containing Fungicide. Am. J. Foren. Med. Path 2021, 42, 62–63. [Google Scholar] [CrossRef]
  11. Celar, F.A.; Kos, K. In vitro compatibility of strain ATCC 74040 with copper fungicides. J. Appl. Entomol. 2020, 144, 388–395. [Google Scholar] [CrossRef]
  12. Giller, K.E.; Witter, E.; McGrath, S.P. Toxicity of heavy metals to microorganisms and microbial processes in agricultural soils: A review. Soil. Biol. Biochem. 1998, 30, 1389–1414. [Google Scholar] [CrossRef]
  13. Knauert, S.; Knauer, K. The role of reactive oxygen species in copper toxicity to two freshwater green algae. J. Phycol. 2008, 44, 311–319. [Google Scholar] [CrossRef]
  14. Tamm, L.; Thuerig, B.; Apostolov, S.; Blogg, H.; Borgo, E.; Corneo, P.E.; Fittje, S.; de Palma, M.; Donko, A.; Experton, C.; et al. Use of Copper-Based Fungicides in Organic Agriculture in Twelve European Countries. Agronomy 2022, 12, 673. [Google Scholar] [CrossRef]
  15. Li, F.; Wei, Y.; Fan, Y.; Li, Z. Degradation Dynamics of Cuaminosul-Fate in Tobacco and Soil. Henan Sci. 2016, 34, 1444–1447. [Google Scholar]
  16. Syrgabek, Y.; Alimzhanova, M. Modern Analytical Methods for the Analysis of Pesticides in Grapes: A Review. Foods 2022, 11, 1623. [Google Scholar] [CrossRef]
  17. González-Curbelo, M.Á.; Varela-Martínez, D.A.; Riaño-Herrera, D.A. Pesticide-residue analysis in soils by the QuEChERS method: A review. Molecules 2022, 27, 4323. [Google Scholar] [CrossRef]
  18. Wahab, S.; Muzammil, K.; Nasir, N.; Khan, M.S.; Ahmad, M.F.; Khalid, M.; Ahmad, W.; Dawria, A.; Reddy, L.K.; Busayli, A.M. Advancement and new trends in analysis of pesticide residues in food: A comprehensive review. Plants 2022, 11, 1106. [Google Scholar] [CrossRef]
  19. Picó, Y.; Blasco, C.; Font, G. Environmental and food applications of LC–tandem mass spectrometry in pesticide-residue analysis: An overview. Mass Spectrom. Rev. 2004, 23, 45–85. [Google Scholar] [CrossRef] [PubMed]
  20. Dong, F.; Chen, X.; Liu, X.; Xu, J.; Li, Y.; Shan, W.; Zheng, Y. Simultaneous determination of five pyrazole fungicides in cereals, vegetables and fruits using liquid chromatography/tandem mass spectrometry. J. Chromatogr. A 2012, 1262, 98–106. [Google Scholar] [CrossRef] [PubMed]
  21. Li, Y.; Dong, F.; Liu, X.; Xu, J.; Li, J.; Kong, Z.; Chen, X.; Liang, X.; Zheng, Y. Simultaneous enantioselective determination of triazole fungicides in soil and water by chiral liquid chromatography/tandem mass spectrometry. J. Chromatogr. A 2012, 1224, 51–60. [Google Scholar] [CrossRef] [PubMed]
  22. Brookhart, A.; Arora, M.; McCullagh, M.; Wilson, I.D.; Plumb, R.S.; Vissers, J.P.C.; Tanna, N. Understanding mobile phase buffer composition and chemical structure effects on electrospray ionization mass spectrometry response. J. Chromatogr. A 2023, 1696, 463966. [Google Scholar] [CrossRef] [PubMed]
  23. Lambert, N.; Miyazaki, S.; Ohira, M.; Tanaka, N.; Felinger, A. Comparison of the kinetic performance of different columns for fast liquid chromatography, emphasizing the contributions of column end structure. J. Chromatogr. A 2016, 1473, 99–108. [Google Scholar] [CrossRef] [PubMed]
  24. Blumberg, L.M. Methodology of quantitative comparison of practically achievable kinetic performance of differently structured liquid chromatography columns. J. Chromatogr. A 2022, 1672, 463039. [Google Scholar] [CrossRef] [PubMed]
  25. Lehotay, S.J.; Son, K.A.; Kwon, H.; Koesukwiwat, U.; Fu, W.S.; Mastovska, K.; Hoh, E.; Leepipatpiboon, N. Comparison of QuEChERS sample preparation methods for the analysis of pesticide residues in fruits and vegetables. J. Chromatogr. A 2010, 1217, 2548–2560. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, Y.B.; Qiao, H.O.; Chen, C.; Wang, Z.L.; Xia, X.D. Determination of nitrofurans metabolites residues in aquatic products by ultra-performance liquid chromatography-tandem mass spectrometry. Food Chem. 2016, 192, 612–617. [Google Scholar] [CrossRef] [PubMed]
  27. Xiu-Ping, Z.; Lin, M.; Lan-Qi, H.; Jian-Bo, C.; Li, Z. The optimization and establishment of QuEChERS-UPLC-MS/MS method for simultaneously detecting various kinds of pesticides residues in fruits and vegetables. J. Chromatogr. B 2017, 1060, 281–290. [Google Scholar] [CrossRef]
  28. Ly, T.K.; Ho, T.D.; Behra, P.; Nhu-Trang, T.T. Determination of 400 pesticide residues in green tea leaves by UPLC-MS/MS and GC-MS/MS combined with QuEChERS extraction and mixed-mode SPE clean-up method. Food Chem. 2020, 326, 126928. [Google Scholar] [CrossRef]
  29. Hu, J.Y.; Li, J.Z. Determination of forchlorfenuron residues in watermelon by solid-phase extraction and high-performance liquid chromatography. J. AOAC Int. 2006, 89, 1635–1640. [Google Scholar] [CrossRef]
  30. Wang, J.; Duan, H.L.; Ma, S.Y.; Zhang, J.; Zhang, Z.Q. Solidification of a Switchable Solvent-Based QuEChERS Method for Detection of 16 Pesticides in Some Fruits and Vegetables. J. Agric. Food Chem. 2019, 67, 8045–8052. [Google Scholar] [CrossRef] [PubMed]
  31. Manirakiza, P.; Covaci, A.; Schepens, P. Improved analytical procedure for determination of chlorinated pesticide residues in human serum using solid phase disc extraction (SPDE), single-step clean-up and gas chromatography. Chromatographia 2002, 55, 353–359. [Google Scholar] [CrossRef]
  32. Han, X.; Huang, J.; Cai, J.; Yu, C.; Jiang, S.; Ye, Q.; Yang, S.; Pan, W.; Zhang, Q.; Wang, Y. Microwave-ultrasonic technique development coupled with natural deep eutectic solvents in anthocyanin extraction from perilla leaves (Perilla frutescens var. Acuta). J. Sci. Food Agric. 2023, 103, 3006–3016. [Google Scholar] [CrossRef] [PubMed]
  33. Jin, X.; Kaw, H.Y.; Zhao, J.; Zou, Y.; He, M.; Li, Z.; Li, D. NLow matrix effect pretreatment method based on gas-liquid micro-extraction technique for determining multi-class pesticides in crops. J. Chromatogr. A 2022, 1675, 463178. [Google Scholar] [CrossRef] [PubMed]
  34. Williams, M.L.; Olomukoro, A.A.; Emmons, R.V.; Goodage, N.H.; Gionfriddo, E. Matrix effects demystified: Strategies for resolving challenges in analytical separations of complex samples. J. Sep. Sci. 2023, 46, e2300571. [Google Scholar] [CrossRef] [PubMed]
  35. Hou, X.; Qiao, T.; Zhao, Y.; Liu, D. Dissipation and safety evaluation of afidopyropen and its metabolite residues in supervised cotton field. Ecotoxicol. Environ. Saf. 2019, 180, 227–233. [Google Scholar] [CrossRef] [PubMed]
  36. Su, Y.; Wang, W.; Hu, J.; Liu, X. Dissipation behavior, residues distribution and dietary risk assessment of tembotrione and its metabolite in maize via QuEChERS using HPLC-MS/MS technique. Ecotoxicol. Environ. Saf. 2020, 191, 110187. [Google Scholar] [CrossRef] [PubMed]
  37. Wojnowski, W.; Tobiszewski, M.; Pena-Pereira, F.; Psillakis, E. AGREEprep—Analytical Greenness Metric for Sample Preparation. TrAC Trends Anal. Chem. 2022, 149, 116553. [Google Scholar] [CrossRef]
  38. Natalia, M.; Wojciech, W.; Justyna, P.; Victoria, S. Blue applicability grade index (BAGI) and software: A new tool for the evaluation of method practicality. Green Chem. 2023, 25, 7598. [Google Scholar]
  39. Jankowska, M.; Lozowicka, B. The processing factors of canning and pasteurization for the most frequently occurring fungicides and insecticides in apples and their application into dietary risk assessment. Food Chem. 2022, 371, 131179. [Google Scholar] [CrossRef]
  40. Li, M.M.; Liu, Y.N.; Fan, B.; Lu, J.; He, Y.; Kong, Z.Q.; Zhu, Y.L.; Jian, Q.; Wang, F.Z. A chemometric processing-factor-based approach to the determination of the fates of five pesticides during apple processing. Lwt-Food Sci. Technol. 2015, 63, 1102–1109. [Google Scholar] [CrossRef]
  41. GB/T 31270.1-2014; Test Guidelines for Environmental Safety Evaluation of Chemical Pesticides Part 4: Soil Adsorption/Deabsorption. Standardization Administration of the People’s Republic of China: Beijing, China, 2014.
  42. Zhao, J.L.; Tan, Z.H.; Wen, Y.; Fan, S.; Liu, C.L. Dissipation of fluazinam in citrus groves and a risk assessment for its dietary intake. J. Sci. Food Agric. 2020, 100, 2052–2056. [Google Scholar] [CrossRef]
Figure 1. The chemical structure of cuaminosulfate.
Figure 1. The chemical structure of cuaminosulfate.
Molecules 29 00794 g001
Figure 2. Chromatograms of cuaminosulfate with different mobile phases: (A) 10 mmol/L ammonium acetate-acetonitrile; (B) 0.1% (v/v) formic acid aqueous-acetonitrile; (C) 0.1% (v/v) ammonia-acetonitrile.
Figure 2. Chromatograms of cuaminosulfate with different mobile phases: (A) 10 mmol/L ammonium acetate-acetonitrile; (B) 0.1% (v/v) formic acid aqueous-acetonitrile; (C) 0.1% (v/v) ammonia-acetonitrile.
Molecules 29 00794 g002
Figure 3. Chromatograms of chromatographic columns: (A) ZORBAX Eclipse Plus C18; (B) Boltimate HILIC; (C) Boltimate XB-NH2; and (D) Acclaim Trinity P1.
Figure 3. Chromatograms of chromatographic columns: (A) ZORBAX Eclipse Plus C18; (B) Boltimate HILIC; (C) Boltimate XB-NH2; and (D) Acclaim Trinity P1.
Molecules 29 00794 g003
Figure 4. Recovery result of single-factor experiments for pretreatment procedure: (A) extraction solvent; (B) liquid–liquid purification; (C) cleanup sorbents.
Figure 4. Recovery result of single-factor experiments for pretreatment procedure: (A) extraction solvent; (B) liquid–liquid purification; (C) cleanup sorbents.
Molecules 29 00794 g004
Figure 5. Response surfaces for Methomyl using the BBD obtained by plotting (A,B) purification time vs. extraction solvent volume, (C,D) cleanup sorbent mass vs. extraction solvent volume, and (E,F) cleanup sorbent mass vs. purification time.
Figure 5. Response surfaces for Methomyl using the BBD obtained by plotting (A,B) purification time vs. extraction solvent volume, (C,D) cleanup sorbent mass vs. extraction solvent volume, and (E,F) cleanup sorbent mass vs. purification time.
Molecules 29 00794 g005
Figure 6. Terminal residues and dissipation dynamics of cuaminosulfate in watermelon and soil.
Figure 6. Terminal residues and dissipation dynamics of cuaminosulfate in watermelon and soil.
Molecules 29 00794 g006
Table 1. Experimental parameters and MS conditions of cuaminosulfate in ESI+ mode.
Table 1. Experimental parameters and MS conditions of cuaminosulfate in ESI+ mode.
MWtR (min)DP (V)Quantification Ion TransitionCE1 (eV)Diagnostic Ion TransitionCE2 (eV)
131.670.5167/65131.9→104.010131.9→122.06
MW: Molecular weight; tR: retention time; DP: declustering potential; CE: collision energy.
Table 2. The calibration curves, determination coefficient (R2), and matrix effect of cuaminosulfate.
Table 2. The calibration curves, determination coefficient (R2), and matrix effect of cuaminosulfate.
MatrixCalibration CurveR2Matrix Effect (%)
0.1% Ammonium hydroxidey = 1,699,198.0640x + 19,444.51770.9991-
Watermelon y = 105,377.7778x + 4980.27780.9998−93.8
Soily = 736,731.1573x +5217.13440.9981−56.6
Table 3. Recovery, relative standard deviation (RSD), and LOQ of cuaminosulfate in watermelon and soil matrix.
Table 3. Recovery, relative standard deviation (RSD), and LOQ of cuaminosulfate in watermelon and soil matrix.
SampleSpiked Level (mg/kg)Intra-Day
Recovery (%)
RSDr (%)Inter-Day
Recovery (%)
RSDR (%)LOQ (mg/kg)
Watermelon0.0595.15.388.36.30.05
0.579.65.493.29.4
180.09.989.44.9
288.36.5101.13.5
Soil0.0593.16.392.47.30.05
0.588.57.487.54.1
580.87.883.94.9
2081.13.4801.36.7
Table 4. Half-life (T1/2) for cuaminosulfate dissipation in different provinces’ field soil conditions.
Table 4. Half-life (T1/2) for cuaminosulfate dissipation in different provinces’ field soil conditions.
FujianGuangxiShanxiAnhui
Collection PointsResidue (mg/kg)Degradation RateResidue (mg/kg)Degradation RateResidue (mg/kg)Degradation RateResidue (mg/kg)Degradation Rate
12Mean(%)12Mean(%)12Mean(%)12Mean(%)
08.49.89.1-12.87.610.2-3.05.04.0-11.08.09.5-
76.78.17.418.99.67.48.516.61.03.02.049.86.07.06.531.5
142.22.02.177.34.17.15.645.32.01.01.562.55.05.05.047.3
212.62.12.474.01.04.02.575.22.03.02.537.81.03.02.078.8
280.51.00.891.62.01.01.585.12.02.02.050.33.02.02.573.7
Regression EquationY = 9.640−0.07386xY = 10.85−0.05616x-Y = 9.581−0.05508x
R20.88560.7615-0.8734
T1/2(d)9.412.3-12.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

He, M.; Wang, Y.; Zhang, L.; Mao, L.; Zhu, L.; Zheng, Y.; Liu, X.; Wu, C. Optimizing Analysis Methods: Rapid and Accurate Determination of Cuaminosulfate Residues with LC-MS/MS Based on Box–Behnken Design Study. Molecules 2024, 29, 794. https://doi.org/10.3390/molecules29040794

AMA Style

He M, Wang Y, Zhang L, Mao L, Zhu L, Zheng Y, Liu X, Wu C. Optimizing Analysis Methods: Rapid and Accurate Determination of Cuaminosulfate Residues with LC-MS/MS Based on Box–Behnken Design Study. Molecules. 2024; 29(4):794. https://doi.org/10.3390/molecules29040794

Chicago/Turabian Style

He, Mingyuan, Yuzhu Wang, Lan Zhang, Liangang Mao, Lizhen Zhu, Yongquan Zheng, Xingang Liu, and Chi Wu. 2024. "Optimizing Analysis Methods: Rapid and Accurate Determination of Cuaminosulfate Residues with LC-MS/MS Based on Box–Behnken Design Study" Molecules 29, no. 4: 794. https://doi.org/10.3390/molecules29040794

Article Metrics

Back to TopTop