A QM/MM investigation of the activation and catalytic mechanism of Fe-only hydrogenases

Fe-only hydrogenases are enzymes that catalyze dihydrogen production or oxidation, due to the presence of an unusual Fe(6)S(6) cluster (the so-called H-cluster) in their active site, which is composed of a Fe(2)S(2) subsite, directly involved in catalysis, and a classical Fe(4)S(4) cubane cluster. Here, we present a hybrid quantum mechanical and molecular mechanical (QM/MM) investigation of the Fe-only hydrogenase from Desulfovibrio desulfuricans, in order to unravel key issues regarding the activation of the enzyme from its completely oxidized inactive state (Hoxinact) and the influence of the protein environment on the structural and catalytic properties of the H-cluster. Our results show that the Fe(2)S(2) subcluster in the Fe(II)Fe(II) redox state - which is experimentally observed for the completely oxidized form of the enzyme - binds a water molecule to one of its metal centers. The computed QM/MM energy values for water binding to the diferrous subsite are in fact over 70 kJ mol(-1); however, the affinity toward water decreases by 1 order of magnitude after a one-electron reduction of H(ox)(inact), thus leading to the release of coordinated water from the H-cluster. The investigation of a catalytic cycle of the Fe-only hydrogenase that implies formation of a terminal hydride ion and a di(thiomethyl)amine (DTMA) molecule acting as an acid/base catalyst indicates that all steps have reasonable reaction energies and that the influence of the protein on the thermodynamic profile of H(2) production catalysis is not negligible. QM/MM results show that the interactions between the Fe(2)S(2) subsite and the protein environment could give place to structural rearrangements of the H-cluster functional for catalysis, provided that the bidentate ligand that bridges the iron atoms in the binuclear subsite is actually a DTMA residue.     

Binding affinities in the SAMPL3 trypsin and host–guest blind tests estimated with the MM/PBSA and LIE methods

We have estimated affinities for the binding of 34 ligands to trypsin and nine guest molecules to three different hosts in the SAMPL3 blind challenge, using the MM/PBSA, MM/GBSA, LIE, continuum LIE, and Glide score methods. For the trypsin challenge, none of the methods were able to accurately predict the experimental results. For the MM/GB(PB)SA and LIE methods, the rankings were essentially random and the mean absolute deviations were much worse than a null hypothesis giving the same affinity to all ligand. Glide scoring gave a Kendall's τ index better than random, but the ranking is still only mediocre, τ = 0.2. However, the range of affinities is small and most of the pairs of ligands have an experimental affinity difference that is not statistically significant. Removing those pairs improves the ranking metric to 0.4-1.0 for all methods except CLIE. Half of the trypsin ligands were non-binders according to the binding assay. The LIE methods could not separate the inactive ligands from the active ones better than a random guess, whereas MM/GBSA and MM/PBSA were slightly better than random (area under the receiver-operating-characteristic curve, AUC = 0.65-0.68), and Glide scoring was even better (AUC = 0.79). For the first host, MM/GBSA and MM/PBSA reproduce the experimental ranking fairly good, with τ = 0.6 and 0.5, respectively, whereas the Glide scoring was considerably worse, with a τ = 0.4, highlighting that the success of the methods is system-dependent.     

Isocyanide in Biochemistry? A Theoretical Investigation of the Electronic Effects and Energetics of Cyanide Ligand Protonation in [FeFe]‐Hydrogenases

The presence of Fe‐bound cyanide ligands in the active site of the proton‐reducing enzymes [FeFe]‐hydrogenases has led to the hypothesis that such Brønsted–Lowry bases could be protonated during the catalytic cycle, thus implying that hydrogen isocyanide (HNC) might have a relevant role in such crucial microbial metabolic paths. We present a hybrid quantum mechanical/molecular mechanical (QM/MM) study of the energetics of CN^? protonation in the enzyme, and of the effects that cyanide protonation can have on [FeFe]‐hydrogenase active sites. A detailed analysis of the electronic properties of the models and of the energy profile associated with H₂ evolution clearly shows that such protonation is dysfunctional for the catalytic process. However, the inclusion of the protein matrix surrounding the active site in our QM/MM models allowed us to demonstrate that the amino acid environment was finely selected through evolution, specifically to lower the Brønsted–Lowry basicity of the cyanide ligands. In fact, the conserved hydrogen‐bonding network formed by these ligands and the neighboring amino acid residues is able to impede CN^? protonation, as shown by the fact that the isocyanide forms of [FeFe]‐hydrogenases do not correspond to stationary points on the enzyme QM/MM potential‐energy surface.     

Interaction strength and shape difference for the h9/2 and h11/2 configurations in163Tm

The strongly shape driving πh_9/2[541]l/2^? configuration with α=+1/2 exhibits some anomalous, and so far unexplained, features concerning the crossing frequency, ?ω_c, the aligned angular momentum, i_x, and interaction strength, at the alignment of the first pair of i_13/2 quasineutrons in several odd-Z rare earth-nuclei. The h_9/2[541]1/2^? and h_11/2[523]7/2^? bands have been studied in the stably deformed rare-earth nucleus¹⁶³Tm to investigate these features. A difference in band crossing frequency of ～ 80 keV between the two bands is found. Rotational bands built on these two configurations have been found to cross in the spin range I=25/2–29/2 ?. Theγ-decay pattern between the two bands is established in the crossing region and analysed in terms of a moderate shape difference between them. A theoretical estimate of the size of the interaction strength between the two bands is presented and compared to the experimental value. The observed band structure in¹⁶³Tm is very similar to that of¹⁶⁷Lu which has 2 protons and 2 neutrons in addition. This observation is discussed in relation to the similarity of the yrast bands of the two even-even “core” nuclei¹⁶²Er and¹⁶⁶Yb, for which theγ-transition energies are identical within ～0.2 keV below the vi_13/2 crossing.