Viscosity-dependent Relaxation Significantly Modulates the Kinetics of CO Recombination in the Truncated Hemoglobin TrHbN from Mycobacterium tuberculosis

Haemoglobins of Mycobacteria: structural features and biological functions

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

Ligand 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.

Role of PheE15 gate in ligand entry and nitric oxide detoxification function of mycobacterium tuberculosis truncated hemoglobin N

The truncated hemoglobin N, HbN, of Mycobacterium tuberculosis is endowed with a potent nitric oxide dioxygenase (NOD) activity that allows it to relieve nitrosative stress and enhance in vivo survival of its host. Despite its small size, the protein matrix of HbN hosts a two-branched tunnel, consisting of orthogonal short and long channels, that connects the heme active site to the protein surface. A novel dual-path mechanism has been suggested to drive migration of O(2) and NO to the distal heme cavity. While oxygen migrates mainly by the short path, a ligand-induced conformational change regulates opening of the long tunnel branch for NO, via a phenylalanine (PheE15) residue that acts as a gate. Site-directed mutagenesis and molecular simulations have been used to examine the gating role played by PheE15 in modulating the NOD function of HbN. Mutants carrying replacement of PheE15 with alanine, isoleucine, tyrosine and tryptophan have similar O(2)/CO association kinetics, but display significant reduction in their NOD function. Molecular simulations substantiated that mutation at the PheE15 gate confers significant changes in the long tunnel, and therefore may affect the migration of ligands. These results support the pivotal role of PheE15 gate in modulating the diffusion of NO via the long tunnel branch in the oxygenated protein, and hence the NOD function of HbN.

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).

The pH dependence of heme pocket hydration and ligand rebinding kinetics in photodissociated carbonmonoxymyoglobin

We monitored the occupancy of a functionally important non-coordinated water molecule in the distal heme pocket of sperm whale myoglobin over the pH range 4.3-9.4. Water occupancy was assessed by using time-resolved spectroscopy to detect the perturbation of the heme visible band absorption spectrum caused by water entry after CO photodissociation ( Goldbeck, R. A., Bhaskaran, S., Ortega, C., Mendoza, J. L., Olson, J. S., Soman, J., Kliger, D. S., and Esquerra, R. M. (2006) Proc. Natl. Acad. Sci. U. S. A. 103, 1254-1259 ). We found that the water occupancy observed during the time interval between ligand photolysis and diffusive recombination decreased by nearly 20% as the pH was lowered below 6. This decrease accounted for most of the concomitant increase in the observed CO bimolecular recombination rate constant, as the lower water occupancy presented a smaller kinetic barrier to CO entry into the pocket at lower pH. These results were consistent with a model in which the distal histidine, which stabilizes the water molecule within the distal pocket by accepting a hydrogen bond, tends to swing out of the pocket upon protonation and destabilize the water occupancy at low pH. Extrapolation of this model to lower pH suggests that the additional increase in ligand association rate constant observed previously in stopped-flow studies at pH 3 may also be due in part to reduced distal water occupancy concomitant with further His64 protonation and coupled protein conformational change.

Conformational dependence of hemoglobin reactivity under high viscosity conditions: the role of solvent slaved dynamics

The concept of protein dynamic states is introduced. This concept is based on (i) protein dynamics being organized hierarchically with respect to solvent slaving and (ii) which tier of dynamics is operative over the time window of a given measurement. The protein dynamic state concept is used to analyze the kinetic phases derived from the recombination of carbon monoxide to sol-gel-encapsulated human adult hemoglobin (HbA) and select recombinant mutants. The temperature-dependent measurements are made under very high viscosity conditions obtained by bathing the samples in an excess of glycerol. The results are consistent with a given tier of solvent slaved dynamics becoming operative at a time delay (with respect to the onset of the measurement) that is primarily solvent- and temperature-dependent. However, the functional consequences of the dynamics are protein- and conformation-specific. The kinetic traces from both equilibrium populations and trapped allosteric intermediates show a consistent progression that exposes the role of both conformation and hydration in the control of reactivity. Iron-zinc symmetric hybrid forms of HbA are used to show the dramatic difference between the kinetic patterns for T state alpha and beta subunits. The overall results support a model for allostery in HbA in which the ligand-binding-induced transition from the deoxy T state to the high -affinity R state proceeds through a progression of T state intermediates.

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.

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.

Proximal influences in two-on-two globins: effect of the Ala69Ser replacement on Synechocystis sp. PCC 6803 hemoglobin

The cyanobacterium Synechocystis sp. PCC 6803 (S6803) expresses a two-on-two globin in which His46 (distal side) and His70 (proximal) function as heme iron axial ligands. His46 can be displaced by O2, CO, and CN-, among others, whereas His70 is not labile under native conditions. The residue preceding the proximal histidine has been implicated in controlling globin axial ligand reactivity; the details of the mechanism, however, are not well understood, and little information exists for bis-histidyl hexacoordinate proteins. In many vertebrate hemoglobins and in the Synechocystis protein, the position is occupied by an alanine, whereas, in myoglobins, it is a serine involved in an intricate hydrogen-bond network. We examined the role of Ala69 in S6803 hemoglobin through the effects of an Ala --> Ser replacement. The substitution resulted in minor structural perturbations, but the response of the holoprotein to temperature-, urea-, and acid-induced denaturation was measurably affected. Enhanced three-state behavior was manifested in the decoupling of heme binding and secondary-structure formation. Urea-gradient gel experiments revealed that the stability of the apoprotein was unchanged by the replacement and that a slight alteration of the folding kinetics occurred in the holoproteins. Cyanide-binding experiments were performed to assess trans effects. The apparent rate constant for association decreased 2-fold upon Ala69Ser replacement. This deceleration was attributed to a change in the lifetime of a state containing a decoordinated His46. The results demonstrated that, as in vertebrate globins and leghemoglobin, proximal influences operate to determine fundamental dynamic and thermodynamic properties of the protein.

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.

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.