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.

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.

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.

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.

Redox sensitive human mitochondrial aconitase and its interaction with frataxin: In vitro and in silico studies confirm that it takes two to tango

Mitochondrial aconitase (ACO2) has been postulated as a redox sensor in the tricarboxylic acid cycle. Its high sensitivity towards reactive oxygen and nitrogen species is due to its particularly labile [4Fe-4S]²⁺ prosthetic group which yields an inactive [3Fe-4S]⁺ cluster upon oxidation. Moreover, ACO2 was found as a main oxidant target during aging and in pathologies where mitochondrial dysfunction is implied. Herein, we report the expression and characterization of recombinant human ACO2 and its interaction with frataxin (FXN), a protein that participates in the de novo biosynthesis of Fe-S clusters. A high yield of pure ACO2 (≥99%, 22 ± 2 U/mg) was obtained and kinetic parameters for citrate, isocitrate, and cis-aconitate were determined. Superoxide, carbonate radical, peroxynitrite, and hydrogen peroxide reacted with ACO2 with second-order rate constants of 10⁸, 10⁸, 10⁵, and 10² M^-1 s^-1, respectively. Temperature-induced unfolding assessed by tryptophan fluorescence of ACO2 resulted in apparent melting temperatures of 51.1 ± 0.5 and 43.6 ± 0.2 °C for [4Fe-4S]²⁺ and [3Fe-4S]⁺ states of ACO2, sustaining lower thermal stability upon cluster oxidation. Differences in protein dynamics produced by the Fe-S cluster redox state were addressed by molecular dynamics simulations. Reactivation of [3Fe-4S]⁺-ACO2 by FXN was verified by activation assays and direct iron-dependent interaction was confirmed by protein-protein interaction ELISA and fluorescence spectroscopic assays. Multimer modeling and protein-protein docking predicted an ACO2-FXN complex where the metal ion binding region of FXN approaches the [3Fe-4S]⁺ cluster, supporting that FXN is a partner for reactivation of ACO2 upon oxidative cluster inactivation.

Mitochondrial Peroxiredoxin 3 Is Rapidly Oxidized and Hyperoxidized by Fatty Acid Hydroperoxides

Human peroxiredoxin 3 (HsPrx3) is a thiol-based peroxidase responsible for the reduction of most hydrogen peroxide and peroxynitrite formed in mitochondria. Mitochondrial disfunction can lead to membrane lipoperoxidation, resulting in the formation of lipid-bound fatty acid hydroperoxides (_LFA-OOHs) which can be released to become free fatty acid hydroperoxides (_fFA-OOHs). Herein, we report that HsPrx3 is oxidized and hyperoxidized by _fFA-OOHs including those derived from arachidonic acid and eicosapentaenoic acid peroxidation at position 15 with remarkably high rate constants of oxidation (>3.5 × 10⁷ M^-1s^-1) and hyperoxidation (~2 × 10⁷ M^-1s^-1). The endoperoxide-hydroperoxide PGG₂, an intermediate in prostanoid synthesis, oxidized HsPrx3 with a similar rate constant, but was less effective in causing hyperoxidation. Biophysical methodologies suggest that HsPrx3 can bind hydrophobic structures. Indeed, molecular dynamic simulations allowed the identification of a hydrophobic patch near the enzyme active site that can allocate the hydroperoxide group of _fFA-OOHs in close proximity to the thiolate in the peroxidatic cysteine. Simulations performed using available and herein reported kinetic data indicate that HsPrx3 should be considered a main target for mitochondrial _fFA-OOHs. Finally, kinetic simulation analysis support that mitochondrial _fFA-OOHs formation fluxes in the range of nM/s are expected to contribute to HsPrx3 hyperoxidation, a modification that has been detected in vivo under physiological and pathological conditions.

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.

Possible molecular basis of the biochemical effects of cysteine-derived persulfides

Persulfides (RSSH/RSS⁻) are species closely related to thiols (RSH/RS⁻) and hydrogen sulfide (H₂S/HS⁻), and can be formed in biological systems in both low and high molecular weight cysteine-containing compounds. They are key intermediates in catabolic and biosynthetic processes, and have been proposed to participate in the transduction of hydrogen sulfide effects. Persulfides are acidic, more acidic than thiols, and the persulfide anions are expected to be the predominant species at neutral pH. The persulfide anion has high nucleophilicity, due in part to the alpha effect, i.e., the increased reactivity of a nucleophile when the neighboring atom has high electron density. In addition, persulfides have electrophilic character, a property that is absent in both thiols and hydrogen sulfide. In this article, the biochemistry of persulfides is described, and the possible ways in which the formation of a persulfide could impact on the properties of the biomolecule involved are discussed.

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.

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.

PIP aquaporin pH sensing is regulated by the length and charge of the C-terminal region

Plant PIP aquaporins play a central role in controlling plant water status. The current structural model for PIP pH-gating states that the main pH sensor is located in loopD and that all the mobile cytosolic elements participate in a complex interaction network that ensures the closed structure. However, the precise participation of the last part of the C-terminal domain (CT) in PIP pH gating remains unknown. This last part has not been resolved in PIP crystal structures and is a key difference between PIP1 and PIP2 paralogues. Here, by a combined experimental and computational approach, we provide data about the role of CT in pH gating of Beta vulgaris PIP. We demonstrate that the length of CT and the positive charge located among its last residues modulate the pH at which the open/closed transition occurs. We also postulate a molecular-based mechanism for the differential pH sensing in PIP homo- or heterotetramers by performing atomistic molecular dynamics simulations (MDS) on complete models of PIP tetramers. Our findings show that the last part of CT can affect the environment of loopD pH sensors in the closed state. Results presented herein contribute to the understanding of how the characteristics of CT in PIP channels play a crucial role in determining the pH at which water transport through these channels is blocked, highlighting the relevance of the differentially conserved very last residues in PIP1 and PIP2 paralogues.

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.

Akt Is S-Palmitoylated: A New Layer of Regulation for Akt

The protein kinase Akt/PKB participates in a great variety of processes, including translation, cell proliferation and survival, as well as malignant transformation and viral infection. In the last few years, novel Akt posttranslational modifications have been found. However, how these modification patterns affect Akt subcellular localization, target specificity and, in general, function is not thoroughly understood. Here, we postulate and experimentally demonstrate by acyl-biotin exchange (ABE) assay and ³H-palmitate metabolic labeling that Akt is S-palmitoylated, a modification related to protein sorting throughout subcellular membranes. Mutating cysteine 344 into serine blocked Akt S-palmitoylation and diminished its phosphorylation at two key sites, T308 and T450. Particularly, we show that palmitoylation-deficient Akt increases its recruitment to cytoplasmic structures that colocalize with lysosomes, a process stimulated during autophagy. Finally, we found that cysteine 344 in Akt1 is important for proper its function, since Akt1-C344S was unable to support adipocyte cell differentiation in vitro. These results add an unexpected new layer to the already complex Akt molecular code, improving our understanding of cell decision-making mechanisms such as cell survival, differentiation and death.

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.

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

Temporal acquisition of the fully folded conformational substate of the Escherichia coli thiol peroxidase by accelerated molecular dynamics simulations.

Multiscale Modeling of Thiol Overoxidation in Peroxiredoxins by Hydrogen Peroxide

In this work, we employ a multiscale quantum-classical mechanics (QM/MM) scheme to investigate the chemical reactivity of sulfenic acids toward hydrogen peroxide, both in aqueous solution and in the protein environment of the peroxiredoxin alkyl hydroperoxide reductase E from Mycobacterium tuberculosis (MtAhpE). The reaction of oxidation of cysteine with hydrogen peroxides, catalyzed by peroxiredoxins, is usually accelerated several orders of magnitude in comparison with the analogous reaction in solution. The resulting cysteine sulfenic acid is then reduced in other steps of the catalytic cycle, recovering the original thiol. However, under some conditions, the sulfenic acid can react with another equivalent of oxidant to form a sulfinic acid. This process is called overoxidation and has been associated with redox signaling. Herein, we employed a multiscale scheme based on density function theory calculations coupled to the classical AMBER force field, developed in our group, to establish the molecular basis of thiol overoxidation by hydrogen peroxide. Our results suggest that residues that play key catalytic roles in the oxidation of MtAhpE are not relevant in the overoxidation process. Indeed, the calculations propose that the process is unfavored by this particular enzyme microenvironment.

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 (H₂S) participates in prokaryotic metabolism and is associated with several physiological functions in mammals. H₂S reacts with oxidized thiol derivatives (i.e. disulfides and sulfenic acids) and thereby forms persulfides, which are plausible transducers of the H₂S-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 H₂S and peroxiredoxins (Prx). MtAhpE-SOH reacted with H₂S, forming a persulfide (MtAhpE-SSH) detectable by mass spectrometry. The rate constant for this reaction was (1.4 ± 0.2) × 10³ 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. H₂S 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 pK_a 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.

From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field

Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.

Chemical Reactivity and Spectroscopy Explored From QM/MM Molecular Dynamics Simulations Using the LIO Code

In this work we present the current advances in the development and the applications of LIO, a lab-made code designed for density functional theory calculations in graphical processing units (GPU), that can be coupled with different classical molecular dynamics engines. This code has been thoroughly optimized to perform efficient molecular dynamics simulations at the QM/MM DFT level, allowing for an exhaustive sampling of the configurational space. Selected examples are presented for the description of chemical reactivity in terms of free energy profiles, and also for the computation of optical properties, such as vibrational and electronic spectra in solvent and protein environments.

PrxQ B from Mycobacterium tuberculosis is a monomeric, thioredoxin-dependent and highly efficient fatty acid hydroperoxide reductase

Mycobacterium tuberculosis (M. tuberculosis) is the intracellular bacterium responsible for tuberculosis disease (TD). Inside the phagosomes of activated macrophages, M. tuberculosis is exposed to cytotoxic hydroperoxides such as hydrogen peroxide, fatty acid hydroperoxides and peroxynitrite. Thus, the characterization of the bacterial antioxidant systems could facilitate novel drug developments. In this work, we characterized the product of the gene Rv1608c from M. tuberculosis, which according to sequence homology had been annotated as a putative peroxiredoxin of the peroxiredoxin Q subfamily (PrxQ B from M. tuberculosis or MtPrxQ B). The protein has been reported to be essential for M. tuberculosis growth in cholesterol-rich medium. We demonstrated the M. tuberculosis thioredoxin B/C-dependent peroxidase activity of MtPrxQ B, which acted as a two-cysteine peroxiredoxin that could function, although less efficiently, using a one-cysteine mechanism. Through steady-state and competition kinetic analysis, we proved that the net forward rate constant of MtPrxQ B reaction was 3 orders of magnitude faster for fatty acid hydroperoxides than for hydrogen peroxide (3×10⁶vs 6×10³M^-1s^-1, respectively), while the rate constant of peroxynitrite reduction was (0.6-1.4) ×10⁶M^-1s^-1 at pH 7.4. The enzyme lacked activity towards cholesterol hydroperoxides solubilized in sodium deoxycholate. Both thioredoxin B and C rapidly reduced the oxidized form of MtPrxQ B, with rates constants of 0.5×10⁶ and 1×10⁶M^-1s^-1, respectively. Our data indicated that MtPrxQ B is monomeric in solution both under reduced and oxidized states. In spite of the similar hydrodynamic behavior the reduced and oxidized forms of the protein showed important structural differences that were reflected in the protein circular dichroism spectra.

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

Nitrosodisulfide S₂NO^- 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 S₂NO^- in water, acetone and acetonitrile, accounting for the UV-Vis signatures and broadening the mechanistic picture of N/S signaling in biochemistry.

Improving Efficiency in SMD Simulations Through a Hybrid Differential Relaxation Algorithm

The fundamental object for studying a (bio)chemical reaction obtained from simulations is the free energy profile, which can be directly related to experimentally determined properties. Although quite accurate hybrid quantum (DFT based)-classical methods are available, achieving statistically accurate and well converged results at a moderate computational cost is still an open challenge. Here, we present and thoroughly test a hybrid differential relaxation algorithm (HyDRA), which allows faster equilibration of the classical environment during the nonequilibrium steering of a (bio)chemical reaction. We show and discuss why (in the context of Jarzynski's Relationship) this method allows obtaining accurate free energy profiles with smaller number of independent trajectories and/or faster pulling speeds, thus reducing the overall computational cost. Moreover, due to the availability and straightforward implementation of the method, we expect that it will foster theoretical studies of key enzymatic processes.

Insights into the mechanism of the reaction between hydrogen sulfide and peroxynitrite

Hydrogen sulfide and peroxynitrite are endogenously generated molecules that participate in biologically relevant pathways. A revision of the kinetic features of the reaction between peroxynitrite and hydrogen sulfide revealed a complex process. The rate constant of peroxynitrite decay, (6.65 ± 0.08) × 10(3) M(-1) s(-1) in 0.05 M sodium phosphate buffer (pH 7.4, 37°C), was affected by the concentration of buffer. Theoretical modeling suggested that, as in the case of thiols, the reaction is initiated by the nucleophilic attack of HS(-) on the peroxide group of ONOOH by a typical bimolecular nucleophilic substitution, yielding HSOH and NO2(-). In contrast to thiols, the reaction then proceeds to the formation of distinct products that absorb near 408 nm. Experiments in the presence of scavengers and carbon dioxide showed that free radicals are unlikely to be involved in the formation of these products. The results are consistent with product formation involving the reactive intermediate HSSH and its fast reaction with a second peroxynitrite molecule. Mass spectrometry and UV-Vis absorption spectra predictions suggest that at least one of the products is HSNO2 or its isomer HSONO.

The extraordinary catalytic ability of peroxiredoxins: a combined experimental and QM/MM study on the fast thiol oxidation step

Peroxiredoxins (Prxs) catalyze the reduction of peroxides, a process of key relevance in a variety of cellular processes. The first step in the catalytic cycle of all Prxs is the oxidation of a cysteine residue to sulfenic acid, which occurs 10(3)-10(7) times faster than in free cysteine. We present an experimental kinetics and hybrid QM/MM investigation to explore the reaction of Prxs with H2O2 using alkyl hydroperoxide reductase E from Mycobacterium tuberculosis as a Prx model. We report for the first time the thermodynamic activation parameters of H2O2 reduction using Prx, which show that protein significantly lowers the activation enthalpy, with an unfavourable entropic effect, compared to the uncatalyzed reaction. The QM/MM simulations show that the remarkable catalytic effects responsible for the fast H2O2 reduction in Prxs are mainly due to an active-site arrangement, which establishes a complex hydrogen bond network activating both reactive species.

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.

Breathing, bubbling, and bending: DNA flexibility from multimicrosecond simulations

Bending of the seemingly stiff DNA double helix is a fundamental physical process for any living organism. Specialized proteins recognize DNA inducing and stabilizing sharp curvatures of the double helix. However, experimental evidence suggests a high protein-independent flexibility of DNA. On the basis of coarse-grained simulations, we propose that DNA experiences thermally induced kinks associated with the spontaneous formation of internal bubbles. Comparison of the protein-induced DNA curvature calculated from the Protein Data Bank with that sampled by our simulations suggests that thermally induced distortions can account for ~80% of the DNA curvature present in experimentally solved structures.

Molecular basis of the mechanism of thiol oxidation by hydrogen peroxide in aqueous solution: challenging the SN2 paradigm

The oxidation of cellular thiol-containing compounds, such as glutathione and protein Cys residues, is considered to play an important role in many biological processes. Among possible oxidants, hydrogen peroxide (H(2)O(2)) is known to be produced in many cell types as a response to a variety of extracellular stimuli and could work as an intracellular messenger. This reaction has been reported to proceed through a S(N)2 mechanism, but despite its importance, the reaction is not completely understood at the atomic level. In this work, we elucidate the reaction mechanism of thiol oxidation by H(2)O(2) for a model methanethiolate system using state of the art hybrid quantum-classical (QM-MM) molecular dynamics simulations. Our results show that the solvent plays a key role in positioning the reactants, that there is a significant charge redistribution in the first stages of the reaction, and that there is a hydrogen transfer process between H(2)O(2) oxygen atoms that occurs after reaching the transition state. These observations challenge the S(N)2 mechanism hypothesis for this reaction. Specifically, our results indicate that the reaction is driven by a tendency of the slightly charged peroxidatic oxygen to become even more negative in the product via an electrophilic attack on the negative sulfur atom. This is inconsistent with the S(N)2 mechanism, which predicts a protonated sulfenic acid and hydroxyl anion as stable intermediates. These intermediates are not found. Instead, the reaction proceeds directly to unprotonated sulfenic acid and water.

A Coarse Grained Model for Atomic-Detailed DNA Simulations with Explicit Electrostatics

Coarse-grain (CG) techniques allow considerable extension of the accessible size and time scales in simulations of biological systems. Although many CG representations are available for the most common biomacromolecules, very few have been reported for nucleic acids. Here, we present a CG model for molecular dynamics simulations of DNA on the multi-microsecond time scale. Our model maps the complexity of each nucleotide onto six effective superatoms keeping the "chemical sense" of specific Watson-Crick recognition. Molecular interactions are evaluated using a classical Hamiltonian with explicit electrostatics calculated under the framework of the generalized Born approach. This CG representation is able to accurately reproduce experimental structures, breathing dynamics, and conformational transitions from the A to the B form in double helical fragments. The model achieves a good qualitative reproduction of temperature-driven melting and its dependence on size, ionic strength, and sequence specificity. Reconstruction of atomistic models from CG trajectories give remarkable agreement with structural, dynamic, and energetic features obtained from fully atomistic simulation, opening the possibility to acquire nearly atomic detail data from CG trajectories.