Ligand migration in nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana

Kinetic and dynamical properties of truncated hemoglobins of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125

Due to the low temperature, the Antarctic marine environment is challenging for protein functioning. Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologs, resulting in enhanced reaction rates at low temperatures. The Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 (PhTAC125) genome is one of the few examples of coexistence of multiple hemoglobin genes encoding, among others, two constitutively transcribed 2/2 hemoglobins (2/2Hbs), also named truncated Hbs (TrHbs), belonging to the Group II (or O), annotated as PSHAa0030 and PSHAa2217. In this work, we describe the ligand binding kinetics and their interrelationship with the dynamical properties of globin Ph-2/2HbO-2217 by combining experimental and computational approaches and implementing a new computational method to retrieve information from molecular dynamic trajectories. We show that our approach allows us to identify docking sites within the protein matrix that are potentially able to transiently accommodate ligands and migration pathways connecting them. Consistently with ligand rebinding studies, our modeling suggests that the distal heme pocket is connected to the solvent through a low energy barrier, while inner cavities play only a minor role in modulating rebinding kinetics.

Probing the Role of Murine Neuroglobin CDloop-D-Helix Unit in CO Ligand Binding and Structural Dynamics

We produced a neuroglobin variant, namely, Ngb CDless, with the excised CDloop- and D-helix, directly joining the C- and E-helices. The CDless variant retained bis-His hexacoordination, and we investigated the role of the CDloop–D-helix unit in controlling the CO binding and structural dynamics by an integrative approach based on X-ray crystallography, rapid mixing, laser flash photolysis, resonance Raman spectroscopy, and molecular dynamics simulations. Rapid mixing and laser flash photolysis showed that ligand affinity was unchanged with respect to the wild-type protein, albeit with increased on and off constants for rate-limiting heme iron hexacoordination by the distal His64. Accordingly, resonance Raman spectroscopy highlighted a more open distal pocket in the CO complex that, in agreement with MD simulations, likely involves His64 swinging inward and outward of the distal heme pocket. Ngb CDless displays a more rigid overall structure with respect to the wild type, abolishing the structural dynamics of the CDloop–D-helix hypothesized to mediate its signaling role, and it retains ligand binding control by distal His64. In conclusion, this mutant may represent a tool to investigate the involvement of CDloop–D-helix in neuroprotective signaling in a cellular or animal model.

Rational Design of a User-Friendly Aptamer/Peptide-Based Device for the Detection of Staphylococcus aureus

The urgent need to develop a detection system for Staphylococcus aureus, one of the most common causes of infection, is prompting research towards novel approaches and devices, with a particular focus on point-of-care analysis. Biosensors are promising systems to achieve this aim. We coupled the selectivity and affinity of aptamers, short nucleic acids sequences able to recognize specific epitopes on bacterial surface, immobilized at high density on a nanostructured zirconium dioxide surface, with the rational design of specifically interacting fluorescent peptides to assemble an easy-to-use detection device. We show that the displacement of fluorescent peptides upon the competitive binding of S. aureus to immobilized aptamers can be detected and quantified through fluorescence loss. This approach could be also applied to the detection of other bacterial species once aptamers interacting with specific antigens will be identified, allowing the development of a platform for easy detection of a pathogen without requiring access to a healthcare environment.

Interaction of human hemoglobin and semi-hemoglobins with the Staphylococcus aureus hemophore IsdB: a kinetic and mechanistic insight

Among multidrug-resistant bacteria, methicillin-resistant Staphylococcus aureus is emerging as one of the most threatening pathogens. S. aureus exploits different mechanisms for its iron supply, but the preferred one is acquisition of organic iron through the expression of hemoglobin (Hb) receptors. One of these, IsdB, belonging to the Isd (Iron-Regulated Surface Determinant) system, was shown to be essential for bacterial growth and virulence. Therefore, interaction of IsdB with Hb represents a promising target for the rational design of a new class of antibacterial molecules. However, despite recent investigations, many structural and mechanistic details of complex formation and heme extraction process are still elusive. By combining site-directed mutagenesis, absorption spectroscopy, surface plasmon resonance and molecular dynamics simulations, we tackled most of the so far unanswered questions: (i) the exact complex stoichiometry, (ii) the microscopic kinetic rates of complex formation, (iii) the IsdB selectivity for binding to, and extracting heme from, α and β subunits of Hb, iv) the role of specific amino acid residues and structural regions in driving complex formation and heme transfer, and (v) the structural/dynamic effect played by the hemophore on Hb.

Cold-Adaptation Signatures in the Ligand Rebinding Kinetics to the Truncated Hemoglobin of the Antarctic Bacterium Pseudoalteromonas haloplanktis TAC125

Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologues, resulting in enhanced reaction rates at low temperatures. In this context, protein-bound water molecules were suggested to play a major role, and their weaker interactions at protein active sites have been associated with cold adaptation. In this work, we tested this hypothesis on truncated hemoglobins (a family of microbial heme-proteins of yet-unclear function) applying molecular dynamics simulations and ligand-rebinding kinetics on a protein from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 in comparison with its thermophilic Thermobifida fusca homologue. The CO rebinding kinetics of the former highlight several geminate phases, with an unusually long-lived geminate intermediate. An articulated tunnel with at least two distinct docking sites was identified by analysis of molecular dynamics simulations and was suggested to be at the origin of the unusual geminate rebinding phase. Water molecules are present in the distal pocket, but their stabilization by TrpG8, TyrB10, and HisCD1 is much weaker than in thermophilic Thermobifida fusca truncated hemoglobin, resulting in a faster geminate rebinding. Our results support the hypothesis that weaker water-molecule interactions at the reaction site are associated with cold adaptation.

The Greenland shark Somniosus microcephalus-Hemoglobins and ligand-binding properties

A large amount of data is currently available on the adaptive mechanisms of polar bony fish hemoglobins, but structural information on those of cartilaginous species is scarce. This study presents the first characterisation of the hemoglobin system of one of the longest-living vertebrate species (392 ± 120 years), the Arctic shark Somniosus microcephalus. Three major hemoglobins are found in its red blood cells and are made of two copies of the same α globin combined with two copies of three very similar β subunits. The three hemoglobins show very similar oxygenation and carbonylation properties, which are unaffected by urea, a very important compound in marine elasmobranch physiology. They display identical electronic absorption and resonance Raman spectra, indicating that their heme-pocket structures are identical or highly similar. The quaternary transition equilibrium between the relaxed (R) and the tense (T) states is more dependent on physiological allosteric effectors than in human hemoglobin, as also demonstrated in polar teleost hemoglobins. Similar to other cartilaginous fishes, we found no evidence for functional differentiation among the three isoforms. The very similar ligand-binding properties suggest that regulatory control of O₂ transport may be at the cellular level and that it may involve changes in the cellular concentrations of allosteric effectors and/or variations of other systemic factors. The hemoglobins of this polar shark have evolved adaptive decreases in O₂ affinity in comparison to temperate sharks.

Characterization of the Heme Pocket Structure and Ligand Binding Kinetics of Non-symbiotic Hemoglobins from the Model Legume Lotus japonicus

Plant hemoglobins (Hbs) are found in nodules of legumes and actinorhizal plants but also in non-symbiotic organs of monocots and dicots. Non-symbiotic Hbs (nsHbs) have been classified into two phylogenetic groups. Class 1 nsHbs show an extremely high O2 affinity and are induced by hypoxia and nitric oxide (NO), whereas class 2 nsHbs have moderate O2 affinity and are induced by cold and cytokinins. The functions of nsHbs are still unclear, but some of them rely on the capacity of hemes to bind diatomic ligands and catalyze the NO dioxygenase (NOD) reaction (oxyferrous Hb + NO → ferric Hb + nitrate). Moreover, NO may nitrosylate Cys residues of proteins. It is therefore important to determine the ligand binding properties of the hemes and the role of Cys residues. Here, we have addressed these issues with the two class 1 nsHbs (LjGlb1-1 and LjGlb1-2) and the single class 2 nsHb (LjGlb2) of Lotus japonicus, which is a model legume used to facilitate the transfer of genetic and biochemical information into crops. We have employed carbon monoxide (CO) as a model ligand and resonance Raman, laser flash photolysis, and stopped-flow spectroscopies to unveil major differences in the heme environments and ligand binding kinetics of the three proteins, which suggest non-redundant functions. In the deoxyferrous state, LjGlb1-1 is partially hexacoordinate, whereas LjGlb1-2 shows complete hexacoordination (behaving like class 2 nsHbs) and LjGlb2 is mostly pentacoordinate (unlike other class 2 nsHbs). LjGlb1-1 binds CO very strongly by stabilizing it through hydrogen bonding, but LjGlb1-2 and LjGlb2 show lower CO stabilization. The changes in CO stabilization would explain the different affinities of the three proteins for gaseous ligands. These affinities are determined by the dissociation rates and follow the order LjGlb1-1 > LjGlb1-2 > LjGlb2. Mutations LjGlb1-1 C78S and LjGlb1-2 C79S caused important alterations in protein dynamics and stability, indicating a structural role of those Cys residues, whereas mutation LjGlb1-1 C8S had a smaller effect. The three proteins and their mutant derivatives exhibited similarly high rates of NO consumption, which were due to NOD activity of the hemes and not to nitrosylation of Cys residues.

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.

Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure-function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1)→Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

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

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

CO rebinding kinetics and molecular dynamics simulations highlight dynamic regulation of internal cavities in human cytoglobin

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

Heterogeneous kinetics of the carbon monoxide association and dissociation reaction to nitrophorin 4 and 7 coincide with structural heterogeneity of the gate-loop

NO is an important signaling molecule in human tissue. However, the mechanisms by which this molecule is controlled and directed are currently little understood. Nitrophorins (NPs) comprise a group of ferriheme proteins originating from blood-sucking insects that are tailored to protect and deliver NO via coordination to and release from the heme iron. Therefore, the kinetics of the association and dissociation reactions were studied in this work using the ferroheme-CO complexes of NP4, NP4(D30N), and NP7 as isoelectronic models for the ferriheme-NO complexes. The kinetic measurements performed by nanosecond laser-flash-photolysis and stopped-flow are accompanied by resonance Raman and FT-IR spectroscopy to characterize the carbonyl species. Careful analysis of the CO rebinding kinetics reveals that in NP4 and, to a larger extent, NP7 internal gas binding cavities are located, which temporarily trap photodissociated ligands. Moreover, changes in the free energy barriers throughout the rebinding and release pathway upon increase of the pH are surprisingly small in case of NP4. Also in case of NP4, a heterogeneous kinetic trace is obtained at pH 7.5, which corresponds to the presence of two carbonyl species in the heme cavity that are seen in vibrational spectroscopy and that are due to the change of the distal heme pocket polarity. Quantification of the two species from FT-IR spectra allowed the fitting of the kinetic traces as two processes, corresponding to the previously reported open and closed conformation of the A-B and G-H loops. With the use of the A-B loop mutant NP4(D30N), it was confirmed that the kinetic heterogeneity is controlled by pH through the disruption of the H-bond between the Asp30 side chain and the Leu130 backbone carbonyl. Overall, this first study on the slow phase of the dynamics of diatomic gas molecule interaction with NPs comprises an important experimental contribution for the understanding of the dynamics involved in the binding/release processes of NO/CO in NPs.

Ligation tunes protein reactivity in an ancient haemoglobin: kinetic evidence for an allosteric mechanism in Methanosarcina acetivorans protoglobin

Protoglobin from Methanosarcina acetivorans (MaPgb) is a dimeric globin with peculiar structural properties such as a completely buried haem and two orthogonal tunnels connecting the distal cavity to the solvent. CO binding to and dissociation from MaPgb occur through a biphasic kinetics. We show that the heterogenous kinetics arises from binding to (and dissociation from) two tertiary conformations in ligation-dependent equilibrium. Ligation favours the species with high binding rate (and low dissociation rate). The equilibrium is shifted towards the species with low binding (and high dissociation) rates for the unliganded molecules. A quantitative model is proposed to describe the observed carbonylation kinetics.

Non-symbiotic haemoglobins-What's happening beyond nitric oxide scavenging?

BACKGROUND AND AIMS

Non-symbiotic haemoglobins have been an active research topic for over 30 years, during which time a considerable portfolio of knowledge has accumulated relative to their chemical and molecular properties, and their presence and mode of induction in plants. While progress has been made towards understanding their physiological role, there remain a number of unanswered questions with respect to their biological function. This review attempts to update recent progress in this area and to introduce a hypothesis as to how non-symbiotic haemoglobins might participate in regulating hormone signal transduction.

PRINCIPAL RESULTS

Advances have been made towards understanding the structural nuances that explain some of the differences in ligand association characteristics of class 1 and class 2 non-symbiotic haemoglobins. Non-symbiotic haemoglobins have been found to function in seed development and germination, flowering, root development and differentiation, abiotic stress responses, pathogen invasion and symbiotic bacterial associations. Microarray analyses under various stress conditions yield uneven results relative to non-symbiotic haemoglobin expression. Increasing evidence of the role of nitric oxide (NO) in hormone responses and the known involvement of non-symbiotic haemoglobins in scavenging NO provide opportunities for fruitful research, particularly at the cellular level.

CONCLUSIONS

Circumstantial evidence suggests that non-symbiotic haemoglobins may have a critical function in the signal transduction pathways of auxin, ethylene, jasmonic acid, salicylic acid, cytokinin and abscisic acid. There is a strong need for research on haemoglobin gene expression at the cellular level relative to hormone signal transduction.

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

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

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

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

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.

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.

Ligand migration and binding in nonsymbiotic hemoglobins of Arabidopsis thaliana

We have studied carbon monoxide (CO) migration and binding in the nonsymbiotic hemoglobins AHb1 and AHb2 of Arabidopsis thaliana using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures. Both proteins have similar amino acid sequences but display pronounced differences in ligand binding properties, at both physiological and cryogenic temperatures. Near neutral pH, the distal HisE7 side chain is close to the heme-bound ligand in the majority of AHb1-CO molecules, as indicated by a low CO stretching frequency at 1921 cm(-1). In this fraction, two CO docking sites can be populated, the primary site B and the secondary site C. When the pH is lowered, a high-frequency stretching band at approximately 1964 cm(-1) grows at the expense of the low-frequency band, indicating that HisE7 protonates and, concomitantly, moves away from the bound ligand. Geminate rebinding barriers are markedly different for the two conformations, and docking site C is not accessible in the low-pH conformation. Rebinding of NO ligands was observed only from site B of AHb1, regardless of conformation. In AHb2, the HisE7 side chain is removed from the bound ligand; rebinding barriers are low, and CO molecules can populate only primary docking site B. These results are interpreted in terms of differences in the active site structures and physiological functions.

Ligand dynamics in heme proteins observed by Fourier transform infrared-temperature derivative spectroscopy

Fourier transform infrared (FTIR) spectroscopy is a powerful tool for the investigation of protein-ligand interactions in heme proteins. Nitric oxide and carbon monoxide are attractive physiologically relevant ligands because their bond stretching vibrations give rise to strong mid-infrared absorption bands that can be measured with exquisite sensitivity and precision using photolysis difference spectroscopy at cryogenic temperatures. These stretching bands are fine-tuned by electrostatic interactions with the environment and, therefore, ligands can be utilized as local probes of structure and dynamics. Bound to the heme iron, the ligand stretching bands are susceptible to changes in the iron-ligand bond and the electric field at the active site. Upon photolysis, the vibrational bands display changes due to ligand relocation to docking sites within the protein, rotational motions of the ligand in these sites and protein conformational changes. Photolysis difference spectra taken over a wide temperature range (3-300K) using specific temperature protocols for sample photodissociation can provide detailed insights into both protein and ligand dynamics. Moreover, temperature-derivative spectroscopy (TDS) has proven to be a particularly powerful technique to study protein-ligand interactions. The FTIR-TDS technique has been extensively applied to studies of carbon monoxide binding to heme proteins, whereas measurements with nitric oxide are still scarce. Here we describe infrared cryo-spectroscopy and present a variety of applications to the study of protein-ligand interactions in heme proteins. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Breathing motions of a respiratory protein revealed by molecular dynamics simulations

Internal cavities, which are central to the biological functions of myoglobin, are exploited by gaseous ligands (e.g., O(2), NO, CO, etc.) to migrate inside the protein matrix. At present, it is not clear whether the ligand makes its own way inside the protein or instead the internal cavities are an intrinsic feature of myoglobin. To address this issue, standard molecular dynamics simulations were performed on horse-heart met-myoglobin with no ligand migrating inside the protein matrix. To reveal intrinsic internal pathways, the use of a statistical approach was applied to the cavity calculation, with special emphasis on the major pathway from the distal pocket to Xe1. Our study points out the remarkable dynamical behavior of Xe4, whose "breathing motions" may facilitate migration of ligands through the distal region. Additionally, our results highlight a two-way path for a ligand to diffuse through the proximal region, possibly allowing an alternative route in case Xe1 is occupied. Finally, our approach has led us to the identification of key residues, such as leucines, that may work as switches between cavities.

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.

Nonsymbiotic hemoglobins and stress tolerance in plants

Hemoglobins (Hbs) are heme containing proteins found in most organisms including animals, bacteria, and plants. Their structure, size, and function are quite diverse among the different organisms. There are three different types of hemoglobins in plants: symbiotic (sHb), nonsymbiotic (nsHb), and truncated hemoglobins (trHb). The nonsymbiotic hemoglobins are divided into: class 1 hemoglobins (nsHb-1s), which have a very high affinity for oxygen: and class 2 hemoglobins (nsHb-2s), which have lower affinity for oxygen, are similar to the sHbs. nsHb-1s are expressed under hypoxia, osmotic stress, nutrient deprivation, cold stress, rhizobial infection, nitric oxide exposure, and fungal infection. Tolerance to stress is very important for the survival of the plant. Hemoglobins are one of many different strategies that plants have evolved to overcome stress conditions and survive. Hbs also react with NO produced under different stress conditions. Class 1 nsHbs are involved in a metabolic pathway involving NO. Those hemoglobins provide an alternative type of respiration to mitochondrial electron transport under limiting oxygen concentrations. Class 1 nsHbs in hypoxic plants act as part of a soluble, terminal, NO dioxygenase system, yielding nitrate from the reaction of oxyHb with NO. The overall reaction sequence, referred to as the nsHb/NO cycle, consumes NADH and maintains ATP levels via an as yet unknown mechanism. Class 2 nsHbs seem to scavenge NO in a similar fashion as class 1 Hbs and are involved in reducing flowering time in Arabidopsis. nsHbs also show peroxidase-like activity and NO metabolism and possibly protect against nitrosative stress in plant-pathogen interaction and in symbiotic interactions. nsHbs may be involved in other stress conditions such as osmotic, nutrient and cold stress together with NO and the function of nsHbs can be in NO metabolism and signal transduction. However, other possible functions cannot be precluded as Hbs have many different functions in other organisms.