Gas sensing is crucial for both prokaryotes and eukaryotes and is primarily performed by heme-based sensors, including H-NOX domains. These systems may provide a new, alternative mode for transporting gaseous molecules in higher organisms, but for the development of such systems, a detailed understanding of the ligand-binding properties is required. Here, we focused on ligand migration within the protein matrix: we performed molecular dynamics simulations on three bacterial (Ka, Ns and Cs) H-NOX proteins and studied the kinetics of CO, NO and O2 diffusion. We compared the response of the protein structure to the presence of ligands, diffusion rate constants, tunnel systems and storage pockets. We found that the rate constant for diffusion decreases in the O2 > NO > CO order in all proteins, and in the Ns > Ks > Cs order if single-gas is considered. Competition between gases seems to seriously influence the residential time of ligands spent in the distal pocket. The channel system is profoundly determined by the overall fold, but the sidechain pattern has a significant role in blocking certain channels by hydrophobic interactions between bulky groups, cation-π interactions or hydrogen bonding triads. The majority of storage pockets are determined by local sidechain composition, although certain functional cavities, such as the distal and proximal pockets are found in all systems. A major guideline for the design of gas transport systems is the need to chemically bind the gas molecule to the protein, possibly joining several proteins with several heme groups together.

Uranium is harmful to human health due to its radiation damage and the ability of uranyl ion (UO2(2+)) to interact with various proteins and disturb their biological functions. Cytochrome b5 (cyt b5) is a highly negatively charged heme protein and plays a key role in mediating cytochrome c (cyt c) signaling in apoptosis by forming a dynamic cyt b5-cyt c complex. In previous molecular modeling study in combination with UV-Vis studies, we found that UO2(2+) is capable of binding to cyt b5 at surface residues, Glu37 and Glu43. In this study, we further investigated the structural consequences of cyt b5 and cyt c, as well as cyt b5-cyt c complex, upon uranyl binding, by fluorescence spectroscopic and circular dichroism techniques. Moreover, we proposed a uranyl binding site for cyt c at surface residues, Glu66 and Glu69, by performing a molecular modeling study. It was shown that uranyl binds to cyt b5 (KD=10 μM), cyt c (KD=87 μM), and cyt b5-cyt c complex (KD=30 μM) with a different affinity, which slightly alters the protein conformation and disturbs the interaction of cyt b5-cyt c complex. Additionally, we investigated the functional consequences of uranyl binding to the protein surface, which decreases the inherent peroxidase activity of cyt c. The information of uranyl-cyt b5/cyt c interactions gained in this study likely provides a clue for the mechanism of uranyl toxicity.