Protein engineering: the potential of remote mutations

Structural Evolution of Bacterial Polyphosphate Degradation Enzyme for Phosphorus Cycling

Living organisms ranging from bacteria to animals have developed their own ways to accumulate and store phosphate during evolution, in particular as the polyphosphate (polyP) granules in bacteria. Degradation of polyP into phosphate is involved in phosphorus cycling, and exopolyphosphatase (PPX) is the key enzyme for polyP degradation in bacteria. Thus, understanding the structure basis of PPX is crucial to reveal the polyP degradation mechanism. Here, it is found that PPX structure varies in the length of ɑ-helical interdomain linker (ɑ-linker) across various bacteria, which is negatively correlated with their enzymatic activity and thermostability - those with shorter ɑ-linkers demonstrate higher polyP degradation ability. Moreover, the artificial DrPPX mutants with shorter ɑ-linker tend to have more compact pockets for polyP binding and stronger subunit interactions, as well as higher enzymatic efficiency (kcat/Km) than that of DrPPX wild type. In Deinococcus-Thermus, the PPXs from thermophilic species possess a shorter ɑ-linker and retain their catalytic ability at high temperatures (70 °C), which may facilitate the thermophilic species to utilize polyP in high-temperature environments. These findings provide insights into the interdomain linker length-dependent evolution of PPXs, which shed light on enzymatic adaption for phosphorus cycling during natural evolution and rational design of enzyme.

Perspectives on improving crop Rubisco by directed evolution

Rubisco catalyses the entry of almost all CO2 into the biosphere and is often the rate-limiting step in plant photosynthesis and growth. Its notoriety as the most abundant protein on Earth stems from the slow and error-prone catalytic properties that require plants, cyanobacteria, algae and photosynthetic bacteria to produce it in high amounts. Efforts to improve the CO2-fixing properties of plant Rubisco has been spurred on by the discovery of more effective isoforms in some algae with the potential to significantly improve crop productivity. Incompatibilities between the protein folding machinery of leaf and algae chloroplasts have, so far, prevented efforts to transplant these more effective Rubisco variants into plants. There is therefore increasing interest in improving Rubisco catalysis by directed (laboratory) evolution. Here we review the advances being made in, and the ongoing challenges with, improving the solubility and/or carboxylation activity of differing non-plant Rubisco lineages. We provide perspectives on new opportunities for the directed evolution of crop Rubiscos and the existing plant transformation capabilities available to evaluate the extent to which Rubisco activity improvements can benefit agricultural productivity.

Improving plastic degrading enzymes via directed evolution

Plastic degrading enzymes have immense potential for use in industrial applications. Protein engineering efforts over the last decade have resulted in considerable enhancement of many properties of these enzymes. Directed evolution, a protein engineering approach that mimics the natural process of evolution in a laboratory, has been particularly useful in overcoming some of the challenges of structure-based protein engineering. For example, directed evolution has been used to improve the catalytic activity and thermostability of polyethylene terephthalate (PET)-degrading enzymes, although its use for the improvement of other desirable properties, such as solvent tolerance, has been less studied. In this review, we aim to identify some of the knowledge gaps and current challenges, and highlight recent studies related to the directed evolution of plastic-degrading enzymes.

Rational peptide design for targeting cancer cell invasion

Cell invasion is an important process in cancer progression and recurrence. Invasion and implantation of cancer cells from their original place to other tissues, by disabling vital organs, challenges the treatment of cancer patients. Given the importance of the matter, many molecular treatments have been developed to inhibit cancer cell invasion. Because of their low production cost and ease of production, peptides are valuable therapeutic molecules for inhibiting cancer cell invasion. In recent years, advances in the field of computational biology have facilitated the design of anti-cancer peptides. In our investigation, using computational biology approaches such as evolutionary analysis, residue scanning, protein-peptide interaction analysis, molecular dynamics, and free energy analysis, our team designed a peptide library with about 100 000 candidates based on A6 (acetyl-KPSSPPEE-amino) sequence which is an anti-invasion peptide. During computational studies, two of the designed peptides that give the highest scores and showed the greatest sequence similarity to A6 were entered into the experimental analysis workflow for further analysis. In experimental analysis steps, the anti-metastatic potency and other therapeutic effects of designed peptides were evaluated using MTT assay, RT-qPCR, zymography analysis, and invasion assay. Our study disclosed that the IK1 (acetyl-RPSFPPEE-amino) peptide, like A6, has great potency to inhibit the invasion of cancer cells.

Nature-inspired Enzyme engineering and sustainable catalysis: biochemical clues from the world of plants and extremophiles

The use of enzymes to accelerate chemical reactions for the synthesis of industrially important products is rapidly gaining popularity. Biocatalysis is an environment-friendly approach as it not only uses non-toxic, biodegradable, and renewable raw materials but also helps to reduce waste generation. In this context, enzymes from organisms living in extreme conditions (extremozymes) have been studied extensively and used in industries (food and pharmaceutical), agriculture, and molecular biology, as they are adapted to catalyze reactions withstanding harsh environmental conditions. Enzyme engineering plays a key role in integrating the structure-function insights from reference enzymes and their utilization for developing improvised catalysts. It helps to transform the enzymes to enhance their activity, stability, substrates-specificity, and substrate-versatility by suitably modifying enzyme structure, thereby creating new variants of the enzyme with improved physical and chemical properties. Here, we have illustrated the relatively less-tapped potentials of plant enzymes in general and their sub-class of extremozymes for industrial applications. Plants are exposed to a wide range of abiotic and biotic stresses due to their sessile nature, for which they have developed various mechanisms, including the production of stress-response enzymes. While extremozymes from microorganisms have been extensively studied, there are clear indications that plants and algae also produce extremophilic enzymes as their survival strategy, which may find industrial applications. Typical plant enzymes, such as ascorbate peroxidase, papain, carbonic anhydrase, glycoside hydrolases and others have been examined in this review with respect to their stress-tolerant features and further improvement via enzyme engineering. Some rare instances of plant-derived enzymes that point to greater exploration for industrial use have also been presented here. The overall implication is to utilize biochemical clues from the plant-based enzymes for robust, efficient, and substrate/reaction conditions-versatile scaffolds or reference leads for enzyme engineering.

Designed active-site library reveals thousands of functional GFP variants

Mutations in a protein active site can lead to dramatic and useful changes in protein activity. The active site, however, is sensitive to mutations due to a high density of molecular interactions, substantially reducing the likelihood of obtaining functional multipoint mutants. We introduce an atomistic and machine-learning-based approach, called high-throughput Functional Libraries (htFuncLib), that designs a sequence space in which mutations form low-energy combinations that mitigate the risk of incompatible interactions. We apply htFuncLib to the GFP chromophore-binding pocket, and, using fluorescence readout, recover >16,000 unique designs encoding as many as eight active-site mutations. Many designs exhibit substantial and useful diversity in functional thermostability (up to 96 °C), fluorescence lifetime, and quantum yield. By eliminating incompatible active-site mutations, htFuncLib generates a large diversity of functional sequences. We envision that htFuncLib will be used in one-shot optimization of activity in enzymes, binders, and other proteins.

Improving thermostability of a PL 5 family alginate lyase with combination of rational design strategies

Alginate lyases with strict substrate specificity possess potential in directed production of alginate oligosaccharides with specific composition. However, their poor thermostability hampered their applications in industry. In this study, an efficient comprehensive strategy including sequence-based analysis, structure-based analysis, and computer-aid ΔΔGfold value calculation was proposed. It was successfully performed on alginate lyase (PMD) with strict poly-β-D-mannuronic acid substrate specificity. Four single-point variants A74V, G75V, A240V, and D250G with increased T_m of 3.94 °C, 5.21 °C, 2.56 °C, and 4.80 °C, respectively, were selected out. After ordered combined mutations, a four-point mutant (M4) was finally generated which displayed remarkable increase on thermostability. The T_m of M4 increased from 42.25 °C to 51.59 °C and its half-life at 50 °C was about 58.9-fold of PMD. Meanwhile, there was no obvious loss of enzyme activity (more than 90% retained). Molecular dynamics simulation analysis insisted that the improvement of thermostability might be attribute to the rigidified region A which might be caused by the newly formed hydrogen bonds and salt bridges introduced by mutations, the lower distance of original hydrogen bonds, and the more compact overall structures.

Complete Depolymerization of PET Wastes by an Evolved PET Hydrolase from Directed Evolution

PETase displays great potential in PET depolymerization. Directed evolution has been limited to engineer PETase due to the lack of high-throughput screening assay. In this study, a novel fluorescence-based high-throughput screening assay employing a newly designed substrate, bis (2-hydroxyethyl) 2-hydroxyterephthalate (termed BHET-OH), was developed for PET hydrolases. The best variant DepoPETase produced 1407-fold more products towards amorphous PET film at 50 °C and showed a 23.3 °C higher T_m value than the PETase WT. DepoPETase enabled complete depolymerization of seven untreated PET wastes and 19.1 g PET waste (0.4 % W_enzyme /W_PET ) in liter-scale reactor, suggesting that it is a potential candidate for industrial PET depolymerization processes. The molecular dynamic simulations revealed that the distal substitutions stabilized the loops around the active sites and transmitted the stabilization effect to the active sites through enhancing inter-loop interactions network.

Role of distal sites in enzyme engineering

The limitations associated with natural enzyme catalysis have triggered the rise of the field of protein engineering. Traditional rational design was based on the analysis of protein structural information and catalytic mechanisms to identify key active sites or ligand binding sites to reshape the substrate pocket. The role and significance of functional sites in the active center have been studied extensively. With a deeper understanding of the structure-catalysis relationship map, the entire protein molecule can be filled with residues that play a substantial role in its structure and function. However, the catalytic mechanism underlying distal mutations remains unclear. The aim of this review was to highlight the criticality of the distal site in enzyme engineering based on the following three aspects: What can distal mutations exert on function from mutability landscape? How do distal sites influence enzyme function? How to predict and design distal mutations? This review provides insights into the catalytic mechanism of enzymes from the global interaction network, knowledge from sequence-structure-dynamics-function relationships, and strategies for distal mutation-based protein engineering.

Tuning Selectivity in CalA Lipase: Beyond Tunnel Engineering

Engineering studies of Candida (Pseudozyma) antarctica lipase A (CalA) have demonstrated the potential of this enzyme in the selective hydrolysis of fatty acid esters of different chain lengths. CalA has been shown to bind substrates preferentially through an acyl-chain binding tunnel accessed via the hydrolytic active site; it has also been shown that selectivity for substrates of longer or shorter chain length can be tuned, for instance by modulating steric hindrance within the tunnel. Here we demonstrate that, whereas the tunnel region is certainly of paramount importance for substrate recognition, residues in distal regions of the enzyme can also modulate substrate selectivity. To this end, we investigate variants that carry one or more substitutions within the substrate tunnel as well as in distal regions. Combining experimental determination of the substrate selectivity using natural and synthetic substrates with computational characterization of protein dynamics and of tunnels, we deconvolute the effect of key substitutions and demonstrate that epistatic interactions contribute to procuring selectivity toward either long-chain or short/medium-chain fatty acid esters. We demonstrate that various mechanisms contribute to the diverse selectivity profiles, ranging from reshaping tunnel morphology and tunnel stabilization to obstructing the main substrate-binding tunnel, highlighting the dynamic nature of the substrate-binding region. This work provides important insights into the versatility of this robust lipase toward diverse applications.

IntEnzyDB: an Integrated Structure-Kinetics Enzymology Database

Data-driven modeling has emerged as a new paradigm for biocatalyst design and discovery. Biocatalytic databases that integrate enzyme structure and function data are in urgent need. Here we describe IntEnzyDB as an integrated structure-kinetics database for facile statistical modeling and machine learning. IntEnzyDB employs a relational database architecture with a flattened data structure, which allows rapid data operation. This architecture also makes it easy for IntEnzyDB to incorporate more types of enzyme function data. IntEnzyDB contains enzyme kinetics and structure data from six enzyme commission classes. Using 1050 enzyme structure-kinetics pairs, we investigated the efficiency-perturbing propensities of mutations that are close or distal to the active site. The statistical results show that efficiency-enhancing mutations are globally encoded and that deleterious mutations are much more likely to occur in close mutations than in distal mutations. Finally, we describe a web interface that allows public users to access enzymology data stored in IntEnzyDB. IntEnzyDB will provide a computational facility for data-driven modeling in biocatalysis and molecular evolution.

Protein scaffolds in human clinics

Fundamental clinical areas such as drug delivery and regenerative medicine require biocompatible materials as mechanically stable scaffolds or as nanoscale drug carriers. Among the wide set of emerging biomaterials, polypeptides offer enticing properties over alternative polymers, including full biocompatibility, biodegradability, precise interactivity, structural stability and conformational and functional versatility, all of them tunable by conventional protein engineering. However, proteins from non-human sources elicit immunotoxicities that might bottleneck further development and narrow their clinical applicability. In this context, selecting human proteins or developing humanized protein versions as building blocks is a strict demand to design non-immunogenic protein materials. We review here the expanding catalogue of human or humanized proteins tailored to execute different levels of scaffolding functions and how they can be engineered as self-assembling materials in form of oligomers, polymers or complex networks. In particular, we emphasize those that are under clinical development, revising their fields of applicability and how they have been adapted to offer, apart from mere mechanical support, highly refined functions and precise molecular interactions.

Strategies to Improve the Biosynthesis of β-Lactam Antibiotics by Penicillin G Acylase: Progress and Prospects

β-Lactam antibiotics are widely used anti-infection drugs that are traditionally synthesized via a chemical process. In recent years, with the growing demand for green alternatives, scientists have turned to enzymatic synthesis. Penicillin G acylase (PGA) is the second most commercially used enzyme worldwide with both hydrolytic and synthetic activities toward antibiotics, which has been used to manufacture the key antibiotic nucleus on an industrial level. However, the large-scale application of PGA-catalyzed antibiotics biosynthesis is still in the experimental stage because of some key limitations, such as low substrate concentration, unsatisfactory yield, and lack of superior biocatalysts. This paper systematically reviews the strategies adopted to improve the biosynthesis of β-lactam antibiotics by adjusting the enzymatic property and manipulating the reaction system in recent 20 years, including mining of enzymes, protein engineering, solvent engineering, in situ product removal, and one-pot reaction cascade. These advances will provide important guidelines for the future use of enzymatic synthesis in the industrial production of β-lactam antibiotics.

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.

Computer-Aided Targeted Mutagenesis of Thermoclostridium caenicola d-Allulose 3-Epimerase for Improved Thermostability

d-Allulose 3-epimerase (DAEase) is a key enzyme in d-allulose bioproduction. DAEase from Thermoclostridium caenicola suffers from poor thermostability, hampering its large-scale applications in industry. In this study, mutants A70P, G107P, F155Y, and D162T with increased melting point temperature (T_m) were obtained by targeted mutagenesis based on the calculation of the free energy of folding. The optimal single-point mutant G107P showed 11.08 h, 5, and 5.70 °C increases in the values of half-life (t_1/2) at 60 °C, the optimum temperature (T_opt), and T_m, respectively. Beneficial mutations were combined by ordered recombination mutagenesis, and the combinational mutant Var3 (G107P/F155Y/D162T/A70P) was generated with ΔT_opt of 10 °C and ΔT_m of 12.25 °C. Its t_1/2 value at 65 °C was more than 140 times higher than that of the wild-type enzyme. Molecular dynamics simulations and homology modeling analysis indicated that the enhanced overall rigidity, increased hydrogen bonds between subunits, and redistributed surface electrostatic charges might be responsible for the improved thermostability of the mutant Var3.

Enhancing computational enzyme design by a maximum entropy strategy

Although computational enzyme design is of great importance, the advances utilizing physics-based approaches have been slow, and further progress is urgently needed. One promising direction is using machine learning, but such strategies have not been established as effective tools for predicting the catalytic power of enzymes. Here, we show that the statistical energy inferred from homologous sequences with the maximum entropy (MaxEnt) principle significantly correlates with enzyme catalysis and stability at the active site region and the more distant region, respectively. This finding decodes enzyme architecture and offers a connection between enzyme evolution and the physical chemistry of enzyme catalysis, and it deepens our understanding of the stability-activity trade-off hypothesis for enzymes. Overall, the strong correlations found here provide a powerful way of guiding enzyme design.

Improving the Thermostability and Catalytic Activity of an Inulosucrase by Rational Engineering for the Biosynthesis of Microbial Inulin

Thermostability and enzymatic activity are two vital indexes determining the application of an enzyme on an industrial scale. A truncated inulosucrase, Laga-ISΔ138-702, from Lactobacillus gasseri showed high catalysis activity. To further enhance its thermostability and activity, multiple sequence alignment (MSA) and rational design based on the modeled structure were performed. Variants A446E, S482A, I614M, and A627S were identified with an improved denaturation temperature (T_m) of more than 1 °C. A combinational mutation method was further carried out to explore the synergistic promotion effects of single-point mutants. Additionally, 33 residues at the N-terminus were truncated to construct mutant M4^N-33. The half-life of M4^N-33 at 55 °C increased by 120 times compared to that of Laga-ISΔ138-702, and the relative activity of M4^N-33 increased up to 152% at the optimal pH and temperature (pH 5.5 and 60 °C). Molecular dynamics (MD) simulations illustrated the decreased b-factor of the surface loop of M4^N-33.

Rate-Perturbing Single Amino Acid Mutation for Hydrolases: A Statistical Profiling

Hydrolases are a critical component for modern chemical, pharmaceutical, and environmental sciences. Identifying mutations that enhance catalytic efficiency presents a roadblock to design and to discover new hydrolases for broad academic and industrial uses. Here, we report the statistical profiling for rate-perturbing mutant hydrolases with a single amino acid substitution. We constructed an integrated structure-kinetics database for hydrolases, IntEnzyDB, which contains 3907 k_cats, 4175 K_Ms, and 2715 Protein Data Bank IDs. IntEnzyDB adopts a relational architecture with a flattened data structure, enabling facile and efficient access to clean and tabulated data for machine learning uses. We conducted statistical analyses on how single amino acids mutations influence the turnover number (i.e., k_cat) and efficiency (i.e., k_cat/K_M), with a particular emphasis on profiling the features for rate-enhancing mutations. The results show that mutation to bulky nonpolar residues with a hydrocarbon chain involves a higher likelihood for rate acceleration than to other types of residues. Linear regression models reveal geometric descriptors of substrate and mutation residues that mediate rate-perturbing outcomes for hydrolases with bulky nonpolar mutations. On the basis of the analyses of the structure-kinetics relationship, we observe that the propensity for rate enhancement is independent of protein sizes. In addition, we observe that distal mutations (i.e., >10 Å from the active site) in hydrolases are significantly more prone to induce efficiency neutrality and avoid efficiency deletion but involve similar propensity for rate enhancement. The studies reveal the statistical features for identifying rate-enhancing mutations in hydrolases, which will potentially guide hydrolase discovery in biocatalysis.

Advance in glycosyltransferases, the important bioparts for production of diversified ginsenosides

Ginsenosides are a series of glycosylated triterpenoids predominantly originated from Panax species with multiple pharmacological activities such as anti-aging, mediatory effect on the immune system and the nervous system. During the biosynthesis of ginsenosides, glycosyltransferases play essential roles by transferring various sugar moieties to the sapogenins in contributing to form structure and bioactivity diversified ginsenosides, which makes them important bioparts for synthetic biology-based production of these valuable ginsenosides. In this review, we summarized the functional elucidated glycosyltransferases responsible for ginsenoside biosynthesis, the advance in the protein engineering of UDP-glycosyltransferases (UGTs) and their application with the aim to provide in-depth understanding on ginsenoside-related UGTs for the production of rare ginsenosides applying synthetic biology-based microbial cell factories in the future.

Substitutions at Nonconserved Rheostat Positions Modulate Function by Rewiring Long-Range, Dynamic Interactions

Amino acid substitutions at nonconserved protein positions can have noncanonical and "long-distance" outcomes on protein function. Such outcomes might arise from changes in the internal protein communication network, which is often accompanied by changes in structural flexibility. To test this, we calculated flexibilities and dynamic coupling for positions in the linker region of the lactose repressor protein. This region contains nonconserved positions for which substitutions alter DNA-binding affinity. We first chose to study 11 substitutions at position 52. In computations, substitutions showed long-range effects on flexibilities of DNA-binding positions, and the degree of flexibility change correlated with experimentally measured changes in DNA binding. Substitutions also altered dynamic coupling to DNA-binding positions in a manner that captured other experimentally determined functional changes. Next, we broadened calculations to consider the dynamic coupling between 17 linker positions and the DNA-binding domain. Experimentally, these linker positions exhibited a wide range of substitution outcomes: Four conserved positions tolerated hardly any substitutions ("toggle"), ten nonconserved positions showed progressive changes from a range of substitutions ("rheostat"), and three nonconserved positions tolerated almost all substitutions ("neutral"). In computations with wild-type lactose repressor protein, the dynamic couplings between the DNA-binding domain and these linker positions showed varied degrees of asymmetry that correlated with the observed toggle/rheostat/neutral substitution outcomes. Thus, we propose that long-range and noncanonical substitutions outcomes at nonconserved positions arise from rewiring long-range communication among functionally important positions. Such calculations might enable predictions for substitution outcomes at a range of nonconserved positions.

Computational Enzyme Engineering Pipelines for Optimized Production of Renewable Chemicals

To enable a sustainable supply of chemicals, novel biotechnological solutions are required that replace the reliance on fossil resources. One potential solution is to utilize tailored biosynthetic modules for the metabolic conversion of CO₂ or organic waste to chemicals and fuel by microorganisms. Currently, it is challenging to commercialize biotechnological processes for renewable chemical biomanufacturing because of a lack of highly active and specific biocatalysts. As experimental methods to engineer biocatalysts are time- and cost-intensive, it is important to establish efficient and reliable computational tools that can speed up the identification or optimization of selective, highly active, and stable enzyme variants for utilization in the biotechnological industry. Here, we review and suggest combinations of effective state-of-the-art software and online tools available for computational enzyme engineering pipelines to optimize metabolic pathways for the biosynthesis of renewable chemicals. Using examples relevant for biotechnology, we explain the underlying principles of enzyme engineering and design and illuminate future directions for automated optimization of biocatalysts for the assembly of synthetic metabolic pathways.

Recent advances in user-friendly computational tools to engineer protein function

Progress in technology and algorithms throughout the past decade has transformed the field of protein design and engineering. Computational approaches have become well-engrained in the processes of tailoring proteins for various biotechnological applications. Many tools and methods are developed and upgraded each year to satisfy the increasing demands and challenges of protein engineering. To help protein engineers and bioinformaticians navigate this emerging wave of dedicated software, we have critically evaluated recent additions to the toolbox regarding their application for semi-rational and rational protein engineering. These newly developed tools identify and prioritize hotspots and analyze the effects of mutations for a variety of properties, comprising ligand binding, protein-protein and protein-nucleic acid interactions, and electrostatic potential. We also discuss notable progress to target elusive protein dynamics and associated properties like ligand-transport processes and allosteric communication. Finally, we discuss several challenges these tools face and provide our perspectives on the further development of readily applicable methods to guide protein engineering efforts.

How a Scorpion Toxin Selectively Captures a Prey Sodium Channel: The Molecular and Evolutionary Basis Uncovered

The growing resistance of insects to chemical pesticides is reducing the effectiveness of conventional methods for pest control and thus, the development of novel insecticidal agents is imperative. Scorpion toxins specific for insect voltage-gated sodium channels (Navs) have been considered as one of the most promising insecticide alternatives due to their host specificity, rapidly evoked toxicity, biodegradability, and the lack of resistance. However, they have not been developed for uses in agriculture and public health, mainly because of a limited understanding of their molecular and evolutionary basis controlling their phylogenetic selectivity. Here, we show that the traditionally defined insect-selective scorpion toxin LqhIT2 specifically captures a prey Nav through a conserved trapping apparatus comprising a three-residue-formed cavity and a structurally adjacent leucine. The former serves as a detector to recognize and bind a highly exposed channel residue conserved in insects and spiders, two major prey items for scorpions; and the latter subsequently seizes the "moving" voltage sensor via hydrophobic interactions to reduce activation energy for channel opening, demonstrating its action in an enzyme-like manner. Based on the established toxin-channel interaction model in combination with toxicity assay, we enlarged the toxic spectrum of LqhIT2 to spiders and certain other arthropods. Furthermore, we found that genetic background-dependent cavity shapes determine the species selectivity of LqhIT2-related toxins. We expect that the discovery of the trapping apparatus will improve our understanding of the evolution and design principle of Nav-targeted toxins from a diversity of arthropod predators and accelerate their uses in pest control.

Hinge-shift mechanism as a protein design principle for the evolution of β-lactamases from substrate promiscuity to specificity

TEM-1 β-lactamase degrades β-lactam antibiotics with a strong preference for penicillins. Sequence reconstruction studies indicate that it evolved from ancestral enzymes that degraded a variety of β-lactam antibiotics with moderate efficiency. This generalist to specialist conversion involved more than 100 mutational changes, but conserved fold and catalytic residues, suggesting a role for dynamics in enzyme evolution. Here, we develop a conformational dynamics computational approach to rationally mold a protein flexibility profile on the basis of a hinge-shift mechanism. By deliberately weighting and altering the conformational dynamics of a putative Precambrian β-lactamase, we engineer enzyme specificity that mimics the modern TEM-1 β-lactamase with only 21 amino acid replacements. Our conformational dynamics design thus re-enacts the evolutionary process and provides a rational allosteric approach for manipulating function while conserving the enzyme active site.

Mutational and biophysical robustness in a prestabilized monobody

The fibronectin type III (FN3) monobody domain is a promising non-antibody scaffold, which features a less complex architecture than an antibody while maintaining analogous binding loops. We previously developed FN3Con, a hyperstable monobody derivative with diagnostic and therapeutic potential. Prestabilization of the scaffold mitigates the stability–function trade-off commonly associated with evolving a protein domain toward biological activity. Here, we aimed to examine if the FN3Con monobody could take on antibody-like binding to therapeutic targets, while retaining its extreme stability. We targeted the first of the Adnectin derivative of monobodies to reach clinical trials, which was engineered by directed evolution for binding to the therapeutic target VEGFR2; however, this function was gained at the expense of large losses in thermostability and increased oligomerization. In order to mitigate these losses, we grafted the binding loops from Adnectin-anti-VEGFR2 (CT-322) onto the prestabilized FN3Con scaffold to produce a domain that successfully bound with high affinity to the therapeutic target VEGFR2. This FN3Con-anti-VEGFR2 construct also maintains high thermostability, including remarkable long-term stability, retaining binding activity after 2 years of storage at 36 °C. Further investigations into buffer excipients doubled the presence of monomeric monobody in accelerated stability trials. These data suggest that loop grafting onto a prestabilized scaffold is a viable strategy for the development of monobody domains with desirable biophysical characteristics and that FN3Con is therefore well-suited to applications such as the evolution of multiple paratopes or shelf-stable diagnostics and therapeutics.

Directed Enzyme Evolution and Encapsulation in Peptide Nanospheres of Quorum Quenching Lactonase as an Antibacterial Treatment against Plant Pathogen

The need to increase agricultural yield has led to an extensive use of antibiotics against plant pathogens, which has resulted in the emergence of resistant strains. Therefore, there is an increasing demand for new methods, preferably with lower chances of developing resistant strains and a lower risk to the environment or public health. Many Gram-negative bacterial pathogens use quorum sensing, a population-density-dependent regulatory mechanism, to monitor the secretion of N-acyl-homoserine lactones (AHLs) and pathogenicity. Therefore, quorum sensing represents an attractive antivirulence target. AHL lactonases hydrolyze AHLs and have potential antibacterial properties; however, their use is limited by thermal instability and durability, or low activity. Here, we demonstrate that an AHL lactonase from the phosphotriesterase-like lactonase family exhibits high activity with the AHL secreted from the plant pathogen Erwinia amylovora and attenuates infection in planta. Using directed enzyme evolution, we were able to increase the enzyme's temperature resistance (T₅₀, the temperature at which 50% of the activity is retained) by 8 °C. Then, by performing enzyme encapsulation in nanospherical capsules composed of tertbutoxycarbonyl-Phe-Phe-OH peptide, the shelf life was extended for more than 5 weeks. Furthermore, the encapsulated and free mutant were able to significantly inhibit up to 70% blossom's infection in the field, achieving the same efficacy as seen with antibiotics commonly used today to treat the plant pathogen. We conclude that specific AHL lactonase can inhibit E. amylovora infection in the field, as it degrades the AHL secreted by this plant pathogen. The combination of directed enzyme evolution and peptide nanostructure encapsulation significantly improved the thermal resistance and shelf life of the enzyme, respectively, increasing its potential in future development as antibacterial treatment.

Evolution of new enzymes by gene duplication and divergence

Thousands of new metabolic and regulatory enzymes have evolved by gene duplication and divergence since the dawn of life. New enzyme activities often originate from promiscuous secondary activities that have become important for fitness due to a change in the environment or a mutation. Mutations that make a promiscuous activity physiologically relevant can occur in the gene encoding the promiscuous enzyme itself, but can also occur elsewhere, resulting in increased expression of the enzyme or decreased competition between the native and novel substrates for the active site. If a newly useful activity is inefficient, gene duplication/amplification will set the stage for divergence of a new enzyme. Even a few mutations can increase the efficiency of a new activity by orders of magnitude. As efficiency increases, amplified gene arrays will shrink to provide two alleles, one encoding the original enzyme and one encoding the new enzyme. Ultimately, genomic rearrangements eliminate co-amplified genes and move newly evolved paralogs to a distant region of the genome.

UMI-linked consensus sequencing enables phylogenetic analysis of directed evolution

The success of protein evolution campaigns is strongly dependent on the sequence context in which mutations are introduced, stemming from pervasive non-additive interactions between a protein's amino acids ('intra-gene epistasis'). Our limited understanding of such epistasis hinders the correct prediction of the functional contributions and adaptive potential of mutations. Here we present a straightforward unique molecular identifier (UMI)-linked consensus sequencing workflow (UMIC-seq) that simplifies mapping of evolutionary trajectories based on full-length sequences. Attaching UMIs to gene variants allows accurate consensus generation for closely related genes with nanopore sequencing. We exemplify the utility of this approach by reconstructing the artificial phylogeny emerging in three rounds of directed evolution of an amine dehydrogenase biocatalyst via ultrahigh throughput droplet screening. Uniquely, we are able to identify lineages and their founding variant, as well as non-additive interactions between mutations within a full gene showing sign epistasis. Access to deep and accurate long reads will facilitate prediction of key beneficial mutations and adaptive potential based on in silico analysis of large sequence datasets.

Ancestral sequence reconstruction produces thermally stable enzymes with mesophilic enzyme-like catalytic properties

Enzymes have high catalytic efficiency and low environmental impact, and are therefore potentially useful tools for various industrial processes. Crucially, however, natural enzymes do not always have the properties required for specific processes. It may be necessary, therefore, to design, engineer, and evolve enzymes with properties that are not found in natural enzymes. In particular, the creation of enzymes that are thermally stable and catalytically active at low temperature is desirable for processes involving both high and low temperatures. In the current study, we designed two ancestral sequences of 3-isopropylmalate dehydrogenase by an ancestral sequence reconstruction technique based on a phylogenetic analysis of extant homologous amino acid sequences. Genes encoding the designed sequences were artificially synthesized and expressed in Escherichia coli. The reconstructed enzymes were found to be slightly more thermally stable than the extant thermophilic homologue from Thermus thermophilus. Moreover, they had considerably higher low-temperature catalytic activity as compared with the T. thermophilus enzyme. Detailed analyses of their temperature-dependent specific activities and kinetic properties showed that the reconstructed enzymes have catalytic properties similar to those of mesophilic homologues. Collectively, our study demonstrates that ancestral sequence reconstruction can produce a thermally stable enzyme with catalytic properties adapted to low-temperature reactions.

The hydrophobicity of an amino acid residue in a flexible loop of KP-43 protease alters activity toward a macromolecule substrate

Abstract

KP-43, a 43-kDa alkaline serine protease, is resistant to chemical oxidants and surfactants, making it suitable for use in laundry detergents. An amino acid residue at position 195, in a unique flexible loop that binds a Ca²⁺ ion, dramatically affects the proteolytic activity and thermal stability of KP-43. In the present study, we obtained 20 variants with substitutions at position 195 and investigated how these residues affect hydrolytic activity toward a macromolecular substrate (casein) and a synthetic tetra-peptide (AAPL). At pH 10, the variant with the highest caseinolytic activity, Tyr195Gln, exhibited 4.4-fold higher activity than the variant with the lowest caseinolytic activity, Tyr195Trp. A significant negative correlation was observed between the hydrophobicity of the residue at position 195 and caseinolytic activity at pH 8–10. At pH 7, the correlation became weak; at pH 6, the correlation reversed to positive. Unlike casein, in the case of hydrolysis of AAPL, no correlation was observed at pH 10 or pH 6. Because the amino acid residue at position 195 is located on the protein surface and considered sufficiently far from the active cleft, the variation in caseinolytic activity between the 20 variants was attributed to changes in interaction efficiency with different states of casein at different pH values. To improve the enzymatic activity, we propose substituting amino acid residues on the protein surface to change the efficiency of interaction with the macromolecular substrates.

Key points

• A single amino acid residue on the protein surface markedly changed enzyme activity.

• The hydrophobicity of the amino acid residue and enzyme activity had a correlation.

• The key amino acid residue for substrate recognition exists on the protein surface.

Electronic supplementary material

The online version of this article (10.1007/s00253-020-10826-2) contains supplementary material, which is available to authorized users.

The 3-ketoacyl-CoA thiolase: an engineered enzyme for carbon chain elongation of chemical compounds

Because of their function of catalyzing the rearrangement of the carbon chains, thiolases have attracted increasing attentions over the past decades. The 3-ketoacyl-CoA thiolase (KAT) is a member of the thiolase, which is capable of catalyzing the Claisen condensation reaction between the two acyl-CoAs, thereby achieving carbon chain elongation. In this way, diverse value-added compounds might be synthesized starting from simple small CoA thioesters. However, most KATs are hampered by low stability and poor substrate specificity, which has hindered the development of large-scale biosynthesis. In this review, the common characteristics in the three-dimensional structure of KATs from different sources are summarized. Moreover, structure-guided rational engineering is discussed as a strategy for enhancing the performance of KATs. Finally, we reviewed the metabolic engineering applications of KATs for producing various energy-storage molecules, such as n-butanol, fatty acids, dicarboxylic acids, and polyhydroxyalkanoates. KEY POINTS: • Summarize the structural characteristics and catalyzation mechanisms of KATs. • Review on the rational engineering to enhance the performance of KATs. • Discuss the applications of KATs for producing energy-storage molecules.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

Practically useful protein-design methods combining phylogenetic and atomistic calculations

Our ability to design new or improved biomolecular activities depends on understanding the sequence-function relationships in proteins. The large size and fold complexity of most proteins, however, obscure these relationships, and protein-optimization methods continue to rely on laborious experimental iterations. Recently, a deeper understanding of the roles of stability-threshold effects and biomolecular epistasis in proteins has led to the development of hybrid methods that combine phylogenetic analysis with atomistic design calculations. These methods enable reliable and even single-step optimization of protein stability, expressibility, and activity in proteins that were considered outside the scope of computational design. Furthermore, ancestral-sequence reconstruction produces insights on missing links in the evolution of enzymes and binders that may be used in protein design. Through the combination of phylogenetic and atomistic calculations, the long-standing goal of general computational methods that can be universally applied to study and optimize proteins finally seems within reach.

One Pot Use of Combilipases for Full Modification of Oils and Fats: Multifunctional and Heterogeneous Substrates

Lipases are among the most utilized enzymes in biocatalysis. In many instances, the main reason for their use is their high specificity or selectivity. However, when full modification of a multifunctional and heterogeneous substrate is pursued, enzyme selectivity and specificity become a problem. This is the case of hydrolysis of oils and fats to produce free fatty acids or their alcoholysis to produce biodiesel, which can be considered cascade reactions. In these cases, to the original heterogeneity of the substrate, the presence of intermediate products, such as diglycerides or monoglycerides, can be an additional drawback. Using these heterogeneous substrates, enzyme specificity can promote that some substrates (initial substrates or intermediate products) may not be recognized as such (in the worst case scenario they may be acting as inhibitors) by the enzyme, causing yields and reaction rates to drop. To solve this situation, a mixture of lipases with different specificity, selectivity and differently affected by the reaction conditions can offer much better results than the use of a single lipase exhibiting a very high initial activity or even the best global reaction course. This mixture of lipases from different sources has been called “combilipases” and is becoming increasingly popular. They include the use of liquid lipase formulations or immobilized lipases. In some instances, the lipases have been coimmobilized. Some discussion is offered regarding the problems that this coimmobilization may give rise to, and some strategies to solve some of these problems are proposed. The use of combilipases in the future may be extended to other processes and enzymes.

A beginner's guide to molecular dynamics simulations and the identification of cross-correlation networks for enzyme engineering

The functional properties of proteins are decided not only by their relatively rigid overall structures, but even more importantly, by their dynamic properties. In a protein, some regions of structure exhibit highly correlated or anti-correlated motions with others, some are highly dynamic but uncorrelated, while other regions are relatively static. The residues with correlated or anti-correlated motions can form a so-called dynamic cross-correlation network, through which information can be transmitted. Such networks have been shown to be critical to allosteric transitions, and ligand binding, and have also been shown to be able to mediate epistatic interactions between mutations. As a result, they are likely to play a significant role in the development of new enzyme engineering strategies. In this chapter, protocols are provided for the assessment of dynamic cross-correlation networks, and for their application in protein engineering. Transketolase from E. coli is used as a model and the software GROMACS is applied for carrying out MD simulations to generate trajectories containing structural ensembles. The trajectory is then used for a dynamic cross correlation analysis using the R package, Bio3D. A matrix of all atom-wise cross-correlation coefficients is finally obtained, which can be displayed in a graphical representation termed a dynamical cross-correlation matrix.

The Role of Conformational Dynamics and Allostery in Modulating Protein Evolution

Advances in sequencing techniques and statistical methods have made it possible not only to predict sequences of ancestral proteins but also to identify thousands of mutations in the human exome, some of which are disease associated. These developments have motivated numerous theories and raised many questions regarding the fundamental principles behind protein evolution, which have been traditionally investigated horizontally using the tip of the phylogenetic tree through comparative studies of extant proteins within a family. In this article, we review a vertical comparison of the modern and resurrected ancestral proteins. We focus mainly on the dynamical properties responsible for a protein's ability to adapt new functions in response to environmental changes. Using the Dynamic Flexibility Index and the Dynamic Coupling Index to quantify the relative flexibility and dynamic coupling at a site-specific, single-amino-acid level, we provide evidence that the migration of hinges, which are often functionally critical rigid sites, is a mechanism through which proteins can rapidly evolve. Additionally, we show that disease-associated mutations in proteins often result in flexibility changes even at positions distal from mutational sites, particularly in the modulation of active site dynamics.

Semi-rational hinge engineering: modulating the conformational transformation of glutamate dehydrogenase for enhanced reductive amination activity towards non-natural substrates

The hinge region was identified to be a promising hotspot for activity engineering of GluDHs, providing a potent alternative for developing high-performance biocatalysts toward valuable optically pure l-amino acid production.

Tailoring Tryptophan Synthase TrpB for Selective Quaternary Carbon Bond Formation

We previously engineered the β-subunit of tryptophan synthase (TrpB), which catalyzes the condensation of l-serine and indole to l-tryptophan, to synthesize a range of noncanonical amino acids from l-serine and indole derivatives or other nucleophiles. Here we employ directed evolution to engineer TrpB to accept 3-substituted oxindoles and form C-C bonds leading to new quaternary stereocenters. Initially, the variants that could use 3-substituted oxindoles preferentially formed N-C bonds on N1 of the substrate. Protecting N1 encouraged evolution toward C-alkylation, which persisted when protection was removed. Six generations of directed evolution resulted in TrpB Pfquat with a 400-fold improvement in activity for alkylation of 3-substituted oxindoles and the ability to selectively form a new, all-carbon quaternary stereocenter at the γ-position of the amino acid products. The enzyme can also alkylate and form all-carbon quaternary stereocenters on structurally similar lactones and ketones, where it exhibits excellent regioselectivity for the tertiary carbon. The configurations of the γ-stereocenters of two of the products were determined via microcrystal electron diffraction (MicroED), and we report the MicroED structure of a small molecule obtained using the Falcon III direct electron detector. Highly thermostable and expressed at >500 mg/L E. coli culture, TrpB Pfquat offers an efficient, sustainable, and selective platform for the construction of diverse noncanonical amino acids bearing all-carbon quaternary stereocenters.