Mónika Galajdaa,
Tímea Fodorb,
Mihály Purgela and
István Fábián*b
aMTA-DE Homogeneous Catalysis and Reaction Mechanisms Research Group, Debrecen, H-4032, Hungary
bDepartment of Inorganic and Analytical Chemistry, University of Debrecen, Debrecen, H-4032, Hungary. E-mail: ifabian@science.unideb.hu
First published on 7th January 2015
The oxidation of pyruvic acid by hypochlorous acid was studied in the pH 1.0–11.0 range. The main path in the redox process is the reaction between the pyruvate ion and HOCl. It was shown that the reaction is second order in both reactants and kHOCl = 2.17 ± 0.04 M−1 s−1 (I = 1.0 M (NaClO4), T = 25.0 °C) with ΔH≠HOCl = 34.6 ± 1.2 kJ mol−1 and ΔS≠HOCl = −120.3 ± 3.6 J mol−1 K−1 for this reaction step. The relatively large negative entropy of activation is consistent with an O atom-transfer mechanism. DFT calculations support this conclusion by predicting a concerted rearrangement of the activated complex which includes the reactants and the solvent water molecule. It was also confirmed that the hydration of pyruvic acid has a significant contribution to the observed kinetic features under acidic conditions. The hydration reaction is a proton catalyzed process and the pseudo-first-order rate constant can be interpreted by considering the protolytic and hydration equilibria of pyruvic acid. The rate constants were determined for the non-catalytic and the proton catalyzed paths: kh0 = 0.24 ± 0.01 s−1 and kh1 = 5.5 ± 0.2 M−1 s−1.
In a broader context, the results also bear relevance on the redox chemistry of the reactants. The reactions of pyruvate ion have extensively been studied in recent years due to the specific role of this species in various biological redox processes.17–20 Furthermore, the simplicity of the molecule allows modelling the oxidation of organic substrates in the presence of various oxidants.21–26 Pyruvic acid is a simple substrate with a carbonyl group next to a carboxylic group. As such, its oxidation in aqueous solution may yield either acetaldehyde or acetic acid. With strong oxidants, such as the hypochlorous acid, the formation of the latter product is expected.
Hypochlorous acid is considered to be a strong disinfecting agent which is broadly used in water treatment technologies27–32 but also plays an important role in oxidative stress in vivo systems.33–40 Due to its outstanding oxidative power, this species can significantly alter the activities and availabilities of biomolecules. Hypochlorous acid may act either as a chlorinating agent or as an O-donor oxidant leading to the formation of very different types of products. It is a key issue how these alternative reaction paths occur under environmental and biological conditions. A survey of the literature reveals that with a few exceptions29,31,41–45 only limited reliable information is available on the intimate details of the oxidation reactions of hypochlorous acid.
In this study we provide a thorough account on the kinetics of the title reaction and elaborate the mechanistic details also using DFT calculations. The results are expected to contribute to better understanding the redox chemistry of the reactants, in particular, the features of HOCl as an O-donor oxidant.
HgO + 2Cl2 + H2O = 2HOCl + Hg2Cl2 | (1) |
Unoxidized chloride was converted into HgCl2 and the mixture of the HOCl/HgCl2/HgO slurry was purged with a steady stream of argon gas into a 0.4 M NaOH solution.
2HOCl = Cl2O + H2O | (2) |
Cl2O + 2NaOH = 2NaOCl + H2O | (3) |
The stock solutions of NaOCl were stored at 5 °C in the dark and were standardized by iodometric titration before use.48 The concentration of excess NaOH in the stock solution was determined by pH titration with standard HClO4.
Other reagents (Na-pyr Aldrich, KI Reanal, NaOH Aldrich) were of reagent grade quality and used without further purification. Samples were prepared using doubly deionized and ultrafiltered water obtained from a MILLI-Q RG (Millipore) water purification system.
The kinetic measurements were made at 25 ± 0.1 °C and 1.0 M ionic strength set with NaClO4 prepared from HClO4 (Reanal) and Na2CO3 (Reanal).49 The concentration of NaClO4 was determined by drying a known volume of the stock solution in an oven at 120 °C and weighing the mass of the recovered solid. The pH was adjusted with perchloric acid or NaOH solutions.
Spectrophotometric measurements were made with a Hewlett-Packard 8543 UV/VIS diode array spectrophotometer (equipped with a built-in magnetic stirrer) using different illumination protocols to explore the possibility of unwanted photoreactions.50 Photochemical interference was not observed in this system. Measurements were made in stoppered quartz cells of 10 mm light path. The temperature of the cell was controlled by a built-in thermoelectric Peltier device.
Conventional fast kinetic measurements were performed in aqueous solution under pseudo-first-order conditions with an Applied Photophysics DX-17MV sequential stopped-flow instrument using 10 mm optical path at 220 and 325 nm. Pseudo-first order rate constants were obtained from the average of five replicate kinetic runs. Quenched stopped-flow measurements were made with a HI-TECH RQF-63 instrument. Data fitting was made with the program packages SCIENTIST51 and ORIGIN52 using non-linear least-squares routines.
1H NMR spectra were recorded by using a Bruker DRX 400 (9.4 T) NMR spectrometer equipped with a Bruker VT-1000 thermo-controller and a BB inverse z gradient probe (5 mm). Pyruvic acid solutions (0.5 M) were prepared in H2O (D2O was added to the samples in capillary). The chemical shifts are reported in ppm, using DSS (4,4-dimethyl-4-silapentane-1-sulfonic acid) as an external standard. The spectra were analyzed with the Bruker WinNMR software package.
CH3COCOOH + HOCl = CH3COOH + CO2 + H+ + Cl− | (4) |
![]() | (5) |
The reaction became faster by decreasing the pH and eventually it could not be monitored spectrophotometrically because the spectral effect became less pronounced. The protolytic equilibrium of the oxidant is shifted toward the acidic form under pH 7.0 and HOCl has relatively small absorbance compared to hypochlorite ion. Thus, the absorbance change was not measurable. An attempt was made to use pyruvate ion in excess and to monitor the reaction at lower wavelengths. However, due to the large absorbance of pyruvate ion, these attempts were not successful. In addition, when the reactant concentrations were increased significantly, the formation of carbon dioxide led to the evolution of small bubbles in the cell which interfered with spectrophotometric monitoring of the reaction.
In order to avoid the noted technical problems, the reaction was studied by using the quenched stopped flow method in the neutral – slightly acidic region. This kinetic technique is more laborious than the standard stopped-flow method and it has sporadically been used for studying inorganic reactions.61,62 The method is based on immediate chemical quenching of the studied reaction, thus it requires the use of a quenching agent, which is capable to remove the reactant(s). The quenching process should be much faster than the reaction studied and should not trigger unknown side reactions. Once the reaction is stopped, the reaction mixture can easily be analysed for the products of either the main reaction and/or the quenching reaction. In the latter case, the concentration of a species formed in equivalent amount form the reactant is determined at a given reaction time. In other words, the amounts of unreacted reactants and product can be determined at various reaction times.
The shortest cumulative dead-time of the HI-TECH RQF-63 unit used in this study is below 4 ms. This time is required to start up a reaction and subsequently freeze it by administering an appropriate quenching agent. Therefore, the instrument can conveniently be used for monitoring reactions with pseudo first order rate constants up to 100 s−1. Pyruvate ion was used in large excess in the quenched stopped-flow experiments and 0.5 M NaI solution was used as quenching agent. Iodide ion instantaneously reduced hypochlorous acid to chloride ion and an equivalent amount of iodine formed. Kinetic traces were recorded by stopping the oxidation reaction at various reaction times and determining the amount of iodine formed by standard potentiometric titration with thiosulfate solution. A representative kinetic curve is plotted in Fig. 1.
![]() | ||
Fig. 1 Typical kinetic trace obtained by the quenched stopped-flow method. C0PYR = 0.500 M, C0HOCl = 3.29 × 10−2 M, pH = 3.88, I = 1.0 M (NaClO4), T = 25.0 °C. Quenching agent: 0.500 M NaI. |
The pseudo-second-order rate constant is shown as a function of pH in Fig. 2. The results obtained from the regular spectrophotometric and the quenched stopped-flow experiments complement each other remarkably well.
![]() | ||
Fig. 2 The pH dependence of the pseudo-second order rate constant for the oxidation of pyruvate ion by hypochlorous acid. Experimental rate constants were obtained by the quenched stopped-flow method (open circles) and conventional kinetic experiments (filled circles). The result of fitting to eqn (7) is also shown (solid line). I = 1.0 M (NaClO4), T = 25.0 °C. |
Pyruvate ion is the dominant form of the substrate above pH 4.0, thus the characteristic pH dependence of the rate constant can be interpreted by considering the fast protolytic equilibrium of hypochlorous acid (eqn (6)).
![]() | (6) |
The pseudo-second-order rate constant is expressed as follows:
![]() | (7) |
Parameter | Value | Unit | Remark |
---|---|---|---|
a This work. | |||
log![]() |
7.40 | Ref. 47 | |
7.40 ± 0.03 | a, pH-metric titration | ||
7.37 ± 0.01 | a, fitting the data to eqn (7) | ||
kOCl | (3.97 ± 0.05) × 10−2 | M−1 s−1 | a |
kHOCl | 2.14 ± 0.02 | M−1 s−1 | a, slightly acidic–alkaline region |
2.20 ± 0.04 | M−1 s−1 | a, acidic region | |
2.17 ± 0.04 | M−1 s−1 | a, average value | |
k′HOCl | (9.1 ± 2.0) × 10−2 | M−1 s−1 | a, acidic region |
ΔH≠HOCl | 34.6 ± 1.2 | kJ mol−1 | a |
ΔS≠HOCl | −120.3 ± 3.6 | J mol−1 K−1 | a |
log![]() |
1.76 ± 0.03 | M−1 | a |
The oxidation by HOCl is the main path of the reaction. The activation parameters for this path were determined by studying the temperature dependence of kobs at pH 2.8 in the 10–45 °C temperature range. Under such conditions, pyruvic acid is fully deprotonated and the contribution of the OCl− path is negligible, i.e. kobs ≈ kHOCl. The activation parameters, ΔS≠ and ΔH≠, were calculated by fitting the kobs values to the Eyring equation (eqn (8)) using a non-linear least squares fitting algorithm (Fig. S2, ESI†). The results are listed in Table 1.
![]() | (8) |
According to literature data, logKp = 1.97 for pyruvate ion,63 thus, it is reasonable to assume that one of the reaction steps is somehow associated with the protonation of this species. In order to explore this possibility, kinetic traces were recorded by acidifying neutral solutions of pyruvate ion in the absence of hypochlorous acid. First order absorbance changes were observed and the first-order rate constants agreed well with that of the faster reaction step in the oxidation experiments at the same pH and pyruvate ion concentration. However, simple proton transfer reactions are diffusion controlled and the protonation of pyruvate ion to pyruvic acid is expected to complete far within the dead-time of the stopped-flow experiments under the conditions applied. It follows that these observations cannot be interpreted by considering only the protonation step, instead, it needs to be assumed that the absorbance change is associated with a slower process triggered by the protonation of pyruvic acid.
A thorough inspection of the kinetic traces confirms this assumption. As shown in Fig. 3, there is a sudden absorbance drop from A0 to A1 within the dead-time of the instrument followed by a slower process associated with an additional absorbance change saturating at A2.
It is well documented that α-ketocarboxylic acids are involved in keto–enol and hydration equilibria.63–66 Scheme 1 shows all the possible equilibrium steps together with the corresponding equilibrium constants for the pyruvate system.63
![]() | ||
Scheme 1 Acid–base, keto–enol and hydration equilibria in aqueous solution of pyruvate ion. Reproduced from ref. 63. |
According to the data in Scheme 1, the enolic form of pyruvic acid and pyruvate ion is present in less than 0.1 and 0.01%, respectively and the hydrated form of pyruvate ion is also negligible. In contrast, the hydration of pyruvic acid is substantial, and the keto and hydrated forms are present in comparable concentrations under acidic conditions. 1H NMR studies support this overall model and the pH dependence of the CH3 peak of pyruvic acid is consistent with the simultaneous presence of the keto and hydrated forms (Fig. S3, ESI†). These features offer a possibility for coherent interpretation of the experimental observations by considering reactions (9) and (10).
![]() | (9) |
![]() | (10) |
![]() | (11) |
The plots of A1 and A2 as a function of pH are consistent with such protolytic steps and allow the calculation of Kp and Kh by fitting A1 and A2 simultaneously to the following expressions:
![]() | (12) |
![]() | (13) |
![]() | ||
Fig. 4 The characteristic absorbance values upon mixing the solution of pyruvate ion with perchloric acid, cf. Fig. 3. The experimental data (dots and squares) were fitted simultaneously to eqn (12) and (13) (solid lines) C0PYR = 5.00 × 10−3 M, I = 1.0 M (NaClO4), T = 25.0 °C. |
Parameter | Value | Unit | Remark |
---|---|---|---|
a This work. | |||
KPYRp | 93.3 | M−1 | Ref. 63 |
32.5 ± 2.0 | M−1 | a | |
Kh | 1.2 ± 0.1 | a | |
K′p | 72.2 | M−1 | a, calculated value, eqn (11) |
57.2 ± 3.4 | M−1 | a, fitted value, eqn (16) | |
kh0 | 0.24 ± 0.01 | s−1 | a |
kh1 | 5.5 ± 0.2 | M−1 s−1 | a |
εpyr | 256.6 ± 3.4 | M−1 cm−1 | a |
εHpyr | 31.8 ± 2.1 | M−1 cm−1 | a |
εHpyr(H2O) | 27.3 ± 2.8 | M−1 cm−1 | a |
The calculated protonation constant of pyruvate ion is about three times smaller than the value reported earlier.63 This is partly due to the difference in the experimental conditions applied, however, other factors should also be considered. The protonation constant was estimated by introducing various approximations in the evaluation of the experimental data in the earlier study. This could easily lead to biased results. In our case, the equilibrium constants were calculated directly from stopped-flow photometric experimental data which are expected to provide more reliable results.
The pseudo-first-order rate constant of hydration increases with increasing [H+] indicating that the reaction of pyruvic acid is acid catalyzed (Fig. 5). The corresponding rate constants can be expressed as follows:
kh = kh0 + kh1[H+]![]() | (14) |
![]() | ||
Fig. 5 The pH dependence of the pseudo-first-order rate constant for the hydration of pyruvic acid: experimental data (dots) and the fitted results on the basis of eqn (15) (solid line). λ = 240 nm, C0PYR = 5.00 × 10−3 M, I = 1.0 M (NaClO4), T = 25.0 °C. |
The protonation equilibrium of pyruvate ion is expected to be diffusion controlled and the pH dependence of the pseudo first order rate constant of hydration, khobs, is given by the following expression:
![]() | (15) |
The forward and reverse rate constants of hydration were calculated by fitting khobs to eqn (15) using fixed values for Kp and Kh. The results are listed in Table 2. In an earlier study, the rate of dehydration of hydrated pyruvic acid was found to be dependent on the pyruvate ion concentration.64 This is in sharp contrast with our finding because the rate constant is independent of pyruvate ion concentration in the 5.0 × 10−4–1.0 × 10−2 M concentration range. Furthermore, according to our results the proton catalysed step is considerably slower and the non-catalytic reaction is much faster than it was reported earlier. It is difficult to identify the main reasons for such differences but several factors should be considered. The technical setup of the kinetic experiments was less sophisticated in the cited study and it is not clear how the relatively fast reactions could be monitored with any acceptable reliability. In addition, the kinetic runs were made with pre-acidified pyruvic acid solutions which may undergo slow decomposition. If the reagent solutions were prepared long before the kinetic experiments, the products of the acidic decomposition could easily interfere with the kinetic studies. Such problems were circumvented in this study by using neutral pyruvate ion solutions as a reactant.
The pH dependence of the rate of the hydration of pyruvic acid can be interpreted by considering the structural features of α-keto carboxylic acids. The acid catalysis confirms that pyruvic acid is involved in the formation of a protonated species along the hydration process. The protonation most likely occurs at the carbonyl group of pyruvic acid, and the formation of an internal hydrogen bond to the carboxylic group may stabilize this intermediate. Such a structure may enhance the possibility of a nucleophile attack by a water molecule at the partially positive carbon atom of the carbonyl group (Scheme 2).
![]() | ||
Scheme 2 The proposed structure of the possible intermediate in the proton catalyzed hydration of pyruvic acid. |
Under acidic conditions, the fast step in the pyruvate ion–HOCl reaction is consistent with the results obtained for the hydration of pyruvic acid. It follows that the slower step in Fig. 6 corresponds to the oxidation of pyruvate ion. The hydration and protonation reactions should be considered as fast pre-equilibria with respect to the redox process. The pseudo first order rate constant is linearly dependent on the concentration of both pyruvate ion and HOCl concentration when they are used in large excess over the other reactant (Fig. S4, ESI†). Thus, the redox reaction remains first order for each reactant under acidic condition. The pH dependence of kobs clearly demonstrates that the reaction between the deprotonated form of the substrate and HOCl is the main reaction path in the redox process and the reaction with Hpyr has minor contribution (Fig. 7).
![]() | ||
Fig. 6 Typical kinetic trace for the oxidation of pyruvic acid by hypochlorous acid under acidic conditions. C0PYR = 1.00 × 10−2 M, C0HOCl = 1.00 × 10−3 M, pH = 1.50, I = 1.0 M (NaClO4), T = 25.0 °C. |
![]() | ||
Fig. 7 The pH dependence of the pseudo-first-order rate constant for the oxidation of pyruvic acid by hypochlorous acid under acidic conditions: experimental data (dots) and the fitted results on the basis of eqn (18) (solid line).). λ = 325 nm, C0PYR = 2.00 × 10−2 M, C0HOCl = 1.00 × 10−3 M, I = 1.0 M (NaClO4), T = 25.0 °C. |
Since the hydration of pyruvic acid is much faster than the redox reaction, the pH dependence shown in Fig. 7 can be interpreted by considering the formal protonation constant of pyruvic acid, K′p. The pseudo-first order rate constant, kobs, was fitted to the following expression:
![]() | (16) |
The pH profile of the second-order rate constant for the oxidation of pyruvate ion by hypochlorous acid in the pH = 1–11 range is shown in Fig. 8.
Recently, we have demonstrated that O-transfer steps are associated with relatively large negative entropies of activation in simple oxidation reactions of oxone.68 This is most likely the consequence of significant geometrical rearrangements associated with the reaction. An O-transfer step would also require massive geometrical changes in the pyruvate ion hypochlorous acid reaction and, therefore, be consistent with a large negative entropy of activation which, in fact, was observed in this system (Table 1). This indicates that a simple O-transfer is probably also operative in the pyruvate ion–hypochlorous acid reaction. The formation of acetic acid requires that the incoming oxygen attacks the carbonyl group and the actual reaction path may include some sort of an intramolecular interaction between the reactants, which is followed by concerted rearrangement of the chemical bonds.
The DFT calculations lend support to this assumption. Energy profiles were calculated for the reaction between pyruvate ion and HOCl. Alternatively, one water molecule was also added to the model checking the possible role of the solvent molecule. The optimized transition states are shown in Fig. 9. Substantial differences between the structural arrangements of the reactants in the transition state were not found in the absence and presence of a water molecule. The calculations predict that the reaction is thermodynamically highly favored (Table S2, ESI†). The oxidation is a concerted process where the bond-lengths of the reactants are changing simultaneously. Thus, the HO–Cl and the Ccarboxyl–Ccarbonyl bonds are elongated somewhat while the Ccarbonyl–OH bond is shortened along the reaction coordinate and the water molecule assists this process via hydrogen bonds.
![]() | ||
Fig. 9 The transition states in the pyruvate ion–HOCl reaction in the absence and presence of one water molecule calculated at the M06/6-311+G(d,p) level. |
As expected, the oxidation involves the attack of the O-atom of HOCl on the carbon atom of the carbonyl group. The results can easily be conceptualized by considering that the carbon atom is partially positive and as such it is prone to an attack by the partially negative moiety of the oxidant. This assumption fully supported by the calculated ESP charges for the transition state (Fig. S5, ESI†). The structure of the transition state suggests a very plausible interpretation of the mechanism of the reaction. The rearrangement of the activated complex can easily lead to simultaneous formation of a C–O bond at the carbonyl site, the breakage of the C–C bond, the release of the chloride ion and a proton transfer from HOCl to the water molecule. This concerted rearrangement of the primary and secondary chemical bonds explains the relatively large entropy of activation observed experimentally. Calculations at various theoretical levels do not predict significant differences in the activation Gibbs free energies (ΔG≠) (Table 3, the corresponding calculated activation parameters are listed in Table S3, ESI†). The added water molecule decreases the ΔG≠ values and makes more exergonic the process. This effect is more pronounced when we use M06 functional and larger basis sets. The calculations predict reasonable structures for the transition state, and explain the stabilization effect of the solvent. However, there is a significant difference between the experimental and calculated activation parameters. Calculations at higher level of theory may resolve this problem by including more water molecules and considering the effects of the excessive network of hydrogen bonds in the solvent. The structures of the transition states and the related imaginary frequencies are given in the ESI (Table S4†).
Functional | Basis set | Pyruvate ion–HOCl | Pyruvate ion–HOCl–H2O |
---|---|---|---|
B3LYP | 6-31+G(d,p) | 133.7 | 129.2 |
6-311+G(d,p) | 124.6 | 121.9 | |
M06 | 6-31+G(d,p) | 127.3 | 113.7 |
6-311+G(d,p) | 119.8 | 101.9 |
Footnote |
† Electronic supplementary information (ESI) available: Figures showing the spectra of the reactants, the temperature dependence of kobs at pH 2.8, the 1H NMR spectra of pyruvic acid as a function of pH, kobs as the function of HOCl concentration at pH 1.0, ESP charges for the transition state; tables on the stoichiometry of the reaction, and DFT results on the Gibbs free energies, calculated activation parameters and coordinates of the transition state. See DOI: 10.1039/c4ra12789g |
This journal is © The Royal Society of Chemistry 2015 |