The superoxide radical switch in the biology of nitric oxide and peroxynitrite.

The free radical nitric oxide (·NO) is a key mediator in different physiological processes such as vasodilation, neurotransmission, inflammation, and cellular immune responses, and thus preserving its bioavailability is essential. In several disease conditions, superoxide radical (O2·-) production increases and leads to the rapid "inactivation" of ·NO by a diffusion-controlled radical termination reaction that yields a potent and short-lived oxidant, peroxynitrite. This reaction not only limits ·NO bioavailability for physiological signal transduction but also can divert and switch the biochemistry of ·NO toward nitrooxidative processes. Indeed, since the early 1990s peroxynitrite (and its secondary derived species) has been linked to the establishment and progression of different acute and chronic human diseases and also to the normal aging process. Here, we revisit an earlier and classical review on the role of peroxynitrite in human physiology and pathology (Pacher P, Beckman J, Liaudet L. Physiol Rev 87: 315-424, 2007) and further integrate, update, and interpret the accumulated evidence over 30 years of research. Innovative tools and approaches for the detection, quantitation, and sub- or extracellular mapping of peroxynitrite and its secondary products (e.g., protein 3-nitrotyrosine) have allowed us to unambiguously connect the complex biochemistry of peroxynitrite with numerous biological outcomes at the physiological and pathological levels. Furthermore, our current knowledge of the ·NO/O2·- and peroxynitrite interplay at the cell, tissue, and organ levels is assisting in the discovery of therapeutic interventions for a variety of human diseases.

Thiol-based chemical probes exhibit antiviral activity against SARS-CoV-2 via allosteric disulfide disruption in the spike glycoprotein.

The development of small-molecules targeting different components of SARS-CoV-2 is a key strategy to complement antibody-based treatments and vaccination campaigns in managing the COVID-19 pandemic. Here, we show that two thiol-based chemical probes that act as reducing agents, P2119 and P2165, inhibit infection by human coronaviruses, including SARS-CoV-2, and decrease the binding of spike glycoprotein to its receptor, the angiotensin-converting enzyme 2 (ACE2). Proteomics and reactive cysteine profiling link the antiviral activity to the reduction of key disulfides, specifically by disruption of the Cys379-Cys432 and Cys391-Cys525 pairs distal to the receptor binding motif in the receptor binding domain (RBD) of the spike glycoprotein. Computational analyses provide insight into conformation changes that occur when these disulfides break or form, consistent with an allosteric role, and indicate that P2119/P2165 target a conserved hydrophobic binding pocket in the RBD with the benzyl thiol-reducing moiety pointed directly toward Cys432. These collective findings establish the vulnerability of human coronaviruses to thiol-based chemical probes and lay the groundwork for developing compounds of this class, as a strategy to inhibit the SARS-CoV-2 infection by shifting the spike glycoprotein redox scaffold.

De novo sequencing and construction of a unique antibody for the recognition of alternative conformations of cytochrome c in cells.

Cytochrome c (cyt c) can undergo reversible conformational changes under biologically relevant conditions. Revealing these alternative cyt c conformers at the cell and tissue level is challenging. A monoclonal antibody (mAb) identifying a key conformational change in cyt c was previously reported, but the hybridoma was rendered nonviable. To resurrect the mAb in a recombinant form, the amino-acid sequences of the heavy and light chains were determined by peptide mapping-mass spectrometry-bioinformatic analysis and used to construct plasmids encoding the full-length chains. The recombinant mAb (R1D3) was shown to perform similarly to the original mAb in antigen-binding assays. The mAb bound to a variety of oxidatively modified cyt c species (e.g., nitrated at Tyr74 or oxidized at Met80), which lose the sixth heme ligation (Fe-Met80); it did not bind to several cyt c phospho- and acetyl-mimetics. Peptide competition assays together with molecular dynamic studies support that R1D3 binds a neoepitope within the loop 40-57. R1D3 was employed to identify alternative conformations of cyt c in cells under oxidant- or senescence-induced challenge as confirmed by immunocytochemistry and immunoaffinity studies. Alternative conformers translocated to the nuclei without causing apoptosis, an observation that was further confirmed after pinocytic loading of oxidatively modified cyt c to B16-F1 cells. Thus, alternative cyt c conformers, known to gain peroxidatic function, may represent redox messengers at the cell nuclei. The availability and properties of R1D3 open avenues of interrogation regarding the presence and biological functions of alternative conformations of cyt c in mammalian cells and tissues.

Catalytic Mechanism of Mycobacterium tuberculosis Methionine Sulfoxide Reductase A.

The oxidation of Met to methionine sulfoxide (MetSO) by oxidants such as hydrogen peroxide, hypochlorite, or peroxynitrite has profound effects on protein function. This modification can be reversed by methionine sulfoxide reductases (msr). In the context of pathogen infection, the reduction of oxidized proteins gains significance due to microbial oxidative damage generated by the immune system. For example, Mycobacterium tuberculosis (Mt) utilizes msrs (MtmsrA and MtmsrB) as part of the repair response to the host-induced oxidative stress. The absence of these enzymes makes Mycobacteria prone to increased susceptibility to cell death, pointing them out as potential therapeutic targets. This study provides a detailed characterization of the catalytic mechanism of MtmsrA using a comprehensive approach, including experimental techniques and theoretical methodologies. Confirming a ping-pong type enzymatic mechanism, we elucidate the catalytic parameters for sulfoxide and thioredoxin substrates (kcat/KM = 2656 ± 525 M-1 s-1 and 1.7 ± 0.8 × 106 M-1 s-1, respectively). Notably, the entropic nature of the activation process thermodynamics, representing ∼85% of the activation free energy at room temperature, is underscored. Furthermore, the current study questions the plausibility of a sulfurane intermediate, which may be a transition-state-like structure, suggesting the involvement of a conserved histidine residue as an acid-base catalyst in the MetSO reduction mechanism. This mechanistic insight not only advances our understanding of Mt antioxidant enzymes but also holds implications for future drug discovery and biotechnological applications.

Revisiting the role of 3-nitrotyrosine residues in the formation of alpha-synuclein oligomers and fibrils.

Nitration of tyrosine residues in alpha-synuclein (a-syn) has been detected in different synucleinopathies, including Parkinson's disease. The potential role of 3-nitrotyrosine formation in a-syn, as an oxidative post-translational modification, is still elusive. In this work, we generated well-characterized tyrosine nitrated a-syn monomers and studied their capability to form oligomers and fibrils. We constructed tyrosine to phenylalanine mutants, containing a single tyrosine residue, a-syn mutant Y(125/133/136)F and Y(39/125/133)F) and assessed the impact in a-syn biophysical properties. Nitrated wild-type a-syn and the Y-F mutants, with one 3-nitrotyrosine residue in either the protein's N-terminal or C-terminal region, showed inhibition of fibril formation but retained the capacity of oligomer formation. The inhibition of a-syn fibrillation occurs even when an important amount of unmodified a-syn is still present. We characterized oligomers from both nitrated and non-nitrated forms of the wild-type protein and the mutant forms obtained. Our results indicate that the formation of 3-nitrotyrosine in a-syn could induce an off-pathway oligomer formation which may have an important impact in the development of synucleinopathies.

Crystal structure of Trypanosoma cruzi heme peroxidase and characterization of its substrate specificity and compound I intermediate.

The protozoan parasite Trypanosoma cruzi is the causative agent of American trypanosomiasis, otherwise known as Chagas disease. To survive in the host, the T. cruzi parasite needs antioxidant defense systems. One of these is a hybrid heme peroxidase, the T. cruzi ascorbate peroxidase-cytochrome c peroxidase enzyme (TcAPx-CcP). TcAPx-CcP has high sequence identity to members of the class I peroxidase family, notably ascorbate peroxidase (APX) and cytochrome c peroxidase (CcP), as well as a mitochondrial peroxidase from Leishmania major (LmP). The aim of this work was to solve the structure and examine the reactivity of the TcAPx-CcP enzyme. Low temperature electron paramagnetic resonance spectra support the formation of an exchange-coupled [Fe(IV)=O Trp233•+] compound I radical species, analogous to that used in CcP and LmP. We demonstrate that TcAPx-CcP is similar in overall structure to APX and CcP, but there are differences in the substrate-binding regions. Furthermore, the electron transfer pathway from cytochrome c to the heme in CcP and LmP is preserved in the TcAPx-CcP structure. Integration of steady state kinetic experiments, molecular dynamic simulations, and bioinformatic analyses indicates that TcAPx-CcP preferentially oxidizes cytochrome c but is still competent for oxidization of ascorbate. The results reveal that TcAPx-CcP is a credible cytochrome c peroxidase, which can also bind and use ascorbate in host cells, where concentrations are in the millimolar range. Thus, kinetically and functionally TcAPx-CcP can be considered a hybrid peroxidase.

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.

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 (H2S). 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 H2S. A pKa 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 SN2 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 pKa 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.

Cardiolipin interactions with cytochrome c increase tyrosine nitration yields and site-specificity.

The interaction between cytochrome c and cardiolipin is a relevant process in the mitochondrial redox homeostasis, playing roles in the mechanism of electron transfer to cytochrome c oxidase and also modulating cytochrome c conformation, reactivity and function. Peroxynitrite is a widespread nitrating agent formed in mitochondria under oxidative stress conditions, and can result in the formation of tyrosine nitrated cytochrome c. Some of the nitro-cytochrome c species undergo conformational changes at physiological pH and increase its peroxidase activity. In this work we evaluated the influence of cardiolipin on peroxynitrite-mediated cytochrome c nitration yields and site-specificity. Our results show that cardiolipin enhances cytochrome c nitration by peroxynitrite and targets it to heme-adjacent Tyr67. Cytochrome c nitration also modifies the affinity of protein with cardiolipin. Using a combination of experimental techniques and computer modeling, it is concluded that structural modifications in the Tyr67 region are responsible for the observed changes in protein-derived radical and tyrosine nitration levels, distribution of nitrated proteoforms and affinity to cardiolipin. Increased nitration of cytochrome c in presence of cardiolipin within mitochondria and the gain of peroxidatic activity could then impact events such as the onset of apoptosis and other processes related to the disruption of mitochondrial redox homeostasis.

Catalysis of Peroxide Reduction by Fast Reacting Protein Thiols.

Life on Earth evolved in the presence of hydrogen peroxide, and other peroxides also emerged before and with the rise of aerobic metabolism. They were considered only as toxic byproducts for many years. Nowadays, peroxides are also regarded as metabolic products that play essential physiological cellular roles. Organisms have developed efficient mechanisms to metabolize peroxides, mostly based on two kinds of redox chemistry, catalases/peroxidases that depend on the heme prosthetic group to afford peroxide reduction and thiol-based peroxidases that support their redox activities on specialized fast reacting cysteine/selenocysteine (Cys/Sec) residues. Among the last group, glutathione peroxidases (GPxs) and peroxiredoxins (Prxs) are the most widespread and abundant families, and they are the leitmotif of this review. After presenting the properties and roles of different peroxides in biology, we discuss the chemical mechanisms of peroxide reduction by low molecular weight thiols, Prxs, GPxs, and other thiol-based peroxidases. Special attention is paid to the catalytic properties of Prxs and also to the importance and comparative outlook of the properties of Sec and its role in GPxs. To finish, we describe and discuss the current views on the activities of thiol-based peroxidases in peroxide-mediated redox signaling processes.

Kinetics of formation and reactivity of the persulfide in the one-cysteine peroxiredoxin from Mycobacterium tuberculosis.

Hydrogen sulfide (H2S) participates in prokaryotic metabolism and is associated with several physiological functions in mammals. H2S reacts with oxidized thiol derivatives (i.e. disulfides and sulfenic acids) and thereby forms persulfides, which are plausible transducers of the H2S-mediated signaling effects. The one-cysteine peroxiredoxin alkyl hydroperoxide reductase E from Mycobacterium tuberculosis (MtAhpE-SH) reacts fast with hydroperoxides, forming a stable sulfenic acid (MtAhpE-SOH), which we chose here as a model to study the interactions between H2S and peroxiredoxins (Prx). MtAhpE-SOH reacted with H2S, forming a persulfide (MtAhpE-SSH) detectable by mass spectrometry. The rate constant for this reaction was (1.4 ± 0.2) × 103 m-1 s-1 (pH 7.4, 25 °C), six times higher than that reported for the reaction with the main low-molecular-weight thiol in M. tuberculosis, mycothiol. H2S was able to complete the catalytic cycle of MtAhpE and, according to kinetic considerations, it could represent an alternative substrate in M. tuberculosis. MtAhpE-SSH reacted 43 times faster than did MtAhpE-SH with the unspecific electrophile 4,4'-dithiodipyridine, a disulfide that exhibits no preferential reactivity with peroxidatic cysteines, but MtAhpE-SSH was less reactive toward specific Prx substrates such as hydrogen peroxide and peroxynitrite. According to molecular dynamics simulations, this loss of specific reactivity could be explained by alterations in the MtAhpE active site. MtAhpE-SSH could transfer its sulfane sulfur to a low-molecular-weight thiol, a process likely facilitated by the low pKa of the leaving thiol MtAhpE-SH, highlighting the possibility that Prx participates in transpersulfidation. The findings of our study contribute to the understanding of persulfide formation and reactivity.

Genetic and transcriptomic evidences suggest ARO10 genes are involved in benzenoid biosynthesis by yeast.

Benzenoids are compounds associated with floral and fruity flavours in flowers, fruits and leaves and present a role in hormonal signalling in plants. These molecules are produced by the phenyl ammonia lyase pathway. However, some yeasts can also synthesize them from aromatic amino acids using an alternative pathway that remains unknown. Hanseniaspora vineae can produce benzenoids at levels up to two orders of magnitude higher than Saccharomyces species, so it is a model microorganism for studying benzenoid biosynthesis pathways in yeast. According to their genomes, several enzymes have been proposed to be involved in a mandelate pathway similar to that described for some prokaryotic cells. Among them, the ARO10 gene product could present benzoylformate decarboxylase activity. This enzyme catalyses the decarboxylation of benzoylformate into benzaldehyde at the end of the mandelate pathway in benzyl alcohol formation. Two homologous genes of ARO10 were found in the two sequenced H. vineae strains. In this study, nine other H. vineae strains were analysed to detect the presence and per cent homology of ARO10 sequences by PCR using specific primers designed for this species. Also, the copy number of the genes was estimated by quantitative PCR. To verify the relation of ARO10 with the production of benzyl alcohol during fermentation, a deletion mutant in the ARO10 gene of Saccharomyces cerevisiae was used. The two HvARO10 paralogues were analysed and compared with other α-ketoacid decarboxylases at the sequence and structural level.

Exploring the conformational transition between the fully folded and locally unfolded substates of Escherichia coli thiol peroxidase.

Thiol peroxidase from Escherichia coli (EcTPx) is a peroxiredoxin that catalyzes the reduction of different hydroperoxides. During the catalytic cycle of EcTPx, the peroxidatic cysteine (CP) is oxidized to a sulfenic acid by peroxide, then the resolving cysteine (CR) condenses with the sulfenic acid of CP to form a disulfide bond, which is finally reduced by thioredoxin. Purified EcTPx as dithiol and disulfide behaves as a monomer under near physiological conditions. Although secondary structure rearrangements are present when comparing different redox states of the enzyme, no significant differences in unfolding free energies are observed under reducing and oxidizing conditions. A conformational change denominated fully folded (FF) to locally unfolded (LU) transition, involving a partial unfolding of αH2 and αH3, must occur to enable the formation of the disulfide bond since the catalytic cysteines are 12 Å apart in the FF conformation of EcTPx. To explore this process, the FF → LU and LU → FF transitions were studied using conventional molecular dynamics simulations and an enhanced conformational sampling technique for different oxidation and protonation states of the active site cysteine residues CP and CR. Our results suggest that the FF → LU transition has a higher associated energy barrier than the refolding LU → FF process in agreement with the relatively low experimental turnover number of EcTPx. Furthermore, in silico designed single-point mutants of αH3 enhanced locally unfolding events, suggesting that the native FF interactions in the active site are not evolutionarily optimized to fully speed-up the conformational transition of wild-type EcTPx.

Kinetics, subcellular localization, and contribution to parasite virulence of a Trypanosoma cruzi hybrid type A heme peroxidase (TcAPx-CcP).

The Trypanosoma cruzi ascorbate peroxidase is, by sequence analysis, a hybrid type A member of class I heme peroxidases [TcAPx-cytochrome c peroxidase (CcP)], suggesting both ascorbate (Asc) and cytochrome c (Cc) peroxidase activity. Here, we show that the enzyme reacts fast with H2O2 (k = 2.9 × 107 M-1⋅s-1) and catalytically decomposes H2O2 using Cc as the reducing substrate with higher efficiency than Asc (kcat/Km = 2.1 × 105 versus 3.5 × 104 M-1⋅s-1, respectively). Visible-absorption spectra of purified recombinant TcAPx-CcP after H2O2 reaction denote the formation of a compound I-like product, characteristic of the generation of a tryptophanyl radical-cation (Trp233•+). Mutation of Trp233 to phenylalanine (W233F) completely abolishes the Cc-dependent peroxidase activity. In addition to Trp233•+, a Cys222-derived radical was identified by electron paramagnetic resonance spin trapping, immunospin trapping, and MS analysis after equimolar H2O2 addition, supporting an alternative electron transfer (ET) pathway from the heme. Molecular dynamics studies revealed that ET between Trp233 and Cys222 is possible and likely to participate in the catalytic cycle. Recognizing the ability of TcAPx-CcP to use alternative reducing substrates, we searched for its subcellular localization in the infective parasite stages (intracellular amastigotes and extracellular trypomastigotes). TcAPx-CcP was found closely associated with mitochondrial membranes and, most interestingly, with the outer leaflet of the plasma membrane, suggesting a role at the host-parasite interface. TcAPx-CcP overexpressers were significantly more infective to macrophages and cardiomyocytes, as well as in the mouse model of Chagas disease, supporting the involvement of TcAPx-CcP in pathogen virulence as part of the parasite antioxidant armamentarium.

Exploring the Catalytic Mechanism of Human Glutamine Synthetase by Computer Simulations.

Glutamine synthetase is an important enzyme that catalyzes the ATP-dependent formation of glutamine from glutamate and ammonia. In mammals, it plays a key role in preventing excitotoxicity in the brain and detoxifying ammonia in the liver. In plants and bacteria, it is fundamental for nitrogen metabolism, being critical for the survival of the organism. In this work, we show how the use of classical molecular dynamics simulations and multiscale quantum mechanics/molecular mechanics simulations allowed us to examine the structural properties and dynamics of human glutamine synthetase (HsGS), as well as the reaction mechanisms involved in the catalytic process with atomic level detail. Our results suggest that glutamine formation proceeds through a two-step mechanism that includes a first step in which the γ-glutamyl phosphate intermediate forms, with a 5 kcal/mol free energy barrier and a -8 kcal/mol reaction free energy, and then a second rate-limiting step involving the ammonia nucleophilic attack, with a free energy barrier of 19 kcal/mol and a reaction free energy of almost zero. A detailed analysis of structural features within each step exposed the relevance of the acid-base equilibrium related to protein residues and substrates in the thermodynamics and kinetics of the reactions. These results provide a comprehensive study of HsGS dynamics and establish the groundwork for further analysis regarding changes in HsGS activity, as occur in natural variants and post-translational modifications.

Nitrosodisulfide [S2NO]- (perthionitrite) is a true intermediate during the "cross-talk" of nitrosyl and sulfide.

Nitrosodisulfide S2NO- is a controversial intermediate in the reactions of S-nitrosothiols with HS- that produce NO and HNO. QM-MM molecular dynamics simulations combined with TD-DFT analysis contribute to a clear identification of S2NO- in water, acetone and acetonitrile, accounting for the UV-Vis signatures and broadening the mechanistic picture of N/S signaling in biochemistry.

Molecular Basis of Hydroperoxide Specificity in Peroxiredoxins: The Case of AhpE from Mycobacterium tuberculosis.

Peroxiredoxins (Prxs) constitute a ubiquitous family of Cys-dependent peroxidases that play essential roles in reducing hydrogen peroxide, peroxynitrite, and organic hydroperoxides in almost all organisms. Members of the Prx subfamilies show differential oxidizing substrate specificities that await explanations at a molecular level. Among them, alkyl hydroperoxide reductases E (AhpE) is a novel subfamily comprising Mycobacterium tuberculosis AhpE and AhpE-like proteins expressed in some bacteria and archaea. We previously reported that MtAhpE reacts ∼10(4) times faster with an arachidonic acid derived hydroperoxide than with hydrogen peroxide, and suggested that this surprisingly high reactivity was related to the presence of a hydrophobic groove at the dimer interface evidenced in the crystallography structure of the enzyme. In this contribution we experimentally confirmed the existence of an exposed hydrophobic patch in MtAhpE. We found that fatty acid hydroperoxide reduction by the enzyme showed positive activation entropy that importantly contributed to catalysis. Computational dynamics indicated that interactions of fatty acid-derived hydroperoxides with the enzyme properly accommodated them inside the active site and modifies enzyme's dynamics. The computed reaction free energy profile obtained via QM/MM simulations is consistent with a greater reactivity in comparison with hydrogen peroxide. This study represents new insights on the understanding of the molecular basis that determines oxidizing substrate selectivity in the peroxiredoxin family, which has not been investigated at an atomic level so far.

Kinetic and structural assessment of the reduction of human 2-Cys peroxiredoxins by thioredoxins.

We have studied the reduction reactions of two cytosolic human peroxiredoxins (Prx) in their disulfide form by three thioredoxins (Trx; two human and one bacterial), with the aim of better understanding the rate and mechanism of those reactions, and their relevance in the context of the catalytic cycle of Prx. We have developed a new methodology based on stopped-flow and intrinsic fluorescence to study the bimolecular reactions, and found rate constants in the range of 105 -106 m-1 s-1 in all cases, showing that there is no marked kinetic preference for the expected Trx partner. By combining experimental findings and molecular dynamics studies, we found that the reactivity of the nucleophilic cysteine (CN ) in the Trx is greatly affected by the formation of the Prx-Trx complex. The protein-protein interaction forces the CN thiolate into an unfavorable hydrophobic microenvironment that reduces its hydration and results in a remarkable acceleration of the thiol-disulfide exchange reactions by more than three orders of magnitude and also produces a measurable shift in the pKa of the CN . This mechanism of activation of the thiol disulfide exchange may help understand the reduction of Prx by alternative reductants involved in redox signaling.

Mechanism of cysteine oxidation by peroxynitrite: An integrated experimental and theoretical study.

Since peroxynitrite was identified as a pathophysiological agent it has been implicated in a great variety of cellular processes. Particularly, peroxynitrite mediated oxidation of cellular thiol-containing compounds such as Cys residues, is a key event which has been extensively studied. Although great advances have been accomplished, the reaction is not completely understood at the atomic level. Aiming to shed light on this subject, we present an integrated kinetic and theoretical study of the oxidation of free Cys by peroxynitrite. We determined pH-independent thermodynamic activation parameters, namely those corresponding to the reaction between the reactive species: Cys thiolate and peroxynitrous acid. We found a pH-independent activation energy of 8.2 ± 0.6 kcal/mol. Simulations were performed using state of the art hybrid quantum-classical (QM-MM) molecular dynamics simulations. Our results are consistent with a SN2 mechanism, with Cys sulfenic acid and nitrite anion as products. The activation barrier is mostly due to the alignment of sulfur's thiolate atom with the oxygen atoms of the peroxide, along with the concomitant charge reorganization and important changes in the solvation profile. This work provides an atomic detailed description of the reaction mechanism and a framework to understand the environment effects on peroxynitrite reactivity with protein thiols.