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

Gas Sensing by Bacterial H-NOX Proteins: An MD Study

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 O₂ 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 O₂ > 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.

Gaseous ligand selectivity of the H-NOX sensor protein from Shewanella oneidensis and comparison to those of other bacterial H-NOXs and soluble guanylyl cyclase

To delineate the commonalities and differences in gaseous ligand discrimination among the heme-based sensors with Heme Nitric oxide/OXygen binding protein (H-NOX) scaffold, the binding kinetic parameters for gaseous ligands NO, CO, and O₂, including K_D, k_on, and k_off, of Shewanella oneidensis H-NOX (So H-NOX) were characterized in detail in this study and compared to those of previously characterized H-NOXs from Clostridium botulinum (Cb H-NOX), Nostoc sp. (Ns H-NOX), Thermoanaerobacter tengcongensis (Tt H-NOX), Vibrio cholera (Vc H-NOX), and human soluble guanylyl cyclase (sGC), an H-NOX analogue. The K_D(NO) and K_D(CO) of each bacterial H-NOX or sGC follow the "sliding scale rule"; the affinities of the bacterial H-NOXs for NO and CO vary in a small range but stronger than those of sGC by at least two orders of magnitude. On the other hand, each bacterial H-NOX exhibits different characters in the stability of its 6c NO complex, reactivity with secondary NO, stability of oxyferrous heme and autoxidation to ferric heme. A facile access channel for gaseous ligands is also identified, implying that ligand access has only minimal effect on gaseous ligand selectivity of H-NOXs or sGC. This comparative study of the binding parameters of the bacterial H-NOXs and sGC provides a basis to guide future new structural and functional studies of each specific heme sensor with the H-NOX protein fold.

Entry of water into the distal heme pocket of soluble guanylate cyclase β1 H-NOX domain alters the ligated CO structure: a resonance Raman and in silico simulation study

Water accessing into the heme pocket and alters the structures of CO–sGC (heme), exhibiting two different v_Fe–CO stretching modes.

A QM/MM study of the reaction mechanism of (R)-hydroxynitrile lyases from Arabidopsis thaliana (AtHNL)

Hydroxynitrile lyases (HNLs) catalyze the conversion of chiral cyanohydrins to hydrocyanic acid (HCN) and aldehyde or ketone. Hydroxynitrile lyase from Arabidopsis thaliana (AtHNL) is the first R-selective HNL enzyme containing an α/β-hydrolases fold. In this article, the catalytic mechanism of AtHNL was theoretically studied by using QM/MM approach based on the recently obtained crystal structure in 2012. Two computational models were constructed, and two possible reaction pathways were considered. In Path A, the calculation results indicate that the proton transfer from the hydroxyl group of cyanohydrin occurs firstly, and then the cleavage of C1-C2 bond and the rotation of the generated cyanide ion (CN(-)) follow, afterwards, CN(-) abstracts a proton from His236 via Ser81. The C1-C2 bond cleavage and the protonation of CN(-) correspond to comparable free energy barriers (12.1 vs. 12.2 kcal mol(-1)), suggesting that both of the two processes contribute a lot to rate-limiting. In Path B, the deprotonation of the hydroxyl group of cyanohydrin and the cleavage of C1-C2 bond take place in a concerted manner, which corresponds to the highest free energy barrier of 13.2 kcal mol(-1). The free energy barriers of Path A and B are very similar and basically agree well with the experimental value of HbHNL, a similar enzyme of AtHNL. Therefore, both of the two pathways are possible. In the reaction, the catalytic triad (His236, Ser81, and Asp208) acts as the general acid/base, and the generated CN(-) is stabilized by the hydroxyl group of Ser81 and the main-chain NH-groups of Ala13 and Phe82.