Search

Magnesium (Mg 2+ ), the second most abundant intracellular cation, plays a crucial role in cellular functions. In this study, we investigate the interaction between Mg 2+ and coenzyme A (CoA), a thiol-containing cofactor central to cellular metabolism also involved in protein modifications. Isothermal titration calorimetry revealed a 1:1 binding stoichiometry between Mg 2+ and free CoA under biologically relevant conditions. Association constants of (537 ± 20) M -1 and (312 ± 7) M -1 were determined at 25 °C and pH 7.2 and 7.8, respectively, suggesting that a significant fraction of CoA is likely bound to Mg 2+ both in the cytosol and in the mitochondrial matrix. Additionally, the process is entropically-driven, and our results support that the origin of the entropy gain is solvent-related. On the other hand, the combination of 1- and 2-dimensional nuclear magnetic resonance spectroscopy with molecular dynamics simulations and unsupervised learning demonstrate a direct coordination between Mg 2+ and the phosphate groups of the 4-phosphopantothenate unit and bound to position 5' of the adenosine ring. Interestingly, the phosphate in position 3' only indirectly contributes to Mg 2+ coordination. Finally, we discuss how the binding of Mg 2+ to CoA perturbates the chemical environment of different CoA atoms, regardless of their apparent proximity to the coordination site, through the modulation of the CoA conformational landscape. This insight holds implications for understanding the impact on both CoA and Mg 2+ functions in physiological and pathological processes., (Copyright © 2024 Elsevier Inc. All rights reserved.)

We present a comprehensive theoretical examination of the structural properties of dianionic polysulfides [S n ] 2- ( n = 2-6), their conjugated monoacids [HS n ] - ( n = 2-6), and a selection of 1e - -oxidized radical anions [S n ] •- ( 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 [S 4 ] •- 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.

Machine learning (ML) methods have reached high accuracy levels for the prediction of in vacuo molecular properties. However, the simulation of large systems solely through ML methods (such as those based on neural network potentials) is still a challenge. In this context, one of the most promising frameworks for integrating ML schemes in the simulation of complex molecular systems are the so-called ML/MM methods. These multiscale approaches combine ML methods with classical force fields (MM), in the same spirit as the successful hybrid quantum mechanics-molecular mechanics methods (QM/MM). The key issue for such ML/MM methods is an adequate description of the coupling between the region of the system described by ML and the region described at the MM level. In the context of QM/MM schemes, the main ingredient of the interaction is electrostatic, and the state of the art is the so-called electrostatic-embedding. In this study, we analyze the quality of simpler mechanical embedding-based approaches, specifically focusing on their application within a ML/MM framework utilizing atomic partial charges derived in vacuo. Taking as reference electrostatic embedding calculations performed at a QM(DFT)/MM level, we explore different atomic charges schemes, as well as a polarization correction computed using atomic polarizabilites. Our benchmark data set comprises a set of about 80k small organic structures from the ANI-1x and ANI-2x databases, solvated in water. The results suggest that the minimal basis iterative stockholder (MBIS) atomic charges yield the best agreement with the reference coupling energy. Remarkable enhancements are achieved by including a simple polarization correction.

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 ( Mt msrA and Mt msrB) 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 Mt msrA 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 ( k cat / K M = 2656 ± 525 M -1 s -1 and 1.7 ± 0.8 × 10 6 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.

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics., (© 2023. The Author(s).)

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Inc.)

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.

The mechanism of the metal centered reduction of metmyoglobin (MbFe III ) by sulfide species (H 2 S/HS - ) under an argon atmosphere has been studied by a combination of spectroscopic, kinetic, and computational methods. Asymmetric S-shaped time-traces for the formation of MbFe II at varying ratios of excess sulfide were observed at pH 5.3 < pH < 8.0 and 25 °C, suggesting an autocatalytic reaction mechanism. An increased rate at more alkaline pHs points to HS - as relevant reactive species for the reduction. The formation of the sulfanyl radical (HS • ) in the slow initial phase was assessed using the spin-trap phenyl N - tert -butyl nitrone. This radical initiates the formation of S-S reactive species as disulfanuidyl/ disulfanudi-idyl radical anions and disulfide (HSSH •- /HSS •2- and HSS - , respectively). The autocatalysis has been ascribed to HSS - , formed after HSSH •- /HSS •2- disproportionation, which behaves as a fast reductant toward the intermediate complex MbFe III (HS - ). We propose a reaction mechanism for the sulfide-mediated reduction of metmyoglobin where only ferric heme iron initiates the oxidation of sulfide species. Beside the chemical interest, this insight into the MbFe III /sulfide reaction under an argon atmosphere is relevant for the interpretation of biochemical aspects of ectopic myoglobins found on hypoxic tissues toward reactive sulfur species.

Coenzyme A (CoA) is a key cellular metabolite which participates in diverse metabolic pathways, regulation of gene expression and the antioxidant defense mechanism. Human NME1 (hNME1), which is a moonlighting protein, was identified as a major CoA-binding protein. Biochemical studies showed that hNME1 is regulated by CoA through both covalent and non-covalent binding, which leads to a decrease in the hNME1 nucleoside diphosphate kinase (NDPK) activity. In this study, we expanded the knowledge on previous findings by focusing on the non-covalent mode of CoA binding to the hNME1. With X-ray crystallography, we solved the CoA bound structure of hNME1 (hNME1-CoA) and determined the stabilization interactions CoA forms within the nucleotide-binding site of hNME1. A hydrophobic patch stabilizing the CoA adenine ring, while salt bridges and hydrogen bonds stabilizing the phosphate groups of CoA were observed. With molecular dynamics studies, we extended our structural analysis by characterizing the hNME1-CoA structure and elucidating possible orientations of the pantetheine tail, which is absent in the X-ray structure due to its flexibility. Crystallographic studies suggested the involvement of arginine 58 and threonine 94 in mediating specific interactions with CoA. Site-directed mutagenesis and CoA-based affinity purifications showed that arginine 58 mutation to glutamate (R58E) and threonine 94 mutation to aspartate (T94D) prevent hNME1 from binding to CoA. Overall, our results reveal a unique mode by which hNME1 binds CoA, which differs significantly from that of ADP binding: the α- and β-phosphates of CoA are oriented away from the nucleotide-binding site, while 3'-phosphate faces catalytic histidine 118 (H118). The interactions formed by the CoA adenine ring and phosphate groups contribute to the specific mode of CoA binding to hNME1.

Persulfides (RSSH/RSS - ) are species closely related to thiols (RSH/RS - ) and hydrogen sulfide (H 2 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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Cuevasanta, Benchoam, Semelak, Möller, Zeida, Trujillo, Alvarez and Estrin.)

Coenzyme A (CoA) is a key cellular metabolite known for its diverse functions in metabolism and regulation of gene expression. CoA was recently shown to play an important antioxidant role under various cellular stress conditions by forming a disulfide bond with proteins, termed CoAlation. Using anti-CoA antibodies and liquid chromatography tandem mass spectrometry (LC-MS/MS) methodologies, CoAlated proteins were identified from various organisms/tissues/cell-lines under stress conditions. In this study, we integrated currently known CoAlated proteins into mammalian and bacterial datasets (CoAlomes), resulting in a total of 2093 CoAlated proteins (2862 CoAlation sites). Functional classification of these proteins showed that CoAlation is widespread among proteins involved in cellular metabolism, stress response and protein synthesis. Using 35 published CoAlated protein structures, we studied the stabilization interactions of each CoA segment (adenosine diphosphate (ADP) moiety and pantetheine tail) within the microenvironment of the modified cysteines. Alternating polar-non-polar residues, positively charged residues and hydrophobic interactions mainly stabilize the pantetheine tail, phosphate groups and the ADP moiety, respectively. A flexible nature of CoA is observed in examined structures, allowing it to adapt its conformation through interactions with residues surrounding the CoAlation site. Based on these findings, we propose three modes of CoA binding to proteins. Overall, this study summarizes currently available knowledge on CoAlated proteins, their functional distribution and CoA-protein stabilization interactions.

Persulfides (RSSH/RSS - ) participate in sulfur trafficking and metabolic processes, and are proposed to mediate the signaling effects of hydrogen sulfide (H 2 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 2 S. A p K 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 N 2 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 p K 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., Competing Interests: Conflict of interest—The authors declare that they have no conflicts of interest with the contents of this article., (© 2020 Benchoam et al.)

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 ( Mt AhpE). 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 Mt AhpE are not relevant in the overoxidation process. Indeed, the calculations propose that the process is unfavored by this particular enzyme microenvironment.

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.