Kinetic modulation in carbonmonoxy derivatives of truncated hemoglobins: the role of distal heme pocket residues and extended apolar tunnel

Myxococcus xanthus truncated globin HbO: in silico analysis and functional characterization

Truncated globins are 20-40 amino acids shorter than full length globins. Till date, globins have been characterized predominantly from bacteria involved in pathogenicity, nitrogen fixation and photosynthesis, where they are implicated in bacterial virulence within the host, protection of nitrogenase from oxygen inactivation and prevention of oxidative damage to the photosynthetic machinery respectively. Myxococcus xanthus, the model myxobacterium, is an obligate aerobe with a multicellular stage in its life cycle where cells encounter oxygen limitation. This work was undertaken to investigate the potential role of the truncated globin in M. xanthus. To examine the role of globins in this unique group of bacteria, the gene coding for a putative truncated globin (HbO) was identified in the genome of M. xanthus DK 1622. The sequence analysis by bioinformatics approaches revealed that HbO from M. xanthus (Mx-HbO) likely adopts a 2-on-2 alpha helical fold of the truncated globins. The gene coding for Mx-HbO was cloned and its expression in E. coli imparted reddish tinge to the cells. The spectral analysis confirmed it to be a functional globin. The expression of Mx-HbO in the heterologous host improved its growth, resulting in the attainment of higher cell density in culture. The transcript of Mx-hbO was induced threefold in the host cells when grown under low aeration condition as compared to the cells grown under high aeration condition. In M. xanthus, an obligate aerobe, where cell growth accompanies swarming, there is a higher density of cells in the middle of the swarm. Our results suggest that Mx-HbO is a functional globin and could facilitate the growth of cells facing oxygen deprivation, the condition prevailing in the middle of the swarm.

Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study

Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.

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.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Ligand migration and hexacoordination in type 1 non-symbiotic rice hemoglobin

Type 1 non-symbiotic rice hemoglobin (rHb1) shows bis-histidyl heme hexacoordination and is capable of binding diatomic ligands reversibly. The biological function is as yet unclear, but the high oxygen affinity makes it unlikely to be involved in oxygen transport. In order to gain insight into possible physiological roles, we have studied CO rebinding kinetics after laser flash photolysis of rHb1 in solution and encapsulated in silica gel. CO rebinding to wt rHb1 in solution occurs through a fast geminate phase with no sign of rebinding from internal docking sites. Encapsulation in silica gel enhances migration to internal cavities. Site-directed mutagenesis of FB10, a residue known to have a key role in the regulation of hexacoordination and ligand affinity, resulted in substantial effects on the rebinding kinetics, partly inhibiting ligand exit to the solvent, enhancing geminate rebinding and enabling ligand migration within the internal cavities. The mutation of HE7, one of the histidyl residues involved in the hexacoordination, prevents hexacoordination, as expected, but also exposes ligand migration through a complex system of cavities. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

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.

Ligand migration through the internal hydrophobic cavities in human neuroglobin

Neuroglobin (Ngb), a member of the globin superfamily, was found in the brain of vertebrates and is suggested to play a neuroprotective function under hypoxic conditions by scavenging nitrogen monoxide (NO) through a dioxygenase activity. In order for such a reaction to efficiently take place and to minimize the release of reactive intermediates in the cytosol, the cosubstrates O(2) and NO and other unstable reaction intermediates should bind sequentially to docking sites in the protein matrix. We have characterized the accessibility of these sites by analyzing the geminate CO rebinding kinetics to the heme moiety observed upon nanosecond flash photolysis of the Ngb-CO complex encapsulated in silica gels. The geminate rebinding phase showed a remarkable complexity, revealing the presence of a system of secondary docking sites where ligands are stored for hundreds of microseconds. Most kinetics steps display little temperature dependence, demonstrating that ligands can easily migrate through the cavities, except for the slowest reaction intermediate, possibly reflecting a structural conformational change reshaping the system of cavities. This conformational change is unrelated with distal His E7 binding to the heme, as it persists for the HE7L mutant. Overall, data are consistent with the presence of a discrete system of docking sites, possibly acting as reservoirs for the putative cosubstrates and for other reactive species involved in the physiologically relevant reaction.

Intrinsic non-symbiotic and truncated haemoglobins and heterologous Vitreoscilla haemoglobin expression in plants

To date, haemoglobins (Hbs) have been shown to exist in all kingdoms of life. The least studied and understood groups are plant non-symbiotic haemoglobins (nsHbs) and the recently found plant truncated Hbs (trHbs). From a biotechnological point of view, the best characterized and almost exclusively applied Hb is the bacterial Vitreoscilla haemoglobin (VHb). In this review, the present state of knowledge of structural features and ligand binding kinetics of plant nsHbs and trHbs and their proposed roles as oxygen carriers, oxygen sensors, and for oxygen storage, in nitric oxide (NO) detoxification, and in peroxidase activity are described. Furthermore, in order to predict the functioning of plant Hbs, their characteristics will be compared with those of the better known bacterial globins. In this context, the effects of heterologous applications of VHb on plants are reviewed. Finally, the challenging future of plant Hb research is discussed.

Effects of active site mutations in haemoglobin I from Lucina pectinata: a molecular dynamic study

Haemoglobin I from Lucina pectinata is a monomeric protein consisting of 142 amino acids. Its active site contains a peculiar arrangement of phenylalanine residues (PheB10, PheCD1 and PheE11) and a distal Gln at position E7. Active site mutations at positions B10, E7 and E11 were performed in deoxy haemoglobin I (HbI), followed by 10 ns molecular dynamic simulations. The results showed that the mutations induced changes in domains far from the active site producing more flexible structures than the native HbI. Distance analyses revealed that the heme pocket amino acids at positions E7 and B10 are extremely sensitive to any heme pocket residue mutation. The high flexibility observed by the E7 position suggests an important role in the ligand binding kinetics in ferrous HbI, while both positions play a major role in the ligand stabilisation processes. Furthermore, our results showed that E11Phe plays a pivotal role in protein stability.

Characterization of ligand migration mechanisms inside hemoglobins from the analysis of geminate rebinding kinetics

The presence of internal hydrophobic cavities and packing defects has been demonstrated for several small globular proteins, including hemoglobins. The reduced thermodynamic stability appears to be compensated for by the capability of controlling ligand diffusion through the protein matrix to the active site, possibly by stocking more than one reactant molecule in selected sites. Photolysis of carbon monoxide complexes of hemoglobins encapsulated in silica gels leads to multiphasic geminate rebinding kinetics at room temperature, reflecting rebinding also from different temporary docking sites inside the protein matrix. A careful analysis of the ligand rebinding kinetics allows the determination of the microscopic rates for the underlying reactions, including those governing the migration to and from the docking sites. This chapter describes the experimental approach used to characterize the ligand rebinding kinetics for heme proteins in silica gels after nanosecond laser flash photolysis and the computational methods necessary to retrieve the kinetic parameters.

The structure and function of plant hemoglobins

Plants, like humans, contain hemoglobin. Three distinct types of hemoglobin exist in plants: symbiotic, non-symbiotic, and truncated hemoglobins. Crystal structures and other structural and biophysical techniques have revealed important knowledge about ligand binding and conformational stabilization in all three types. In symbiotic hemoglobins (leghemoglobins), ligand binding regulatory mechanisms have been shown to differ dramatically from myoglobin and red blood cell hemoglobin. In the non-symbiotic hemoglobins found in all plants, crystal structures and vibrational spectroscopy have revealed the nature of the structural transition between the hexacoordinate and ligand-bound states. In truncated hemoglobins, the abbreviated globin is porous, providing tunnels that may assist in ligand binding, and the bound ligand is stabilized by more than one distal pocket residue. Research has implicated these plant hemoglobins in a number of possible functions differing among hemoglobin types, and possibly between plant species.

Identification of ligand-binding pathways in truncated hemoglobins using locally enhanced sampling molecular dynamics

This chapter reviews the use of a locally enhanced sampling molecular dynamics (LESMD) for the study of ligand binding in truncated hemoglobins. The method, however, can be applied to any protein-ligand system. Truncated hemoglobins appear to have a tunnel(s) potentially used by the ligand to bind. These structural features give some indication of how the ligand moves through the protein to bind but do not give the complete picture. The LESMD method has been used to investigate the pathways of ligand binding to group I truncated hemoglobins from the eubacteria Mycobacterium tuberculosis, the ciliated protozoa Paramecium caudatum, and the unicellular alga Chlamydomonas eugametos.

Ligand migration in nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana

AHb1 is a hexacoordinated type 1 nonsymbiotic hemoglobin recently discovered in Arabidopsis thaliana. To gain insight into the ligand migration inside the protein, we studied the CO rebinding kinetics of AHb1 encapsulated in silica gels, in the presence of glycerol. The CO rebinding kinetics after nanosecond laser flash photolysis exhibits complex ligand migration patterns, consistent with the existence of discrete docking sites in which ligands can temporarily be stored before rebinding to the heme at different times. This finding may be of relevance to the physiological NO dioxygenase activity of this protein, which requires sequential binding of two substrates, NO and O2, to the heme.

Ligand recombination and a hierarchy of solvent slaved dynamics: the origin of kinetic phases in hemeproteins

Ligand recombination studies play a central role both for characterizing different hemeproteins and their conformational states but also for probing fundamental biophysical processes. Consequently, there is great importance to providing a foundation from which one can understand the physical processes that give rise to and modulate the large range of kinetic patterns associated with ligand recombination in myoglobins and hemoglobins. In this work, an overview of cryogenic and solution phase recombination phenomena for COMb is first reviewed and then a new paradigm is presented for analyzing the temperature and viscosity dependent features of kinetic traces in terms of multiple phases that reflect which tier(s) of solvent slaved protein dynamics is (are) operative on the photoproduct population during the time course of the measurement. This approach allows for facile inclusion of both ligand diffusion among accessible cavities and conformational relaxation effects. The concepts are illustrated using kinetic traces and MEM populations derived from the CO recombination process for wild type and mutant myoglobins either in sol-gel matrices bathed in glycerol or in trehalose-derived glassy matrices.

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.

Dynamical Regulation of Ligand Migration by a Gate-Opening Molecular Switch in Truncated Hemoglobin-N fromMycobacterium tuberculosis

Truncated hemoglobin-N is believed to constitute a defense mechanism of Mycobacterium tuberculosis against NO produced by macrophages, which is converted to the harmless nitrate anion. This process is catalyzed very efficiently, as the enzyme activity is limited by ligand diffusion. By using extended molecular dynamics simulations we explore the mechanism that regulates ligand diffusion and, particularly, the role played by residues that assist binding of O2 to the heme group. Our data strongly support the hypothesis that the access of NO to the heme cavity is dynamically regulated by the TyrB10-GlnE11 pair, which acts as a molecular switch that controls opening of the ligand diffusion tunnel. Binding of O2 to the heme group triggers local conformational changes in the TyrB10-GlnE11 pair, which favor opening of the PheE15 gate residue through global changes in the essential motions of the protein skeleton. The complex pattern of conformational changes triggered upon O2 binding is drastically altered in the GlnE11-->Ala and TyrB10-->Phe mutants, which justifies the poor enzymatic activity observed experimentally for the TyrB10-->Phe form. The results support a molecular mechanism evolved to ensure access of NO to the heme cavity in the oxygenated form of the protein, which should warrant survival of the microorganism under stress conditions.

Protein fold and structure in the truncated (2/2) globin family

Analysis of amino acids sequences and protein folds has recently unraveled the structural bases and details of several proteins from the recently discovered "truncated hemoglobin" family. The analysis here presented, in agreement with previous surveys, shows that truncated hemoglobins can be classified in three main groups, based on their structural properties. Crystallographic analyses have shown that all three groups adopt a 2-on-2 alpha-helical sandwich fold, resulting from apparent editing of the classical 3-on-3 alpha-helical sandwich of vertebrate and invertebrate conventional globins. Specific structural features distinguish each of the three groups. Among these, a protein matrix tunnel system is typical of group I, a Trp residue at the G8 topological site is conserved in groups II and III, and TyrB10 is almost invariant through the three groups. A strongly intertwined network of hydrogen bonds stabilizes the heme bound ligand, despite variability of the heme distal residues observed in the different proteins considered. Details of ligand recognition in the three groups are discussed at the light of residue conservation and of differing ligand diffusion pathways to the heme. Based on structural analyses of the family-specific fold, we endorse a recent proposal of leaving the "truncated hemoglobins" term, that does not represent properly the observed 2-on-2 alpha-helical sandwich fold, and adopting the simple "2/2Hb" term to concisely address this protein family.

Crystallographic evidence for dioxygen interactions with iron proteins

The interaction of dioxygen with iron plays a key role in many important biological processes, such as dioxygen transport in the bloodstream and the reduction of dioxygen by iron in respiration. However, the catalytic mechanisms employed, for example in ligand oxidation, are not fully understood at the current time despite intensive biochemical, spectroscopic and structural studies. This review outlines the structural evidence obtained by X-ray crystallographic methods for the nature of the interactions between dioxygen and the metal in iron-containing proteins. Proteins involved in iron transport or electron transfer are not included.

The reactivity with CO of AHb1 and AHb2 from Arabidopsis thaliana is controlled by the distal HisE7 and internal hydrophobic cavities

The nonsymbiotic hemoglobins, AHb1 and AHb2, have recently been isolated from Arabidopsis thaliana. Using steady-state and time-resolved spectroscopic methods, we show that Fe2+ AHb1 contains a mixture of penta- and hexacoordinated heme, while Fe2+ AHb2 is fully hexacoordinated. In the CO complexes, polar interactions and H-bonds with the ligand are stronger for AHb1 than for AHb2. The ligand binding kinetics are substantially different, reflecting the distribution between the penta- and hexacoordinated species, and indicate that protein dynamics and ligand migration pathways are very specific for each of the two proteins. In particular, a very small, non-exponential geminate rebinding observed in AHb1 suggests that the distal heme cavity is connected with the exterior by a relatively open channel. The large, temperature-dependent geminate rebinding observed for AHb2 implies a major role of protein dynamics in the ligand migration from the distal cavity to the solvent. The structures of AHb1 and AHb2, modeled on the basis of the homologous rice hemoglobin, exhibit a different cavity system that is fully compatible with the observed ligand binding kinetics. Overall, these kinetic and structural data are consistent with the putative NO-dioxygenase activity previously attributed to AHb1, whereas the role of AHb2 remains elusive.

Time-resolved methods in Biophysics. 2. Monitoring haem proteins at work with nanosecond laser flash photolysis

Haem proteins have long been the most studied proteins in biophysics, and have become paradigms for the characterization of fundamental biomolecular processes as ligand binding and regulatory conformational transitions. The presence of the haem prosthetic group, the absorbance spectrum of which has a ligation sensitive region conveniently located in the UV-visible range, has offered a powerful and sensitive tool for the investigation of molecular functions. The central Fe atom is capable of reversibly binding diatomic ligands, including O(2), CO, and NO. The Fe-ligand bond is photolabile, and a reactive unligated state can be transiently generated with a pulsed laser. The photodissociated ligands quickly rebind to the haem and the process can be monitored by transient absorbance methods. The ligand rebinding kinetics reflects protein dynamics and ligand migration within the protein inner cavities. The characterization of these processes was done in the past mainly by low temperature experiments. The use of silica gels to trap proteins allows the characterization of internal ligand dynamics at room temperature. In order to show the potential of the laser flash photolysis techniques, combined with modern numerical analysis methods, we report experiments conducted on two non-symbiotic haemoglobins from Arabidopsis thaliana. The comparison between time courses recorded on haemoglobins in solution and encapsulated in silica gels allows for the highlighting of different interplays between protein dynamics and ligand migration.

Ligand-induced dynamical regulation of NO conversion in Mycobacterium tuberculosis truncated hemoglobin-N

Mycobacterium tuberculosis, the causative agent of human tuberculosis, is forced into latency by nitric oxide produced by macrophages during infection. In response to nitrosative stress M. tuberculosis has evolved a defense mechanism that relies on the oxygenated form of "truncated hemoglobin" N (trHbN), formally acting as NO-dioxygenase, yielding the harmless nitrate ion. X-ray crystal structures have shown that trHbN hosts a two-branched protein matrix tunnel system, proposed to control diatomic ligand migration to the heme, as the rate-limiting step in NO conversion to nitrate. Extended molecular dynamics simulations (0.1 micros), employed here to characterize the factors controlling diatomic ligand diffusion through the apolar tunnel system, suggest that O2 migration in deoxy-trHbN is restricted to a short branch of the tunnel, and that O2 binding to the heme drives conformational and dynamical fluctuations promoting NO migration through the long tunnel branch. The simulation results suggest that trHbN has evolved a dual-path mechanism for migration of O2 and NO to the heme, to achieve the most efficient NO detoxification.

Modulation of reactivity and conformation within the T-quaternary state of human hemoglobin: the combined use of mutagenesis and sol-gel encapsulation

A range of conformationally distinct functional states within the T quaternary state of hemoglobin are accessed and probed using a combination of mutagenesis and sol-gel encapsulation that greatly slow or eliminate the T --> R transition. Visible and UV resonance Raman spectroscopy are used to probe the proximal strain at the heme and the status of the alpha(1)beta(2) interface, respectively, whereas CO geminate and bimolecular recombination traces in conjunction with MEM (maximum entropy method) analysis of kinetic populations are used to identify functionally distinct T-state populations. The mutants used in this study are Hb(Nbeta102A) and the alpha99-alpha99 cross-linked derivative of Hb(Wbeta37E). The former mutant, which binds oxygen noncooperatively with very low affinity, is used to access low-affinity ligated T-state conformations, whereas the latter mutant is used to access the high-affinity end of the distribution of T-state conformations. A pattern emerges within the T state in which ligand reactivity increases as both the proximal strain and the alpha(1)beta(2) interface interactions are progressively lessened after ligand binding to the deoxy T-state species. The ligation and effector-dependent interplay between the heme environment and the stability of the Trp beta37 cluster in the hinge region of the alpha(1)beta(2) interface appears to determine the distribution of the ligated T-state species generated upon ligand binding. A qualitative model is presented, suggesting that different T quaternary structures modulate the stability of different alphabeta dimer conformations within the tetramer.

Modeling truncated hemoglobin vibrational dynamics

We present a study on the near equilibrium dynamics of two small proteins in the family of truncated hemoglobins, developed under the framework of a Gaussian network approach. Effective beta carbon atoms are taken into account besides Calphas for all residues but glycines in the coarse-graining procedure, without leading to an increase in the degrees of freedom (beta Gaussian Model). Normalized covariance matrix and deformation along slowest modes with collective character are analyzed, pointing out anticorrelations between functionally relevant sites for the proteins under study. In particular, we underline the functional motions of an extended tunnel-cavity system running inside the protein matrix, which provide a pathway for small ligands binding with the iron in the heme group. We give a rough estimate of the order of magnitude of the relaxation times of the slowest two overdamped modes and compare results with previous studies on globins.

The Position 68(E11) Side Chain in Myoglobin Regulates Ligand Capture, Bond Formation with Heme Iron, and Internal Movement into the Xenon Cavities

After photodissociation, ligand rebinding to myoglobin exhibits complex kinetic patterns associated with multiple first-order geminate recombination processes occurring within the protein and a simpler bimolecular phase representing second-order ligand rebinding from the solvent. A smooth transition from cryogenic-like to solution phase properties can be obtained by using a combination of sol-gel encapsulation, addition of glycerol as a bathing medium, and temperature tuning (-15 --> 65 degrees C). This approach was applied to a series of double mutants, myoglobin CO (H64L/V68X, where X = Ala, Val, Leu, Asn, and Phe), which were designed to examine the contributions of the position 68(E11) side chain to the appearance and disappearance of internal rebinding phases in the absence of steric and polar interactions with the distal histidine. Based on the effects of viscosity, temperature, and the stereochemistry of the E11 side chain, the three major phases, B --> A, C --> A, and D --> A, can be assigned, respectively, to ligand rebinding from the following: (i) the distal heme pocket, (ii) the xenon cavities prior to large amplitude side chain conformational relaxation, and (iii) the xenon cavities after significant conformational relaxation of the position 68(E11) side chain. The relative amplitudes of the B --> A and C --> A phases depend markedly on the size and shape of the E11 side chain, which regulates sterically both ligand return to the heme iron atom and ligand migration to the xenon cavities. The internal xenon cavities provide a transient docking site that allows side chain relaxations and the entry of water into the vacated distal pocket, which in turn slows ligand recombination markedly.

A hierarchy of functionally important relaxations within myoglobin based on solvent effects, mutations and kinetic model

Geminate CO rebinding in myoglobin is studied for two viscous solvents, trehalose and sol-gel (bathed in 100% glycerol) at several temperatures. Mutations in key distal hemepocket residues are used to eliminate or enhance specific relaxation modes. The time-resolved data are analyzed with a modified Agmon-Hopfield model which is capable of providing excellent fits in cases where a single relaxation mode is dominant. Using this approach, we determine the relaxation rate constants of specific functionally important modes, obtaining also their Arrhenius activation energies. We find a hierarchy of distal pocket modes controlling the rebinding kinetics. The "heme access mode" (HAM) is responsible for the major slow-down in rebinding. It is a solvent-coupled cooperative mode which restricts ligand return from the xenon cavities. Bulky side-chains, like those His64 and Trp29 (in the L29W mutant), operate like overdamped pendulums which move over and block the binding site. They may be either unslaved (His64) or moderately slaved (Trp29) to the solvent. Small side-chain relaxations, most notably of leucines, are revealed in some mutants (V68L, V68A). They are conjectured to facilitate inter-cavity ligand motion. When all relaxations are arrested (H64L in trehalose), we observe pure inhomogeneous kinetics with no temperature dependence, suggesting that proximal relaxation is not a factor on the investigated timescale.

Tyrosine B10 and heme-ligand interactions of Lucina pectinata hemoglobin II: control of heme reactivity

The distal pocket of hemoglobin II (HbII) from Lucina pectinata is characterized by the presence of a GlnE7 and a TyrB10. To elucidate the functional properties of HbII, biophysical studies were conducted on HbII and a HbI PheB10Tyr site-directed mutant. The pH titration data at neutral conditions showed visible bands at 486, 541, 577 and 605 nm for both proteins. This suggests the possible existence of a conformational equilibrium between an open and closed configuration due to the interactions of the TyrB10, ligand, and heme iron. The kinetic behavior for the reaction of both ferric proteins with H2O2 indicates that the rate for the formation of the ferryl intermediates species varies with pH, suggesting that the reaction is strongly dependent on the conformational states. At basic pH values, the barrier for the reaction increases as the tyrosine adopts a closed conformation and the ferric hydroxyl replaces the met-aquo species. The existence of these conformers is further supported by resonance Raman (RR) data, which indicate that in a neutral environment, the ferric HbII species is present as a possible mixture of coordination and spin states, with values at 1558 and 1580 cm(-1) for the nu2 marker, and 1479, 1492, and 1503 cm(-1) for the nu3 mode. Moreover, the presence of the A3 and A(o) conformers at 1924 and 1964 cm(-1) in the HbII-CO infrared spectra confirms the existence of an open and closed conformation due to the orientation of the TyrB10 with respect to the heme active center.

Spectroscopic and Functional Characterization of T State Hemoglobin Conformations Encapsulated in Silica Gels

Oxygen binding curves of sol-gel-encapsulated deoxy human adult hemoglobin (HbA) have previously revealed two distinct noncooperative populations with oxygen binding affinities approximately 1000 and 100 times lower than that of the high-affinity R state. The two populations which have been termed the low-affinity (LA) and high-affinity (HA) T states can be selectively stabilized using two different encapsulation protocols for deoxy-HbA. The present study seeks to understand the factors giving rise to these different affinity states. Visible and UV resonance Raman spectroscopies are used to characterize the conformational properties of both the deoxy and deoxy-turned-carbonmonoxy (CO) derivatives of HbA derived from the two encapsulation protocols. The geminate and bimolecular recombination of CO to the photodissociated CO derivatives is used to characterize the functional properties of the slowly evolving encapsulated populations. The results show that the initial deoxy-HbA populations are conformationally indistinguishable with respect to encapsulation protocol. The addition of CO to sol-gel-encapsulated deoxy-HbA triggers a detectable progression of conformational and functional changes. Visible resonance Raman spectra of the CO photoproduct reveal a progression of changes of the iron-proximal histidine stretching frequencies: 215, 222, 227, and 230 cm(-1). The low and high values correspond to the initial deoxy T state and liganded R (R(2)) state species, respectively. The 222 and 227 cm(-1) species are generated using encapsulation protocols that give rise to what are termed the LA and HA T states, respectively. The UV resonance Raman spectra of these and related species indicate that the progression from deoxy T to LA to HA is associated with a progressive loosening of T state constraints within the hinge and switch regions of the alpha(1)beta(2) interface. The time scale for the progression is determined by a balance between the ligation-initiated evolution toward high-affinity conformations and factors such as allosteric effectors, gel matrix, and added glycerol that slow ligand-binding-induced relaxation. Thus, it appears that the encapsulation protocol-dependent rate of ligand-binding-induced relaxation determines the functional properties of the initially encapsulated deoxy-HbA population.

Cyanide Binding to Hexacoordinate Cyanobacterial Hemoglobins: Hydrogen-Bonding Network and Heme Pocket Rearrangement in Ferric H117A Synechocystis Hemoglobin

The truncated hemoglobin (Hb) from the cyanobacterium Synechocystis sp. PCC 6803 is a bis-histidyl hexacoordinate complex in the absence of exogenous ligands. This protein can form a covalent cross-link between His117 in the H-helix and the heme 2-vinyl group. Cross-linking, the physiological importance of which has not been established, is avoided with the His117Ala substitution. In the present work, H117A Hb was used to explore exogenous ligand binding to the heme group. NMR and thermal denaturation data showed that the replacement was of little consequence to the structural and thermodynamic properties of ferric Synechocystis Hb. It did, however, decelerate the association of cyanide ions with the heme iron. Full complexation required hours, instead of minutes, of incubation at optical and NMR concentrations. At neutral pH and in the presence of excess cyanide, binding occurred with a first-order dependence on cyanide concentration, eliminating distal histidine decoordination as the rate-limiting step. The cyanide complex of the H117A variant was characterized for the conformational changes occurring as the histidine on the distal side, His46 (E10), was displaced. Extensive rearrangement allowed Tyr22 (B10) to insert in the heme pocket and Gln43 (E7) and Gln47 (E11) to come in contact with it. H-bond formation to the bound cyanide was identified in solution with the use of (1)H(2)O/(2)H(2)O mixtures. Cyanide binding also resulted in a change in the ratio of heme orientational isomers, in a likely manifestation of heme environment reshaping. Similar observations were made with the related Synechococcus sp. PCC 7002 H117A Hb, except that cyanide binding was rapid in this protein. In both cases, the (15)N chemical shift of bound cyanide was reminiscent of that in peroxidases and the orientation of the proximal histidine was as in other truncated Hbs. The ensemble of the data provided insight into the structural cooperativity of the heme pocket scaffold and pointed to the reactive 117 site of Synechocystis Hb as a potential determinant of biophysical and, perhaps, functional properties.

Neuroglobin and cytoglobin: Two new entries in the hemoglobin superfamily*

Neuroglobin (Ngb) and cytoglobin (Cygb) are two newly discovered intracellular members of the vertebrate hemoglobin (Hb) family. Ngb, predominantly expressed in nerve cells, is of ancient evolutionary origin and is homologous to nerve-globins of invertebrates. Cygb, present in many different tissues, shares common ancestry with myoglobin (Mb) and can be traced to early vertebrate evolution. Ngb and Cygb display the classical three-on-three α-helical globin fold and are endowed with a hexa-coordinate heme Fe atom, in both their ferrous and ferric forms, having the heme distal HisE7 residue as the endogenous sixth ligand. Reversible intramolecular hexa- to penta-coordination of the heme Fe atom modulates Ngb and Cygb ligand-binding properties. In Ngb and Cygb, ligand migration to/from the heme distal site may be assisted by protein/matrix tunnel cavity systems. The physiological roles of Ngb and Cygb are poorly understood. Ngb may protect neuronal cells from hypoxic-ischemic insults, may act as oxidative stress-responsive sensor protein, and may sustain NO/O(2) scavenging and/or reactive oxygen species (ROS) detoxification. Cygb, located in the cytoplasm of fibroblasts, chondroblasts, osteoblasts, and hepatic stellate cells, has been hypothesized to be involved in collagen synthesis. In neurons, Cygb, located in both cytoplasm and nucleus, may provide O(2) for enzymatic reactions, and may be involved in a ROS (NO)-signaling pathway(s). Here, we review current knowledge on Ngb and Cygb in terms of their structure, function, and evolutionary links to the well-known human HbA and Mb.

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.

Crystallographic Analysis of Synechocystis Cyanoglobin Reveals the Structural Changes Accompanying Ligand Binding in a Hexacoordinate Hemoglobin

The crystal structures of cyanide and azide-bound forms of the truncated hemoglobin from Synechocystis are presented at 1.8 angstroms resolution. A comparison with the structure of the endogenously liganded protein reveals a conformational shift unprecedented in hemoglobins, and provides the first picture of a hexacoordinate hemoglobin in both the bis-histidyl and the exogenously coordinated states. The structural changes between the different conformations are confined to two regions of the protein; the B helix, and the E helix, including the EF loop. A molecular "hinge" controlling movement of the E helix is observed in the EF loop, which is composed of three principal structural elements: Arg64, the heme-d-propionate, and a three-residue extension of the F helix. Additional features of the structural transition between the two protein conformations are discussed as they relate to the complex ligand-binding behavior observed in hexacoordinate hemoglobins, and the potential physiological function of this class of proteins.

Mapping protein matrix cavities in human cytoglobin through Xe atom binding

Cytoglobin is the fourth recognized globin type, almost ubiquitously distributed in human tissues; its function is still poorly understood. Cytoglobin displays a core region of about 150 residues, structurally related to hemoglobin and myoglobin, and two extra segments, about 20 residues each, at the N- and C-termini. The core region hosts a large apolar cavity, held to provide a ligand diffusion pathway to/from the heme, and/or ligand temporary docking sites. Here we report the crystal structure (2.4A resolution, R-factor 19.1%) of a human cytoglobin mutant bearing the CysB2(38) --> Ser and CysE9(83) --> Ser substitutions (CYGB*), treated under pressurized xenon. Three Xe atoms bind to the heme distal site region of CYGB* mapping the protein matrix apolar cavity. Despite the conserved globin fold, the cavity found in CYGB* is structured differently from those recognized to play a functional role in myoglobin, neuroglobin, truncated hemoglobins, and Cerebratulus lacteus mini-hemoglobin.

Heme-Ligand Tunneling in Group I Truncated Hemoglobins

Truncated hemoglobins (trHbs) are small hemoproteins forming a separate cluster within the hemoglobin superfamily; their functional roles in bacteria, plants, and unicellular eukaryotes are marginally understood. Crystallographic investigations have shown that the trHb fold (a two-on-two alpha-helical sandwich related to the globin fold) hosts a protein matrix tunnel system offering a potential path for ligand diffusion to the heme distal site. The tunnel topology is conserved in group I trHbs, although with modulation of its size/structure. Here, we present a crystallographic investigation on trHbs from Mycobacterium tuberculosis, Chlamydomonas eugametos, and Paramecium caudatum, showing that treatment of trHb crystals under xenon pressure leads to binding of xenon atoms at specific (conserved) sites along the protein matrix tunnel. The crystallographic results are in keeping with data from molecular dynamics simulations, where a dioxygen molecule is left free to diffuse within the protein matrix. Modulation of xenon binding over four main sites is related to the structural properties of the tunnel system in the three trHbs and may be connected to their functional roles. In a parallel crystallographic investigation on M. tuberculosis trHbN, we show that butyl isocyanide also binds within the apolar tunnel, in excellent agreement with concepts derived from the xenon binding experiments. These results, together with recent data on atypical CO rebinding kinetics to group I trHbs, underline the potential role of the tunnel system in supporting diffusion, but also accumulation in multiple copies, of low polarity ligands/molecules within group I trHbs.

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.