Design of a switchable eliminase

Enzyme Enhancement Through Computational Stability Design Targeting NMR-Determined Catalytic Hotspots

Enzymes are the quintessential green catalysts, but realizing their full potential for biotechnology typically requires improvement of their biomolecular properties. Catalysis enhancement, however, is often accompanied by impaired stability. Here, we show how the interplay between activity and stability in enzyme optimization can be efficiently addressed by coupling two recently proposed methodologies for guiding directed evolution. We first identify catalytic hotspots from chemical shift perturbations induced by transition-state-analogue binding and then use computational/phylogenetic design (FuncLib) to predict stabilizing combinations of mutations at sets of such hotspots. We test this approach on a previously designed de novo Kemp eliminase, which is already highly optimized in terms of both activity and stability. Most tested variants displayed substantially increased denaturation temperatures and purification yields. Notably, our most efficient engineered variant shows a ∼3-fold enhancement in activity (kcat ∼ 1700 s-1, kcat/KM ∼ 4.3 × 105 M-1 s-1) from an already heavily optimized starting variant, resulting in the most proficient proton-abstraction Kemp eliminase designed to date, with a catalytic efficiency on a par with naturally occurring enzymes. Molecular simulations pinpoint the origin of this catalytic enhancement as being due to the progressive elimination of a catalytically inefficient substrate conformation that is present in the original design. Remarkably, interaction network analysis identifies a significant fraction of catalytic hotspots, thus providing a computational tool which we show to be useful even for natural-enzyme engineering. Overall, our work showcases the power of dynamically guided enzyme engineering as a design principle for obtaining novel biocatalysts with tailored physicochemical properties, toward even anthropogenic reactions.

Emergence of specific binding and catalysis from a designed generalist binding protein

The evolution of binding and catalysis played a central role in the emergence of life. While natural proteins have finely tuned affinities for their primary ligands, they also bind weakly and promiscuously to other molecules, which serve as starting points for stepwise, incremental evolution of entirely new specificities. Thus, modern proteins emerged from the joint exploration of sequence and structural space. The ability of natural proteins to bind promiscuously to small molecule fragments has been widely evaluated using methods including crystallographic fragment screening. However, this approach had not been applied to de novo proteins. Here, we apply this method to explore the promiscuity of a de novo small molecule-binding protein ABLE. As in Nature, we found ABLE was capable of forming weak complexes, which were found to be excellent starting points for evolving entirely new functions, including a binder of a turn-on fluorophore and a highly efficient and specific Kemp eliminase enzyme. This work shows how Nature and protein designers can take advantage of promiscuous binding interactions to evolve new proteins with specialized functions.

Evolution Enhances Kemp Eliminase Activity by Optimizing Oxyanion Stabilization and Conformational Flexibility

The base-promoted Kemp elimination reaction has been used as a model system for enzyme design. Among the multiple computationally designed and evolved Kemp eliminases generated along the years, the HG3-to-HG3.17 evolutionary trajectory is particularly interesting due to the high catalytic efficiency of HG3.17 and the debated role of glutamine 50 (Gln50) as potential oxyanion stabilizer. This study aims to elucidate the structural and dynamic changes along the evolutionary pathway from HG3 to HG3.17 that contribute to improved catalytic efficiency. In particular, we evaluate key variants along the HG3 evolutionary trajectory via molecular dynamics simulations coupled to non-covalent interactions and water analysis. Our computational study indicates that HG3.17 can adopt a catalytically competent conformation promoted by a water-mediated network of non-covalent interactions, in which aspartate 127 (Asp127) is properly positioned for proton abstraction and Gln50 and to some extent mutation cysteine 84 (Cys84) contribute to oxyanion stabilization. We find that HG3.17 exhibits a rather high flexibility of Gln50, which is regulated by the conformation adopted by the active site residue tryptophan 44 (Trp44). This interplay between Gln50 and Trp44 positioning induced by distal active site mutations affects the water-mediated network of non-covalent interactions, Gln50 preorganization, and water content of the active site pocket.

Fishing for Catalysis: Experimental Approaches to Narrowing Search Space in Directed Evolution of Enzymes

Directed evolution has transformed protein engineering offering a path to rapid improvement of protein properties. Yet, in practice it is limited by the hyper-astronomic protein sequence search space, and approaches to identify mutagenic hot spots, i.e., locations where mutations are most likely to have a productive impact, are needed. In this perspective, we categorize and discuss recent progress in the experimental approaches (broadly defined as structural, bioinformatic, and dynamic) to hot spot identification. Recent successes in harnessing protein dynamics and machine learning approaches provide new opportunities for the field and will undoubtedly help directed evolution reach its full potential.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

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.

NMR-guided directed evolution

Directed evolution is a powerful tool for improving existing properties and imparting completely new functionalities to proteins1-4. Nonetheless, its potential in even small proteins is inherently limited by the astronomical number of possible amino acid sequences. Sampling the complete sequence space of a 100-residue protein would require testing of 20100 combinations, which is beyond any existing experimental approach. In practice, selective modification of relatively few residues is sufficient for efficient improvement, functional enhancement and repurposing of existing proteins5. Moreover, computational methods have been developed to predict the locations and, in certain cases, identities of potentially productive mutations6-9. Importantly, all current approaches for prediction of hot spots and productive mutations rely heavily on structural information and/or bioinformatics, which is not always available for proteins of interest. Moreover, they offer a limited ability to identify beneficial mutations far from the active site, even though such changes may markedly improve the catalytic properties of an enzyme10. Machine learning methods have recently showed promise in predicting productive mutations11, but they frequently require large, high-quality training datasets, which are difficult to obtain in directed evolution experiments. Here we show that mutagenic hot spots in enzymes can be identified using NMR spectroscopy. In a proof-of-concept study, we converted myoglobin, a non-enzymatic oxygen storage protein, into a highly efficient Kemp eliminase using only three mutations. The observed levels of catalytic efficiency exceed those of proteins designed using current approaches and are similar with those of natural enzymes for the reactions that they are evolved to catalyse. Given the simplicity of this experimental approach, which requires no a priori structural or bioinformatic knowledge, we expect it to be widely applicable and to enable the full potential of directed enzyme evolution.

Active site center redesign increases protein stability preserving catalysis in thioredoxin

The stabilization of natural proteins is a long-standing desired goal in protein engineering. Optimizing the hydrophobicity of the protein core often results in extensive stability enhancements. However, the presence of totally or partially buried catalytic charged residues, essential for protein function, has limited the applicability of this strategy. Here, focusing on the thioredoxin, we aimed to augment protein stability by removing buried charged residues in the active site without loss of catalytic activity. To this end, we performed a charged-to-hydrophobic substitution of a buried and functional group, resulting in a significant stability increase yet abolishing catalytic activity. Then, to simulate the catalytic role of the buried ionizable group, we designed a combinatorial library of variants targeting a set of seven surface residues adjacent to the active site. Notably, more than 50% of the library variants restored, to some extent, the catalytic activity. The combination of experimental study of 2% of the library with the prediction of the whole mutational space by partial least squares regression revealed that a single point mutation at the protein surface is sufficient to fully restore the catalytic activity without thermostability cost. As a result, we engineered one of the highest thermal stabilities reported for a protein with a natural occurring fold (137°C). Further, our hyperstable variant preserves the catalytic activity both in vitro and in vivo.

Phosphorylation promotes the endonuclease-like activity of human centrin 2

Centrin is a member of the EF-hand superfamily of calcium-binding proteins, which is involved in the nucleotide excision repair (NER). Reversible phosphorylation of centrin is an important regulatory mechanism in vivo and is closely related to many physiological processes. To explore the possible role of centrin in NER, the endonuclease-like activity of human centrin 2 (HsCen2) regulated by phosphorylation in the absence or presence of Tb³⁺ was investigated by spectroscopy techniques, gel electrophoresis, and molecular docking simulation in 10 mM Hepes, pH 7.4. The results showed that phosphorylation weakened the binding of Tb³⁺ to HsCen2 and enhanced the binding of DNA to HsCen2. Phosphorylation improves the endonuclease-like activity of HsCen2. In addition, Tb³⁺ is favorable for DNA binding and endonuclease-like activity of HsCen2 before and after phosphorylation. These results provide clear insights into the effects of phosphorylation on the properties of HsCen2 and offer important clues for further exploration of how phosphorylation affects protein-driven functions.

Phosphorylation weakened the binding of Tb³⁺ to HsCen2, enhanced the binding of DNA to HsCen2; and improves the endonuclease-like activity of HsCen2; Additionally, the endonuclease-like activity of HsCen2 or HsCen2p is regulated up by Tb³⁺-binding.

Tuning the Catalytic Activity of Synthetic Enzyme KE15 with DNA

Efficiency improvement of synthetic enzymes through scaffold modifications suffers from limitations in terms of effectiveness, cost, and potential devastating consequences for protein structural stability. Here, we propose an alternative to scaffold modification, within electrostatic preorganization theory, where the enzyme's greater environment is designed to support the evolution of the reaction in the active site. We demonstrate the feasibility of such an approach by placing a (polar) DNA fragment in the surroundings of the Kemp eliminase enzyme KE15 (structure from Houk's group) and computing the resulting change in catalytic activity. We find that the introduction of a DNA fragment magnifies the contribution of protein residues to the stabilization of the transition state, estimated from electric field calculations with polarizable molecular dynamics. Our randomly generated test systems reveal a 2.0 kcal/mol reduction in activation energy, suggesting that even more significant catalytic improvements could be made by optimizing DNA size, sequence, and orientation with respect to the enzyme, validating our approach.

Rational design of allosteric switchable catalysts

Allosteric regulation, in many cases, involves switching the activities of natural enzymes, which further affects the enzymatic network and cell signaling in the living systems. The research on the construction of allosteric switchable catalysts has attracted broad interests, aiming to control the progress and asymmetry of catalytic reactions, expand the chemical biology toolbox, substitute unstable natural enzymes in the biological detection and biosensors, and fabricate the biomimetic cascade reactions. Thus, in this review, we summarize the recent outstanding works in switchable catalysts based on the allosterism of single molecules, supramolecular complexes, and self-assemblies. The concept of allosterism was extended from natural proteins to polymers, organic molecules, and supramolecular systems. In terms of the difference between these building scaffolds, a variety of design methods that tailor biological and synthetic molecules into controllable catalysts were introduced with emphasis.

Structural characterization and Kemp eliminase activity of the Mycobacterium smegmatis Ketosteroid Isomerase

The Kemp elimination reaction, involving the ring-opening of benzoxazole and its derivatives under the action of natural enzymes or chemical catalysts, has been of interest to researchers since its discovery. Because this reaction does not exist in all currently known metabolic pathways, the computational design of Kemp eliminases has provided valuable insights into principles of enzymatic catalysis. However, it was discovered that the naturally occurring promiscuous enzymes ydbC, xapA and ketosteroid isomerase also can catalyze Kemp elimination. Here, we report the crystal structure of ketosteroid isomerase (KSI) from Mycobacterium smegmatis MC2 155. MsKSI crystallizes in the P2₁2₁2₁ space group with two molecules in an asymmetric unit, and ultracentrifugation data confirms that it forms a stable dimer in solution, consistent with the 1.9 Å-resolution structure. Our assays confirm that MsKSI accelerates the Kemp elimination of 5-nitrobenzoxazole (5NBI) with an optimal pH of 5.5. A 2.35 Å resolution crystal structure of the MsKSI-5NBI complex reveals that the substrate 5NBI is bound in the active pocket of the enzyme composed of hydrophobic residues. In addition, the Glu127 residue is proposed to play an important role as a general base in proton transfer and breaking weak O-N bonds to open the five-membered ring. This work provides a starting point for exploring the artificial modification of MsKSI using the natural enzyme as the backbone.

How directed evolution reshapes the energy landscape in an enzyme to boost catalysis

The advent of biocatalysts designed computationally and optimized by laboratory evolution provides an opportunity to explore molecular strategies for augmenting catalytic function. Applying a suite of nuclear magnetic resonance, crystallography, and stopped-flow techniques to an enzyme designed for an elementary proton transfer reaction, we show how directed evolution gradually altered the conformational ensemble of the protein scaffold to populate a narrow, highly active conformational ensemble and accelerate this transformation by nearly nine orders of magnitude. Mutations acquired during optimization enabled global conformational changes, including high-energy backbone rearrangements, that cooperatively organized the catalytic base and oxyanion stabilizer, thus perfecting transition-state stabilization. The development of protein catalysts for many chemical transformations could be facilitated by explicitly sampling conformational substates during design and specifically stabilizing productive substates over all unproductive conformations.

Superior Proton-Transfer Catalytic Promiscuity of Cytochrome c in Self-Organized Media

Evolutionarily elderly proteins commonly feature greater catalytic promiscuity. Cytochrome c is among the first set of proteins in evolution to have known prospects in electron transport and peroxidative properties. Here, we report that cyt c is also a proficient proton-transfer catalyst and enhances the Kemp elimination (KE; model reaction to show proton transfer catalytic property) by ∼750-fold on self-organized systems like micelles and vesicles. The self-organized systems mimic the mitochondrial environment in vitro for cyt c. Using an array of biophysical and biochemical mutational assays, both acid-base and redox mechanistic pathways have been explored. The histidine moiety close to hemin group (His18) is mainly responsible for proton abstraction to promote the concerted E2 pathway for KE catalysis when cyt c is in its oxidized form; this has also been confirmed by a H18A mutant of cyt c. However, the redox pathway is predominant under reducing conditions in the presence of dithiothreitol over the pH range 6-7.4. Interestingly, we found almost 750-fold enhanced KE catalysis by cyt c compared to aqueous buffer. Overall, in addition to providing mechanistic insights, the data reveal an unprecedented catalytic property of cyt c that could be of high importance in an evolutionary perspective considering its role in delineating the phylogenic tree and also towards generating programmable designer biocatalysts.

Enhancing a de novo enzyme activity by computationally-focused ultra-low-throughput screening

Directed evolution has revolutionized protein engineering. Still, enzyme optimization by random library screening remains sluggish, in large part due to futile probing of mutations that are catalytically neutral and/or impair stability and folding. FuncLib is a novel approach which uses phylogenetic analysis and Rosetta design to rank enzyme variants with multiple mutations, on the basis of predicted stability. Here, we use it to target the active site region of a minimalist-designed, de novo Kemp eliminase. The similarity between the Michaelis complex and transition state for the enzymatic reaction makes this system particularly challenging to optimize. Yet, experimental screening of a small number of active-site variants at the top of the predicted stability ranking leads to catalytic efficiencies and turnover numbers (∼2 × 104 M-1 s-1 and ∼102 s-1) for this anthropogenic reaction that compare favorably to those of modern natural enzymes. This result illustrates the promise of FuncLib as a powerful tool with which to speed up directed evolution, even on scaffolds that were not originally evolved for those functions, by guiding screening to regions of the sequence space that encode stable and catalytically diverse enzymes. Empirical valence bond calculations reproduce the experimental activation energies for the optimized eliminases to within ∼2 kcal mol-1 and indicate that the enhanced activity is linked to better geometric preorganization of the active site. This raises the possibility of further enhancing the stability-guidance of FuncLib by computational predictions of catalytic activity, as a generalized approach for computational enzyme design.

Rational molecular design for improving digestive enzyme resistance of beta-glucosidase from Trichoderma viride based on inhibition of bound state formation

Beta-glucosidase (BGL1) is widely used in animal feed industries. However, degradation caused by digestive enzymes in the intestine hampers its application. Improving the resistance of feed enzymes against proteases is crucial in livestock farming. To improve the resistance of beta-glucosidase against pepsin and trypsin, a rational molecular design based on the inhibition of bound-state formation and secondary design was developed. The strategy includes: (1) prediction of the interaction surface of the pepsin-BGL1 complex structure, (2) prediction of key amino acids affecting the formation of the complex, (3) optimization of pepsin-resistant mutants by structural evaluation, (4) secondary molecular design based on pepsin-resistant mutants, and optimization of pepsin and trypsin-resistant mutants. Two BGL1 protein mutants (BGL1^Q627C and BGL1^{Q627C/R543H/R646W}) were constructed, and then mutated and wild-type BGL1s were expressed in Pichia pastoris. The half-life of BGL1^Q627C and BGL1^{Q627C/R543H/R646W} were 1.36 and 1.51 times that of the wild type upon pepsin exposure, respectively. For trypsin resistance, the half-life were 0.93 and 1.53 times that of the wild type, respectively. Compare to those of the wild type, most of the basic enzymatic properties of both mutants were not significantly changed except for increased Michaelis constants. The rational design method can be used as a guide for modifying other feed enzymes.

Minimalist de novo Design of Protein Catalysts

The field of protein design has grown enormously in the past few decades. In this review we discuss the minimalist approach to design of artificial enzymes, in which protein sequences are created with the minimum number of elements for folding and function. This method relies on identifying starting points in catalytically inert scaffolds for active site installation. The progress of the field from the original helical assemblies of the 1980s to the more complex structures of the present day is discussed, highlighting the variety of catalytic reactions which have been achieved using these methods. We outline the strengths and weaknesses of the minimalist approaches, describe representative design cases and put it in the general context of the de novo design of proteins.

Kemp Eliminases of the AlleyCat Family Possess High Substrate Promiscuity

Minimalist enzymes designed to catalyze model reactions provide useful starting points for creating catalysts for practically important chemical transformations. We have shown that Kemp eliminases of the AlleyCat family facilitate conversion of leflunomide (an immunosupressor pro-drug) to its active form teriflunomide with outstanding rate enhancement (nearly four orders of magnitude) and catalytic proficiency (more than seven orders of magnitude) without any additional optimization. This remarkable activity is achieved by properly positioning the substrate in close proximity to the catalytic glutamate with very high pKa.

Design of catalytic polypeptides and proteins

The first part of this review article lists examples of complete, empirical de novo design that made important contributions to the development of the field and initiated challenging projects. The second part of this article deals with computational design of novel enzymes in native protein scaffolds; active designs were refined through random and site-directed mutagenesis producing artificial enzymes with nearly native enzyme- like activities against a number of non-natural substrates. Combining aspects of de novo design and biological evolution of nature's enzymes has started and will accelerate the development of novel enzyme activities.

Strategies for designing non-natural enzymes and binders

The design of tailor-made enzymes is a major goal in biochemical research that can result in wide-range applications and will lead to a better understanding of how proteins fold and function. In this review we highlight recent advances in enzyme and small molecule binder design. A focus is placed on novel strategies for the design of scaffolds, developments in computational methods, and recent applications of these techniques on receptors, sensors, and enzymes. Further, the integration of computational and experimental methodologies is discussed. The outlined examples of designed enzymes and binders for various purposes highlight the importance of this topic and underline the need for tailor-made proteins.

A Designed Enzyme Promotes Selective Post-translational Acylation

A computationally designed, allosterically regulated catalyst (CaM M144H) produced by substituting a single residue in calmodulin, a non-enzymatic protein, is capable of efficient and site selective post-translational acylation of lysines in peptides with highly diverse sequences. Calmodulin's binding partners are involved in regulating a large number of cellular processes; this new chemical-biology tool will help to identify them and provide structural insight into their interactions with calmodulin.

Protein Environment: A Crucial Triggering Factor in Josephin Domain Aggregation: The Role of 2,2,2-Trifluoroethanol

The protein ataxin-3 contains a polyglutamine stretch that triggers amyloid aggregation when it is expanded beyond a critical threshold. This results in the onset of the spinocerebellar ataxia type 3. The protein consists of the globular N-terminal Josephin domain and a disordered C-terminal tail where the polyglutamine stretch is located. Expanded ataxin-3 aggregates via a two-stage mechanism: first, Josephin domain self-association, then polyQ fibrillation. This highlights the intrinsic amyloidogenic potential of Josephin domain. Therefore, much effort has been put into investigating its aggregation mechanism(s). A key issue regards the conformational requirements for triggering amyloid aggregation, as it is believed that, generally, misfolding should precede aggregation. Here, we have assayed the effect of 2,2,2-trifluoroethanol, a co-solvent capable of stabilizing secondary structures, especially α-helices. By combining biophysical methods and molecular dynamics, we demonstrated that both secondary and tertiary JD structures are virtually unchanged in the presence of up to 5% 2,2,2-trifluoroethanol. Despite the preservation of JD structure, 1% of 2,2,2-trifluoroethanol suffices to exacerbate the intrinsic aggregation propensity of this domain, by slightly decreasing its conformational stability. These results indicate that in the case of JD, conformational fluctuations might suffice to promote a transition towards an aggregated state without the need for extensive unfolding, and highlights the important role played by the environment on the aggregation of this globular domain.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Role of conformational dynamics in the evolution of novel enzyme function

The free energy landscape concept that describes enzymes as an ensemble of differently populated conformational sub-states in dynamic equilibrium is key for evaluating enzyme activity, enantioselectivity, and specificity. Mutations introduced in the enzyme sequence can alter the populations of the pre-existing conformational states, thus strongly modifying the enzyme ability to accommodate alternative substrates, revert its enantiopreferences, and even increase the activity for some residual promiscuous reactions. In this feature article, we present an overview of the current experimental and computational strategies to explore the conformational free energy landscape of enzymes. We provide a series of recent publications that highlight the key role of conformational dynamics for the enzyme evolution towards new functions and substrates, and provide some perspectives on how conformational dynamism should be considered in future computational enzyme design protocols.

Rational and Semirational Protein Design

This mini review gives an overview over different design approaches and methodologies applied in rational and semirational enzyme engineering. The underlying principles for engineering novel activities, enantioselectivity, substrate specificity, stability, and pH optimum are summarized.

Computational Design of Protein Function

The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.

De novo active sites for resurrected Precambrian enzymes

Protein engineering studies often suggest the emergence of completely new enzyme functionalities to be highly improbable. However, enzymes likely catalysed many different reactions already in the last universal common ancestor. Mechanisms for the emergence of completely new active sites must therefore either plausibly exist or at least have existed at the primordial protein stage. Here, we use resurrected Precambrian proteins as scaffolds for protein engineering and demonstrate that a new active site can be generated through a single hydrophobic-to-ionizable amino acid replacement that generates a partially buried group with perturbed physico-chemical properties. We provide experimental and computational evidence that conformational flexibility can assist the emergence and subsequent evolution of new active sites by improving substrate and transition-state binding, through the sampling of many potentially productive conformations. Our results suggest a mechanism for the emergence of primordial enzymes and highlight the potential of ancestral reconstruction as a tool for protein engineering.

Artificial Metalloenzymes: Reaction Scope and Optimization Strategies

The incorporation of a synthetic, catalytically competent metallocofactor into a protein scaffold to generate an artificial metalloenzyme (ArM) has been explored since the late 1970's. Progress in the ensuing years was limited by the tools available for both organometallic synthesis and protein engineering. Advances in both of these areas, combined with increased appreciation of the potential benefits of combining attractive features of both homogeneous catalysis and enzymatic catalysis, led to a resurgence of interest in ArMs starting in the early 2000's. Perhaps the most intriguing of potential ArM properties is their ability to endow homogeneous catalysts with a genetic memory. Indeed, incorporating a homogeneous catalyst into a genetically encoded scaffold offers the opportunity to improve ArM performance by directed evolution. This capability could, in turn, lead to improvements in ArM efficiency similar to those obtained for natural enzymes, providing systems suitable for practical applications and greater insight into the role of second coordination sphere interactions in organometallic catalysis. Since its renaissance in the early 2000's, different aspects of artificial metalloenzymes have been extensively reviewed and highlighted. Our intent is to provide a comprehensive overview of all work in the field up to December 2016, organized according to reaction class. Because of the wide range of non-natural reactions catalyzed by ArMs, this was done using a functional-group transformation classification. The review begins with a summary of the proteins and the anchoring strategies used to date for the creation of ArMs, followed by a historical perspective. Then follows a summary of the reactions catalyzed by ArMs and a concluding critical outlook. This analysis allows for comparison of similar reactions catalyzed by ArMs constructed using different metallocofactor anchoring strategies, cofactors, protein scaffolds, and mutagenesis strategies. These data will be used to construct a searchable Web site on ArMs that will be updated regularly by the authors.

Control of Heme Coordination and Catalytic Activity by Conformational Changes in Peptide-Amphiphile Assemblies

Self-assembling peptide materials have gained significant attention, due to well-demonstrated applications, but they are functionally underutilized. To advance their utility, we use noncovalent interactions to incorporate the biological cofactor heme-B for catalysis. Heme-proteins achieve differing functions through structural and coordinative variations. Here, we replicate this phenomenon by highlighting changes in heme reactivity as a function of coordination, sequence, and morphology (micelles versus fibers) in a series of simple peptide amphiphiles with the sequence c16-xyL₃K₃-CO₂H where c16 is a palmitoyl moiety and xy represents the heme binding region: AA, AH, HH, and MH. The morphology of this peptide series is characterized using transmission electron and atomic force microscopies as well as dynamic light scattering. Within this small library of peptide constructs, we show that three spectroscopically (UV/visible and electron paramagnetic resonance) distinct heme environments were generated: noncoordinated/embedded high-spin, five-coordinate high-spin, and six-coordinate low-spin. The resulting material's functional dependence on sequence and supramolecular morphology is highlighted 2-fold. First, the heme active site binds carbon monoxide in both micelles and fibers, demonstrating that the heme active site in both morphologies is accessible to small molecules for catalysis. Second, peroxidase activity was observed in heme-containing micelles yet was significantly reduced in heme-containing fibers. We briefly discuss the implications these findings have in the production of functional, self-assembling peptide materials.

Exploring the Development of Ground-State Destabilization and Transition-State Stabilization in Two Directed Evolution Paths of Kemp Eliminases

Computer-aided enzyme design presents a major challenge since in most cases it has not resulted in an impressive catalytic power. The reasons for the problems with computational design include the use of nonquantitative approaches, but they may also reflect other difficulties that are not completely obvious. Thus, it is very useful to try to learn from the trend in directed evolution experiments. Here we explore the nature of the refinement of Kemp eliminases by directed evolution, trying to gain an understanding of related requirements from computational design. The observed trend in the directed evolution refinement of KE07 and HG3 are reproduced, showing that in the case of KE07 the directed evolution leads to ground-state destabilization, whereas in the case of HG3 the directed evolution leads to transition-state stabilization. The nature of the different paths of the directed evolution is examined and discussed. The present study seems to indicate that computer-aided enzyme design may require more than calculations of the effect of single mutations and should be extended to calculations of the effect of simultaneous multiple mutations (that make a few residues preorganized effectively). However, the analysis of two known evolution paths can still be accomplished using the relevant sequences and structures. Thus, by comparing two directed evolution paths of Kemp eliminases we reached the important conclusion that the more effective path leads to transition-state stabilization.

Insights Into How Heme Reduction Potentials Modulate Enzymatic Activities of a Myoglobin-based Functional Oxidase

Heme-copper oxidase (HCO) is a class of respiratory enzymes that use a heme-copper center to catalyze O₂ reduction to H₂ O. While heme reduction potential (E°') of different HCO types has been found to vary >500 mV, its impact on HCO activity remains poorly understood. Here, we use a set of myoglobin-based functional HCO models to investigate the mechanism by which heme E°' modulates oxidase activity. Rapid stopped-flow kinetic measurements show that increasing heme E°' by ca. 210 mV results in increases in electron transfer (ET) rates by 30-fold, rate of O₂ binding by 12-fold, O₂ dissociation by 35-fold, while decreasing O₂ affinity by 3-fold. Theoretical calculations reveal that E°' modulation has significant implications on electronic charge of both heme iron and O₂ , resulting in increased O₂ dissociation and reduced O₂ affinity at high E°' values. Overall, this work suggests that fine-tuning E°' in HCOs and other heme enzymes can modulate their substrate affinity, ET rate and enzymatic activity.

A redox-mediated Kemp eliminase

The acid/base-catalysed Kemp elimination of 5-nitro-benzisoxazole forming 2-cyano-4-nitrophenol has long served as a design platform of enzymes with non-natural reactions, providing new mechanistic insights in protein science. Here we describe an alternative concept based on redox catalysis by P450-BM3, leading to the same Kemp product via a fundamentally different mechanism. QM/MM computations show that it involves coordination of the substrate's N-atom to haem-Fe(II) with electron transfer and concomitant N-O heterolysis liberating an intermediate having a nitrogen radical moiety Fe(III)-N· and a phenoxyl anion. Product formation occurs by bond rotation and H-transfer. Two rationally chosen point mutations cause a notable increase in activity. The results shed light on the prevailing mechanistic uncertainties in human P450-catalysed metabolism of the immunomodulatory drug leflunomide, which likewise undergoes redox-mediated Kemp elimination by P450-BM3. Other isoxazole-based pharmaceuticals are probably also metabolized by a redox mechanism. Our work provides a basis for designing future artificial enzymes.

Incorporating an allosteric regulatory site in an antibody through backbone design

Allosteric regulation underlies living cells' ability to sense changes in nutrient and signaling-molecule concentrations, but the ability to computationally design allosteric regulation into non-allosteric proteins has been elusive. Allosteric-site design is complicated by the requirement to encode the relative stabilities of active and inactive conformations of the same protein in the presence and absence of both ligand and effector. To address this challenge, we used Rosetta to design the backbone of the flexible heavy-chain complementarity-determining region 3 (HCDR3), and used geometric matching and sequence optimization to place a Zn²⁺ -coordination site in a fluorescein-binding antibody. We predicted that due to HCDR3's flexibility, the fluorescein-binding pocket would configure properly only upon Zn²⁺ application. We found that regulation by Zn²⁺ was reversible and sensitive to the divalent ion's identity, and came at the cost of reduced antibody stability and fluorescein-binding affinity. Fluorescein bound at an order of magnitude higher affinity in the presence of Zn²⁺ than in its absence, and the increase in fluorescein affinity was due almost entirely to faster fluorescein on-rate, suggesting that Zn²⁺ preorganized the antibody for fluorescein binding. Mutation analysis demonstrated the extreme sensitivity of Zn²⁺ regulation on the atomic details in and around the metal-coordination site. The designed antibody could serve to study how allosteric regulation evolved from non-allosteric binding proteins, and suggests a way to designing molecular sensors for environmental and biomedical targets.

Why copper is preferred over iron for oxygen activation and reduction in haem-copper oxidases

Haem-copper oxidase (HCO) catalyses the natural reduction of oxygen to water using a haem-copper centre. Despite decades of research on HCOs, the role of non-haem metal and the reason for nature's choice of copper over other metals such as iron remains unclear. Here, we use a biosynthetic model of HCO in myoglobin that selectively binds different non-haem metals to demonstrate 30-fold and 11-fold enhancements in the oxidase activity of Cu- and Fe-bound HCO mimics, respectively, as compared with Zn-bound mimics. Detailed electrochemical, kinetic and vibrational spectroscopic studies, in tandem with theoretical density functional theory calculations, demonstrate that the non-haem metal not only donates electrons to oxygen but also activates it for efficient O-O bond cleavage. Furthermore, the higher redox potential of copper and the enhanced weakening of the O-O bond from the higher electron density in the d orbital of copper are central to its higher oxidase activity over iron. This work resolves a long-standing question in bioenergetics, and renders a chemical-biological basis for the design of future oxygen-reduction catalysts.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Kemp Eliminase Activity of Ketosteroid Isomerase

Kemp eliminases represent the most successful class of computationally designed enzymes, with rate accelerations of up to 109-fold relative to the rate of the same reaction in aqueous solution. Nevertheless, several other systems such as micelles, catalytic antibodies, and cavitands are known to accelerate the Kemp elimination by several orders of magnitude. We found that the naturally occurring enzyme ketosteroid isomerase (KSI) also catalyzes the Kemp elimination. Surprisingly, mutations of D38, the residue that acts as a general base for its natural substrate, produced variants that catalyze the Kemp elimination up to 7000-fold better than wild-type KSI does, and some of these variants accelerate the Kemp elimination more than the computationally designed Kemp eliminases. Analysis of the D38N general base KSI variant suggests that a different active site carboxylate residue, D99, performs the proton abstraction. Docking simulations and analysis of inhibition by active site binders suggest that the Kemp elimination takes place in the active site of KSI and that KSI uses the same catalytic strategies of the computationally designed enzymes. In agreement with prior observations, our results strengthen the conclusion that significant rate accelerations of the Kemp elimination can be achieved with very few, nonspecific interactions with the substrate if a suitable catalytic base is present in a hydrophobic environment. Computational design can fulfill these requirements, and the design of more complex and precise environments represents the next level of challenges for protein design.