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

Electronic Structures of Radical-Pair-Forming Cofactors in a Heliobacterial Reaction Center

Photosynthetic reaction centers (RCs) are membrane proteins converting photonic excitations into electric gradients. The heliobacterial RCs (HbRCs) are assumed to be the precursors of all known RCs, making them a compelling subject for investigating structural and functional relationships. A comprehensive picture of the electronic structure of the HbRCs is still missing. In this work, the combination of selective isotope labelling of 13C and 15N nuclei and the utilization of photo-CIDNP MAS NMR (photochemically induced dynamic nuclear polarization magic-angle spinning nuclear magnetic resonance) allows for highly enhanced signals from the radical-pair-forming cofactors. The remarkable magnetic-field dependence of the solid-state photo-CIDNP effect allows for observation of positive signals of the electron donor cofactor at 4.7 T, which is interpreted in terms of a dominant contribution of the differential relaxation (DR) mechanism. Conversely, at 9.4 T, the emissive signals mainly originate from the electron acceptor, due to the strong activation of the three-spin mixing (TSM) mechanism. Consequently, we have utilized two-dimensional homonuclear photo-CIDNP MAS NMR at both 4.7 T and 9.4 T. These findings from experimental investigations are corroborated by calculations based on density functional theory (DFT). This allows us to present a comprehensive investigation of the electronic structure of the cofactors involved in electron transfer (ET).

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.

Structural Determinants of Redox Conduction Favor Robustness over Tunability in Microbial Cytochrome Nanowires

Structural determinants of a 103-fold variation in electrical conductivity for helical homopolymers of tetra-, hexa-, and octa-heme cytochromes (named Omc- E, S, and Z, respectively) from Geobacter sulfurreducens are investigated with the Pathways model for electron tunneling, classical molecular dynamics, and hybrid quantum/classical molecular mechanics. Thermally averaged electronic couplings for through-space heme-to-heme electron transfer in the "nanowires" computed with density functional theory are ≤0.015 eV. Pathways analyses also indicate that couplings match within a factor of 5 for all "nanowires", but some alternative tunneling routes are found involving covalent protein backbone bonds (Omc- S and Z) or propionic acid-ligating His H-bonds on adjacent hemes (OmcZ). Reorganization energies computed from electrostatic vertical energy gaps or a version of the Marcus continuum expression parameterized on the total (donor + acceptor) solvent-accessible surface area typically agree within 20% and fall within the range 0.48-0.98 eV. Reaction free energies in all three "nanowires" are ≤|0.28| eV, even though Coulombic interactions primarily tune the site redox energies by 0.7-1.2 eV. Given the conserved energetic parameters, redox conductivity differs by < 103-fold among the cytochrome "nanowires". Redox currents do not exceed 3.0 × 10-3 pA at a physiologically relevant 0.1 V bias, with the slowest electron transfers being on a (μs) timescale much faster than typical (ms) enzymatic turnovers. Thus, the "nanowires" are proposed to be functionally robust to variations in structure that provide a habitat-customized protein interface. The 30 pA to 30 nA variation in conductivity previously reported from atomic force microscopy experiments is not intrinsic to the structures and/or does not result from the physiologically relevant redox conduction mechanism.

Structural Determinants of Redox Conduction Favor Robustness over Tunability in Microbial Cytochrome Nanowires

Helical homopolymers of multiheme cytochromes catalyze biogeochemically significant electron transfers with a reported 10 3 -fold variation in conductivity. Herein, classical molecular dynamics and hybrid quantum/classical molecular mechanics are used to elucidate the structural determinants of the redox potentials and conductivities of the tetra-, hexa-, and octaheme outer-membrane cytochromes E, S, and Z, respectively, from Geobacter sulfurreducens . Second-sphere electrostatic interactions acting on minimally polarized heme centers are found to regulate redox potentials over a computed 0.5-V range. However, the energetics of redox conduction are largely robust to the structural diversity: Single-step electronic couplings (⟨H mn ⟩), reaction free energies , and reorganization energies (λ mn ) are always respectively <|0.026|, <|0.26|, and between 0.5 - 1.0 eV. With these conserved parameter ranges, redox conductivity differed by less than a factor of 10 among the 'nanowires' and is sufficient to meet the demands of cellular respiration if 10 2 - 10 3 'nanowires' are expressed. The 'nanowires' are proposed to be differentiated by the protein packaging to interface with a great variety of environments, and not by conductivity, because the rate-limiting electron transfers are elsewhere in the respiratory process. Conducting-probe atomic force microscopy measurements that find conductivities 10 3 -10 6 -fold more than cellular demands are suggested to report on functionality that is either not used or not accessible under physiological conditions. The experimentally measured difference in conductivity between Omc- S and Z is suggested to not be an intrinsic feature of the CryoEM-resolved structures.

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.

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.

The redox potential of a heme cofactor in Nitrosomonas europaea cytochrome c peroxidase: a polarizable QM/MM study

Redox reactions are crucial to biological processes that protect organisms against oxidative stress. Metalloenzymes, such as peroxidases which reduce excess reactive oxygen species into water, play a key role in detoxification mechanisms. Here we present the results of a polarizable QM/MM study of the reduction potential of the electron transfer heme in the cytochrome c peroxidase of Nitrosomonas europaea. We have found that environment polarization does not substantially affect the computed value of the redox potential. Particular attention has been given to analyzing the role of electrostatic interactions within the protein environment and the solvent on tuning the redox potential of the heme co-factor. We have found that the electrostatic interactions predominantly explain the fluctuations of the vertical ionization/attachment energies of the heme for the sampled configurations, and that the long range electrostatic interactions (up to 40 Å) contribute substantially to the absolute values of the vertical energy gaps.

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/.

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.

Understanding the Reactivity of Mn-Oxo Porphyrins for Substrate Hydroxylation: Theoretical Predictions and Experimental Evidence Reconciled

The Mn-oxo porphyrin (MnOP) mechanism for substrate hydroxylation is computationally studied with the aim to better understand reactivity in these systems. Theoretical studies suggest Mn(V)OP species to be very reactive intermediates with thermally accessible reaction barriers represented by low-spin/high-spin-crossover occurring in the Mn(V)OP oxidant, and kinetics for selected Mn(V)OP species indeed find high reactivity. On the other hand, MnOP complexes lead to modest yields in hydroxylation reactions of several different substrates, implying low rate constants and high reaction barriers. The resolution of this inconsistency is very important to understand the reactivity of Mn-oxo porphyrins and to improve the catalytic conditions. In this work we use the toluene hydroxylation by the Mn(V)OP(H₂O)⁺ complex as a case study to gain deep insight into the reaction mechanism. Minimum energy crossing point (MECP) results on the H-abstraction process from toluene indicate a first crossover from a singlet to a triplet spin state of the Mn(V)OP(H₂O)⁺ species with a thermally accessible barrier, followed by a very facile H-abstraction by the triplet complex. Issues concerning (i) the validation of the level of the density functional theory employed (BP86) to describe the singlet-triplet energy gap in the Mn(V)OP(H₂O)⁺ system versus highly accurate DMRG-CASPT2/CC calculations, and (ii) the influence of the axial ligand (X = none, Cl^-, CH₃CN, OH^-, and O^2-) on MnOP reactivity, which models the different experimental conditions, are addressed. The ligand trans influence mainly controls the reactivity through the singlet-triplet energy gap modulation, with the porphyrin ruffling distortion also finely tuning it. Finally, a stepwise model for the H-abstraction process is proposed which allows a direct comparison between the calculated and experimentally measured Gibbs free activation energy barriers ( Zhang et al. J. Am. Chem. Soc. 2005 , 127 , 6573 - 6582 ). The low yields in catalysis are shown not to be due to low reactivity of Mn(V).

A Density Functional Theory-Based Scheme to Compute the Redox Potential of a Transition Metal Complex: Applications to Heme Compound

We estimated the redox potential of a model heme compound by using the combination of our density functionals with a computational scheme, which corrects the solvation energy to the normal solvent model. Among many density functionals, the LC-BOP12 functional gave the smallest mean absolute error of 0.16 V in the test molecular sets. The application of these methods revealed that the redox potential of a model heme can be controlled within 200 mV by changing the protonation state and even within 20 mV by the flipping of the ligand histidine. In addition, the redox potential depends on the inverse of the dielectric constant, which controls the surroundings. The computational results also imply that a system with a low dielectric constant avoids the charged molecule by controlling either the redox potential or the protonation system.