Computation of Rate Constants for Diffusion of Small Ligands to and from Buried Protein Active Sites. Computational Approaches for Studying Enzyme Mechanism Part B. Elsevier; 2016. pp. 299-326. [DOI: 10.1016/bs.mie.2016.05.039]

Gas Sensing by Bacterial H-NOX Proteins: An MD Study

Gas sensing is crucial for both prokaryotes and eukaryotes and is primarily performed by heme-based sensors, including H-NOX domains. These systems may provide a new, alternative mode for transporting gaseous molecules in higher organisms, but for the development of such systems, a detailed understanding of the ligand-binding properties is required. Here, we focused on ligand migration within the protein matrix: we performed molecular dynamics simulations on three bacterial (Ka, Ns and Cs) H-NOX proteins and studied the kinetics of CO, NO and O₂ diffusion. We compared the response of the protein structure to the presence of ligands, diffusion rate constants, tunnel systems and storage pockets. We found that the rate constant for diffusion decreases in the O₂ > NO > CO order in all proteins, and in the Ns > Ks > Cs order if single-gas is considered. Competition between gases seems to seriously influence the residential time of ligands spent in the distal pocket. The channel system is profoundly determined by the overall fold, but the sidechain pattern has a significant role in blocking certain channels by hydrophobic interactions between bulky groups, cation-π interactions or hydrogen bonding triads. The majority of storage pockets are determined by local sidechain composition, although certain functional cavities, such as the distal and proximal pockets are found in all systems. A major guideline for the design of gas transport systems is the need to chemically bind the gas molecule to the protein, possibly joining several proteins with several heme groups together.

O2 sensitivity and H2 production activity of hydrogenases-A review

Hydrogenases are metalloproteins capable of catalyzing the interconversion between molecular hydrogen and protons and electrons. The iron-sulfur clusters within the enzyme enable rapid relay of electrons which are either consumed or generated at the active site. Their unparalleled catalytic efficiency has attracted attention, especially for potential use in H₂ production and/or fuel cell technologies. However, there are limitations to using hydrogenases, especially due to their high O₂ sensitivity. The subclass, called [FeFe] hydrogenases, are particularly more vulnerable to O₂ but proficient in H₂ production. In this review, we provide an overview of mechanistic and protein engineering studies focused on understanding and enhancing O₂ tolerance of the enzyme. The emphasis is on ongoing studies that attempt to overcome O₂ sensitivity of the enzyme while it catalyzes H₂ production in an aerobic environment. We also discuss pioneering attempts to utilize the enzyme in biological H₂ production and other industrial processes, as well as our own perspective on future applications.

Electrochemical Investigations of Hydrogenases and Other Enzymes That Produce and Use Solar Fuels

Many enzymes that produce or transform small molecules such as O₂, H₂, and CO₂ embed inorganic cofactors based on transition metals. Their active site, where the chemical reaction occurs, is buried in and protected by the protein matrix, and connected to the solvent in several ways: chains of redox cofactors mediate long-range electron transfer; static or dynamic tunnels guide the substrate, product and inhibitors; amino acids and water molecules transfer protons. The catalytic mechanism of these enzymes is therefore delocalized over the protein and involves many different steps, some of which determine the response of the enzyme under conditions of stress (extreme redox conditions, presence of inhibitors, light), the catalytic rates in the two directions of the reaction and their ratio (the "catalytic bias"). Understanding all the steps in the catalytic cycle, including those that occur on sites of the protein that are remote from the active site, requires a combination of biochemical, structural, spectroscopic, theoretical, and kinetic methods. Here we argue that kinetics should be used to the fullest extent, by extracting quantitative information from the comparison of data and kinetic models and by exploring the combination of experimental kinetics and theoretical chemistry. In studies of these catalytic mechanisms, direct electrochemistry, the technique which we use and contribute to develop, has become unescapable. It simply consists in monitoring the changes in activity of an enzyme that is wired to an electrode by recording an electric current. We have described kinetic models that can be used to make sense of these data and to learn about various aspects of the mechanism that are difficult to probe using more conventional methods: long-range electron transfer, diffusion along gas channels, redox-driven (in)activations, active site chemistry and photoreactivity under conditions of turnover. In this Account, we highlight a few results that illustrate our approach. We describe how electrochemistry can be used to monitor substrate and inhibitor diffusion along the gas channels of hydrogenases and we discuss how the kinetics of intramolecular diffusion relates to global properties such as resistance to oxygen and catalytic bias. The kinetics and/or thermodynamics of intramolecular electron transfer may also affect the catalytic bias, the catalytic potentials on either side of the equilibrium potential, and the overpotentials for catalysis (defined as the difference between the catalytic potentials and the open circuit potential). This is understood by modeling the shape of the steady-state catalytic response of the enzyme. Other determinants of the catalytic rate, such as domain motions, have been probed by examining the transient catalytic response recorded at fast scan rates. Last, we show that combining electrochemical investigations and MD, DFT, and TD-DFT calculations is an original way of probing the reactivity of the H-cluster of hydrogenase, in particular its reactions with CO, O₂, and light. This approach contrasts with the usual strategy which aims at stabilizing species that are presumed to be catalytic intermediates, and determining their structure using spectroscopic or structural methods.

First Principles Calculation of the Reaction Rates for Ligand Binding to Myoglobin: The Cases of NO and CO

Ligand binding by proteins is among the most fundamental processes in nature. Among these processes the binding of small gas molecules, such as O₂ , CO and NO to heme proteins has traditionally received vivid interest, which was further boosted by their recently recognized significant role in gas sensing in the body. At the heart of the binding of these ligands to the heme group is the spinforbidden reaction between high-spin iron(II) and the ligand yielding a low-spin adduct. We use computational means to address the complete mechanism of CO and NO binding by myoglobin. Considering that it involves several steps occurring on different time scales, molecular dynamics simulations were performed to address the diffusion of the ligand through the enzyme, and DFT calculations in combination with statistical rate calculation to investigate the spin-forbidden reaction. The calculations yielded rate constants in qualitative agreement with experiments and revealed that the bottleneck of NO and CO binding is different; for NO, diffusion was found to be rate-limiting, whereas for CO, the spin-forbidden step is the slowest.

Quantum chemical approaches to [NiFe] hydrogenase

The mechanism by which [NiFe] hydrogenase catalyses the oxidation of molecular hydrogen is a significant yet challenging topic in bioinorganic chemistry. With far-reaching applications in renewable energy and carbon mitigation, significant effort has been invested in the study of these complexes. In particular, computational approaches offer a unique perspective on how this enzyme functions at an electronic and atomistic level. In this article, we discuss state-of-the art quantum chemical methods and how they have helped deepen our comprehension of [NiFe] hydrogenase. We outline the key strategies that can be used to compute the (i) geometry, (ii) electronic structure, (iii) thermodynamics and (iv) kinetic properties associated with the enzymatic activity of [NiFe] hydrogenase and other bioinorganic complexes.