Heme Distortion Research Articles

Heme pocket modulates protein conformation and diguanylate cyclase activity of a tetrameric globin coupled sensor

Bacteria use the second messenger cyclic dimeric guanosine monophosphate (c-di-GMP) to control biofilm formation and other key phenotypes in response to environmental signals. Changes in oxygen levels can alter c-di-GMP signaling through a family of proteins termed globin coupled sensors (GCS) that contain diguanylate cyclase domains. Previous studies have found that GCS diguanylate cyclase activity is controlled by ligand binding to the heme within the globin domain, with oxygen binding resulting in the greatest increase in catalytic activity. Herein, we present evidence that heme-edge residues control O2-dependent signaling in PccGCS, a GCS protein from Pectobacterium carotovorum, by modulating heme distortion. Using enzyme kinetics, resonance Raman spectroscopy, small angle X-ray scattering, and multi-wavelength analytical ultracentrifugation, we have developed an integrated model of the full-length PccGCS tetramer and have identified conformational changes associated with ligand binding, heme conformation, and cyclase activity. Taken together, these studies provide new insights into the mechanism by which O2 binding modulates activity of diguanylate cyclase-containing GCS proteins.

Computational Exploration of Minimum Energy Reaction Pathway of N2O Formation from Intermediate I of P450nor Using an Active Center Model.

P450nor is a heme-containing enzyme that catalyzes the conversion of nitric oxide (NO) to nitrous oxide (N2O). Its catalytic mechanism has attracted attention in chemistry, biology, and environmental engineering. The catalytic cycle of P450nor is proposed to consist of three major steps. The reaction mechanism for the last step, N2O generation, remains unknown. In this study, the reaction pathway of the N2O generation from the intermediate I was explored with the B3LYP calculations using an active center model after the examination of the validity of the model. In the validation, we compared the heme distortions between P450nor and other oxidoreductases, suggesting a small effect of protein environment on the N2O generation reaction in P450nor. We then evaluated the electrostatic environment effect of P450nor on the hydride affinity to the active site with quantum mechanics/molecular mechanics (QM/MM) calculations, confirming that the affinity was unchanged with or without the protein environment. The active center model for P450nor showed that the N2O generation process in the enzymatic reaction undergoes a reasonable barrier height without protein environment. Consequently, our findings strongly suggest that the N2O generation reaction from the intermediate I depends sorely on the intrinsic reactivity of the heme cofactor bound on cysteine residue.

Random Forest Classifier of Heme Proteins Using Porphyrin Distortions and Axial Ligands of Heme

We achieved success in the construction of a random forest classifier that is capable of accurately classifying heme proteins for oxygen storage/transport function and redox function, with a high level of precision (>83%), by utilizing heme porphyrin distortions and axial ligands of heme as features. Our findings suggest that the heme distortions and the axial ligands are significantly correlated with oxidoreductase functions and oxygen storage/transport functions, respectively, all of which are commonly observed in heme proteins. We constructed a random forest classifier using heme porphyrin distortions and axial ligands as features to successfully classify and predict two distinct classes of heme protein functions with remarkable accuracy. Our findings revealed that heme porphyrin distortion played a significant role in the classification of oxidoreductases, while axial ligand was found to be crucial for the classification of oxygen-binding proteins.

Regulatory mechanism of formaldehyde release in heme degradation catalyzed by Staphylococcus aureus IsdG

IsdG-type enzymes catalyze the noncanonical degradation of heme to iron, staphylobilin (SB), and formaldehyde (HCHO), presumably by binding heme in an unusually distorted conformation. Their unique mechanism has been elucidated for MhuD from Mycobacterium tuberculosis, revealing an unusual ring opening of hydroxyheme by dioxygenation. A similar mechanism has been postulated for other IsdG enzymes; however, MhuD, which is special as an IsdG-type enzyme, retains a formyl group in the linearized tetrapyrrole. Recent reports on Staphylococcus aureus IsdG have suggested the formation of SB retaining a formyl group (formyl-SB), but its identification is preliminary. Furthermore, the reaction properties of formyl-SB and the mechanism of HCHO release remain unclear. In this study, the complex reaction of S.aureus IsdG was reexamined to elucidate its mechanism, including the identification of reaction products and their control mechanisms. Depending on the reaction conditions, IsdG produced both SB and formyl-SB as the main product, the latter of which was isolated and characterized by MS and NMR measurements. The formyl-SB product was generated upon the reaction between hydroxyheme-IsdG and O2 without reduction, indicating the dioxygenation mechanism as found for MhuD. Under reducing conditions, hydroxyheme-IsdG was converted also to SB and HCHO by activating another O2 molecule. These results provide the first overview of the complicated IsdG reaction. The heme distortion in the IsdG-type enzymes is shown to generally promote ring cleavage by dioxygenation. The presence or absence of HCHO release can be influenced by many factors, and the direct identification of S.aureus heme catabolites is of interest.

The H-NOX protein structure adapts to different mechanisms in sensors interacting with nitric oxide.

Some classes of bacteria within phyla possess protein sensors identified as homologous to the heme domain of soluble guanylate cyclase, the mammalian NO-receptor. Named H-NOX domain (Heme-Nitric Oxide or OXygen-binding), their heme binds nitric oxide (NO) and O2 for some of them. The signaling pathways where these proteins act as NO or O2 sensors appear various and are fully established for only some species. Here, we investigated the reactivity of H-NOX from bacterial species toward NO with a mechanistic point of view using time-resolved spectroscopy. The present data show that H-NOXs modulate the dynamics of NO as a function of temperature, but in different ranges, changing its affinity by changing the probability of NO rebinding after dissociation in the picosecond time scale. This fundamental mechanism provides a means to adapt the heme structural response to the environment. In one particular H-NOX sensor the heme distortion induced by NO binding is relaxed in an ultrafast manner (∼15 ps) after NO dissociation, contrarily to other H-NOX proteins, providing another sensing mechanism through the H-NOX domain. Overall, our study links molecular dynamics with functional mechanism and adaptation.

Elucidation of the Correlation between Heme Distortion and Tertiary Structure of the Heme-Binding Pocket Using a Convolutional Neural Network.

Heme proteins serve diverse and pivotal biological functions. Therefore, clarifying the mechanisms of these diverse functions of heme is a crucial scientific topic. Distortion of heme porphyrin is one of the key factors regulating the chemical properties of heme. Here, we constructed convolutional neural network models for predicting heme distortion from the tertiary structure of the heme-binding pocket to examine their correlation. For saddling, ruffling, doming, and waving distortions, the experimental structure and predicted values were closely correlated. Furthermore, we assessed the correlation between the cavity shape and molecular structure of heme and demonstrated that hemes in protein pockets with similar structures exhibit near-identical structures, indicating the regulation of heme distortion through the protein environment. These findings indicate that the tertiary structure of the heme-binding pocket is one of the factors regulating the distortion of heme porphyrin, thereby controlling the chemical properties of heme relevant to the protein function; this implies a structure–function correlation in heme proteins.

Global Analysis of Heme Proteins Elucidates the Correlation between Heme Distortion and the Heme-Binding Pocket.

Heme proteins play diverse and important biological roles, from electron transfer and chemical catalysis to oxygen transport and/or storage. Although the distortion of heme porphyrin correlates with the physical properties of heme, such as the redox potential and oxygen affinity, the relationship between heme distortion and the heme protein environment is unclear. Here, we tested the hypothesis that the protein environment of the heme-binding pocket determines heme distortion (conformation). We analyzed the correlations between the amino acid composition of the heme-binding pocket and the magnitude of heme distortion along 12 vibrational modes using machine learning. A correlation was detected in the three lowest vibrational modes. Analysis of heme distortions in nearly the same environments of the heme-binding pocket supported this notion. Our analyses indicate that the heme-binding pocket environment is a major factor impacting the distortion of heme porphyrin along the three lowest vibrational modes. In addition, statistical analysis of the distortion of heme porphyrin revealed that the peaks of distributions of the ruffling and breathing distortions are shifted from 0 (the equilibrium structure). Both the ruffling and breathing distortions are correlated with the redox potential of heme, so that heme molecules with these distortions have a lower redox potential than planar molecules. These findings explain the structure-function relationship of heme.

Analysis of Fluctuation in the Heme-Binding Pocket and Heme Distortion in Hemoglobin and Myoglobin.

Heme is located in the active site of proteins and has diverse and important biological functions, such as electron transfer and oxygen transport and/or storage. The distortion of heme porphyrin is considered an important factor for the diverse functions of heme because it correlates with the physical properties of heme, such as oxygen affinity and redox potential. Therefore, clarification of the relationship between heme distortion and the protein environment is crucial in protein science. Here, we analyzed the fluctuation in heme distortion in the protein environment for hemoglobin and myoglobin using molecular dynamics (MD) simulations and quantum mechanical (QM) calculations as well as statistical analysis of the protein structures of hemoglobin and myoglobin stored in Protein Data Bank. Our computation and statistical analysis showed that the protein environment for hemoglobin and myoglobin prominently affects the doming distortion of heme porphyrin, which correlates with its oxygen affinity, and that the magnitude of distortion is different between hemoglobin and myoglobin. These results suggest that heme distortion is affected by its protein environment and fluctuates around its fitted conformation, leading to physical properties that are appropriate for protein functions.

Interactive Regulation between Aliphatic Hydroxylation and Aromatic Hydroxylation of Thaxtomin D in TxtC: A Theoretical Investigation.

TxtC is an unusual bifunctional cytochrome P450 that is able to perform sequential aliphatic and aromatic hydroxylation of the diketopiperazine substrate thaxtomin D in two distinct sites to produce thaxtomin A. Though the X-ray structure of TxtC complexed with thaxtomin D revealed a binding mode for its aromatic hydroxylation, the preferential hydroxylation site is aliphatic C14. It is thus intriguing to unravel how TxtC accomplishes such two-step catalytic hydroxylation on distinct aliphatic and aromatic carbons and why the aliphatic site is preferred in the hydroxylation step. In this work, by employing molecular docking and molecular dynamics (MD) simulation, we revealed that thaxtomin D could adopt two different conformations in the TxtC active site, which were equal in energy with either the aromatic C20-H or aliphatic C14-H pointing toward the active Cpd I oxyferryl moiety. Further ONIOM calculations indicated that the energy barrier for the rate-limiting hydroxylation step on the aliphatic C14 site was 9.6 kcal/mol more favorable than that on the aromatic C20 site. The hydroxyl group on the monohydroxylated intermediate thaxtomin B C14 site formed hydrogen bonds with Ser280 and Thr385, which induced the l-Phe moiety to rotate around the Cβ-Cγ bond of the 4-nitrotryptophan moiety. Thus, it adopted an energetically favorable conformation with aromatic C20 adjacent to the oxyferryl moiety. In addition, the hydroxyl group induced solvent water molecules to enter the active site, which propelled thaxtomin B toward the heme plane and resulted in heme distortion. Based on this geometrical layout, the rate-limiting aromatic hydroxylation energy barrier decreased to 15.4 kcal/mol, which was comparable to that of the thaxtomin D aliphatic hydroxylation process. Our calculations indicated that heme distortion lowered the energy level of the lowest Cpd I α-vacant orbital, which promoted electron transfer in the rate-limiting thaxtomin B aromatic hydroxylation step in TxtC.

Ultrafast dynamics of heme distortion in the O2-sensor of a thermophilic anaerobe bacterium

Heme-Nitric oxide and Oxygen binding protein domains (H-NOX) are found in signaling pathways of both prokaryotes and eukaryotes and share sequence homology with soluble guanylate cyclase, the mammalian NO receptor. In bacteria, H-NOX is associated with kinase or methyl accepting chemotaxis domains. In the O2-sensor of the strict anaerobe Caldanaerobacter tengcongensis (Ct H-NOX) the heme appears highly distorted after O2 binding, but the role of heme distortion in allosteric transitions was not yet evidenced. Here, we measure the dynamics of the heme distortion triggered by the dissociation of diatomics from Ct H-NOX using transient electronic absorption spectroscopy in the picosecond to millisecond time range. We obtained a spectroscopic signature of the heme flattening upon O2 dissociation. The heme distortion is immediately (<1 ps) released after O2 dissociation to produce a relaxed state. This heme conformational change occurs with different proportions depending on diatomics as follows: CO < NO < O2. Our time-resolved data demonstrate that the primary structural event of allostery is the heme distortion in the Ct H-NOX sensor, contrastingly with hemoglobin and the human NO receptor, in which the primary structural events are respectively the motion of the proximal histidine and the rupture of the iron-histidine bond.

Horizontal and Vertical Push Effects in Saddled Zinc Porphyrin Complexes: Implications for Heme Distortion.

A heme oxygen binding behavior was described through a unique geometric and electronic comparison of zinc porphyrin complexes. In this work, a charge transfer model for saddled metalloporphyrin complexes outlined the push effects of the ring nonplanarity and axial imidazole, and the pull effect of the axial dioxygen. The origin and role of the horizontal (ring nonplanarity) push effect and its relationship to the vertical (axial ligand) push/pull effect and its contribution to dioxygen binding were considered from the perspectives of crystal structures, theoretical calculations, and bathochromic shifts. Single-point energy and molecular orbital calculations starting from crystal structures were used to obtain the electronic structures of zinc porphyrin complexes. This study not only revealed that the electronic behavior of metalloporphyrins is driven by ring nonplanarity and axial ligation but also afforded new insight into the oxygen carrier mechanism in heme.

Structural Insights into Oxygen-Dependent Signal Transduction within Globin Coupled Sensors.

In order to respond to external stimuli, bacteria have evolved sensor proteins linking external signals to intracellular outputs that can then regulate downstream pathways and phenotypes. Globin coupled sensor proteins (GCSs) serve to link environmental O2 levels to cellular processes by coupling a heme-containing sensor globin domain to a catalytic output domain. However, the mechanism by which O2 binding activates these proteins is currently unknown. To provide insights into the signaling mechanism, two distinct dimeric complexes of the isolated globin domain of the GCS from Bordetella pertussis ( BpeGlobin) were solved via X-ray crystallography in which differences in ligand-bound states were observed. Both monomers of one dimer contain Fe(II)-O2 states, while the other dimer consists of the Fe(III)-H2O and Fe(II)-O2 states. These data provide the first molecular insights into the heme pocket conformation of the active Fe(II)-O2 form of these enzymes. In addition, heme distortion modes and heme-protein interactions were found to correlate with the ligation state of the globin, suggesting that these conformational changes play a role in O2-dependent signaling. Fourier transform infrared spectroscopy (FTIR) of the full-length GCS from B. pertussis ( BpeGReg) and the closely related GCS from Pectobacterium carotovorum ssp. carotovorum ( PccGCS) confirmed the importance of an ordered water within the heme pocket and two distal residues (Tyr43 and Ser68) as hydrogen-bond donors. Taken together, this work provides mechanistic insights into BpeGReg O2 sensing and the signaling mechanisms of diguanylate cyclase-containing GCS proteins.

Statistical and quantum-chemical analysis of the effect of heme porphyrin distortion in heme proteins: Differences between oxidoreductases and oxygen carrier proteins

We studied the effects of heme distortions on the chemical properties in two protein classes, oxidoreductases and oxygen carrier proteins, in the Protein Data Bank by a combination of statistics and quantum chemistry. A statistical analysis of heme distortion provided the characteristic distortions. They were reasonably correlated with the redox potentials and oxygen affinities evaluated by quantum chemical calculations, implying that heme proteins utilize these distortions to enhance their functions. The combination of statistical analysis and quantum chemical calculations is useful for the elucidation of the structure-function relationships.

Geometric deconstruction of core and electron activation of a π-system in a series of deformed porphyrins: mimics of heme.

The predominant distortion of heme is responsible for its electronic activity, catalytic ability and spectral properties. In this work, altogether 12 new X-ray structures of saddled, waved and ruffled porphyrins are reported. Three types of deformed porphyrins as mimics of heme were evaluated and analyzed by geometric deconstruction, spectral comparison, and electrochemical tracking, which shows a unique relationship of deformation fashions and distortion degree to the geometry of the core and electron transfer ability of rings in these enzyme containing porphyrins. These mimics can adjust their core geometry for changing the structures of potential metals; while for rings themselves, they can also regulate the electron activity by switching the HOMO of the large π systems. These deformed porphyrins can be used as ideal mimics for heme. These findings help us to understand the principle and contribution of these deformations to electron transfer in catalytic oxidation and photoreactions. The nonplanar mimics have been synthesized through a modular synthetic approach under Adler-Longo or Lindsey condensation conditions.

Density functional study of porphyrin distortion effects on redox potential of heme.

Heme is involved in various biochemical roles in hemoproteins. In the present study, the effect of heme distortion on the redox potential was systematically investigated with density functional calculations. We focused on the ruffled and saddled distortions of heme, which correspond to the two lowest-frequency normal modes. Our computations demonstrated that the ruffled distortion tended to reduce the redox potential of heme and that the transition of the electronic configuration occurred from (dxz , dyz )3 (dxy )2 to (dxz , dyz )4 (dxy )1 . In contrast, the saddled distortion had a tendency toward an increase in the redox potential, and no transition of the electronic configuration occurred. In experiments, these tendencies were found in the relationship between with the heme distortions and the redox potentials in cytochrome c3 . © 2017 Wiley Periodicals, Inc.

Measurement of Heme Ruffling Changes in MhuD Using UV-vis Spectroscopy.

For decades it has been known that an out-of-plane ruffling distortion of heme perturbs its UV-vis absorption (Abs) spectrum, but whether increased ruffling induces a red or blue shift of the Soret band has remained a topic of debate. This debate has been resolved by the spectroscopic and computational characterization of Mycobacterium tuberculosis MhuD presented here, an enzyme that converts heme, oxygen, and reducing equivalents to nonheme iron and mycobilin. W66F and W66A MhuD have been characterized using (1)H nuclear magnetic resonance, Abs, and magnetic circular dichroism spectroscopies, and the data have been used to develop an experimentally validated theoretical model of ruffled, ferric heme. The PBE density functional theory (DFT) model that has been developed accurately reproduces the observed spectral changes from wild type enzyme, and the underlying quantum mechanical origins of these ruffling-induced changes were revealed by analyzing the PBE DFT description of the electronic structure. Small amounts of heme ruffling have no influence on the energy of the Q-band and blue-shift the Soret band due to symmetry-allowed mixing of the Fe 3dxy and porphyrin a2u orbitals. Larger amounts of ruffling red-shift both the Q and Soret bands due to disruption of π-bonding within the porphyrin ring.

Crystal Structure of CYP106A2 in Substrate-Free and Substrate-Bound Form.

CYP106A2 from Bacillus megaterium ATCC 13368 is known as a bacterial steroid hydroxylase that is also capable of hydroxylating a variety of terpenoids. To analyze the substrate specificity of this enzyme further, different resin acids of the abietane and pimarane types were tested with regard to binding and conversion. Product formation could be shown for all tested substrates. Spectroscopic studies revealed type I binding spectra for isopimaric acid, but dehydroabietic acid did not induce a high-spin shift of the enzyme. Interestingly, binding of abietic acid resulted in a type II difference spectrum typical for nitrogenous inhibitors. Co-crystallization of CYP106A2 with abietic acid and structure determination revealed bending of the heme cofactor when abietic acid was bound in the active site. Quantum chemical calculations strongly suggest that this heme distortion is the cause of the unusual spectroscopic characteristics.

Effects of protein structure on iron-polypeptide vibrational dynamic coupling in cytochrome c.

Cytochrome c (Cyt c) has a heme covalently bound to the polypeptide via a Cys-X-X-Cys-His (CXXCH) linker that is located in the interface region for protein–protein interactions. To determine whether the polypeptide matrix influences iron vibrational dynamics, nuclear resonance vibrational spectroscopy (NRVS) measurements were performed on 57Fe-labeled ferric Hydrogenobacter thermophilus cytochrome c-552, and variants M13V, M13V/K22M, and A7F, which have structural modifications that alter the composition or environment of the CXXCH pentapeptide loop. Simulations of the NRVS data indicate that the 150–325 cm–1 region is dominated by NHis–Fe–SMet axial ligand and polypeptide motions, while the 325–400 cm–1 region shows dominant contributions from ν(Fe–NPyr) (Pyr = pyrrole) and other heme-based modes. Diagnostic spectral signatures that directly relate to structural features of the heme active site are identified using a quantum chemistry-centered normal coordinate analysis (QCC-NCA). In particular, spectral features that directly correlate with CXXCH loop stiffness, the strength of the Fe–His interaction, and the degree of heme distortion are identified. Cumulative results from our investigation suggest that compared to the wild type (wt), variants M13V and M13V/K22M have a more rigid CXXCH pentapeptide segment, a stronger Fe–NHis interaction, and a more ruffled heme. Conversely, the A7F variant has a more planar heme and a weaker Fe–NHis bond. These results are correlated to the observed changes in reduction potential between wt protein and the variants studied here. Implications of these results for Cyt c biogenesis and electron transfer are also discussed.

Fractional transfer of a free unpaired electron to overcome energy barriers in the formation of Fe(4+) from Fe(3+) during the core contraction of macrocycles: implication for heme distortion.

The free unpaired electron in Fe(3+) ions cannot be directly removed, and needs a transfer pathway with at least four steps to overcome the high energy barriers to form Fe(4+) ions. Fine changes in the electronic structure of Fe(3+) ions on spin conversion were identified through a deeper analysis of the diffraction, spectral and electrochemical data for six non-planar iron porphyrins. Fe(3+) ions can form four d electron tautomers as the compression of the central ion is increased. This indicates that the Fe(3+) ion undergoes a multistep electron transfer where the total energy gap of electron transfer is split into several smaller gaps to form high-valent Fe(4+) ions. We find that the interchange of these four electron tautomers is clearly related to the core size of the macrocycle in the current series. The large energy barrier to produce iron(iv) complexes is overcome through a gradient effect of multiple energy levels. In addition, a possible porphyrin Fe(3+)˙ radical may be formed from its stable isoelectronic form, porphyrin Fe(3+), under strong core contraction. These results indicate the important role of heme distortion in its catalytic oxidation functions.

Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study.

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or β hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe–histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or β chains, prior to the quaternary R–T and T–R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe–His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem.2002, 98, 149). Experimentally, the νFeHis evolution is faster for β than for α chains, and pump–probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 μs for β but not for α hemes, an interval previously shown to be the first step in the R–T transition. In the α chain νFeHis dropped sharply at 20 μs, the final step in the R–T transition. The time courses are fully consistent with recent computational mapping of the R–T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A.2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.