Implementing a web-based introductory bioinformatics course for non-bioinformaticians that incorporates practical exercises.

A recent scientific discipline, bioinformatics, defined as using informatics for the study of biological problems, is now a requirement for the study of biological sciences. Bioinformatics has become such a powerful and popular discipline that several academic institutions have created programs in this field, allowing students to become specialized. However, biology students who are not involved in a bioinformatics program also need a solid toolbox of bioinformatics software and skills. Therefore, we have developed a completely online bioinformatics course for non-bioinformaticians, entitled "BIF-1901 Introduction à la bio-informatique et à ses outils (Introduction to bioinformatics and bioinformatics tools)," given by the Department of Biochemistry, Microbiology, and Bioinformatics of Université Laval (Quebec City, Canada). This course requires neither a bioinformatics background nor specific skills in informatics. The underlying main goal was to produce a completely online up-to-date bioinformatics course, including practical exercises, with an intuitive pedagogical framework. The course, BIF-1901, was conceived to cover the three fundamental aspects of bioinformatics: (1) informatics, (2) biological sequence analysis, and (3) structural bioinformatics. This article discusses the content of the modules, the evaluations, the pedagogical framework, and the challenges inherent to a multidisciplinary, fully online course. © 2017 by The International Union of Biochemistry and Molecular Biology, 46(1):31-38, 2018.

Mechanism of the Nitric Oxide Dioxygenase Reaction of Mycobacterium tuberculosis Hemoglobin N.

Many globins convert •NO to innocuous NO3- through their nitric oxide dioxygenase (NOD) activity. Mycobacterium tuberculosis fights the oxidative and nitrosative stress imposed by its host (the toxic effects of O2•- and •NO species and their OONO- and •NO2 derivatives) through the action of truncated hemoglobin N (trHbN), which catalyzes the NOD reaction with one of the highest rates among globins. The general NOD mechanism comprises the following steps: binding of O2 to the heme, diffusion of •NO into the heme pocket and formation of peroxynitrite (OONO-), isomerization of OONO-, and release of NO3-. Using quantum mechanics/molecular mechanics free-energy calculations, we show that the NOD reaction in trHbN follows a mechanism in which heme-bound OONO- undergoes homolytic cleavage to give FeIV═O2- and the •NO2 radical but that these potentially harmful intermediates are short-lived and caged by the heme pocket residues. In particular, the simulations show that Tyr33(B10) side chain is shielded from FeIV═O2- and •NO2 (and protected from irreversible oxidation and nitration) by forming stable hydrogen bonds with Gln58(E11) side chain and Leu54(E7) backbone. Aromatic residues Phe46(CD1), Phe32(B9), and Tyr33(B10) promote NO3- dissociation via C-H···O bonding and provide stabilizing interactions for the anion along its egress route.

The N-terminal pre-A region of Mycobacterium tuberculosis 2/2HbN promotes NO-dioxygenase activity.

UNLABELLED: A unique defense mechanisms by which Mycobacterium tuberculosis protects itself from nitrosative stress is based on the O2 -dependent NO-dioxygenase (NOD) activity of truncated hemoglobin 2/2HbN (Mt2/2HbN). The NOD activity largely depends on the efficiency of ligand migration to the heme cavity through a two-tunnel (long and short) system; recently, it was also correlated with the presence at the Mt2/2HbN N-terminus of a short pre-A region, not conserved in most 2/2HbNs, whose deletion results in a drastic reduction of NO scavenging. In the present study, we report the crystal structure of Mt2/2HbN-ΔpreA, lacking the pre-A region, at a resolution of 1.53 Å. We show that removal of the pre-A region results in long range effects on the protein C-terminus, promoting the assembly of a stable dimer, both in the crystals and in solution. In the Mt2/2HbN-ΔpreA dimer, access of heme ligands to the short tunnel is hindered. Molecular dynamics simulations show that the long tunnel branch is the only accessible pathway for O2 -ligand migration to/from the heme, and that the gating residue Phe(62)E15 partly restricts the diameter of the tunnel. Accordingly, kinetic measurements indicate that the kon value for peroxynitrite isomerization by Mt2/2HbN-ΔpreA-Fe(III) is four-fold lower relative to the full-length protein, and that NO scavenging by Mt2/2HbN-ΔpreA-Fe(II)-O2 is reduced by 35-fold. Therefore, we speculate that Mt2/2HbN evolved to host the pre-A region as a mechanism for preventing dimerization, thus reinforcing the survival of the microorganism against the reactive nitrosative stress in macrophages. DATABASE: Coordinates and structure factors have been deposited in the Protein Data Bank under accession number 5AB8.

Molecular model of hemoglobin N from Mycobacterium tuberculosis bound to lipid bilayers: a combined spectroscopic and computational study.

A singular aspect of the 2-on-2 hemoglobin structures of groups I and II is the presence of tunnels linking the protein surface to the distal heme pocket, supporting the storage and the diffusion of small apolar ligands to/from the buried active site. As the solubility of apolar ligands is greater in biological membranes than in solution, the association of these proteins with biological membranes may improve the efficiency of ligand capture. As very little is known on this subject, we have investigated the interactions between hemoglobin N (HbN), a group I 2-on-2 hemoglobin from the pathogenic Mycobacterium tuberculosis (Mtb), and biological membranes using both experimental techniques and MD simulations. HbN has a potent nitric oxide dioxygenase activity (HbN-Fe²âº-O2 + •NO + H2O → HbN-Fe³âº-OH2 + NO3⁻) that is thought to protect the aerobic respiration of Mtb from inhibition by •NO. Three different membrane compositions were chosen for the studies, representative of the mycobacterial plasma membrane and the mammalian cell membranes. Both the experimental and the modeling results agreed with each other and allow for a detailed molecular description of HbN in association with membranes of different compositions. The results indicated that HbN is a peripheral protein, and the association with the membranes occurred via the pre-A, G, and H helices. In addition, HbN would be allowed to modulate the binding to the membranes via electrostatic interactions between the lipid membranes and the Asp100 residue. In its membrane-bound form the short tunnel of HbN is oriented toward the membrane interior and the other tunnels point toward the solvent. Such protein orientation would facilitate the uptake of nonpolar substrates from the membrane and the release of products to the solvent. It is interesting to note that the pre-A, G, and H helices are conserved among HbN from a few other Mycobacteria.

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations.

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe²âº-O2 + (•)NO → trHbN-Fe³âº-OH2 + NO3⁻) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (•)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (•)NO and O2 to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.

Structural characterization of a group II 2/2 hemoglobin from the plant pathogen Agrobacterium tumefaciens.

Within the 2/2 hemoglobin sub-family, no group II 2/2Hbs from proteobacteria have been so far studied. Here we present the first structural characterization of a group II 2/2Hb from the soil and phytopathogenic bacterium Agrobacterium tumefaciens (At-2/2HbO). The crystal structure of ferric At-2/2HbO (reported at 2.1Å resolution) shows the location of specific/unique heme distal site residues (e.g., His(42)CD1, a residue distinctive of proteobacteria group II 2/2Hbs) that surround a heme-liganded water molecule. A highly intertwined hydrogen-bonded network, involving residues Tyr(26)B10, His(42)CD1, Ser(49)E7, Trp(93)G8, and three distal site water molecules, stabilizes the heme-bound ligand. Such a structural organization suggests a path for diatomic ligand diffusion to/from the heme. Neither a similar distal site structuring effect nor the presence of distal site water molecules has been so far observed in group I and group III 2/2Hbs, thus adding new distinctive information to the complex picture of currently available 2/2Hb structural and functional data. This article is part of a Special Issue entitled: Protein Structure and Function in the Crystalline State.

Theoretical investigations of nitric oxide channeling in Mycobacterium tuberculosis truncated hemoglobin N.

Mycobacterium tuberculosis group I truncated hemoglobin trHbN catalyzes the oxidation of nitric oxide (NO) to nitrate with a second-order rate constant k approximately 745 microM(-1) s(-1) at 23 degrees C (nitric oxide dioxygenase reaction). It was proposed that this high efficiency is associated with the presence of hydrophobic tunnels inside trHbN structure that allow substrate diffusion to the distal heme pocket. In this work, we investigated the mechanisms of NO diffusion within trHbN tunnels in the context of the nitric oxide dioxygenase reaction using two independent approaches. Molecular dynamics simulations of trHbN were performed in the presence of explicit NO molecules. Successful NO diffusion from the bulk solvent to the distal heme pocket was observed in all simulations performed. The simulations revealed that NO interacts with trHbN at specific surface sites, composed of hydrophobic residues located at tunnel entrances. The entry and the internal diffusion of NO inside trHbN were performed using the Long, Short, and EH tunnels reported earlier. The Short tunnel was preferentially used by NO to reach the distal heme pocket. This preference is ascribed to its hydrophobic funnel-shape entrance, covering a large area extending far from the tunnel entrance. This funnel-shape entrance triggers the frequent formation of solvent-excluded cavities capable of hosting up to three NO molecules, thereby accelerating NO capture and entry. The importance of hydrophobicity of entrances for NO capture is highlighted by a comparison with a polar mutant for which residues at entrances were mutated with polar residues. A complete map of NO diffusion pathways inside trHbN matrix was calculated, and NO molecules were found to diffuse from Xe cavity to Xe cavity. This scheme was in perfect agreement with the three-dimensional free-energy distribution calculated using implicit ligand sampling. The trajectories showed that NO significantly alters the dynamics of the key amino acids of Phe(62)(E15), a residue proposed to act as a gate controlling ligand traffic inside the Long tunnel, and also of Ile(119)(H11), at the entrance of the Short tunnel. It is noteworthy that NO diffusion inside trHbN tunnels is much faster than that reported previously for myoglobin. The results presented in this work shed light on the diffusion mechanism of apolar gaseous substrates inside protein matrix.

Structural characterization of the tunnels of Mycobacterium tuberculosis truncated hemoglobin N from molecular dynamics simulations.

The structure of oxygenated trHbN from Mycobacterium tuberculosis shows an extended heme distal hydrogen-bond network that includes Tyr33(B10), Gln58(E11), and the bound O(2). In addition, trHbN structure shows a network of hydrophobic cavities organized in two orthogonal branches. In the present work, the structure and the dynamics of oxygenated and deoxygenated trHbN in explicit water was investigated from 100 ns molecular dynamics (MD) simulations. Results show that, depending on the presence or the absence of a coordinated O(2), the Tyr33(B10) and Gln58(E11) side chains adopt two different configurations in concert with hydrogen bond network rearrangement. In addition, our data indicate that Tyr33(B10) and Gln58(E11) control the dynamics of Phe62(E15). In deoxy-trHbN, Phe62(E15) is restricted to one conformation. Upon O(2) binding, the conformation of Gln58(E11) changes and residue Phe62(E15) fluctuates between two conformations. We also conducted a systematic study of trHbN tunnels by analyzing thousands of MD snapshots with CAVER. The results show that tunnel formation is the result of the dynamic reshaping of short-lived hydrophobic cavities. The analyses indicate that the presence of these cavities is likely linked to the rigid structure of trHbN and also reveal two tunnels, EH and BE, that link the protein surface to the buried distal heme pocket and not present in the crystallographic structure. The cavities are sufficiently large to accomodate and store ligands. Tunnel dynamics in trHbN was found to be controlled by the side-chain conformation of the Tyr33(B10), Gln58(E11), and Phe62(E15) residues. Importantly, in contrast to recently published works, our extensive systematic studies show that the presence or absence of a coordinated dioxygen does not control the opening of the long tunnel but rather the opening of the EH tunnel. In addition, the data lead to new and distinctly different conclusion on the impact of the Phe62(E15) residue on trHbN tunnels. We propose that the EH and the long tunnels are used for apolar ligands storage. The trajectories bring important new structural insights related to trHbN function and to ligand diffusion in proteins.

Ligand binding to truncated hemoglobin N from Mycobacterium tuberculosis is strongly modulated by the interplay between the distal heme pocket residues and internal water.

The survival of Mycobacterium tuberculosis requires detoxification of host *NO. Oxygenated Mycobacterium tuberculosis truncated hemoglobin N catalyzes the rapid oxidation of nitric oxide to innocuous nitrate with a second-order rate constant (k'(NOD) approximately 745 x 10(6) m(-1) x s(-1)), which is approximately 15-fold faster than the reaction of horse heart myoglobin. We ask what aspects of structure and/or dynamics give rise to this enhanced reactivity. A first step is to expose what controls ligand/substrate binding to the heme. We present evidence that the main barrier to ligand binding to deoxy-truncated hemoglobin N (deoxy-trHbN) is the displacement of a distal cavity water molecule, which is mainly stabilized by residue Tyr(B10) but not coordinated to the heme iron. As observed in the Tyr(B10)/Gln(E11) apolar mutants, once this kinetic barrier is lowered, CO and O(2) binding is very rapid with rates approaching 1-2 x 10(9) m(-1) x s(-1). These large values almost certainly represent the upper limit for ligand binding to a heme protein and also indicate that the iron atom in trHbN is highly reactive. Kinetic measurements on the photoproduct of the *NO derivative of met-trHbN, where both the *NO and water can be directly followed, revealed that water rebinding is quite fast (approximately 1.49 x 10(8) s(-1)) and is responsible for the low geminate yield in trHbN. Molecular dynamics simulations, performed with trHbN and its distal mutants, indicated that in the absence of a distal water molecule, ligand access to the heme iron is not hindered. They also showed that a water molecule is stabilized next to the heme iron through hydrogen-bonding with Tyr(B10) and Gln(E11).

Ferrous Campylobacter jejuni truncated hemoglobin P displays an extremely high reactivity for cyanide - a comparative study.

Campylobacter jejuni hosts two hemoglobins (Hbs). The Camplylobacter jejuni single-domain Hb (called Cgb) is homologous to the globin domain of flavohemoglobin, and it has been proposed to protect the bacterium against nitrosative stress. The second Hb is called Ctb (hereafter Cj-trHbP), belongs to truncated Hb group III, and has been hypothesized to be involved in O(2) chemistry. Here, the kinetics and thermodynamics of cyanide binding to ferric and ferrous Cj-trHbP [Cj-trHbP(III) and Cj-trHbP(II), respectively] are reported and analyzed in parallel with those of related heme proteins, with particular reference to those from Mycobacterium tuberculosis. The affinity of cyanide for Cj-trHbP(II) is higher than that reported for any known (in)vertebrate globin by more than three orders of magnitude (K = 1.2 x 10(-6) m). This can be fully attributed to the highest (ever observed for a ferrous Hb) cyanide-binding association rate constant (k(on) = 3.3 x 10(3) m(-1).s(-1)), even though the binding process displays a rate-limiting step (k(max) = 9.1 s(-1)). Cj-trHbP(III) shows a very high affinity for cyanide (L = 5.8 x 10(-9) m); however, cyanide association kinetics are independent of cyanide concentration, displaying a rate-limiting step (l(max) = 2.0 x 10(-3) s(-1)). Values of the first-order rate constant for cyanide dissociation from Cj-trHbP(II)-cyanide and Cj-trHbP(III)-cyanide (k(off) =5.0 x 10(-3) s(-1) and l(off) > or = 1 x 10(-4) s(-1), respectively) are similar to those reported for (in)vertebrate globins. The very high affinity of cyanide for Cj-trHbP(II), reminiscent of that of horseradish peroxidase(II), suggests that this globin may participate in cyanide detoxification.

The roles of Tyr(CD1) and Trp(G8) in Mycobacterium tuberculosis truncated hemoglobin O in ligand binding and on the heme distal site architecture.

The crystal structure of the cyano-met form of Mt-trHbO revealed two unusual distal residues Y(CD1) and W(G8) forming a hydrogen-bond network with the heme-bound ligand [Milani, M., et al. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 5766-5771]. W(G8) is an invariant residue in group II and group III trHbs and has no counterpart in other globins. A previous study reported that changing Y(CD1) for a Phe causes a significant increase in the O2 combination rate, but almost no change in the O2 dissociation rate [Ouellet, H., et al. (2003) Biochemistry 42, 5764-5774]. Here we investigated the role of the W(G8) in ligand binding by using resonance Raman spectroscopy, stopped-flow spectrophotometry, and X-ray crystallography. For this purpose, W(G8) was changed, by site-directed mutagenesis, to a Phe in both the wild-type protein and the mutant Y(CD1)F to create the single mutant W(G8)F and the double mutant Y(CD1)F/W(G8)F, respectively. Resonance Raman results suggest that W(G8) interacts with the heme-bound O2 and CO, as evidenced by the increase of the Fe-O2 stretching mode from 559 to 564 cm-1 and by the lower frequency of the Fe-CO stretching modes (514 and 497 cm-1) compared to that of the wild-type protein. Mutation of W(G8) to Phe indicates that this residue controls ligand binding, as evidenced by a dramatic increase of the combination rates of both O2 and CO. Also, the rate of O2 dissociation showed a 90-1000-fold increase in the W(G8)F and Y(CD1)F/W(G8)F mutants, that is in sharp contrast with the values obtained for the other distal mutants Y(B10)F and Y(CD1)F [Ouellet, H., et al. (2003) Biochemistry 42, 5764-5774]. Taken together, these data indicate a pivotal role for the W(G8) residue in O2 binding and stabilization.

Mycobacterial truncated hemoglobins: from genes to functions.

Infections caused by bacteria belonging to genus Mycobacterium are among the most challenging threats for human health. The ability of mycobacteria to persist in vivo in the presence of reactive nitrogen and oxygen species implies the presence in these bacteria of effective detoxification mechanisms. Mycobacterial truncated hemoglobins (trHbs) have recently been implicated in scavenging of reactive nitrogen species. Individual members from each trHb family (N, O, and P) can be present in the same mycobacterial species. The distinct features of the heme active site structure combined with different ligand binding properties and in vivo expression patterns of mycobacterial trHbs suggest that these globins may accomplish diverse functions. Here, recent genomic, structural and biochemical information on mycobacterial trHbs is reviewed, with the aim of providing further insights into the role of these globins in mycobacterial physiology.

Reaction of Mycobacterium tuberculosis truncated hemoglobin O with hydrogen peroxide: evidence for peroxidatic activity and formation of protein-based radicals.

In this work, we investigated the reaction of ferric Mycobacterium tuberculosis truncated hemoglobin O (trHbO) with hydrogen peroxide. Stopped-flow spectrophotometric experiments under single turnover conditions showed that trHbO reacts with H(2)O(2) to give transient intermediate(s), among which is an oxyferryl heme, different from a typical peroxidase Compound I (oxyferryl heme pi-cation radical). EPR spectroscopy indicated evidence for both tryptophanyl and tyrosyl radicals, whereas redox titrations demonstrated that the peroxide-treated protein product retains 2 oxidizing eq. We propose that Compound I formed transiently is reduced with concomitant oxidation of Trp(G8) to give the detected oxoferryl heme and a radical on Trp(G8) (detected by EPR of the trHbO Tyr(CD1)Phe mutant). In the wild-type protein, the Trp(G8) radical is in turn reduced rapidly by Tyr(CD1). In a second cycle, Trp(G8) may be reoxidized by the ferryl heme to yield ferric heme and two protein radicals. In turn, these migrate to form tyrosyl radicals on Tyr(55) and Tyr(115), which lead, in the absence of a reducing substrate, to oligomerization of the protein. Steady-state kinetics in the presence of H(2)O(2) and the one-electron donor 2,2'-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) indicated that trHbO has peroxidase activity, in accord with the presence of typical peroxidase intermediates. These findings suggest an oxidation/reduction function for trHbO and, by analogy, for other Group II trHbs.

Structural determinants in the group III truncated hemoglobin from Campylobacter jejuni.

Truncated hemoglobins (trHbs) constitute a distinct lineage in the globin superfamily, distantly related in size and fold to myoglobin and monomeric hemoglobins. Their phylogenetic analyses revealed that three groups (I, II, and III) compose the trHb family. Group I and II trHbs adopt a simplified globin fold, essentially composed of a 2-on-2 alpha-helical sandwich, wrapped around the heme group. So far no structural data have been reported for group III trHbs. Here we report the three-dimensional structure of the group III trHbP from the eubacterium Campylobacter jejuni. The 2.15-A resolution crystal structure of C. jejuni trHbP (cyano-met form) shows that the 2-on-2 trHb fold is substantially conserved in the trHb group III, despite the absence of the Gly-based sequence motifs that were considered necessary for the attainment of the trHb specific fold. The heme crevice presents important structural modifications in the C-E region and in the FG helical hinge, with novel surface clefts at the proximal heme site. Contrary to what has been observed for group I and II trHbs, no protein matrix tunnel/cavity system is evident in C. jejuni trHbP. A gating movement of His(E7) side chain (found in two alternate conformations in the crystal structure) may be instrumental for ligand entry to the heme distal site. Sequence conservation allows extrapolating part of the structural results here reported to the whole trHb group III.

Ligand interactions in the distal heme pocket of Mycobacterium tuberculosis truncated hemoglobin N: roles of TyrB10 and GlnE11 residues.

The crystallographic structure of oxygenated trHbN from Mycobacterium tuberculosis showed an extended heme distal site hydrogen-bonding network that includes Y(B10), Q(E11), and the bound O(2) (Milani, M., et al. (2001) EMBO J. 20, 3902-3909). In the present work, we analyze the effects that substitutions at the B10 and E11 positions exert on the heme and its coordinated ligands, using steady-state resonance Raman spectroscopy, absorption spectroscopy and X-ray crystallography. Our results show that (1) residues Y(B10) and Q(E11) control the binding and the ionization state of the heme-bound water molecules in ferric trHbN and are important in keeping the sixth coordination position vacant in deoxy trHbN; (2) residue Q(E11) plays a role in maintaining the integrity of the proximal Fe-His bond in deoxy trHbN; (3) in wild-type oxy-trHbN, the size and hydrogen-bonding capability of residue E11 is important to sustain proper interaction between Y(B10) and the heme-bound O(2); (4) CO-trHbN is in a conformational equilibrium, where either the Y(B10) or the Q(E11) residue interacts with the heme-bound CO; and (5) Y(B10) and Q(E11) residues control the conformation (and likely the dynamics) of the protein matrix tunnel gating residue F(E15). These findings suggest that the functional processes of ligand binding and diffusion are controlled in trHbN through the dynamic interaction of residues Y(B10), Q(E11), F(E15), and the heme ligand.

Three globin lineages belonging to two structural classes in genomes from the three kingdoms of life.

Although most globins, including the N-terminal domains within chimeric proteins such as flavohemoglobins and globin-coupled sensors, exhibit a 3/3 helical sandwich structure, many bacterial, plant, and ciliate globins have a 2/2 helical sandwich structure. We carried out a comprehensive survey of globins in the genomes from the three kingdoms of life. Bayesian phylogenetic trees based on manually aligned sequences indicate the possibility of past horizontal globin gene transfers from bacteria to eukaryotes. blastp searches revealed the presence of 3/3 single-domain globins related to the globin domains of the bacterial and fungal flavohemoglobins in many bacteria, a red alga, and a diatom. Iterated psi-blast searches based on groups of globin sequences found that only the single-domain globins and flavohemoglobins recognize the eukaryote 3/3 globins, including vertebrate neuroglobins, alpha- and beta-globins, and cytoglobins. The 2/2 globins recognize the flavohemoglobins, as do the globin coupled sensors and the closely related single-domain protoglobins. However, the 2/2 globins and the globin-coupled sensors do not recognize each other. Thus, all globins appear to be distributed among three lineages: (i) the 3/3 plant and metazoan globins, single-domain globins, and flavohemoglobins; (ii) the bacterial 3/3 globin-coupled sensors and protoglobins; and (iii) the bacterial, plant, and ciliate 2/2 globins. The three lineages may have evolved from an ancestral 3/3 or 2/2 globin. Furthermore, it appears likely that the predominant functions of globins are enzymatic and that oxygen transport is a specialized development that accompanied the evolution of metazoans.

Structural bases for heme binding and diatomic ligand recognition in truncated hemoglobins.

Truncated hemoglobins (trHbs) are low-molecular-weight oxygen-binding heme-proteins distributed in eubacteria, cyanobacteria, unicellular eukaryotes, and in higher plants, constituting a distinct group within the hemoglobin (Hb) superfamily. TrHbs display amino acid sequences 20-40 residues shorter than classical (non)vertebrate Hbs and myoglobins, to which they are scarcely related by sequence similarity. The trHb tertiary structure is based on a 2-on-2 alpha-helical sandwich, which represents a striking editing of the highly conserved 3-on-3 alpha-helical globin fold, achieved through deletion/truncation of alpha-helices and specific residue substitutions. Despite their 'minimal' polypeptide chain span, trHbs display an inner tunnel/cavity system held to support ligand diffusion to/from the heme distal pocket, accumulation of heme ligands within the protein matrix, and/or multiligand reactions. Moreover, trHbs bind and effectively stabilize the heme and recognize diatomic ligands (i.e., O2, CO, NO, and cyanide), albeit with varying thermodynamic and kinetic parameters. Here, structural bases for heme binding and diatomic ligand recognition by trHbs are reviewed.

Viscosity-dependent relaxation significantly modulates the kinetics of CO recombination in the truncated hemoglobin TrHbN from Mycobacterium tuberculosis.

Kinetic traces were generated for the nanosecond and slower rebinding of photodissociated CO to trHbN in solution and in porous sol-gel matrices as a function of viscosity, conformation, and mutation. TrHbN is one of the two truncated hemoglobins from Mycobacterium tuberculosis. The kinetic traces were analyzed in terms of three distinct phases. These three phases are ascribed to rebinding: (i) from the distal heme pocket, (ii) from the adjacent apolar tunnel prior to conformational relaxation, and (iii) from the apolar tunnel subsequent to conformational relaxation. The fractional content of each of these phases was shown to be a function of the viscosity and, in the case of the sol-gel-encapsulated samples, sample preparation history. The observed kinetic patterns support a model consisting of the following elements: (i) the viscosity and conformation-sensitive dynamics of the Tyr(B10) side chain facilitate diffusion of the dissociated ligand from the distal heme pocket into the adjacent tunnel; (ii) the distal heme pocket architecture determines ligand access from the tunnel back to the heme iron; (iii) the distal heme pocket architecture is governed by a ligand-dependent hydrogen bonding network that limits the range of accessible side chain positions; and (iv) the apolar tunnel linking the heme site to the solvent biases the competition between water and ligand for occupancy of the vacated polar distal heme pocket greatly toward the nonpolar ligand. Implications of these finding with respect to biological function are discussed.

Cyanide binding to truncated hemoglobins: a crystallographic and kinetic study.

Cyanide is one of the few diatomic ligands able to interact with the ferric and ferrous heme-Fe atom. Here, the X-ray crystal structure of the cyanide derivative of ferric Mycobacterium tuberculosis truncated hemoglobin-N (M. tuberculosis trHbN) has been determined at 2.0 A (R-general = 17.8% and R-free = 23.5%), and analyzed in parallel with those of M. tuberculosis truncated hemoglobin-O (M. tuberculosis trHbO), Chlamydomonas eugametos truncated hemoglobin (C. eugametos trHb), and sperm whale myoglobin, generally taken as a molecular model. Cyanide binding to M. tuberculosis trHbN is stabilized directly by residue TyrB10(33), which may assist the deprotonation of the incoming ligand and the protonation of the outcoming cyanide. In M. tuberculosis trHbO and in C. eugametos trHb the ligand is stabilized by the distal pocket residues TyrCD1(36) and TrpG8(88), and by the TyrB10(20) - GlnE7(41) - GlnE11(45) triad, respectively. Moreover, kinetics for cyanide binding to ferric M. tuberculosis trHbN and trHbO and C. eugametos trHb, for ligand dissociation from the ferrous trHbs, and for the reduction of the heme-Fe(III)-cyanide complex have been determined, at pH 7.0 and 20.0 degrees C. Despite the different heme distal site structures and ligand interactions, values of the rate constant for cyanide binding to ferric (non)vertebrate heme proteins are similar, being influenced mainly by the presence in the heme pocket of proton acceptor group(s), whose function is to assist the deprotonation of the incoming ligand (i.e., HCN). On the other hand, values of the rate constant for the reduction of the heme-Fe(III)-cyanide (non)vertebrate globins span over several orders of magnitude, reflecting the different ability of the heme proteins considered to give productive complex(es) with dithionite or its reducing species SO(2)(-). Furthermore, values of the rate constant for ligand dissociation from heme-Fe(II)-cyanide (non)vertebrate heme proteins are very different, reflecting the different nature and geometry of the heme distal residue(s) hydrogen-bonded to the heme-bound cyanide.