A microscopic model for gas diffusion dynamics in a [NiFe]-hydrogenase

Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase

[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 the 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.

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.

Exploring the gas access routes in a [NiFeSe] hydrogenase using crystals pressurized with krypton and oxygen

Hydrogenases are metalloenzymes that catalyse both H₂ evolution and uptake. They are gas-processing enzymes with deeply buried active sites, so the gases diffuse through channels that connect the active site to the protein surface. The [NiFeSe] hydrogenases are a special class of hydrogenases containing a selenocysteine as a nickel ligand; they are more catalytically active and less O₂-sensitive than standard [NiFe] hydrogenases. Characterisation of the channel system of hydrogenases is important to understand how the inhibitor oxygen reaches the active site to cause oxidative damage. To this end, crystals of Desulfovibrio vulgaris Hildenborough [NiFeSe] hydrogenase were pressurized with krypton and oxygen, and a method for tracking labile O₂ molecules was developed, for mapping a hydrophobic channel system similar to that of the [NiFe] enzymes as the major route for gas diffusion.

Exploring Conformational Change of Adenylate Kinase by Replica Exchange Molecular Dynamic Simulation

Replica exchange molecular dynamics (REMD) simulation is a popular enhanced sampling method that is widely used for exploring the atomic mechanism of protein conformational change. However, the requirement of huge computational resources for REMD, especially with the explicit solvent model, largely limits its application. In this study, the availability and efficiency of a variant of velocity-scaling REMD (vsREMD) was assessed with adenylate kinase as an example. Although vsREMD achieved results consistent with those from conventional REMD and experimental studies, the number of replicas required for vsREMD (30) was much less than that for conventional REMD (80) to achieve a similar acceptance rate (∼0.2), demonstrating high efficiency of vsREMD to characterize the protein conformational change and associated free-energy profile. Thus, vsREMD is a highly efficient approach for studying the large-scale conformational change of protein systems.

Kinetics of CO2 diffusion in human carbonic anhydrase: a study using molecular dynamics simulations and the Markov-state model

Molecular dynamics (MD) simulations, in combination with the Markov-state model (MSM), were applied to probe CO₂ diffusion from an aqueous solution into the active site of human carbonic anhydrase II (hCA-II), an enzyme useful for enhanced CO₂ capture and utilization. The diffusion process in the hydrophobic pocket of hCA-II was illustrated in terms of a two-dimensional free-energy landscape. We found that CO₂ diffusion in hCA-II is a rate-limiting step in the CO₂ diffusion-binding-reaction process. The equilibrium distribution of CO₂ shows its preferential accumulation within a hydrophobic domain in the protein core region. An analysis of the committors and reactive fluxes indicates that the main pathway for CO₂ diffusion into the active site of hCA-II is through a binding pocket where residue Gln¹³⁶ contributes to the maximal flux. The simulation results offer a new perspective on the CO₂ hydration kinetics and useful insights toward the development of novel biochemical processes for more efficient CO₂ sequestration and utilization.

Ribonucleotide Reductase Requires Subunit Switching in Hypoxia to Maintain DNA Replication

Cells exposed to hypoxia experience replication stress but do not accumulate DNA damage, suggesting sustained DNA replication. Ribonucleotide reductase (RNR) is the only enzyme capable of de novo synthesis of deoxyribonucleotide triphosphates (dNTPs). However, oxygen is an essential cofactor for mammalian RNR (RRM1/RRM2 and RRM1/RRM2B), leading us to question the source of dNTPs in hypoxia. Here, we show that the RRM1/RRM2B enzyme is capable of retaining activity in hypoxia and therefore is favored over RRM1/RRM2 in order to preserve ongoing replication and avoid the accumulation of DNA damage. We found two distinct mechanisms by which RRM2B maintains hypoxic activity and identified responsible residues in RRM2B. The importance of RRM2B in the response to tumor hypoxia is further illustrated by correlation of its expression with a hypoxic signature in patient samples and its roles in tumor growth and radioresistance. Our data provide mechanistic insight into RNR biology, highlighting RRM2B as a hypoxic-specific, anti-cancer therapeutic target.

•

RRM2B is induced in response to hypoxia in both cell models and patient datasets

•

RRM2B retains activity in hypoxic conditions and is the favored RNR subunit in hypoxia

•

Loss of RRM2B has detrimental consequences for cell fate, specifically in hypoxia

•

RRM2B depletion enhanced hypoxic-specific apoptosis and increased radiosensitivity

Structural Plasticity in Globins: Role of Protein Dynamics in Defining Ligand Migration Pathways

Globins are a family of proteins characterized by the presence of the heme prosthetic group and involved in variety of biological functions in the cell. Due to their biological relevance and widespread distribution in all kingdoms of life, intense research efforts have been devoted to disclosing the relationships between structural features, protein dynamics, and function. Particular attention has been paid to the impact of differences in amino acid sequence on the topological features of docking sites and cavities and to the influence of conformational flexibility in facilitating the migration of small ligands through these cavities. Often, tunnels are carved in the interior of globins, and ligand exchange is regulated by gating residues. Understanding the subtle intricacies that relate the differences in sequence with the structural and dynamical features of globins with the ultimate aim of rationalizing the thermodynamics and kinetics of ligand binding continues to be a major challenge in the field. Due to the evolution of computational techniques, significant advances into our understanding of these questions have been made. In this review we focus our attention on the analysis of the ligand migration pathways as well as the function of the structural cavities and tunnels in a series of representative globins, emphasizing the synergy between experimental and theoretical approaches to gain a comprehensive knowledge into the molecular mechanisms of this diverse family of 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.

Aerobic damage to [FeFe]-hydrogenases: activation barriers for the chemical attachment of O2

[FeFe]-hydrogenases are the best natural hydrogen-producing enzymes but their biotechnological exploitation is hampered by their extreme oxygen sensitivity. The free energy profile for the chemical attachment of O2 to the enzyme active site was investigated by using a range-separated density functional re-parametrized to reproduce high-level ab initio data. An activation free-energy barrier of 13 kcal mol(-1) was obtained for chemical bond formation between the di-iron active site and O2, a value in good agreement with experimental inactivation rates. The oxygen binding can be viewed as an inner-sphere electron-transfer process that is strongly influenced by Coulombic interactions with the proximal cubane cluster and the protein environment. The implications of these results for future mutation studies with the aim of increasing the oxygen tolerance of this enzyme are discussed.

Substrate channel in nitrogenase revealed by a molecular dynamics approach

Mo-dependent nitrogenase catalyzes the biological reduction of N2 to two NH3 molecules at FeMo-cofactor buried deep inside the MoFe protein. Access of substrates, such as N2, to the active site is likely restricted by the surrounding protein, requiring substrate channels that lead from the surface to the active site. Earlier studies on crystallographic structures of the MoFe protein have suggested three putative substrate channels. Here, we have utilized submicrosecond atomistic molecular dynamics simulations to allow the nitrogenase MoFe protein to explore its conformational space in an aqueous solution at physiological ionic strength, revealing a putative substrate channel. The viability of this observed channel was tested by examining the free energy of passage of N2 from the surface through the channel to FeMo-cofactor, resulting in the discovery of a very low energy barrier. These studies point to a viable substrate channel in nitrogenase that appears during thermal motions of the protein in an aqueous environment and that approaches a face of FeMo-cofactor earlier implicated in substrate binding.

Uncovering a dynamically formed substrate access tunnel in carbon monoxide dehydrogenase/acetyl-CoA synthase

The transport of small ligands to active sites of proteins is the basis of vital processes in biology such as enzymatic catalysis and cell signaling, but also of more destructive ones including enzyme inhibition and oxidative damage. Here, we show how a diffusion-reaction model solved by means of molecular dynamics and density functional theory calculations provides novel insight into the transport of small ligands in proteins. In particular, we unravel the existence of an elusive, dynamically formed gas channel, which CO2 takes to diffuse from the solvent to the active site (C-cluster) of the bifunctional multisubunit enzyme complex carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS). Two cavities forming this channel are temporarily created by protein fluctuations and are not apparent in the X-ray structures. The ligand transport is controlled by two residues at the end of this tunnel, His113 and His116, and occurs on the same time scale on which chemical binding to the active site takes place (0.1-1 ms), resulting in an overall binding rate on the second time scale. We find that upon reduction of CO2 to CO, the newly formed Fe-hydroxy ligand greatly strengthens the hydrogen-bond network, preventing CO from exiting the protein through the same way that CO2 takes to enter the protein. This is the basis for directional transport of CO from the production site (C-cluster of CODH subunit) to the utilization site (A-cluster of ACS subunit). In view of these results, a general picture emerges of how large proteins guide small ligands toward their active sites.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Molecular basis of [FeFe]-hydrogenase function: an insight into the complex interplay between protein and catalytic cofactor

The precise electrochemical features of metal cofactors that convey the functions of redox enzymes are essentially determined by the specific interaction pattern between cofactor and enclosing protein environment. However, while biophysical techniques allow a detailed understanding of the features characterizing the cofactor itself, knowledge about the contribution of the protein part is much harder to obtain. [FeFe]-hydrogenases are an interesting class of enzymes that catalyze both, H2 oxidation and the reduction of protons to molecular hydrogen with significant efficiency. The active site of these proteins consists of an unusual prosthetic group (H-cluster) with six iron and six sulfur atoms. While H-cluster architecture and catalytic states during the different steps of H2 turnover have been thoroughly investigated during the last 20 years, possible functional contributions from the polypeptide framework were only assumed according to the level of conservancy and X-ray structure analyses. Due to the recent development of simpler and more efficient expression systems the role of single amino acids can now be experimentally investigated. This article summarizes, compares and categorizes the results of recent investigations based on site directed and random mutagenesis according to their informative value about structure function relationships in [FeFe]-hydrogenases. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Structural features of [NiFeSe] and [NiFe] hydrogenases determining their different properties: a computational approach

Hydrogenases are metalloenzymes that catalyze the reversible reaction H(2)<->2H(+) + 2e(-), being potentially useful in H(2) production or oxidation. [NiFeSe] hydrogenases are a particularly interesting subgroup of the [NiFe] class that exhibit tolerance to O(2) inhibition and produce more H(2) than standard [NiFe] hydrogenases. However, the molecular determinants responsible for these properties remain unknown. Hydrophobic pathways for H(2) diffusion have been identified in [NiFe] hydrogenases, as have proton transfer pathways, but they have never been studied in [NiFeSe] hydrogenases. Our aim was, for the first time, to characterize the H(2) and proton pathways in a [NiFeSe] hydrogenase and compare them with those in a standard [NiFe] hydrogenase. We performed molecular dynamics simulations of H(2) diffusion in the [NiFeSe] hydrogenase from Desulfomicrobium baculatum and extended previous simulations of the [NiFe] hydrogenase from Desulfovibrio gigas (Teixeira et al. in Biophys J 91:2035-2045, 2006). The comparison showed that H(2) density near the active site is much higher in [NiFeSe] hydrogenase, which appears to have an alternative route for the access of H(2) to the active site. We have also determined a possible proton transfer pathway in the [NiFeSe] hydrogenase from D. baculatum using continuum electrostatics and Monte Carlo simulation and compared it with the proton pathway we found in the [NiFe] hydrogenase from D. gigas (Teixeira et al. in Proteins 70:1010-1022, 2008). The residues constituting both proton transfer pathways are considerably different, although in the same region of the protein. These results support the hypothesis that some of the special properties of [NiFeSe] hydrogenases could be related to differences in the H(2) and proton pathways.

Mechanistic insight into the blocking of CO diffusion in [NiFe]-hydrogenase mutants through multiscale simulation

[NiFe]-hydrogenases are fascinating biological catalysts with potential application in biofuel cells. However, a severe problem in practical application is the strong sensitivity of hydrogenase to gaseous inhibitor molecules such as CO and O(2). Recently, a number of successful protein engineering studies have been reported that aimed at lowering the access of diatomic inhibitors to the active site pocket, but the molecular mechanism conferring increased resistance remained unclear. Here we use a multiscale simulation approach combining molecular dynamics with a master equation formalism to explain the steady drop in CO diffusion rate observed for the mutants V74M L122A, V74M L122M, and V74M of Desulfovibrio fructosovorans [NiFe]-hydrogenase. We find that diffusion in these variants is controlled by two gates, one between residues 74 and 476 and the other between residues 74 and 122. The existence of two control points in different locations explains why the reduction in the experimental diffusion rate does not simply correlate with the width of the main gas channel. We also find that in the more effective mutation (V74M) CO molecules are still able to reach the active site through transitions that are gated by the microsecond dihedral motions of the side chain of R476 and the thermal fluctuations of the width of the gas channel defined by M74 and L122. Reflecting on the molecular information gained from simulation, we discuss future mutation experiments that could further lower the diffusion rates of small ligands inhibiting [NiFe]-hydrogenase.

Kinetics and thermodynamics of gas diffusion in a NiFe hydrogenase

We have investigated O₂ and H₂ transport across a NiFe hydrogenase at the atomic scale by means of computational methods. The Wild Type protein has been compared with the V74Q mutant. Two distinct methodologies have been applied to study the gas access to the active site. Temperature locally enhanced sampling simulations have emphasized the importance of protein dynamics on gas diffusion. The O₂ diffusion free energy profiles, obtained by umbrella sampling, are in agreement with the known kinetic data and show that in the V74Q mutant, the inhibition process is lowered from both a kinetic and a thermodynamic point of view.

Regioselectivity of H cluster oxidation

The H(2)-evolving potential of [FeFe] hydrogenases is severely limited by the oxygen sensitivity of this class of enzymes. Recent experimental studies on hydrogenase from C. reinhardtii point to O(2)-induced structural changes in the [Fe(4)S(4)] subsite of the H cluster. Here, we investigate the mechanistic basis of this observation by means of density functional theory. Unexpectedly, we find that the isolated H cluster shows a pathological catalytic activity for the formation of reactive oxygen species such as O(2)(-) and HO(2)(-). After protonation of O(2)(-), an OOH radical may coordinate to the Fe atoms of the cubane, whereas H(2)O(2) specifically reacts with the S atoms of the cubane-coordinating cysteine residues. Both pathways are accompanied by significant structural distortions that compromise cluster integrity and thus catalytic activity. These results explain the experimental observation that O(2)-induced inhibition is accompanied by distortions of the [Fe(4)S(4)] moiety and account for the irreversibility of this process.