Protonation states of intermediates in the reaction mechanism of [NiFe] hydrogenase studied by computational methods.

TY - JOUR

T1 - Protonation states of intermediates in the reaction mechanism of [NiFe] hydrogenase studied by computational methods.

AU - Dong, Geng

AU - Ryde, Ulf

PY - 2016

Y1 - 2016

N2 - The [NiFe] hydrogenases catalyse the reversible conversion of H2 to protons and electrons. The active site consists of a Fe ion with one carbon monoxide, two cyanide, and two cysteine (Cys) ligands. The latter two bridge to a Ni ion, which has two additional terminal Cys ligands. It has been suggested that one of the Cys residues is protonated during the reaction mechanism. We have used combined quantum mechanical and molecular mechanics (QM/MM) geometry optimisations, large QM calculations with 817 atoms, and QM/MM free energy simulations, using the TPSS and B3LYP methods with basis sets extrapolated to the quadruple zeta level to determine which of the four Cys residues is more favourable to protonate for four putative states in the reaction mechanism, Ni-SIa, Ni-R, Ni-C, and Ni-L. The calculations show that for all states, the terminal Cys-546 residue is most easily protonated by 14-51 kJ/mol, owing to a more favourable hydrogen-bond pattern around this residue in the protein.

AB - The [NiFe] hydrogenases catalyse the reversible conversion of H2 to protons and electrons. The active site consists of a Fe ion with one carbon monoxide, two cyanide, and two cysteine (Cys) ligands. The latter two bridge to a Ni ion, which has two additional terminal Cys ligands. It has been suggested that one of the Cys residues is protonated during the reaction mechanism. We have used combined quantum mechanical and molecular mechanics (QM/MM) geometry optimisations, large QM calculations with 817 atoms, and QM/MM free energy simulations, using the TPSS and B3LYP methods with basis sets extrapolated to the quadruple zeta level to determine which of the four Cys residues is more favourable to protonate for four putative states in the reaction mechanism, Ni-SIa, Ni-R, Ni-C, and Ni-L. The calculations show that for all states, the terminal Cys-546 residue is most easily protonated by 14-51 kJ/mol, owing to a more favourable hydrogen-bond pattern around this residue in the protein.

U2 - 10.1007/s00775-016-1348-9

DO - 10.1007/s00775-016-1348-9

M3 - Article

C2 - 26940957

VL - 21

SP - 383

EP - 394

JO - Journal of Biological Inorganic Chemistry

JF - Journal of Biological Inorganic Chemistry

SN - 1432-1327

IS - 3

ER -