Dynamics of nitric oxide controlled by protein complex in bacterial system

Spatial transcriptome-guided multi-scale framework connects P. aeruginosa metabolic states to oxidative stress biofilm microenvironment

With the generation of spatially resolved transcriptomics of microbial biofilms, computational tools can be used to integrate this data to elucidate the multi-scale mechanisms controlling heterogeneous biofilm metabolism. This work presents a Multi-scale model of Metabolism In Cellular Systems (MiMICS) which is a computational framework that couples a genome-scale metabolic network reconstruction (GENRE) with Hybrid Automata Library (HAL), an existing agent-based model and reaction-diffusion model platform. A key feature of MiMICS is the ability to incorporate multiple -omics-guided metabolic models, which can represent unique metabolic states that yield different metabolic parameter values passed to the extracellular models. We used MiMICS to simulate Pseudomonas aeruginosa regulation of denitrification and oxidative stress metabolism in hypoxic and nitric oxide (NO) biofilm microenvironments. Integration of P. aeruginosa PA14 biofilm spatial transcriptomic data into a P. aeruginosa PA14 GENRE generated four PA14 metabolic model states that were input into MiMICS. Characteristic of aerobic, denitrification, and oxidative stress metabolism, the four metabolic model states predicted different oxygen, nitrate, and NO exchange fluxes that were passed as inputs to update the agent's local metabolite concentrations in the extracellular reaction-diffusion model. Individual bacterial agents chose a PA14 metabolic model state based on a combination of stochastic rules, and agents sensing local oxygen and NO. Transcriptome-guided MiMICS predictions suggested microscale denitrification and oxidative stress metabolic heterogeneity emerged due to local variability in the NO biofilm microenvironment. MiMICS accurately predicted the biofilm's spatial relationships between denitrification, oxidative stress, and central carbon metabolism. As simulated cells responded to extracellular NO, MiMICS revealed dynamics of cell populations heterogeneously upregulating reactions in the denitrification pathway, which may function to maintain NO levels within non-toxic ranges. We demonstrated that MiMICS is a valuable computational tool to incorporate multiple -omics-guided metabolic models to mechanistically map heterogeneous microbial metabolic states to the biofilm microenvironment.

Are Protein Cavities and Pockets Commonly Used by Redox Active Signalling Molecules?

It has been well known for a long time that inert gases, such as xenon (Xe), have significant biological effects. As these atoms are extremely unlikely to partake in direct chemical reactions with biomolecules such as proteins, lipids, and nucleic acids, there must be some other mode of action to account for the effects reported. It has been shown that the topology of proteins allows for cavities and hydrophobic pockets, and it is via an interaction with such protein structures that inert gases are thought to have their action. Recently, it has been mooted that the relatively inert gas molecular hydrogen (H₂) may also have its effects via such a mechanism, influencing protein structures and actions. H₂ is thought to also act via interaction with redox active compounds, particularly the hydroxyl radical (^·OH) and peroxynitrite (ONOO^-), but not nitric oxide (NO^·), superoxide anions (O₂^·-) or hydrogen peroxide (H₂O₂). However, instead of having a direct interaction with H₂, is there any evidence that these redox compounds can also interact with Xe pockets and cavities in proteins, either having an independent effect on proteins or interfering with the action of inert gases? This suggestion will be explored here.

A 2.2 Å cryoEM structure of a quinol-dependent NO Reductase shows close similarity to respiratory oxidases

Quinol-dependent nitric oxide reductases (qNORs) are considered members of the respiratory heme-copper oxidase superfamily, are unique to bacteria, and are commonly found in pathogenic bacteria where they play a role in combating the host immune response. qNORs are also essential enzymes in the denitrification pathway, catalysing the reduction of nitric oxide to nitrous oxide. Here, we determine a 2.2 Å cryoEM structure of qNOR from Alcaligenes xylosoxidans, an opportunistic pathogen and a denitrifying bacterium of importance in the nitrogen cycle. This high-resolution structure provides insight into electron, substrate, and proton pathways, and provides evidence that the quinol binding site not only contains the conserved His and Asp residues but also possesses a critical Arg (Arg720) observed in cytochrome bo3, a respiratory quinol oxidase.

Trapping of a Mononitrosyl Nonheme Intermediate of Nitric Oxide Reductase by Cryo-Photolysis of Caged Nitric Oxide

Characterization of short-lived reaction intermediates is essential for elucidating the mechanism of the reaction catalyzed by metalloenzymes. Here, we demonstrated that the photolysis of a caged compound under cryogenic temperature followed by thermal annealing is an invaluable technique for trapping of short-lived reaction intermediates of metalloenzymes through the study of membrane-integrated nitric oxide reductase (NOR) that catalyzes reductive coupling of two NO molecules to N₂O at its heme/nonheme Fe_B binuclear center. Although NO produced by the photolysis of caged NO did not react with NOR under cryogenic temperature, annealing to ∼160 K allowed NO to diffuse and react with NOR, which was evident from the appearance of EPR signals assignable to the S = 3/2 state. This indicates that the nonheme Fe_B-NO species can be trapped as the intermediate. Time-resolved IR spectroscopy with the use of the photolysis of caged NO as a reaction trigger showed that the intermediate formed at 10 μs gave the NO stretching frequency at 1683 cm^-1 typical of nonheme Fe-NO, confirming that the combination of the cryo-photolysis of caged NO and annealing enabled us to trap the reaction intermediate. Thus, the cryo-photolysis of the caged compound has great potential for the characterization of short-lived reaction intermediates.

Sequence-specific destabilization of azurin by tetramethylguanidinium-dipeptide ionic liquids

The thermal unfolding of the copper redox protein azurin was studied in the presence of four different dipeptide-based ionic liquids (ILs) utilizing tetramethylguanidinium as the cation. The four dipeptides have different sequences including the amino acids Ser and Asp: TMG-AspAsp, TMG-SerSer, TMG-SerAsp, and TMG-AspSer. Thermal unfolding curves generated from temperature-dependent fluorescence spectroscopy experiments showed that TMG-AspAsp and TMG-SerSer have minor destabilizing effects on the protein while TMG-AspSer and TMG-SerAsp strongly destabilize azurin. Red-shifted fluorescence signatures in the 25 °C correlate with the observed protein destabilization in the solutions with TMG-AspSer and TMG-SerAsp. These signals could correspond to interactions between the Asp residue in the dipeptide and the azurin Trp residue in the unfolded state. These results, supported by appropriate control experiments, suggest that dipeptide sequence-specific interactions lead to selective protein destabilization and motivate further studies of TMG-dipeptide ILs.

•

Four different dipeptide-based ionic liquids were prepared with TMG as cation.

•

The dipeptides included all four sequences containing Asp and Ser.

•

The ionic liquids destabilizing effects on the protein azurin were measured.

•

Different dipeptide sequences have different effects on the protein stability.

CarbonIron Electron Transport Channels in Porphyrin-Graphene Complex for ppb-Level Room Temperature NO Gas Sensing

It is a great challenge to develop efficient room-temperature sensing materials and sensors for nitric oxide (NO) gas, which is a biomarker molecule used in the monitoring of inflammatory respiratory diseases. Herein, Hemin (Fe (III)-protoporphyrin IX) is introduced into the nitrogen-doped reduced graphene oxide (N-rGO) to obtain a novel sensing material HNG-ethanol. Detailed XPS spectra and DFT calculations confirm the formation of carbon-iron bonds in HNG-ethanol during synthesis process, which act as electron transport channels from graphene to Hemin. Owing to this unique chemical structure, HNG-ethanol exhibits superior gas sensing properties toward NO gas (R_a /R_g  = 3.05, 20 ppm) with a practical limit of detection (LOD) of 500 ppb and reliable repeatability (over 5 cycles). The HNG-ethanol sensor also possesses high selectivity against other exhaled gases, high humidity resistance, and stability (less than 3% decrease over 30 days). In addition, a deep understanding of the gas sensing mechanisms is proposed for the first time in this work, which is instructive to the community for fabricating sensing materials based on graphene-iron derivatives in the future.

Potential of Time-Resolved Serial Femtosecond Crystallography Using High Repetition Rate XFEL Sources

This perspective review describes emerging techniques and future opportunities for time-resolved serial femtosecond crystallography (TR-SFX) experiments using high repetition rate XFEL sources. High repetition rate sources are becoming more available with the European XFEL in operation and the recently upgraded LCLS-II will be available in the near future. One efficient use of these facilities for TR-SFX relies on pump–probe experiments using a laser to trigger a reaction of light-responsive proteins or mix-and-inject experiments for light-unresponsive proteins. With the view to widen the application of TR-SFX, the promising field of photocaged compounds is under development, which allows the very fast laser triggering of reactions that is no longer limited to naturally light-responsive samples. In addition to reaction triggering, a key concern when performing an SFX experiment is efficient sample usage, which is a main focus of new high repetition rate-compatible sample delivery methods.

The Biologically Relevant Coordination Chemistry of Iron and Nitric Oxide: Electronic Structure and Reactivity

Nitric oxide (NO) is an important signaling molecule that is involved in a wide range of physiological and pathological events in biology. Metal coordination chemistry, especially with iron, is at the heart of many biological transformations involving NO. A series of heme proteins, nitric oxide synthases (NOS), soluble guanylate cyclase (sGC), and nitrophorins, are responsible for the biosynthesis, sensing, and transport of NO. Alternatively, NO can be generated from nitrite by heme- and copper-containing nitrite reductases (NIRs). The NO-bearing small molecules such as nitrosothiols and dinitrosyl iron complexes (DNICs) can serve as an alternative vehicle for NO storage and transport. Once NO is formed, the rich reaction chemistry of NO leads to a wide variety of biological activities including reduction of NO by heme or non-heme iron-containing NO reductases and protein post-translational modifications by DNICs. Much of our understanding of the reactivity of metal sites in biology with NO and the mechanisms of these transformations has come from the elucidation of the geometric and electronic structures and chemical reactivity of synthetic model systems, in synergy with biochemical and biophysical studies on the relevant proteins themselves. This review focuses on recent advancements from studies on proteins and model complexes that not only have improved our understanding of the biological roles of NO but also have provided foundations for biomedical research and for bio-inspired catalyst design in energy science.

Thermodynamic destabilization of azurin by four different tetramethylguanidinium amino acid ionic liquids

The thermal unfolding of the copper redox protein azurin was studied in the presence of four different amino acid-based ionic liquids (ILs), all of which have tetramethylguanidium as cation. The anionic amino acid includes two with alcohol side chains, serine and threonine, and two with carboxylic acids, aspartate and glutamate. Control experiments showed that amino acids alone do not significantly change protein stability and pH changes anticipated by the amino acid nature have only minor effects on the protein. With the ILs, the protein is destabilized and the melting temperature is decreased. The two ILs with alcohol side chains strongly destabilize the protein while the two ILs with acid side chains have weaker effects. Unfolding enthalpy (ΔH_unf^°) and entropy (ΔS_unf^°) values, derived from fits of the unfolding data, show that some ILs increase ΔH_unf^°while others do not significantly change this value. All ILs, however, increase ΔS_unf^°. MD simulations of both the folded and unfolded protein conformations in the presence of the ILs provide insight into the different IL-protein interactions and how they affect the ΔH_unf^° values. The simulations also confirm that the ILs increase the unfolded state entropies which can explain the increased ΔS_unf^° values.

Crystal structure of a photosynthetic LH1-RC in complex with its electron donor HiPIP

Photosynthetic electron transfers occur through multiple components ranging from small soluble proteins to large integral membrane protein complexes. Co-crystallization of a bacterial photosynthetic electron transfer complex that employs weak hydrophobic interactions was achieved by using high-molar-ratio mixtures of a soluble donor protein (high-potential iron-sulfur protein, HiPIP) with a membrane-embedded acceptor protein (reaction center, RC) at acidic pH. The structure of the co-complex offers a snapshot of a transient bioenergetic event and revealed a molecular basis for thermodynamically unfavorable interprotein electron tunneling. HiPIP binds to the surface of the tetraheme cytochrome subunit in the light-harvesting (LH1) complex-associated RC in close proximity to the low-potential heme-1 group. The binding interface between the two proteins is primarily formed by uncharged residues and is characterized by hydrophobic features. This co-crystal structure provides a model for the detailed study of long-range trans-protein electron tunneling pathways in biological systems.

The high potential iron-sulfur (HiPIP) proteins are direct electron donors to the light-harvesting-reaction center complexes (LH1-RC) in photosynthetic β- and γ-Proteobacteria. Here, the authors present the 2.9 Å crystal structure of the HiPIP-bound LH1-RC complex from the thermophilic purple sulfur bacterium Thermochromatium tepidum and discuss mechanistic implications for the electron transfer pathway.

Current Knowledge and Future Directions in Developing Strategies to Combat Pseudomonas aeruginosa Infection

In the face of growing antimicrobial resistance, there is an urgent need for the development of effective strategies to target Pseudomonas aeruginosa. This metabolically versatile bacterium can cause a wide range of severe opportunistic infections in patients with serious underlying medical conditions, such as those with burns, surgical wounds or people with cystic fibrosis. Many of the key adaptations that arise in this organism during infection are centered on core metabolism and virulence factor synthesis. Interfering with these processes may provide a new strategy to combat infection which could be combined with conventional antibiotics. This review will provide an overview of the most recent work that has advanced our understanding of P. aeruginosa infection. Strategies that exploit this recent knowledge to combat infection will be highlighted alongside potential alternative therapeutic options and their limitations.

Biological and Bioinspired Inorganic N-N Bond-Forming Reactions

The metallobiochemistry underlying the formation of the inorganic N-N-bond-containing molecules nitrous oxide (N2O), dinitrogen (N2), and hydrazine (N2H4) is essential to the lifestyles of diverse organisms. Similar reactions hold promise as means to use N-based fuels as alternative carbon-free energy sources. This review discusses research efforts to understand the mechanisms underlying biological N-N bond formation in primary metabolism and how the associated reactions are tied to energy transduction and organismal survival. These efforts comprise studies of both natural and engineered metalloenzymes as well as synthetic model complexes.

Functional interactions between nitrite reductase and nitric oxide reductase from Paracoccus denitrificans

Denitrification is a microbial pathway that constitutes an important part of the nitrogen cycle on earth. Denitrifying organisms use nitrate as a terminal electron acceptor and reduce it stepwise to nitrogen gas, a process that produces the toxic nitric oxide (NO) molecule as an intermediate. In this work, we have investigated the possible functional interaction between the enzyme that produces NO; the cd₁ nitrite reductase (cd₁NiR) and the enzyme that reduces NO; the c-type nitric oxide reductase (cNOR), from the model soil bacterium P. denitrificans. Such an interaction was observed previously between purified components from P. aeruginosa and could help channeling the NO (directly from the site of formation to the side of reduction), in order to protect the cell from this toxic intermediate. We find that electron donation to cNOR is inhibited in the presence of cd₁NiR, presumably because cd₁NiR binds cNOR at the same location as the electron donor. We further find that the presence of cNOR influences the dimerization of cd₁NiR. Overall, although we find no evidence for a high-affinity, constant interaction between the two enzymes, our data supports transient interactions between cd₁NiR and cNOR that influence enzymatic properties of cNOR and oligomerization properties of cd₁NiR. We speculate that this could be of particular importance in vivo during metabolic switches between aerobic and denitrifying conditions.

Bioinspired oxidation of oximes to nitric oxide with dioxygen by a nonheme iron(II) complex

The ability of two iron(II) complexes, [(Tp^Ph2)Fe^II(benzilate)] (1) and [(Tp^Ph2)(Fe^II)₂(NPP)₃] (2) (Tp^Ph2 = hydrotris(3,5-diphenylpyrazol-1-yl)borate, NPP-H = α-isonitrosopropiophenone), of a monoanionic facial N3 ligand in the O₂-dependent oxidation of oximes is reported. The mononuclear complex 1 reacts with dioxygen to decarboxylate the iron-coordinated benzilate. The oximate-bridged dinuclear complex (2), which contains a high-spin (Tp^Ph2)Fe^II unit and a low-spin iron(II)-oximate unit, activates dioxygen at the high-spin iron(II) center. Both the complexes exhibit the oxidative transformation of oximes to the corresponding carbonyl compounds with the incorporation of one oxygen atom from dioxygen. In the oxidation process, the oxime units are converted to nitric oxide (NO) or nitroxyl (HNO). The iron(II)-benzilate complex (1) reacts with oximes to afford HNO, whereas the iron(II)-oximate complex (2) generates NO. The results described here suggest that the oxidative transformation of oximes to NO/HNO follows different pathways depending upon the nature of co-ligand/reductant.Graphic abstract.

Dimeric structures of quinol-dependent nitric oxide reductases (qNORs) revealed by cryo-electron microscopy

Quinol-dependent nitric oxide reductases (qNORs) are membrane-integrated, iron-containing enzymes of the denitrification pathway, which catalyze the reduction of nitric oxide (NO) to the major ozone destroying gas nitrous oxide (N₂O). Cryo-electron microscopy structures of active qNOR from Alcaligenes xylosoxidans and an activity-enhancing mutant have been determined to be at local resolutions of 3.7 and 3.2 Å, respectively. They unexpectedly reveal a dimeric conformation (also confirmed for qNOR from Neisseria meningitidis) and define the active-site configuration, with a clear water channel from the cytoplasm. Structure-based mutagenesis has identified key residues involved in proton transport and substrate delivery to the active site of qNORs. The proton supply direction differs from cytochrome c-dependent NOR (cNOR), where water molecules from the cytoplasm serve as a proton source similar to those from cytochrome c oxidase.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Changes in the abundance and community composition of different nitrogen cycling groups in response to fumigation with 1,3-dichloropropene

The fumigant 1,3-dichloropropene (1,3-D) is widely-used to control pathogenic bacteria, fungi, nematodes and insects in soil before a crop is planted. Although fumigants in general have been reported to have a 'fertilizer effect' in the soil by increasing nitrogen availability, little is known of how a specific fumigant such as 1,3-D affects available nitrogen. This study used real-time quantitative PCR (qPCR) and 16S rRNA gene amplicon sequencing techniques to investigate the effects of 1,3-D on microorganisms involved in nitrogen cycling that were present in 2 soils: Jiangxi lateritic red soil and Beijing fluvo-aquic soil. The fumigant 1,3-D temporarily decreased the abundance of 11 functional genes involved in nitrogen-fixing, nitrification and denitrification in both soil types. Different nitrogen cycling groups recovered to the unfumigated level in various incubation phases. Microorganisms containing nifH, nxrB, napA and qnorB genes were most vulnerable to 1,3-D fumigation. However, a stronger and longer inhibition effect of 1,3-D on these 11 functional genes was observed in Jiangxi soil than in Beijing soil. At the same time, the abundance of nifH, AOBamoA, nirS, qnorB and cnorB genes was significantly increased 59 days after 1,3-D fumigation. Fumigation with 1,3-D significantly reduced the nitrogen-fixing bacteria Azospirillum and Paenibacillus; the nitrifiers Nitrosomonas and Nitrospira; and the denitrifiers Pseudomonas, Paracoccus and Sphingomonas. Conversely, fumigation with 1,3-D increased the nitrogen-fixing bacteria Bradyrhizobium and Rhizobium; the nitrification bacteria Nitrosospira and Nitrolancea; and the denitrification bacteria Sphingobium, Alcanivorax, Bacillus, Streptomyces and Aeromonas. Fumigation with 1,3-D therefore caused significant shifts in the species composition and number of microbes directly involved in nitrogen cycling in the short-term. These results contribute toward a better understanding of the impact of 1,3-D fumigation on various types of soil nitrogen-cycling groups.

Responses of Nitrogen-Cycling Microorganisms to Dazomet Fumigation

The influence of soil fumigation on microorganisms involved in transforming nitrogen remains little understood, despite the use of fumigants for many decades to control soil-borne pathogens and plant-parasitic nematodes. We used real-time PCR (quantitative polymerase chain reaction) and 16S rRNA gene amplicon sequencing techniques to monitor changes in the diversity and community structure of microorganisms associated with nitrogen transfer after the soil was fumigated with dazomet (DZ). We also examined nitrous oxide (N₂O) emissions from these microorganisms present in fumigated fluvo-aquic soil and lateritic red soil. Fumigation with DZ significantly reduced the abundance of 16S rRNA and nitrogen cycling functional genes (nifH, AOA amoA, AOB amoA, nxrB, narG, napA, nirK, nirS, cnorB, qnorB, and nosZ). At the same time, N₂O production rates increased between 9.9 and 30 times after fumigation. N₂O emissions were significantly correlated with NH 4 + , dissolved amino acids and microbial biomass nitrogen, but uncorrelated with functional gene abundance. Diversity indices showed that DZ temporarily stimulated bacterial diversity as well as caused a significant change in bacterial community composition. For example, DZ significantly decreased populations of N₂-fixing bacteria Mesorhizobium and Paenibacillus, nitrifiers Nitrosomonas, and the denitrifiers Bacillus, Pseudomonas, and Paracoccus. The soil microbial community had the ability to recover to similar population levels recorded in unfumigated soils when the inhibitory effects of DZ fumigation were no longer evident. The microbial recovery rate, however, depended on the physicochemical properties of the soil. These results provided useful information for environmental safety assessments of DZ in China, for improving our understanding of the N-cycling pathways in fumigated soils, and for determining the potential responses of different N-cycling groups after fumigation.

Surface-Enhanced Infrared Absorption Spectroscopy of Bacterial Nitric Oxide Reductase under Electrochemical Control Using a Vibrational Probe of Carbon Monoxide

Nitric oxide reductases (NORs) reduce nitric oxide to nitrous oxide in the denitrification pathway of the global nitrogen cycle. NORs contain four iron cofactors and the NO reduction occurs at the heme b₃/nonheme Fe_B binuclear active site. The determination of reduction potentials of the iron cofactors will help us elucidate the enzymatic reaction mechanism. However, previous reports on these potentials remain controversial. Herein, we performed electrochemical and surface-enhanced infrared absorption (SEIRA) spectroscopic measurements of Pseudomonas aeruginosa NOR immobilized on gold electrodes. Cyclic voltammograms exhibited two reduction peaks at -0.11 and -0.44 V vs SHE, and a SEIRA spectrum using a vibrational probe of CO showed a characteristic band at 1972 cm^-1 at -0.4 V vs SHE, which was assigned to νCO of heme b₃-CO. Our results suggest that the reduction of heme b₃ initiates the enzymatic NO reduction.

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.