Aerobic Damage to [FeFe]-Hydrogenases: Activation Barriers for the Chemical Attachment of O2

A Dynamic Water Channel Affects O2 Stability in [FeFe]-Hydrogenases

[FeFe]-hydrogenases are capable of reducing protons at a high rate. However, molecular oxygen (O2 ) induces the degradation of their catalytic cofactor, the H-cluster, which consists of a cubane [4Fe4S] subcluster (4FeH ) and a unique diiron moiety (2FeH ). Previous attempts to prevent O2 -induced damage have focused on enhancing the protein's sieving effect for O2 by blocking the hydrophobic gas channels that connect the protein surface and the 2FeH . In this study, we aimed to block an O2 diffusion pathway and shield 4FeH instead. Molecular dynamics (MD) simulations identified a novel water channel (WH ) surrounding the H-cluster. As this hydrophilic path may be accessible for O2 molecules we applied site-directed mutagenesis targeting amino acids along WH in proximity to 4FeH to block O2 diffusion. Protein film electrochemistry experiments demonstrate increased O2 stabilities for variants G302S and S357T, and MD simulations based on high-resolution crystal structures confirmed an enhanced local sieving effect for O2 in the environment of the 4FeH in both cases. The results strongly suggest that, in wild type proteins, O2 diffuses from the 4FeH to the 2FeH . These results reveal new strategies for improving the O2 stability of [FeFe]-hydrogenases by focusing on the O2 diffusion network near the active site.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O₂ and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Binding of Reactive Oxygen Species at Fe-S Cubane Clusters

Reactive oxygen species (ROS) play an important role in the biochemistry of the cell and occur in degenerative processes as well as in signal transduction. Iron-sulfur proteins are particularly oxygen-sensitive and their inorganic cofactors frequently undergo ROS-induced decomposition reactions. As experimental knowledge about these processes is still incomplete we present here a quantum chemical study of the relative energetics for the binding of the most relevant ROS to [Fe4S4] clusters. We find that cubane clusters with one uncoordinated Fe atom (as found, for instance, in aconitase) bind all oxygen derivatives considered, whereas activation of triplet O2 to singlet O2 is required for binding to valence-saturated iron centers in these clusters. The radicals NO and OH feature the most exothermic binding energies to Fe atoms. Direct sulfoxidation of coordinating cysteine residues is only possible by OH or H2O2 as attacking agents. The thermodynamic picture of ROS binding to iron-sulfur clusters established here can serve as a starting point for studying reactivity-modulating effects of the cluster-embedding protein environment on ROS-induced decomposition of iron-sulfur proteins.

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

[FeFe]-hydrogenase abundance and diversity along a vertical redox gradient in Great Salt Lake, USA

The use of [FeFe]-hydrogenase enzymes for the biotechnological production of H₂ or other reduced products has been limited by their sensitivity to oxygen (O₂). Here, we apply a PCR-directed approach to determine the distribution, abundance, and diversity of hydA gene fragments along co-varying salinity and O₂ gradients in a vertical water column of Great Salt Lake (GSL), UT. The distribution of hydA was constrained to water column transects that had high salt and relatively low O₂ concentrations. Recovered HydA deduced amino acid sequences were enriched in hydrophilic amino acids relative to HydA from less saline environments. In addition, they harbored interesting variations in the amino acid environment of the complex H-cluster metalloenzyme active site and putative gas transfer channels that may be important for both H₂ transfer and O₂ susceptibility. A phylogenetic framework was created to infer the accessory cluster composition and quaternary structure of recovered HydA protein sequences based on phylogenetic relationships and the gene contexts of known complete HydA sequences. Numerous recovered HydA are predicted to harbor multiple N- and C-terminal accessory iron-sulfur cluster binding domains and are likely to exist as multisubunit complexes. This study indicates an important role for [FeFe]-hydrogenases in the functioning of the GSL ecosystem and provides new target genes and variants for use in identifying O₂ tolerant enzymes for biotechnological applications.