Role of Pre-A Motif in Nitric Oxide Scavenging by Truncated Hemoglobin, HbN, of Mycobacterium tuberculosis

Mycobacterium tuberculosis truncated hemoglobin, HbN, is endowed with a potent nitric-oxide dioxygenase activity and has been found to relieve nitrosative stress and enhance in vivo survival of a heterologous host, Salmonella enterica Typhimurium, within the macrophages. These findings implicate involvement of HbN in the defense of M. tuberculosis against nitrosative stress. The protein carries a tunnel system composed of a short and a long tunnel branch that has been proposed to facilitate diatomic ligand migration to the heme and an unusual Pre-A motif at the N terminus, which does not contribute significantly to the structural integrity of the protein, as it protrudes out of the compact globin fold. Strikingly, deletion of Pre-A region from the M. tuberculosis HbN drastically reduces its ability to scavenge nitric oxide (NO), whereas its insertion at the N terminus of Pre-A lacking HbN of Mycobacterium smegmatis improved its nitric-oxide dioxygenase activity. Titration of the oxygenated adduct of HbN and its mutants with NO indicated that the stoichiometric oxidation of protein is severalfold slower when the Pre-A region is deleted in HbN. Molecular dynamics simulations show that the excision of Pre-A motif results in distinct changes in the protein dynamics, which cause the gate of the tunnel long branch to be trapped into a closed conformation, thus impeding migration of diatomic ligands toward the heme active site. The present study, thus, unequivocally demonstrates vital function of Pre-A region in NO scavenging and unravels its unique role by which HbN might attain its efficient NO-detoxification ability.Unlike many pathogens that are overtly harmful to their host, Mycobacterium tuberculosis can persist for years within humans in a clinically latent state. The success of tubercle bacillus to establish long term persistent infection within the human host lies in its ability to survive and resist hazardous environment of its intracellular niche. Early infection events involve entry and multiplication within the bacteriostatic environment of macrophages (1, 2), where an inducible nitric-oxide synthase generates copious amounts of nitric oxide (NO), which plays an important role in the host defense against microbial pathogens (3, 4) and restricts their growth and survival by inhibiting key enzymes such as the terminal respiratory oxidases (5) and iron-sulfur centers of key enzymes such as aconitase (6, 7). In addition, NO combines with superoxide produced by respiring cells to form the highly oxidizing agent peroxynitrite (8). An efficient NO scavenging system, therefore, is required by microbial pathogens to cope with NO poisoning during their intracellular regime.M. tuberculosis has evolved resistance mechanisms by which toxic effects of NO and nitrosative stress can be evaded. One of the unique defense mechanisms by which it can protect itself from reactive nitrogen species relies on the oxygenated form of a group I truncated hemoglobin, HbN, which very efficiently converts NO into harmless nitrate (9–11). Studies conducted in our laboratory have shown that HbN carries a potent O₂-dependent NO dioxygenase (NOD)³ activity and protects its host from nitrosative stress (11) and microbicidal activities of macrophages (12). It has also been demonstrated that disruption of glbN gene, encoding HbN, in Mycobacterium bovis bacillus Calmette-Guérin causes a dramatic reduction in the NO-consuming activity of stationary phase cells, resulting in marked NO-induced inhibition of aerobic respiration relative to the wild type cells (10). Overall, these studies strongly support the NO-scavenging and detoxification roles of HbN, which may be vital for in vivo survival and pathogenicity of M. tuberculosis.X-ray crystallographic and computational studies have highlighted the relationship between structure and NOD function of key residues in M. tuberculosis HbN. The protein hosts a two-branched protein matrix tunnel system (13, 14) and has evolved a novel dual-path mechanism to drive migration of O₂ and NO to the distal heme cavity where a ligand-induced conformational change regulates the opening of the tunnel branch via phenylalanine (Phe(E15)) residue that acts as a gate for ligand entry (15). Extended molecular dynamics (MD) simulations showed that the essential motions of the protein backbone and the positioning of Phe(E15) gate changes significantly after O₂ binding, which facilitates access of a diatomic ligand (NO) to the active site to carry out NO detoxification (16, 17). Finally, a novel pathway has been identified for the egression of nitrate anion (18).HbN homologs are less represented in living organisms as compared with its group II molecular relative HbO. A variety of physiological functions are displayed by HbN-type truncated hemoglobins (trHbs). For example, Chlamydomonas eugametos HbN has been implicated in photosynthesis (19, 20), whereas the probable function of Paramecium caudatum HbN is O₂ supply to mitochondria (21), and HbN of the cyanobacterium Nostoc commune may be part of a microaerobically functioning electron transfer system (22). At present, M. tuberculosis HbN is the only type I truncated hemoglobin for which a potent NOD activity has been ascribed. Recent genome data from various mycobacterial species indicate that HbN-type trHbs exhibit 70–80% sequence similarity (21), which a priori suggests a common globin fold and a similar biological function. Nevertheless, the HbN counterpart in Mycobacterium smegmatis displays severalfold lower NOD activity (23) despite having conserved features for the truncated globin fold and close homology with M. tuberculosis HbN. In fact, residues suggested to play a crucial role for ligand migration through the protein matrix are fully conserved (11). The critical question that emerges now is how M. tuberculosis HbN is able to achieve an efficient NO detoxification ability and what aspects of its structure and/or dynamics contribute to this property.This study aims to understand the molecular basis of the differences in the NO-metabolizing ability of HbN from M. tuberculosis and M. smegmatis. The results demonstrate a vital contribution of Pre-A helical region in regulation of NOD activity, unraveling a unique mechanism that contributes to modulate the highly efficient NO detoxification activity of HbN.     

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.     

Truncated Hemoglobin,HbN, Is Post-translationally Modified in Mycobacterium tuberculosis and Modulates Host-Pathogen Interactions during Intracellular Infection

Mycobacterium tuberculosis (Mtb) is a phenomenally successful human pathogen having evolved mechanisms that allow it to survive within the hazardous environment of macrophages and establish long term, persistent infection in the host against the control of cell-mediated immunity. One such mechanism is mediated by the truncated hemoglobin, HbN, of Mtb that displays a potent O₂-dependent nitric oxide dioxygenase activity and protects its host from the toxicity of macrophage-generated nitric oxide (NO). Here we demonstrate for the first time that HbN is post-translationally modified by glycosylation in Mtb and remains localized on the cell membrane and the cell wall. The glycan linkage in the HbN was identified as mannose. The elevated expression of HbN in Mtb and M. smegmatis facilitated their entry within the macrophages as compared with isogenic control cells, and mutation in the glycan linkage of HbN disrupted this effect. Additionally, HbN-expressing cells exhibited higher survival within the THP-1 and mouse peritoneal macrophages, simultaneously increasing the intracellular level of proinflammatory cytokines IL-6 and TNF-α and suppressing the expression of co-stimulatory surface markers CD80 and CD86. These results, thus, suggest the involvement of HbN in modulating the host-pathogen interactions and immune system of the host apart from protecting the bacilli from nitrosative stress inside the activated macrophages, consequently driving cells toward increased infectivity and intracellular survival.     

Oxygen binding and NO scavenging properties of truncated hemoglobin, HbN, of Mycobacterium smegmatis

Unraveling of microbial genome data has indicated that two distantly related truncated hemoglobins (trHbs), HbN and HbO, might occur in many species of slow-growing pathogenic mycobacteria. Involvement of HbN in bacterial defense against NO toxicity and nitrosative stress has been proposed. A gene, encoding a putative HbN homolog with conserved features of typical trHbs, has been identified within the genome sequence of fast-growing mycobacterium, Mycobacterium smegmatis. Sequence analysis of M. smegmatis HbN indicated that it is relatively smaller in size and lacks N-terminal pre-A region, carrying 12-residue polar sequence motif that is present in HbN of M. tuberculosis. HbN encoding gene of M. smegmatis was expressed in E. coli as a 12.8kD homodimeric heme protein that binds oxygen reversibly with high affinity (P50 approximately 0.081 mm Hg) and autooxidizes faster than M. tuberculosis HbN. The circular dichroism spectra indicate that HbN of M. smegmatis and M. tuberculosis are structurally similar. Interestingly, an hmp mutant of E. coli, unable to metabolize nitric oxide, exhibited very low NO uptake activity in the presence of M. smegmatis HbN as compared to HbN of M. tuberculosis. On the basis of cellular heme content, specific nitric oxide dioxygenase (NOD) activity of M. smegmatis HbN was nearly one-third of that from M. tuberculosis. Additionally, the hmp mutant of E. coli, carrying M. smegmatis HbN, exhibited nearly 10-fold lower cell survival under nitrosative stress and nitrite derived reactive nitrogen species as compared to the isogenic strain harboring HbN of M. tuberculosis. Taken together, these results suggest that NO metabolizing activity and protection provided by M. smegmatis HbN against toxicity of NO and reactive nitrogen is significantly lower than HbN of M. tuberculosis. The lower efficiency of M. smegmatis HbN for NO detoxification as compared to M. tuberculosis HbN might be related to different level of NO exposure and nitrosative stress faced by these mycobacteria during their cellular metabolism.     

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 ≈ 745 μM⁻¹ s⁻¹ at 23°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⁶²(E15), a residue proposed to act as a gate controlling ligand traffic inside the Long tunnel, and also of Ile¹¹⁹(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.     

Mechanistic Insight into the Enzymatic Reduction of Truncated Hemoglobin N of Mycobacterium tuberculosis: ROLE OF THE CD LOOP AND PRE-A MOTIF IN ELECTRON CYCLING*

Many pathogenic microorganisms have evolved hemoglobin-mediated nitric oxide (NO) detoxification mechanisms, where a globin domain in conjunction with a partner reductase catalyzes the conversion of toxic NO to innocuous nitrate. The truncated hemoglobin HbN of Mycobacterium tuberculosis displays a potent NO dioxygenase activity despite lacking a reductase domain. The mechanism by which HbN recycles itself during NO dioxygenation and the reductase that participates in this process are currently unknown. This study demonstrates that the NADH-ferredoxin/flavodoxin system is a fairly efficient partner for electron transfer to HbN with an observed reduction rate of 6.2 μm/min⁻¹, which is nearly 3- and 5-fold faster than reported for Vitreoscilla hemoglobin and myoglobin, respectively. Structural docking of the HbN with Escherichia coli NADH-flavodoxin reductase (FdR) together with site-directed mutagenesis revealed that the CD loop of the HbN forms contacts with the reductase, and that Gly⁴⁸ may have a vital role. The donor to acceptor electron coupling parameters calculated using the semiempirical pathway method amounts to an average of about 6.4 10⁻⁵ eV, which is lower than the value obtained for E. coli flavoHb (8.0 10⁻⁴ eV), but still supports the feasibility of an efficient electron transfer. The deletion of Pre-A abrogated the heme iron reduction by FdR in the HbN, thus signifying its involvement during intermolecular interactions of the HbN and FdR. The present study, thus, unravels a novel role of the CD loop and Pre-A motif in assisting the interactions of the HbN with the reductase and the electron cycling, which may be vital for its NO-scavenging function.     

Theoretical investigation on the diatomic ligand migration process and ligand binding properties in non‐O2‐binding H‐NOX domain

The Nostoc sp (Ns) H‐NOX (heme‐nitric oxide or OXygen‐binding) domain shares 35% sequence identity with soluble guanylate cyclase (sGC) and exhibits similar ligand binding property with the sGC. Previously, our molecular dynamic (MD) simulation work identified that there exists a Y‐shaped tunnel system hosted in the Ns H‐NOX interior, which servers for ligand migration. The tunnels were then confirmed by Winter et al. [PNAS 2011;108(43):E 881–889] recently using x‐ray crystallography with xenon pressured conditions. In this work, to further investigate how the protein matrix of Ns H‐NOX modulates the ligand migration process and how the distal residue composition affects the ligand binding prosperities, the free energy profiles for nitric oxide (NO), carbon monooxide (CO), and O₂ migration are explored using the steered MDs simulation and the ligand binding energies are calculated using QM/MM schemes. The potential of mean force profiles suggest that the longer branch of the tunnel would be the most favorable route for NO migration and a second NO trapping site other than the distal heme pocket along this route in the Ns H‐NOX was identified. On the contrary, CO and O₂ would prefer to diffuse via the shorter branch of the tunnel. The QM/MM (quantum mechanics/molecular mechanics) calculations suggest that the hydrophobic distal pocket of Ns H‐NOX would provide an approximately vacuum environment and the ligand discrimination would be determined by the intrinsic binding properties of the diatomic gas ligand to the heme group. Proteins 2013; 81:1363–1376. © 2013 Wiley Periodicals, Inc.     

Nitric oxide scavenging and detoxification by the Mycobacterium tuberculosis haemoglobin,HbN in Escherichia coli

Nitric oxide (NO), generated in large amounts within the macrophages, controls and restricts the growth of internalized human pathogen, Mycobacterium tuberculosis H37Rv. The molecular mechanism by which tubercle bacilli survive within macrophages is currently of intense interest. In this work, we have demonstrated that dimeric haemoglobin, HbN, from M. tuberculosis exhibits distinct nitric oxide dioxygenase (NOD) activity and protects growth and cellular respiration of heterologous hosts, Escherichia coli and Mycobacterium smegmatis, from the toxic effect of exogenous NO and the NO-releasing compounds. A flavohaemoglobin (HMP)-deficient mutant of E. coli, unable to metabolize NO, acquired an oxygen-dependent NO consumption activity in the presence of HbN. On the basis of cellular haem content, the specific NOD activity of HbN was nearly 35-fold higher than the single-domain Vitreoscilla haemoglobin (VHb) but was sevenfold lower than the two-domain flavohaemoglobin. HbN-dependent NO consumption was sustained with repeated addition of NO, demonstrating that HbN is catalytically reduced within E. coli. Aerobic growth and respiration of a flavohaemoglobin (HMP) mutant of E. coli was inhibited in the presence of exogenous NO but remained insensitive to NO inhibition when these cells produced HbN, VHb or flavohaemoglobin. M. smegmatis, carrying a native HbN very similar to M. tuberculosis HbN, exhibited a 7.5-fold increase in NO uptake when exposed to gaseous NO, suggesting NO-induced NOD activity in these cells. In addition, expression of plasmid-encoded HbN of M. tuberculosis in M. smegmatis resulted in 100-fold higher NO consumption activity than the isogenic control cells. These results provide strong experimental evidence in support of NO scavenging and detoxification function for the M. tuberculosis HbN. The catalytic NO scavenging by HbN may be highly advantageous for the survival of tubercle bacilli during infection and pathogenesis.     

Gates and Binding Pockets for Nitric Oxide with Cytochrome c′, According to Molecular Dynamics

Random‐acceleration molecular‐dynamics (RAMD) simulations with models of homodimeric 6‐ligated distal‐NO and 5‐ligated proximal‐NO cytochrome c′ complexes, in TIP3 H₂O, showed two distinct, non‐intercommunicating worlds. In the framework of a long cavity formed by four protein helices with heme at one extremity, NO was observed to follow different pathways with the two complexes to reach the solvent. With the 6‐ligated complex, NO was observed to progress by exploiting protein internal channels created by thermal fluctuations, and be temporarily trapped into binding pockets before reaching the preferred gate at the heme end of the cavity. In contrast, with the 5‐ligated complex, NO was observed to surface the solvent‐exposed helix 7, up to a gate at the other extremity of the protein, only occasionally finding an earlier, direct way out toward the solvent. That only bulk NO gets involved in forming the 5‐ligated proximal‐NO complex is in agreement with previous experimental observations, while the occurrence of binding pockets suggests that also reservoir NO might play a role with the distal‐NO complex.     

Bacterial denitrifying nitric oxide reductases and aerobic respiratory terminal oxidases use similar delivery pathways for their molecular substrates

The superfamily of heme?copper oxidoreductases (HCOs) include both NO and O₂ reductases. Nitric oxide reductases (NORs) are bacterial membrane enzymes that catalyze an intermediate step of denitrification by reducing nitric oxide (NO) to nitrous oxide (N₂O). They are structurally similar to heme?copper oxygen reductases (HCOs), which reduce O₂ to water. The experimentally observed apparent bimolecular rate constant of NO delivery to the deeply buried catalytic site of NORs was previously reported to approach the diffusion-controlled limit (10⁸–10⁹?M^?1?s^?1). Using the crystal structure of cytochrome-c dependent NOR (cNOR) from Pseudomonas aeruginosa, we employed several protocols of molecular dynamics (MD) simulation, which include flooding simulations of NO molecules, implicit ligand sampling and umbrella sampling simulations, to elucidate how NO in solution accesses the catalytic site of this cNOR. The results show that NO partitions into the membrane, enters the enzyme from the lipid bilayer and diffuses to the catalytic site via a hydrophobic tunnel that is resolved in the crystal structures. This is similar to what has been found for O₂ diffusion through the closely related O₂ reductases. The apparent second order rate constant approximated using the simulation data is ~5?×?10⁸?M^?1?s^?1, which is optimized by the dynamics of the amino acid side chains lining in the tunnel. It is concluded that both NO and O₂ reductases utilize well defined hydrophobic tunnels to assure that substrate diffusion to the buried catalytic sites is not rate limiting under physiological conditions.     

Catalase evolved to concentrate H2O2 at its active site

Catalase is a homo-tetrameric enzyme that has its heme active site deeply buried inside the protein. Its only substrate, hydrogen peroxide (H₂O₂), reaches the heme through a 45 Å-long channel. Large-subunit catalases, but not small-subunit catalases, have a loop (gate loop) that interrupts the major channel. Two accesses lead to a gate that opens the final section of the channel to the heme; gates from the R-related subunits are interconnected. Using molecular dynamic simulations of the Neurospora crassa catalase-1 tetramer in a box of water (48,600 molecules) or 6 M H₂O₂, it is shown that the number of H₂O₂ molecules augments at the surface of the protein and in the accesses to the gate and the final section of the channel. Increase in H₂O₂ is due to the prevalence and distribution of amino acids that have an increased residency for H₂O₂ (mainly histidine, proline and charged residues), which are localized at the protein surface and the accesses to the gate. In the section of the channel from the heme to the gate, turnover rate of water molecules was faster than for H₂O₂ and increased residence sites for water and H₂O₂ were determined. In the presence of H₂O₂, the exclusion of water molecules from a specific site suggests a mechanism that could contend with the competing activity of water, allowing for catalase high kinetic efficiency.     

Nitrosylation Mechanisms of Mycobacterium tuberculosis and Campylobacter jejuni Truncated Hemoglobins N,O, and P

Truncated hemoglobins (trHbs) are widely distributed in bacteria and plants and have been found in some unicellular eukaryotes. Phylogenetic analysis based on protein sequences shows that trHbs branch into three groups, designated N (or I), O (or II), and P (or III). Most trHbs are involved in the O₂/NO chemistry and/or oxidation/reduction function, permitting the survival of the microorganism in the host. Here, a detailed comparative analysis of kinetics and/or thermodynamics of (i) ferrous Mycobacterium tubertulosis trHbs N and O (Mt-trHbN and Mt-trHbO, respectively), and Campylobacter jejuni trHb (Cj-trHbP) nitrosylation, (ii) nitrite-mediated nitrosylation of ferrous Mt-trHbN, Mt-trHbO, and Cj-trHbP, and (iii) NO-based reductive nitrosylation of ferric Mt-trHbN, Mt-trHbO, and Cj-trHbP is reported. Ferrous and ferric Mt-trHbN and Cj-trHbP display a very high reactivity towards NO; however, the conversion of nitrite to NO is facilitated primarily by ferrous Mt-trHbN. Values of kinetic and/or thermodynamic parameters reflect specific trHb structural features, such as the ligand diffusion pathways to/from the heme, the heme distal pocket structure and polarity, and the ligand stabilization mechanisms. In particular, the high reactivity of Mt-trHbN and Cj-trHbP reflects the great ligand accessibility to the heme center by two protein matrix tunnels and the E7-path, respectively, and the penta-coordination of the heme-Fe atom. In contrast, the heme-Fe atom of Mt-trHbO the ligand accessibility to the heme center of Mt-trHbO needs large conformational readjustments, thus limiting the heme-based reactivity. These results agree with different roles of Mt-trHbN, Mt-trHbO, and Cj-trHbP in vivo.     

NO binding induced conformational changes in a truncated hemoglobin from Mycobacterium tuberculosis

The resonance Raman spectra of the NO-bound ferric derivatives of wild-type HbN and the B10 Tyr --> Phe mutant of HbN, a hemoglobin from Mycobacterium tuberculosis, were examined with both Soret and UV excitation. The Fe-N-O stretching and bending modes of the NO derivative of the wild-type protein were tentatively assigned at 591 and 579 cm(-1), respectively. Upon B10 mutation, the Fe-NO stretching mode was slightly enhanced and the bending mode diminished in amplitude. In addition, the N-O stretching mode shifted from 1914 to 1908 cm(-1). These data suggest that the B10 Tyr forms an H-bond(s) with the heme-bound NO and causes it to bend in the wild-type protein. To further investigate the interaction between the B10 Tyr and the heme-bound NO, we examined the UV Raman spectrum of the B10 Tyr by subtracting the B10 mutant spectrum from the wild-type spectrum. It was found that, upon NO binding to the ferric protein, the Y(8a) mode of the B10 Tyr shifted from 1616 to 1622 cm(-1), confirming a direct interaction between the B10 Tyr and the heme-bound NO. Furthermore, the Y(8a) mode of the other two Tyr residues at positions 16 and 72 that are remote from the heme was also affected by NO binding, suggesting that NO binding to the distal site of the heme triggers a large-scale conformational change that propagates through the pre-F helix loop to the E and B helices. This large-scale conformational change triggered by NO binding may play an important role in regulating the ligand binding properties and/or the chemical reactivity of HbN.     

How does a valine residue that modulates heme-NO binding kinetics in inducible NO synthase regulate enzyme catalysis?

Nitric oxide (NO) release from nitric oxide synthases (NOSs) depends on the dissociation of a ferric heme-NO product complex (Fe^IIINO) that forms immediately after NO is made in the heme pocket. The NOS-like enzyme of Bacillus subtilis (bsNOS) has 10-20 fold slower Fe^IIINO dissociation rate (k_d) and NO association rate (k_on) compared to mammalian NOS counterparts. We previously showed that an Ile for Val substitution at the opening of the heme pocket in bsNOS contributes to these differences. The complementary mutation in mouse inducible NOS oxygenase domain (Val346Ile) decreased the NO k_on and k_d by 8 and 3-fold, respectively, compared to wild-type iNOSoxy, and also slowed the reductive processing of the heme-O₂ catalytic intermediate. To investigate how these changes affect steady-state catalytic behaviors, we generated and characterized the V346I mutant of full-length inducible NOS (iNOS). The mutant exhibited a 4-5 fold lower NO synthesis activity, an apparent uncoupled NADPH consumption, and formation of a heme-NO complex during catalysis that was no longer sensitive to solution NO scavenging. We found that these altered catalytic behaviors were not due to changes in the heme reduction rate or in the stability of the enzyme heme-O₂ intermediate, but instead were due to the slower NO k_on and k_d and a slower oxidation rate of the enzyme ferrous heme-NO complex. Computer simulations that utilized the measured kinetic values confirmed this interpretation, and revealed that the V346I iNOS has an enhanced NADPH-dependent NO dioxygenase activity that converts almost 1 NO to nitrate for every NO that the enzyme releases into solution. Together, our results highlight the importance of heme pocket geometry in tuning the NO release versus NO dioxygenase activities of iNOS.     

Nitrosylation of c heme in cd1-nitrite reductase is enhanced during catalysis

The reduction of nitrite into nitric oxide (NO) in denitrifying bacteria is catalyzed by nitrite reductase. In several species, this enzyme is a heme-containing protein with one c heme and one d₁ heme per monomer (cd₁NiR), encoded by the nirS gene.     

Nitric Oxide Dynamics in Truncated Hemoglobin: Docking Sites, Migration Pathways, and Vibrational Spectroscopy from Molecular Dynamics Simulations

Atomistic simulations of nitric oxide (NO) dynamics and migration in the trHbN of Mycobacterium tuberculosis are reported. From extensive molecular dynamics simulations (48 ns in total), the structural and energetic properties of the ligand docking sites in the protein have been characterized and a connectivity network between the ligand docking sites has been built. Several novel migration and exit pathways are found and are analyzed in detail. The interplay between a hydrogen-bonding network involving residues Tyr³³ and Gln⁵⁸ and the bound O₂ ligand is discussed and the role of Phe⁶² residue in ligand migration is examined. It is found that Phe⁶² is directly involved in controlling ligand migration. This is reminiscent of His⁶⁴ in myoglobin, which also plays a central role in CO migration pathways. Finally, infrared spectra of the NO molecule in different ligand docking sites of the protein are calculated. The pocket-specific spectra are typically blue-shifted by 5-10 cm⁻¹, which should be detectable in future spectroscopic experiments.     

Determination of Ligand Pathways in Globins: APOLAR TUNNELS VERSUS POLAR GATES*

Although molecular dynamics simulations suggest multiple interior pathways for O₂ entry into and exit from globins, most experiments indicate well defined single pathways. In 2001, we highlighted the effects of large-to-small amino acid replacements on rates for ligand entry and exit onto the three-dimensional structure of sperm whale myoglobin. The resultant map argued strongly for ligand movement through a short channel from the heme iron to solvent that is gated by the distal histidine (His-64(E7)) near the solvent edge of the porphyrin ring. In this work, we have applied the same mutagenesis mapping strategy to the neuronal mini-hemoglobin from Cerebratulus lacteus (CerHb), which has a large internal tunnel from the heme iron to the C-terminal ends of the E and H helices, a direction that is 180° opposite to the E7 channel. Detailed comparisons of the new CerHb map with expanded results for Mb show unambiguously that the dominant (>90%) ligand pathway in CerHb is through the internal tunnel, and the major (>75%) ligand pathway in Mb is through the E7 gate. These results demonstrate that: 1) mutagenesis mapping can identify internal pathways when they exist; 2) molecular dynamics simulations need to be refined to address discrepancies with experimental observations; and 3) alternative pathways have evolved in globins to meet specific physiological demands.     

Dynamics comparison of two myoglobins with a distinct heme active site

Sperm whale myoglobin (swMb) is a well-studied heme protein, both experimentally and theoretically. Comparatively, little attention has been paid to another member of Mb family, Aplysia limacina myoglobin (apMb). swMb and apMb have the same overall structure and perform the same biological function, i.e., O₂ carrier, while using a distinct heme active site. To provide insights into the structure-function relationship for these two Mbs, we herein made a comparison in terms of their dynamics properties using molecular dynamics simulation. We analyzed the overall structure and protein motions, as well as the intramolecular contacts, namely salt-bridges and hydrogen bonds, especially the interactions between the heme propionate groups and the polypeptide chain. The internal cavities in apMb were also compared to the well-known xenon and other cavities in swMb. Based on current simulations, we propose a unique “arginine gate” for apMb, which has a similar function to the histidine gate observed for swMb in previous studies. This study provides valuable information for understanding the homology of heme proteins, and also aids in rational design of structural and functional heme proteins by alternating the heme active site.     

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

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