Characterization of the bottlenecks and pathways for inhibitor dissociation from [NiFe] hydrogenase

Farzin Sohraby,Ariane Nunes Alves
DOI: https://doi.org/10.26434/chemrxiv-2023-0gs2b-v2
2024-04-24
Abstract:[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O2 and CO. One strategy to engineer [NiFe] hydrogenases and achieve O2 and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site. In this work, we characterized the unbinding pathways of CO in complex with the wild type and 10 different mutants of [NiFe] hydrogenase from Desulfovibrio fructosovorans using τ-Random Accelerated Molecular Dynamics (τRAMD) to enhance the sampling of unbinding events. The ranking provided by the relative residence times computed with τRAMD is in agreement with experiments. Extensive data analysis of the simulations revealed that, from the two bottlenecks proposed in previous studies for the transit of gas molecules (residues 74 and 122, and residues 74 and 476), only one of them (residues 74 and 122) effectively modulates diffusion and residence times for CO. We also computed pathway probabilities for the unbinding of CO, O2 and H2 from the wild type [NiFe] hydrogenase and we observed that, while the most probable pathways are the same, the secondary pathways are different. We propose that introducing mutations to block the most probable paths, in combination with mutations to open the main secondary path used by H2, can be a feasible strategy to achieve CO and O2 resistance in the [NiFe] hydrogenase from Desulfovibrio fructosovorans.
Chemistry
What problem does this paper attempt to address?
The paper primarily investigates the dynamic process of inhibitor (such as CO and O2) desorption from the [NiFe] hydrogenase (an important biocatalyst), particularly how to enhance its tolerance to inhibitors through mutations. The τ-RAMD (τ-random accelerated molecular dynamics) method was employed to simulate and analyze the pathways and timescales of inhibitor dissociation from the active site of the enzyme. By simulating the wild-type and 10 different mutants of [NiFe] hydrogenase, the researchers found that the bottleneck region between amino acid residues 74 and 122 plays a crucial role in the diffusion and residence time of CO. Although these mutations can delay the binding of inhibitors, they do not generate enzymes with complete resistance to inhibitors. By analyzing the dissociation pathways of different gas molecules, they proposed a strategy that could achieve resistance to CO and O2 by blocking the main pathway and opening up secondary pathways. This work provides a theoretical basis for the design of [NiFe] hydrogenases with greater inhibitor tolerance.