Kinetics, Thermodynamics and Mechanism of Enzymatic Degradation of Zearalenone in Degummed Corn Oil

The kinetics and thermodynamics of the enzymatic degradation of zearalenone (ZEN) in degummed corn oil were investigated by analyzing the impacts of temperature, pH, ZEN hydrolase dosage and ZEN concentration on the initial reaction rate. The kinetic study found that the maximum reaction rate was 0.97 μmol × kg−1 min−1, the Michaelis constant (Km) was 11,476 μmol × kg−1 and the Michaelis equation was V = 0.97[S]/(11,476 + [S]). The thermodynamic study showed that the activation energy (Ea) was 70.37 kJ·mol−1, the activation enthalpy change of the reaction (ΔH) > 0, the free energy of activation (ΔG) > 0 and the activation entropy change (ΔS) < 0, indicating the reaction could not be spontaneous. The reaction mechanism of ZEN was studied by a hybrid quadrupole orbitrap mass spectrometer. It was found that ZEN first generated the intermediate G/L/D/W-ZEN+H2O, followed by generating the intermediate W-ZEN-H2O under the action of a degrading enzyme. Then, the lactone bond was opened to produce C18H24O6, and finally the decarboxylation product C17H24O4 formed automatically.


Introduction
Zearalenone (ZEN), a non-steroidal estrogen mycotoxin, is mainly produced by Fusarium roseum, Fusarium graminearum, Fusarium tricinum and other Fusarium microorganisms during their secondary metabolism and is the most widely contaminated mycotoxin in the world [1,2]. Corn is easily polluted by ZEN in the process of planting and harvesting [3]. As ZEN is a benzoic acid lactone compound and exhibits a certain lipophilicity, it is very easily transferred into crude corn oil during corn oil production. The highest content of ZEN in corn crude oil was more than 8000 µg × kg −1 [4]. Zearalenone has a similar phenolic dihydroxyacetone structure to β-Estrogens, so it can competitively bind to estrogen receptors in the body after entering humans and animals, thus causing reproductive toxicity, genotoxicity, carcinogenicity and immunotoxicity [5][6][7].
Rapid and effective detoxification methods have attracted increasing attention due to the huge economic losses and health risks caused by ZEN pollution. At present, the detoxification methods of ZEN can be divided into three categories: physical, chemical and biological methods. However, the physical and chemical methods have defects such as incomplete detoxification, the loss of some small molecular nutrients, secondary pollutants and low detoxification efficiency [8][9][10][11][12]. Biological enzymatic detoxification has become a research hotspot because of its irreversible reaction process, mild reaction conditions, no secondary pollutants, low price and many other advantages [12]. Zearalenone hydrolase (ZHD) has the advantages of a clear detoxification mechanism, high degradation activity and no medium required for reaction, showing good application prospects [13].
In 2002, Takahashi Ando et al. carried out a thorough study on the mechanism of ZEN degradation catalyzed by lactone hydrolase ZHD101 (as shown in Figure 1) and found that lactone hydrolase can catalyze the cleavage of lactone bonds in ZEN molecules, opening the lactone ring, producing a molecule of carboxyl and a molecule of hydroxyl, and then spontaneously decarboxylating to form a completely non-toxic degradation product (1-(3,5-dihydroxybenzene-10 -hydroxy-1 -trans undecene-6 -one)). A molecule of carbon dioxide is released at the same time [14]. The catalytic mechanism of the enzyme is mainly the result of the joint action of amino acid residues. During the catalytic process, intermediate products bound with amino acids will be formed [21]. degradation activity and no medium required for reaction, showing good application prospects [13]. In 2002, Takahashi Ando et al. successfully expressed the lactone hydrolase encoding gene ZHD101 from Clonostachys rosea IFO7063 in Escherichia coli [14]. Since then, ZHD101 has gradually become a research hotspot and has been successfully expressed in a variety of hosts including Escherichia coli [14], Pichia pastoris [15], Saccharomyces cerevisiae [16] and Lactobacillus reuteri [17]. To obtain ZEN lactone hydrolase with better enzymology properties, some other microbial lactone hydrolases, such as ZHD607 [15] from Phylophora americana and Neurospora crassa ZENC [18], zhd518 [19] from Rhinocladiella machinenziei, CbZHD [20] from Cladophialophora bantiana and ZHD607 [15] from Phylophora americana, have been gradually excavated and studied.
In 2002, Takahashi Ando et al. carried out a thorough study on the mechanism of ZEN degradation catalyzed by lactone hydrolase ZHD101 (as shown in Figure 1) and found that lactone hydrolase can catalyze the cleavage of lactone bonds in ZEN molecules, opening the lactone ring, producing a molecule of carboxyl and a molecule of hydroxyl, and then spontaneously decarboxylating to form a completely non-toxic degradation product (1-(3,5-dihydroxybenzene-10′-hydroxy-1′-trans undecene-6′-one)). A molecule of carbon dioxide is released at the same time [14]. The catalytic mechanism of the enzyme is mainly the result of the joint action of amino acid residues. During the catalytic process, intermediate products bound with amino acids will be formed [21]. It has been reported that ZEN removal from corn oil by enzyme is an environmentally friendly and efficient method [22]. However, the kinetics and mechanism of reaction in crude corn oil have not been clarified. In this study, the kinetics, thermodynamics and mechanism of the enzymatic degradation of zearalenone in degummed corn oil (DCO) were investigated. The structure of the degradation products was identified to explore its reaction mechanism, so as to provide a theoretical basis for the removal of ZEN from corn oil by ZHD.

Effects of pH on Initial Reaction Rate
As shown in Figure 2, the effects of different pH on the initial reaction rate were investigated at 35 °C and the ZHD dosage was 30 mg × kg −1 . It has been reported that ZEN removal from corn oil by enzyme is an environmentally friendly and efficient method [22]. However, the kinetics and mechanism of reaction in crude corn oil have not been clarified. In this study, the kinetics, thermodynamics and mechanism of the enzymatic degradation of zearalenone in degummed corn oil (DCO) were investigated. The structure of the degradation products was identified to explore its reaction mechanism, so as to provide a theoretical basis for the removal of ZEN from corn oil by ZHD.

Effects of pH on Initial Reaction Rate
As shown in Figure 2, the effects of different pH on the initial reaction rate were investigated at 35 • C and the ZHD dosage was 30 mg × kg −1 .
It can be seen from Figure 2 that the initial reaction rate changes in a bell shape with the increase in pH. When the pH was in the range of 7.5-8.5, the initial reaction rate was the highest, which is related to the dissociation state of the active group of the enzyme in different acid-base systems. The combination of the active group of the enzyme and ZEN is in the best state when the pH is 7.5-8.5, making the reaction rate faster. When the pH is higher or lower than the value, the action of the acid or alkali causes the dissociation of the active group of the enzyme to change into an inactive state, thus reducing the reaction rate [23].  It can be seen from Figure 2 that the initial reaction rate changes in a bell shape with the increase in pH. When the pH was in the range of 7.5-8.5, the initial reaction rate was the highest, which is related to the dissociation state of the active group of the enzyme in different acid-base systems. The combination of the active group of the enzyme and ZEN is in the best state when the pH is 7.5-8.5, making the reaction rate faster. When the pH is higher or lower than the value, the action of the acid or alkali causes the dissociation of the active group of the enzyme to change into an inactive state, thus reducing the reaction rate [23].

Effects of ZHD Dosage on Initial Reaction Rate
As shown in Figure 3, the effect of different ZHD dosages on the initial reaction rate was investigated at 35 °C and pH 8.0.

Effects of ZHD Dosage on Initial Reaction Rate
As shown in Figure 3, the effect of different ZHD dosages on the initial reaction rate was investigated at 35 • C and pH 8.0. It can be seen from Figure 2 that the initial reaction rate changes in a bell shape with the increase in pH. When the pH was in the range of 7.5-8.5, the initial reaction rate was the highest, which is related to the dissociation state of the active group of the enzyme in different acid-base systems. The combination of the active group of the enzyme and ZEN is in the best state when the pH is 7.5-8.5, making the reaction rate faster. When the pH is higher or lower than the value, the action of the acid or alkali causes the dissociation of the active group of the enzyme to change into an inactive state, thus reducing the reaction rate [23].

Effects of ZHD Dosage on Initial Reaction Rate
As shown in Figure 3, the effect of different ZHD dosages on the initial reaction rate was investigated at 35 °C and pH 8.0.  It can be seen from Figure 3 that the initial reaction rate increases linearly with the increase in the ZHD dosage and then tends to slow down. During enzyme catalysis, the enzyme firstly combines with ZEN to form an intermediate complex. When the enzyme concentration is low, the reaction rate is slow. With the increase in the enzyme concentration, the rate increases in proportion. When the enzyme concentration is large, the enzyme changes from fully saturated to excessive. Under the condition of constant substrate concentration, the reaction rate does not change much [24].

Effects of ZEN Concentration on Initial Reaction Rate
As shown in Figure 4, the effect of different ZEN concentrations on the initial reaction rate was investigated at 35 • C and pH 8.0 and the ZHD dosage was 30 mg × kg −1 . changes from fully saturated to excessive. Under the condition of constant substrate concentration, the reaction rate does not change much [24].

Effects of ZEN Concentration on Initial Reaction Rate
As shown in Figure 4, the effect of different ZEN concentrations on the initial reaction rate was investigated at 35 °C and pH 8.0 and the ZHD dosage was 30 mg × kg −1 . Zearalenone degradation by ZHD is a hydrolysis reaction, and the reaction substrate is ZEN and water [25]. However, since the concentration of water changes little during the reaction process, it can be treated as a constant, so the reaction can be regarded as a single substrate enzyme catalyzed reaction. It can be seen from Figure 4 that when the enzyme concentration is constant, the initial reaction rate catalyzed by the enzyme increases in a positive proportion with the increase in the substrate concentration, showing a first-order reaction.

Enzyme Reaction Kinetics
The Lineweaver Burk method is used for drawing Figure 4, and the results are shown in Figure 5. Zearalenone degradation by ZHD is a hydrolysis reaction, and the reaction substrate is ZEN and water [25]. However, since the concentration of water changes little during the reaction process, it can be treated as a constant, so the reaction can be regarded as a single substrate enzyme catalyzed reaction. It can be seen from Figure 4 that when the enzyme concentration is constant, the initial reaction rate catalyzed by the enzyme increases in a positive proportion with the increase in the substrate concentration, showing a first-order reaction.

Enzyme Reaction Kinetics
The Lineweaver Burk method is used for drawing Figure 4, and the results are shown in Figure 5. According to the Lineweaver Burk method, the intersection point with the Y axis is Vm −1 , and the intersection point with the X axis is −Km −1 . Fitting the obtained curve, the equation (y=11,626x + 0.9355) is obtained, so that Vm is 1.07 μmol × kg −1 × min −1 and Km is 12,428 μmol × kg −1 . Thus, the Michaelis equation of this reaction is V = 1.07[S]/(12,428 + [S]).
The Wilkinson statistical method includes two steps of nonlinear multiplication and Taylor expansion as shown in Tables 1 and 2. Figure 5 is solved by the Wilkinson statistical method. As can be seen from Table 1 Vm0 and Km0 are the estimated values of the maximum reaction rate and Michaelis constant, respectively. According to the Lineweaver Burk method, the intersection point with the Y axis is Vm −1 , and the intersection point with the X axis is −Km −1 . Fitting the obtained curve, the equation (y = 11,626x + 0.9355) is obtained, so that Vm is 1.07 µmol × kg −1 × min −1 and Km is 12,428 µmol × kg −1 . Thus, the Michaelis equation of this reaction is V = 1.07[S]/(12,428 + [S]).
The Wilkinson statistical method includes two steps of nonlinear multiplication and Taylor expansion as shown in Tables 1 and 2. Figure 5 is solved by the Wilkinson statistical method. As can be seen from Table 1 Vm 0 and Km 0 are the estimated values of the maximum reaction rate and Michaelis constant, respectively.
The f and f' in Table 2 are calculated according to the following formula: According to Table 2, the following results can be obtained that The final Michaelis equation calculated by the above statistical method is There is a certain difference between Vm and Km calculated by the Lineweaver buck method and the Wilkinson statistical method. The Wilkinson statistical method is considered to be reliable [26]. Therefore, the Michaelis equation obtained by the Wilkinson statistical method is finally selected as the final result.
There is a certain gap between Km and V max obtained in this study and the data reported by other scholars [18]. Km reflects the affinity between the enzyme and substrate. The larger the Km, the smaller the affinity between the enzyme and the substrate, and the higher the required substrate concentration when the reaction rate reaches 1/2 of the maximum reaction rate [27]. The substrate ZEN in this study is dissolved in the oil phase, while the enzyme is dissolved in the water phase. The two phases' immiscibility requires a large substrate concentration to increase the contact probability between the enzyme and the substrate. Therefore, the Km obtained in this study is larger than in others.

Effects of Temperature on Initial Reaction Rate
As shown in Figure 6, the effect of different temperatures on the initial reaction rate was investigated at Ph 7.5, ZHD dosage 30 mg × kg −1 . a large substrate concentration to increase the contact probability between the enzyme and the substrate. Therefore, the Km obtained in this study is larger than in others.

Effects of Temperature on Initial Reaction Rate
As shown in Figure 6, the effect of different temperatures on the initial reaction rate was investigated at Ph 7.5, ZHD dosage 30 mg × kg −1 . It can be seen from Figure 6 that the initial reaction rate increases first and then decreases with the increase in temperature. This is because the influence of temperature on the enzyme reaction rate includes two aspects. On the one hand, the material thermodynamic movement intensifies with increasing temperatures, and the enzyme reaction rate increases because the contact probability between substrate molecules increases. On the other hand, the enzyme protein denaturates gradually with the increase in temperature, and the enzyme activity decreases, resulting in a decrease in the reaction rate [28]. The temperature when the rate reaches the highest point is the optimal reaction temperature of the enzyme. It can be seen from Figure 2 that the enzyme reached the optimum reaction It can be seen from Figure 6 that the initial reaction rate increases first and then decreases with the increase in temperature. This is because the influence of temperature on the enzyme reaction rate includes two aspects. On the one hand, the material thermodynamic movement intensifies with increasing temperatures, and the enzyme reaction rate increases because the contact probability between substrate molecules increases. On the other hand, the enzyme protein denaturates gradually with the increase in temperature, and the enzyme activity decreases, resulting in a decrease in the reaction rate [28]. The temperature when the rate reaches the highest point is the optimal reaction temperature of the enzyme. It can be seen from Figure 2 that the enzyme reached the optimum reaction temperature at about 40 • C. The result was consistent with the conclusion reported by Wang [19].

Determination of Thermodynamic Parameters
It can be seen from Figure 6 that, in the range of 25-40 • C, the initial reaction rate is linear with the increase in temperature. As shown in Figure 7, the Arrhenius plot was applied to calculate the activation energy (Ea). The activation energy (Ea) is 70.37 kJ × mol −1 .

Determination of Thermodynamic Parameters
It can be seen from Figure 6 that, in the range of 25-40 °C, the initial reaction rate is linear with the increase in temperature. As shown in Figure 7, the Arrhenius plot was applied to calculate the activation energy (Ea). The activation energy (Ea) is 70.37 kJ × mol −1 .
At present, there is no data report on the Ea of ZHD, but the Ea of lactone hydrolase for other use is 56.8kJ mol −1 , which is little lower than that calculated in this study [29]. The slightly higher activation energy (Ea) indicates that the formation of the enzyme substrate complex requires more energy. This feature makes the industrial application of ZHD require more energy supplement, and this feature should be taken into account in the process application. As shown in Table 3, the thermodynamic parameters for ZEN degradation by ZHD were calculated according to Equations (4)-(6).  At present, there is no data report on the Ea of ZHD, but the Ea of lactone hydrolase for other use is 56.8 kJ mol −1 , which is little lower than that calculated in this study [29]. The slightly higher activation energy (Ea) indicates that the formation of the enzyme substrate complex requires more energy. This feature makes the industrial application of ZHD require more energy supplement, and this feature should be taken into account in the process application.
As shown in Table 3, the thermodynamic parameters for ZEN degradation by ZHD were calculated according to Equations (4)-(6). Table 3. Thermodynamic parameters for ZEN degradation by ZHD. It can be seen from Table 3 that ∆H and ∆G are more than 0, and ∆S is less than 0, indicating that the reaction cannot be spontaneous. The enthalpy change of activation usually represents the energy transferred from the outside to make the active site of the enzyme fit with the substrate. If the value is greater than 0, it indicates that the reaction is endothermic [30]. ∆S < 0 indicates that the reaction is a process of entropy reduction [31], which is because ZEN is dissolved in oil and needs to be exchanged to the reaction interface, thus losing some degrees of freedom. ∆G > 0 indicates that the reaction is not spontaneous [32], and it needs to increase the oil-water contact area by providing energy.

Mechanism of ZEN Hydrolase
As shown in Figure 8, the ZEN in the DCO before the reaction was analyzed by LC-MS. It can be seen from Figure 8  fragment ion (m/z 301.14) under mass spectrometry conditions. Francesca et al. [33] reported that the above two ions are characteristic ions of zearalenone LC-MS, and the results are consistent with the literature.
The chromatogram of the negative ion scanning results of the enzymolysis reaction products is shown in Figure 9, wherein Figure 9A shows the chromatographic peak of ZEN (5.43 min); Figure 9B,C show that ZEN in DCO generates product B (2.18 min) and product C (5.67 min), respectively, under the condition of enzymatic hydrolysis.
The mass spectra of B and C in Figure 9 are shown in Figures 10A and 10B It can be seen from Figure 10 that the molecular weights of the two degradation products of ZEN are 335.15 and 291.16, respectively. The molecular weight of A is 18 higher than that of ZEN, which is the molecular weight of water. This indicates that the hydrolysis of the lactone bond of ZEN opens and becomes carboxyl and hydroxyl groups. The molecular weight of B (43.99) was less than that of A, which is the molecular weight of carbon dioxide, indicating that product A (water-added product) further reacts to remove the carboxyl group generated by hydrolysis to obtain product B (decarboxylation product).
The chromatogram of the positive ion scanning results of the enzymolysis reaction products is shown in Figure 11. It can be seen from Figure 8 that the chromatographic peak with a retention time of 8.52 min represents ZEN. The chromatographic identification process with 6.30 min and 8.13 min retention times is shown in Figures 12 and 13.

Mechanism of ZEN Hydrolase
As shown in Figure 8, the ZEN in the DCO before the reaction was analyzed by LC-MS. It can be seen from Figure 8 that ZEN has a peak retention time of 8.51 min in ESI+ ion mode, and can form the [C18H22O5+H] + parent ion (m/z 319.15) and [C18H22O5-H2O] + fragment ion (m/z 301.14) under mass spectrometry conditions. Francesca et al. [33] reported that the above two ions are characteristic ions of zearalenone LC-MS, and the results are consistent with the literature.
The chromatogram of the negative ion scanning results of the enzymolysis reaction products is shown in Figure 9, wherein Figure 9A shows the chromatographic peak of ZEN (5.43 min); Figure 9B,C show that ZEN in DCO generates product B (2.18 min) and product C (5.67 min), respectively, under the condition of enzymatic hydrolysis. The mass spectra of B and C in Figure 9 are shown in Figure 10A and Figure 10B It can be seen from Figure 10 that the molecular weights of the two degradation prod-  The chromatogram of the positive ion scanning results of the enzymolysis reaction products is shown in Figure 11. It can be seen from Figure 8 that the chromatographic peak with a retention time of 8.52 min represents ZEN. The chromatographic identification process with 6.30 min and 8.13 min retention times is shown in Figures 12 and 13. As shown in Figure 12, the mass spectrum with a retention time of 6.30 min implies that the parent ion of the substance is double-charged m/z 808.61. Further, the difference The chromatogram of the positive ion scanning results of the enzymolysis reaction products is shown in Figure 11. It can be seen from Figure 8 that the chromatographic peak with a retention time of 8.52 min represents ZEN. The chromatographic identification process with 6.30 min and 8.13 min retention times is shown in Figures 12 and 13. As shown in Figure 12, the mass spectrum with a retention time of 6.30 min implies that the parent ion of the substance is double-charged m/z 808.61. Further, the difference is calculated by the fragment ions (m/z 752.06, m/z 695.52, m/z 622.97) to get exactly three charge (m/z = 809, z = 2). Finally, based on the research results of Qi et al. [21], it is found that Gly32, Leu132 and Trp183 form strong hydrogen bonds in different parts of ZEN, especially the hydrogen bond force between water molecules and a variety of amino acids near the phenol hydroxyl group. It can be inferred that the substance is G/L/D/W-ZEN+H2O, and it is determined as the intermediate product I of the enzymatic hydrolysis reaction.  As shown in Figure 13, the mass spectrum with a retention time of 8.13 min implies that the parent ion of the substance is m/z 486.27. Further, by calculating the difference with the ZEN common characteristic fragment ion (m/z 301.14), we can get exactly one amino acid residue, tryptophan (Trp, W) in turn, with a neutral loss of Δ185. It may be that Gly32, Leu132 and Trp183 form strong hydrogen bonds in different parts of ZEN, especially the hydrogen bond force between water molecules and a variety of amino acids near the phenol hydroxyl group. It can be inferred that the substance is G/L/D/W-ZEN+H2O, and it is determined as the intermediate product I of the enzymatic hydrolysis reaction.  As shown in Figure 13, the mass spectrum with a retention time of 8.13 min implies that the parent ion of the substance is m/z 486.27. Further, by calculating the difference with the ZEN common characteristic fragment ion (m/z 301.14), we can get exactly one amino acid residue, tryptophan (Trp, W) in turn, with a neutral loss of Δ185. It may be As shown in Figure 12, the mass spectrum with a retention time of 6.30 min implies that the parent ion of the substance is double-charged m/z 808.61. Further, the difference is calculated by the fragment ions (m/z 752.06, m/z 695.52, m/z 622.97) to get exactly three amino acid residues, namely glycine (Gly, G), leucine (Leu, L) and tryptophan (Trp, W). The above amino acid residues correspond to neutral losses ∆57, ∆113 and ∆186, respectively. Further calculation of the molecular weight shows that the substance also contains aspartic acid residues (Asp, D). As aspartic acid belongs to an acidic amino acid, it is easy to lose a water molecule (H 2 O) and add an additional positive charge under the condition of collision-induced dissociation (CID), so the whole molecular system shows a double charge (m/z = 809, z = 2). Finally, based on the research results of Qi et al. [21], it is found that Gly32, Leu132 and Trp183 form strong hydrogen bonds in different parts of ZEN, especially the hydrogen bond force between water molecules and a variety of amino acids near the phenol hydroxyl group. It can be inferred that the substance is G/L/D/W-ZEN+H 2 O, and it is determined as the intermediate product I of the enzymatic hydrolysis reaction.
As shown in Figure 13, the mass spectrum with a retention time of 8.13 min implies that the parent ion of the substance is m/z 486.27. Further, by calculating the difference with the ZEN common characteristic fragment ion (m/z 301.14), we can get exactly one amino acid residue, tryptophan (Trp, W) in turn, with a neutral loss of ∆185. It may be that Trp contributes a proton H to combine with -OH of ZEN to lose molecule water. Finally, we can infer that the substance is W-ZEN-H 2 O, and it is determined as the intermediate product II of the enzymatic hydrolysis reaction.
It is reported that the structure of ZEN hydrolase can be divided into two parts: the folding center domain and the cap structure of the hydrolase [34]. A large groove is formed at the junction of the two structures, which is proved to be the substrate binding site by the substrate complex structure. The benzene ring part of the substrate ZEN is mainly fixed by hydrogen bonds, and the lactone ring part is mainly bound to the active center by hydrophobic force.
The pathway of ZEN degradation by ZEN hydrolase is shown in Figure 14. that Trp contributes a proton H to combine with -OH of ZEN to lose molecule water. Finally, we can infer that the substance is W-ZEN-H2O, and it is determined as the intermediate product II of the enzymatic hydrolysis reaction. It is reported that the structure of ZEN hydrolase can be divided into two parts: the folding center domain and the cap structure of the hydrolase [34]. A large groove is formed at the junction of the two structures, which is proved to be the substrate binding site by the substrate complex structure. The benzene ring part of the substrate ZEN is mainly fixed by hydrogen bonds, and the lactone ring part is mainly bound to the active center by hydrophobic force.
The pathway of ZEN degradation by ZEN hydrolase is shown in Figure 14. It can be inferred from Figure 14 that, under the action of ZEN hydrolase, a water molecule binds to the 4-hydroxy group of the benzene ring of ZEN to form the intermediate product Ⅰ, so that ZEN and the enzyme form a hydrogen bond fixation. After that, glycine, leucine and aspartic acid are removed, tryptophan is retained, and intramolecular dehydration forms intermediate product Ⅱ. Then, two water molecule bond intermediate product Ⅱ tryptophan is removed, and the lactone bond is opened to form water-added product C18H24O6. Finally, the water-added product loses a carboxyl group automatically to form the decarboxylation product C17H24O4.

Conclusions
The kinetics and mechanism of the enzymatic degradation of ZEN in DCO were investigated, and the following conclusions were obtained. It can be inferred from Figure 14 that, under the action of ZEN hydrolase, a water molecule binds to the 4-hydroxy group of the benzene ring of ZEN to form the intermediate product I, so that ZEN and the enzyme form a hydrogen bond fixation. After that, glycine, leucine and aspartic acid are removed, tryptophan is retained, and intramolecular dehydration forms intermediate product II. Then, two water molecule bond intermediate product II tryptophan is removed, and the lactone bond is opened to form water-added product C 18 H 24 O 6 . Finally, the water-added product loses a carboxyl group automatically to form the decarboxylation product C 17 H 24 O 4 .

Conclusions
The kinetics and mechanism of the enzymatic degradation of ZEN in DCO were investigated, and the following conclusions were obtained.
(1) The initial reaction rate increases firstly and then decreases with the increases in temperature and pH. It rises in a straight line and then tends to slow down with the increase in the dosage of enzyme added. Moreover, with the increase in substrate concentration, the initial reaction rate presents a positive proportion, which is a first-order reaction.
(3) Thermodynamic study shows that the activation energy (Ea) is 70.37 kJ·mol, the activation enthalpy change (∆H) > 0, the free energy (∆G) > 0, the activation entropy change (∆S) < 0 and the reaction cannot be spontaneous.
(4) The reaction mechanism of ZEN was studied by liquid chromatography mass spectrometry. It was found that ZEN first produced the intermediate G/L/D/W-ZEN+H 2 O, and then produced the intermediate W-ZEN-H 2 O under the action of the degrading enzyme. Then, the lactone bond was opened to produce the water-added product C 18 H 24 O 6 , and finally the decarboxylation product C 17 H 24 O 4 formed automatically.

Materials
Degummed corn oil (DCO) was obtained from Hebei Yufeng bio-engineering Co., Ltd. in Xingtai, China, and its contents of ZEN were 554 µg × kg −1 . Zearalenone hydrolase was obtained from Rongvezyme Preparation Co., Ltd., Zhuhai, China; the enzyme activity was 243 U × g −1 , which was measured according to the methods of Wang [19]. Zearalenone standard was purchased from Sigma-Aldrich (St. Louis, MO, USA). Tris, hydrochloric acid and other solvents and chemicals were all of analytical grade (Sinopharm Holdings Ltd., Shanghai, China).

Detection of ZEN Content
The concentrations of ZEN in the corn oil before and after processing were determined by high-performance liquid chromatography with a fluorescence detector (Waters 2695, Waters Co., Ltd., Milford, MA, USA) in accordance with GB 5009.209-2016 (National Standard of the People's Republic China, 2016).
Separation was conducted by a C 18 column (4 µm particle size, 150 mm × 4.6 mm, Waters Co., Ltd., Milford, MA, USA) at 25 • C, and the injection volume was set as 100 µL. A mobile phase consisting of acetonitrile-water-methanol [46:46:8, v/v/v] was used at a flow rate of 1 mL × min −1 .

Determination of the Initial Reaction Rate of ZEN Hydrolase
Degummed corn oil (500 g), 0.1 mol/L Tris HCI buffer (3 mL, pH 7.5) and 15 mg ZEN hydrolase were added to a beaker. Then, the sample was added to a high-speed shear machine (IKA ULTRA-TURRAX UTL2000, IKA Co., Ltd., Staufen, Germany) to mix the samples at 5000 r/min for 2 min. After that, the mixed sample (100 g) was put into the enzyme reactor at a certain temperature. After 10 min, the temperature was raised to 95 • C to terminate the reaction. Then, the sample was separated by a centrifuge at 10,000 r/min for 10 min. The upper oil layer was dried at 0.09 mPa and 80 • C. Finally, the content of the ZEN of the sample was determined by HPLC.
The initial reaction rate was calculated by the equation where V is the initial reaction rate (µmol × kg −1 × min −1 ); M 0 is the content of ZEN in DCO after the reaction (µg × kg −1 ); M 1 is the content of ZEN in DCO before the reaction (µg × kg −1 ); t is the reaction time (10 min); M is the molar mass of ZEN (318 g × mol −1 ).

Determination of Kinetic Constants and Thermodynamic Parameters
To study the effects of pH on the initial reaction rate of ZEN degradation, different buffers were used for pH 6.5-9.0.
To study the effect of ZHD dosage on the initial reaction rate of ZEN degradation, different ZHD dosages were used for 10-60 mg × kg −1 .
To study the effect of the concentration of ZEN on the initial reaction rate of ZEN degradation, different concentrations of ZEN were used for 500-7000 µg × kg −1 . Lineweaver-Burk plot and Wilkinson statistical method were applied to calculate the Michaelis constant (Km) and maximum velocity (Vmax).
The effect of reaction temperature on the initial reaction rate of ZEN degradation was investigated at a range of temperatures between 25 and 50 • C. The activation energy (Ea) was calculated from the slope of the Arrhenius plot according to the equation: lnk = −Ea/(RT) (4) where k, rate of activation at T; R, gas constant (8.314 J × mol −1 × K −1 ); T, the absolute temperature in Kelvin. The thermodynamic parameters for substrate hydrolysis were calculated according to equation [35].

Statistical Analyses
All experiments were carried out in triplicate independently and SPSS (version 19.0, SPSS, Inc., Chicago, IL, USA) was used for data analysis. Results are presented as mean ± standard deviation of triplicate measurements.