The Importance of the Scaffold for de Novo Enzymes: A Case Study with Kemp Eliminase

Tuning Acid-Base Chemistry at an Electrified Gold/Water Interface

Acid-base reactions are ubiquitous in solution chemistry, as well as in electrochemistry. However, macroscopic concepts derived in solutions, such as pKa and pH, differ significantly at electrified metal-aqueous interfaces due to specific solvation and applied voltage. Here, we measure the pKa values of an amino acid, glycine, at a gold/water interface under a varying applied voltage by means of spectroscopic titration. With the help of simulations, we propose a general model to understand potential-dependent shifts in pKa values in terms of local hydrophobicity and electric fields. These parameters can be tuned by adjusting the metal surface and applied voltage, respectively, offering promising, but still unexplored, paths to regulate reactivity. Our results change the focus with respect to common interpretations based on, for example, apparent local pH effects and open interesting perspectives for electrochemical reaction steering.

Electrostatics as a Guiding Principle in Understanding and Designing Enzymes

Enzyme design faces challenges related to the implementation of the basic principles that govern the catalytic activity in natural enzymes. In this work, we revisit basic electrostatic concepts that have been shown to explain the origin of enzymatic efficiency like preorganization and reorganization. Using magnitudes such as the electrostatic potential and the electric field generated by the protein, we explain how these concepts work in different enzymes and how they can be used to rationalize the consequences of point mutations. We also discuss examples of protein design in which electrostatic effects have been implemented. For the near future, molecular simulations, coupled with the use of machine learning methods, can be used to implement electrostatics as a guiding principle for enzyme designs.

Substrate Positioning Dynamics Involves a Non-Electrostatic Component to Mediate Catalysis

Substrate positioning dynamics (SPD) orients the substrate in the active site, thereby influencing catalytic efficiency. However, it remains unknown whether SPD effects originate primarily from electrostatic perturbation inside the enzyme or can independently mediate catalysis with a significant non-electrostatic component. In this work, we investigated how the non-electrostatic component of SPD affects transition state (TS) stabilization. Using high-throughput enzyme modeling, we selected Kemp eliminase variants with similar electrostatics inside the enzyme but significantly different SPD. The kinetic parameters of these mutants were experimentally characterized. We observed a valley-shaped, two-segment linear correlation between the TS stabilization free energy (converted from kinetic parameters) and substrate positioning index (a metric to quantify SPD). The energy varies by approximately 2 kcal/mol. Favorable SPD was observed for the distal mutant R154W, increasing the proportion of reactive conformations and leading to the lowest activation free energy. These results indicate the substantial contribution of the non-electrostatic component of SPD to enzyme catalytic efficiency.

Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions

In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

An algorithm to find the optimal oriented external electrostatic field for annihilating a reaction barrier in a polarizable molecular system

The use of oriented external electric fields (OEEFs) to promote and control chemical reactivity has motivated many theoretical and computational studies in the last decade to model the action of OEEFs on a molecular system and its effects on chemical processes. Given a reaction, a central goal in this research area is to predict the optimal OEEF (oOEEF) required to annihilate the reaction energy barrier with the smallest possible field strength. Here, we present a model rooted in catastrophe and optimum control theories that allows us to find the oOEEF for a given reaction valley in the potential energy surface (PES). In this model, the effective (or perturbed) PES of a polarizable molecular system is constructed by adding to the original, non-perturbed, PES a term accounting for the interaction of the OEEF with the intrinsic electric dipole and polarizability of the molecular system, so called the polarizable molecular electric dipole (PMED) model. We demonstrate that the oOEEF can be established by locating a point in the original PES with unique topological properties: the optimal barrier breakdown or bond-breaking point (oBBP). The essential feature of the oBBP structure is the fact that this point maintains its topological properties for all the applied OEEFs, also for the unperturbed PES, thus becoming much more relevant than the commonly used minima and transition state structures. The PMED model proposed here has been implemented in an open access package and is shown to successfully predict the oOEEF for two processes: an isomerization reaction of a cumulene derivative and the Huisgen cycloaddition reaction.

EnzyHTP Computational Directed Evolution with Adaptive Resource Allocation

Directed evolution facilitates enzyme engineering via iterative rounds of mutagenesis. Despite the wide applications of high-throughput screening, building "smart libraries" to effectively identify beneficial variants remains a major challenge in the community. Here, we developed a new computational directed evolution protocol based on EnzyHTP, a software that we have previously reported to automate enzyme modeling. To enhance the throughput efficiency, we implemented an adaptive resource allocation strategy that dynamically allocates different types of computing resources (e.g., GPU/CPU) based on the specific need of an enzyme modeling subtask in the workflow. We implemented the strategy as a Python library and tested the library using fluoroacetate dehalogenase as a model enzyme. The results show that compared to fixed resource allocation where both CPU and GPU are on-call for use during the entire workflow, applying adaptive resource allocation can save 87% CPU hours and 14% GPU hours. Furthermore, we constructed a computational directed evolution protocol under the framework of adaptive resource allocation. The workflow was tested against two rounds of mutational screening in the directed evolution experiments of Kemp eliminase (KE07) with a total of 184 mutants. Using folding stability and electrostatic stabilization energy as computational readout, we identified all four experimentally observed target variants. Enabled by the workflow, the entire computation task (i.e., 18.4 μs MD and 18,400 QM single-point calculations) completes in 3 days of wall-clock time using ∼30 GPUs and ∼1000 CPUs.

Evaluating the Transition State Stabilization/Destabilization Effects of the Electric Fields from Scaffold Residues by a QM/MM Approach

The protein scaffolds of enzymes not only provide structural support for the catalytic center but also exert preorganized electric fields for electrostatic catalysis. In recent years, uniform oriented external electric fields (OEEFs) have been widely applied to enzymatic reactions to mimic the electrostatic effects of the environment. However, the electric fields exerted by individual residues in proteins may be quite heterogeneous across the active site, with varying directions and strengths at different positions of the active site. Here, we propose a QM/MM-based approach to evaluate the effects of the electric fields exerted by individual residues in the protein scaffold. In particular, the heterogeneity of the residue electric fields and the effect of the native protein environment can be properly accounted for by this QM/MM approach. A case study of the O-O heterolysis reaction in the catalytic cycle of TyrH shows that (1) for scaffold residues that are relatively far from the active site, the heterogeneity of the residue electric field in the active site is not very significant and the electrostatic stabilization/destabilization due to each residue can be well approximated with the interaction energy between a uniform electric field and the QM region dipole; (2) for scaffold residues near the active site, the residue electric fields can be highly heterogeneous along the breaking O-O bond. In such a case, approximating the residue electric fields as uniform fields may misrepresent the overall electrostatic effect of the residue. The present QM/MM approach can be applied to evaluate the residues' electrostatic impact on enzymatic reactions, which also can be useful in computational optimization of electric fields to boost the enzyme catalysis.

Geometry of Charge Density as a Reporter on the Role of the Protein Scaffold in Enzymatic Catalysis: Electrostatic Preorganization and Beyond

Enzymes host active sites inside protein macromolecules, which have diverse, often incredibly complex, and atom-expensive structures. It is an outstanding question what the role of these expensive scaffolds might be in enzymatic catalysis. Answering this question is essential to both enzymology and the design of artificial enzymes with proficiencies that will match those of the best natural enzymes. Protein rigidifying the active site, contrasted with the dynamics and vibrational motion promoting the reaction, as well as long-range electrostatics (also known as electrostatic preorganization) were all proposed as central contributions of the scaffold to the catalysis. Here, we show that all these effects inevitably produce changes in the quantum mechanical electron density in the active site, which in turn defines the reactivity. The phenomena are therefore fundamentally inseparable. The geometry of the electron density-a scalar field characterized by a number of mathematical features such as critical points-is a rigorous and convenient descriptor of enzymatic catalysis and a reporter on the role of the protein. We show how this geometry can be analyzed, linked to the reaction barriers, and report in particular on intramolecular electric fields in enzymes. We illustrate these tools on the studies of electrostatic preorganization in several representative enzyme classes, both natural and artificial. We highlight the forward-looking aspects of the approach.

A catastrophe theory-based model for optimal control of chemical reactions by means of oriented electric fields

The effect of oriented external electric fields (OEEF) on chemical reactivity has been studied theoretically and computationally in the last decades. A central goal in this research area is to predict the orientation and the smallest amplitude electric field that renders a barrierless chemical process with the smallest possible strength. Recently, a model to find the optimal electric field has been proposed and described (Bofill JM et al., J. Chem. Theory Comput. 18:935, 2022). We here proof that this model is based on catastrophe and optimum control theories. Based on both theories a technical treatment of the model is given and applied to a two-dimensional generic example that provides insight into its nature and capability. Finally, the model is applied to determine the optimal OEEF for the trans-to-cis isomerization of a [3]cumulene derivative.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

Natural Evolution Provides Strong Hints about Laboratory Evolution of Designer Enzymes

Laboratory evolution combined with computational enzyme design provides the opportunity to generate novel biocatalysts. Nevertheless, it has been challenging to understand how laboratory evolution optimizes designer enzymes by introducing seemingly random mutations. A typical enzyme optimized with laboratory evolution is the abiological Kemp eliminase, initially designed by grafting active site residues into a natural protein scaffold. Here, we relate the catalytic power of laboratory-evolved Kemp eliminases to the statistical energy ([Formula: see text]) inferred from their natural homologous sequences using the maximum entropy model. The [Formula: see text] of designs generated by directed evolution is correlated with enhanced activity and reduced stability, thus displaying a stability-activity trade-off. In contrast, the [Formula: see text] for mutants in catalytic-active remote regions (in which remote residues are important for catalysis) is strongly anticorrelated with the activity. These findings provide an insight into the role of protein scaffolds in the adaption to new enzymatic functions. It also indicates that the valley in the [Formula: see text] landscape can guide enzyme design for abiological catalysis. Overall, the connection between laboratory and natural evolution contributes to understanding what is optimized in the laboratory and how new enzymatic function emerges in nature, and provides guidance for computational enzyme design.

Computational Modeling of Supramolecular Metallo-organic Cages-Challenges and Opportunities

Self-assembled metallo-organic cages have emerged as promising biomimetic platforms that can encapsulate whole substrates akin to an enzyme active site. Extensive experimental work has enabled access to a variety of structures, with a few notable examples showing catalytic behavior. However, computational investigations of metallo-organic cages are scarce, not least due to the challenges associated with their modeling and the lack of accurate and efficient protocols to evaluate these systems. In this review, we discuss key molecular principles governing the design of functional metallo-organic cages, from the assembly of building blocks through binding and catalysis. For each of these processes, computational protocols will be reviewed, considering their inherent strengths and weaknesses. We will demonstrate that while each approach may have its own specific pitfalls, they can be a powerful tool for rationalizing experimental observables and to guide synthetic efforts. To illustrate this point, we present several examples where modeling has helped to elucidate fundamental principles behind molecular recognition and reactivity. We highlight the importance of combining computational and experimental efforts to speed up supramolecular catalyst design while reducing time and resources.

Molecular Dynamics-Derived Descriptor Informs the Impact of Mutation on the Catalytic Turnover Number in Lactonase Across Substrates

Molecular dynamics simulations have been extensively employed to reveal the roles of protein dynamics in mediating enzyme catalysis. However, simulation-derived predictive descriptors that inform the impacts of mutations on catalytic turnover numbers remain largely unexplored. In this work, we report the identification of molecular modeling-derived descriptors to predict mutation effect on the turnover number of lactonase SsoPox with both native and non-native substrates. The study consists of 10 enzyme-substrate complexes resulting from a combination of five enzyme variants with two substrates. For each complex, we derived 15 descriptors from molecular dynamics simulations and applied principal component analysis to rank the predictive capability of the descriptors. A top-ranked descriptor was identified, which is the solvent-accessible surface area (SASA) ratio of the substrate to the active site pocket. A uniform volcano-shaped plot was observed in the distribution of experimental activation free energy against the SASA ratio. To achieve efficient lactonase hydrolysis, a non-native substrate-bound enzyme variant needs to involve a similar range of the SASA ratio to the native substrate-bound wild-type enzyme. The descriptor reflects how well the enzyme active site pocket accommodates a substrate for reaction, which has the potential of guiding optimization of enzyme reaction turnover for non-native chemical transformations.

Controlling Chemical Reactivity with Optimally Oriented Electric Fields: A Generalization of the Newton Trajectory Method

The use of oriented external electric fields (OEEF) as a tool to accelerate chemical reactions has recently attracted much interest. A new model to calculate the optimal OEEF of the least intensity to induce a barrierless chemical reaction path is presented. A suitable ansatz is provided by defining an effective potential energy surface (PES), which considers the unperturbed or original PES of the molecular reactive system and the action of a constant OEEF on the overall dipole moment of system. Based on a generalization of the Newton Trajectories (NT) method, it is demonstrated that the optimal OEEF can be determined upon locating a special point of the potential energy surface (PES), the so-called "optimal bond-breaking point" (optimal BBP), for which two different algorithms are proposed. At this point, the gradient of the original or unperturbed PES is an eigenvector of zero eigenvalue of the Hessian matrix of the effective PES. A thorough discussion of the geometrical aspects of the optimal BBP and the optimal OEEF is provided using a two-dimensional model, and numerical calculations of the optimal OEEF for a S_N2 reaction and the 1,3-dipolar retrocycloaddition of isoxazole to fulminic acid plus acetylene reaction serve as a proof of concept. The knowledge of the orientation of optimal OEEF provides a practical way to reduce the effective barrier of a given chemical process.

Linear Combination of Atomic Dipoles to Calculate the Bond and Molecular Dipole Moments of Molecules and Molecular Liquids

We report a linear combination of atomic dipole (LCAD) method for calculating the bond dipole moments of molecules. We show that the LCAD method reproduces the known molecular dipole moments of small to large molecules with a small error with respect to experimental and benchmark ab initio calculations, and molecular dipole distributions of bulk water that agree with maximally localized Wannier functions. The bond dipole moments derived from LCAD are also chemically interpretable in terms of the trend in bond ionicity in going from neutral to charged molecules. Moreover, the LCAD method accurately captures the influence of electric fields, supported by the correct trend in the change of the dipole moment under a uniform external electric field. The better grounding of bond dipole calculations indicates that it should also serve as a useful approach to bond dipole-field models used in catalysis or to reconstruct the small dipole of a H-terminated graphene flake.

The Interplay of Electrostatics and Chemical Positioning in the Evolution of Antibiotic Resistance in TEM β-Lactamases

The interplay of enzyme active site electrostatics and chemical positioning is important for understanding the origin(s) of enzyme catalysis and the design of novel catalysts. We reconstruct the evolutionary trajectory of TEM-1 β-lactamase to TEM-52 toward extended-spectrum activity to better understand the emergence of antibiotic resistance and to provide insights into the structure-function paradigm and noncovalent interactions involved in catalysis. Utilizing a detailed kinetic analysis and the vibrational Stark effect, we quantify the changes in rates and electric fields in the Michaelis and acyl-enzyme complexes for penicillin G and cefotaxime to ascertain the evolutionary role of electric fields to modulate function. These data are combined with MD simulations to interpret and quantify the substrate-dependent structural changes during evolution. We observe that this evolutionary trajectory utilizes a large preorganized electric field and substrate-dependent chemical positioning to facilitate catalysis. This governs the evolvability, substrate promiscuity, and protein fitness landscape in TEM β-lactamase antibiotic resistance.

An electrostatic duel: subtle differences in the catalytic performance of monoamine oxidase A and B isoenzymes elucidated at the residue level using quantum computations

The origin of the immense catalytic power of enzymes remains one of the biggest unresolved questions in biochemistry, with electrostatics being one of the main contenders. Herein, we report results that not only confirm that electrostatics is the driving force behind enzyme catalysis, but also that it is capable of tuning subtle differences in the catalytic performance between structurally similar enzymes, as demonstrated using the example of isoenzymes, monoamine oxidases A and B. Using our own computationally efficient multiscale model [A. Prah, et al., ACS Catal., 2019, 9, 1231] we analyzed the rate-limiting step of the reaction between phenylethylamine and both isoenzymes and deduced that the electrostatic environment provided by isoenzyme B has a perceivably higher catalytic influence on all the considered parameters of the reaction (energy barrier, charge transfer, dipole moment, and HOMO-LUMO gap). This is in full agreement with the available experimental kinetic data and with our own simulations of the reaction in question. In-depth analysis of individual amino acid contributions of both isoenzymes to the barrier (based on the interaction between the electric field provided by the enzyme and the dipole moment of the reacting moiety) shows that the majority of the difference between the isoenzymes can be attributed to a small number of sizable differences between the aligned amino acid pairs, whereas in most of the pairs the difference in contribution to the barrier is vanishingly small. These results suggest that electrostatics largely controls the substrate selectivity of enzymes and validates our approach as being capable of discerning fine nuances in the selectivity of structurally related isoenzymes.

Evolution of dynamical networks enhances catalysis in a designer enzyme

Activation heat capacity is emerging as a crucial factor in enzyme thermoadaptation, as shown by the non-Arrhenius behaviour of many natural enzymes. However, its physical origin and relationship to the evolution of catalytic activity remain uncertain. Here we show that directed evolution of a computationally designed Kemp eliminase reshapes protein dynamics, which gives rise to an activation heat capacity absent in the original design. These changes buttress transition-state stabilization. Extensive molecular dynamics simulations show that evolution results in the closure of solvent-exposed loops and a better packing of the active site. Remarkably, this gives rise to a correlated dynamical network that involves the transition state and large parts of the protein. This network tightens the transition-state ensemble, which induces a negative activation heat capacity and non-linearity in the activity-temperature dependence. Our results have implications for understanding enzyme evolution and suggest that selectively targeting the conformational dynamics of the transition-state ensemble by design and evolution will expedite the creation of novel enzymes.

Designing better enzymes: Insights from directed evolution

De novo enzymes can be created by computational design and directed evolution. Here, we review recent insights into the origins of catalytic power in evolved designer enzymes to pinpoint opportunities for next-generation designs: Evolution precisely organizes active sites, introduces catalytic H-bonding networks, invokes electrostatic catalysis, and creates dynamical networks embedding the active site in a reactive protein scaffold. Such insights foster our fundamental knowledge of enzyme catalysis and fuel the future design of tailor-made enzymes.

Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts

As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.

Catalytic Effect of Electric Fields on the Kemp Elimination Reactions with Neutral Bases

The effect of solvent reaction fields and oriented electric fields on the Kemp elimination reaction between methylamine or imidazole and 5-nitrobenzisoxazole has been theoretically studied. The Kemp reaction is the most widely used for the design of new enzymes. Our results, using the SMD continuous model for solvents, are in quite good agreement with the experimental fact that the rate of the analogous reaction with butylamine is one order of magnitude smaller in water than in acetonitrile. In the case of external electric fields, our results show that they can increase or decrease the energy barrier depending on the magnitude and orientation of the field. A duly oriented electric field may have a notable catalytic effect on the reaction. So, external electric fields and reaction fields due to the medium can contribute to the design of new enzymes. Several factors that must be taken into account to increase the catalytic effect are discussed.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

Bottom-Up Nonempirical Approach To Reducing Search Space in Enzyme Design Guided by Catalytic Fields

Currently developed protocols of theozyme design still lead to biocatalysts with much lower catalytic activity than enzymes existing in nature, and, so far, the only avenue of improvement was the in vitro laboratory-directed evolution (LDE) experiments. In this paper, we propose a different strategy based on "reversed" methodology of mutation prediction. Instead of common "top-down" approach, requiring numerous assumptions and vast computational effort, we argue for a "bottom-up" approach that is based on the catalytic fields derived directly from transition state and reactant complex wave functions. This enables direct one-step determination of the general quantitative angular characteristics of optimal catalytic site and simultaneously encompasses both the transition-state stabilization (TSS) and ground-state destabilization (GSD) effects. We further extend the static catalytic field approach by introducing a library of atomic multipoles for amino acid side-chain rotamers, which, together with the catalytic field, allow one to determine the optimal side-chain orientations of charged amino acids constituting the elusive structure of a preorganized catalytic environment. Obtained qualitative agreement with experimental LDE data for Kemp eliminase KE07 mutants validates the proposed procedure, yielding, in addition, a detailed insight into possible dynamic and epistatic effects.

Globally optimal catalytic fields for a Diels-Alder reaction

In a previous paper [M. Dittner and B. Hartke, J. Chem. Theory Comput. 14, 3547 (2018)], we introduced a preliminary version of our GOCAT (globally optimal catalyst) concept in which electrostatic catalysts are designed for arbitrary reactions by global optimization of distributed point charges that surround the reaction. In this first version, a pre-defined reaction path was kept fixed. This unrealistic assumption allowed for only small catalytic effects. In the present work, we extend our GOCAT framework by a sophisticated and robust on-the-fly reaction path optimization, plus further concomitant algorithm adaptions. This allows smaller and larger excursions from a pre-defined reaction path under the influence of the GOCAT point-charge surrounding, all the way to drastic mechanistic changes. In contrast to the restricted first GOCAT version, this new version is able to address real-life catalysis. We demonstrate this by applying it to the electrostatic catalysis of a prototypical Diels-Alder reaction. Without using any prior information, this procedure re-discovers theoretically and experimentally established features of electrostatic catalysis of this very reaction, including a field-dependent transition from the synchronous, concerted textbook mechanism to a zwitterionic two-step mechanism, and diastereomeric discrimination by suitable electric field components.

Interplay of water and a supramolecular capsule for catalysis of reductive elimination reaction from gold

Supramolecular assemblies have gained tremendous attention due to their ability to catalyze reactions with the efficiencies of natural enzymes. Using ab initio molecular dynamics, we identify the origin of the catalysis by the supramolecular capsule Ga₄L₆^12- on the reductive elimination reaction from gold complexes and assess their similarity to natural enzymes. By comparing the free energies of the reactants and transition states for the catalyzed and uncatalyzed reactions, we determine that an encapsulated water molecule generates electric fields that contributes the most to the reduction in the activation free energy. Although this is unlike the biomimetic scenario of catalysis through direct host-guest interactions, the electric fields from the nanocage also supports the transition state to complete the reductive elimination reaction with greater catalytic efficiency. However it is also shown that the nanocage poorly organizes the interfacial water, which in turn creates electric fields that misalign with the breaking bonds of the substrate, thus identifying new opportunities for catalytic design improvements in nanocage assemblies.

Kemp Elimination Reaction Catalyzed by Electric Fields

The Kemp elimination reaction is the most widely used in the de novo design of new enzymes. The effect of two different kinds of electric fields in the reactions of acetate as a base with benzisoxazole and 5-nitrobenzisoxazole as substrates have been theoretically studied. The effect of the solvent reaction field has been calculated using the SMD continuum model for several solvents; we have shown that solvents inhibit both reactions, the decrease of the reaction rate being larger as far as the dielectric constant is increased. The diminution of the reaction rate is especially remarkable between aprotic organic solvents and protic solvents as water, the electrostatic term of the hydrogen bonds being the main factor for the large inhibitory effect of water. The presence of an external electric field oriented in the direction of the charge transfer (z axis) increases it and, so, the reaction rate. In the reaction of the nitro compound, if the electric field is oriented in an orthogonal direction (x axis) the charge transfer to the NO₂ group is favored and there is a subsequent increase of the reaction rate. However, this increase is smaller than the one produced by the field in the z axis. It is worthwhile mentioning that one of the main effects of external electric fields of intermediate intensity is the reorientation of the reactants. Finally, the implications of our results in the de novo design of enzymes are discussed.

Quantum Mechanical Description of Electrostatics Provides a Unified Picture of Catalytic Action Across Methyltransferases

Methyl transferases (MTases) are a well-studied class of enzymes for which competing enzymatic enhancement mechanisms have been suggested, ranging from structural methyl group CH···X hydrogen bonds (HBs) to electrostatic- and charge-transfer-driven stabilization of the transition state (TS). We identified all Class I MTases for which reasonable resolution (<2.0 Å) crystal structures could be used to form catalytically competent ternary complexes for multiscale (i.e., quantum-mechanical/molecular-mechanical or QM/MM) simulation of the S_N2 methyl transfer reaction coordinate. The four Class I MTases studied have both distinct functions (e.g., protein repair or biosynthesis) and substrate nucleophiles (i.e., C, N, or O). While CH···X HBs stabilize all reactant complexes, no universal TS stabilization role is found for these interactions in MTases. A consistent picture is instead obtained through analysis of charge transfer and electrostatics, wherein much of cofactor-substrate charge separation is maintained in the TS region, and electrostatic potential is correlated with substrate nucleophilicity (i.e., intrinsic reactivity).

Globally Optimal Catalytic Fields - Inverse Design of Abstract Embeddings for Maximum Reaction Rate Acceleration

The search for, and understanding of, good catalysts for chemical reactions is a central issue for chemists. Here, we present first steps toward developing a general computational framework to better support this task. This framework combines efficient, unbiased global optimization techniques with an abstract representation of the catalytic environment, to shrink the search space. To analyze the resulting catalytic embeddings, we employ dimensionality reduction and clustering techniques. This not only provides an inverse design approach to new catalytic embeddings but also illuminates the actual interactions behind catalytic effects. All this is illustrated here with a strictly electrostatic model for the environment and with two versions of a selected example reaction. We close with detailed discussions of future improvements of our framework.