Engineered chimeras reveal the structural basis of hexacoordination in globins: A case study of neuroglobin and myoglobin

Background

Myoglobin (Mb) and neuroglobin (Ngb) are representative members of pentacoordinated and bis-histidyl, hexacoordinated globins. In spite of their low sequence identity, they show surprisingly similar three-dimensional folds. The ability of Ngb to form a hexacoordinated bis-histidyl complex with the distal HisE7 has a strong impact on ligand affinity. The factors governing such different behaviors have not been completely understood yet, even though they are extremely relevant to establish structure–function relationships within the globin superfamily.

Methods

In this work we generated chimeric proteins by swapping a previously identified regulatory segment – the CD region – and evaluated comparatively the structural and functional properties of the resulting proteins by molecular dynamics simulations, and spectroscopic and kinetic investigations.

Results

Our results show that chimeric proteins display heme coordination properties displaced towards those expected for the corresponding CD region. In particular, in the absence of exogenous ligands, chimeric Mb is found as a partially hexacoordinated bis-histidyl species, whereas chimeric Ngb shows a lower equilibrium constant for forming a hexacoordinated bis-histidyl species.

Conclusions

While these results confirm the regulatory role of the CD region for bis-histidyl hexacoordination, they also suggest that additional sources contribute to fine tune the equilibrium.General significanceGlobins constitute a ubiquitous group of heme proteins widely found in all kingdoms of life. These findings raise challenging questions regarding the structure–function relationships in these proteins, as bis-histidyl hexacoordination emerges as a novel regulatory mechanism of the physiological function of globins.     

X-Ray crystallographic structural characteristics of Arabidopsis hemoglobin I and their functional implications

Genome of the model dicot flowering plant, Arabidopsis thaliana, a popular tool for understanding molecular biology of plant physiology, encodes all three classes of plant hemoglobins that differ in their sequence, ligand binding and spectral properties. As such these globins are of considerable attention. Crystal structures of few members of plant class I nonsymbiotic hemoglobin have been described earlier. Here we report the crystal structure of Arabidopsis class I hemoglobin (AHb1) to 2.2 ? and compare its key features with the structures of similar nonsymbiotic hemoglobin from other species. Crystal structure of AHb1 is homologous to the related members with similar globin fold and heme pocket architecture. The structure is homodimeric in the asymmetric unit with both distal and proximal histidines coordinating to the heme iron atom. Residues lining the dimeric interface are also conserved in AHb1 with the exception of additional electrostatic interaction between H112 and E113 of each subunit and that involving Y119 through two water molecules. In addition, differences in heme pocket non-covalent interactions, a novel Ser residue at F7 position, Xe binding site variability, internal cavity topology differences, CD loop conformation and stability and other such properties might explain kinetic variability in AHb1. Detailed cavity analysis of AHb1 showed the presence of a novel long tunnel connecting the distal pockets of both the monomers. Presence of such tunnel, along with conformational heterogeneity observed in the two chains, might suggest cooperative ligand binding and support its role in NO scavenging. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Different roles of protein dynamics and ligand migration in non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana

The ligand rebinding kinetics after photolysis of the CO complexes of Arabidopsis thaliana hemoglobins AHb1 and AHb2 in solution show very different amplitudes in the geminate phase, reflecting different migration pathways of the photodissociated ligand in the system of internal cavities accessible from the heme. The dependence of the geminate phase on CO concentration, temperature, encapsulation in silica gels and presence of glycerol confirms a remarkable difference in the internal structure of the two proteins and a dramatically different role of protein dynamics in regulating the reactivity with CO. This finding strongly supports the idea that they have distinct physiological functions.     

The quantum mechanically mixed-spin state in a non-symbiotic plant hemoglobin: the effect of distal mutation on AHb1 from Arabidopsis thaliana

Non-symbiotic hemoglobins are hexacoordinated heme proteins found in all plants. To gain insight into the importance of the heme hexacoordination and the coordinated distal histidine in general for the possible physiological functions of these proteins, the distal His(E7) of Arabidopsis thaliana hemoglobin (AHb1) was substituted by a leucine residue. The heme properties of the wild-type and mutant proteins have been characterized by electronic absorption, resonance Raman and electron paramagnetic resonance spectroscopic studies at room and low temperatures. Significant differences between the wild-type and mutant proteins have been detected. The most striking is the formation of an uncommon quantum mechanically mixed-spin heme species in the mutant. This is the first observation of such a spin state in a plant hemoglobin. The proportion of this species, which at room temperature coexists with a minor pentacoordinated high-spin form, increases markedly at low temperature.     

Oxygen and nitric oxide rebinding kinetics in nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana

Type 1 nonsymbiotic hemoglobin from Arabidopsis thaliana (AHb1) shows a partial bis-histidyl hexacoordination but can reversibly bind diatomic ligands. The physiological function is still unclear, but the high oxygen affinity rules out a function related to O2 sensing, carrying, or storing. To gain insight into its possible functional roles, we have investigated its O2 and NO rebinding kinetics after laser flash photolysis. The rate constants of the rebinding from the primary docking site for both O2 and NO are higher than CO, with lower photolysis yields. Moreover, the amplitude of the geminate phase increases and, as for CO, the numerical analysis of the experimental curves suggests the existence of an internal pathway leading, with high rate, to an additional docking site. However, the accessibility to this site seems to be strongly ligand-dependent, being lower for O2 and higher for NO.     

Kinetics and computational studies of ligand migration in nitrophorin 7 and its Δ1–3 mutant

Nitrophorins (NPs) are nitric oxide (NO)-carrying heme proteins found in the saliva of the blood-sucking insect Rhodnius prolixus. Though NP7 exhibits a large sequence resemblance with other NPs, two major differential features are the ability to interact with negatively charged cell surfaces and the presence of a specific N-terminus composed of three extra residues (Leu1-Pro2-Gly3). The aim of this study is to examine the influence of the N-terminus on the ligand binding, and the topological features of inner cavities in closed and open states of NP7, which can be associated to the protein structure at low and high pH, respectively. Laser flash photolysis measurements of the CO rebinding kinetics to NP7 and its variant NP7(Δ1–3), which lacks the three extra residues at the N-terminus, exhibit a similar pattern and support the existence of a common kinetic mechanism for ligand migration and binding. This is supported by the existence of a common topology of inner cavities, which consists of two docking sites in the heme pocket and a secondary site at the back of the protein. The ligand exchange between these cavities is facilitated by an additional site, which can be transiently occupied by the ligand in NP7, although it is absent in NP4. These features provide a basis to explain the enhanced internal gas hosting capacity found experimentally in NP7 and the absence of ligand rebinding from secondary sites in NP4. The current data allow us to speculate that the processes of docking to cell surfaces and NO release may be interconnected in NP7, thereby efficiently releasing NO into a target cell. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Oxygen binding to Arabidopsis thaliana AHb2 nonsymbiotic hemoglobin: evidence for a role in oxygen transport

Nonsymbiotic hemoglobins AHb1 and AHb2 discovered in Arabidopsis thaliana are likely to carry out distinct physiological roles, in consideration of their differences in sequence, structure, expression pattern, and tissue localization. Despite a relatively fast autoxidation in the presence of O(2) , we were able to collect O(2) -binding curves for AHb2 in the presence of a reduction enzymatic system. AHb2 binds O(2) noncooperatively with a p50 of 0.021 ± 0.003 Torr, a value consistent with a recently proposed role in O(2) transport. The analysis of the internal cavities derived from the structures sampled in molecular dynamics simulations confirms strong differences with AHb1, proposed to work as a NO deoxygenase in vivo. Overall, our results are consistent with a role for AHb2 as an oxygen carrier, as recently proposed on the basis of experiments on AHb2-overexpressing mutants of A. thaliana.     

Electron transfer function versus oxygen delivery: a comparative study for several hexacoordinated globins across the animal kingdom

Caenorhabditis elegans globin GLB-26 (expressed from gene T22C1.2) has been studied in comparison with human neuroglobin (Ngb) and cytoglobin (Cygb) for its electron transfer properties. GLB-26 exhibits no reversible binding for O(2) and a relatively low CO affinity compared to myoglobin-like globins. These differences arise from its mechanism of gaseous ligand binding since the heme iron of GLB-26 is strongly hexacoordinated in the absence of external ligands; the replacement of this internal ligand, probably the E7 distal histidine, is required before binding of CO or O(2) as for Ngb and Cygb. Interestingly the ferrous bis-histidyl GLB-26 and Ngb, another strongly hexacoordinated globin, can transfer an electron to cytochrome c (Cyt-c) at a high bimolecular rate, comparable to those of inter-protein electron transfer in mitochondria. In addition, GLB-26 displays an unexpectedly rapid oxidation of the ferrous His-Fe-His complex without O(2) actually binding to the iron atom, since the heme is oxidized by O(2) faster than the time for distal histidine dissociation. These efficient mechanisms for electron transfer could indicate a family of hexacoordinated globin which are functionally different from that of pentacoordinated globins.     

Crystal structure of a nonsymbiotic plant hemoglobin

BACKGROUND: Nonsymbiotic hemoglobins (nsHbs) form a new class of plant proteins that is distinct genetically and structurally from leghemoglobins. They are found ubiquitously in plants and are expressed in low concentrations in a variety of tissues including roots and leaves. Their function involves a biochemical response to growth under limited O(2) conditions. RESULTS: The first X-ray crystal structure of a member of this class of proteins, riceHb1, has been determined to 2.4 A resolution using a combination of phasing techniques. The active site of ferric riceHb1 differs significantly from those of traditional hemoglobins and myoglobins. The proximal and distal histidine sidechains coordinate directly to the heme iron, forming a hemichrome with spectral properties similar to those of cytochrome b(5). The crystal structure also shows that riceHb1 is a dimer with a novel interface formed by close contacts between the G helix and the region between the B and C helices of the partner subunit. CONCLUSIONS: The bis-histidyl heme coordination found in riceHb1 is unusual for a protein that binds O(2) reversibly. However, the distal His73 is rapidly displaced by ferrous ligands, and the overall O(2) affinity is ultra-high (K(D) approximately 1 nM). Our crystallographic model suggests that ligand binding occurs by an upward and outward movement of the E helix, concomitant dissociation of the distal histidine, possible repacking of the CD corner and folding of the D helix. Although the functional relevance of quaternary structure in nsHbs is unclear, the role of two conserved residues in stabilizing the dimer interface has been identified.     

Unraveling the molecular basis for ligand binding in truncated hemoglobins: The trHbO Bacillus subtilis case

Truncated hemoglobins (trHbs) are heme proteins present in bacteria, unicellular eukaryotes, and higher plants. Their tertiary structure consists in a 2‐over‐2 helical sandwich, which display typically an inner tunnel/cavity system for ligand migration and/or storage. The microorganism Bacillus subtilis contains a peculiar trHb, which does not show an evident tunnel/cavity system connecting the protein active site with the solvent, and exhibits anyway a very high oxygen association rate. Moreover, resonant Raman results of CO bound protein, showed that a complex hydrogen bond network exists in the distal cavity, making it difficult to assign unambiguously the residues involved in the stabilization of the bound ligand. To understand these experimental results with atomistic detail, we performed classical molecular dynamics simulations of the oxy, carboxy, and deoxy proteins. The free energy profiles for ligand migration suggest that there is a key residue, GlnE11, that presents an alternate conformation, in which a wide ligand migration tunnel is formed, consistently with the kinetic data. This tunnel is topologically related to the one found in group I trHbs. On the other hand, the results for the CO and O₂ bound protein show that GlnE11 is directly involved in the stabilization of the cordinated ligand, playing a similar role as TyrB10 and TrpG8 in other trHbs. Our results not only reconcile the structural data with the kinetic information, but also provide additional insight into the general behaviour of trHbs. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

HMGB binding to DNA: single and double box motifs

High mobility group (HMG) proteins are nuclear proteins believed to significantly affect DNA interactions by altering nucleic acid flexibility. Group B (HMGB) proteins contain HMG box domains known to bind to the DNA minor groove without sequence specificity, slightly intercalating base pairs and inducing a strong bend in the DNA helical axis. A dual-beam optical tweezers system is used to extend double-stranded DNA (dsDNA) in the absence as well as presence of a single box derivative of human HMGB2 [HMGB2(box A)] and a double box derivative of rat HMGB1 [HMGB1(box A+box B)]. The single box domain is observed to reduce the persistence length of the double helix, generating sharp DNA bends with an average bending angle of 99 ± 9° and, at very high concentrations, stabilizing dsDNA against denaturation. The double box protein contains two consecutive HMG box domains joined by a flexible tether. This protein also reduces the DNA persistence length, induces an average bending angle of 77 ± 7°, and stabilizes dsDNA at significantly lower concentrations. These results suggest that single and double box proteins increase DNA flexibility and stability, albeit both effects are achieved at much lower protein concentrations for the double box. In addition, at low concentrations, the single box protein can alter DNA flexibility without stabilizing dsDNA, whereas stabilization at higher concentrations is likely achieved through a cooperative binding mode.     

CO Rebinding Kinetics and Molecular Dynamics Simulations Highlight Dynamic Regulation of Internal Cavities in Human Cytoglobin

Cytoglobin (Cygb) was recently discovered in the human genome and localized in different tissues. It was suggested to play tissue-specific protective roles, spanning from scavenging of reactive oxygen species in neurons to supplying oxygen to enzymes in fibroblasts. To shed light on the functioning of such versatile machinery, we have studied the processes supporting transport of gaseous heme ligands in Cygb. Carbon monoxide rebinding shows a complex kinetic pattern with several distinct reaction intermediates, reflecting rebinding from temporary docking sites, second order recombination, and formation (and dissociation) of a bis-histidyl heme hexacoordinated reaction intermediate. Ligand exit to the solvent occurs through distinct pathways, some of which exploit temporary docking sites. The remarkable change in energetic barriers, linked to heme bis-histidyl hexacoordination by HisE7, may be responsible for active regulation of the flux of reactants and products to and from the reaction site on the distal side of the heme. A substantial change in both protein dynamics and inner cavities is observed upon transition from the CO-liganded to the pentacoordinated and bis-histidyl hexacoordinated species, which could be exploited as a signalling state. These findings are consistent with the expected versatility of the molecular activity of this protein.     

Solution structure and backbone dynamics of human Raf-1 kinase inhibitor protein

Human Raf-1 kinase inhibitor protein (hRKIP) is a small multi-functional protein of 187 residues. It contains a conserved pocket, which binds a wide range of ligands from various small molecules to distinct proteins. To provide a structural basis for the ligand diversity of RKIP, we herein determined the solution structure of hRKIP, and analyzed its structural dynamics. In solution, hRKIP mainly comprises two antiparallel β sheets, two α helices and two 3₁₀ helices. NMR dynamic analysis reveals that the overall structure of hRKIP is rigid, but its C-terminal helix which is close to the ligand-binding site is mobile. In addition, residues around the ligand-binding pocket exhibit significant conformational exchange on the μs–ms timescale. Conformational flexibility may allow the ligand-binding pocket and the C-terminal helix to adopt various conformations to interact with different substrates. This work may shed light on the underlying molecular mechanisms of how hRKIP recognizes and binds diverse substrate ligands.     

Arabidopsis nonsymbiotic hemoglobin AHb1 modulates nitric oxide bioactivity

Nitric oxide (NO) is a widespread signaling molecule, and numerous targets of its action exist in plants. Whereas the activity of NO in erythrocytes, microorganisms, and invertebrates has been shown to be regulated by several hemoglobins, the function of plant hemoglobins in NO detoxification has not yet been elucidated. Here, we show that Arabidopsis thaliana nonsymbiotic hemoglobin AHb1 scavenges NO through production of S-nitrosohemoglobin and reduces NO emission under hypoxic stress, indicating its role in NO detoxification. However, AHb1 does not affect NO-mediated hypersensitive cell death in response to avirulent Pseudomonas syringae, suggesting that it is not involved in the removal of NO bursts originated from acute responses when NO mediates crucial defense signaling functions.     

Relative Domain Folding and Stability of a Membrane Transport Protein

There is a limited understanding of the folding of multidomain membrane proteins. Lactose permease (LacY) of Escherichia coli is an archetypal member of the major facilitator superfamily of membrane transport proteins, which contain two domains of six transmembrane helices each. We exploit chemical denaturation to determine the unfolding free energy of LacY and employ Trp residues as site-specific thermodynamic probes. Single Trp LacY mutants are created with the individual Trps situated at mirror image positions on the two LacY domains. The changes in Trp fluorescence induced by urea denaturation are used to construct denaturation curves from which unfolding free energies can be determined. The majority of the single Trp tracers report the same stability and an unfolding free energy of approximately + 2 kcal mol^− 1. There is one exception; the fluorescence of W33 at the cytoplasmic end of helix I on the N domain is unaffected by urea. In contrast, the equivalent position on the first helix, VII, of the C-terminal domain exhibits wild-type stability, with the single Trp tracer at position 243 on helix VII reporting an unfolding free energy of + 2 kcal mol^− 1. This indicates that the region of the N domain of LacY at position 33 on helix I has enhanced stability to urea, when compared the corresponding location at the start of the C domain. We also find evidence for a potential network of stabilising interactions across the domain interface, which reduces accessibility to the hydrophilic substrate binding pocket between the two domains.     

Probing hydrogen peroxide oxidation kinetics of wild-type Synechocystis catalase-peroxidase (KatG) and selected variants

Catalase-peroxidases (KatGs) are unique bifunctional heme peroxidases that exhibit peroxidase and substantial catalase activities. Nevertheless, the reaction pathway of hydrogen peroxide dismutation, including the electronic structure of the redox intermediate that actually oxidizes H₂O₂, is not clearly defined. Several mutant proteins with diminished overall catalase but wild-type-like peroxidase activity have been described in the last years. However, understanding of decrease in overall catalatic activity needs discrimination between reduction and oxidation reactions of hydrogen peroxide. Here, by using sequential-mixing stopped-flow spectroscopy, we have investigated the kinetics of the transition of KatG compound I (produced by peroxoacetic acid) to its ferric state by trapping the latter as cyanide complex. Apparent bimolecular rate constants (pH 6.5, 20 °C) for wild-type KatG and the variants Trp122Phe (lacks KatG-typical distal adduct), Asp152Ser (controls substrate access to the heme cavity) and Glu253Gln (channel entrance) are reported to be 1.2 × 10⁴ M^− 1 s^− 1, 30 M^− 1 s^− 1, 3.4 × 10³ M^− 1 s^− 1, and 8.6 × 10³ M^− 1 s^− 1, respectively. These findings are discussed with respect to steady-state kinetic data and proposed reaction mechanism(s) for KatG. Assets and drawbacks of the presented method are discussed.     

The comparative study on the solution structures of the oxidized bovine microsomal cytochrome b5 and mutant V45H

A comparative study on the solution structures of bovine microsomal cytochrome b5 (Tb5) and the mutant V45H has been achieved by 1D and 2D 1H-NMR spectroscopy to clarify the differences in the solution conformations between these two proteins. The results reveal that the global folding of the V45H mutant in solution is unchanged, but the subtle changes exist in the orientation of the axial ligand His39, and heme vinyl groups. The side chain of His45 in V45H mutant extends to the outer edge of the heme pocket leaving a cavity at the site originally occupied by the inner methyl group of Val45 residue. In addition, the imidazole ring of axial ligand His39 rotates counterclockwise by approximately 3 degrees around the His-Fe-His axis, and the 4-heme vinyl group turns to the space vacated by the removed side chain due to the mutation. Furthermore, the helix III of the heme pocket undergoes outward displacement, while the linkage between helix II and III is shifted leftward. These observations are not only consistent with the pattern of the pseudocontact shifts of the heme protons, but also well account for the lower stability of V45H mutant against heat and urea.     

Low-mass molecular dynamics simulation: A simple and generic technique to enhance configurational sampling

CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal–isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277 K and 1 atm with the first folding event occurring as early as 66.1 ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal–isobaric MD simulations performed on commodity computers—an important step forward in quantitative biology.