Oxidation of amino acids by peracetic acid: Reaction kinetics, pathways and theoretical calculations.

Peracetic acid (PAA) is a sanitizer with increasing use in food, medical and water treatment industries. Amino acids are important components in targeted foods for PAA treatment and ubiquitous in natural waterbodies and wastewater effluents as the primary form of dissolved organic nitrogen. To better understand the possible reactions, this work investigated the reaction kinetics and transformation pathways of selected amino acids towards PAA. Experimental results demonstrated that most amino acids showed sluggish reactivity to PAA except cysteine (CYS), methionine (MET), and histidine (HIS). CYS showed the highest reactivity with a very rapid reaction rate. Reactions of MET and HIS with PAA followed second-order kinetics with rate constants of 4.6 ± 0.2, and 1.8 ± 0.1 M-1⋅s-1 at pH 7, respectively. The reactions were faster at pH 5 and 7 than at pH 9 due to PAA speciation. Low concentrations of H2O2 coexistent with PAA contributed little to the oxidation of amino acids. The primary oxidation products of amino acids with PAA were [O] addition compounds on the reactive sites at thiol, thioether and imidazole groups. Theoretical calculations were applied to predict the reactivity and regioselectivity of PAA electrophilic attacks on amino acids and improved mechanistic understanding. As an oxidative disinfectant, the reaction of PAA with organics to form byproducts is inevitable; however, this study shows that PAA exhibits lower and more selective reactivity towards biomolecules such as amino acids than other common disinfectants, causing less concern of toxic disinfection byproducts. This attribute may allow greater stability and more targeted actions of PAA in various applications.

Introduction

Peracetic acid (PAA), as a wide-spectrum antimicrobial agent, has been used in water disinfection and food sanitization since early 1940s (Rossoni and Gaylarde, 2000; Kitis, 2004). Commercial PAA solution is a colorless, quaternary equilibrium mixture of PAA, hydrogen peroxide (H2O2), acetic acid and water according to Equation (1) below:

The oxidation potential of PAA (1.76 V) is close to H2O2 (1.80 V), higher than aqueous chlorine (1.48 V) and chlorine dioxide (1.28 V) (Luukkonen and Pehkonen, 2017). PAA has been reported to be a more efficient disinfectant than chlorine and chlorine dioxide in inactivating pathogenic and indicator microbes (Salgot et al., 2002; Kitis, 2004; Dell'Erba et al., 2007). The advantages of PAA disinfection includes stability, low cost, and harmless disinfection byproducts (DBPs) (Dell'Erba et al., 2007; Voukkali and Zorpas, 2015). Thus, as a disinfectant or sterilant, PAA has a wide application in various industries including water treatment, food processing, medical, and textile (Kitis, 2004; Falsanisi et al., 2006). A number of pilot- or full-scale trials of using PAA as an alternative to chlorine for wastewater disinfection have been conducted in Europe and North America (Veschetti et al., 2003; Koivunen and Heinonen-Tanski, 2005; González et al., 2012; Eckert, 2013). In the U.S., the Environmental Protection Agency has authorized the use of PAA in fruits, vegetables and meat processing as well as in wastewater and storm water disinfection (EPA, 1998, 1999, 2012).

The traditional chlorine disinfection in water treatment and food production may cause the formation of various chlorinated DBPs, some of which are potentially carcinogenic (Fukayama et al., 1986). Since the increasing use of PAA, the possible reaction of organic matter with PAA and formed byproducts have drawn great attention. An earlier study reported that PAA failed to degrade several target pharmaceuticals in wastewater at low doses (Hey et al., 2012). On the other hand, Cai et al. (2017) demonstrated that the combination of PAA with UV could lead to rapid degradation of pharmaceuticals with investigation of involved reactive species. Zhang et al. (2017) showed that PAA alone could degrade several β-lactam antibiotics in surface water and wastewater relatively quickly; however, the reaction rates were much slower in comparison to those with common oxidants such as O3, chlorine, and Fe(VI). Meanwhile, the detected transformation products of β-lactams were mostly sulfoxide products instead of mineralized compounds, indicating the modest oxidizing ability of PAA. In general, the above studies were consistent with an early research that reported PAA solution was virtually unaffected in the presence of organic matter (Rodgers et al., 2004).

Amino acids are among the primary components of dissolved organic nitrogen in the aquatic environment. They are ubiquitous in the environment at concentrations up to 10 g⋅L−1 in surface waters (Trehy et al., 1986) and higher in wastewaters. Amino acids are also widely found in meats, fruits and vegetables in the form of proteins (Wu, 2009). Thus, amino acids have a great chance to encounter disinfectants during water treatment or food processing. The transformation of amino acids in chlorination processes has been well studied, with corresponding aldehydes and nitriles as the primary byproducts as well as toxic chlorinated products such as dihaloacetonitriles and chloral hydrate (Trehy et al., 1986; Pattison and Davies, 2001). In contrast, the potential reactions between amino acids and PAA during water treatment remain unclear and have been scarcely investigated. It is important to evaluate the transformation of amino acids and whether harmful byproducts are formed during PAA reactions in order to improve the applications of PAA in water treatment.

Some previous studies indicated that the amino acid composition in blood meal or organic matter may change after PAA treatment (Schnitzer and Hindle, 1980; Kerkaert et al., 2011; Hicks et al., 2015). Schnitzer and Hindle (1980) reported that PAA treatment may decrease the amino acid-N component in humic materials. Literature also suggested that oxidation by PAA may result in increase of hydrophilicity of the proteins through crosslinking destruction, whereas hypochlorous acid oxidation commonly results in protein aggregation (Kerkaert et al., 2011; Hicks et al., 2015). Extensive literature search indicates that there has been little information regarding the reaction kinetics and mechanism of the oxidation of amino acids by PAA to date. Thus, this study was motivated to address this knowledge gap.

The specific objectives of this work were to: (1) determine the reactivities of a range of amino acids towards PAA; (2) assess the impact of background H2O2 and solution pH on the oxidation reaction of amino acids with PAA; and (3) elucidate the transformation products and reaction mechanism via product characterization and theoretical analysis of amino acids' reactive sites. All the experiments were conducted using PAA dosage and pH conditions relevant to those common in wastewater treatment or food processing with PAA applications.

Materials and methods

Chemicals

Analytical-grade L-arginine (ARG), L-histidine (HIS), L-aspartic acid (ASP), and L-glycine (GLY) were purchased from Acros Organics. Other amino acids including L-cysteine (CYS), L-proline (PRO), L-methionine (MET), L-tyrosine (TYR), and L-glutamic acid (GLU) were obtained from MP Biomedicals at >98% purity. The structures of the selected amino acids are shown in Fig. S1. PAA solution (∼39% PAA, 6% H2O2 and ∼45% acetic acid, w/w), sodium dihydrogen phosphate, disodium hydrogen phosphate, sodium thiosulfate, and sodium borohydride were purchased from Sigma Aldrich (St. Louis, MO, USA). Hydrogen peroxide solution (∼30% H2O2 in water, w/w) was purchased from Fisher Scientific (Waltham, MA, USA). Reagent-grade deionized water (≥18 mΩ-cm) was generated from a Milli-Q integral water purification system (Millipore, Billerica, MA, USA) and used to prepare working solutions. The working solution of each amino acid at 10 mM was freshly prepared and used within 24 h. Fresh PAA working solutions at 100 mM were prepared every week by diluting the PAA stock. The H2O2 working solution was prepared at 11 mM. All stock and working solutions were stored at 4–5 °C before use.

PAA oxidation experiments

PAA oxidation experiments were conducted in 200-mL glass beakers with magnetic stirring at room temperature. Reaction solutions were maintained at desired pH using phosphate buffer (10 mM). Reactions were initiated by spiking an appropriate amount of PAA or H2O2 working solution into the buffered solution containing 10 μM amino acids. Aliquots were periodically taken and injected into HPLC vials containing an excess amount of sodium thiosulfate or sodium borohydride relative to the PAA concentration. The quenched samples were analyzed within 12 h after a derivatization process as described below. The reactions with H2O2 were conducted similarly by adding an appropriate amount of H2O2 working solution to the buffered reaction solution containing 10 μM selected amino acids.

Analytical methods

Titration methods were applied to determine the concentrations of PAA and H2O2, respectively, in the PAA stock solution as previously described (Cai et al., 2017). Iodometric titration was used to determine the combined concentration of PAA and H2O2, and potassium permanganate titration under acidic pH was used to determine the H2O2 concentration in the stock. The latter method was also used to determine the H2O2 concentration in pure H2O2 stock solution. The lower concentrations of PAA in the PAA working solution and in the oxidation experiments were quantified by the standard N,N-diethyl-p-phenylenediamine (DPD) titration method (APHA et al., 1998) as described previously (Cai et al., 2017).

The concentrations of residual amino acids in the reactions were determined using an Agilent 1100 high performance liquid chromatography (HPLC) equipped with a fluorescence detector after derivatization (detailed in Supporting Information (SI) Text S1). Monobromobimane (mBBr) was used as a derivatization reagent for CYS determination, while o-phthalaldehyde (OPA) was used for the other amino acids. After derivatization, 20 μL of sample was separated on an Agilent Zorbax SB-C18 column (2.1 × 150 mm, 5 μm) with a flow rate of 0.3 mL/min. The mobile phase was a mixture of (A) sodium acetate buffer (5 mM) at pH 5.9, and (B) acetonitrile. The eluent composition changed linearly from initially 10% A to 50% at 7 min and then was so maintained for 3 min, followed by re-equilibration to the starting condition. The excitation and emission wavelengths for fluorescence detection were set respectively at 395 and 475 nm for the mBBr derivatization method, and 340 and 455 nm for the OPA derivatization method.

The oxidation products of selected amino acids by PAA were analyzed using an Agilent 1100 Series LC/MSD system equipped with an Agilent Zorbax SB-C18 column (2.1 × 150 mm, 5 μm) with a flow rate at 0.3 mL/min. The mobile phase consisted of (A) deionized water with formic acid (0.1%, v/v), and (B) methanol. The gradient of the mobile phase B was kept at 5% in the first 5 min, then ramped to 90% at 10 min and maintained for 3 min, and then re-equilibrated to the starting condition. The injection volume was 10 μL. The mass spectrometer was set at negative electrospray ionization mode (ESI-) with a 70 V fragmentation voltage and 4000 V capillary voltage. Nitrogen was used both as the drying gas at 350 °C (6.0 L/min) and nebulizer gas (25 psi). A full-scan mode from 50 to 500 m/z was applied to examine the transformation products.

Computational method

All calculations were performed using Gaussian 03 software (Frisch et al., 2003). Geometry optimization calculations were executed without any constraints using the B3LYP method with 6-31G(d,p) basis set (Hariharan and Pople, 1973; Lee et al., 1988; Becke, 1993; Raghavachari, 2000). Condensed Fukui function (CFF), an important concept in conceptual density functional theory, was computed based on natural population analysis (NPA) charge to predict the regioselectivity of the oxidation of amino acids by PAA (Parr and Yang, 1984; Olah et al., 2002). All of the above calculations were performed in the gas phase.

Results and discussion

Evaluation of PAA reactivity to different amino acids

To assess the possible reaction between amino acids and PAA, a total of nine essential amino acids, which were widely detected in vegetables, fruits, and meats, were selected to study their reactivity towards PAA oxidation. As shown in Fig. 1, most amino acids are sluggish to PAA oxidation except CYS, MET and HIS. CYS was completely oxidized by PAA in 30 min, suggesting its high reactivity towards PAA oxidation. The high reactivity of CYS may be caused by its thiol group, which could be easily oxidized to its disulfide form (RS-SR) (Nakamoto and Bardwell, 2004; Rehder and Borges, 2010). MET and HIS were also observed with >25% transformation under the same condition. This is consistent with an early study that reported the sulfur-containing amino acids in meat would greatly decreased after PAA treatment, while most of the other amino acids showed little changes (Hicks et al., 2015). Unlike PAA treatment, chlorination was found with high reactivity towards all amino acids. The presence of amino acids during chlorination processes may greatly reduce the disinfection efficiency and produce DBPs such as trihalomethanes and halonitriles, which may have adverse effects on human health (Na and Olson, 2007). On the basis that the amounts of HIS, CYS, and MET in fruits and vegetables are typically low (Brückner and Westhauser, 2003) and the low reactivities of other amino acids to PAA, PAA oxidation probably will not change the organic composition of the treated food or wastewater significantly and much less byproducts would be formed, in comparison with other common oxidants such as chlorine.

Fig. 1

Losses of nine selected amino acids by PAA (Amino acids: 10 μM, PAA: 131 μM, incubation time: 30 min; room temperature). The order of nine amino acids in Fig. 1 is set based on their amounts in strawberry samples (from highest to lowest) (Gebhardt et al., 2008). Error bars represent standard deviation from duplicate samples.

Reaction kinetics of selected amino acids with PAA

Owing to their higher reactivity, CYS, MET and HIS were selected to further study their reaction kinetics with PAA. Notably, under the employed reaction conditions ([CYS]0 = 10 μM, [PAA]0 = 131 μM (i.e., 10 ppm), pH = 7.0), CYS was completely oxidized by PAA within 1 min. The reaction was too fast to quantify the rate constant using the batch reactor set-up. The reactions of MET and HIS with PAA were not as rapid and followed first-order kinetics with respect to the amino acid concentration (Fig. S2a). In addition, the pseudo-first-order rate constant (kobs, min−1) was found to increase linearly with PAA concentration (131–1310 μM) (Fig. S2b). Hence, the oxidation of amino acids by PAA can be expressed by the second-order kinetics (see Equations (2), (3)), with the apparent second-order rate constant (kapp, M−1⋅s−1) obtained by dividing the kobs value by the PAA concentration.where kobs is the observed first-order rate constant, kapp is the apparent second-order rate constant for the overall reaction, [Amino acid] and [PAA] are the concentration of amino acid and PAA, respectively.

The kapp values for MET and HIS with PAA were determined to be 4.6 ± 0.2, and 1.8 ± 0.1 M−1⋅s−1, respectively. To be noted, the kapp value for CYS towards PAA should be much larger, at least greater than 5.8 × 102 M−1⋅s−1, which is near the upper limit that could be measured by the batch reaction set-up in this study. The difference in kapp among amino acids towards PAA is likely caused by the side chain groups, which will be fully discussed later based on the products identification and reactive sites evaluation. Meanwhile, Table 1 indicates that the reaction rates of amino acids with other common oxidants are usually several orders of magnitude higher than those with PAA. This is similar to our earlier study that showed PAA has a slower reaction rate towards β-lactam antibiotics than other water treatment oxidants (Zhang et al., 2017). As a promising disinfectant, PAA shows limited ability to degrade organic biomolecules and micropollutants. This attribute may be beneficial for the further application of PAA in various industries, as PAA will not be easily consumed by organic matter in water and will not generate chlorinated DBPs or other harmful byproducts commonly formed from chlorine disinfection (Acero et al., 2010).

Table 1

Values for the second-order reaction rate constants (M−1⋅s−1) for CYS, MET, and HIS with chemical oxidants at pH 7.0–7.4.a

	HOCl	HBrO	O3	PAA
CYS	3.0 × 107	1.2 × 107	4.4 × 109	>5.8 × 102
MET	3.8 × 107	3.6 × 106	4.0 × 106	4.6
HIS	1.0 × 105	3.0 × 106	5.3 × 106	1.8

Rate constants were from experiments in this study for PAA (pH 7.0) and from literature (Pryor et al., 1984; Pattison and Davies, 2001, 2004) for HClO, HBrO, and O3 (pH 7.4).

Effect of H2O2 on oxidation of amino acids by PAA

H2O2 is an important component in commercial PAA solutions; however, the role of H2O2 in the oxidation of amino acids by PAA solutions is unknown. In this work, CYS, MET and HIS were reacted with H2O2 only to study the impact of H2O2 on amino acids oxidation. It was found that the concentration of background H2O2 (about 33 μM) in the PAA solution was too low to cause any significant transformation of amino acids at pH 7. Thus, the concentration of H2O2 was increased to 131 μM (the same as the PAA concentration, as some commercial PAA solutions may have a lower PAA/H2O2 molar ratio (Shah et al., 2015)) in order to assess the reaction kinetics. Fig. 2 shows that 131 μM of H2O2 could only oxidize 15% of CYS, 5.5% of MET, and 3.4% of HIS in 30 min. Meanwhile, under the same concentration of PAA, PAA could oxidize 100%, 65%, and 28% of the amino acids, respectively. The second-order rate constants of CYS and MET with H2O2 at pH 7 were about 5.7 × 10−1 and 8.9 × 10−2 M−1⋅s−1, respectively. These kinetic results are comparable to the earlier research that have calculated the second-order rate constants of CYS and MET in the oxidation by H2O2 to be around 2.7 × 10−1 M−1⋅s−1 (pH 5.7) and 1.6 × 10−2 (pH 5.5–6.0), respectively (Boonvisut et al., 1982; Chu et al., 2016). Evidently, the oxidation of amino acids by H2O2 is much slower than that by PAA, indicating that PAA is the primary oxidant in PAA solution dominating the transformation of amino acids. This finding is consistent with our previous study that it was PAA, not H2O2, that induced the oxidation reaction of β-lactams, and there was no obvious synergistic or inhibitory effect between PAA and H2O2 (Zhang et al., 2017). Earlier literature also has concluded that the antimicrobial effect of PAA is far greater than H2O2 (Lubello et al., 2002; Kitis, 2004). Thus, in the applications of PAA solutions, it is mainly PAA that is responsible for antibacterial and compound oxidation.

Fig. 2

Time course of transformation of CYS, MET, and HIS with PAA solution (solid dots) and H2O2 solution (valid dots). Experimental conditions: [amino acids]0 = 10 μM, [PAA]0 = 131 μM, [H2O2]0 = 131 μM, pH = 7, room temperature.

Effect of pH on PAA oxidation of amino acids

PAA and amino acids are both ionizable compounds. The primary species of PAA and amino acids under different pH conditions may vary according to the acid dissociation constants, which may cause different reaction behaviors in the oxidation of amino acids by PAA. As shown in Fig. 3, the pseudo first-order rate constants of MET and HIS by PAA oxidation at pH 5 and 7 were comparable, but decreased significantly at pH 9.

Fig. 3

Pseudo first-order rate constants of MET, and HIS with PAA at different pHs. [amino acids]0 = 10 μM, [PAA]0 = 131 μM, room temperature.

The pKa of PAA is around 8.2 (Luukkonen and Pehkonen, 2017) and the acid-base dissociation could result in varied distribution of neutral (PAAH) and charged (PAA−) species. According to the speciation, neutral PAAH is dominant (over 94%) at pH 5 and 7, while anionic PAA− becomes dominant (over 89%) at pH 9. Speciation of amino acids may be another factor. For MET and HIS, the acidity constants of the carboxylic acid group on the α-carbon are 2.28 and 1.82, and the acidity constants of the protonated α-amino group are 9.21, and 9.17, respectively. HIS has an additional acidity constant for its imidazole ring —NH at 6.00 (Linde, 1991). The species distributions of PAA, MET and HIS at pH 5, 7, and 9 are depicted in Table S1 and Fig. S3.

From pH 5 to 7, the amount (99.9–92.1%) of PAAH changed little and could explain why the rate constants were comparable at pH 5 and 7 as it was PAAH having the strong oxidizing ability rather than PAA− (Zhang et al., 2017). At pH 9, the lower distribution (10.4%) of PAAH could lead to the lower rate constants for oxidation of amino acids. For amino acids' speciation, MET had the same neutral, zwitterion species dominant at pH 5 and 7, while HIS's imidazole —NH group changed from being protonated to non-protonated from pH 5 to 7 (Fig. S3). The very similar rate constants of HIS at pH 5 and 7 suggest that the impact of HIS species on HIS oxidation by PAAH was negligible. Thus, it can be concluded that the rate constants for PAA oxidation are lower at alkaline solution due to the lower distribution of PAAH species, while the speciation of the targeted amino acids have much less influence on the reaction rate.

Reaction products and pathways

The transformation products of CYS, MET, and HIS by PAA oxidation were assessed and the results are presented in Table 2, with proposed transformation pathways in Fig. 4. The MS results showed that disulfide (CYS dimer), sulfinic acid (CYSO2H), and sulfonic acid (CYSO3H) were the three major products of oxidation of CYS by PAA. The formation of CYSOH was not detected, which was likely due to its high instability (Penn et al., 1978). The above compounds are well-known CYS oxidation products (2M-1, M+16, M+32, and M+48, where M = molecular weight of parent CYS) under various conditions because of the reactive thiol functional group (Milne and Zika, 1993; Greer et al., 1996; Hawkins and Davies, 2001; Mailloux et al., 2014). In this study, PAA likely attacks the highly nucleophilic thiol group, which may result in formation of CYSOH. Compounds with the —SOH group were found to be extremely difficult to detect and quantify due to the much lower pKa than that of the —SH group (Enami et al., 2009; Hugo et al., 2009; McGrath et al., 2010). CYSOH may couple with another CYS to form the CYS dimer through intermolecular disulfide bridging. CYSOH can also undergo further oxidation to generate CYSO2H and CYSO3H (Fig. 4a).

Table 2

MS identified products of CYS, MET, and HIS with PAA oxidation.a

	m/z	Molecular weight	Retention time (min)	Peak Areab (×105)
CYS	120	121(M)	1.60	–
	152	M+32	1.69	0.06
	168	M+48	1.56	0.56
	239	2M −1	1.58	0.15
MET	148	149 (M)	2.31	1.66
	164	M+16	1.69	0.71
HIS	154	155 (M)	2.20	4.00
	186	M+32	2.12	4.76

Initial concentration of amino acids and PAA were 10 and 131 μM, respectively.

Peak area values for amino acids and corresponding products at 1 min for CYS, 30 min for MET and HIS.

Fig. 4

Proposed reaction pathways for (a) CYS, (b) MET, and (c) HIS oxidation by PAA.

Ion chromatography was employed to determine whether oxidation of CYS by PAA could generate small ions. The results suggested that CYSO3H was the final product as only a new peak of CYSO3 – emerged at 12.9 min, while no SO4 2– (standard R.T. = 11.2 min) or other small organic carboxylic acids were detected (Fig. S4). This is rather different from other chemical oxidants, such as O3 or HClO, which were able to break bonds and finally form ions including aliphatic acids, SO42−, and/or NO3− (Mudd et al., 1969; Wenk et al., 2013).

The reaction pathway for MET oxidation by PAA is rather simple that MET was expected to be oxidized to generate the M+16 intermediate (i.e. MET sulfoxide) (Fig. 4b). Previous research also observed MET sulfoxide as the primary product during oxidation of MET by H2O2 (Keck, 1996). MET sulfoxide is believed to contribute to biological aging as it increases with age in body tissues. Oxidation of MET may cause the accumulation of damaged, nonfunctional proteins (Shringarpure and Davies, 2002; Stadtman et al., 2005).

The M+32 (m/z 186) is the predominant intermediate (based on peak area) of HIS oxidation by PAA, while the dimeric product (m/z 323) may be formed via nucleophilic addition of m/z 186 by another HIS at the C9 site followed by loss of a water (Fig. 4c). Similar reaction products were also detected in the oxidation of HIS by 1O2, •OH, or in other oxidation systems (Tomita et al., 1969; Agon et al., 2006; Amano et al., 2014; Liu et al., 2014).

Based on the results above, it can be concluded that PAA is a weak oxidant that the corresponding oxidation products are mostly substrates with one or more [O] addition, which is consistent with previous observations (Zhang et al., 2017). The thiol group of CYS was oxidized to sulfinic and sulfonic acids, as well as to generate CYS dimers. The thioether sulfur of MET is oxidized to sulfoxide by PAA. [O] addition on the imidazole ring was also demonstrated to be the primary reaction pathway for HIS oxidation by PAA. Moreover, PAA appears to be only reactive to organics with specific functional groups, such as thiol, thioether sulfur, as well as specific groups with high electron distribution.

Theoretical prediction of reactive sites on selected amino acids

PAA is a well-known electrophilic oxidant, as the molecular orbitals shown in Fig. S5 well demonstrates the electron-deficient property around the peroxyl group of PAA. According to the products identified and general oxidation reaction mechanism, PAA is likely to accept electrons from amino acids and give out an oxygen atom. The condensed Fukui function (CFF), commonly used in prediction of reactive sites of electrophilic, nucleophilic, and radical attacks (Parr and Yang, 1984), was used to predict the reactive sites on selected amino acids towards PAA attacks. The detailed description for the CFF concept and the calculation method are presented in Text S2.

As presented in Fig. 5, the ability of losing electron was quantified and visualized on every site of amino acids by the calculated Fukui index (). The on sulfur atom in CYS and MET is much higher than that on the other sites, suggesting the strong affinity of sulfur-containing organics towards PAA oxidation, which is in accordance with the fact that most of the known organic compounds that react fast with PAA are those with sulfur atom(s) (Chipiso et al., 2016; Zhang et al., 2017). Moreover, the corresponding products are mostly sulfoxides and other following products. Meanwhile, the distribution on HIS molecule is less differential, with values of 0.103, 0.095, 0.116, and 0.148 at C1, C5, C7, and C9, respectively. Thus, the oxygen atom from PAA likely prefers addition onto C9 or C7 of HIS as we proposed in Fig. 4c. Then, the [O] addition intermediate coupled with another HIS molecule and formed a dimer product. Based on the products identified and theoretical calculations, we proposed that the sites with higher CFF values are more likely to be attacked by PAA, forming corresponding [O] addition products and other derivatives.

Fig. 5

Visualized Fukui index () on selected amino acids by electrophilic PAA attacks. (a) CYS, (b) MET, (c) HIS, (d) Fukui electrophilic attack index () distribution data. The diameters of blue circles represent relative values on the specific site.

Since the PAA-induced direct oxidation was considered to be electrophilic attacks, the reactivity of amino acids towards PAA oxidation is likely related to their nucleophilicities. Thus for the comparison of different amino acids, the nucleophilicity index was used, which measures the energy stabilization when an optimal electronic charge transfer from the system to the environment occurs (De Vleeschouwer et al., 2007). To describe the nucleophilic character of a reactive site within a molecule, a local nucleophilicity index can be proposed:where N is the global nucleophilic index, which is calculated based on the electronic chemical potential and chemical hardness (see SI Text S3); and displays its maximum value, i.e. at the active site (k) of the nucleophile.

Note that the local nucleophilicity index contains a global contribution N as a factor of the genuine local nucleophilic index . Table 3 summarizes the values of electrostatic descriptors, including the global and local nucleophilicity indexes, for 12 amino acids. CYS shows the S5 sulfur atom as the preferential site for electrophilic attack on the basis of the Fukui function,  = 0.439, and the local nucleophilicity,  = 0.850 eV. MET presents a local nucleophilicity value of  = 0.768 eV at sulfur atom S6. Note that HIS, having multiple reactive sites due to the close Fukui function values at C7 and C9, presents a local nucleophilicity value of  = 0.653 eV with the sum of two reactive sites. The local nucleophilicity index of CYS, MET and HIS are much higher than those of the other amino acids, suggesting that they are easier to undergo electrophilic attacks on the specific sites. This is consistent with our experimental results that CYS, MET and HIS were oxidized more easily than the other amino acids by PAA (Fig. 1). Moreover, the Nk values of three more amino acids (leucine (LEU), threonine (THR), and tryptophan (TRP)) with typical alkyl, hydroxyl and indole side chains were also calculated, and their reactivities towards PAA oxidation can thus be predicted. The results in Table 3 indicate that LEU and THR may show less reactivity, while TRP could be a potential substrate for PAA oxidation due to its side chain indole. According to our work and published literature, compounds containing sulfur atoms or rings with nitrogen atoms are always found to be highly nucleophilic, and more likely to undergo PAA oxidation (Chipiso et al., 2016; Zhang et al., 2017). However, descriptor could only qualitatively predict the reactivities of amino acids, and more parameters and factors should be considered to establish a better model.

Table 3

Electrostatic descriptors for 12 amino acids at B3LYP/6-31g(d,p) level. The units for these descriptors are all eV except .

Amino acids	Electron affinity (A)	Ionization potential (I)	Electronic chemical potential (μ)	Chemical hardness (η)	Global nucleophilicity index (N)	fk−	Local nucleophilicity index (N_k)
ARG	−1.338	8.048	−3.355	9.386	1.668	0.271	0.452
ASP	−2.208	9.228	−3.510	11.437	1.856	0.110	0.204
CYS	−2.130	8.874	−3.372	11.004	1.936	0.439	0.850
GLU	−1.880	9.114	−3.617	10.995	1.681	0.105	0.176
GLY	−2.846	9.478	−3.316	12.324	2.242	0.161	0.361
HIS	−2.330	8.060	−2.865	10.390	2.532	0.258a	0.653
MET	−1.417	8.089	−3.336	9.506	1.708	0.450	0.768
PRO	−2.717	8.460	−2.872	11.177	2.710	0.125	0.339
TYR	−1.636	7.855	−3.110	9.491	1.963	0.152	0.298
LEU	−2.728	9.282	−3.277	12.009	2.237	0.121	0.271
THR	−2.639	9.341	−3.351	11.980	2.133	0.118	0.252
TRP	−1.269	7.577	−3.154	8.846	1.778	0.285a	0.507

The sum of the largest two values was applied due to the wide distribution of reactive sites on HIS molecule.

Conclusions

This work demonstrates that PAA has low reactivity with amino acids except those with reactive side chains including thiol group (CYS), thioether group (MET), and imidazole ring (HIS). The background H2O2 in the PAA solution has a negligible effect on the reactions of amino acids with PAA. PAA oxidation of MET and HIS is much slower at pH 9 than at pH 5 and 7 due to the speciation of PAA. The combined results from experimental measurements and theoretical calculations show that PAA is prone to attack amino acids with strong nucleophilicity, delivering [O] to the electron-rich sites. Overall, it can be concluded that the widely used PAA solution is much less reactive towards organics than other common disinfectants, thus resulting in longer stability and fewer byproducts formation. Moreover, PAA treatment primarily generates [O] addition products which are less harmful compared to other DBPs of concerns commonly found with disinfection oxidants such as chlorine. From the perspective of disinfectant stability and DBPs formation, the results of this study suggest that PAA is a more stable and safer oxidant, which could be suitable in food processing, wastewater treatment as well as other industries in the future.

Conflict of interest

The authors declare no conflict of interest.