Abstract
Introduction
Direct thiol–disulfide interchange is the rate determining step in the folding process of proteins that have to form structural disulfide bonds. Although spontaneous thiol–disulfide interchange is slow (kinetically not competent on the folding timescale), enzyme catalysis accelerate the reactions in vivo. Chemically, similar mechanisms operate in the endoplasmic reticulum (ER; Fig. 1) (15) and mitochondrial inter membrane space (IMS; Fig. 2) (45) of eukaryotes and in the periplasm of prokaryotes (Fig. 3) (26, 68), although via different enzyme pathways. In all three cellular compartments, the new disulfide bond in the folding polypeptide is introduced by a series of inter- and intramolecular thiol–disulfide interchange reactions via oxidoreductases, and the oxidizing equivalents are supplied by molecular oxygen.



Protein thiols are also major targets of reactive oxygen species in vivo and oxidative conversion of thiols into disulfides and their subsequent reduction (e.g., see Figs. 4 and 5, respectively), is an alternative mechanism for thiol–disulfide exchange. In fact, the thiol–disulfide pool is thought to be primarily responsible for intracellular redox homeostasis and these reactions are important both for antioxidant defense and redox regulation of cell signaling (108). Therefore, delicate enzymatic pathways exist to control the direction of the electron-flow via cascades of one- and two-electron redox reactions (79).


It is now generally accepted that in cellular systems, both thiol–disulfide interchange and thiol-oxidation/reduction reactions are nonequilibrium dynamic processes, which are kinetically, not thermodynamically controlled (51, 53, 57). In other words, redox potentials and equilibrium constants only highlight whether a reaction is favorable, but partitioning of particular pathways depend on relative rates. Enzymes play essential roles in these processes via fine tuning reaction activation energies, which determine the final outcome of the oxidative stimuli or the position of structural disulfides in native proteins. This review is focused on the kinetics and mechanisms of direct substitution (the classical thiol–disulfide interchange) and redox reactions that result in thiol–disulfide exchange. Basic kinetic and mechanistic principles as well as parameters that influence reaction rates are discussed from a chemical perspective with relevant biological examples.
The Classical Thiol–Disulfide Interchange Reaction
From the chemical point of view, the reaction is a nucleophilic substitution of a thiol in disulfides with another thiol (Scheme 1, Reaction 1).

A number of studies have shown that the reaction proceeds with a simple SN2 type nucleophilic substitution mechanism. The SN2 mechanism is consistent with a one-step reaction via a single transition state complex with no intermediate formation (Fig. 6). The nucleophile is the deprotonated thiolate anion, which attacks the reacting sulfur of the disulfide moiety. Several pieces of theoretical and experimental evidence suggest a linear trisulfide-like transition state, with the negative charge being delocalized, but most abundant on the attacking and leaving sulfurs (10, 30, 32, 63).

In their seminal work, Fava et al. have reported that the rate equation for Reaction 1 exhibits first order dependency on both the thiol and disulfide concentrations (30), which is indicative of a bimolecular reaction. The observed base catalysis was consistent with the deprotonated thiolate being a better nucleophile compared to the thiol. In addition, their data allowed the exclusion of an alternative mechanism, where base hydrolysis of the disulfide would result in a sulfenic acid derivative, followed by condensation with the attacking thiol to give the product disulfide.
The noncatalyzed reaction is slow (in the order of k=0.1−10 M −1s−1 at pH 7, see Table 1) (30, 56, 93, 100), but oxidoreductase enzymes [such as thioredoxin (Trx), glutaredoxin (Grx), the bacterial disulfide bond protein family A-D (DsbA-D), protein disulfide isomerase (PDI), Mia40, ER oxidoreductin (Ero1), essential for respiration and vegetative growth sulfhydryl oxidase, Erv1, see Figs. 1 –3 and 5] are designed to accelerate the process up to k=104–106 M −1s−1 (for an abbreviated list of rate constants, see Table 1). A general feature of these enzymes is that they contain a Trx fold domain with a redox active CXXC motif at the active site, where in the reduced forms, the N-terminal Cys serves as the nucleophile. Therefore, the first bimolecular thiol–disulfide exchange reaction of a reduced oxidoreductase enzyme with its substrate results in an intermediate mixed disulfide species between the central disulfide sulfur of the substrate and the N-terminal Cys of the oxidoreductase. This is followed by an intramolecular nucleophilic attack of the C-terminal Cys on the N-terminal Cys of the oxidoreductase to give the oxidized CXXC motif and release the reduced substrate thiol. An exception is Mia40 (responsible for disulfide bond formation in the IMS see Fig. 2), which is not a member of the Trx family and has a redox active N-terminal CPC and two structural CX9C motifs. Another distinctive feature of Mia40 is that it forms the intermediate mixed disulfides via the C-terminal Cys of its CPC (12).
DTT, dithiothreitol; GSSG, oxidized glutathione; iP, phosphate buffer; GSH, reduced glutathione; nDsbD, N-terminal domain of DsbD; Grx, glutaredoxin; Trx, thioredoxin; I27S24-S55, I27 domain of human cardiac titin; HPLC, high-performance liquid chromatography, GOR, glutathione oxidoreductase; CoA, coenzyme A; CoAS-SG, glutathionylated CoA; GrxS-SG, gultathionylated GrX.
Factors That Influence the Rates of Thiol–Disulfide Interchange Reactions
In this section, a number of factors will be discussed, which were proposed to have rate enhancing or inhibitory effects on thiol–disulfide interchange reactions. Most reports consider the factors (i.e., the pKa and the nucleophilicity) that influence the reactivity of the attacking thiol. However, it is important to emphasize that from the kinetic point of view, factors that affect the stability of the leaving group or the electrophilicity of the central disulfide sulfur, will have just as important roles in decreasing the activation barrier.
Box 1: A Few General Kinetic Considerations for Studying Thiol–Disulfide Interchange Reactions
The rate equation and rate determining concentrations
For a true bimolecular thiol–disulfide interchange reaction (Scheme 1, Reaction 1), the following rate equation applies:
As shown by Equation 12, the rate law of a bimolecular thiol–disulfide interchange reaction exhibits first order dependency on both reactant concentrations. Under pseudo-first order conditions (one of the reactants is in >10-fold excess over the other), the observed kinetic traces will fit to an exponential equation. However, if the concentrations of the two reactants are comparable (which is often the case under biological conditions), then the data has to be treated with a more complicated second order equation.
Another important implication of this phenomenon is that under pseudo-first order conditions, the half life of the reaction will only be dependent on the concentration of the reactant in excess. In contrast, under second order conditions, the half life will depend on both reactant concentrations.
Reversible reactions
The above considerations apply when the reaction essentially proceeds to completion, and hence, is practically irreversible. However, thiol–disulfide exchange reactions in many cases approach equilibrium with substantial amounts of reactants remaining (such as in the case of the oxidation of DsbA or DsbC with DsbB or DsbD, respectively).
In this case, the rate law will have the following form:
and therefore, fitting the data with a single exponential (when, for example, one of the reactants is in excess over the other) does not result
The effect of pH
The kinetics of thiol–disulfide interchange reactions are pH dependent, because of the much stronger nucleophilicity of the thiolate compared to the thiol. The difference in the nucleophilicities is so dramatic that in many cases the reaction proceeds via the thiolate even when it is a minor species (i.e., at pH values far below the thiol pKa). A model that takes into account the different protonated forms of the attacking thiol is shown below:
From this model, the following rate equation can be derived:
Protonation–deprotonation reactions are usually fast in aqueous media, and hence, it is generally accepted that Reaction 15 can be treated using the pre-equilibrium assumption. In addition, if we assume that the reaction only proceeds via Reaction 17 under the applied experimental conditions and take into account that [RSH]tot=[RSH] + [RS-], then the rate equation will be:
Hence, the reactant concentration-independent apparent rate constant (
where k 17 is independent of pH. Practically this means that the rate of the reaction will increase with the pH up to the point where most of the attacking thiol will be deprotonated. Note that in polyprotic molecules, such as GSH, the deprotonated form of the sulfur center will be determined by microscopic acid dissociation constants (see Acid dissociation constants of reacting functional groups section), which further complicates the kinetic analyses (74, 78, 79).
Acid dissociation constants of reacting functional groups
Equation 20 in Box1 implies that at physiological pH, the rate of thiol–disulfide interchange reactions will be inversely dependent on the attacking thiol's pKa, because the deprotonated thiolate form is a better nucleophile than the thiol. However, due to the linear-free energy relationship, the nucleophilicity of the thiolate will increase with its pKa, which translates in an increase in k
17, representing an opposite effect on the overall reaction rate (93, 100, 107). The relationship between k
17 and the pKa of the reacting thiol is described by a Brønsted equation:
where βnuc is the Brønsted coefficient that was established to be ∼0.5 for the reaction of various different thiols with oxidized glutathione (GSSG) (100). Based on this Brønsted relationship, Figure 7 shows the predicted

The rate of Reaction 17 is also influenced by the pKa values of the central and leaving groups, which is described by a more comprehensive Brønsted equation:
where βnuc, βc, and βlg are Brønsted coefficients and pKa nuc, pKa c, and pKa lg are acid dissociation constants for the attacking, central and leaving groups (pKa RSH, pKa R′SH, and pKa R′′SH for Reaction 17), respectively.
The measured Brønsted coefficients (βnuc=0.59, βc=−0.40, and βlg=−0.59) indicate that the pKa values of the attacking and leaving thiols will have a larger influence compared to the central thiol. The negative βlg corroborates Wilson and coworkers' observation (107) that the rate decreases with the increase in the leaving thiol pKa, representing an opposite effect compared to that of the attacking thiol's pKa. This is explained by the increase in the stability of the leaving thiolate with the decrease of pKa lg, which results in a drop of the transition state energy. Therefore, the aforementioned extremely low pKa values of DsbA, DsbC, or Grx1 Cys residues will have an accelerating effect on the reduction of the mixed disulfide bonds (which form as intermediates during the reduction of: DsbA by DsbB, protein disulfides by DsbC, or glutathionylated proteins by Grx1) via leaving group stabilization in the transition state structure (see Fig. 6).
Microscopic versus macroscopic acid dissociation constants
The protonation state of a functional group in a polyprotic acid is characterized by microscopic acid dissociation constants, a phenomenon that is often overlooked in kinetic studies.
For glutathione, the two carboxylic groups are both deprotonated under physiological conditions, and therefore we can assume that their pKa values do not influence the kinetics of its thiol–disulfide exchange reactions. However, the terminal amino and sulfhydryl groups are mostly in their protonated forms at pH 7 resulting in two labile protons for reduced glutathione (GSH), see Equation 23 on Scheme 2a.

Ka 1 and Ka 2 (in Scheme 2a) are macroscopic acid dissociation constants that determine the overall protonation state of GSH. However, these equilibrium constants cannot be assigned to individual functional groups, because HL2− (Scheme 2a) represents a mixture of two protonation isomers [for a recent review of microscopic pKa's of Cys and GSH see (79)], where the proton is either attached to the thiol or the amino group of GSH (see Scheme 2b). Therefore, microscopic acid dissociation constants are required to calculate the concentration of thiolate functional groups. It is particularly important in kinetic studies because the rate of Reaction 1 will be more influenced by the protonation state of the thiol group than by the overall protonation state of GSH. However, protonation at a different site could not only influence thiol pKa values, but additional charge effects could also contribute to the overall reactivity of the molecule (79).
The ionizability of protein thiols is even more challenging to measure. These are influenced by many factors, including charged states, orientations, inductive effects of neighboring functional groups, and solvent accessibility (90). In addition, proteins are dynamic entities and changing the pH often results in small, but sometimes substantial conformational changes. Therefore, in a protein titration data set (which is in most cases the method of measuring the pKa), all of the above-mentioned variables make it hard to evaluate the sole contribution of the deprotonation of a particular functional group. This gives a different chemical meaning for protein pKa values compared to those of small molecules. On the other hand, it is important to understand how changes in pH affect the reactivities of protein thiols. Therefore, this parameter is best studied by measuring the pH profile of the kinetics of the reaction with its biologically relevant reaction partner [for more details see Ref. (79)]. For example, the pKa value of the C-terminal Cys in Trx is likely to be different when the enzyme is in its reduced state compared to when it is engaged in a mixed disulfide intermediate (between the N-terminal Trx Cys and the substrate) (19). Indeed, it has been proposed that due to a conformational change in the Trx–arsenate reductase mixed disulfide complex, new H-bonding interactions result in a transient drop in the pKa of the C-terminal Trx Cys, activating it for a nucleophilic attack on the N-terminal Cys to release the substrate (22, 91).
Steric factors
Steric factors play particularly important roles in the kinetics of reactions that proceed via an associative SN2 type mechanism, because of a crowded transition state structure. A linear conformation was suggested for the attacking, central, and leaving groups for the best orientation of their molecular orbitals see, for example, (10, 30, 33). Bulky functional groups can hinder the access of the attacking thiol and increase the activation energy. On the other hand, steric interactions can also stabilize the binding of particular substrates, which is an important factor in substrate specificity of many enzymes. A nice example for a steric hindrance-mediated kinetic insulation is the periplasmic disulfide shuttle system (Fig. 3). Although the pKa values of reduced DsbA and DsbC are very similar [pKa=3.5(82) and 4.1(99), respectively], in the periplasm, DsbA is predominantly oxidized, while DsbC is reduced. These functional redox statuses of DsbA and DsbC are maintained by DsbB and DsbD, respectively (see Fig. 3). Since all of these proteins are in the same cellular compartment and thermodynamically, the nonfunctional reduction of DsbA by DsbD or oxidation of DsbB by DsbC would be favorable, they are under strict kinetic control: DsbD reduces DsbC (k=4×106 M−1s−1)>3 orders of magnitude faster than DsbA (k=9×102 M−1s−1) and similarly DsbB oxidizes DsbA (k=3×106 M −1s−1)>3 orders of magnitude faster than DsbC (k=2×103 M −1s−1) (92), which were proposed to be due to steric factor-mediated substrate specificity (26, 49).
In addition, steric factors can introduce strain on a disulfide bond, and accelerate its interchange reactions by making it more labile. In DsbA, for example, the tense conformation of the oxidized form contributes to its oxidizing properties (112). It is also an important driving force of the folding process, where stabilization of the structured intermediate disulfide moieties are key in the overall kinetics and metastable disulfides are preferentially reshuffled (81).
An additional factor in thiol–disulfide exchange reactions is the formation of energetically favorable cyclic disulfides. This, together with the fact that cyclic disulfide formation is (i) entropy reduction facilitated and (ii) occurs at a larger effective thiol concentration compared to intermolecular reactions (due to the close proximity of the other thiol group), is the reason why dithiothreitol (DTT) is preferred over mercaptoethanol to reduce protein disulfides. In more detail, during the reduction of protein disulfides, the facilitated ring closure step results in a shift in the equilibrium toward the formation of oxidized DTT (a cyclic disulfide) and reduced protein thiols, while in the case of mercaptoethanol (even at large concentrations of the reducing agent), substantial amounts of mixed disulfide species can form. It is noteworthy that a ring opening reaction (the back reaction of the reduction of proteins with dithiols) would be faster for a more strained cyclic disulfide (as mentioned above for DsbA). However, investigating the kinetics of the reactions of many different dithiols with DTNB, Whitesides et al. reported a faster intramolecular interchange reaction (i.e., ring closure) step compared to the initial attack of the dithiol on the disulfide (in all studied cases even at very low reactant concentrations) (104). This would suggest that, even in cases when more strained intramolecular disulfides are generated, the equilibrium is shifted toward the formation of the cyclic disulfide species (when the dithiol is used in excess).
Similar principles would apply for the dithiol oxidoreductases, which also contributes to the thermodynamics and kinetics of their thiol–disulfide interchange reactions. In general, the rate of reactions of dithiols with disulfides and the reducing potentials of dithiols depend on: (i) the strain in the resulting cyclic disulfide, (ii) the conformational freedom of the dithiol, (iii) the entropy factor (which suggests a more unfavorable reaction for the formation of larger loops), and (iv) the steric proximity of the two thiol functional groups (particularly in proteins). These factors all contribute to the large differences that were observed in the equilibrium and rate constants for the reduction of disulfides by dithiols (e.g., orders of magnitude larger equilibrium constants were observed for the reactions of Trx compared to DsbA) (48, 115).
Charges, coulombic and H-bonding interactions with neighboring functional groups
In general, neighboring functional groups that stabilize the charge distribution of the transition state accelerate, and interactions that lower the energetics of the reactant complexes inhibit Reaction 1.
The aminoacid sequence of peptides was shown to have a profound effect on the rate of thiol–disulfide exchange reactions with particular influence of ionic aminoacids in close proximity to the thiolate (111). In proteins, due to their tertiary and quaternary structures, charged aminoacids that are sequentially well separated can also come close to the reaction center and affect reactivity. Positive or partial positive charges (electron withdrawing functional groups) stabilize the thiolate or the leaving group and, hence, catalyze the reactions, while negative charges have opposite effects.
Ion pairing, dipole or H-bonding interactions can also contribute to the nucleophilicity of the reacting thiolate. Electron-donating functional groups will increase its nucleophilicity, but H-bond donation will have an opposite effect (18, 33). For example, a recent report suggested an accelerating effect for the reaction of peroxiredoxins (Prx) with peroxides via a shift of H-bonds from the thiolate to the substrate to increase the nucleophilicity of the thiolate (33). Note that this latter reaction is an oxidation reaction, but the nucleophilicity of the thiolate in redox reactions is similarly an important rate-enhancing factor (see the Factors that influence the rate of thiol oxidation section).
For a direct thiol–disulfide interchange reaction, a hydrophobic environment is catalytic (32, 95). In Mia40, for example, two CX9C connected helices create a hydrophobic pocket, which serves as a substrate-binding groove (94). Similarly, a large hydrophobic region behind the active site of PDI was demonstrated to have an important role in its activity. In addition to preferential binding of unfolded substrates and assisting conformational changes, the hydrophobic environment was shown to lower the activation energy of the thiol–disulfide interchange reaction (32, 59). Furthermore, SN2 reactions in general favor an aprotic environment, partly because polar protic solvents tend to have a larger stabilizing effect on the reactants than on the transition state complex, which increases the activation energy. Thiol–disulfide interchange reactions indeed were found to be ∼3 orders of magnitude faster in aprotic solvents than in water (95), which is consistent with a more delocalized negative charge in the transition state, providing experimental credence for the theoretical calculations.
Inter- versus intramolecular reactions
The thiol–disulfide interchange can occur between two different reactants or in an intramolecular fashion. Intramolecular reactions can gain considerable rate enhancement compared to intermolecular reactions due to several different reasons (14, 64): (i) In an intermolecular reaction, a wide distribution of rotamers exist and statistical factors (as well as concentration) determine the probability of the conformation that undergoes reaction. On the other hand, in intramolecular reactions, the orientation of the reacting functional groups could be in the most favorable position for the molecular orbitals to overlap. (ii) In intramolecular reactions, interactions with neighboring residues could lock the reacting functional groups in a conformation that closely resembles the transition state and, thereby, considerably lower the activation energy. (iii) Intramolecular reactions are entropy reduction-facilitated via the loss of the translational entropy, which accompanies the bringing together of the reactants. For example, a highly favorable thiol–disulfide exchange between the C- and N-terminal sites of DsbD was proposed (k=1×105 M−1s−1), although the kinetics were investigated in an intermolecular fashion via the purified C-terminal (cDsbD) and N-terminal domains of DsbD (nDsbD; see Table 1). Therefore, the true intramolecular disulfide relay in DsbD is likely to be much faster due to the above-mentioned factors and unlikely to be rate limiting in the overall reduction of DsbC (although the direct reduction of DsbC by the nDsbD is fast k=3.9×106 M−1s−1) (92).
In addition, in intramolecular enzymatic reactions, the pKa-s of the reacting and leaving thiol groups could be optimized for the thiol–disulfide exchange via conformational changes as proposed for Trx [see the effect of pKa section and (19, 91)].
Grx represents a unique exception in the Trx-fold family of proteins, because the intramolecular disulfide formation between the two Cys in their CXXC motif is not favorable during the reduction of glutathionylated proteins (Fig. 5b). Grx recycle glutathionylated protein thiols (37, 50, 89, 114) via a Grx-GSH intermediate (see Fig. 5), which is subsequently reduced in an intermolecular reaction with another GSH to give GSSG and reduced Grx (31). Hence, the intramolecular Grx disulfide formation (which occurs occasionally as a side reaction) is outcompeted by the intermolecular reaction with GSH. Therefore, this is a nice example how preferential binding of a substrate (via steric, coulombic, H-bonding, etc. interactions) not only results in specificity, but has catalytic effect on the reaction (in this case, up to the extent to outcompete an entropy reduction facilitated intramolecular reaction) (50). Another factor that may contribute to the favorable intermolecular reduction of Grx-GSH is the stabilization of the bound GSH moiety (114) to inhibit the attack of the C-terminal Grx Cys via steric isolation. On the other hand, it is important to note that Grx also reduce protein disulfides (Fig. 5c). These reactions were proposed to occur via the so-called dithiol mechanism, which is different from the reduction of glutathionylated proteins. In the first step, a mixed disulfide forms as a result of a nucleophilic attack of the N-terminal Grx Cys on the protein disulfide moiety. This mixed disulfide is reduced in a Trx-like mechanism, that is, by the C-terminal Cys of the Grx CXXC motif. The Grx disulfide is recycled by GSH, where glutathionylation of the N-terminal Cys is followed by an intermolecular thiol–disulfide exchange with another GSH to give GSSG (just like in the case of the reduction of glutathionylated proteins) (65).
Mechanical force
Mechanical force can also trigger the thiol–disulfide exchange by activating the disulfide bond. For example, it was shown that in the adhesion ligand von Willebrand factor physiologically relevant fluid shear stresses promoted the thiol–disulfide exchange, which increased platelet binding (20).
Using force clamp spectroscopy, a series of publications from the Fernandez group allowed deeper insight into the molecular bases of thiol–disulfide exchange reactions (3, 63, 105). They have shown that the rate of Reaction 1 is exponentially dependent on the applied force when pulling the disulfide moiety apart. The rate is appropriately described by an Arrhenius equation (29):
where A is the pre-exponential factor, Ea is the activation energy of the reaction, k B is the Boltzmann constant, T is the temperature, and Δx r is the distance to the transition state along the reaction coordinate. The fitted value of Δx r holds information on the disulfide bond elongation although theoretical calculations suggest that this parameter should be handled with care (62). Nevertheless, the fact that the reduction of a disulfide by different thiol nucleophiles resulted in very similar Δx r values suggests that these reactions exhibit similar energetic and structural properties (2).
The kinetics and mechanisms of disulfide reduction by Trx was also investigated with this technique. The method shed light on the catalytic mechanism and suggested a force-dependent reorientation of the substrate at the active site to achieve a previously proposed linear -S-S-S- like transition state structure (105).
More recently, these investigators also managed to measure the rates of intramolecular thiol–disulfide exchange reactions (3). In the presence of two available thiol groups within the same protein, they could investigate their reactivities separately with the disulfide moiety. The obtained four-fold difference in the reactivities of the two thiol groups could be well explained by the previously discussed leaving group pKa effects. Importantly, this investigation for the first time allowed insight into the kinetics of thiol–disulfide isomerization reactions.
Prevention of nonfavorable side reactions
In complex biological systems, unwanted side reactions could inhibit the enzymatic catalysis. A good example is the IMS disulfide relay system, which is inhibited in the absence of GSH or the Zn2+ binding protein (Hot13), although via different mechanisms. Accumulation of intermediate mixed disulfides between Mia40 and its substrate slows the process down significantly. The counterproductive catalytic role of GSH (a reducing agent catalyzing an oxidation process) was explained by reduction of this mixed disulfide intermediate to recycle Mia40 in its reactive form (13). On the other hand, Zn2+ was shown to inhibit Mia40 reduction by Erv1 (67) and the Cys-rich Hot13 protein was proposed to catalyze the overall protein-folding process in the IMS via binding Zn2+ (70), and thereby accelerating Erv1-mediated Mia40 reduction (see Fig. 2).
Oxidation of Thiols
In many respects, thiol oxidation is similar to the thiol–disulfide interchange and by a variety of secondary reactions, oxidation of one thiol can lead to disulfide formation on another. Thiol oxidation to disulfide is at least a two-step process involving intermediate species. Oxidation can occur via one- and two-electron redox pathways depending on the nature of the oxidant.
The thiol-mediated redox state of the cell is predominantly controlled by the Trx and glutathione systems. Both these systems use NADPH as a fuel of their reducing capacity, although via different mechanisms. Trx are kept reduced by thioredoxin reductase (TrxR), which transfers the electrons from NADPH, while GSH is kept reduced by the glutathione oxidoreductase/NADPH couple. GSH and Trx can transfer their reducing equivalents to many different enzymes, including Prx, methionine sulfoxide reductase and ribonucleotide reductase for Trx; or glutathione peroxidase (Gpx or Orp1), glutathione S-transferase (GST), and Grx for GSH. Hence, they are involved in a variety of cellular events.
Two-electron oxidation pathways
Oxidation of Cys by two-electron oxidants, in most cases, occurs via formation of a Cys sulfenic acid (RSOH) derivative (Reactions 2–6 on Scheme 1). RSOH often forms directly in a bimolecular reaction with the oxidant (e.g., with peroxides; Reaction 2a). Alternatively, it can be a secondary hydrolytic intermediate in the oxidation of thiols with, for example, hypohalous acids (HOX, X=Cl, Br or I) or hypothiocyanite (OSCN−) (74, 76, 96) that predominantly react via X+, not OH+, transfer (34, 74, 101) (Reactions 2b–3). Sulfenic acids are generally reactive; they are rapidly reduced by thiols to form the corresponding disulfide species (Reaction 4).
Box 2: A Few General Kinetic Considerations for Studying Thiol Oxidation Reactions
The rate equation and rate determining concentrations
Two-electron thiol oxidation to disulfide [Scheme 1 model (
This implies that all subsequent reaction steps are fast and the steady-state approximation applies for the change in the concentration of the intermediates over time. Even in this case, the relationship between the observed rate constant and the true rate constant of the initial bimolecular reaction will depend on the observable (i.e., whether we follow the decay of the thiol or the formation of the disulfide). If we follow the decay of the thiol, then the observed rate constant needs to be corrected with a proportionality constant to account for the fast loss of a second thiol in Reaction 4, and hence, the rate law will have the following form:
In some cases, the situation can be more complicated and depending on the mechanism and conditions, the proportionality constant can be even larger (73). In these cases, a more comprehensive kinetic analysis is required to establish the concentration dependencies of the rate equation. Another important point is that if we follow the formation of the disulfide, we need to make sure that the rate determining step under the applied conditions does not change.
The effect of pH
The pH dependency of thiol oxidation reactions are generally more complicated compared to thiol–disulfide exchange reactions. For example, the pH profile of the rate law for Reaction 2 will not only depend on the protonation status of the thiol, but that of the oxidant's as well. Protonation of the oxidant accelerates the reaction because the protonated form is a better electrophile, which is an opposite effect to the aforementioned favorable deprotonation of the thiol. Therefore, the plot of k obs as a function of pH will have a bell-shaped profile, where the pH at the maximum will correspond to the average of the pKa values of the reactants [e.g., (78)]. Sometimes, alternative pathways via the less reactive species (e.g., the deprotonated form of the oxidant at higher pH) also exhibit competitive rates and, therefore, need to be taken into account (74, 78). Similarly to the kinetics of thiol–disulfide interchange reactions, microscopic acid dissociation constants (see Acid dissociation constants of reacting functional groups section) are required to describe speciation of thiolate derivatives of polyprotic molecules (such as Cys or GSH) (74, 78, 79). Furthermore, the rate of Reaction 4 is also pH dependent. As for the initial oxidation step (i.e., Reaction 2), the favored species are the thiolate (nucleophile) and the protonated RSOH (electrophile), and hence, a bell-shaped pH dependency is expected for this reaction too. In addition, the two consecutive reactions could exhibit pH dependencies to a different extent, which may result in a change in the rate determining step with the pH (74).
Studying fast reactions by competition kinetics
Under physiological conditions, some reactions are so fast that they cannot be followed directly. For a non-enzymatic reaction, the rate can be slowed down by a dramatic move away from pH 7 (either toward acidic or basic conditions), because of the above-mentioned reasons. However, most proteins are only stable at around pH 7, and pH-induced conformational transitions can result in changes in the mechanism. Therefore, to obtain kinetic information on fast protein oxidation reactions, the competition kinetics approach is often used [see e.g., (23, 77, 83, 87)]. This is based on the introduction into the reaction mixture of another reactant that competes for the oxidant. By changing the relative concentrations of the protein and the new reactant and taking into account the second order rate constant for the oxidation of the new reactant, the rate constant of interest (i.e., for the oxidation of the protein) can be calculated from product distribution. This approach is widely accepted, but there are potential pitfalls that are needed to be taken into account during the analysis: (i) Cross reactions of intermediates and products can change product distribution and confine the kinetic analysis. (ii) Not efficient mixing (i.e., slower mixing compared to reaction rates) can cause locally large concentrations of reactants due to inhomogenity, which could also change product ratios. (iii) For such competition kinetics measurements, the proportionality constant cannot be determined due to the lack of information on secondary reactions.
Peroxides react relatively slowly with most Cys residues with second order rate constants ranging from k=1 to a few hundred M −1s−1 (108) (e.g., see Table 2). However, thiol peroxidases contain Cys residues (Cysp; called the peroxidative Cys) that exhibit exceptional reactivity with peroxides having several orders of magnitude larger second order rate constants k=106–108 M −1s−1 [such as Prx; the Escherichia coli transcription factor (OxyR); or Gpx; see Table 2] (8, 23, 36, 61, 87). Based on a crude modeling exercise on the kinetics of cellular thiol oxidation by H2O2 in a cytosol-like environment, we recently put forward some concepts in relation to peroxide-mediated redox signaling (79, 108): (i) Prx1 and 2 and Gpx1 are the most favored targets of H2O2 and if they are effectively recycled they should consume the majority of the cytosolic peroxide. Due to the high reactivity of Prx and Gpx (and abundance in the case of Prx), direct oxidation of GSH and other important redox-regulated proteins (like GAPDH or protein tyrosine phosphatase) are prevented up till the concentration of H2O2 reaches tens of μM-s, and hence, the redox signal that involves H2O2 is most likely transmitted via Prx and/or Gpx.
OxyR, the Escherichia coli transcription factor; Prx, peroxiredoxins; PBS, phosphate-buffered saline; DTPA, diethylene triamine penta acetic acid; HRP, horseradish peroxidase.
The second part of the aforementioned simplest oxidation-mediated thiol–disulfide conversion mechanism is an electrophilic–nucleophilic condensation reaction between RSOH and another thiol to give the disulfide (Reaction 4, where RSOH is the electrophile and RSH is the nucleophile). It is significantly more challenging to investigate the kinetics of this reaction because RSOH derivatives are in most cases very reactive and do not exhibit characteristic spectral features (36, 74, 78). Often the production of a RSOH is rate limiting and much slower than its condensation with another thiol (66). In these reactions, the steady-state approximation applies to RSOH concentrations, which makes it difficult to acquire kinetic information for Reactions 4 (or 5) in small molecules (6, 7). However, in reactions of RSH with more reactive oxidants, such as the HOX [at pH=7 k HOCl=∼107, k HOBr=∼107 (4, 74, 84, 85)], Reaction 4 becomes rate determining and could be followed by stopped-flow spectroscopy at high pH (which is necessary to slow down the reaction to the stopped-flow timescale, see Box 2 The effect of pH section) (74). Extrapolation of the data to pH 7 gave an estimate of 105 M−1s−1 for the lower limit of the second order rate constant for the reaction of the parent CySOH with Cys under physiological conditions. (74). In contrast, in some proteins, steric hindrance restricts formation of a protein–protein disulfide and stabilizes the RSOH derivative. Notwithstanding, the RSOH derivative of human serum albumin exhibits surprisingly low reactivities even with small molecular weight thiols, having second order rate constants in the range of 2–100 M −1s−1 (102). In these cases, steric hindrance should be less prevalent and the proposed explanation is that the reacting Cys (CyS34) is sitting in a hydrophobic pocket and not readily accessible for the solvent and the reactants. (27).
RSOH can also act as a nucleophile, for example, in its self-condensation or disproportionation reactions to give thiosulfinate ester (RS-S(=O)R′) or sulfinic acid (RSO2H) derivatives, respectively. Both reactions were observed during the oxidation of free Cys with hypohalites even at >10-fold excess of thiol over the oxidant (74, 78). (RS-S(=O)R′) react with thiols to give a disulfide and RSOH (which is then converted to disulfide by another thiol) (39, 40, 73, 78). Therefore, these reactions can be considered as intermediate steps during thiol–disulfide conversion. However, it is important to note that in the presence of more than one type of thiol, intermediate formation of RS-S(=O)R′ can result in different product distribution, that is, different ratios of the mixed disulfides could form via the different pathways. On the other hand, RSO2H does not react with thiols in aqueous solutions and, therefore, diverts the reactions from the disulfide formation path. Notwithstanding, reduction of RSO2H derivatives of Prx back to the corresponding thiols by sulfiredoxin (Srx) occurs by a mechanism that involves the formation of a mixed (RS-S(=O)R′) intermediate (17, 54).
Thiol–disulfide conversion via radical-mediated oxidation pathways
One-electron oxidation of Cys thiols give thiyl radicals (Reactions 7 and 8) (103). The reaction proceeds via hydrogen abstraction from the thiol group or electron transfer through the thiolate species (9). The biologically important second order rate constants of these reactions range from 10–1000 M −1s−1 for superoxide [still debated due to the initiation of a chain reaction, see Refs. (52, 109, 117)] through 108 M −1s−1 (1) for carbon centered radicals up to the diffusion-controlled limit as with the hydroxyl radical. Rapid recombination of two thiyl radicals can give a closed shell disulfide species, which would complete the thiol–disulfide conversion (Reaction 9). However, although the reaction is fast [2k=1.5×109 M −1s−1 for GSH (46)], it is less likely to occur under physiological conditions because thiyl radicals exhibit high reactivities with other more abundant molecules. In biological systems, the most favorable is probably its reaction with thiolates (k>109 M −1s−1) (5) to give the disulfide anion radical (RSSR′•−; Reaction 10). RSSR′•− are powerful reducing radicals and react favorably with oxygen (which was predicted to be their predominant reaction under physiological conditions) to give superoxide and the corresponding disulfide (Reaction 11). Alternatively, RS. can react with oxygen (k ∼ 108–109 M −1s−1) in a reversible reaction to give a thiyl peroxyl radical RSOO•), although this reaction is predicted to be less prevalent in biological systems. The fate of RSOO• depends on the conditions, but many reaction pathways result in Cys derivatives with higher oxidation states (79).
In biological systems, one- and two-electron redox reactions often accompany each other. For example, in the disulfide relay systems, the oxidizing or reducing equivalents for oxidoreductases are supplied by molecular oxygen and NADPH, respectively, via one-electron redox reactions (see Figs. 1 –3 and 5). The switch between the one- and two-electron transport pathways is mediated via intramolecular flavin adenine dinucleotide–disulfide interactions in flavin containing oxidoreductases, such as Erv1, Ero1, Grx, or TrxR (11, 43, 106, 118).
Although one-electron-mediated pathways are prevalent in biological systems, discussion of factors that influence their kinetics is outside the scope of this review.
Factors That Influence the Rates of Thiol Oxidation
Many of the previously mentioned determinants of the rate of thiol–disulfide interchange, including the pKa and the nucleophilicity of the thiol (33), have similar influence on the kinetics of thiol oxidation reactions. In addition, the pKa of the oxidant in reactions 2 and 4 will also affect its electrophilicity (see Box 2) and, thus, have a marked effect on the kinetics [see e.g., (71, 72, 75, 76)].
As with the thiol–disulfide interchange, factors that enhance the stability of the leaving group will also accelerate oxidation reactions. For example, during the reduction of different peroxide species by thiols, the effect of the reacting and leaving group pKa values on the kinetics can also be described by Brønsted equations (33).
Reaction 2 was also suggested to proceed via an SN2-type mechanism with a linear conformation of the thiol sulfur and the oxidant's attacking and leaving groups (33). [Although a recent insightful computational study challenges this proposal (116)]. Therefore, similarly to the thiol–disulfide exchange, larger oxidants could react slower with protein thiols due to steric hindrance. This effect might be responsible for the observed drop of the reactivities of Prx 2 and 3 with the size of the peroxide oxidant (86).
H-bond donating functional groups have profound effects on thiol oxidation reactions from many perspectives. For example, in a recent study, high-level theoretical calculations (at the G4 level) showed that the rate of the reaction of a thiol with a peroxide species is markedly accelerated via H-bond donating interactions at the sulfur as well as both oxygen centers in a cooperative manner (77). H-bond donation to the reacting oxygen has an accelerating effect via increasing its electrophilicity, and H-bond donation to the leaving oxygen stabilizes the transition state (as well as the leaving group) and decreases the activation energy. However, it is a bit less intuitive to envisage how H-bond donation to the sulfur center accelerates the reaction, because in general, H-bonding interactions to a nucleophile decrease its nucleophilic character (18). The effect was best explained by a more pronounced stabilization of the transition state structure compared to the reactant complexes, which further lowers the activation barrier of the reaction (77). On the protein level, these H-bonding effects are demonstrated by the example of the peroxidase activity of Prx, where two highly conserved Arg residues (∼3.5 and 6.5 Å away from the reacting sulfur) were proposed to serve as H-bond donors to the sulfur center as well as the two peroxide oxygens (see Fig. 8). Comprehensive kinetic analysis with a series of mutant enzymes revealed the importance of these two Arg residues in the peroxidase activities of Prx2 and Prx3. Point mutation of either of the Arg residues resulted in five orders of magnitude drop in the second order rate constants of the reactions of both isoforms with H2O2. When both Arg residues were mutated (in Prx2 as well as Prx3), the reactivities decreased further by two orders of magnitude resulting in similar reactivities of the double mutants as free Cys. Note that similar rate-enhancing properties of neighboring functional groups via H-bonding interactions were proposed for Prx (although via different mechanisms) in two different studies based on structural considerations of the active sites (44, 80). Structural studies also suggest that the Arg that is in the close proximity of the Cp (Arg 123 for Prx 3, see Fig. 8) plays an important role in activating the reacting peroxide oxygen, but a different orientation of the peroxide was proposed (44). Hence, further kinetic and structural analysis is required to establish which functional groups are involved in H-bond donation to the leaving oxygen in the different Prx isoforms.

In OxyR, an Arg residue (Arg 266) in the close proximity of the reacting thiolate serves a similar role (21). For both Prx and OxyR, disulfide formation is accompanied by a conformational change, which takes the Cys residues away from the corresponding Arg to stabilize the disulfide moiety (i.e., the close proximity of Arg would lower the pKa and make the disulfide-bound Cysp, a better leaving group).
Side reactions can also inhibit oxidation-mediated thiol–disulfide exchange processes. For example, during catalytic reduction of peroxides by the Prx/Trx/TrxR/NADPH system, at high peroxide concentrations, RSO2H can form, which takes the protein out of the catalytic cycle (110, 113) (see Fig. 4). Although RSO2H are eventually recycled by Srx, this process is significantly slower and, hence, inhibitory to the overall reduction process of oxidized Prx.
Interconversion of Substitution and Redox-Mediated Thiol–Disulfide Exchange Pathways
A beautiful example for the interconversion of oxidative and substitution-mediated thiol–disulfide exchange is the regulation of hydroperoxide homeostasis in S. cerevisiae. The peroxidase Orp1 reduces the peroxide to give an Orp1-SOH derivative, which reacts with the Yes-associated protein 1 (Yap1) to give a mixed disulfide intermediate. The mixed disulfide is subsequently reduced in an intramolecular thiol–disulfide interchange reaction to release Orp1 and give the oxidized Yap1 (a monomer with an intramolecular disulfide), which is its activated form. The signal regulatory role of Trx in this process also occurs via thiol–disulfide interchange reactions with the intramolecular Orp1 or Yap1 disulfide to recycle them in their active forms (24).
Recent kinetic considerations suggested that Prx will also serve as cellular peroxide sentinels (79, 108) and might transmit the redox signal by initiation of a cascade of subsequent reactions potentially via similar thiol–disulfide-mediated mechanisms (35, 79, 88, 108): Prx senses the signal (H2O2) in a rapid bimolecular reaction to give Prx-SOH. Prx-SOH then could react with another protein thiol to form a mixed disulfide, which could further engage in inter- and intramolecular thiol–disulfide exchange reactions and transmit the signal. Alternatively, the oxidized Prx dimer (as in Fig. 4) can transmit the oxidizing equivalents to other proteins via thiol–disulfide interchange. An example for this latter mechanism is Trx oxidation by the Prx dimer (as in Fig. 5a), which potentially contributes to redox regulation of cell signaling, because Trx is involved as a thiol reducing agent in many signaling pathways.
Another example for the overlap of the oxidative and redox-mediated thiol–disulfide conversion is the capture of the H2O2 by Prx4 that is produced by Ero1 in the ER. Oxidized Prx4 can subsequently oxidize PDI (to give its active form) via thiol–disulfide interchange, which provides a more efficient utilization of oxidizing equivalents in the protein-folding process (see Fig. 1) (119).
Conclusions and Future Directions
Thiol–disulfide interchange and thiol oxidation reactions are constantly occurring in cellular systems. It is now widely appreciated that the dynamics of these reactions are governed via kinetic rather than thermodynamic forces. Due to the promiscuous chemical properties of sulfur, the reactions can proceed via diverse pathways and give a plethora of different products and intermediate species. The lack of distinct spectroscopic signatures of the various thiol derivatives, as well as the large scale of chemical factors that influence the reaction rates poses challenges for kinetic investigations. The situation is further complicated by the fact that in cellular systems, some of the reaction cascades travel through different cellular compartments (see e.g., Figs. 1 –3). Therefore, electron shuttling or diffusion/transport of reactants across membranes can become rate determining for the overall process. The above factors represent serious challenges and limitations to kinetic investigations that are aimed to better understand biological events. However, recent technological advances in kinetic, product identification/quantification, imaging, and other biochemical methodologies provide additional tools to carry out more comprehensive mechanistic studies. Therefore, a systematic rigorous chemical kinetics approach, which is routinely used in investigations of small molecule reactions, is becoming available to study more complicated enzymatic systems, which is expected to allow a deeper insight into cellular functions.
Footnotes
Acknowledgments
The financial support of FP7-PEOPLE-2010-RG (Marie Curie International Reintegration Grant; grant No.: PIRG08-GA-2010-277006) is greatly acknowledged. I am grateful to Prof. Christine Winterbourn for critical reading and for my colleagues for proofreading the manuscript.
