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

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.

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.

Mechanism of Peroxynitrite Interaction with Ferric M. tuberculosis Nitrobindin: A Computational Study

Nitrobindins (Nbs) are all-β-barrel heme proteins present along the evolutionary ladder. They display a highly solvent-exposed ferric heme group with the iron atom being coordinated by the proximal His residue and a water molecule at the distal position. Ferric nitrobindins (Nb(III)) play a role in the conversion of toxic peroxynitrite (ONOO-) to harmless nitrate, with the value of the second-order rate constant being similar to those of most heme proteins. The value of the second-order rate constant of Nbs increases as the pH decreases; this suggests that Nb(III) preferentially reacts with peroxynitrous acid (ONOOH), although ONOO- is more nucleophilic. In this work, we shed light on the molecular basis of the ONOO- and ONOOH reactivity of ferric Mycobacterium tuberculosis Nb (Mt-Nb(III)) by dissecting the ligand migration toward the active site, the water molecule release, and the ligand binding process by computer simulations. Classical molecular dynamics simulations were performed by employing a steered molecular dynamics approach and the Jarzynski equality to obtain ligand migration free energy profiles for both ONOO- and ONOOH. Our results indicate that ONOO- and ONOOH migration is almost unhindered, consistent with the exposed metal center of Mt-Nb(III). To further analyze the ligand binding process, we computed potential energy profiles for the displacement of the Fe(III)-coordinated water molecule using a hybrid QM/MM scheme at the DFT level and a nudged elastic band approach. These results indicate that ONOO- exhibits a much larger barrier for ligand displacement than ONOOH, suggesting that water displacement is assisted by protonation of the leaving group by the incoming ONOOH.

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.

Nitric oxide binding to ferrous nitrobindins: A computer simulation investigation

Nitrobindins (Nbs) represent an evolutionary conserved all-β-barrel heme-proteins displaying a highly solvent-exposed heme-Fe(III) atom, coordinated by a proximal His residue. Interestingly, even if the distal side is exposed to the solvent, the value of the second order rate constants for ligand binding to the ferrous derivative is almost one order of magnitude lower than those reported for myoglobins (Mbs). Noteworthy, nitric oxide binding to the sixth coordination position of the heme-Fe(II)-atom causes the cleavage or the severe weakening of the proximal His-Fe(II) bond. Here, we provide a computer simulation investigation to shed light on the molecular basis of ligand binding kinetics, by dissecting the ligand binding process into the ligand migration and the bond formation steps. Classical molecular dynamics simulations were performed employing a steered molecular dynamics approach and the Jarzinski equality to obtain ligand migration free energy profiles. The formation of the heme-Fe(II)-NO bond took into consideration the iron atom displacement from the heme plane. The ligand migration is almost unhindered, and the low rate constant for NO binding is due to the large displacement of the Fe(II) atom with respect to the heme plane responsible for the barrier for the Fe(II)-NO bond formation. In addition, we investigated the weakening and breaking of the proximal His-Fe(II) bond, observed experimentally upon NO binding, by means of a combination of classical molecular dynamics simulations and quantum-classical (QM-MM) optimizations. In both human and M. tuberculosis Nbs, a stable alternative conformation of the proximal His residue interacting with a network of water molecules was observed.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Binding mechanism of disulfide species to ferric hemeproteins: The case of metmyoglobin

The interactions of the heme iron of hemeproteins with sulfide and disulfide compounds are of potential interest as physiological signaling processes. While the interaction with hydrogen sulfide has been described computationally and experimentally, the reaction with disulfide, and specifically the molecular mechanism for ligand binding has not been studied in detail. In this work, we study the association process for disulfane and its conjugate base disulfanide at different pH conditions. Additionally, by means of advanced sampling techniques based on multiple steered molecular dynamics, we provide free energy profiles for ligand migration for both acid/base species, showing a similar behavior to the previously reported for the related H2S/HS¯ pair. Finally, we studied the ligand interchange reaction (H2O/H2S, HS¯ and H2O/HSSH, HSS¯) by means of hybrid quantum mechanics-molecular mechanics calculations. We show that the anionic species are able to displace more efficiently the H2O bound to the iron, and that the H-bond network in the distal cavity can help the neutral species to perform the reaction. Altogether, we provide a molecular explanation for the experimental information and show that the global association process depends on a fine balance between the migration towards the active site and the ligand interchange reaction.

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.

Best Practices on QM/MM Simulations of Biological Systems

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

Tailoring Cooperative Emission in Molecules: Superradiance and Subradiance from First-Principles Simulations

Cooperative optical effects provide a pathway to both the amplification (superradiance) and the suppression (subradiance) of photon emission from electronically excited states. These captivating phenomena offer a rich variety of possibilities for photonic technologies aimed at electromagnetic energy manipulation, including lasers and high-speed emitting devices in the case of superradiance or optical energy storage in that of subradiance. The employment of molecules as the building pieces in these developments requires a precise understanding of the roles of separation, orientation, spatial distribution, and applied fields, which remains challenging for theory and experiments. These questions are addressed here through ab initio quantum dynamics simulations of collective emission on the basis of a novel semiclassical formalism and time-dependent density functional theory. By establishing the configurations leading to decoherence and how the fine-tuning of a pulse can accumulate or release optical energy in H₂ arrays, this report provides fundamental insight toward the design of real superradiant and subradiant devices.

Non-Aromatic Fluorescence in Biological Matter: The Exception or the Rule?

While in the vast majority of cases fluorescence in biological matter has been attributed to aromatic or conjugated groups, peptides associated with neurodegenerative diseases, such as Alzheimer's, Parkinson's, or Huntington's, have been recently shown to display an intrinsic visible fluorescence even in the absence of aromatic residues. This has called the attention of researchers from many different fields, trying to understand the origin of this peculiar behavior and, at the same time, motivating the search for novel strategies to control the optical properties of new biophotonic materials. Today, after nearly 15 years of its discovery, there is a growing consensus about the mechanism underlying this phenomenon, namely, that electronic interactions between non-optically active molecules can result in supramolecular assemblies that are fluorescent. Despite this progress, many aspects of this phenomenon remain uncharted territory. In this Perspective, we lay down the state-of-the-art in the field highlighting the open questions from both experimental and theoretical fronts in this fascinating emerging area of non-aromatic fluorescence.

Mycobacterium tuberculosis DosS binds H2S through its Fe3+ heme iron to regulate the DosR dormancy regulon

Mycobacterium tuberculosis (Mtb) senses and responds to host-derived gasotransmitters NO and CO via heme-containing sensor kinases DosS and DosT and the response regulator DosR. Hydrogen sulfide (H2S) is an important signaling molecule in mammals, but its role in Mtb physiology is unclear. We have previously shown that exogenous H2S can modulate expression of genes in the Dos dormancy regulon via an unknown mechanism(s). Here, we test the hypothesis that Mtb senses and responds to H2S via the DosS/T/R system. Using UV-Vis and EPR spectroscopy, we show that H2S binds directly to the ferric (Fe3+) heme of DosS (KDapp = 5.30 μM) but not the ferrous (Fe2+) form. No interaction with DosT(Fe2+-O2) was detected. We found that the binding of sulfide can slowly reduce the DosS heme iron to the ferrous form. Steered Molecular Dynamics simulations show that H2S, and not the charged HS- species, can enter the DosS heme pocket. We also show that H2S increases DosS autokinase activity and subsequent phosphorylation of DosR, and H2S-mediated increases in Dos regulon gene expression is lost in Mtb lacking DosS. Finally, we demonstrate that physiological levels of H2S in macrophages can induce DosR regulon genes via DosS. Overall, these data reveal a novel mechanism whereby Mtb senses and responds to a third host gasotransmitter, H2S, via DosS(Fe3+). These findings highlight the remarkable plasticity of DosS and establish a new paradigm for how bacteria can sense multiple gasotransmitters through a single heme sensor kinase.

Short hydrogen bonds enhance nonaromatic protein-related fluorescence

Intrinsic fluorescence of nonaromatic amino acids is a puzzling phenomenon with an enormous potential in biophotonic applications. The physical origins of this effect, however, remain elusive. Herein, we demonstrate how specific hydrogen bond networks can modulate fluorescence. We highlight the key role played by short hydrogen bonds, present in the protein structure, on the ensuing fluorescence. We provide detailed experimental and molecular evidence to explain these unusual nonaromatic optical properties. Our findings should benefit the design of novel optically active biomaterials for applications in biosensing and imaging.

Fluorescence in biological systems is usually associated with the presence of aromatic groups. Here, by employing a combined experimental and computational approach, we show that specific hydrogen bond networks can significantly affect fluorescence. In particular, we reveal that the single amino acid L-glutamine, by undergoing a chemical transformation leading to the formation of a short hydrogen bond, displays optical properties that are significantly enhanced compared with L-glutamine itself. Ab initio molecular dynamics simulations highlight that these short hydrogen bonds prevent the appearance of a conical intersection between the excited and the ground states and thereby significantly decrease nonradiative transition probabilities. Our findings open the door to the design of new photoactive materials with biophotonic applications.

Electron transfer pathways from quantum dynamics simulations

This work explores the possibility of simulating an electron transfer process between a donor and an acceptor in real time using time-dependent density functional theory electron dynamics. To achieve this objective, a central issue to resolve is the definition of the initial state. This must be a non-equilibrium electronic state able to trigger the charge transfer dynamics; here, two schemes are proposed to prepare such states. One is based on the combination of the density matrices of the donor and acceptor converged separately with appropriate charges (for example, -1 for the donor and +1 for the acceptor). The second approach relied on constrained DFT to localize the charge on each fragment. With these schemes, electron transfer processes are simulated in different model systems of increasing complexity: an atomic hydrogen dimer, a polyacetylene chain, and the active site of the T. cruzi hybrid type A heme peroxidase, for which two possible electron transfer paths have been postulated. For the latter system, the present methodology applied in a hybrid Quantum Mechanics - Molecular Mechanics framework allows us to establish the relative probabilities of each path and provides insight into the inhibition of the electron transfer provoked by the substitution of tryptophan by phenylalanine in the W233F mutant.

Dissipative Equation of Motion for Electromagnetic Radiation in Quantum Dynamics

The dynamical description of the radiative decay of an electronically excited state in realistic many-particle systems is an unresolved challenge. In the present investigation electromagnetic radiation of the charge density is approximated as the power dissipated by a classical dipole, to cast the emission in closed form as a unitary single-electron theory. This results in a formalism of unprecedented efficiency, critical for ab initio modeling, which exhibits at the same time remarkable properties: it quantitatively predicts decay rates, natural broadening, and absorption intensities. Exquisitely accurate excitation lifetimes are obtained from time-dependent DFT simulations for C^{2+}, B^{+}, and Be, of 0.565, 0.831, and 1.97 ns, respectively, in accord with experimental values of 0.57±0.02, 0.86±0.07, and 1.77-2.5 ns. Hence, the present development expands the frontiers of quantum dynamics, bringing within reach first-principles simulations of a wealth of photophysical phenomena, from fluorescence to time-resolved spectroscopies.

Fluorescence Quantum Yields in Complex Environments from QM-MM TDDFT Simulations: The Case of Indole in Different Solvents

Fluorescence is commonly exploited to probe microscopic properties. An important example is tryptophan in protein environments, where variations in fluorescence quantum yield, and in absorption and emission maxima, are used as indicators of changes in the environment. Modeling the fluorescence quantum yield requires the determination of both radiative and nonradiative decay constants, both on the potential energy surface of the excited fluorophore. Furthermore, the inclusion of complex environments implies their accurate representation as well as extensive configurational sampling. In this work, we present and test various methodologies based on time-dependent density functional theory (TDDFT) and quantum mechanics/molecular mechanics (QM/MM) dynamics that take all of these requirements into account to provide a quantitative prediction of the effect of the environment on the fluorescence quantum yield of indole, a tryptophan fluorophore. This investigation paves the way for applications to the realistic spectroscopic characterization of the local protein environment of tryptophan from computer simulations.

Vibronic Dynamics of Photodissociating ICN from Simulations of Ultrafast X-Ray Absorption Spectroscopy

Ultrafast UV-pump/soft-X-ray-probe spectroscopy is a subject of great interest since it can provide detailed information about dynamical photochemical processes with ultrafast resolution and atomic specificity. Here, we focus on the photodissociation of ICN in the 1 Π1 excited state, with emphasis on the transient response in the soft-X-ray spectral region as described by the ab initio spectral lineshape averaged over the nuclear wavepacket probability density. We find that the carbon K-edge spectral region reveals a rich transient response that provides direct insights into the dynamics of frontier orbitals during the I-CN bond cleavage process. The simulated UV-pump/soft-X-ray-probe spectra exhibit detailed dynamical information, including a time-domain signature for coherent vibration associated with the photogenerated CN fragment.

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.

Reaction Path Analysis from Potential Energy Contributions Using Forces: An Accessible Estimator of Reaction Coordinate Adequacy

The calculation of potential energy and free-energy profiles along complex chemical reactions or rare event processes is of great interest because of their importance for many areas in chemistry, molecular biology, and material science. One typical way to generate these profiles is to add a bias potential to modify the energy surface, which can act on a selected degree of freedom in the system. However, in these cases, the quality of the result is strongly dependent on the selection of the degree of freedom over which this bias potential acts. The present work introduces a simple method for the analysis of the degree of freedom selected to describe a chemical process. The proposed methodology is based on the decomposition of contributions to the potential energy profiles by the integration of forces along a reaction path, which allows evaluating the different contributions to the energy change. This could be useful for discriminating the contributions to the energy arising from different regions of the system, which is particularly useful in systems with complex environments that must be represented using hybrid quantum mechanics/molecular mechanics schemes. Furthermore, this methodology allows in generating a quick and simple analysis of the degree of freedom which is used to describe the potential energy profile associated with the reactive process. This is computationally more accessible than the corresponding free-energy profile and can therefore be used as a simple estimator of reaction coordinate adequacy.

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.

Structure of the Human ACP-ISD11 Heterodimer

In recent years, the mammalian mitochondrial protein complex for iron-sulfur cluster assembly has been the focus of important studies. This is partly because of its high degree of relevance in cell metabolism and because mutations of the involved proteins are the cause of several human diseases. Cysteine desulfurase NFS1 is the key enzyme of the complex. At present, it is well-known that the active form of NFS1 is stabilized by the small protein ISD11. In this work, the structure of the human mitochondrial ACP-ISD11 heterodimer was determined at 2.0 Å resolution. ACP-ISD11 forms a cooperative unit stabilized by several ionic interactions, hydrogen bonds, and apolar interactions. The 4'-phosphopantetheine-acyl chain, which is covalently bound to ACP, interacts with several residues of ISD11, modulating together with ACP the foldability of ISD11. Recombinant human ACP-ISD11 was able to interact with the NFS1 desulfurase, thus yielding an active enzyme, and the NFS1/ACP-ISD11 core complex was activated by frataxin and ISCU proteins. Internal motions of ACP-ISD11 were studied by molecular dynamics simulations, showing the persistence of the interactions between both protein chains. The conformation of the dimer is similar to that found in the context of the (NFS1/ACP-ISD11)₂ supercomplex core, which contains the Escherichia coli ACP instead of the human variant. This fact suggests a sequential mechanism for supercomplex consolidation, in which the ACP-ISD11 complex may fold independently and, after that, the NFS1 dimer would be stabilized.