Computational Understanding of the Selectivities in Metalloenzymes

Quantum chemical studies of the reaction mechanisms of enzymatic CO2 conversion

Enzymatic capture and conversion of carbon dioxide (CO2) into value-added chemicals are of great interest in the field of biocatalysis and have a positive impact on climate change. The quantum chemical methods, recognized as valuable tools for studying reaction mechanisms, have been widely employed in investigating the reaction mechanisms of the enzymes involved in CO2 utilization. In this perspective, we review the mechanistic studies of representative enzymes that are either currently used or have the potential for converting CO2, utilizing the quantum chemical cluster approach and the quantum mechanical/molecular mechanical (QM/MM) method. We begin by summarizing current trends in enzymatic CO2 conversion, followed by a brief description of the computational details of quantum chemical methods. Then, a series of representative examples of the computational modeling of biocatalytic CO2 conversion are presented, including the reduction of CO2 to C1 species (carbon monoxide and formate), and the fixation of CO2 to form aliphatic and aromatic carboxylic acids. The microscopic views of reaction mechanisms obtained from these studies are helpful in guiding the rational design of current enzymes and the discovery of novel enzymes with enhanced performance in converting CO2. Additionally, they provide key information for the de novo design of new-to-nature enzymes. To conclude, we present a perspective on the potential combination of machine learning with quantum description in the study of enzymatic conversion of CO2.

The mechanistic insights into different aspects of promiscuity in metalloenzymes

Enzymes are nature's ultimate machinery to catalyze complex reactions. Though enzymes are evolved to catalyze specific reactions, they also show significant promiscuity in reactions and substrate selection. Metalloenzymes contain a metal ion or metal cofactor in their active site, which is crucial in their catalytic activity. Depending on the metal and its coordination environment, the metal ion or cofactor may function as a Lewis acid or base and a redox center and thus can catalyze a plethora of natural reactions. In fact, the versatility in the oxidation state of the metal ions provides metalloenzymes with a high level of catalytic adaptability and promiscuity. In this chapter, we discuss different aspects of promiscuity in metalloenzymes by using several recent experimental and theoretical works as case studies. We start our discussion by introducing the concept of promiscuity and then we delve into the mechanistic insight into promiscuity at the molecular level.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

The Quantum Chemical Cluster Approach in Biocatalysis

ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.

Wastewater toxicity removal: Integrated chemical and effect-based monitoring of full-scale conventional activated sludge and membrane bioreactor plants

The literature is currently lacking effect-based monitoring studies targeted at evaluating the performance of full-scale membrane bioreactor plants. In this research, a monitoring campaign was performed at a full-scale wastewater treatment facility with two parallel lines (traditional activated sludge and membrane bioreactor). Beside the standard parameters (COD, nitrogen, phosphorus, and metals), 6 polynuclear aromatic hydrocarbons, 29 insecticides, 2 herbicides, and 3 endocrine disrupting compounds were measured. A multi-tiered battery of bioassays complemented the investigation, targeting different toxic modes of action and employing various biological systems (uni/multicellular, prokaryotes/eukaryotes, trophic level occupation). A traffic light scoring approach was proposed to quickly visualize the impact of treatment on overall toxicity that occurred after the exposure to raw and concentrated wastewater. Analysis of the effluents of the CAS and MBR lines show very good performance of the two systems for removal of organic micropollutants and metals. The most noticeable differences between CAS and MBR occurred in the concentration of suspended solids; chemical analyses did not show major differences. On the other hand, bioassays demonstrated better performance for the MBR. Both treatment lines complied with the Italian law's "ecotoxicity standard for effluent discharge in surface water". Yet, residual biological activity was still detected, demonstrating the adequacy and sensitivity of the toxicological tools, which, by their inherent nature, allow the overall effects of complex mixtures to be taken into account.

Radical SAM Enzymes and Metallocofactor Assembly: A Structural Point of View

This Review focuses on the structure-function relationship of radical S-adenosyl-l-methionine (SAM) enzymes involved in the assembly of metallocofactors corresponding to the active sites of [FeFe]-hydrogenase and nitrogenase [MoFe]-protein. It does not claim to correspond to an extensive review on the assembly machineries of these enzyme active sites, for which many good reviews are already available, but instead deals with the contribution of structural data to the understanding of their chemical mechanism (Buren et al. Chem. Rev.2020, 142 ( (25), ) 11006-11012; Britt et al. Chem. Sci.2020, 11 ( (38), ), 10313-10323). Hence, we will present the history and current knowledge about the radical SAM maturases HydE, HydG, and NifB as well as what, in our opinion, should be done in the near future to overcome the existing barriers in our understanding of this fascinating chemistry that intertwine organic radicals and organometallic complexes.

Successes and challenges in multiscale modelling of artificial metalloenzymes: the case study of POP-Rh2 cyclopropanase

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. On top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries require both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allows taking a step further. As a showcase, the enantioselective profiles of the most outstanding variants of an artificial Rh₂-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers (Nat. Commun., 2015, 6, 7789 and Nat. Chem., 2018, 10, 318-324) have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, the latter being the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take the form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

Revealing the Mechanism of Isethionate Sulfite-Lyase by QM/MM Calculations

Isethionate sulfite-lyase (IseG) is a recently characterized glycyl radical enzyme (GRE) that catalyzes radical-mediated C-S bond cleavage of isethionate to produce acetaldehyde and sulfite. Herein, we use quantum mechanical/molecular mechanical (QM/MM) calculations to investigate the detailed catalytic reaction mechanism of IseG. Our calculations indicate that a previously proposed direct 1,2-elimination mechanism is disfavored. Instead, we suggest a new 1,2-migration mechanism for this enzymatic reaction: a key stepwise 1,2-SO₃^- radical migration occurs after the catalytically active cysteinyl radical grabs a hydrogen atom from isethionate, followed by hydrogen atom transfer from cysteine to a 1-hydroxylethane-1-sulfonate radical intermediate. Finally, the elimination of sulfite from 1-hydroxylethane-1-sulfonate to result in the final product is likely to occur outside the enzyme. Glu468 in the active site is found to help orient the substrate rather than grabbing a proton from the hydroxyl group of the substrate. Our findings help reveal the mechanisms of radical-mediated C-S bond cleavage of organosulfonates catalyzed by GREs and expand the understanding of radical-based enzymatic catalysis.

Redox Isomerism in the S3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory

The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S3 state. Our results enabled us to evaluate conflicting models for the S3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Differences in the Nature of the Phosphoryl Transfer Transition State in Protein Phosphatase 1 and Alkaline Phosphatase: Insights from QM Cluster Models

Quantum mechanical (QM) cluster models are used to probe effects on the catalytic properties of protein phosphatase 1 (PP1) and alkaline phosphatase (AP) due to metal ions and active site residues. The calculations suggest that the phosphoryl transfer transition states in PP1 are synchronous in nature with a significant degree of P-Olg cleavage, while those in AP are tighter with a modest degree of P-Olg cleavage and a range of P-Onuc formation. Similar to observations made in our recent work, a significant degree of cross talk between the forming and breaking P-O bonds complicates the interpretation of the Brønsted relation, especially in regard to AP for which the computed βlg/βEQ,lg value does not correlate with the degree of P-Olg cleavage regardless of the metal ions in the active site. By comparison, the correlation between βlg/βEQ,lg and the P-Olg bond order is more applicable to PP1, which generally exhibits less variation in the transition state than AP. Results for computational models with swapped metal ions between PP1 and AP suggest that the metal ions modulate both the nature of the transition state and the degrees of sensitivity of the transition state to the leaving group. In the reactant state, the degree of the scissile bond polarization is also different in the two enzymes, although this difference appears to be largely determined by the active site residues rather than the metal ions. Therefore, both the identity of the metal ion and the positioning of polar or charged residues in the active site contribute to the distinct catalytic characteristics of these enzymes. Several discrepancies observed between the QM cluster results and the available experimental data highlight the need for further QM/MM method developments for the quantitative analysis of metalloenzymes that contain open-shell transition metal ions.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

Formation of Unstable and very Reactive Chemical Species Catalyzed by Metalloenzymes: A Mechanistic Overview

Nature has tailored a wide range of metalloenzymes that play a vast array of functions in all living organisms and from which their survival and evolution depends on. These enzymes catalyze some of the most important biological processes in nature, such as photosynthesis, respiration, water oxidation, molecular oxygen reduction, and nitrogen fixation. They are also among the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions of temperature, pH, and pressure. In the absence of these enzymes, these reactions would proceed very slowly, if at all, suggesting that these enzymes made the way for the emergence of life as we know today. In this review, the structure and catalytic mechanism of a selection of diverse metalloenzymes that are involved in the production of highly reactive and unstable species, such as hydroxide anions, hydrides, radical species, and superoxide molecules are analyzed. The formation of such reaction intermediates is very difficult to occur under biological conditions and only a rationalized selection of a particular metal ion, coordinated to a very specific group of ligands, and immersed in specific proteins allows these reactions to proceed. Interestingly, different metal coordination spheres can be used to produce the same reactive and unstable species, although through a different chemistry. A selection of hand-picked examples of different metalloenzymes illustrating this diversity is provided and the participation of different metal ions in similar reactions (but involving different mechanism) is discussed.

Conversion mechanism of enoyl thioesters into acyl thioesters catalyzed by 2-enoyl-thioester reductases from Candida Tropicalis

Enoyl thioester reductase from Candida tropicalis (Etr1p) catalyzes the NADPH-dependent conversion of enoyl thioesters into acyl thioesters, which are essential in fatty acid and second metabolite biosynthesis.