Multiscale Modeling of Thiol Overoxidation in Peroxiredoxins by Hydrogen Peroxide

Inorganic Polysulfides in Solution: Structural Properties and Conformational Isomerism

We present a comprehensive theoretical examination of the structural properties of dianionic polysulfides [Sn]2- (n = 2-6), their conjugated monoacids [HSn]- (n = 2-6), and a selection of 1e--oxidized radical anions [Sn]•- (n = 2-4), in aqueous and dimethyl sulfoxide (DMSO) solutions. We investigated the structures and stabilities of various conformational isomers within these families of compounds by employing Quantum Mechanics-Molecular Mechanics (QM-MM) Molecular Dynamics (MD) simulations. The explicit inclusion of solvent molecules in the calculations revealed stable conformational structures that were previously unreported and might have appreciable concentrations in real systems. The interconversions between the isomeric structures proceed on the order of hundreds of picoseconds and are energetically similar to the isomerization processes in substituted cyclohexanes. We also conducted a detailed analysis of the stability of different isomers of the radical anion [S4]•- in solution. Our findings highlight the significant influence of the solvent on the isomerizations, a result that could be particularly relevant for enhancing the performance of metal-sulfur batteries.

Glutathione and peroxisome redox homeostasis

Despite intensive research on peroxisome biochemistry, the role of glutathione in peroxisomal redox homeostasis has remained a matter of speculation for many years, and only recently has this issue started to be experimentally addressed. Here, we summarize and compare data from several organisms on the peroxisome-glutathione topic. It is clear from this comparison that the repertoire of glutathione-utilizing enzymes in peroxisomes of different organisms varies widely. In addition, the available data suggest that the kinetic connectivity between the cytosolic and peroxisomal pools of glutathione may also be different in different organisms, with some possessing a peroxisomal membrane that is promptly permeable to glutathione whereas in others this may not be the case. However, regardless of the differences, the picture that emerges from all these data is that glutathione is a crucial component of the antioxidative system that operates inside peroxisomes in all organisms.

Minimum Free Energy Pathways of Reactive Processes with Nudged Elastic Bands

The determination of minimum free energy pathways (MFEP) is one of the most widely used strategies to study reactive processes. For chemical reactions in complex environments, the combination of quantum mechanics (QM) with a molecular mechanics (MM) representation is usually necessary in a hybrid QM/MM framework. However, even within the QM/MM approximation, the affordable sampling of the phase space is, in general, quite restricted. To reduce drastically the computational cost of the simulations, several methods such as umbrella sampling require performing a priori a selection of a reaction coordinate. The quality of the computed results, in an affordable computational time, is intimately related to the reaction coordinate election which is, in general, a nontrivial task. In this work, we provide an approach to model reactive processes in complex environments that does not require the a priori selection of a reaction coordinate. The proposed methodology combines QM/MM simulations with an extrapolation of the nudged elastic bands (NEB) method to the free energy surface (FENEB). We present and apply our own FENEB scheme to optimize MFEP in different reactive processes, using QM/MM frameworks at semiempirical and density functional theory levels. Our implementation is based on performing the FENEB optimization by uncoupling the optimization of the band in a perpendicular and tangential direction. In each step, a full optimization with the spring force is performed, which guarantees that the images remain evenly distributed. The robustness of the method and the influence of sampling on the quality of the optimized MFEP and its associated free energy barrier are studied. We show that the FENEB method provides a good estimation of the reaction barrier even with relatively short simulation times, supporting that its combination with QM/MM frameworks provides an adequate tool to study chemical processes in complex environments.

Structural and functional fine mapping of cysteines in mammalian glutaredoxin reveal their differential oxidation susceptibility

Protein-S-glutathionylation is a post-translational modification involving the conjugation of glutathione to protein thiols, which can modulate the activity and structure of key cellular proteins. Glutaredoxins (GLRX) are oxidoreductases that regulate this process by performing deglutathionylation. However, GLRX has five cysteines that are potentially vulnerable to oxidative modification, which is associated with GLRX aggregation and loss of activity. To date, GLRX cysteines that are oxidatively modified and their relative susceptibilities remain unknown. We utilized molecular modeling approaches, activity assays using recombinant GLRX, coupled with site-directed mutagenesis of each cysteine both individually and in combination to address the oxidizibility of GLRX cysteines. These approaches reveal that C8 and C83 are targets for S-glutathionylation and oxidation by hydrogen peroxide in vitro. In silico modeling and experimental validation confirm a prominent role of C8 for dimer formation and aggregation. Lastly, combinatorial mutation of C8, C26, and C83 results in increased activity of GLRX and resistance to oxidative inactivation and aggregation. Results from these integrated computational and experimental studies provide insights into the relative oxidizability of GLRX's cysteines and have implications for the use of GLRX as a therapeutic in settings of dysregulated protein glutathionylation.

Novel Lennard-Jones Parameters for Cysteine and Selenocysteine in the AMBER Force Field

Cysteine is a common amino acid with a thiol group that plays a pivotal role in a variety of scenarios in redox biochemistry. In contrast, selenocysteine, the 21st amino acid, is only present in 25 human proteins. Classical force-field parameters for cysteine and selenocysteine are still scarce. In this context, we present a methodology to obtain Lennard-Jones parameters for cysteine and selenocysteine in different physiologically relevant oxidation and protonation states. The new force field parameters obtained in this work are available at https://github.com/MALBECC/AMBER-parameters-database. The parameters were adjusted to reproduce water radial distribution functions obtained by density functional theory ab initio molecular dynamics. We validated the results by evaluating the impact of the choice of parameters on the structure and dynamics in classical molecular dynamics simulations of representative proteins containing catalytic cysteine/selenocysteine residues. There are significant changes in protein structure and dynamics depending on the parameters choice, specifically affecting the residues close to the catalytic sites.

Peroxiredoxins wear many hats: Factors that fashion their peroxide sensing personalities

Peroxiredoxins (Prdxs) sense and assess peroxide levels, and signal through protein interactions. Understanding the role of the multiple structural and post-translational modification (PTM) layers that tunes the peroxiredoxin specificities is still a challenge. In this review, we give a tabulated overview on what is known about human and bacterial peroxiredoxins with a focus on structure, PTMs, and protein-protein interactions. Armed with numerous cellular and atomic level experimental techniques, we look at the future and ask ourselves what is still needed to give us a clearer view on the cellular operating power of Prdxs in both stress and non-stress conditions.

Sulfenic acid in human serum albumin: Reaction with thiols, oxidation and spontaneous decay

Human serum albumin (HSA) contains 17 disulfides and only one reduced cysteine, Cys34, which can be oxidized to a relatively stable sulfenic acid (HSA-SOH). This derivative has been previously detected and quantified. However, its properties are poorly understood. Herein, HSA-SOH formation from the exposure of HSA to hydrogen peroxide was confirmed using the sulfenic acid probe bicyclo [6.1.0]nonyne-biotin (BCN-Bio1), and by direct detection by whole protein mass spectrometry. The decay pathways of HSA-SOH were studied. HSA-SOH reacted with a thiol leading to the formation of a mixed disulfide. The reaction occurred through a concerted or direct displacement mechanism (S_N2) with the thiolate (RS^-) as nucleophile towards HSA-SOH. The net charge of the thiolate affected the value of the rate constant. In the presence of hydrogen peroxide, HSA-SOH was further oxidized to sulfinic acid (HSA-SO₂H) and sulfonic acid (HSA-SO₃H). The rate constants of these reactions were estimated. Lastly, HSA-SOH spontaneously decayed in solution. Mass spectrometry experiments suggested that the decay product is a sulfenylamide (HSA-SN(R')R″). Chromatofocusing analysis showed that the overoxidation with hydrogen peroxide predominates at alkaline pH whereas the spontaneous decay predominates at acidic pH. The present findings provide insights into the reactivity and fate of the sulfenic acid in albumin, which are also of relevance to numerous sulfenic acid-mediated processes in redox biology and catalysis.

Acidity and nucleophilic reactivity of glutathione persulfide

Persulfides (RSSH/RSS^-) participate in sulfur trafficking and metabolic processes, and are proposed to mediate the signaling effects of hydrogen sulfide (H₂S). Despite their growing relevance, their chemical properties are poorly understood. Herein, we studied experimentally and computationally the formation, acidity, and nucleophilicity of glutathione persulfide (GSSH/GSS^-), the derivative of the abundant cellular thiol glutathione (GSH). We characterized the kinetics and equilibrium of GSSH formation from glutathione disulfide and H₂S. A pK_a of 5.45 for GSSH was determined, which is 3.49 units below that of GSH. The reactions of GSSH with the physiologically relevant electrophiles peroxynitrite and hydrogen peroxide, and with the probe monobromobimane, were studied and compared with those of thiols. These reactions occurred through S_N2 mechanisms. At neutral pH, GSSH reacted faster than GSH because of increased availability of the anion and, depending on the electrophile, increased reactivity. In addition, GSS^- presented higher nucleophilicity with respect to a thiolate with similar basicity. This can be interpreted in terms of the so-called α effect, i.e. the increased reactivity of a nucleophile when the atom adjacent to the nucleophilic atom has high electron density. The magnitude of the α effect correlated with the Brønsted nucleophilic factor, β_nuc, for the reactions with thiolates and with the ability of the leaving group. Our study constitutes the first determination of the pK_a of a biological persulfide and the first examination of the α effect in sulfur nucleophiles, and sheds light on the chemical basis of the biological properties of persulfides.