Assessment of Embedding Schemes in a Hybrid Machine Learning/Classical Potentials (ML/MM) Approach.

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.

Inorganic Polysulfides in Solution: Structural Properties and Conformational Isomerism.

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

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.

The carbonyl-lock mechanism underlying non-aromatic fluorescence in biological matter.

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.

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.

Autocatalytic Mechanism in the Anaerobic Reduction of Metmyoglobin by Sulfide Species.

The mechanism of the metal centered reduction of metmyoglobin (MbFeIII) by sulfide species (H2S/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 MbFeII 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 MbFeIII(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 MbFeIII/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.

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

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

A Unique Mode of Coenzyme A Binding to the Nucleotide Binding Pocket of Human Metastasis Suppressor NME1.

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.

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 (H2S/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.

Profiling the Site of Protein CoAlation and Coenzyme A Stabilization Interactions.

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.

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.

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.

Theoretical investigation of the mechanism of nitroxyl decomposition in aqueous solution.

Nitroxyl (HNO) is a species that has been proposed recently to play different roles in nitrosative stress processes. HNO decomposition in aqueous solution leading to N2O is a fast reaction that competes with many biochemical reactions in which HNO may be involved. Since molecular determinants of this reaction are still not fully understood, we present in this work an exhaustive analysis of the mechanism in terms of electronic-structure calculations as well as state of the art hybrid quantum mechanics/molecular mechanics molecular dynamics simulations. We characterized the reaction mechanism and computed free energy profiles for the reaction steps using an umbrella sampling procedure. We propose a first dimerization step followed by an acid-base equilibria. Afterwards, the product is formed from two main pathways involving cis-hyponitrous acid (cis-HONNOH) and its conjugate basis as intermediate. Our calculations show preference for the anionic pathway under physiological conditions and allow us to rationalize the results in terms of a molecular description of specific interactions with the solvent. These interactions turn out to be determinant in the stabilization of transition states and, thereby, modifying the free energy barriers. We predict a strong pH-dependence of the overall kinetics of N2O formation, related with the fraction of reactive species available in solution. Finally, we suggest experimental procedures which could validate this mechanism.