Competing elimination and substitution reactions of simple acyclic disulfides

Disulfides form persulfides at alkaline pH leading to potential overestimations in the cold cyanolysis method

It is well established that proteins and peptides can release sulfur under alkaline treatment, mainly through the β-elimination of disulfides with the concomitant formation of persulfides and dehydroalanine derivatives. In this study, we evaluated the formation of glutathione persulfide (GSSH/GSS-) by exposure of glutathione disulfide (GSSG) to alkaline conditions. The kinetics of the reaction between GSSG and HO- was investigated by UV-Vis absorbance, reaction with 5,5'-dithio-bis-(2-nitrobenzoic acid) (DTNB), and cold cyanolysis, obtaining an apparent second-order rate constant of ∼10-3 M-1 s-1 at 25 °C. The formation of GSSH and the dehydroalanine derivative was confirmed by HPLC and/or mass spectrometry. However, the mixtures did not equilibrate in a timescale of hours, and additional species, including thiol and diverse sulfane sulfur compounds were also formed, probably through further reactions of the persulfide. Cold cyanolysis is frequently used to quantify persulfides, since it measures sulfane sulfur. This method involves a step in which the sample to be analyzed is incubated with cyanide at alkaline pH. When cold cyanolysis was applied to samples containing GSSG, sulfane sulfur products that were not present in the original sample were measured. Thus, our results reveal the risk of overestimating the amount of sulfane sulfur compounds in samples that contain disulfides due to their decay to persulfides and other sulfane sulfur compounds at alkaline pH. Overall, our study highlights that the β-elimination of disulfides is a potential source of persulfides, although we do not recommend the preparation of GSSH from incubation of GSSG in alkali. Our study also highlights the importance of being cautious when doing and interpreting cold cyanolysis experiments.

Does Steric Hindrance Actually Govern the Competition between Bimolecular Substitution and Elimination Reactions?

Bimolecular nucleophilic substitution (SN2) and elimination (E2) reactions are prototypical examples of competing reaction mechanisms, with fundamental implications in modern chemical synthesis. Steric hindrance (SH) is often considered to be one of the dominant factors determining the most favorable reaction out of the SN2 and E2 pathways. However, the picture provided by classical chemical intuition is inevitably grounded on poorly defined bases. In this work, we try to shed light on the aforementioned problem through the analysis and comparison of the evolution of the steric energy (EST), settled within the IQA scheme and experienced along both reaction mechanisms. For such a purpose, the substitution and elimination reactions of a collection of alkyl bromides (R-Br) with the hydroxide anion (OH-) were studied in the gas phase at the M06-2X/aug-cc-pVDZ level of theory. The results show that, generally, EST recovers the appealing trends already anticipated by chemical intuition and organic chemistry, supporting the role that SH is classically claimed to play in the competition between SN2 and E2 reactions.

Gas-Phase Reactions of Dimethyl Disulfide with Aliphatic Carbanions - A Mass Spectrometry and Computational Study

Ion-molecule reactions of Me₂S₂ with a wide range of aliphatic carbanions differing by structure and proton affinity values have been studied in the gas phase using mass spectrometry techniques and DFT calculations. The analysis of the spectra shows a variety of product ions formed via different reaction mechanisms, depending on the structure and proton affinity of the carbanion. Product ions of thiophilic reaction (m/z 47), S_N2 (m/z 79), and E2 elimination - addition sequence of reactions (m/z 93) can be observed. Primary products of thiophilic reaction can undergo subsequent S_N2 and proton transfer reactions. Gibbs free energy profiles calculated for experimentally observed reactions using PBE0/6-311+G(2d,p) method show good agreement with experimental results. Graphical Abstract ᅟ.

Reactivity of hydropersulfides toward the hydroxyl radical unraveled: disulfide bond cleavage, hydrogen atom transfer, and proton-coupled electron transfer

Hydropersulfides (RSSH) are highly reactive as nucleophiles and hydrogen atom transfer reagents. These chemical properties are believed to be key for them to act as antioxidants in cells. The reaction involving the radical species and the disulfide bond (S-S) in RSSH, a known redox-active group, however, has been scarcely studied, resulting in an incomplete understanding of the chemical nature of RSSH. We have performed a high-level theoretical investigation on the reactions of the hydroxyl radical (˙OH) toward a set of RSSH (R = -H, -CH₃, -NH₂, -C(O)OH, -CN, and -NO₂). The results show that S-S cleavage and H-atom abstraction are the two competing channels. The electron inductive effect of R induces selective ˙OH substitution at one sulfur atom upon S-S cleavage, forming RSOH and ˙SH for the electron donating groups (EDGs), whereas producing HSOH and ˙SR for the electron withdrawing groups (EWGs). The H-Atom abstraction by ˙OH follows a classical hydrogen atom transfer (hat) mechanism, producing RSS˙ and H₂O. Surprisingly, a proton-coupled electron transfer (pcet) process also occurs for R being an EDG. Although for RSSH having EWGs hat is the leading channel, S-S cleavage can be competitive or even dominant for the EDGs. The overall reactivity of RSSH toward ˙OH attack is greatly enhanced with the presence of an EDG, with CH₃SSH being the most reactive species found in this study (overall rate constant: 4.55 × 10¹² M^-1 s^-1). Our results highlight the complexity in RSSH reaction chemistry, the extent of which is closely modulated by the inductive effect of the substituents in the case of the oxidation by hydroxyl radicals.

Thiol redox biochemistry: insights from computer simulations

Thiol redox chemical reactions play a key role in a variety of physiological processes, mainly due to the presence of low-molecular-weight thiols and cysteine residues in proteins involved in catalysis and regulation. Specifically, the subtle sensitivity of thiol reactivity to the environment makes the use of simulation techniques extremely valuable for obtaining microscopic insights. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of selected relevant issues in thiol redox biochemistry, such as investigations on (1) the protonation state of cysteine in protein, (2) two-electron oxidation of thiols by hydroperoxides, chloramines, and hypochlorous acid, (3) mechanistic and kinetics aspects of the de novo formation of disulfide bonds and thiol-disulfide exchange, (4) formation of sulfenamides, (5) formation of nitrosothiols and transnitrosation reactions, and (6) one-electron oxidation pathways.

Dramatic substituent effects on the mechanisms of nucleophilic attack on Se-S bridges

The reactions of XSeSX, XSeSY, and YSeSX (X, Y = CH3, NH2, OH, F) with F(-) and CN(-) nucleophiles have been investigated by means of B3PW91/6-311+G(2df,p) and G4 calculations. In systems where the two substituents are not identical (XSeSY), the more stable of the two possible isomers corresponds to those in which the most electronegative substituent is attached to Se. Nucleophilic attack takes place at Se, independent of the nature of the nucleophile, with the only exception being XSeSF (X = CH3 , NH2 , OH), in which case the attack occurs at S. In agreement with recent results for disulfide and diselenide linkages, the mechanisms leading to Se-S bond cleavage are not always the more favorable ones because for highly electronegative substituents the most favorable process is fission of the chalcogen-substituent bond. These dissimilarities in the observed reactivity pattern as a function of the electronegativity of the substituents are due to the fact that the σ-type Se-S antibonding orbital, which for low-electronegative substituents is the lowest unnoccupied molecular orbital (LUMO), becomes strongly destabilized when the electronegativity of the substituent increases, and is replaced by an antibonding π-type Se-X (or S-X) orbital. In contrast, however, with what has been found for disulfide and diselenide derivatives, the observed reactivity does not change with the nature of the nucleophile. The activation strain model provides interesting insight into these processes, showing that in most cases the activation barriers are the consequence of subtle differences in the strain or in the interaction energies.

Mechanism of S(H)2 reactions of disulfides: frontside vs backside, stepwise vs concerted

Density functional theory calculations indicate that the S(H)2 reactions of disulfides with alkyl or aryl radicals take place via concerted backside displacement. The activation energies for reactions of Me* with RSSR (R = Me, Et, (i)Pr, (t)Bu) increase with the size of R, since larger R groups prevent the formation of an ideal geometry for SOMO-LUMO overlap. Frontside transition states can also be located, but these lie at least 11 kcal mol(-1) above the corresponding backside transition states.

Frontside versus Backside S(N)2 substitution at group 14 atoms: origin of reaction barriers and reasons for their absence

We have theoretically studied the gas-phase nucleophilic substitution at group-14 atoms (S(N)2@A) in the model reactions of Cl(-)+AH(3)Cl (A=C, Si, Ge, Sn, and Pb) using relativistic density functional theory (DFT) at ZORA-OLYP/TZ2P. Firstly, we wish to explore and understand how the reaction coordinate zeta, and potential energy surfaces (PES) along zeta, vary as the center of nucleophilic attack changes from carbon to the heavier group-14 atoms. Secondly, a comparison between the more common backside reaction (S(N)2-b) and the frontside pathway (S(N)2-f) is performed. The S(N)2-b reaction is found to have a central barrier for A=C, but none for the other group-14 atoms, A=Si-Pb. Relativistic effects destabilize reactant complexes and transition species by up to 10 kcal mol(-1) (for S(N)2-f@Pb), but they do not change relative heights of barriers. We also address the nature of the transformation in the frontside S(N)2-f reactions in terms of turnstile rotation versus Berry-pseudorotation mechanism.

Mechanism of thiolate-disulfide interchange reactions in biochemistry

Both density functional theory (DFT) (B3LYP) and CCSD ab initio calculations were employed in a theoretical investigation of the mechanism of thiolate-disulfide exchange reactions. The reaction pathway for degenerate thiolate-disulfide exchange reactions with dimethyl disulfide has been shown to proceed through a SN2-like transition structure that is very close in energy to the corresponding trisulfur anionic intermediate ([delta-S-S-Sdelta(-)]). When relatively small substituents are involved, the level of theory must be increased to CCSD to make this rather subtle mechanistic distinction. With the more sterically hindered exchange reaction involving t-butyl mercaptide and di-t-butyl disulfide, the potential energy surface exhibits a distinct preference for the S(N)2 displacement pathway with an activation barrier of 9.8 kcal/mol. When corrections for solvent polarity are included (COSMO), an S(N)2 mechanism is also implicated in both polar and nonpolar solvents. DFT studies on thiolate-disulfide exchange, when the substituent is a model peptide, strongly support the intermediacy of a trisulfur intermediate that lies 10.7 kcal/mol below isolated reactants. A well depth of this magnitude should provide a sufficient lifetime of the intermediate to accommodate the requisite conformational adjustments that accompanies formation of the new disulfide bond.

Mechanism of SN2 disulfide bond cleavage by phosphorus nucleophiles. Implications for biochemical disulfide reducing agents

The B3LYP variant of DFT has been used to study the mechanism of S-S bond scission in dimethyl disulfide by a phosphorus nucleophile, trimethylphospine (TMP). The reaction is highly endothermic in the gas phase and requires significant external stabilization of the charged products. DFT calculations (B3LYP) were performed with explicit (water molecules added) and implicit solvent corrections (COSMO model). The transition structures for this SN2 displacement reaction in a number of model systems have been located and fully characterized. The reaction barriers calculated with different approaches for different systems are quite close (around 11 kcal/mol). Remarkably, the calculations suggest that the reaction is almost barrierless with respect to the preorganized reaction complex and that most of the activation energy is required to rearrange the disulfide and TMP to its most effective orientation for the SMe group transfer way. Different reactivities of different phosphorus nucleophiles were suggested to be the result of steric effects, as manifested largely by varying amounts of hindrance to solvation of the initial product phosphonium ion. These data indicate that the gas-phase addition of a phosphine to the disulfide moiety will most likely form a phosphonium cation-thiolate anion salt, in the presence of four or more water molecules, that provide sufficient H-bonding stabilization to allow displacement of the thiolate anion, a normal uncomplicated SN2 transition state is to be expected.

DFT study of the glutathione peroxidase-like activity of phenylselenol incorporating solvent-assisted proton exchange

Modeling of the glutathione peroxidase-like activity of phenylselenol has been accomplished using density-functional theory and solvent-assisted proton exchange (SAPE). SAPE is a modeling technique intended to mimic solvent participation in proton transfer associated with chemical reaction. Within this method, explicit water molecules incorporated into the gas-phase model allow relay of a proton through the water molecules from the site of protonation in the reactant to that in the product. The activation barriers obtained by SAPE for the three steps of the GPx-like mechanism of PhSeH fall within the limits expected for a catalytic system at physiological temperatures (DeltaG(1)++ = 19.1 kcal/mol; DeltaG(2)++= 6.6 kcal/mol; G(3)++ = 21.7 kcal/mol) and are significantly lower than studies which require direct proton transfer. The size of the SAPE network is also considered for the model of the reduction of the selenenic acid, step 2 of the GPx-like cycle. Use of a four-water network better accommodates the reaction pathway and reduces the activation barrier by 5 kcal/mol over the two-water model.