Role of Heme Distortion on Oxygen Affinity in Heme Proteins: The Protoglobin Case

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restraints for nonprotein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in chlorophyll proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C2 symmetry of photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and cryogenic electron microscopy PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to a 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restraints play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested from cyanobacteria to algae. These results highlight both the opportunities and limitations inherent in structure-based analyses of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

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.

MicroED Structure of a Protoglobin Reactive Carbene Intermediate

Microcrystal electron diffraction (MicroED) is an emerging technique that has shown great potential for describing new chemical and biological molecular structures. Several important structures of small molecules, natural products, and peptides have been determined using ab initio methods. However, only a couple of novel protein structures have thus far been derived by MicroED. Taking advantage of recent technological advances, including higher acceleration voltage and using a low-noise detector in counting mode, we have determined the first structure of an Aeropyrum pernix protoglobin (ApePgb) variant by MicroED using an AlphaFold2 model for phasing. The structure revealed that mutations introduced during directed evolution enhance carbene transfer activity by reorienting an α helix of ApePgb into a dynamic loop, making the catalytic active site more readily accessible. After exposing the tiny crystals to the substrate, we also trapped the reactive iron-carbenoid intermediate involved in this engineered ApePgb's new-to-nature activity, a challenging carbene transfer from a diazirine via a putative metallo-carbene. The bound structure discloses how an enlarged active site pocket stabilizes the carbene bound to the heme iron and, presumably, the transition state for the formation of this key intermediate. This work demonstrates that improved MicroED technology and the advancement in protein structure prediction now enable investigation of structures that was previously beyond reach.

Interaction of nitrite with ferric protoglobin from Methanosarcina acetivorans - an interesting model for spectroscopic studies of the haem-ligand interaction

Protoglobin from Methanosarcina acetivorans (MaPgb) is a dimeric globin belonging to the same lineage of the globin superfamily as globin-coupled sensors. A putative role in the scavenging of reactive nitrogen and oxygen species has been suggested as a possible adaptation mechanism of the host organism to different gaseous environments in the course of evolution. A combination of optical absorption, electronic circular dichroism (ECD), resonance Raman (rRaman), and electron paramagnetic resonance (EPR) reveal the unusual in vitro reaction of ferric MaPgb with nitrite. In contrast to other globins, a large excess of nitrite did not induce the formation of a nitriglobin form in MaPgb. Surprisingly, the addition of nitrite in mildly acidic pH led to the formation of a stable nitric-oxide ligated ferric form of the protein (MaPgb-NO). Furthermore, the 300-700 nm ECD spectrum of ferric MaPgb is for the first time reported and discussed, showing strong differences in the Soret and Q ellipticity compared to ferric myoglobin, in line with the unusually strongly ruffled haem group of MaPgb and the related quantum-mechanical admixture of the S = 5/2 and S = 3/2 state of its ferric form. The Soret and Q ellipticity change strongly upon formation of MaPgb-NO, revealing a significant effect of the nitric-oxide ligation on the haem group and pocket. The related changes in the asymmetric pyrrole half-ring stretching vibration modes observed in the rRaman spectra give experimental support to earlier theoretical models, in which an important role of the in-plane breathing modes of the haem was predicted for the stabilization of the binding of diatomic gases to MaPgb.

How Nanoparticles Modify Adsorbed Proteins: Impact of Silica Nanoparticles on the Hemoglobin Active Site

The adsorption of proteins on surfaces has been studied for a long time, but the relationship between the structural and functional properties of the adsorbed protein and the adsorption mechanism remains unclear. Using hemoglobin adsorbed on silica nanoparticles, we have previously shown that hemoglobin's affinity towards oxygen increases with adsorption. Nevertheless, it was also shown that there were no significant changes in the quaternary and secondary structures. In order to understand the change in activity, we decided in this work to focus on the active sites of hemoglobin, the heme and its iron. After measuring adsorption isotherms of porcine hemoglobin on Ludox silica nanoparticles, we analyzed the structural modifications of adsorbed hemoglobin by X-ray absorption spectroscopy and circular dichroism spectra in the Soret region. It was found that upon adsorption, there were modifications in the heme pocket environment due to changes in the angles of the heme vinyl functions. These alterations can explain the greater affinity observed.

Prediction of Protein Function from Tertiary Structure of the Active Site in Heme Proteins by Convolutional Neural Network

Structure-function relationships in proteins have been one of the crucial scientific topics in recent research. Heme proteins have diverse and pivotal biological functions. Therefore, clarifying their structure-function correlation is significant to understand their functional mechanism and is informative for various fields of science. In this study, we constructed convolutional neural network models for predicting protein functions from the tertiary structures of heme-binding sites (active sites) of heme proteins to examine the structure-function correlation. As a result, we succeeded in the classification of oxygen-binding protein (OB), oxidoreductase (OR), proteins with both functions (OB-OR), and electron transport protein (ET) with high accuracy. Although the misclassification rate for OR and ET was high, the rates between OB and ET and between OB and OR were almost zero, indicating that the prediction model works well between protein groups with quite different functions. However, predicting the function of proteins modified with amino acid mutation(s) remains a challenge. Our findings indicate a structure-function correlation in the active site of heme proteins. This study is expected to be applied to the prediction of more detailed protein functions such as catalytic reactions.

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.

Biocatalytic Carbene Transfer Using Diazirines

Biocatalytic carbene transfer from diazo compounds is a versatile strategy in asymmetric synthesis. However, the limited pool of stable diazo compounds constrains the variety of accessible products. To overcome this restriction, we have engineered variants of Aeropyrum pernix protoglobin (ApePgb) that use diazirines as carbene precursors. While the enhanced stability of diazirines relative to their diazo isomers enables access to a diverse array of carbenes, they have previously resisted catalytic activation. Our engineered ApePgb variants represent the first example of catalysts for selective carbene transfer from these species at room temperature. The structure of an ApePgb variant, determined by microcrystal electron diffraction (MicroED), reveals that evolution has enhanced access to the heme active site to facilitate this new-to-nature catalysis. Using readily prepared aryl diazirines as model substrates, we demonstrate the application of these highly stable carbene precursors in biocatalytic cyclopropanation, N-H insertion, and Si-H insertion reactions.

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.

PyDISH: database and analysis tools for heme porphyrin distortion in heme proteins

Heme participates in a wide range of biological functions such as oxygen transport, electron transport, oxygen reduction, transcriptional regulation and so on. While the mechanism of each function has been investigated for many heme proteins, the origin of the diversity of the heme functions is still unclear and a crucial scientific issue. We have constructed a database of heme proteins, named Python-based database and analyzer for DIStortion of Heme porphyrin (PyDISH), which also contains some analysis tools. The aim of PyDISH is to integrate the information on the structures of hemes and heme proteins and the functions of heme proteins. This database will provide the structure-function relationships focusing on heme porphyrin distortion and lead to the elucidation of the origin of the functional diversity of heme proteins. In addition, the insights obtained from the database can be used for the design of protein function. PyDISH contains the structural data of more than 13 000 hemes extracted from the Protein Data Bank, including heme porphyrin distortion, axial ligands coordinating to the heme and the orientation of the propionate sidechains of heme. PyDISH also has information about the protein domains, including Uniprot ID, protein fold by CATH ID, organism, coordination distance and so on. The analytical tools implemented in PyDISH allow users to not only browse and download the data but also analyze the structures of heme porphyrin by using the analytical tools implemented in PyDISH. PyDISH users will be able to utilize the obtained results for the design of protein function. Database URL: http://pydish.bio.info.hiroshima-cu.ac.jp/.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

Iron transitions during activation of allosteric heme proteins in cell signaling

Allosteric heme proteins can fulfill a very large number of different functions thanks to the remarkable chemical versatility of heme through the entire living kingdom. Their efficacy resides in the ability of heme to transmit both iron coordination changes and iron redox state changes to the protein structure. Besides the properties of iron, proteins may impose a particular heme geometry leading to distortion, which allows selection or modulation of the electronic properties of heme. This review focusses on the mechanisms of allosteric protein activation triggered by heme coordination changes following diatomic binding to proteins as diverse as the human NO-receptor, cytochromes, NO-transporters and sensors, and a heme-activated potassium channel. It describes at the molecular level the chemical capabilities of heme to achieve very different tasks and emphasizes how the properties of heme are determined by the protein structure. Particularly, this reviews aims at giving an overview of the exquisite adaptability of heme, from bacteria to mammals.

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.

Distinctive structural properties of THB11, a pentacoordinate Chlamydomonas reinhardtii truncated hemoglobin with N- and C-terminal extensions

Hemoglobins (Hbs) utilize heme b as a cofactor and are found in all kingdoms of life. The current knowledge reveals an enormous variability of Hb primary sequences, resulting in topological, biochemical and physiological individuality. As Hbs appear to modulate their reactivities through specific combinations of structural features, predicting the characteristics of a given Hb is still hardly possible. The unicellular green alga Chlamydomonas reinhardtii contains 12 genes encoding diverse Hbs of the truncated lineage, several of which possess extended N- or C-termini of unknown function. Studies on some of the Chlamydomonas Hbs revealed yet unpredictable structural and biochemical variations, which, along with a different expression of their genes, suggest diverse physiological roles. Chlamydomonas thus represents a promising system to analyze the diversification of Hb structure, biochemistry and physiology. Here, we report the crystal structure, resolved to 1.75 Å, of the heme-binding domain of cyanomet THB11 (Cre16.g662750), one of the pentacoordinate algal Hbs, which offer a free Fe-coordination site in the reduced state. The overall fold of THB11 is conserved, but individual features such as a kink in helix E, a tilted heme plane and a clustering of methionine residues at a putative tunnel exit appear to be unique. Both N- and C-termini promote the formation of oligomer mixtures, and the absence of the C terminus results in reduced nitrite reduction rates. This work widens the structural and biochemical knowledge on the 2/2Hb family and suggests that the N- and C-terminal extensions of the Chlamydomonas 2/2Hbs modulate their reactivity by intermolecular interactions.

Diboron Porphyrins: The Raman Signature of the In-Plane Tetragonal Elongation of the Macrocycle

We describe an unusual in-plane type of porphyrin core distortion, tetragonal elongation (TE), observed experimentally in diboron porphyrins. The vibrational spectra of several of these complexes exhibit shifts that we have assigned to this TE distortion by comparing experimental spectra with DFT computational findings. The influence of TE in porphyrin systems was isolated using DFT analysis of the well-known model compounds Ni(II)porphine and Zn(II)porphine, with the macrocycle ring constrained to eliminate the influence of out-of-plane (OOP) distortions. A significant down-shift in frequencies was observed for porphyrin normal vibrational modes, particularly the in-plane A_1g/B_1g modes that are dominated by contributions from stretching and bending of C_α-C_m coordinates. In contrast, TE had little effect on the v(Pyr_halfring) and δ(Pyr_def) modes, though the lowered symmetry of the system resulted in significant splitting of the B_2u and B_3u modes. The impact of the TE distortion upon the diboron porphyrin vibrational spectrum was probed experimentally using Raman spectroscopy of B₂O₂(BCl₃)₂(TTP), B₂OF₂(TTP), and B₂OPhOH₂(TTP) (TTP = 5,10,15,20-(tetra- p-tolyl)porphyrin). Comparing the experimentally obtained spectral signatures to the computational findings allowed us to assign the large shifts observed for the v₂ and v₃ modes to the TE distortion in diboron porphyrins.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

Conformational control of cofactors in nature - the influence of protein-induced macrocycle distortion on the biological function of tetrapyrroles

Tetrapyrrole-containing proteins are one of the most fundamental classes of enzymes in nature and it remains an open question to give a chemical rationale for the multitude of biological reactions that can be catalyzed by these pigment-protein complexes. There are many fundamental processes where the same (i.e., chemically identical) porphyrin cofactor is involved in chemically quite distinct reactions. For example, heme is the active cofactor for oxygen transport and storage (hemoglobin, myoglobin) and for the incorporation of molecular oxygen in organic substrates (cytochrome P450). It is involved in the terminal oxidation (cytochrome c oxidase) and the metabolism of H2O2 (catalases and peroxidases) and catalyzes various electron transfer reactions in cytochromes. Likewise, in photosynthesis the same chlorophyll cofactor may function as a reaction center pigment (charge separation) or as an accessory pigment (exciton transfer) in light harvesting complexes (e.g., chlorophyll a). Whilst differences in the apoprotein sequences alone cannot explain the often drastic differences in physicochemical properties encountered for the same cofactor in diverse protein complexes, a critical factor for all biological functions must be the close structural interplay between bound cofactors and the respective apoprotein in addition to factors such as hydrogen bonding or electronic effects. Here, we explore how nature can use the same chemical molecule as a cofactor for chemically distinct reactions using the concept of conformational flexibility of tetrapyrroles. The multifaceted roles of tetrapyrroles are discussed in the context of the current knowledge on distorted porphyrins. Contemporary analytical methods now allow a more quantitative look at cofactors in protein complexes and the development of the field is illustrated by case studies on hemeproteins and photosynthetic complexes. Specific tetrapyrrole conformations are now used to prepare bioengineered designer proteins with specific catalytic or photochemical properties.

Ultrafast Dynamics of Carboxy-Hemoglobin: Two-Dimensional Infrared Spectroscopy Experiments and Simulations

This Letter presents a comparison between experimental and simulated 2D mid-infrared spectra of carboxy-hemoglobin in the spectral region of the carbon monoxide stretching mode. The simulations rely on a fluctuating potential energy surface that includes both the effect of heme and the protein surroundings computed from molecular dynamics simulations. A very good agreement between theory and experiment is obtained with no adjustable parameters. The simulations show that the effect of the distal histidine through the hydrogen bond is strong and is directly responsible for the slow decay of the frequency-frequency correlation function on a 10 ps time scale. This study confirms that fluctuations in carboxy-hemoglobin are more inhomogeneous than those in the more frequently studied carboxy-myoglobin. The comparison between simulations and experiments brings valuable information on the complex relation between protein structure and spectral diffusion.

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.

Nitrite-reductase and peroxynitrite isomerization activities of Methanosarcina acetivorans protoglobin

Within the globin superfamily, protoglobins (Pgb) belong phylogenetically to the same cluster of two-domain globin-coupled sensors and single-domain sensor globins. Multiple functional roles have been postulated for Methanosarcina acetivorans Pgb (Ma-Pgb), since the detoxification of reactive nitrogen and oxygen species might co-exist with enzymatic activity(ies) to facilitate the conversion of CO to methane. Here, the nitrite-reductase and peroxynitrite isomerization activities of the CysE20Ser mutant of Ma-Pgb (Ma-Pgb*) are reported and analyzed in parallel with those of related heme-proteins. Kinetics of nitrite-reductase activity of ferrous Ma-Pgb* (Ma-Pgb*-Fe(II)) is biphasic and values of the second-order rate constant for the reduction of NO₂^– to NO and the concomitant formation of nitrosylated Ma-Pgb*-Fe(II) (Ma-Pgb*-Fe(II)-NO) are k_app1 = 9.6±0.2 M^–1 s^–1 and k_app2 = 1.2±0.1 M^–1 s^–1 (at pH 7.4 and 20°C). The k_app1 and k_app2 values increase by about one order of magnitude for each pH unit decrease, between pH 8.3 and 6.2, indicating that the reaction requires one proton. On the other hand, kinetics of peroxynitrite isomerization catalyzed by ferric Ma-Pgb* (Ma-Pgb*-Fe(III)) is monophasic and values of the second order rate constant for peroxynitrite isomerization by Ma-Pgb*-Fe(III) and of the first order rate constant for the spontaneous conversion of peroxynitrite to nitrate are h_app = 3.8×10⁴ M^–1 s^–1 and h₀ = 2.8×10^–1 s^–1 (at pH 7.4 and 20°C). The pH-dependence of h_on and h₀ values reflects the acid-base equilibrium of peroxynitrite (pK_a = 6.7 and 6.9, respectively; at 20°C), indicating that HOONO is the species that reacts preferentially with the heme-Fe(III) atom. These results highlight the potential role of Pgbs in the biosynthesis and scavenging of reactive nitrogen and oxygen species.

Protoglobin: structure and ligand-binding properties

Protoglobin is the first globin identified in Archaea; its biological role is still unknown, although it can bind O2, CO and NO reversibly in vitro. The X-ray structure of Methanosarcina acetivorans protoglobin revealed several peculiar structural features. Its tertiary structure can be considered as an expanded version of the canonical globin fold, characterised by the presence of a pre-A helix (named Z) and a 20-residue N-terminal extension. Other unusual trends are a large distortion of the haem moiety, and its complete burial in the protein matrix due to the extended CE and FG loops and the 20-residue N-terminal loop. Access of diatomic ligands to the haem has been proposed to be granted by two tunnels, which are mainly defined by helices B/G (tunnel 1) and B/E (tunnel 2), and whose spatial orientation and topology give rise to an almost orthogonal two-tunnel system unprecedented in other globins. At a quaternary level, protoglobin forms a tight dimer, mostly based on the inter-molecular four-helix bundle built by the G- and H-helices, similar to that found in globin-coupled sensor proteins, which share with protoglobin a common phylogenetic origin. Such unique structural properties, together with an unusually low O2 dissociation rate and a selectivity ratio for O2/CO binding that favours O2 ligation, make protoglobin a peculiar case for gaining insight into structure to function relationships within the globin superfamily. While recent structural and biochemical data have given answers to important questions, the functional issue is still unclear and it is expected to represent the major focus of future investigations.

Porphyrin π-stacking in a heme protein scaffold tunes gas ligand affinity

The role of π-stacking in controlling redox and ligand binding properties of porphyrins has been of interest for many years. The recent discovery of H-NOX domains has provided a model system to investigate the role of porphyrin π-stacking within a heme protein scaffold. Removal of a phenylalanine-porphyrin π-stack dramatically increased O2, NO, and CO affinities and caused changes in redox potential (~40mV) without any structural changes. These results suggest that small changes in redox potential affect ligand affinity and that π-stacking may provide a novel route to engineer heme protein properties for new functions.

Structure and haem-distal site plasticity in Methanosarcina acetivorans protoglobin

Protoglobin from Methanosarcina acetivorans C2A (MaPgb), a strictly anaerobic methanogenic Archaea, is a dimeric haem-protein whose biological role is still unknown. As other globins, protoglobin can bind O₂, CO and NO reversibly in vitro, but it displays specific functional and structural properties within members of the hemoglobin superfamily. CO binding to and dissociation from the haem occurs through biphasic kinetics, which arise from binding to (and dissociation from) two distinct tertiary states in a ligation-dependent equilibrium. From the structural viewpoint, protoglobin-specific loops and a N-terminal extension of 20 residues completely bury the haem within the protein matrix. Thus, access of small ligand molecules to the haem is granted by two apolar tunnels, not common to other globins, which reach the haem distal site from locations at the B/G and B/E helix interfaces. Here, the roles played by residues Trp(60)B9, Tyr(61)B10 and Phe(93)E11 in ligand recognition and stabilization are analyzed, through crystallographic investigations on the ferric protein and on selected mutants. Specifically, protein structures are reported for protoglobin complexes with cyanide, with azide (also in the presence of Xenon), and with more bulky ligands, such as imidazole and nicotinamide. Values of the rate constant for cyanide dissociation from ferric MaPgb-cyanide complexes have been correlated to hydrogen bonds provided by Trp(60)B9 and Tyr(61)B10 that stabilize the haem-Fe(III)-bound cyanide. We show that protoglobin can strikingly reshape, in a ligand-dependent way, the haem distal site, where Phe(93)E11 acts as ligand sensor and controls accessibility to the haem through the tunnel system by modifying the conformation of Trp(60)B9.

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.

Marked difference in the electronic structure of cyanide-ligated ferric protoglobins and myoglobin due to heme ruffling

Electron paramagnetic resonance experiments reveal a significant difference between the principal g values (and hence ligand-field parameters) of the ferric cyanide-ligated form of different variants of the protoglobin of Methanosarcina acetivorans (MaPgb) and of horse heart myoglobin (hhMb). The largest principal g value of the ferric cyanide-ligated MaPgb variants is found to be significantly lower than for any of the other globins reported so far. This is at least partially caused by the strong heme distortions as proven by the determination of the hyperfine interaction of the heme nitrogens and mesoprotons. Furthermore, the experiments confirm recent theoretical predictions [Forti, F.; Boechi, L., Bikiel, D., Martí, M.A.; Nardini, M.; Bolognesi, M.; Viappiani, C.; Estrin, D.; Luque, F. J. J. Phys. Chem. B 2011, 115, 13771-13780] that Phe(G8)145 plays a crucial role in the ligand modulation in MaPgb. Finally, the influence of the N-terminal 20 amino-acid chain on the heme pocket in these protoglobins is also proven.

Steric bulkiness of pyrrole substituents and the out-of-plane deformations of porphyrins: nickel(II) octaisopropylporphyrin and its meso-nitro derivative

Sterically bulky substituents at the β-carbons of the pyrrole rings of porphyrins are sufficient to cause large out-of-plane porphyrin distortions even in the absence of substituent groups at the meso carbons. It is well established that substituents at the meso-positions only or at both the β-pyrrole and the meso-positions are sufficiently bulky to result in large non-planar distortions of the macrocycle. However, no systematic studies of the effects of bulky β-pyrrole substituents alone have been reported. Herein, molecular simulations and X-ray crystallography of nickel(II) 2,3,7,8,12,13,17,18-octa(isopropyl)porphyrin reveal that large out-of-plane distortions (>1.5 Å) are induced by the steric repulsion of the β-isopropyl groups but fail to lead to a single strongly energetically favored conformer. The molecular simulations indicate that multiple conformers differing in the orientation of the isopropyl groups and the macrocycle conformation coexist in solution and this is confirmed by resonance Raman spectroscopy. Large downshifts in the structure-sensitive lines result from the non-planar distortion, and line broadenings result from structural heterogeneity. The heterogeneity originates from tradeoffs between energy contributions from steric repulsion and macrocycle distortion. Simulations for 5-nitro-2,3,7,8,12,13,17,18-octa(isopropyl)porphyrin suggest two possible orientations of the nitro group with respect to the macrocycle mean plane — one nearly vertical (as in the crystal structure) and another that is nearly parallel. INDO/S semiempirical calculations indicate an orbital of the NO2 group resides between the porphyrin frontier orbitals with significant mixing of the nitro and porphyrin orbitals.KEYWORDS: porphyrin, non-planar, resonance Raman, X-ray crystallography, crystal structure, isopropyl, nitro, conformer, molecular mechanics, molecular simulations, density functional theory, steric crowding, conformational heterogeneity.

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.

Ligation tunes protein reactivity in an ancient haemoglobin: kinetic evidence for an allosteric mechanism in Methanosarcina acetivorans protoglobin

Protoglobin from Methanosarcina acetivorans (MaPgb) is a dimeric globin with peculiar structural properties such as a completely buried haem and two orthogonal tunnels connecting the distal cavity to the solvent. CO binding to and dissociation from MaPgb occur through a biphasic kinetics. We show that the heterogenous kinetics arises from binding to (and dissociation from) two tertiary conformations in ligation-dependent equilibrium. Ligation favours the species with high binding rate (and low dissociation rate). The equilibrium is shifted towards the species with low binding (and high dissociation) rates for the unliganded molecules. A quantitative model is proposed to describe the observed carbonylation kinetics.

Revealing substituent effects on the electronic structure and planarity of Ni-porphyrins

Using density functional theory, we have studied the effects on structural and electronic consequences (including HOMO-LUMO energy gaps, vertical ionization potentials (IP_v), and vertical electron affinities (EA_v)) of the following two factors: (a) meso- and β-substituents acting as inductive donors (CH₃), inductive acceptors that are electron-donating through resonance (Br), inductive electron acceptors (CF₃), and resonance enabled acceptors (NO₂); and (b) complete replacement of pyrrole nitrogens with P-atoms. The principal results of the study are: (1) For the bare Ni-porphyrin, the solvents were found not to affect the HOMO-LUMO gaps but to change the IP_v and EA_v noticeably. (2) In the series CH₃ → Br → CF₃ → NO₂ the HOMO-LUMO energy gaps, IP_v, and EA_v increase for both meso- and β-substituents. The ruffling distortion of the porphyrin core is retained, and becomes stronger for the two acceptor groups. In general, effects of meso-substituents on the ruffling distortion of the porphyrin core is more pronounced. (3) Most significantly, complete replacement of pyrrole nitrogens in the NiP with phosphorus atoms produces the species, NiP(P)₄, with the structural and electronic features drastically different from the original NiP. This implies that NiP(P)₄ can possess interesting and unusual novel properties, including aromaticity and reactivity, leading to its various beneficial potential applications. Furthermore, NiP(P)₄ high stability both in the gas phase and different solvents was shown, implying the feasibility of its synthesis.

Modulation of ligand-field parameters by heme ruffling in cytochromes c revealed by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectra of variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht c-552) and Pseudomonas aeruginosa cytochrome c(551) (Pa c-551) are analyzed to determine the effect of heme ruffling on ligand-field parameters. Mutations introduced at positions 13 and 22 in Ht c-552 were previously demonstrated to influence hydrogen bonding in the proximal heme pocket and to tune reduction potential (E(m)) over a range of 80 mV [Michel, L. V.; Ye, T.; Bowman, S. E. J.; Levin, B. D.; Hahn, M. A.; Russell, B. S.; Elliott, S. J.; Bren, K. L. Biochemistry 2007, 46, 11753-11760]. These mutations are shown here to also increase heme ruffling as E(m) decreases. The primary effect on electronic structure of increasing heme ruffling is found to be a decrease in the axial ligand-field term Δ/λ, which is proposed to arise from an increase in the energy of the d(xy) orbital. Mutations at position 7, previously demonstrated to influence heme ruffling in Pa c-551 and Ht c-552, are utilized to test this correlation between molecular and electronic structure. In conclusion, the structure of the proximal heme pocket of cytochromes c is shown to play a role in determining heme conformation and electronic structure.