Binding and docking interactions of NO, CO and O₂in heme proteins as probed by density functional theory

Mycobacterial and Human Ferrous Nitrobindins: Spectroscopic and Reactivity Properties

Structural and functional properties of ferrous Mycobacterium tuberculosis (Mt-Nb) and human (Hs-Nb) nitrobindins (Nbs) were investigated. At pH 7.0 and 25.0 °C, the unliganded Fe(II) species is penta-coordinated and unlike most other hemoproteins no pH-dependence of its coordination was detected over the pH range between 2.2 and 7.0. Further, despite a very open distal side of the heme pocket (as also indicated by the vanishingly small geminate recombination of CO for both Nbs), which exposes the heme pocket to the bulk solvent, their reactivity toward ligands, such as CO and NO, is significantly slower than in most hemoproteins, envisaging either a proximal barrier for ligand binding and/or crowding of H₂O molecules in the distal side of the heme pocket which impairs ligand binding to the heme Fe-atom. On the other hand, liganded species display already at pH 7.0 and 25 °C a severe weakening (in the case of CO) and a cleavage (in the case of NO) of the proximal Fe-His bond, suggesting that the ligand-linked movement of the Fe(II) atom onto the heme plane brings about a marked lengthening of the proximal Fe-imidazole bond, eventually leading to its rupture. This structural evidence is accompanied by a marked enhancement of both ligands dissociation rate constants. As a whole, these data clearly indicate that structural–functional relationships in Nbs strongly differ from what observed in mammalian and truncated hemoproteins, suggesting that Nbs play a functional role clearly distinct from other eukaryotic and prokaryotic hemoproteins.

Spin Interconversion of Heme-Peroxo-Copper Complexes Facilitated by Intramolecular Hydrogen-Bonding Interactions

Synthetic peroxo-bridged high-spin (HS) heme-(μ-η2:η1-O22-)-Cu(L) complexes incorporating (as part of the copper ligand) intramolecular hydrogen-bond (H-bond) capabilities and/or steric effects are herein demonstrated to affect the complex's electronic and geometric structure, notably impacting the spin state. An H-bonding interaction with the peroxo core favors a low-spin (LS) heme-(μ-η1:η1-O22-)-Cu(L) structure, resulting in a reversible temperature-dependent interconversion of spin state (5 coordinate HS to 6 coordinate LS). The LS state dominates at low temperatures, even in the absence of a strong trans-axial heme ligand. Lewis base addition inhibits the H-bond facilitated spin interconversion by competition for the H-bond donor, illustrating the precise H-bonding interaction required to induce spin-crossover (SCO). Resonance Raman spectroscopy (rR) shows that the H-bonding pendant interacts with the bridging peroxide ligand to stabilize the LS but not the HS state. The H-bond (to the Cu-bound O atom) acts to weaken the O-O bond and strengthen the Fe-O bond, exhibiting ν(M-O) and ν(O-O) values comparable to analogous known LS complexes with a strong donating trans-axial ligand, 1,5-dicyclohexylimidazole, (DCHIm)heme-(μ-η1:η1-O22-)-Cu(L). Variable-temperature (-90 to -130 °C) UV-vis and 2H NMR spectroscopies confirm the SCO process and implicate the involvement of solvent binding. Examining a case of solvent binding without SCO, thermodynamic parameters were obtained from a van't Hoff analysis, accounting for its contribution in SCO. Taken together, these data provide evidence for the H-bond group facilitating a core geometry change and allowing solvent to bind, stabilizing a LS state. The rR data, complemented by DFT analysis, reveal a stronger H-bonding interaction with the peroxo core in the LS compared to the HS complexes, which enthalpically favors the LS state. These insights enhance our fundamental understanding of secondary coordination sphere influences in metalloenzymes.

Local Energy Decomposition of Open-Shell Molecular Systems in the Domain-Based Local Pair Natural Orbital Coupled Cluster Framework

Local energy decomposition (LED) analysis decomposes the interaction energy between two fragments calculated at the domain-based local pair natural orbital CCSD(T) (DLPNO-CCSD(T)) level of theory into a series of chemically meaningful contributions and has found widespread applications in the study of noncovalent interactions. Herein, an extension of this scheme that allows for the analysis of interaction energies of open-shell molecular systems calculated at the UHF-DLPNO-CCSD(T) level is presented. The new scheme is illustrated through applications to the CH₂···X (X = He, Ne, Ar, Kr, and water) and heme···CO interactions in the low-lying singlet and triplet spin states. The results are used to discuss the mechanism that governs the change in the singlet–triplet energy gap of methylene and heme upon adduct formation.

The reductive phase of Rhodobacter sphaeroides cytochrome c oxidase disentangled by CO ligation

Cytochrome c oxidase (CcO) is a membrane protein of the respiratory chain that catalytically reduces molecular oxygen (O2) to water while translocating protons across the membrane. The enzyme hosts two copper and two heme iron moieties (heme a/heme a3). The atomic details of the sequential steps that go along with this redox-driven proton translocation are a matter of debate. Particularly for the reductive phase of CcO that precedes oxygen binding experimental data are scarce. Here, we use CcO under anaerobic conditions where carbon monoxide (CO) is bound to heme a3 which in tandem with CuB forms the binuclear center (BNC). Fourier-transform infrared (FTIR) absorption spectroscopy is combined with electro-chemistry to probe different redox and protonation states populated by variation of the external electrostatic potential. With this approach, the redox behavior of heme a and the BNC could be separated and the corresponding redox potentials were determined. We also infer the protonation of one of the propionate side chains of heme a3 to correlate with the oxidation of heme a. Experimental changes in the local electric field surrounding CO bound to heme a3 are determined by their vibrational Stark effect and agree well with electrostatic computations. The comparison of experimental and computational results indicates that changes of the heme a3/CuB redox state are coupled to proton transfer towards heme a3. The latter supports the role of the heme a3 propionate D as proton loading site.

ns-μs Time-Resolved Step-Scan FTIR of ba₃ Oxidoreductase from Thermus thermophilus: Protonic Connectivity of w941-w946-w927

Time-resolved step-scan FTIR spectroscopy has been employed to probe the dynamics of the ba₃ oxidoreductase from Thermus thermophilus in the ns-μs time range and in the pH/pD 6–9 range. The data revealed a pH/pD sensitivity of the D372 residue and of the ring-A propionate of heme a₃. Based on the observed transient changes a model in which the protonic connectivity of w941-w946-927 to the D372 and the ring-A propionate of heme a₃ is described.

Nanosecond ligand migration and functional protein relaxation in ba3 oxidoreductase: Structures of the B0, B1 and B2 intermediate states

Nanosecond time-resolved step-scan FTIR spectroscopy (nTRS (2) -FTIR) has been applied to literally probe the active site of the carbon monoxide (CO)-bound thermophilic ba3 heme-copper oxidoreductase as it executes its function. The nTRS (2) - snapshots of the photolysed heme a3 Fe-CO/CuB species captured a "transition state" whose side chains prevent the photolysed CO to enter the docking cavity. There are three sets of ba3 photoproduct bands of docked CO with different orientation exhibiting different kinetics. The trajectories of the "docked" CO at 2122, 2129 and 2137cm(-1) is referred to in the literature as B2, B1 and B0 intermediate states, respectively. The present data provided direct evidence for the role of water in controlling ligand orientation in an intracavity protein environment.

Docking and migration of carbon monoxide in nitrogenase: the case for gated pockets from infrared spectroscopy and molecular dynamics

Evidence of a CO docking site near the FeMo cofactor in nitrogenase has been obtained by Fourier transform infrared spectroscopy-monitored low-temperature photolysis. We investigated the possible migration paths for CO from this docking site using molecular dynamics calculations. The simulations support the notion of a gas channel with multiple internal pockets from the active site to the protein exterior. Travel between pockets is gated by the motion of protein residues. Implications for the mechanism of nitrogenase reactions with CO and N2 are discussed.

Using spectroscopic tools to probe porphyrin deformation and porphyrin-protein interactions

The reactivity and functionality of heme proteins are to a significant extent determined by the conformation of their functional heme groups and by the interaction of axial ligands with their protein environment. This review focuses on experimental methods and theoretical concepts for elucidating symmetry lowering perturbations of the heme induced by the protein environment of the heme pocket. First, we discuss a variety of methods which can be used to probe the electric field at the heme, including spectral hole burning as well as low temperature absorption and room temperature circular dichroism spectroscopy. Second, we show how heme deformations can be described as superposition of deformations along normal coordinates, thereby using the irreducible representations of the D4h point group as a classification tool. Finally, resonance Raman spectroscopy is introduced as a tool to probe the deformations of metalloprophyrins in solution and in protein matrices by measuring and comparing intensities and depolarization properties rather than wavenumber positions.

Role of the distal hydrogen-bonding network in regulating oxygen affinity in the truncated hemoglobin III from Campylobacter jejuni

Oxygen affinity in heme-containing proteins is determined by a number of factors, such as the nature and conformation of the distal residues that stabilize the heme bound-oxygen via hydrogen-bonding interactions. The truncated hemoglobin III from Campylobacter jejuni (Ctb) contains three potential hydrogen-bond donors in the distal site: TyrB10, TrpG8, and HisE7. Previous studies suggested that Ctb exhibits an extremely slow oxygen dissociation rate due to an interlaced hydrogen-bonding network involving the three distal residues. Here we have studied the structural and kinetic properties of the G8(WF) mutant of Ctb and employed state-of-the-art computer simulation methods to investigate the properties of the O(2) adduct of the G8(WF) mutant, with respect to those of the wild-type protein and the previously studied E7(HL) and/or B10(YF) mutants. Our data indicate that the unique oxygen binding properties of Ctb are determined by the interplay of hydrogen-bonding interactions between the heme-bound ligand and the surrounding TyrB10, TrpG8, and HisE7 residues.

Variation and analysis of second-sphere interactions and axial histidinate character in c-type cytochromes

The electron-donating properties of the axial His ligand to heme iron in cytochromes c (cyts c) are found to be correlated with the midpoint reduction potential (E(m)) in variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht cyt c(552)) in which mutations have been made in and near the Cys-X-X-Cys-His (CXXCH) heme-binding motif. To probe the strength of the His-Fe(III) interaction, we have measured (13)C nuclear magnetic resonance (NMR) chemical shifts for (13)CN(-) bound to heme iron trans to the axial His in Ht Fe(III) cyt c(552) variants. We observe a linear relationship between these (13)C chemical shifts and E(m), indicating that the His-Fe(III) bond strength correlates with E(m). To probe a conserved hydrogen bonding interaction between the axial His Hdelta1 and the backbone carbonyl of a Pro residue, we measured the pK(a) of the axial His Hdelta1 proton (pK(a(2))), which we propose to relate to the His-Fe(III) interaction, reduction potential, and local electrostatic effects. The observed linear relationship between the axial His (13)Cbeta chemical shift and E(m) is proposed to reflect histidinate (anionic) character of the ligand. A linear relationship also is seen between the average heme methyl (1)H chemical shift and E(m) which may reflect variation in axial His electron-donating properties or in the ruffling distortion of the heme plane. In summary, chemical shifts of the axial His and exogenous CN(-) bound trans to His are shown to be sensitive probes of the His-Fe(III) interaction in variants of Ht cyt c(552) and display trends that correlate with E(m).