Virág Bogdándi,
Gábor Lente* and
István Fábián
Department of Inorganic and Analytical Chemistry, University of Debrecen, Debrecen, Hungary. E-mail: lenteg@science.unideb.hu
First published on 30th July 2015
The stoichiometry and kinetics of the reaction between isoniazid (isonicotinic hydrazide, INH) and the hypochlorite ion was studied in a basic aqueous medium (pH > 10) by spectrophotometric, hydrogen-1 nuclear magnetic resonance, electrospray mass spectrometric and stopped-flow methods. At oxidant excess, the sole organic product was confirmed to be isonicotinic acid. At isoniazid excess, a coupled condensation product formed from isoniazid and isonicotinic acid was also detected in relatively high amounts. The stoichiometric ratio of the two reactants (INH/OCl−) was 0.46 or 0.67 depending on which of the reagents was in excess initially. Under all conditions, the rate determining step of the overall reaction was confirmed to be the rapid bimolecular reaction between the neutral from of isoniazid and the anionic hypochlorite ion. The second order rate constant of the process was estimated to be (3.4 ± 0.7) × 104 M−1 s−1 at 25.0 °C.
Isoniazid was introduced into clinical practice in the 1950s. Since then, a lot of efforts have been devoted to study the molecular background of its specific antimycobacterial activity and the possible emergence of resistance to it. It seems now clear that the oxidation of isoniazid plays a major role both in its mechanism of action and metabolism in the body.2–16 The activation of isoniazid is carried out by a catalase-peroxidase enzyme: it involves coupling it to a NADH,4,6–8 and the resulting complex is a strong inhibitor of the enoyl-acyl carrier protein reductase enzyme,9 and consequently stops the synthesis of mycolic acid, which is required for building the bacterial cell wall. Because of this mechanism, there is still active research aimed at understanding the oxidation patterns of isoniazid under a variety of conditions.14–16
Hypochlorous acid and its conjugate base, hypochlorite ion, are ubiquitous in biological systems. This simple inorganic oxidant is produced by neutrophils to fight invading microbes and also has a role in inflammations.17–22 It is the product of the myeloperoxidase-catalyzed oxidation of chloride ion in the presence hydrogen peroxide17,22 and may react with a whole range of organic and inorganic reactants in the body. The kinetics of these processes has an obviously important role in governing the typical concentrations of available HOCl/OCl−.20,21,23,24 An additional, and industrially very important interest in the reactions of hypochlorite ion stems from the use of chlorination in water treatment technologies.23–28
In this paper, we report a detailed experimental study of the reaction between isoniazid and hypochlorite ion. In addition to the mentioned biological relevance, there are various points of chemical interest in this process as well. A flow injection determination of isoniazid was reported earlier based on chemiluminescence after reaction with electrogenerated hypochlorite ion.29–31 The oxidation methods of acid hydrazides are also of interest in organic chemistry typically as a tool in peptide synthesis or as a de-protection strategy to obtain carboxylic acids.32–38
At this point, it should be mentioned that there is an earlier kinetic study published on the reaction between isoniazid and hypochlorite ion in basic medium.39 In that work, the process was followed through sampling and offline analysis by iodometric titration. Reportedly, typical reaction lifetimes were several minutes or even hours. Our own observations contrast this earlier report39 very sharply. Three independent methods (spectrophotometry, NMR and mass spectrometry) showed that the reaction between isoniazid and hypochlorite ion is complete within 5 s. Furthermore, the earlier study reported a retarding effect of added chloride ion,39 which is unreasonable for hypochlorite ion reactions in the basic range as chloride ion does not participate in any reaction under these conditions.23,24 The earlier work39 assumed that hypochlorite ion can be selectively monitored in this system by iodometric titration. In fact, literature results show that isoniazid reacts with iodine readily.40 Therefore, it is impossible to determine the concentration of hypochlorite ion by this method in the presence of excess isoniazid. Consequently, the conclusions reported in that earlier study39 must be entirely flawed. One might speculate that a completely erroneous signal was followed by iodometry in the earlier work,39 which might have arisen because of the known influence of ambient oxygen during iodometric titrations at longer times. However, because sufficiently detailed description of the experimental procedures is not reported in the article,39 it is very difficult to assess what the exact origin of the misleading conclusion was.
The objective of this work was to study the reaction of isoniazid with hypochlorite ion through establishing the stoichiometry and the rate equation of the process. Another aim was to interpret this rate equation on a molecular basis. The reaction itself was expected to be a potential source of active isoniazid loss during tuberculosis treatment, when the body produces significant amount of hypochlorite ion to fight the infection, therefore contributing to the interpretation of the well-known need for unusually high therapeutic doses.1–5 In addition, characterizing the process was expected to yield more insights into the metabolism of isoniazid.
Chloride ion free solutions of sodium hypochlorite were prepared as described earlier in the literature.23–25,41 The concentrations of NaOCl solutions were determined by iodometric titration,23–25,41 whereas the concentration of excess NaOH in the same solution was determined by acid-base titration with perchloric acid as detailed in a previous paper.24 In all kinetic measurements, constant ionic strength was maintained by adding an appropriate amount of NaClO4. High purity sodium perchlorate was specially prepared for this purpose by a literature method.42 Solutions of sodium hydroxide were prepared from solid NaOH and their concentrations were standardized by routine acid–base titration against primary standard potassium hydrogen phthalate (Reanal). NaCl (Sigma-Aldrich) was used as an additive in some experiments. For NMR experiments, D2O and NaOD (40% in D2O) were purchased from Sigma-Aldrich and used without further purification.
ESI-MS measurements were carried out with a Bruker micrOTOFQ mass spectrometer equipped with a quadrupole and time-of-flight (Q-TOF) analyzer in positive ion mode. The nebulizer pressure was set to 6 psi, N2 was used as drying gas at 4 l min−1, the temperature was maintained at 160 °C, and the applied capillary voltage was −3.0 kV. The typical error of mass measurement was below 0.005 in the m/z range from 80 to 600. Spectra were accumulated and recorded by a digitizer at a rate of 2 GHz.
Stopped-flow experiments showed that the redox reaction between the two reagents was quite fast: spectral changes were over within 5 s in all cases. Fig. 1 shows typical spectral changes in the UV range during the course of the process. Under these particular conditions, the absorbance change is complete within 2 s. The ratio of the initial reactants in this experiment is close to stoichiometric as it is confirmed later. Both hypochlorite ion and isoniazid have absorptions in this spectral region, so the initial spectrum reflects the combination of these absorbance contributions. The final spectrum is dominated by the organic final product as chloride ion does not have light absorption in this spectral region. Further preliminary spectrophotometric experiments showed that it is possible to monitor the concentration change of isoniazid with high selectivity at 320 nm. Therefore, this wavelength was selected for later kinetic studies.
The results were evaluated using a spectrophotometric titration plot, in which the absorbance was corrected for dilution and plotted as a function of the molar ratios of the two reactants.48,49 The dilution is defined as η = (V + Vini)/Vini, where Vini is the initial volume of the titration (2.00 cm3 in our case), whereas V is the total volume of the added titrant (up to 1.50 cm3 in our experimental setup). When ηA is plotted as a function of the molar ratio of the reactants (as in Fig. 2), a combination of two straight lines is expected for an irreversible process, and the intersection of the two straight lines gives the stoichiometric ratio. This general statement is proved by a detailed derivation that can be found in the ESI. Fig. 2 presents the results of such an analysis of our titration data. Two sets of points are shown as the titration was carried out in both ways. At first, isoniazid was titrated with sodium hypochlorite and the break point occurred at [OCl−]/[INH] = 1.49 ± 0.03. When a solution of NaOCl was titrated with an INH solution in the second case, the break point was observed at [INH]/[OCl−] = 0.46 ± 0.01. It should be added that these plots depend on the wavelength used, but the position of the break point does not. For Fig. 2, the wavelength selection was made to ensure that the break point is determined with the highest accuracy possible. Fig. S1 and S2 in the ESI† show the entire data set recorded during the two titrations, whereas Fig. S3† shows a plot analogous to Fig. 2 but using different wavelengths. Fig. S4† also proves that the absorbance measured in the kinetic experiment shown in Fig. 1 at 260 nm correlates linearly with the absorbance measured at 320 nm. Therefore, it is clear that these two wavelengths carry the same kinetic information. This is also true for any other wavelengths where meaningful absorbance change can be measured.
In a simple process, it would be expected that the break point of the first titration occurs at the reciprocal of the break point of the second titration (about 2.2). This is clearly not the case here. Yet the plots themselves are combinations of straight lines as expected for a simple reaction. The significant difference between the two titrations is that in the first case (where the break point occurs at 1.49), isoniazid is titrated with hypochlorite ion gradually. Therefore, the process goes under isoniazid excess after the addition of each aliquot. In the other case, hypochlorite is the excess reagent when the reaction proceeds. A simple interpretation of our observations is that the product (or product mixture) is different in these two cases.
Earlier literature results about the oxidation of hydrazide derivatives show that three different oxidized products and three different stoichiometries might be plausible in this case.32–38 These are as follows:
IN–CONHNH2 + 2OCl− + OH− → IN–COO− + N2 + 2H2O + 2Cl− | (1) |
2IN–CONHNH2 + OCl− → 2IN–CONH2 + N2 + H2O + Cl− | (2) |
2IN–CONHNH2 + 2OCl− → IN–CONH–NHCO–IN + N2 + 2H2O + 2Cl− | (3) |
Eqn (1) is the four-electron oxidation of isoniazid producing isonicotinic acid (INA, the anion is shown because of the basic medium). In eqn (2), a formally one-electron oxidation produces isonicotinic amide (INM). Eqn (3) shows a coupled product, 2-(4-pyridinylcarbonyl)hydrazide (IND), as a result of a formally two-electron oxidation. It seems that under conditions where OCl− is in excess during the process, eqn (1) with a molar ratio of 1:
2 INH
:
OCl− dominates and isonicotinic acid may be the only organic product. Under reverse conditions, however, the lower (ca. 1
:
1.5) INH
:
OCl− implies that at least one of the products of eqn (2) and (3) must be formed in significant quantities. As the independent experiments (described in the next paragraphs) never confirmed any formation of isonicotinic amide, it is concluded that the coupled product shown in eqn (3) is responsible for this difference. At oxidant excess, there is some slight deviation from the 0.50 molar ratio expected based on eqn (1), which is 8% and, given the estimated standard errors, probably significant. This slight deviation might be due to the fact that some of the nitrogen atoms of isoniazid are oxidized to nitrate ion:
IN–CONHNH2 + 7OCl− + 3OH− → IN–COO− + 2NO3− + 3H2O + 7Cl− | (4) |
An INH/OCl− ratio of 0.46 would imply about 9:
2 N2
:
NO3− final molar ratio in the nitrogen-containing products at oxidant excess, which means that nitrate ion would be a very minor by-product.
As isonicotinic acid (INA) and isonicotinic amide (INM) are available independently, some experiments were carried out with these substrates directly. UV-spectrophotometric experiments revealed no change at all for an hour when isonicotinic acid (INA) was mixed with sodium hypochlorite in aqueous solution (Fig. S5†). On the other hand, isonicotinic amide (INM) underwent significant change in a few minutes when mixed with hypochlorite ion in aqueous solution (Fig. S6†). Therefore, INM is not a plausible product of isoniazid oxidation under these conditions.
Further efforts to identify the products were done by ESI-MS spectrometry. Isoniazid and sodium hypochlorite were mixed in a dilute solution under conditions when INH was in excess. The ESI-MS spectrum of the resulting mixture is shown in Fig. 3. The spectrum confirms the presence of excess isoniazid (the peak at 138.065 corresponds to INH + H+, whereas the peak at 160.047 matches INH + Na+). Two peaks attributable to isonicotinic acid are also visible (146.021: INA + Na+; 168.002: INA − H+ + 2Na+) confirming the identity of this product. As the experiment was done at an excess of isoniazid, formation of the coupled product IND is possible, and is actually confirmed by four peaks at the higher mass end of the spectrum (243.084: IND + H+; 254.078: 2IND + H+ + Na+ (2 + ion); 265.064: IND + Na+; 287.046: IND − H+ + 2Na+). Signals attributable to isonicotinic amide INM are missing from the spectrum.
![]() | ||
Fig. 3 ESI-MS spectrum of the product mixture at isoniazid excess. The peak marked with * is a known instrument background peak that does not belong to the sample. |
Further investigation of the products was carried out by 1H-NMR spectroscopy. Some of the results are shown in the ESI as Fig. S7.† The 1H-NMR spectrum of isoniazid in D2O in basic medium shows two doublets: the first at 7.56 ppm, the second at 8.40 ppm. When an excess of NaOCl is added, the structure of the spectrum remains the same, but the doublets appear at 8.42 ppm and 8.90 ppm and the coupling constant increases. When solid isonicotinic acid is added to this sample, no new peaks appear in the spectrum. Therefore, isonicotinic acid is the sole organic product of the oxidation at an excess of hypochlorite ion, which agrees well with the results of the spectrophotometric titration.
In summary, when hypochlorite ion is used in excess, isoniazid is oxidized only to isonicotinic acid. At isoniazid excess, a mixture of isonicotinic acid and the coupled molecule IND is formed. It should be added that the relative amount of IND seems to be dependent on the conditions. A higher absolute concentration of isoniazid seems to favor the formation of IND. Unfortunately, kinetic studies were not feasible in the region where IND is a significant product as the formation of nitrogen bubbles interfered with the spectrophotometric detection.
As pseudo-first order conditions could not be applied, the evaluation was based on determining the initial rate50,51 of absorbance change at 320 nm. This wavelength proved to be selective to isoniazid under the experimental conditions used in this work. The dependence of initial rates on the concentration of isoniazid is shown in Fig. 5 at two different pH values. In both cases, the points define a reasonably good straight line without an intercept. It follows that the process is first-order with respect to the organic reactant. A similar concentration dependence for hypochlorite ion is shown in Fig. 6. Again, the rates are directly proportional to [OCl−] at both hydroxide ion concentrations. These observations prove that the order of the reaction is 1 with respect to both reagents and the overall rate law is second order:
v = kapp[INH][OCl−] | (5) |
![]() | ||
Fig. 5 Initial rate as a function of isoniazid concentration in the oxidation of isoniazid with hypochlorite ion. [OCl−] = 1.0 mM, T = 25.0 °C, l = 1.00 cm. |
![]() | ||
Fig. 6 Initial rate as a function of hypochlorite ion concentration in the oxidation of isoniazid with hypochlorite ion. [INH] = 0.50 mM, T = 25.0 °C, l = 1.00 cm. |
In this equation, kapp is a pH-dependent apparent rate constant that also incorporates the molar absorptivity of isoniazid because v is defined as the initial rate of absorbance change. It should be emphasized that [INH] here means the total concentration of all different protonated forms of isoniazid.
The effect of chloride ion was also tested as it is a product of the reaction. As shown in Fig. S8 of the ESI,† addition of Cl− prior to initiating the reaction does not have any effect on the initial rate and does not influence the kinetic curves, either. This clearly confirms the expectations, as chloride ion usually only has effects on the reactions of hypochlorous acid in acidic medium, where the equilibria involving the formation of chlorine (Cl2) are relevant.23–25,41
As shown in Fig. 4, a change in the concentration of hydroxide ion strongly influences the shape of the detected traces. As [OH−] increases, the initial absorbance of the kinetic curves increases, but the initial rate of absorbance change decreases. These two phenomena can be interpreted by a single acid–base equilibrium. In addition, it is also to be noted that the final absorbance does not decrease to zero. The pH-dependence of the final absorbance was essentially the same as the pH dependence of the initial absorbance, so the final absorbance reflects that some unreacted isoniazid remains after reaction. The initial ratio of reactants was 1:
2 (INH
:
OCl−), so this observation gives some support to the assumption that the experimentally determined stoichiometry (0.46 INH
:
OCl−) is somewhat different from that shown in eqn (1).
The change in the initial absorbance values is a consequence of a reaction that occurs within the dead time of the stopped-flow instrument (about 1 ms). This is clearly not the redox reaction, which is slower. It was also shown that the spectrum of isoniazid (without the addition of any oxidant) depends on the pH in this range (Fig. S9 in the ESI†). Isoniazid is known to take part in two distinct acid–base equilibria in the measurable pH range (2.0–12.0).52 The first equilibrium (i.e. the one with the lowest pK) corresponds to the protonation of the pyridine moiety and is characterized by a pK of 3.53 at 25.0 °C in aqueous solution containing 20% of methanol.52 Thus, this protolytic step is only significant in the acidic pH range. The second step is a deprotonation of the amidic nitrogen of isoniazid, and is characterized by the following equilibrium:
IN–CONHNH2 = IN–CON−NH2 + H+ | (6) |
The corresponding acid dissociation constant is defined as follows:
![]() | (7) |
It should be pointed out that the notation [IN–CONHNH2] was introduced here to mark the concentration of the neutral form of isoniazid, in contrast with [INH], which means the total concentration of all protonated forms. The literature value of the pK is 11.14,52 which is within our investigated pH range. Therefore, the acid–base equilibrium shown in eqn (6) could be responsible for the observations. The forms IN–CONHNH2 and IN–CON−NH2 have different spectra, this causes the pH-dependence of the initial absorbance values. Using standard considerations, the following formula can be derived for the initial absorbance:
![]() | (8) |
![]() | (9) |
The dependence of the initial absorbance on [OH−]T (where [OH−]T is the analytical concentration of sodium hydroxide) is shown in Fig. 7. The value of pKw, the ionic product of water, was fixed at 13.77, which is the literature value for 1.0 M NaClO4 medium.53 In the initial attempts to fit the experimental data, AlowpH was found to be very close to zero and its standard deviation was larger than the fitted value. It was concluded that AlowpH is not a significant term and its value was fixed at zero. The remaining two parameters were estimated by least squares fitting as follows: AhighpH = 1.87 ± 0.02 and pK = 11.16 ± 0.03. The results show that the form IN–CON−NH2 absorbs strongly, but IN–CONHNH2 does not have significant absorption at the measurement wavelength (320 nm). The fitted pK is in excellent agreement with the literature value of 11.14 despite the fact that the two values were determined in slightly different media.
![]() | ||
Fig. 7 Initial absorbance as a function of pH in the oxidation of isoniazid with hypochlorite ion. [INH] = 0.50 mM, [OCl−] = 1.0 mM, T = 25.0 °C, l = 1.00 cm. |
The pH-dependence of the initial rate can be interpreted by considering two different factors. First, the apparent molar absorbance of isoniazid changes as a result of equilibrium 6. A simple way to correct for this effect is to divide the initial rate by the initial absorbance. Second, the rate of isoniazid consumption changes with shifting pH as the reactivities of the two different forms are not the same. The usual derivation, which is based on considering reaction 6 as a pre-equilibrium, gives the following results:
![]() | (10) |
In this formula, rlowpH is characteristic of the reaction of the IN–CON−NH2 form, whereas rhighpH belongs to the IN–CON−NH− form. Fig. 8 shows the a plot of v0/At=0 as a function of [OH−]T (see eqn (9) and (10) regarding the calculation of [OH−]). The first attempts at fitting gave very low values for rhighpH, which were again smaller than the estimated standard deviation. Therefore, the contribution of this term to the overall rate is never significant and its value was fixed at zero in the final evaluation. The other two parameters were estimated as follows: rlowpH = 29.4 ± 1.4 s−1 and pK = 11.00 ± 0.04 (using 13.77 as pKw again53). The estimated pK is in acceptable agreement with the one estimated from the initial absorbances and also with the literature value.
![]() | ||
Fig. 8 Absorbance-corrected initial rate as a function of pH in the oxidation of isoniazid with hypochlorite ion. [INH] = 0.50 mM, [OCl−] = 1.0 mM, T = 25.0 °C, l = 1.00 cm. |
These results imply that hypochlorite ion only reacts with the neutral form of isoniazid. The rate determining step is a simple second order process:
![]() | (11) |
The value of the bimolecular rate constant (k) can be calculated from the initial rate of absorbance change at 320 nm (v0) from any of the experiments using the following formula:
![]() | (12) |
Based on the total kinetic data set, k = (3.4 ± 0.7) × 104 M−1 s−1 was estimated. The detailed calculations are given in Table S1 of the ESI.†
At this point, it should be noted that all of the previous results were derived from measuring initial rates. For a complete kinetic description, it needs to be proved that the rate equation obtained in this way is suitable for describing the full kinetic traces. The kinetic curves of the overall second order process between reactants R1 and R2 under non-flooding conditions can be described as follows:50,51
![]() | (13) |
When the initial ratio of the reactants is exactly stoichiometric ([R1]0 = [R2]0/μ), the time dependence of the concentration of R1 is given as follows:50,51
![]() | (14) |
For any other initial conditions, the formula can be given as:
![]() | (15) |
The absorbance can be calculated from Beer's law using εR1, εR2 and εPi for the molar absorptivities of R1, R2 and products Pi, in order:
A = εR1[R1] + εR2[R2] + ΣεPi[Pi] | (16) |
The time dependence of the absorbance for the case described by eqn (14) is therefore:
![]() | (17) |
The parameters P1, P2 and P3 are related to the other quantities through the following equations:
![]() | (18) |
When one of the two reactants is in excess (i.e. eqn (15) is valid), the following formula is valid for the absorbance:50,51
![]() | (19) |
P1, P2, P3 and P4 are parameters in these equations that are related to the initial concentration of the reactants and the molar absorptivities of the reactants and products as follows:
![]() | (20) |
For the measurements where the initial hydroxide ion concentration was high enough (≥10 mM) to be practically constant throughout the entire experiment, these equations were used to fit the full kinetic traces measured. As shown by Fig. 9 and Fig. S10 in the ESI,† the fits were excellent. Some of the numerical results of these fittings are also displayed in the ESI as Table S2.† The excellent fit of these integrated expression shows that the rate equation established based on the initial rates describes the entire studied process, no secondary reactions emerge at longer reaction times. However, these integrated expressions can only be used for curves for which the initial hydroxide ion concentration is at least 10 mM, otherwise [OH−] also changes significantly during the process, and the full integration of rate law (cf. eqn (10)) is not possible analytically any more.
![]() | ||
Fig. 9 Example of fitting experimental data to eqn (19). [INH] = 0.50 mM, [OCl−] = 1.15 mM, [OH−]T = 100 mM, T = 25.0 °C, l = 1.00 cm. The grey line represent the fitted curve. The letter A represents absorbance. |
Footnote |
† Electronic supplementary information (ESI) available: Figures referred to in the main text of the article. See DOI: 10.1039/c5ra13540k |
This journal is © The Royal Society of Chemistry 2015 |