Heme pocket structural properties of a bacterial truncated hemoglobin from Thermobifida fusca

Proximal ligand tunes active site structure and reactivity in bacterial L. monocytogenes coproheme ferrochelatase

Ferrochelatases catalyze the insertion of ferrous iron into the porphyrin during the heme b biosynthesis pathway, which is fundamental for both prokaryotes and eukaryotes. Interestingly, in the active site of ferrochelatases, the proximal ligand coordinating the porphyrin iron of the product is not conserved, and its catalytic role is still unclear. Here we compare the L. monocytogenes bacterial coproporphyrin ferrochelatase native enzyme together with selected variants, where the proximal Tyr residue was replaced by a His (i.e. the most common ligand in heme proteins), a Met or a Phe (as in human and actinobacterial ferrochelatases, respectively), in their Fe(III), Fe(II) and Fe(II)-CO adduct forms. The study of the active site structure and the activity of the proteins in solution has been performed by UV-vis electronic absorption and resonance Raman spectroscopies, biochemical characterization, and classical MD simulations. All the mutations alter the H-bond interactions between the iron porphyrin propionate groups and the protein, and induce effects on the activity, depending on the polarity of the proximal ligand. The overall results confirm that the weak or non-existing coordination of the porphyrin iron by the proximal residue is essential for the binding of the substrate and the release of the final product.

Truncated (2/2) hemoglobin: Unconventional structures and functional roles in vivo and in human pathogenesis

Truncated hemoglobins (trHbs) build a sub-class of the globin family, found in eubacteria, cyanobacteria, unicellular eukaryotes, and in higher plants; among these, selected human pathogens are found. The trHb fold is based on a 2/2 α-helical sandwich, consisting of a simplified and reduced-size version of the classical 3/3 α-helical sandwich of vertebrate and invertebrate globins. Phylogenetic analysis indicates that trHbs further branch into three groups: group I (or trHbN), group II (or trHbO), and group III (or trHbP), each group being characterized by specific structural features. Among these, a protein matrix tunnel, or a cavity system implicated in diatomic ligand diffusion through the protein matrix, is typical of group I and group II, respectively. In general, a highly intertwined network of hydrogen bonds stabilizes the heme bound ligand, despite variability of the heme distal residues in the different trHb groups. Notably, some organisms display genes from more than one trHb group, suggesting that trHbN, trHbO, and trHbP may support different functions in vivo, such as detoxification of reactive nitrogen and oxygen species, respiration, oxygen storage/sensoring, thus aiding survival of an invading microorganism. Here, structural features and proposed functions of trHbs from human pathogens are reviewed.

Conformational Flexibility Drives Cold Adaptation in Pseudoalteromonas haloplanktis TAC125 Globins

Significance: Temperature is one of the most important drivers in shaping protein adaptations. Many biochemical and physiological processes are influenced by temperature. Proteins and enzymes from organisms living at low temperature are less stable in comparison to high-temperature adapted proteins. The lower stability is generally due to greater conformational flexibility. Recent Advances: Adaptive changes in the structure of cold-adapted proteins may occur at subunit interfaces, distant from the active site, thus producing energy changes associated with conformational transitions transmitted to the active site by allosteric modulation, valid also for monomeric proteins in which tertiary structural changes may play an essential role. Critical Issues: Despite efforts, the current experimental and computational methods still fail to produce general principles on protein evolution, since many changes are protein and species dependent. Environmental constraints or other biological cellular signals may override the ancestral information included in the structure of the protein, thus introducing inaccuracy in estimates and predictions on the evolutionary adaptations of proteins in response to cold adaptation. Future Directions: In this review, we describe the studies and approaches used to investigate stability and flexibility in the cold-adapted globins of the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125. In fact, future research directions will be prescient on more detailed investigation of cold-adapted proteins and the role of fluctuations between different conformational states.

Active Sites of O2-Evolving Chlorite Dismutases Probed by Halides and Hydroxides and New Iron-Ligand Vibrational Correlations

O₂-evolving chlorite dismutases (Clds) fall into two subfamilies, which efficiently convert ClO₂^- to O₂ and Cl^-. The Cld from Dechloromonas aromatica (DaCld) represents the chlorite-decomposing homopentameric enzymes found in perchlorate- and chlorate-respiring bacteria. The Cld from the Gram-negative human pathogen Klebsiella pneumoniae (KpCld) is representative of the second subfamily, comprising homodimeric enzymes having truncated N-termini. Here steric and nonbonding properties of the DaCld and KpCld active sites have been probed via kinetic, thermodynamic, and spectroscopic behaviors of their fluorides, chlorides, and hydroxides. Cooperative binding of Cl^- to KpCld drives formation of a hexacoordinate, high-spin aqua heme, whereas DaCld remains pentacoordinate and high-spin under analogous conditions. Fluoride coordinates to the heme iron in KpCld and DaCld, exhibiting ν(Fe^III-F) bands at 385 and 390 cm^-1, respectively. Correlation of these frequencies with their CT1 energies reveals strong H-bond donation to the F^- ligand, indicating that atoms directly coordinated to heme iron are accessible to distal H-bond donation. New vibrational frequency correlations between either ν(Fe^III-F) or ν(Fe^III-OH) and ν(Fe^II-His) of Clds and other heme proteins are reported. These correlations orthogonalize proximal and distal effects on the bonding between iron and exogenous π-donor ligands. The axial Fe-X vibrations and the relationships between them illuminate both similarities and differences in the H-bonding and electrostatic properties of the distal and proximal heme environments in pentameric and dimeric Clds. Moreover, they provide general insight into the structural basis of reactivity toward substrates in heme-dependent enzymes and their mechanistic intermediates, especially those containing the ferryl moiety.

Structural and Functional Significance of the N- and C-Terminal Appendages in Arabidopsis Truncated Hemoglobin

Plant hemoglobins constitute three distinct groups: symbiotic, nonsymbiotic, and truncated hemoglobins. Structural investigation of symbiotic and nonsymbiotic (class I) hemoglobins revealed the presence of a vertebrate-like 3/3 globin fold in these proteins. In contrast, plant truncated hemoglobins are similar to bacterial truncated hemoglobins with a putative 2/2 α-helical globin fold. While multiple structures have been reported for plant hemoglobins of the first two categories, for plant truncated globins only one structure has been reported of late. Here, we report yet another crystal structure of the truncated hemoglobin from Arabidopsis thaliana (AHb3) with two water molecules in the heme pocket, of which one is distinctly coordinated to the heme iron, unlike the only available crystal structure of AHb3 with a hydroxyl ligand. AHb3 was monomeric in its crystallographic asymmetric unit; however, dimer was evident in the crystallographic symmetry, and the globin indeed existed as a stable dimer in solution. The tertiary structure of the protein exhibited a bacterial-like 2/2 α-helical globin fold with an additional N-terminal α-helical extension and disordered C-termini. To address the role of these extended termini in AHb3, which is yet unknown, N- and C-terminal deletion mutants were created and characterized and molecular dynamics simulations performed. The C-terminal deletion had an insignificant effect on most properties but perturbed the dimeric equilibrium of AHb3 and significantly influenced azide binding kinetics in the ferric state. These results along with the disordered nature of the C-terminus indicated its putative role in intramolecular or intermolecular interactions probably regulating protein-ligand and protein-protein interactions. While the N-terminal deletion did not change the overall globin fold, stability, or ligand binding kinetics, it seemed to have influenced coordination at the heme iron, the hydration status of the active site, and the quaternary structure of AHb3. Evidence indicated that the N-terminus is the predominant factor regulating the quaternary interaction appropriate to physiological requirements, dynamics of the side chains in the heme pocket, and tunnel organization in the protein matrix.

Structural flexibility of the heme cavity in the cold-adapted truncated hemoglobin from the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125

Truncated hemoglobins build one of the three branches of the globin protein superfamily. They display a characteristic two-on-two α-helical sandwich fold and are clustered into three groups (I, II and III) based on distinct structural features. Truncated hemoglobins are present in eubacteria, cyanobacteria, protozoa and plants. Here we present a structural, spectroscopic and molecular dynamics characterization of a group-II truncated hemoglobin, encoded by the PSHAa0030 gene from Pseudoalteromonas haloplanktis TAC125 (Ph-2/2HbO), a cold-adapted Antarctic marine bacterium hosting one flavohemoglobin and three distinct truncated hemoglobins. The Ph-2/2HbO aquo-met crystal structure (at 2.21 Å resolution) shows typical features of group-II truncated hemoglobins, namely the two-on-two α-helical sandwich fold, a helix Φ preceding the proximal helix F, and a heme distal-site hydrogen-bonded network that includes water molecules and several distal-site residues, including His(58)CD1. Analysis of Ph-2/2HbO by electron paramagnetic resonance, resonance Raman and electronic absorption spectra, under varied solution conditions, shows that Ph-2/2HbO can access diverse heme ligation states. Among these, detection of a low-spin heme hexa-coordinated species suggests that residue Tyr(42)B10 can undergo large conformational changes in order to act as the sixth heme-Fe ligand. Altogether, the results show that Ph-2/2HbO maintains the general structural features of group-II truncated hemoglobins but displays enhanced conformational flexibility in the proximity of the heme cavity, a property probably related to the functional challenges, such as low temperature, high O2 concentration and low kinetic energy of molecules, experienced by organisms living in the Antarctic environment.

Functional and Spectroscopic Characterization of Chlamydomonas reinhardtii Truncated Hemoglobins

The single-cell green alga Chlamydomonas reinhardtii harbors twelve truncated hemoglobins (Cr-TrHbs). Cr-TrHb1-1 and Cr-TrHb1-8 have been postulated to be parts of the nitrogen assimilation pathway, and of a NO-dependent signaling pathway, respectively. Here, spectroscopic and reactivity properties of Cr-TrHb1-1, Cr-TrHb1-2, and Cr-TrHb1-4, all belonging to clsss 1 (previously known as group N or group I), are reported. The ferric form of Cr-TrHb1-1, Cr-TrHb1-2, and Cr-TrHb1-4 displays a stable 6cLS heme-Fe atom, whereas the hexa-coordination of the ferrous derivative appears less strongly stabilized. Accordingly, kinetics of azide binding to ferric Cr-TrHb1-1, Cr-TrHb1-2, and Cr-TrHb1-4 are independent of the ligand concentration. Conversely, kinetics of CO or NO₂⁻ binding to ferrous Cr-TrHb1-1, Cr-TrHb1-2, and Cr-TrHb1-4 are ligand-dependent at low CO or NO₂⁻ concentrations, tending to level off at high ligand concentrations, suggesting the presence of a rate-limiting step. In agreement with the different heme-Fe environments, the pH-dependent kinetics for CO and NO₂−binding to ferrous Cr-TrHb1-1, Cr-TrHb1-2, and Cr-TrHb1-4 are characterized by different ligand-linked protonation events. This raises the question of whether the simultaneous presence in C. reinhardtii of multiple TrHb1s may be related to different regulatory roles.

Bridging Theory and Experiment to Address Structural Properties of Truncated Haemoglobins: Insights from Thermobifida fusca HbO

In this chapter, we will discuss the paradigmatic case of Thermobifida fusca (Tf-trHb) HbO in its ferrous and ferric states and its behaviour towards a battery of possible ligands. This choice was dictated by the fact that it has been one of the most extensively studied truncated haemoglobins, both in terms of spectroscopic and molecular dynamics studies. Tf-trHb typifies the structural properties of group II trHbs, as the active site is characterized by a highly polar distal environment in which TrpG8, TyrCD1, and TyrB10 provide three potential H-bond donors in the distal cavity capable of stabilizing the incoming ligands. The role of these residues in key topological positions, and their interplay with the iron-bound ligands, has been addressed in studies carried out on the CO, F(-), OH(-), CN(-), and HS(-) adducts formed with the wild-type protein and a combinatorial set of mutants, in which the distal polar residues, TrpG8, TyrCD1, and TyrB10, have been singly, doubly, or triply replaced by a Phe residue. In this context, such a complete analysis provides an excellent benchmark for the investigation of the relationship between protein structure and function, allowing one to translate physicochemical properties of the active site into the observed functional behaviour. Tf-trHb will be compared with other members of the group II trHbs and, more generally, with members of the other trHb subgroups.

Role of local structure and dynamics of small ligand migration in proteins: a study of a mutated truncated hemoprotein from Thermobifida fusca by time resolved MIR spectroscopy

Carbon monoxide recombination dynamics in a mutant of the truncated hemoglobin from Thermobida fusca (3F-Tf-trHb) has been analyzed by means of ultrafast Visible-pump/MidIR-probe spectroscopy and compared with that of the wild-type protein. In 3F-Tf-trHb, three topologically relevant amino acids, responsible for the ligand stabilization through the formation of a H-bond network (TyrB10 TyrCD1 and TrpG8), have been replaced by Phe residues. X-ray diffraction data show that Phe residues in positions B10 and G8 maintain the same rotameric arrangements as Tyr and Trp in the wild-type protein, while Phe in position CD1 displays significant rotameric heterogeneity. Photodissociation of the ligand has been induced by exciting the sample with 550 nm pump pulses and the CO rebinding has been monitored in two mid-IR regions respectively corresponding to the ν(CO) stretching vibration of the iron-bound CO (1880-1980 cm(-1)) and of the dissociated free CO (2050-2200 cm(-1)). In both the mutant and wild-type protein, a significant amount of geminate CO rebinding is observed on a subnanosecond time scale. Despite the absence of the distal pocket hydrogen-bonding network, the kinetics of geminate rebinding in 3F-Tf-trHb is very similar to the wild-type, showing how the reactivity of dissociated CO toward the heme is primarily regulated by the effective volume and flexibility of the distal pocket and by caging effects exerted on the free CO on the analyzed time scale.

Molecular basis of thermal stability in truncated (2/2) hemoglobins

BACKGROUND

Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures.

METHODS

We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position.

RESULTS

The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible.

CONCLUSIONS

This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability.

GENERAL SIGNIFICANCE

These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

The diversity of 2/2 (truncated) globins

Small size globins that have been defined as 'truncated haemoglobins' or as '2/2 haemoglobins' have increasingly been discovered in microorganisms since the early 1990s. Analysis of amino acid sequences allowed to distinguish three groups that collect proteins with specific and common structural properties. All three groups display 3D structures that are based on four main α-helices, which are a subset of the conventional eight-helices globin fold. Specific features, such as the presence of protein matrix tunnels that are held to promote diffusion of functional ligands to/from the haem, distinguish members of the three groups. Haem distal sites vary for their accessibility, local structures, polarity, and ligand stabilization mechanisms, suggesting functional roles that are related to O2/NO chemistry. In a few cases, such activities have been proven in vitro and in vivo through deletion mutants. The issue of 2/2 haemoglobin varied biological functions throughout the three groups remains however fully open.

Haemoglobins of Mycobacteria: structural features and biological functions

The genus Mycobacterium is comprised of Gram-positive bacteria occupying a wide range of natural habitats and includes species that range from severe intracellular pathogens to economically useful and harmless microbes. The recent upsurge in the availability of microbial genome data has shown that genes encoding haemoglobin-like proteins are ubiquitous among Mycobacteria and that multiple haemoglobins (Hbs) of different classes may be present in pathogenic and non-pathogenic species. The occurrence of truncated haemoglobins (trHbs) and flavohaemoglobins (flavoHbs) showing distinct haem active site structures and ligand-binding properties suggests that these Hbs may be playing diverse functions in the cellular metabolism of Mycobacteria. TrHbs and flavoHbs from some of the severe human pathogens such as Mycobacterium tuberculosis and Mycobacterium leprae have been studied recently and their roles in effective detoxification of reactive nitrogen and oxygen species, electron cycling, modulation of redox state of the cell and facilitation of aerobic respiration have been proposed. This multiplicity in the function of Hbs may aid these pathogens to cope with various environmental stresses and survive during their intracellular regime. This chapter provides recent updates on genomic, structural and functional aspects of Mycobacterial Hbs to address their role in Mycobacteria.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

H-bonding networks of the distal residues and water molecules in the active site of Thermobifida fusca hemoglobin

The ferric form of truncated hemoglobin II from Thermobifida fusca (Tf-trHb) and its triple mutant WG8F-YB10F-YCD1F at neutral and alkaline pH, and in the presence of CN(-) have been characterized by resonance Raman spectroscopy, electron paramagnetic resonance spectroscopy, and molecular dynamics simulations. Tf-trHb contains three polar residues in the distal site, namely TrpG8, TyrCD1 and TyrB10. Whereas TrpG8 can act as a potential hydrogen-bond donor, the tyrosines can act as donors or acceptors. Ligand binding in heme-containing proteins is determined by a number of factors, including the nature and conformation of the distal residues and their capability to stabilize the heme-bound ligand via hydrogen-bonding and electrostatic interactions. Since both the RR Fe-OH(-) and Fe-CN(-) frequencies are very sensitive to the distal environment, detailed information on structural variations has been obtained. The hydroxyl ligand binds only the WT protein giving rise to two different conformers. In form 1 the anion is stabilized by H-bonds with TrpG8, TyrCD1 and a water molecule, in turn H-bonded to TyrB10. In form 2, H-bonding with TyrCD1 is mediated by a water molecule. Unlike the OH(-) ligand, CN(-) binds both WT and the triple mutant giving rise to two forms with similar spectroscopic characteristics. The overall results clearly indicate that H-bonding interactions both with distal residues and water molecules are important structural determinants in the active site of Tf-trHb. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Small ligand-globin interactions: reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Carbon monoxide recombination dynamics in truncated hemoglobins studied with visible-pump midIR-probe spectroscopy

Carbon monoxide recombination dynamics upon photodissociation with visible light has been characterized by means of ultrafast visible-pump/MidIR probe spectroscopy for the truncated hemoglobins from Thermobifida fusca and Bacillus subtilis. Photodissociation has been induced by exciting the sample at two different wavelengths: 400 nm, corresponding to the heme absorption in the B-band, and 550 nm, in the Q-bands. The bleached iron-CO coordination band located at 1850-1950 cm(-1) and the free CO absorption band in the region 2050-2200 cm(-1) have been observed by probe pulses tuned in the appropriate infrared region. The kinetic traces measured at 1850-1950 cm(-1) reveal multiexponential subnanosecond dynamics that have been interpreted as arising from fast geminate recombination of the photolyzed CO. A compared analysis of the crystal structure of the two proteins reveals a similar structure of their distal heme pocket, which contains conserved polar and aromatic amino acid residues closely interacting with the iron ligand. Although fast geminate recombination is observed in both proteins, several kinetic differences can be evidenced, which can be interpreted in terms of a different structural flexibility of the corresponding heme distal pockets. The analysis of the free CO band-shape and of its dynamic evolution brings out novel features about the nature of the docking site inside the protein cavity.

Following ligand migration pathways from picoseconds to milliseconds in type II truncated hemoglobin from Thermobifida fusca

CO recombination kinetics has been investigated in the type II truncated hemoglobin from Thermobifida fusca (Tf-trHb) over more than 10 time decades (from 1 ps to ∼100 ms) by combining femtosecond transient absorption, nanosecond laser flash photolysis and optoacoustic spectroscopy. Photolysis is followed by a rapid geminate recombination with a time constant of ∼2 ns representing almost 60% of the overall reaction. An additional, small amplitude geminate recombination was identified at ∼100 ns. Finally, CO pressure dependent measurements brought out the presence of two transient species in the second order rebinding phase, with time constants ranging from ∼3 to ∼100 ms. The available experimental evidence suggests that the two transients are due to the presence of two conformations which do not interconvert within the time frame of the experiment. Computational studies revealed that the plasticity of protein structure is able to define a branched pathway connecting the ligand binding site and the solvent. This allowed to build a kinetic model capable of describing the complete time course of the CO rebinding kinetics to Tf-trHb.

Ligand- and proton-linked conformational changes of the ferrous 2/2 hemoglobin of Pseudoalteromonas haloplanktis TAC125

The spectroscopic and ligand-binding properties of a 2/2 globin from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 have been studied in the ferrous state. It displays two major conformations characterized by CO-association rates that differ by a factor of 20, with relative fractions that depend on pH. A dynamic equilibrium is found between the two conformations, as indicated by an enhanced slower phase when lower CO levels were used to allow a longer time to facilitate the transition. The deoxy form, in the absence of external ligands, is a mixture of a predominant six-coordinate low spin form and a five-coordinate high-spin state; the proportion of low spin increasing at alkaline pH. In addition, at temperatures above the physiological temperature of 1 °C, an enhanced tendency of the protein to oxidize is observed.

Structural analysis in nonsymbiotic hemoglobins: what can we learn from inner cavities?

Plants contain three classes of hemoglobins which are not associated with nitrogen fixing bacteria, and have been accordingly termed nonsymbiotic hemoglobins. The function of nonsymbiotic hemoglobins is as yet mostly unknown. A NO dioxygenase activity has been proposed and demonstrated for some of them in vitro. In this context, a sound molecular mechanism that relates the structure with the biological activity is crucial to suggest a given physiological role. Insight into such a mechanism is now facilitated by recent progress made in both experimental and computational techniques. These studies have highlighted a number of key structural features implicated in the function of nonsymbiotic hemoglobins. The bis-histidyl hexacoordination of the heme in both its ferric and ferrous states provides a powerful and general tool to modulate reactivity, protein dynamics, and shape of the cavities. In addition, the specific arrangement of distal cavity residues provides effective protection against autoxidation. Inspection of the static crystal structures available for both liganded and unliganded states seems unsufficient to explain the function of these proteins. Function appears to be intimately linked with protein flexibility, which influences the dynamical behavior of inner cavities, capable of delivering apolar reactants to the reaction site, and removing charged reaction products. In this mini review, we demonstrate how the integration of information derived from experimental assays and computational studies is valuable and can shed light into the linkage between structural plasticity of nonsymbiotic hemoglobins and their biological role.