Handicap-Recover Evolution Leads to a Chemically Versatile, Nucleophile-Permissive Protease

Ultrahigh-Throughput Directed Evolution of a Metal-Free α/β-Hydrolase with a Cys-His-Asp Triad into an Efficient Phosphotriesterase

Finding new mechanistic solutions for biocatalytic challenges is key in the evolutionary adaptation of enzymes, as well as in devising new catalysts. The recent release of man-made substances into the environment provides a dynamic testing ground for observing biocatalytic innovation at play. Phosphate triesters, used as pesticides, have only recently been introduced into the environment, where they have no natural counterpart. Enzymes have rapidly evolved to hydrolyze phosphate triesters in response to this challenge, converging onto the same mechanistic solution, which requires bivalent cations as a cofactor for catalysis. In contrast, the previously identified metagenomic promiscuous hydrolase P91, a homologue of acetylcholinesterase, achieves slow phosphotriester hydrolysis mediated by a metal-independent Cys-His-Asp triad. Here, we probe the evolvability of this new catalytic motif by subjecting P91 to directed evolution. By combining a focused library approach with the ultrahigh throughput of droplet microfluidics, we increase P91's activity by a factor of ≈360 (to a k_cat/K_M of ≈7 × 10⁵ M^-1 s^-1) in only two rounds of evolution, rivaling the catalytic efficiencies of naturally evolved, metal-dependent phosphotriesterases. Unlike its homologue acetylcholinesterase, P91 does not suffer suicide inhibition; instead, fast dephosphorylation rates make the formation of the covalent adduct rather than its hydrolysis rate-limiting. This step is improved by directed evolution, with intermediate formation accelerated by 2 orders of magnitude. Combining focused, combinatorial libraries with the ultrahigh throughput of droplet microfluidics can be leveraged to identify and enhance mechanistic strategies that have not reached high efficiency in nature, resulting in alternative reagents with novel catalytic machineries.

Between Protein Fold and Nucleophile Identity: Multiscale Modeling of the TEV Protease Enzyme-Substrate Complex

The cysteine protease from the tobacco etch virus (TEVp) is a well-known and widely utilized enzyme. TEVp's chymotrypsin-like fold is generally associated with serine catalytic triads that differ in terms of a reaction mechanism from the most well-studied papain-like cysteine proteases. The question of what dominates the TEVp mechanism, nucleophile identity, or structural composition has never been previously addressed. Here, we use enhanced sampling multiscale modeling to uncover that TEVp combines the features of two worlds in such a way that potentially hampers its activity. We show that TEVp cysteine is strictly in the anionic form in a free enzyme similar to papain. Peptide binding shifts the equilibrium toward the nucleophile's protonated form, characteristic of chymotrypsin-like proteases, although the cysteinyl anion form is still present and interconversion is rapid. This way cysteine protonation generates enzyme states that are a diversion from the most effective course of action, with only 13.2% of Michaelis complex sub-states able to initiate the reaction. As a result, we propose an updated view on the reaction mechanism catalyzed by TEVp. We also demonstrate that AlphaFold is able to construct protease-substrate complexes with high accuracy. We propose that our findings open a way for its industrious use in enzymological tasks. Unique features of TEVp discovered in this work open a discussion on the evolutionary history and trade-offs of optimizing serine triad-associated folds to cysteine as a nucleophile.

Construction, Investigation and Application of TEV Protease Variants with Improved Oxidative Stability

Tobacco etch virus protease (TEVp) is a useful tool for removing fusion tags, but wild-type TEVp is less stable under oxidized redox state. In this work, we introduced and combined C19S, C110S and C130S into TEVp variants containing T17S, L56V, N68D, I77V and S135G to improve protein solubility, and S219V to inhibit self-proteolysis. The solubility and cleavage activity of the constructed variants in Escherichia coli strains including BL21(DE3), BL21(DE3)pLys, Rossetta(DE3) and Origami(DE3) under the same induction conditions were analyzed and compared. The desirable soluble amounts, activity, and oxidative stability were identified to be reluctantly favored in the TEVp. Unlike C19S, C110S and C130S hardly impacted on decreasing protein solubility in the BL21(DE3), but they contributed to improved tolerance to the oxidative redox state in vivo and in vitro. After two fusion proteins were cleaved by purified TEVp protein containing double mutations under the oxidized redox state, the refolded disulfide-rich bovine enterokinase catalytic domain or maize peroxidase with enhanced yields were released from the regenerated amorphous cellulose via affinity absorption of the cellulose-binding module as the affinity tag.

Discovery of an ultraspecific triuret hydrolase (TrtA) establishes the triuret biodegradation pathway

Triuret (carbonyldiurea) is an impurity found in industrial urea fertilizer (<0.1% w/w) that is applied, worldwide, around 300 million pounds each year on agricultural lands. In addition to anthropogenic sources, endogenous triuret has been identified in amoeba and human urine, the latter being diagnostic for hypokalemia. The present study is the first to describe the metabolic breakdown of triuret, which funnels into biuret metabolism. We identified the gene responsible for triuret decomposition (trtA) in bacterial genomes, clustered with biuH, which encodes biuret hydrolase and has close protein sequence homology. TrtA is a member of the isochorismatase-like hydrolase (IHL) protein family, similarly to BiuH, and has a catalytic efficiency (k_cat/K_M) of 6 x 10⁵ M⁻¹s⁻¹, a K_M for triuret of 20 μM, and exquisite substrate specificity. Indeed, TrtA has four orders of magnitude less activity with biuret. Crystal structures of TrtA in apo and holo form were solved and compared with the BiuH structure. The high substrate selectivity was found to be conveyed by second shell residues around each active site. Mutagenesis of residues conserved in TrtA to the alternate consensus found in BiuHs revealed residues critical to triuret hydrolase activity but no single mutant evolved more biuret activity, and likely a combination of mutations is required to interconvert between TrtA, BiuH functions. TrtA-mediated triuret metabolism is relatively rare in recorded genomes (1–2%), but is largely found in plant-associated, nodulating, and endophytic bacteria. This study suggests functions for triuret hydrolase in certain eukaryotic intermediary processes and prokaryotic intermediary or biodegradative metabolism

Differences in protein structural regions that impact functional specificity in GT2 family β-glucan synthases

Most cell wall and secreted β-glucans are synthesised by the CAZy Glycosyltransferase 2 family (www.cazy.org), with different members catalysing the formation of (1,4)-β-, (1,3)-β-, or both (1,4)- and (1,3)-β-glucosidic linkages. Given the distinct physicochemical properties of each of the resultant β-glucans (cellulose, curdlan, and mixed linkage glucan, respectively) are crucial to their biological and biotechnological functions, there is a desire to understand the molecular evolution of synthesis and how linkage specificity is determined. With structural studies hamstrung by the instability of these proteins to solubilisation, we have utilised in silico techniques and the crystal structure for a bacterial cellulose synthase to further understand how these enzymes have evolved distinct functions. Sequence and phylogenetic analyses were performed to determine amino acid conservation, both family-wide and within each sub-family. Further structural analysis centred on comparison of a bacterial curdlan synthase homology model with the bacterial cellulose synthase crystal structure, with molecular dynamics simulations performed with their respective β-glucan products bound in the trans-membrane channel. Key residues that differentially interact with the different β-glucan chains and have sub-family-specific conservation were found to reside at the entrance of the trans-membrane channel. The linkage-specific catalytic activity of these enzymes and hence the type of β-glucan chain built is thus likely determined by the different interactions between the proteins and the first few glucose residues in the channel, which in turn dictates the position of the acceptor glucose. The sequence-function relationships for the bacterial β-glucan synthases pave the way for extending this understanding to other kingdoms, such as plants.

Molecular basis for the production of cyclic peptides by plant asparaginyl endopeptidases

Asparaginyl endopeptidases (AEPs) are proteases that have crucial roles in plant defense and seed storage protein maturation. Select plant AEPs, however, do not function as proteases but as transpeptidases (ligases) catalyzing the intra-molecular ligation of peptide termini, which leads to peptide cyclization. These ligase-type AEPs have potential biotechnological applications ranging from in vitro peptide engineering to plant molecular farming, but the structural features enabling these enzymes to catalyze peptide ligation/cyclization rather than proteolysis are currently unknown. Here, we compare the sequences, structures, and functions of diverse plant AEPs by combining molecular modeling, sequence space analysis, and functional testing in planta. We find that changes within the substrate-binding pocket and an adjacent loop, here named the “marker of ligase activity”, together play a key role for AEP ligase efficiency. Identification of these structural determinants may facilitate the discovery of more ligase-type AEPs and the engineering of AEPs with tailored catalytic properties.

Asparaginyl endopeptidases (AEPs) are plant proteases that can also function as ligases, catalyzing the production of cyclic plant peptides. Here, the authors identify structural features that govern AEP ligase activity, providing insights to aid the discovery and engineering of ligase-type AEPs.

Identification of (S)-selective transaminases for the asymmetric synthesis of bulky chiral amines

The use of transaminases to access pharmaceutically relevant chiral amines is an attractive alternative to transition-metal-catalysed asymmetric chemical synthesis. However, one major challenge is their limited substrate scope. Here we report the creation of highly active and stereoselective transaminases starting from fold class I. The transaminases were developed by extensive protein engineering followed by optimization of the identified motif. The resulting enzymes exhibited up to 8,900-fold higher activity than the starting scaffold and are highly stereoselective (up to >99.9% enantiomeric excess) in the asymmetric synthesis of a set of chiral amines bearing bulky substituents. These enzymes should therefore be suitable for use in the synthesis of a wide array of potential intermediates for pharmaceuticals. We also show that the motif can be engineered into other protein scaffolds with sequence identities as low as 70%, and as such should have a broad impact in the field of biocatalytic synthesis and enzyme engineering.

Reverse evolution leads to genotypic incompatibility despite functional and active site convergence

Understanding the extent to which enzyme evolution is reversible can shed light on the fundamental relationship between protein sequence, structure, and function. Here, we perform an experimental test of evolutionary reversibility using directed evolution from a phosphotriesterase to an arylesterase, and back, and examine the underlying molecular basis. We find that wild-type phosphotriesterase function could be restored (>10⁴-fold activity increase), but via an alternative set of mutations. The enzyme active site converged towards its original state, indicating evolutionary constraints imposed by catalytic requirements. We reveal that extensive epistasis prevents reversions and necessitates fixation of new mutations, leading to a functionally identical sequence. Many amino acid exchanges between the new and original enzyme are not tolerated, implying sequence incompatibility. Therefore, the evolution was phenotypically reversible but genotypically irreversible. Our study illustrates that the enzyme's adaptive landscape is highly rugged, and different functional sequences may constitute separate fitness peaks.

DOI:http://dx.doi.org/10.7554/eLife.06492.001

Enzymes in bacteria and other organisms are built following instructions contained within each cell's DNA. Changes in the DNA, that is to say, mutations, can alter the shape and activity of the enzymes that are produced, which can ultimately affect the ability of the organism to survive and reproduce. Mutations that are beneficial to the organism are more likely to be passed on to future generations, which can lead to populations changing over time.

The DNA sequences that an organism carries are referred to as its ‘genotype’ and the resulting physical characteristics of the organism are known as its ‘phenotype’. Studies of evolution tend to focus on how particular species or molecules become more different over time. However, one area that remains controversial is whether it is possible for evolution to be reversed so that an organism or molecule returns to a previous form.

An enzyme called PTE is said to have phosphotriesterase activity because it catalyzes this particular type of chemical reaction. Recently, a group of researchers used a method called ‘directed evolution’ to demonstrate that it is possible for PTE to evolve in a way that means it loses its phosphotriesterase activity and becomes able to catalyze a different type of chemical reaction. Here, Kaltenbach et al.—including some of the researchers from the previous work—investigated whether it was possible to use the same method to reverse this evolution and restore the enzyme's original activity.

The experiments show that reverse evolution is possible as phosphotriesterase activity was restored to the PTE enzyme from the previous study. However, although the phenotype of the final enzyme matched that of the original PTE enzyme, the genotypes did not match as the DNA sequences of the genes that encode these enzymes differ. The DNA does not revert to its original sequence because the effect of individual mutations on the phenotype depends on what other mutations are present. For example, as the enzyme evolved its new activity, additional mutations accumulated that did not alter enzyme activity. During the reverse evolution experiment, some of these mutations could have started to exert influence on the phenotype so that different mutations were required to restore the phosphotriesterase activity.

In the future, Kaltenbach et al.'s findings may aid efforts to engineer artificial enzymes for use in medicine or industry.

DOI:http://dx.doi.org/10.7554/eLife.06492.002