Abstract
The superfamily of heme‑copper oxidoreductases (HCOs) include both NO and O(2) 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(2)O). They are structurally similar to heme‑copper oxygen reductases (HCOs), which reduce O(2) 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(8)-10(9) 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(2) diffusion through the closely related O(2) reductases. The apparent second order rate constant approximated using the simulation data is ~5 × 10(8) 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(2) reductases utilize well defined hydrophobic tunnels to assure that substrate diffusion to the buried catalytic sites is not rate limiting under physiological conditions.