Evidence of a Thermodynamic Ramp for Hole Hopping to Protect a Redox Enzyme from Oxidative Damage

Tuning the peroxidase activity of artificial P450 peroxygenase by engineering redox-sensitive residues

Cytochrome P450 monooxygenases (P450s) are well recognized as versatile bio-oxidation catalysts. However, the catalytic functions of P450s are highly dependent on NAD(P)H and redox partner proteins. Our group has recently reported the use of a dual-functional small molecule (DFSM) for generating peroxygenase activity of P450BM3, a long-chain fatty acid hydroxylase from Bacillus megaterium. The DFSM-facilitated P450BM3 peroxygenase system exhibited excellent peroxygenation activity and regio-/enantioselectivity for various organic substrates, such as styrenes, thioanisole, small alkanes, and alkylbenzenes. Very recently, we demonstrated that the DFSM-facilitated P450BM3 peroxygenase could be switched to a peroxidase by engineering the redox-sensitive tyrosine residues in P450BM3. Given the great potential of P450 peroxidase for C-H oxyfunctionalization, we herein report scrutiny of the effect of mutating redox-sensitive residues on peroxidase activity by deeply screening all redox-sensitive residues of P450BM3, namely methionines, tryptophans, cysteines, and phenylalanines. As a result, six beneficial mutations at positions M212, F81, M112, F173, M177, and F77 were screened out from 78 constructed mutants, and significantly enhanced the peroxidase activity of P450BM3 in the presence of Im-C6-Phe, a typical DFSM molecule. Further combination of the beneficial mutations resulted in a more than 100-fold improvement in peroxidase activity compared with that of the combined parent enzyme and DFSM, comparable to or better than most natural peroxidases. In addition, mutations of redox-sensitive residues even dramatically increased, by more than 300-fold, the peroxidase activity of the starting F87A enzyme in the absence of the DFSM, despite the far lower apparent catalytic turnover number compared with the DFSM-P450 system. This study provides new insights and a potential strategy for regulating the catalytic promiscuity of P450 enzymes for multiple functional oxidations.

Functionalization of Glucose Oxidase in Organic Solvent: Towards Direct Electrical Communication across Enzyme-Electrode Interface

Enzymatic biosensors based on glucose oxidase has been proven to be one of the effective strategies for the detection of glucose and contributed to health improvements. Therefore, research and debates to date are ongoing in an attempt to find the most effective way to detect this analyte using this enzyme as the recognition center. The 3rd generation biosensors using direct electron transfer (DET) type enzymes are a great way towards practical devices. In this work, we developed a simple method for the functionalization of glucose oxidase with redoxable ferrocene groups in chloroform. The enzyme retained its activity after storage in this organic solvent and after the functionalization procedures. This enzyme functionalization strategy was employed to develop the biosensing monolayer-based platforms for the detection of glucose utilizing the quasi-DET mechanism. As a result of an electrochemical regeneration of the catalytic center, the formation of harmful H₂O₂ is minimized during enzymatic electrocatalysis.

Insights into the Thermodynamics and Kinetics of Amino-Acid Radicals in Proteins

Some oxidoreductase enzymes use redox-active tyrosine, tryptophan, cysteine, and/or glycine residues as one-electron, high-potential redox (radical) cofactors. Amino-acid radical cofactors typically perform one of four tasks-they work in concert with a metallocofactor to carry out a multielectron redox process, serve as storage sites for oxidizing equivalents, activate the substrate molecules, or move oxidizing equivalents over long distances. It is challenging to experimentally resolve the thermodynamic and kinetic redox properties of a single-amino-acid residue. The inherently reactive and highly oxidizing properties of amino-acid radicals increase the experimental barriers further still. This review describes a family of stable and well-structured model proteins that was made specifically to study tyrosine and tryptophan oxidation-reduction. The so-called α₃X model protein system was combined with very-high-potential protein film voltammetry, transient absorption spectroscopy, and theoretical methods to gain a comprehensive description of the thermodynamic and kinetic properties of protein tyrosine and tryptophan radicals.

Functional and protective hole hopping in metalloenzymes

Electrons can tunnel through proteins in microseconds with a modest release of free energy over distances in the 15 to 20 Å range. To span greater distances, or to move faster, multiple charge transfers (hops) are required. When one of the reactants is a strong oxidant, it is convenient to consider the movement of a positively charged "hole" in a direction opposite to that of the electron. Hole hopping along chains of tryptophan (Trp) and tyrosine (Tyr) residues is a critical function in several metalloenzymes that generate high-potential intermediates by reactions with O2 or H2O2, or by activation with visible light. Examination of the protein structural database revealed that Tyr/Trp chains are common protein structural elements, particularly among enzymes that react with O2 and H2O2. In many cases these chains may serve a protective role in metalloenzymes by deactivating high-potential reactive intermediates formed in uncoupled catalytic turnover.

Oxalate decarboxylase uses electron hole hopping for catalysis

The hexameric low-pH stress response enzyme oxalate decarboxylase catalyzes the decarboxylation of the oxalate mono-anion in the soil bacterium Bacillus subtilis. A single protein subunit contains two Mn-binding cupin domains, and catalysis depends on Mn(III) at the N-terminal site. The present study suggests a mechanistic function for the C-terminal Mn as an electron hole donor for the N-terminal Mn. The resulting spatial separation of the radical intermediates directs the chemistry toward decarboxylation of the substrate. A π-stacked tryptophan pair (W96/W274) links two neighboring protein subunits together, thus reducing the Mn-to-Mn distance from 25.9 Å (intrasubunit) to 21.5 Å (intersubunit). Here, we used theoretical analysis of electron hole-hopping paths through redox-active sites in the enzyme combined with site-directed mutagenesis and X-ray crystallography to demonstrate that this tryptophan pair supports effective electron hole hopping between the C-terminal Mn of one subunit and the N-terminal Mn of the other subunit through two short hops of ∼8.5 Å. Replacement of W96, W274, or both with phenylalanine led to a large reduction in catalytic efficiency, whereas replacement with tyrosine led to recovery of most of this activity. W96F and W96Y mutants share the wildtype tertiary structure. Two additional hole-hopping networks were identified leading from the Mn ions to the protein surface, potentially protecting the enzyme from high Mn oxidation states during turnover. Our findings strongly suggest that multistep hole-hopping transport between the two Mn ions is required for enzymatic function, adding to the growing examples of proteins that employ aromatic residues as hopping stations.

Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design

The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.

Cooperative protein-solvent tuning of proton transfer energetics: carbonic anhydrase as a case study

We investigate the coupling between the proton transfer (PT) energetics and the protein-solvent dynamics using the intra-molecular PT in wild type (wt) human carbonic anhydrase II and its ten-fold faster mutant Y7F/N67Q as a test case. We calculate the energy variation upon PT, and from that we also calculate the PT reaction free energy, making use of a hybrid quantum mechanics/molecular dynamics approach. In agreement with the experimental data, we obtain that the reaction free energy is basically the same in the two systems. Yet, we show that the instantaneous PT energy is on average lower in the mutant possibly contributing to the faster PT rate. Analysis of the contribution to the PT energetics of the solvent and of each protein residue, also not in the vicinity of the active site, provides evidence for electrostatic tuning of the PT energy arising from the combined effect of the solvent and the protein environment. These findings open up a way to the more general task of the rational design of mutants with either enhanced or reduced PT rate.

Tuning Proton Transfer Thermodynamics in SARS-Cov-2 Main Protease: Implications for Catalysis and Inhibitor Design

In this comutational work a hybrid quantum mechanics/molecular mechanics approach, the MD-PMM approach, is used to investigate the proton transfer reaction the activates the catalytic activity of SARS-CoV-2 main protease. The proton transfer thermodynamics is investigated for the apo ensyme (i.e., without any bound substrate or inhibitor) and in the presence of a inhibitor, N3, which was previously shown to covalently bind SARS-CoV-2 main protease.

Fully Atomistic Multiscale Approach for pKa Prediction

The ionization state of titratable amino acids strongly affects proteins structure and functioning in a large number of biological processes. It is therefore essential to be able to characterize the pK_a of ionizable groups inside proteins and to understand its microscopic determinants in order to gain insights into many functional properties of proteins. A big effort has been devoted to the development of theoretical approaches for the prediction of deprotonation free energies, yet the accurate theoretical/computational calculation of pK_a values is recognized as a current challenge. A methodology based on a hybrid quantum/classical approach is here proposed for the computation of deprotonation free energies. The method is applied to calculate the pK_a of formic acid, methylammonium, and methanethiol, providing results in good agreement with the corresponding experimental estimates. The pK_a is also calculated for aspartic acid and lysine as single residues in solution and for three aspartic/glutamic acids inside a well-characterized protein: hen egg white lysozyme. While for small molecules the method is able to deal with multiple protonation states of all titratable groups, this becomes computationally very expensive for proteins. The calculated pK_a values for the single amino acids (except for the zwitterionic aspartic acid) and inside the protein display a systematic shift with respect to the experimental values that suggests that the fine balance between hydrophobic and polar interactions might be not accurately reproduced by the usual classical force-fields, thus affecting the computation of deprotonation free energies. The calculated pK_a shifts inside the protein are in good agreement with the corresponding experimental ones (within 1 pK_a unit), well reproducing the pK_a changes due to the protein environment even in the case of large pK_a shifts.