Distal regulation of heme binding of heme oxygenase-1 mediated by conformational fluctuations

Comparison of the Backbone Dynamics of Dehaloperoxidase-Hemoglobin Isoenzymes

Dehaloperoxidase (DHP) is a multifunctional hemeprotein with a functional switch generally regulated by the chemical class of the substrate. Its two isoforms, DHP-A and DHP-B, differ by only five amino acids and have an almost identical protein fold. However, the catalytic efficiency of DHP-B for oxidation by a peroxidase mechanism ranges from 2- to 6-fold greater than that of DHP-A depending on the conditions. X-ray crystallography has shown that many substrates and ligands have nearly identical binding in the two isoenzymes, suggesting that the difference in catalytic efficiency could be due to differences in the conformational dynamics. We compared the backbone dynamics of the DHP isoenzymes at pH 7 through heteronuclear relaxation dynamics at 11.75, 16.45, and 19.97 T in combination with four 300 ns MD simulations. While the overall dynamics of the isoenzymes are similar, there are specific local differences in functional regions of each protein. In DHP-A, Phe35 undergoes a slow chemical exchange between two conformational states likely coupled to a swinging motion of Tyr34. Moreover, Asn37 undergoes fast chemical exchange in DHP-A. Given that Phe35 and Asn37 are adjacent to Tyr34 and Tyr38, it is possible that their dynamics modulate the formation and migration of the active tyrosyl radicals in DHP-A at pH 7. Another significant difference is that both distal and proximal histidines have a 15-18% smaller S2 value in DHP-B, thus their greater flexibility could account for the higher catalytic activity. The distal histidine grants substrate access to the distal pocket. The greater flexibility of the proximal histidine could also accelerate H2O2 activation at the heme Fe by increased coupling of an amino acid charge relay to stabilize the ferryl Fe(IV) oxidation state in a Poulos-Kraut "push-pull"-type peroxidase mechanism.

Crystal structure of higher plant heme oxygenase-1 and its mechanism of interaction with ferredoxin

Heme oxygenase (HO) converts heme to carbon monoxide, biliverdin, and free iron, products that are essential in cellular redox signaling and iron recycling. In higher plants, HO is also involved in the biosynthesis of photoreceptor pigment precursors. Despite many common enzymatic reactions, the amino acid sequence identity between plant-type and other HOs is exceptionally low (∼19.5%), and amino acids that are catalytically important in mammalian HO are not conserved in plant-type HOs. Structural characterization of plant-type HO is limited to spectroscopic characterization by electron spin resonance, and it remains unclear how the structure of plant-type HO differs from that of other HOs. Here, we have solved the crystal structure of Glycine max (soybean) HO-1 (GmHO-1) at a resolution of 1.06 Å and carried out the isothermal titration calorimetry measurements and NMR spectroscopic studies of its interaction with ferredoxin, the plant-specific electron donor. The high-resolution X-ray structure of GmHO-1 reveals several novel structural components: an additional irregularly structured region, a new water tunnel from the active site to the surface, and a hydrogen-bonding network unique to plant-type HOs. Structurally important features in other HOs, such as His ligation to the bound heme, are conserved in GmHO-1. Based on combined data from X-ray crystallography, isothermal titration calorimetry, and NMR measurements, we propose the evolutionary fine-tuning of plant-type HOs for ferredoxin dependency in order to allow adaptation to dynamic pH changes on the stroma side of the thylakoid membrane in chloroplast without losing enzymatic activity under conditions of fluctuating light.

Recent Advances in the Understanding of the Reaction Chemistries of the Heme Catabolizing Enzymes HO and BVR Based on High Resolution Protein Structures

In mammals, catabolism of the heme group is indispensable for life. Heme is first cleaved by the enzyme Heme Oxygenase (HO) to the linear tetrapyrrole Biliverdin IXα (BV), and BV is then converted into bilirubin by Biliverdin Reductase (BVR). HO utilizes three Oxygen molecules (O2) and seven electrons supplied by NADPH-cytochrome P450 oxidoreductase (CPR) to open the heme ring and BVR reduces BV through the use of NAD(P)H. Structural studies of HOs, including substrate-bound, reaction intermediate-bound, and several specific inhibitor-bound forms, reveal details explaining substrate binding to HO and mechanisms underlying-specific HO reaction progression. Cryo-trapped structures and a time-resolved spectroscopic study examining photolysis of the bond between the distal ligand and heme iron demonstrate how CO, produced during the HO reaction, dissociates from the reaction site with a corresponding conformational change in HO. The complex structure containing HO and CPR provides details of how electrons are transferred to the heme-HO complex. Although the tertiary structure of BVR and its complex with NAD+ was determined more than 10 years ago, the catalytic residues and the reaction mechanism of BVR remain unknown. A recent crystallographic study examining cyanobacterial BVR in complex with NADP+ and substrate BV provided some clarification regarding these issues. Two BV molecules are bound to BVR in a stacked manner, and one BV may assist in the reductive catalysis of the other BV. In this review, recent advances illustrated by biochemical, spectroscopic, and crystallographic studies detailing the chemistry underlying the molecular mechanism of HO and BVR reactions are presented.

Structural, Photophysical, and Photochemical Characterization of Zinc Protoporphyrin IX in a Dimeric Variant of an Iron Storage Protein: Insights into the Mechanism of Photosensitized H2 Generation

Some of us have previously reported the preparation of a dimeric form of the iron storage protein, bacterioferritin (Bfr), in which the native heme b is substituted with the photosensitizer, Zn(II)-protoporphyrin IX (ZnPP-Bfr dimer). We further showed that the ZnPP-Bfr dimer can serve as a photosensitizer for platinum-catalyzed H₂ generation in aqueous solution without the usually added electron relay between photosensitizer and platinum ( Clark , E. R. , Inorg. Chem. 2017 , 56 , 4584 - 4593 ). We proposed reductive or oxidative quenching pathways involving the ZnPP anion radical (ZnPP^•-) or the ZnPP cation radical, (ZnPP^•+), respectively. The present report describes structural, photophysical, and photochemical properties of the ZnPP in the ZnPP-Bfr dimer. X-ray absorption spectroscopic studies at 10 K showed a mixture of five- and six-coordinated Zn centers with axial coordination by one long Zn-SγMet distance of ∼2.8 Å and ∼40% having an additional shorter Zn-S distance of ∼2.4 Å, in addition to the expected 4 nitrogen atom coordination from the porphyrin. The ZnPP in ZnPP-Bfr dimer was prone to photosensitized oxidation to ZnPP^•+. The ZnPP^•+ was rapidly reduced by ascorbic acid, which we previously determined was essential for photosensitized H₂ production in this system. These results are consistent with an oxidative quenching pathway involving electron transfer from ³ZnPP* to platinum, which may be assisted by a flexible ZnPP axial coordination sphere. However, the low quantum yield for H₂ production (∼1%) in this system could make reductive quenching difficult to detect, and can, therefore, not be completely ruled out. The ZnPP-Bfr dimer provides a simple but versatile framework for mechanistic assessment and optimization of porphyrin-photosensitized H₂ generation without an electron relay between porphyrin and the platinum catalyst.

Overview of Relaxation Dispersion NMR Spectroscopy to Study Protein Dynamics and Protein-Ligand Interactions

Proteins and nucleic acids are central to all biological processes. NMR spectroscopy has proven to be excellent for studying the dynamics of these macromolecules over various timescales. Relaxation rates and heteronuclear nuclear Overhauser-effect values can resolve motion on pico- to nanosecond timescales, residual dipolar couplings provide information on submicro- to millisecond timescales, and even slower dynamics over seconds to hours can be resolved by hydrogen-exchange experiments. Relaxation dispersion experiments are especially valuable because they resolve motion on micro- to millisecond timescales, encompassing biomolecular motions associated with ligand binding, enzymatic catalysis, and domain-domain opening. These experiments provide structural, kinetic, and thermodynamic information on "invisible" excited conformational states. Relaxation dispersion can be applied not only to single biomolecules but also to protein-ligand complexes to study the kinetics and thermodynamics of association and dissociation. We review recent developments in relaxation dispersion methodology, outline the R_1ρ relaxation dispersion experiment, and discuss application to biomolecular interactions. © 2018 by John Wiley & Sons, Inc.

Prediction of the interaction of metallic moieties with proteins: An update for protein-ligand docking techniques

In this article, we present a new approach to expand the range of application of protein-ligand docking methods in the prediction of the interaction of coordination complexes (i.e., metallodrugs, natural and artificial cofactors, etc.) with proteins. To do so, we assume that, from a pure computational point of view, hydrogen bond functions could be an adequate model for the coordination bonds as both share directionality and polarity aspects. In this model, docking of metalloligands can be performed without using any geometrical constraints or energy restraints. The hard work consists in generating the convenient atom types and scoring functions. To test this approach, we applied our model to 39 high-quality X-ray structures with transition and main group metal complexes bound via a unique coordination bond to a protein. This concept was implemented in the protein-ligand docking program GOLD. The results are in very good agreement with the experimental structures: the percentage for which the RMSD of the simulated pose is smaller than the X-ray spectra resolution is 92.3% and the mean value of RMSD is < 1.0 Å. Such results also show the viability of the method to predict metal complexes-proteins interactions when the X-ray structure is not available. This work could be the first step for novel applicability of docking techniques in medicinal and bioinorganic chemistry and appears generalizable enough to be implemented in most protein-ligand docking programs nowadays available. © 2017 Wiley Periodicals, Inc.

Chronic obstructive pulmonary disease candidate gene prioritization based on metabolic networks and functional information

Chronic obstructive pulmonary disease (COPD) is a multi-factor disease, in which metabolic disturbances played important roles. In this paper, functional information was integrated into a COPD-related metabolic network to assess similarity between genes. Then a gene prioritization method was applied to the COPD-related metabolic network to prioritize COPD candidate genes. The gene prioritization method was superior to ToppGene and ToppNet in both literature validation and functional enrichment analysis. Top-ranked genes prioritized from the metabolic perspective with functional information could promote the better understanding about the molecular mechanism of this disease. Top 100 genes might be potential markers for diagnostic and effective therapies.

Elucidation of potential sites for antibody engineering by fluctuation editing

Target-specific monoclonal antibodies can be routinely acquired, but the sequences of naturally acquired antibodies are not always affinity-matured and methods that increase antigen affinity are desirable. Most biophysical studies have focused on the complementary determining region (CDR), which directly contacts the antigen; however, it remains difficult to increase the affinity as much as desired. While strategies to alter the CDR to increase antibody affinity are abundant, those that target non-CDR regions are scarce. Here we describe a new method, designated fluctuation editing, which identifies potential mutation sites and engineers a high-affinity antibody based on conformational fluctuations observed by NMR relaxation dispersion. Our data show that relaxation dispersion detects important fluctuating residues that are not located in the CDR and that increase antigen–antibody affinity by point mutation. The affinity-increased mutants are shown to fluctuate less in their free form and to form a more packed structure in their antigen-bound form.

Structural and mutational analyses of the Leptospira interrogans virulence-related heme oxygenase provide insights into its catalytic mechanism

Heme oxygenase from Leptospira interrogans is an important virulence factor. During catalysis, redox equivalents are provided to this enzyme by the plastidic-type ferredoxin-NADP⁺ reductase also found in L. interrogans. This process may have evolved to aid this bacterial pathogen to obtain heme-iron from their host and enable successful colonization. Herein we report the crystal structure of the heme oxygenase-heme complex at 1.73 Å resolution. The structure reveals several distinctive features related to its function. A hydrogen bonded network of structural water molecules that extends from the catalytic site to the protein surface was cleared observed. A depression on the surface appears to be the H⁺ network entrance from the aqueous environment to the catalytic site for O₂ activation, a key step in the heme oxygenase reaction. We have performed a mutational analysis of the F157, located at the above-mentioned depression. The mutant enzymes were unable to carry out the complete degradation of heme to biliverdin since the reaction was arrested at the verdoheme stage. We also observed that the stability of the oxyferrous complex, the efficiency of heme hydroxylation and the subsequent conversion to verdoheme was adversely affected. These findings underscore a long-range communication between the outer fringes of the hydrogen-bonded network of structural waters and the heme active site during catalysis. Finally, by analyzing the crystal structures of ferredoxin-NADP⁺ reductase and heme oxygenase, we propose a model for the productive association of these proteins.

Quantitative analysis of protein-ligand interactions by NMR

Protein-ligand interactions have been commonly studied through static structures of the protein-ligand complex. Recently, however, there has been increasing interest in investigating the dynamics of protein-ligand interactions both for fundamental understanding of the underlying mechanisms and for drug development. NMR is a versatile and powerful tool, especially because it provides site-specific quantitative information. NMR has widely been used to determine the dissociation constant (KD), in particular, for relatively weak interactions. The simplest NMR method is a chemical-shift titration experiment, in which the chemical-shift changes of a protein in response to ligand titration are measured. There are other quantitative NMR methods, but they mostly apply only to interactions in the fast-exchange regime. These methods derive the dissociation constant from population-averaged NMR quantities of the free and bound states of a protein or ligand. In contrast, the recent advent of new relaxation-based experiments, including R2 relaxation dispersion and ZZ-exchange, has enabled us to obtain kinetic information on protein-ligand interactions in the intermediate- and slow-exchange regimes. Based on R2 dispersion or ZZ-exchange, methods that can determine the association rate, kon, dissociation rate, koff, and KD have been developed. In these approaches, R2 dispersion or ZZ-exchange curves are measured for multiple samples with different protein and/or ligand concentration ratios, and the relaxation data are fitted to theoretical kinetic models. It is critical to choose an appropriate kinetic model, such as the two- or three-state exchange model, to derive the correct kinetic information. The R2 dispersion and ZZ-exchange methods are suitable for the analysis of protein-ligand interactions with a micromolar or sub-micromolar dissociation constant but not for very weak interactions, which are typical in very fast exchange. This contrasts with the NMR methods that are used to analyze population-averaged NMR quantities. Essentially, to apply NMR successfully, both the type of experiment and equation to fit the data must be carefully and specifically chosen for the protein-ligand interaction under analysis. In this review, we first explain the exchange regimes and kinetic models of protein-ligand interactions, and then describe the NMR methods that quantitatively analyze these specific interactions.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

Extracting protein dynamics information from overlapped NMR signals using relaxation dispersion difference NMR spectroscopy

Protein dynamics plays important roles in many biological events, such as ligand binding and enzyme reactions. NMR is mostly used for investigating such protein dynamics in a site-specific manner. Recently, NMR has been actively applied to large proteins and intrinsically disordered proteins, which are attractive research targets. However, signal overlap, which is often observed for such proteins, hampers accurate analysis of NMR data. In this study, we have developed a new methodology called relaxation dispersion difference that can extract conformational exchange parameters from overlapped NMR signals measured using relaxation dispersion spectroscopy. In relaxation dispersion measurements, the signal intensities of fluctuating residues vary according to the Carr-Purcell-Meiboon-Gill pulsing interval, whereas those of non-fluctuating residues are constant. Therefore, subtraction of each relaxation dispersion spectrum from that with the highest signal intensities, measured at the shortest pulsing interval, leaves only the signals of the fluctuating residues. This is the principle of the relaxation dispersion difference method. This new method enabled us to extract exchange parameters from overlapped signals of heme oxygenase-1, which is a relatively large protein. The results indicate that the structural flexibility of a kink in the heme-binding site is important for efficient heme binding. Relaxation dispersion difference requires neither selectively labeled samples nor modification of pulse programs; thus it will have wide applications in protein dynamics analysis.

Utility of heme analogues to intentionally modify heme-globin interactions in myoglobin

Myoglobin reconstitution with various synthetic heme analogues was reviewed to follow the consequences of modified heme-globin interactions. Utility of dimethyl sulfoxide as the solvent for water-insoluble hemes was emphasized. Proton NMR spectroscopy revealed that loose heme-globin contacts in the heme pocket eventually caused the dynamic heme rotation around the iron-histidine bond. The full rotational rate was estimated to be about 1400 s(-1) at room temperature for 1,4,5,8-tetramethylhemin. The X-ray analysis of the myoglobin containing iron porphine, the smallest heme without any side chains, showed that the original globin fold was well conserved despite the serious disruption of native heme-globin contacts. Comparison between the two myoglobins with static and rotatory prosthetic groups indicated that the oxygen and carbon monoxide binding profiles were almost unaffected by the heme motion. On the other hand, altered tetrapyrrole array of porphyrin dramatically changed the dissociation constant of oxygen from 0.0005 mm Hg of porphycene-myoglobin to ∞ in oxypyriporphyrin-myoglobin. Heme-globin interactions in myoglobin were also monitored with circular dichroism spectroscopy. The observation on several reconstituted protein revealed an unrecognized role of the propionate groups in protoheme. Shortening of heme 6,7-propionates to carboxylates resulted in almost complete disappearance of the positive circular dichroism band in the Soret region. The theoretical analysis suggested that the disappeared circular dichroism band reflected the cancellation effects between different conformers of the carboxyl groups directly attached to heme periphery. The above techniques were proposed to be applicable to other hemoproteins to create new biocatalysts. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.