Luciana Capece

Evidence for a ferryl intermediate in a heme-based dioxygenase

In contrast to the wide spectrum of cytochrome P450 monooxygenases, there are only 2 heme-based dioxygenases in humans: tryptophan dioxygenase (hTDO) and indoleamine 2,3-dioxygenase (hIDO). hTDO and hIDO catalyze the same oxidative ring cleavage reaction of L-tryptophan to N-formyl kynurenine, the initial and rate-limiting step of the kynurenine pathway. Despite immense interest, the mechanism by which the 2 enzymes execute the dioxygenase reaction remains elusive. Here, we report experimental evidence for a key ferryl intermediate of hIDO that supports a mechanism in which the 2 atoms of dioxygen are inserted into the substrate via a consecutive 2-step reaction. This finding introduces a paradigm shift in our understanding of the heme-based dioxygenase chemistry, which was previously believed to proceed via simultaneous incorporation of both atoms of dioxygen into the substrate. The ferryl intermediate is not observable during the hTDO reaction, highlighting the structural differences between the 2 dioxygenases, as well as the importance of stereoelectronic factors in modulating the reactions.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

The first step of the dioxygenation reaction carried out by tryptophan dioxygenase and indoleamine 2,3-dioxygenase as revealed by quantum mechanical/molecular mechanical studies

Tryptophan dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO) are two heme-containing enzymes which catalyze the conversion of l-tryptophan to N-formylkynurenine (NFK). In mammals, TDO is mostly expressed in liver and is involved in controlling homeostatic serum tryptophan concentrations, whereas IDO is ubiquitous and is involved in modulating immune responses. Previous studies suggested that the first step of the dioxygenase reaction involves the deprotonation of the indoleamine group of the substrate by an evolutionarily conserved distal histidine residue in TDO and the heme-bound dioxygen in IDO. Here, we used classical molecular dynamics and hybrid quantum mechanical/molecular mechanical methods to evaluate the base-catalyzed mechanism. Our data suggest that the deprotonation of the indoleamine group of the substrate by either histidine in TDO or heme-bound dioxygen in IDO is not energetically favorable. Instead, the dioxygenase reaction can be initiated by a direct attack of heme-bound dioxygen on the C2=C3 bond of the indole ring, leading to a protein-stabilized 2,3-alkylperoxide transition state and a ferryl epoxide intermediate, which subsequently recombine to generate NFK. The novel sequential two-step oxygen addition mechanism is fully supported by our recent resonance Raman data that allowed identification of the ferryl intermediate (Lewis-Ballester et al. in Proc Natl Acad Sci USA 106:17371–17376, 2009). The results reveal the subtle differences between the TDO and IDO reactions and highlight the importance of protein matrix in modulating stereoelectronic factors for oxygen activation and the stabilization of both transition and intermediate states.

Complete Reaction Mechanism of Indoleamine 2,3-Dioxygenase as Revealed by QM/MM Simulations

Indoleamine 2,3-dioxygenase (IDO) and tryptophan dioxygenase (TDO) are two heme proteins that catalyze the oxidation reaction of tryptophan (Trp) to N-formylkynurenine (NFK). Human IDO (hIDO) has recently been recognized as a potent anticancer drug target, a fact that triggered intense research on the reaction and inhibition mechanisms of hIDO. Our recent studies revealed that the dioxygenase reaction catalyzed by hIDO and TDO is initiated by addition of the ferric iron-bound superoxide to the C(2)═C(3) bond of Trp to form a ferryl and Trp-epoxide intermediate, via a 2-indolenylperoxo radical transition state. The data demonstrate that the two atoms of dioxygen are inserted into the substrate in a stepwise fashion, challenging the paradigm of heme-based dioxygenase chemistry. In the current study, we used QM/MM methods to decipher the mechanism by which the second ferryl oxygen is inserted into the Trp-epoxide to form the NFK product in hIDO. Our results show that the most energetically favored pathway involves proton transfer from Trp-NH(3)(+) to the epoxide oxygen, triggering epoxide ring opening and a concerted nucleophilic attack of the ferryl oxygen to the C(2) of Trp that leads to a metastable reaction intermediate. This intermediate subsequently converts to NFK, following C(2)-C(3) bond cleavage and the associated back proton transfer from the oxygen to the amino group of Trp. A comparative study with Xantomonas campestris TDO (xcTDO) indicates that the reaction follows a similar pathway, although subtle differences distinguishing the two enzyme reactions are evident. The results underscore the importance of the NH(3)(+) group of Trp in the two-step ferryl-based mechanism of hIDO and xcTDO, by acting as an acid catalyst to facilitate the epoxide ring-opening reaction and ferryl oxygen addition to the indole ring.

High pressure reveals structural determinants for globin hexacoordination: Neuroglobin and myoglobin cases

The influence of pressure on the equilibrium between five‐(5c) and six‐coordination (6c) forms in neuroglobin (Ngb) and myoglobin (Mb) has been examined by means of molecular dynamics (MD) simulations at normal and high pressure. The results show that the main effect of high pressure is to reduce the protein mobility without altering the structure in a significant manner. Moreover, our data suggest that the equilibrium between 5c and 6c states in globins is largely controlled by the structure and dynamics of the C‐D region. Finally, in agreement with the available experimental data, the free energy profiles obtained from steered MD for both proteins indicate that high pressure enhances hexacoordination. In Ngb, the shift in equilibrium is mainly related to an increase in the 6c→5c transition barrier, whereas in Mb such a shift is primarily due to a destabilization of the 5c state. Proteins 2009.

Dynamical Characterization of the Heme NO Oxygen Binding (HNOX) Domain. Insight into Soluble Guanylate Cyclase Allosteric Transition

Since the discovery of soluble guanylate cyclase (sGC) as the mammalian receptor for nitric oxide (NO), numerous studies have been performed in order to understand how sGC transduces the NO signal. However, the structural basis of sGC activation is still not completely elucidated. Spectroscopic and kinetic studies showed that the key step in the activation mechanism was the NO-induced breaking of the iron proximal histidine bond in the so-called 6c-NO to 5c-NO transition. The main breakthrough in the understanding of sGC activation mechanism came, however, from the elucidation of crystal structures for two different prokaryotic heme NO oxygen (HNOX) domains, which are homologues to the sGC heme domain. In this work we present computer simulation results of Thermoanaerobacter tencogensis HNOX that complement these structural studies, yielding molecular explanations to several poorly understood properties of these proteins. Specifically, our results explain the differential ligand binding patterns of the HNOX domains according to the nature of proximal and distal residues. We also show that the natural dynamics of these proteins is intimately related with the proposed conformational dependent activation process, which involves mainly the alphaFbeta1 loop and the alphaA-alphaC distal subdomain. The results from the sGC models also support this view and suggest a key role for the alphaFbeta1 loop in the iron proximal histidine bond breaking process and, therefore, in the sGC activation mechanism.

Protein dynamics and ligand migration interplay as studied by computer simulation.

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Oxygen affinity controlled by dynamical distal conformations: The soybean leghemoglobin and the Paramecium caudatum hemoglobin cases

The binding of diatomic ligands, such as O2, NO, and CO, to heme proteins is a process intimately related with their function. In this work, we analyzed by means of a combination of classical Molecular Dynamics (MD) and Hybrid Quantum–Classical (QM/MM) techniques the existence of multiple conformations in the distal site of heme proteins and their influence on oxygen affinity regulation. We considered two representative examples: soybean leghemoglobin (Lba) and Paramecium caudatum truncated hemoglobin (PcHb). The results presented in this work provide a molecular interpretation for the kinetic, structural, and mutational data that cannot be obtained by assuming a single distal conformation. Proteins 2007.

Small ligand–globin interactions: Reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Substrate Stereo-specificity in Tryptophan dioxygenase and Indoleamine 2,3- dioxygenase

The first and rate‐limiting step of the kynurenine pathway, in which tryptophan (Trp) is converted to N‐formylkynurenine is catalyzed by two heme‐containing proteins, Indoleamine 2,3‐dioxygenase (IDO), and Tryptophan 2,3‐dioxygenase (TDO). In mammals, TDO is found exclusively in liver tissue, IDO is found ubiquitously in all tissues. IDO has become increasingly popular in pharmaceutical research as it was found to be involved in many physiological situations, including immune escape of cancer. More importantly, small‐molecule inhibitors of IDO are currently utilized in cancer therapy. One of the main concerns for the design of human IDO (hIDO) inhibitors is that they should be selective enough to avoid inhibition of TDO. In this work, we have used a combination of classical molecular dynamics (MD) and hybrid quantum‐classical (QM/MM) methodologies to establish the structural basis that determine the differences in (a) the interactions of TDO and IDO with small ligands (CO/O2) and (b) the substrate stereo‐specificity in hIDO and TDO. Our results indicate that the differences in small ligand bound structures of IDO and TDO arise from slight differences in the structure of the bound substrate complex. The results also show that substrate stereo‐specificity of TDO is achieved by the perfect fit of L‐Trp, but not D‐Trp, which exhibits weaker interactions with the protein matrix. For hIDO, the presence of multiple stable binding conformations for L/D‐Trp reveal the existence of a large and dynamic active site. Taken together, our data allow determination of key interactions useful for the future design of more potent hIDO‐selective inhibitors. Proteins 2010;