Ligand preference and orientation in b- and c-type heme-binding proteins

Structural Basis for Monomer-Dimer Transition of Dri1 Upon Heme Binding

Domain related to iron (DRI) contains approximately 90 residues and is involved in iron and heme metabolism. Recent discoveries have annotated Dri1, a DRI-only protein from the cyanobacterium Synechocystis, as a regulator of succinate dehydrogenase in a b-type heme-dependent manner or as a c-type heme oxygenase. Here, we report high-resolution structures of Dri1 in complex with b-type and c-type hemes, respectively. Bis-His-ligated heme is located in the middle of the dimeric Dri1 complex with heme b, as well as in the complex of monomeric Dri1 with c-type heme, but distinct heme binding modes are revealed. Structural analyses suggest that Dri1 may participate in the succinate dehydrogenase activity and/or the metabolism of cytochromes.

Diversity and evolution of nitric oxide reduction in bacteria and archaea

Nitrous oxide is a potent greenhouse gas whose production is catalyzed by nitric oxide reductase (NOR) members of the heme-copper oxidoreductase (HCO) enzyme superfamily. We identified several previously uncharacterized HCO families, four of which (eNOR, sNOR, gNOR, and nNOR) appear to perform NO reduction. These families have novel active-site structures and several have conserved proton channels, suggesting that they might be able to couple NO reduction to energy conservation. We isolated and biochemically characterized a member of the eNOR family from the bacterium Rhodothermus marinus and found that it performs NO reduction. These recently identified NORs exhibited broad phylogenetic and environmental distributions, greatly expanding the diversity of microbes in nature capable of NO reduction. Phylogenetic analyses further demonstrated that NORs evolved multiple times independently from oxygen reductases, supporting the view that complete denitrification evolved after aerobic respiration.

A Chemical Proteomic Map of Heme-Protein Interactions

Heme is an essential cofactor for many human proteins as well as the primary transporter of oxygen in blood. Recent studies have also established heme as a signaling molecule, imparting its effects through binding with protein partners rather than through reactivity of its metal center. However, the comprehensive annotation of such heme-binding proteins in the human proteome remains incomplete. Here, we describe a strategy which utilizes a heme-based photoaffinity probe integrated with quantitative proteomics to map heme-protein interactions across the proteome. In these studies, we identified 350+ unique heme-protein interactions, the vast majority of which were heretofore unknown and consist of targets from diverse functional classes, including transporters, receptors, enzymes, transcription factors, and chaperones. Among these proteins is the immune-related interleukin receptor-associated kinase 1 (IRAK1), where we provide preliminary evidence that heme agonizes its catalytic activity. Our findings should improve the current understanding of heme's regulation as well as its signaling functions and facilitate new insights of its roles in human disease.

A noncanonical cytochrome c stimulates calcium binding by PilY1 for type IVa pili formation

Type IVa pili (T4aP) are versatile bacterial cell surface structures that undergo extension/adhesion/retraction cycles powered by the cell envelope-spanning T4aP machine. In this machine, a complex composed of four minor pilins and PilY1 primes T4aP extension and is also present at the pilus tip mediating adhesion. Similar to many several other bacteria, Myxococcus xanthus contains multiple minor pilins/PilY1 sets that are incompletely understood. Here, we report that minor pilins and PilY1 (PilY1.1) of cluster_1 form priming and tip complexes contingent on calcium and a noncanonical cytochrome c (TfcP) with an unusual His/Cys heme ligation. We provide evidence that TfcP is unlikely to participate in electron transport and instead stimulates calcium binding by PilY1.1 at low-calcium concentrations, thereby stabilizing PilY1.1 and enabling T4aP function in a broader range of calcium concentrations. These results not only identify a previously undescribed function of cytochromes c but also illustrate how incorporation of an accessory factor expands the environmental range under which the T4aP system 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.

Protein-Embedded Metalloporphyrin Arrays Templated by Circularly Permuted Tobacco Mosaic Virus Coat Proteins

Bioenergetic processes in nature have relied on networks of cofactors for harvesting, storing, and transforming the energy from sunlight into chemical bonds. Models mimicking the structural arrangement and functional crosstalk of the cofactor arrays are important tools to understand the basic science of natural systems and to provide guidance for non-natural functional biomaterials. Here, we report an artificial multiheme system based on a circular permutant of the tobacco mosaic virus coat protein (cpTMV). The double disk assembly of cpTMV presents a gap region sandwiched by the two C₂-symmetrically related disks. Non-native bis-his coordination sites formed by the mutation of the residues in this gap region were computationally screened and experimentally tested. A cpTMV mutant Q101H was identified to create a circular assembly of 17 protein-embedded hemes. Biophysical characterization using X-ray crystallography, cyclic voltammetry, and electron paramagnetic resonance (EPR) suggested both structural and functional similarity to natural multiheme cytochrome c proteins. This protein framework offers many further engineering opportunities for tuning the redox properties of the cofactors and incorporating non-native components bearing varied porphyrin structures and metal centers. Emulating the electron transfer pathways in nature using a tunable artificial system can contribute to the development of photocatalytic materials and bioelectronics.

Identifying metal binding amino acids based on backbone geometries as a tool for metalloprotein engineering

Metal cofactors within proteins perform a versatile set of essential cellular functions. In order to take advantage of the diverse functionality of metalloproteins, researchers have been working to design or modify metal binding sites in proteins to rationally tune the function or activity of the metal cofactor. This study has performed an analysis on the backbone atom geometries of metal-binding amino acids among 10 different metal binding sites within the entire protein data bank. A set of 13 geometric parameters (features) was identified that is capable of predicting the presence of a metal cofactor in the protein structure with overall accuracies of up to 97% given only the relative positions of their backbone atoms. The decision tree machine-learning algorithm used can quickly analyze an entire protein structure for the presence of sets of primary metal coordination spheres upon mutagenesis, independent of their original amino acid identities. The methodology was designed for application in the field of metalloprotein engineering. A cluster analysis using the data set was also performed and demonstrated that the features chosen are useful for identifying clusters of structurally similar metal-binding sites.

A novel catalytic heme cofactor in SfmD with a single thioether bond and a bis-His ligand set revealed by a de novo crystal structural and spectroscopic study

SfmD is a heme-dependent enzyme in the biosynthetic pathway of saframycin A. Here, we present a 1.78 Å resolution de novo crystal structure of SfmD, which unveils a novel heme cofactor attached to the protein with an unusual Hx n HxxxC motif (n ∼ 38). This heme cofactor is unique in two respects. It contains a single thioether bond in a cysteine-vinyl link with Cys317, and the ferric heme has two axial protein ligands, i.e., His274 and His313. We demonstrated that SfmD heme is catalytically active and can utilize dioxygen and ascorbate for a single-oxygen insertion into 3-methyl-l-tyrosine. Catalytic assays using ascorbate derivatives revealed the functional groups of ascorbate essential to its function as a cosubstrate. Abolishing the thioether linkage through mutation of Cys317 resulted in catalytically inactive SfmD variants. EPR and optical data revealed that the heme center undergoes a substantial conformational change with one axial histidine ligand dissociating from the iron ion in response to substrate 3-methyl-l-tyrosine binding or chemical reduction by a reducing agent, such as the cosubstrate ascorbate. The labile axial ligand was identified as His274 through redox-linked structural determinations. Together, identifying an unusual heme cofactor with a previously unknown heme-binding motif for a monooxygenase activity and the structural similarity of SfmD to the members of the heme-based tryptophan dioxygenase superfamily will broaden understanding of heme chemistry.

Superior Proton-Transfer Catalytic Promiscuity of Cytochrome c in Self-Organized Media

Evolutionarily elderly proteins commonly feature greater catalytic promiscuity. Cytochrome c is among the first set of proteins in evolution to have known prospects in electron transport and peroxidative properties. Here, we report that cyt c is also a proficient proton-transfer catalyst and enhances the Kemp elimination (KE; model reaction to show proton transfer catalytic property) by ∼750-fold on self-organized systems like micelles and vesicles. The self-organized systems mimic the mitochondrial environment in vitro for cyt c. Using an array of biophysical and biochemical mutational assays, both acid-base and redox mechanistic pathways have been explored. The histidine moiety close to hemin group (His18) is mainly responsible for proton abstraction to promote the concerted E2 pathway for KE catalysis when cyt c is in its oxidized form; this has also been confirmed by a H18A mutant of cyt c. However, the redox pathway is predominant under reducing conditions in the presence of dithiothreitol over the pH range 6-7.4. Interestingly, we found almost 750-fold enhanced KE catalysis by cyt c compared to aqueous buffer. Overall, in addition to providing mechanistic insights, the data reveal an unprecedented catalytic property of cyt c that could be of high importance in an evolutionary perspective considering its role in delineating the phylogenic tree and also towards generating programmable designer biocatalysts.

The pigment binding behaviour of water-soluble chlorophyll protein (WSCP)

Water-soluble chlorophyll proteins (WSCPs) are homotetrameric proteins that bind four chlorophyll (Chl) molecules in identical binding sites, which makes WSCPs a good model to study protein-pigment interactions. In a previous study, we described preferential binding of Chl a or Chl b in various WSCP versions. Chl b binding is preferred when a hydrogen bond can be formed between the C7 formyl of the chlorin macrocycle and the protein, whereas Chl a is preferred when Chl b binding is sterically unfavorable. Here, we determined the binding affinities and kinetics of various WSCP versions not only for Chl a/b, but also for chlorophyllide (Chlide) a/b and pheophytin (Pheo) a/b. Altered KD values are responsible for the Chl a/b selectivity in WSCP whereas differences in the reaction kinetics are neglectable in explaining different Chl a/b preferences. WSCP binds both Chlide and Pheo with a lower affinity than Chl, which indicates the importance of the phytol chain and the central Mg2+ ion as interaction sites between WSCP and pigment. Pheophorbide (Pheoide), lacking both the phytol chain and the central Mg2+ ion, can only be bound as Pheoide b to a WSCP that has a higher affinity for Chl b than Chl a, which underlines the impact of the C7 formyl-protein interaction. Moreover, WSCP was able to bind protochlorophyllide and Mg-protoporphyrin IX, which suggests that neither the size of the π electron system of the macrocycle nor the presence of a fifth ring at the macrocycle notably affect the binding to WSCP. WSCP also binds heme to form a tetrameric complex, suggesting that heme is bound in the Chl-binding site.

High-affinity binding and catalytic activity of His/Tyr-based sequences: Extending heme-regulatory motifs beyond CP

BACKGROUND & MOTIVATION

Peptides and proteins can interact with heme through His, Tyr, or Cys in heme-regulatory motifs (HRMs). The Cys-Pro dipeptide is a well investigated HRM, but for His and Tyr such a distinct motif is currently unknown. In addition, many heme-peptide complexes, such as heme-amyloid β, can display a peroxidase-like activity, albeit there is little understanding of how the local primary and secondary coordination environment influences catalytic activity. We thus systematically evaluated a series of His- and Tyr-based peptides to identify sequence features for high-affinity heme binding and their impact on the catalytic activity of heme.

METHODS

We employed solid-phase peptide synthesis to produce 58 nonapeptides, which were investigated by UV/vis, resonance Raman, and 2D NMR spectroscopy. A chromogenic assay was used to determine the catalytic activity of the heme-peptide complexes.

RESULTS

Heme-binding affinity and binding mode were found to be dependent on the coordinating amino acid and spacer length between multiple potential coordination sites in a motif. In particular, HXH and HXXXH motifs showed strong heme binding. Analysis of the peroxidase-like activity revealed that some of these peptides and also HXXXY motifs enhance the catalytic activity of heme significantly.

CONCLUSIONS

We identify HXH, HXXXH, and HXXXY as potential new HRMs with functional properties. Several peptides displayed a strikingly high peroxidase-like activity.

GENERAL SIGNIFICANCE

The identification of HRMs allows to discover yet unknown heme-regulated proteins, and consequently, enhances our current understanding of pathologies involving labile heme.

HeMoQuest: a webserver for qualitative prediction of transient heme binding to protein motifs

Background

The notion of heme as a regulator of many physiological processes via transient binding to proteins is one that is recently being acknowledged. The broad spectrum of the effects of heme makes it important to identify further heme-regulated proteins to understand physiological and pathological processes. Moreover, several proteins were shown to be functionally regulated by interaction with heme, yet, for some of them the heme-binding site(s) remain unknown. The presented application HeMoQuest enables identification and qualitative evaluation of such heme-binding motifs from protein sequences.

Results

We present HeMoQuest, an online interface (http://bit.ly/hemoquest) to algorithms that provide the user with two distinct qualitative benefits. First, our implementation rapidly detects transient heme binding to nonapeptide motifs from protein sequences provided as input. Additionally, the potential of each predicted motif to bind heme is qualitatively gauged by assigning binding affinities predicted by an ensemble learning implementation, trained on experimentally determined binding affinity data. Extensive testing of our implementation on both existing and new manually curated datasets reveal that our method produces an unprecedented level of accuracy (92%) in identifying those residues assigned “heme binding” in all of the datasets used. Next, the machine learning implementation for the prediction and qualitative assignment of binding affinities to the predicted motifs achieved 71% accuracy on our data.

Conclusions

Heme plays a crucial role as a regulatory molecule exerting functional consequences via transient binding to surfaces of target proteins. HeMoQuest is designed to address this imperative need for a computational approach that enables rapid detection of heme-binding motifs from protein datasets. While most existing implementations attempt to predict sites of permanent heme binding, this application is to the best of our knowledge, the first of its kind to address the significance of predicting transient heme binding to proteins.

Functionalization of Bacterial Microcompartment Shell Proteins With Covalently Attached Heme

Heme is a versatile redox cofactor that has considerable potential for synthetic biology and bioelectronic applications. The capacity to functionalize non-heme-binding proteins with covalently bound heme moieties in vivo could expand the variety of bioelectronic materials, particularly if hemes could be attached at defined locations so as to facilitate position-sensitive processes like electron transfer. In this study, we utilized the cytochrome maturation system I to develop a simple approach that enables incorporation of hemes into the backbone of target proteins in vivo. We tested our methodology by targeting the self-assembling bacterial microcompartment shell proteins, and inserting functional hemes at multiple locations in the protein backbone. We found substitution of three amino acids on the target proteins promoted heme attachment with high occupancy. Spectroscopic measurements suggested these modified proteins covalently bind low-spin hemes, with relative low redox midpoint potentials (about -210 mV vs. SHE). Heme-modified shell proteins partially retained their self-assembly properties, including the capacity to hexamerize, and form inter-hexamer attachments. Heme-bound shell proteins demonstrated the capacity to integrate into higher-order shell assemblies, however, the structural features of these macromolecular complexes was sometimes altered. Altogether, we report a versatile strategy for generating electron-conductive cytochromes from structurally-defined proteins, and provide design considerations on how heme incorporation may interface with native assembly properties in engineered proteins.

HEMEsPred: Structure-Based Ligand-Specific Heme Binding Residues Prediction by Using Fast-Adaptive Ensemble Learning Scheme

Heme is an essential biomolecule that widely exists in numerous extant organisms. Accurately identifying heme binding residues (HEMEs) is of great importance in disease progression and drug development. In this study, a novel predictor named HEMEsPred was proposed for predicting HEMEs. First, several sequence- and structure-based features, including amino acid composition, motifs, surface preferences, and secondary structure, were collected to construct feature matrices. Second, a novel fast-adaptive ensemble learning scheme was designed to overcome the serious class-imbalance problem as well as to enhance the prediction performance. Third, we further developed ligand-specific models considering that different heme ligands varied significantly in their roles, sizes, and distributions. Statistical test proved the effectiveness of ligand-specific models. Experimental results on benchmark datasets demonstrated good robustness of our proposed method. Furthermore, our method also showed good generalization capability and outperformed many state-of-art predictors on two independent testing datasets. HEMEsPred web server was available at http://www.inforstation.com/HEMEsPred/ for free academic use.

In-Cell Enzymology To Probe His-Heme Ligation in Heme Oxygenase Catalysis

Heme oxygenase (HO) is a ubiquitous enzyme with key roles in inflammation, cell signaling, heme disposal, and iron acquisition. HO catalyzes the oxidative conversion of heme to biliverdin (BV) using a conserved histidine to coordinate the iron atom of bound heme. This His-heme interaction has been regarded as being essential for enzyme activity, because His-to-Ala mutants fail to convert heme to biliverdin in vitro. We probed a panel of proximal His mutants of cyanobacterial, human, and plant HO enzymes using a live-cell activity assay based on heterologous co-expression in Escherichia coli of each HO mutant and a fluorescent biliverdin biosensor. In contrast to in vitro studies with purified proteins, we observed that multiple HO mutants retained significant activity within the intracellular environment of bacteria. X-ray crystallographic structures of human HO1 H25R with bound heme and additional functional studies suggest that HO mutant activity inside these cells does not involve heme ligation by a proximal amino acid. Our study reveals unexpected plasticity in the active site binding interactions with heme that can support HO activity within cells, suggests important contributions by the surrounding active site environment to HO catalysis, and can guide efforts to understand the evolution and divergence of HO function.

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.

Heme-bound SiaA from Streptococcus pyogenes: Effects of mutations and oxidation state on protein stability

The protein SiaA (HtsA) is part of a heme uptake pathway in Streptococcus pyogenes. In this report, we present the heme binding of the alanine mutants of the axial histidine (H229A) and methionine (M79A) ligands, as well as a lysine (K61A) and cysteine (C58A) located near the heme propionates (based on homology modeling) and a control mutant (C47A). pH titrations gave pKa values ranging from 9.0 to 9.5, close to the value of 9.7 for WT SiaA. Resonance Raman spectra of the mutants suggested that the ferric heme environment may be distinct from the wild-type; spectra of the ferrous states were similar. The midpoint reduction potential of the K61A mutant was determined by spectroelectrochemical titration to be 61±3mV vs. SHE, similar to the wild-type protein (68±3mV). The addition of guanidine hydrochloride showed two processes for protein denaturation, consistent with heme loss from protein forms differing by the orientation of the heme in the binding pocket (the half-life for the slower process ranged from less than half a day to two days). The ease of protein unfolding was related to the strength of interaction of the residues with the heme. We hypothesize that kinetically facile but only partial unfolding, followed by a very slow approach to the completely unfolded state, may be a fundamental attribute of heme trafficking proteins. Small motions to release/transfer the heme accompanied by resistance to extensive unfolding may preserve the three dimensional form of the protein for further uptake and release.

Tuning of Hemes b Equilibrium Redox Potential Is Not Required for Cross-Membrane Electron Transfer

In biological energy conversion, cross-membrane electron transfer often involves an assembly of two hemes b. The hemes display a large difference in redox midpoint potentials (ΔE_{m_}b), which in several proteins is assumed to facilitate cross-membrane electron transfer and overcome a barrier of membrane potential. Here we challenge this assumption reporting on heme b ligand mutants of cytochrome bc₁ in which, for the first time in transmembrane cytochrome, one natural histidine has been replaced by lysine without loss of the native low spin type of heme iron. With these mutants we show that ΔE_{m_}b can be markedly increased, and the redox potential of one of the hemes can stay above the level of quinone pool, or ΔE_{m_}b can be markedly decreased to the point that two hemes are almost isopotential, yet the enzyme retains catalytically competent electron transfer between quinone binding sites and remains functional in vivo. This reveals that cytochrome bc₁ can accommodate large changes in ΔE_{m_}b without hampering catalysis, as long as these changes do not impose overly endergonic steps on downhill electron transfer from substrate to product. We propose that hemes b in this cytochrome and in other membranous cytochromes b act as electronic connectors for the catalytic sites with no fine tuning in ΔE_{m_}b required for efficient cross-membrane electron transfer. We link this concept with a natural flexibility in occurrence of several thermodynamic configurations of the direction of electron flow and the direction of the gradient of potential in relation to the vector of the electric membrane potential.

SCMHBP: prediction and analysis of heme binding proteins using propensity scores of dipeptides

Background

Heme binding proteins (HBPs) are metalloproteins that contain a heme ligand (an iron-porphyrin complex) as the prosthetic group. Several computational methods have been proposed to predict heme binding residues and thereby to understand the interactions between heme and its host proteins. However, few in silico methods for identifying HBPs have been proposed.

Results

This work proposes a scoring card method (SCM) based method (named SCMHBP) for predicting and analyzing HBPs from sequences. A balanced dataset of 747 HBPs (selected using a Gene Ontology term GO:0020037) and 747 non-HBPs (selected from 91,414 putative non-HBPs) with an identity of 25% was firstly established. Consequently, a set of scores that quantified the propensity of amino acids and dipeptides to be HBPs is estimated using SCM to maximize the predictive accuracy of SCMHBP. Finally, the informative physicochemical properties of 20 amino acids are identified by utilizing the estimated propensity scores to be used to categorize HBPs. The training and mean test accuracies of SCMHBP applied to three independent test datasets are 85.90% and 71.57%, respectively. SCMHBP performs well relative to comparison with such methods as support vector machine (SVM), decision tree J48, and Bayes classifiers. The putative non-HBPs with high sequence propensity scores are potential HBPs, which can be further validated by experimental confirmation. The propensity scores of individual amino acids and dipeptides are examined to elucidate the interactions between heme and its host proteins. The following characteristics of HBPs are derived from the propensity scores: 1) aromatic side chains are important to the effectiveness of specific HBP functions; 2) a hydrophobic environment is important in the interaction between heme and binding sites; and 3) the whole HBP has low flexibility whereas the heme binding residues are relatively flexible.

Conclusions

SCMHBP yields knowledge that improves our understanding of HBPs rather than merely improves the prediction accuracy in predicting HBPs.

Redox state dependence of axial ligand dynamics in Nitrosomonas europaea cytochrome c552

Analysis of NMR spectra reveals that the heme axial Met ligand orientation and dynamics in Nitrosomonas europaea cytochrome c552 (Ne cyt c) are dependent on the heme redox state. In the oxidized state, the heme axial Met is fluxional, interconverting between two conformers related to each other by inversion through the Met δS atom. In the reduced state, there is no evidence of fluxionality, with the Met occupying one conformation similar to that seen in the homologous Pseudomonas aeruginosa cytochrome c551. Comparison of the observed and calculated pseudocontact shifts for oxidized Ne cyt c using the reduced protein structure as a reference structure reveals a redox-dependent change in the structure of the loop bearing the axial Met (loop 3). Analysis of nuclear Overhauser effects (NOEs) and existing structural data provides further support for the redox state dependence of the loop 3 structure. Implications for electron transfer function are discussed.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Direct tests of enzymatic heme degradation by the malaria parasite Plasmodium falciparum

Malaria parasites generate vast quantities of heme during blood stage infection via hemoglobin digestion and limited de novo biosynthesis, but it remains unclear if parasites metabolize heme for utilization or disposal. Recent in vitro experiments with a heme oxygenase (HO)-like protein from Plasmodium falciparum suggested that parasites may enzymatically degrade some heme to the canonical HO product, biliverdin (BV), or its downstream metabolite, bilirubin (BR). To directly test for BV and BR production by P. falciparum parasites, we DMSO-extracted equal numbers of infected and uninfected erythrocytes and developed a sensitive LC-MS/MS assay to quantify these tetrapyrroles. We found comparable low levels of BV and BR in both samples, suggesting the absence of HO activity in parasites. We further tested live parasites by targeted expression of a fluorescent BV-binding protein within the parasite cytosol, mitochondrion, and plant-like plastid. This probe could detect exogenously added BV but gave no signal indicative of endogenous BV production within parasites. Finally, we recombinantly expressed and tested the proposed heme degrading activity of the HO-like protein, PfHO. Although PfHO bound heme and protoporphyrin IX with modest affinity, it did not catalyze heme degradation in vivo within bacteria or in vitro in UV absorbance and HPLC assays. These observations are consistent with PfHO's lack of a heme-coordinating His residue and suggest an alternative function within parasites. We conclude that P. falciparum parasites lack a canonical HO pathway for heme degradation and thus rely fully on alternative mechanisms for heme detoxification and iron acquisition during blood stage infection.

Prediction of heme binding residues from protein sequences with integrative sequence profiles

Background

The heme-protein interactions are essential for various biological processes such as electron transfer, catalysis, signal transduction and the control of gene expression. The knowledge of heme binding residues can provide crucial clues to understand these activities and aid in functional annotation, however, insufficient work has been done on the research of heme binding residues from protein sequence information.

Methods

We propose a sequence-based approach for accurate prediction of heme binding residues by a novel integrative sequence profile coupling position specific scoring matrices with heme specific physicochemical properties. In order to select the informative physicochemical properties, we design an intuitive feature selection scheme by combining a greedy strategy with correlation analysis.

Results

Our integrative sequence profile approach for prediction of heme binding residues outperforms the conventional methods using amino acid and evolutionary information on the 5-fold cross validation and the independent tests.

Conclusions

The novel feature of an integrative sequence profile achieves good performance using a reduced set of feature vector elements.

Heme cytotoxicity and the pathogenesis of immune-mediated inflammatory diseases

Heme, iron (Fe) protoporphyrin IX, functions as a prosthetic group in a range of hemoproteins essential to support life under aerobic conditions. The Fe contained within the prosthetic heme groups of these hemoproteins can catalyze the production of reactive oxygen species. Presumably for this reason, heme must be sequestered within those hemoproteins, thereby shielding the reactivity of its Fe-heme. However, under pathologic conditions associated with oxidative stress, some hemoproteins can release their prosthetic heme groups. While this heme is not necessarily damaging per se, it becomes highly cytotoxic in the presence of a range of inflammatory mediators such as tumor necrosis factor. This can lead to tissue damage and, as such, exacerbate the pathologic outcome of several immune-mediated inflammatory conditions. Presumably, targeting “free heme” may be used as a therapeutic intervention against these diseases.

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.

Computational prediction of heme-binding residues by exploiting residue interaction network

Computational identification of heme-binding residues is beneficial for predicting and designing novel heme proteins. Here we proposed a novel method for heme-binding residue prediction by exploiting topological properties of these residues in the residue interaction networks derived from three-dimensional structures. Comprehensive analysis showed that key residues located in heme-binding regions are generally associated with the nodes with higher degree, closeness and betweenness, but lower clustering coefficient in the network. HemeNet, a support vector machine (SVM) based predictor, was developed to identify heme-binding residues by combining topological features with existing sequence and structural features. The results showed that incorporation of network-based features significantly improved the prediction performance. We also compared the residue interaction networks of heme proteins before and after heme binding and found that the topological features can well characterize the heme-binding sites of apo structures as well as those of holo structures, which led to reliable performance improvement as we applied HemeNet to predicting the binding residues of proteins in the heme-free state. HemeNet web server is freely accessible at http://mleg.cse.sc.edu/hemeNet/.

Design principles for chlorophyll-binding sites in helical proteins

The cyclic tetrapyrroles, viz. chlorophylls (Chl), their bacterial analogs bacteriochlorophylls, and hemes are ubiquitous cofactors of biological catalysis that are involved in a multitude of reactions. One systematic approach for understanding how Nature achieves functional diversity with only this handful of cofactors is by designing de novo simple and robust protein scaffolds with heme and/or (bacterio)chlorophyll [(B)Chls]-binding sites. This strategy is currently mostly implemented for heme-binding proteins. To gain more insight into the factors that determine heme-/(B)Chl-binding selectivity, we explored the geometric parameters of (B)Chl-binding sites in a nonredundant subset of natural (B)Chl protein structures. Comparing our analysis to the study of a nonredundant database of heme-binding helical histidines by Negron et al. (Proteins 2009;74:400-416), we found a preference for the m-rotamer in (B)Chl-binding helical histidines, in contrast to the preferred t-rotamer in heme-binding helical histidines. This may be used for the design of specific heme- or (B)Chl-binding sites in water-soluble helical bundles, because the rotamer type defines the positioning of the bound cofactor with respect to the helix interface and thus the protein-binding site. Consensus sequences for (B)Chl binding were identified by combining a computational and database-derived approach and shown to be significantly different from the consensus sequences recommended by Negron et al. (Proteins 2009;74:400-416) for heme-binding helical proteins. The insights gained in this work on helix- (B)Chls-binding pockets provide useful guidelines for the construction of reasonable (B)Chl-binding protein templates that can be optimized by computational tools.

HemeBIND: a novel method for heme binding residue prediction by combining structural and sequence information

BACKGROUND

Accurate prediction of binding residues involved in the interactions between proteins and small ligands is one of the major challenges in structural bioinformatics. Heme is an essential and commonly used ligand that plays critical roles in electron transfer, catalysis, signal transduction and gene expression. Although much effort has been devoted to the development of various generic algorithms for ligand binding site prediction over the last decade, no algorithm has been specifically designed to complement experimental techniques for identification of heme binding residues. Consequently, an urgent need is to develop a computational method for recognizing these important residues.

RESULTS

Here we introduced an efficient algorithm HemeBIND for predicting heme binding residues by integrating structural and sequence information. We systematically investigated the characteristics of binding interfaces based on a non-redundant dataset of heme-protein complexes. It was found that several sequence and structural attributes such as evolutionary conservation, solvent accessibility, depth and protrusion clearly illustrate the differences between heme binding and non-binding residues. These features can then be separately used or combined to build the structure-based classifiers using support vector machine (SVM). The results showed that the information contained in these features is largely complementary and their combination achieved the best performance. To further improve the performance, an attempt has been made to develop a post-processing procedure to reduce the number of false positives. In addition, we built a sequence-based classifier based on SVM and sequence profile as an alternative when only sequence information can be used. Finally, we employed a voting method to combine the outputs of structure-based and sequence-based classifiers, which demonstrated remarkably better performance than the individual classifier alone.

CONCLUSIONS

HemeBIND is the first specialized algorithm used to predict binding residues in protein structures for heme ligands. Extensive experiments indicated that both the structure-based and sequence-based methods have effectively identified heme binding residues while the complementary relationship between them can result in a significant improvement in prediction performance. The value of our method is highlighted through the development of HemeBIND web server that is freely accessible at http://mleg.cse.sc.edu/hemeBIND/.

Crystal structures of CO and NO adducts of MauG in complex with pre-methylamine dehydrogenase: implications for the mechanism of dioxygen activation

MauG is a diheme enzyme responsible for the post-translational formation of the catalytic tryptophan tryptophylquinone (TTQ) cofactor in methylamine dehydrogenase (MADH). MauG can utilize hydrogen peroxide, or molecular oxygen and reducing equivalents, to complete this reaction via a catalytic bis-Fe(IV) intermediate. Crystal structures of diferrous, Fe(II)-CO, and Fe(II)-NO forms of MauG in complex with its preMADH substrate have been determined and compared to one another as well as to the structure of the resting diferric MauG-preMADH complex. CO and NO each bind exclusively to the 5-coordinate high-spin heme with no change in ligation of the 6-coordinate low-spin heme. These structures reveal likely roles for amino acid residues in the distal pocket of the high-spin heme in oxygen binding and activation. Glu113 is implicated in the protonation of heme-bound diatomic oxygen intermediates in promoting cleavage of the O-O bond. Pro107 is shown to change conformation on the binding of each ligand and may play a steric role in oxygen activation by positioning the distal oxygen near Glu113. Gln103 is in a position to provide a hydrogen bond to the Fe(IV)═O moiety that may account for the unusual stability of this species in MauG.

Heme proteins--diversity in structural characteristics, function, and folding

The characteristics of heme prosthetic groups and their binding sites have been analyzed in detail in a data set of nonhomologous heme proteins. Variations in the shape, volume, and chemical composition of the binding site, in the mode of heme binding and in the number and nature of heme-protein interactions are found to result in significantly different heme environments in proteins with different functions in biology. Differences are also seen in the properties of the apo states of the proteins. The apo states of proteins that bind heme permanently in their functional form show some disorder, ranging from local unfolding in the heme binding pocket to complete unfolding to give a random coil. In contrast, proteins that bind heme transiently are fully folded in their apo and holo states, presumably allowing both apo and holo forms to remain biologically active resisting aggregation or proteolysis. The principles identified here provide a framework for the design of de novo proteins that will exhibit tight heme ligand binding and for the identification of the function of structural genomic target proteins with heme ligands.

Multiconformation continuum electrostatics analysis of the effects of a buried Asp introduced near heme a in Rhodobacter sphaeroides cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces O(2) to water via a series of proton-coupled electron transfers, generating a transmembrane electrochemical gradient. Coupling electron and proton transfer requires changing the pK(a) values of buried residues at each stage in the reaction cycle. Heme a is a key cofactor in the CcO electron transfer chain. Mutation of Ser44 to Asp has been reported [Mills, D. A., et al. (2008) Biochemistry 47, 11499-11509], changing the hydrogen bond acceptor from His102, the heme a axial ligand in Rhodobactor sphaeroides CcO. This adds an acidic residue to the CcO interior. The electrochemical behavior of heme a in wild-type and S44D CcO is compared using the continuum electrostatics program MCCE. The introduced, deeply buried Asp remains ionized at physiological pH only when the nearby heme is oxidized. Heme a reduction is now calculated to be strongly coupled to Asp proton binding, while with Ser44, it is weakly coupled to small protonation shifts at multiple sites, increasing the pH dependence in the mutant. At pH 7, the partially ionized Asp 44 is calculated to lower the heme redox potential by 50 mV as expected given the thermodynamics of coupled electron and proton transfers. This highlights an curious finding in the experimental results where a low Asp pK(a) is found together with a stabilized reduced heme. The stabilization of a heme oxidation in a model complex by a hydrogen bond to the axial His ligand calculated with continuum electrostatics and with density functional theory were in good agreement.

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

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

Understanding properties of cofactors in proteins: redox potentials of synthetic cytochromes b

Haehnel et al. synthesized 399 different artificial cytochrome b (aCb) models. They consist of a template-assisted four-helix bundle with one embedded heme group. Their redox potentials were measured and cover the range from -148 to -89 mV. No crystal structures of these aCb are available. Therefore, we use the chemical composition and general structural principles to generate atomic coordinates of 31 of these aCb mutants, which are chosen to cover the whole interval of redox potentials. We start by modeling the coordinates of one aCb from scratch. Its structure remains stable after energy minimization and during molecular dynamics simulation over 2 ns. Based on this structure, coordinates of the other 30 aCb mutants are modeled. The calculated redox potentials for these 31 aCb agree within 10 mV with the experimental values in terms of root mean square deviation. Analysis of the dependence of heme redox potential on protein environment shows that the shifts in redox potentials relative to the model systems in water are due to the low-dielectric medium of the protein and the protonation states of the heme propionic acid groups, which are influenced by the surrounding amino acids. Alternatively, we perform a blind prediction of the same redox potentials using an empirical approach based on a linear scoring function and reach a similar accuracy. Both methods are useful to understand and predict heme redox potentials. Based on the modeled structure we can understand the detailed structural differences between aCb mutants that give rise to shifts in heme redox potential. On the other hand, one can explore the correlation between sequence variations and aCb redox potentials more directly and on much larger scale using the empirical prediction scheme, which--thanks to its simplicity--is much faster.

Analysis of the electrochemistry of hemes with E(m)s spanning 800 mV

The free energy of heme reduction in different proteins is found to vary over more than 18 kcal/mol. It is a challenge to determine how proteins manage to achieve this enormous range of E(m)s with a single type of redox cofactor. Proteins containing 141 unique hemes of a-, b-, and c-type, with bis-His, His-Met, and aquo-His ligation were calculated using Multi-Conformation Continuum Electrostatics (MCCE). The experimental E(m)s range over 800 mV from -350 mV in cytochrome c(3) to 450 mV in cytochrome c peroxidase (vs. SHE). The quantitative analysis of the factors that modulate heme electrochemistry includes the interactions of the heme with its ligands, the solvent, the protein backbone, and sidechains. MCCE calculated E(m)s are in good agreement with measured values. Using no free parameters the slope of the line comparing calculated and experimental E(m)s is 0.73 (R(2) = 0.90), showing the method accounts for 73% of the observed E(m) range. Adding a +160 mV correction to the His-Met c-type hemes yields a slope of 0.97 (R(2) = 0.93). With the correction 65% of the hemes have an absolute error smaller than 60 mV and 92% are within 120 mV. The overview of heme proteins with known structures and E(m)s shows both the lowest and highest potential hemes are c-type, whereas the b-type hemes are found in the middle E(m) range. In solution, bis-His ligation lowers the E(m) by approximately 205 mV relative to hemes with His-Met ligands. The bis-His, aquo-His, and His-Met ligated b-type hemes all cluster about E(m)s which are approximately 200 mV more positive in protein than in water. In contrast, the low potential bis-His c-type hemes are shifted little from in solution, whereas the high potential His-Met c-type hemes are raised by approximately 300 mV from solution. The analysis shows that no single type of interaction can be identified as the most important in setting heme electrochemistry in proteins. For example, the loss of solvation (reaction field) energy, which raises the E(m), has been suggested to be a major factor in tuning in situ E(m)s. However, the calculated solvation energy vs. experimental E(m) shows a slope of 0.2 and R(2) of 0.5 thus correlates weakly with E(m)s. All other individual interactions show even less correlation with E(m). However the sum of these terms does reproduce the range of observed E(m)s. Therefore, different proteins use different aspects of their structures to modulate the in situ heme electrochemistry. This study also shows that the calculated E(m)s are relatively insensitive to different heme partial charges and to the protein dielectric constant used in the simulation.

p3d--Python module for structural bioinformatics

Background

High-throughput bioinformatic analysis tools are needed to mine the large amount of structural data via knowledge based approaches. The development of such tools requires a robust interface to access the structural data in an easy way. For this the Python scripting language is the optimal choice since its philosophy is to write an understandable source code.

Results

p3d is an object oriented Python module that adds a simple yet powerful interface to the Python interpreter to process and analyse three dimensional protein structure files (PDB files). p3d's strength arises from the combination of a) very fast spatial access to the structural data due to the implementation of a binary space partitioning (BSP) tree, b) set theory and c) functions that allow to combine a and b and that use human readable language in the search queries rather than complex computer language. All these factors combined facilitate the rapid development of bioinformatic tools that can perform quick and complex analyses of protein structures.

Conclusion

p3d is the perfect tool to quickly develop tools for structural bioinformatics using the Python scripting language.