Biochemistry. Loop grafting and the origins of enzyme species

Remote loop evolution reveals a complex biological function for chitinase enzymes beyond the active site

Loops are small secondary structural elements that play a crucial role in the emergence of new enzyme functions. However, the evolutionary molecular mechanisms how proteins acquire these loop elements and obtain new function is poorly understood. To address this question, we study glycoside hydrolase family 19 (GH19) chitinase-an essential enzyme family for pathogen degradation in plants. By revealing the evolutionary history and loops appearance of GH19 chitinase, we discover that one loop which is remote from the catalytic site, is necessary to acquire the new antifungal activity. We demonstrate that this remote loop directly accesses the fungal cell wall, and surprisingly, it needs to adopt a defined structure supported by long-range intramolecular interactions to perform its function. Our findings prove that nature applies this strategy at the molecular level to achieve a complex biological function while maintaining the original activity in the catalytic pocket, suggesting an alternative way to design new enzyme function.

Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers

The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.

Engineering Substrate Promiscuity of Nucleoside Phosphorylase Via an Insertions-Deletions Strategy

Nucleoside phosphorylases (NPs) are the key enzymes in the nucleoside metabolism pathway and are widely employed for the synthesis of nucleoside analogs, which are difficult to access via conventional synthetic methods. NPs are generally classified as purine nucleoside phosphorylase (PNP) and pyrimidine or uridine nucleoside phosphorylase (PyNP/UP), based on their substrate preference. Here, based on the evolutionary information on the NP-I family, we adopted an insertions-deletions (InDels) strategy to engineer the substrate promiscuity of nucleoside phosphorylase AmPNPΔS2V102 K, which exhibits both PNP and UP activities from a trimeric PNP (AmPNP) of Aneurinibacillus migulanus. Furthermore, the AmPNPΔS2V102 K exerted phosphorylation activities toward arabinose nucleoside, fluorosyl nucleoside, and dideoxyribose, thereby broadening the unnatural-ribose nucleoside substrate spectrum of AmPNP. Finally, six purine nucleoside analogues were successfully synthesized, using the engineered AmPNPΔS2V102 K instead of the traditional "two-enzymes PNP/UP" approach. These results provide deep insights into the catalytic mechanisms of the PNP and demonstrate the benefits of using the InDels strategy to achieve substrate promiscuity in an enzyme, as well as broadening the substrate spectrum of the enzyme.

Sequence - dynamics - function relationships in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) are crucial regulators of cellular signaling. Their activity is regulated by the motion of a conserved loop, the WPD-loop, from a catalytically inactive open to a catalytically active closed conformation. WPD-loop motion optimally positions a catalytically critical residue into the active site, and is directly linked to the turnover number of these enzymes. Crystal structures of chimeric PTPs constructed by grafting parts of the WPD-loop sequence of PTP1B onto the scaffold of YopH showed WPD-loops in a wide-open conformation never previously observed in either parent enzyme. This wide-open conformation has, however, been observed upon binding of small molecule inhibitors to other PTPs, suggesting the potential of targeting it for drug discovery efforts. Here, we have performed simulations of both enzymes and show that there are negligible energetic differences in the chemical step of catalysis, but significant differences in the dynamical properties of the WPD-loop. Detailed interaction network analysis provides insight into the molecular basis for this population shift to a wide-open conformation. Taken together, our study provides insight into the links between loop dynamics and chemistry in these YopH variants specifically, and how WPD-loop dynamic can be engineered through modification of the internal protein interaction network.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

AlphaFold, allosteric, and orthosteric drug discovery: Ways forward

Drug discovery is arguably a highly challenging and significant interdisciplinary aim. The stunning success of the artificial intelligence-powered AlphaFold, whose latest version is buttressed by an innovative machine-learning approach that integrates physical and biological knowledge about protein structures, raised drug discovery hopes that unsurprisingly, have not come to bear. Even though accurate, the models are rigid, including the drug pockets. AlphaFold's mixed performance poses the question of how its power can be harnessed in drug discovery. Here we discuss possible ways of going forward wielding its strengths, while bearing in mind what AlphaFold can and cannot do. For kinases and receptors, an input enriched in active (ON) state models can better AlphaFold's chance of rational drug design success.

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

Insights into the importance of WPD-loop sequence for activity and structure in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) possess a conserved mobile catalytic loop, the WPD-loop, which brings an aspartic acid into the active site where it acts as an acid/base catalyst. Prior experimental and computational studies, focused on the human enzyme PTP1B and the PTP from Yersinia pestis, YopH, suggested that loop conformational dynamics are important in regulating both catalysis and evolvability. We have generated a chimeric protein in which the WPD-loop of YopH is transposed into PTP1B, and eight chimeras that systematically restored the loop sequence back to native PTP1B. Of these, four chimeras were soluble and were subjected to detailed biochemical and structural characterization, and a computational analysis of their WPD-loop dynamics. The chimeras maintain backbone structural integrity, with somewhat slower rates than either wild-type parent, and show differences in the pH dependency of catalysis, and changes in the effect of Mg²⁺. The chimeric proteins' WPD-loops differ significantly in their relative stability and rigidity. The time required for interconversion, coupled with electrostatic effects revealed by simulations, likely accounts for the activity differences between chimeras, and relative to the native enzymes. Our results further the understanding of connections between enzyme activity and the dynamics of catalytically important groups, particularly the effects of non-catalytic residues on key conformational equilibria.

Insertions and deletions in protein evolution and engineering

Protein evolution or engineering studies are traditionally focused on amino acid substitutions and the way these contribute to fitness. Meanwhile, the insertion and deletion of amino acids is often overlooked, despite being one of the most common sources of genetic variation. Recent methodological advances and successful engineering stories have demonstrated that the time is ripe for greater emphasis on these mutations and their understudied effects. This review highlights the evolutionary importance and biotechnological relevance of insertions and deletions (indels). We provide a comprehensive overview of approaches that can be employed to include indels in random, (semi)-rational or computational protein engineering pipelines. Furthermore, we discuss the tolerance to indels at the structural level, address how domain indels can link the function of unrelated proteins, and feature studies that illustrate the surprising and intriguing potential of frameshift mutations.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

LoopGrafter: a web tool for transplanting dynamical loops for protein engineering

The transplantation of loops between structurally related proteins is a compelling method to improve the activity, specificity and stability of enzymes. However, despite the interest of loop regions in protein engineering, the available methods of loop-based rational protein design are scarce. One particular difficulty related to loop engineering is the unique dynamism that enables them to exert allosteric control over the catalytic function of enzymes. Thus, when engaging in a transplantation effort, such dynamics in the context of protein structure need consideration. A second practical challenge is identifying successful excision points for the transplantation or grafting. Here, we present LoopGrafter (https://loschmidt.chemi.muni.cz/loopgrafter/), a web server that specifically guides in the loop grafting process between structurally related proteins. The server provides a step-by-step interactive procedure in which the user can successively identify loops in the two input proteins, calculate their geometries, assess their similarities and dynamics, and select a number of loops to be transplanted. All possible different chimeric proteins derived from any existing recombination point are calculated, and 3D models for each of them are constructed and energetically evaluated. The obtained results can be interactively visualized in a user-friendly graphical interface and downloaded for detailed structural analyses.

Characterization and in planta validation of a CHI4 chitinase from Cajanus platycarpus (Benth.) Maesen for its efficacy against pod borer, Helicoverpa armigera (Hübner)

BACKGROUND

Pigeonpea, Cajanus cajan is one of the economically important legume food crops and a major source of dietary proteins. Management of pod borer, Helicoverpa armigera has been prominent among crop improvement programs. Lack of resistance sources in the cultivated germplasm and crossing incompatibility with pod borer-resistant wild relatives have prompted biotechnological interventions. Identification and exploitation of genes from pigeonpea wild relatives in host plant resistance towards the pod borer assumes pertinence. Dynamic transcriptome analysis of the wild relative vis a vis cultivated pigeonpea identified a CHI4 chitinase as one of the putative insect resistance genes.

RESULTS

The study presents variations in important amino acids in CHI4 chitinases from C. cajan and its wild relative C. platycarpus. Comparative protein modeling and docking analysis of the two proteins demonstrated differences in substrate binding efficacy of the chitinase from C. platycarpus which resulted in a minimum binding energy of -8.7 kcal mol^-1 . Furthermore, we successfully evaluated the insecticidal activity of the chitinase from C. platycarpus against H. armigera challenge through heterologous expression in tobacco. Molecular characterization of transgenic plants confirmed that their efficacy against H. armigera was a result of the integration of CHI4 from C. platycarpus.

CONCLUSION

Docking analysis demonstrated effective substrate interaction as a possible reason for efficacy against pod borer in the chitinase from C. platycarpus. This was authenticated by successful overexpression and bioefficacy assessment against H. armigera in tobacco. The CHI4 gene from C. platycarpus can be useful in the mitigation of H. armegira in pigeonpea as well as in other crops. © 2021 Society of Chemical Industry.

Active-site loop variations adjust activity and selectivity of the cumene dioxygenase

Active-site loops play essential roles in various catalytically important enzyme properties like activity, selectivity, and substrate scope. However, their high flexibility and diversity makes them challenging to incorporate into rational enzyme engineering strategies. Here, we report the engineering of hot-spots in loops of the cumene dioxygenase from Pseudomonas fluorescens IP01 with high impact on activity, regio- and enantioselectivity. Libraries based on alanine scan, sequence alignments, and deletions along with a novel insertion approach result in up to 16-fold increases in activity and the formation of novel products and enantiomers. CAVER analysis suggests possible increases in the active pocket volume and formation of new active-site tunnels, suggesting additional degrees of freedom of the substrate in the pocket. The combination of identified hot-spots with the Linker In Loop Insertion approach proves to be a valuable addition to future loop engineering approaches for enhanced biocatalysts.

Primary and promiscuous functions coexist during evolutionary innovation through whole protein domain acquisitions

Molecular examples of evolutionary innovation are scarce and generally involve point mutations. Innovation can occur through larger rearrangements, but here experimental data is extremely limited. Integron integrases innovated from double-strand- toward single-strand-DNA recombination through the acquisition of the I2 α-helix. To investigate how this transition was possible, we have evolved integrase IntI1 to what should correspond to an early innovation state by selecting for its ancestral activity. Using synonymous alleles to enlarge sequence space exploration, we have retrieved 13 mutations affecting both I2 and the multimerization domains of IntI1. We circumvented epistasis constraints among them using a combinatorial library that revealed their individual and collective fitness effects. We obtained up to 10⁴-fold increases in ancestral activity with various asymmetrical trade-offs in single-strand-DNA recombination. We show that high levels of primary and promiscuous functions could have initially coexisted following I2 acquisition, paving the way for a gradual evolution toward innovation.

Rosetta-Enabled Structural Prediction of Permissive Loop Insertion Sites in Proteins

While loop motifs frequently play a major role in protein function, our understanding of how to rationally engineer proteins with novel loop domains remains limited. In the absence of rational approaches, the incorporation of loop domains often destabilizes proteins, thereby requiring massive screening and selection to identify sites that can accommodate loop insertion. We developed a computational strategy for rapidly scanning the entire structure of a scaffold protein to determine the impact of loop insertion at all possible amino acid positions. This approach is based on the Rosetta kinematic loop modeling protocol and was demonstrated by identifying sites in lipase that were permissive to insertion of the LAP peptide. Interestingly, the identification of permissive sites was dependent on the contribution of the residues in the near-loop environment on the Rosetta score and did not correlate with conventional structural features (e.g., B-factors). As evidence of this, several insertion sites (e.g., following residues 17, 47-49, and 108), which were predicted and confirmed to be permissive, interrupted helices, while others (e.g., following residues 43, 67, 116, 119, and 121), which are situated in loop regions, were nonpermissive. This approach was further shown to be predictive for β-glucosidase and human phosphatase and tensin homologue (PTEN), and to facilitate the engineering of insertion sites through in silico mutagenesis. By enabling the design of loop-containing protein libraries with high probabilities of soluble expression, this approach has broad implications in many areas of protein engineering, including antibody design, improving enzyme activity, and protein modification.

Remodeling enzyme active sites by stepwise loop insertion

The remolding active site loops via residue insertion/deletion as well as substitution is thought to play a key role in enzyme divergent evolution. However, enzyme engineering by residue insertion in active site loops often severely perturbs the protein structural integrity and causes protein misfolding and activity loss. We have designed a stepwise loop insertion strategy (StLois), in which a pair of randomized residues is introduced in a stepwise manner, efficiently collating mutational fitness effects. The strategy of StLois constitutes three key steps. First, the target regions should be identified through structural and functional analysis on the counterpart enzymes. Second, pair residues can be introduced in loop regions through insertion with NNK codon degeneracy. Third, the best hit used as a template for the next round mutagenesis. The residue insertion process can repeat as many times as necessary. By using the StLois method, we have evolved the substrate preference of a lactonase to phosphotriesterase. In this chapter, we describe the detailed StLois technique, which efficiently expands the residue in the loop region and remolds the architecture of enzyme active site for novel catalytic properties.

Engineering the specificity of Streptococcus pyogenes sortase A by loop grafting

Sortases are a group of enzymes displayed on the cell‐wall of Gram‐positive bacteria. They are responsible for the attachment of virulence factors onto the peptidoglycan in a transpeptidation reaction through recognition of a pentapeptide substrate. Most housekeeping sortases recognize one specific pentapeptide motif; however, Streptococcus pyogenes sortase A (SpSrtA WT) recognizes LPETG, LPETA and LPKLG motifs. Here, we examined SpSrtA's flexible substrate specificity by investigating the role of the β7/β8 loop in determining substrate specificity. We exchanged the β7/β8 loop in SpSrtA with corresponding β7/β8 loops from Staphylococcus aureus (SaSrtA WT) and Bacillus anthracis (BaSrtA WT). While the BaSrtA‐derived variant showed no enzymatic activity toward either LPETG or LPETA substrates, the activity of the SaSrtA‐derived mutant toward the LPETA substrate was completely abolished. Instead, the mutant had an improved activity toward LPETG, the preferred substrate of SaSrtA WT.

Structure of the plastic-degrading Ideonella sakaiensis MHETase bound to a substrate

The extreme durability of polyethylene terephthalate (PET) debris has rendered it a long-term environmental burden. At the same time, current recycling efforts still lack sustainability. Two recently discovered bacterial enzymes that specifically degrade PET represent a promising solution. First, Ideonella sakaiensis PETase, a structurally well-characterized consensus α/β-hydrolase fold enzyme, converts PET to mono-(2-hydroxyethyl) terephthalate (MHET). MHETase, the second key enzyme, hydrolyzes MHET to the PET educts terephthalate and ethylene glycol. Here, we report the crystal structures of active ligand-free MHETase and MHETase bound to a nonhydrolyzable MHET analog. MHETase, which is reminiscent of feruloyl esterases, possesses a classic α/β-hydrolase domain and a lid domain conferring substrate specificity. In the light of structure-based mapping of the active site, activity assays, mutagenesis studies and a first structure-guided alteration of substrate specificity towards bis-(2-hydroxyethyl) terephthalate (BHET) reported here, we anticipate MHETase to be a valuable resource to further advance enzymatic plastic degradation.

Highly active enzymes by automated combinatorial backbone assembly and sequence design

Automated design of enzymes with wild-type-like catalytic properties has been a long-standing but elusive goal. Here, we present a general, automated method for enzyme design through combinatorial backbone assembly. Starting from a set of homologous yet structurally diverse enzyme structures, the method assembles new backbone combinations and uses Rosetta to optimize the amino acid sequence, while conserving key catalytic residues. We apply this method to two unrelated enzyme families with TIM-barrel folds, glycoside hydrolase 10 (GH10) xylanases and phosphotriesterase-like lactonases (PLLs), designing 43 and 34 proteins, respectively. Twenty-one GH10 and seven PLL designs are active, including designs derived from templates with <25% sequence identity. Moreover, four designs are as active as natural enzymes in these families. Atomic accuracy in a high-activity GH10 design is further confirmed by crystallographic analysis. Thus, combinatorial-backbone assembly and design may be used to generate stable, active, and structurally diverse enzymes with altered selectivity or activity.

Computationally designed enzymes often show lower activity or stability than their natural counterparts. Here, the authors present an evolution-inspired method for automated enzyme design, creating stable enzymes with accurate active site architectures and wild-type-like activities.

Stepwise Loop Insertion Strategy for Active Site Remodeling to Generate Novel Enzyme Functions

The remodeling of active sites to generate novel biocatalysts is an attractive and challenging task. We developed a stepwise loop insertion strategy (StLois), in which randomized residue pairs are inserted into active site loops. The phosphotriesterase-like lactonase from Geobacillus kaustophilus (GkaP-PLL) was used to investigate StLois's potential for changing enzyme function. By inserting six residues into active site loop 7, the best variant ML7-B6 demonstrated a 16-fold further increase in catalytic efficiency toward ethyl-paraoxon compared with its initial template, that is a 609-fold higher, >10⁷ fold substrate specificity shift relative to that of wild-type lactonase. The remodeled variants displayed 760-fold greater organophosphate hydrolysis activity toward the organophosphates parathion, diazinon, and chlorpyrifos. Structure and docking computations support the source of notably inverted enzyme specificity. Considering the fundamental importance of active site loops, the strategy has potential for the rapid generation of novel enzyme functions by loop remodeling.

LoopX: A Graphical User Interface-Based Database for Comprehensive Analysis and Comparative Evaluation of Loops from Protein Structures

Due to their crucial role in function, folding, and stability, protein loops are being targeted for grafting/designing to create novel or alter existing functionality and improve stability and foldability. With a view to facilitate a thorough analysis and effectual search options for extracting and comparing loops for sequence and structural compatibility, we developed, LoopX a comprehensively compiled library of sequence and conformational features of ∼700,000 loops from protein structures. The database equipped with a graphical user interface is empowered with diverse query tools and search algorithms, with various rendering options to visualize the sequence- and structural-level information along with hydrogen bonding patterns, backbone φ, ψ dihedral angles of both the target and candidate loops. Two new features (i) conservation of the polar/nonpolar environment and (ii) conservation of sequence and conformation of specific residues within the loops have also been incorporated in the search and retrieval of compatible loops for a chosen target loop. Thus, the LoopX server not only serves as a database and visualization tool for sequence and structural analysis of protein loops but also aids in extracting and comparing candidate loops for a given target loop based on user-defined search options.

Multitarget selection of catalytic antibodies with β-lactamase activity using phage display

β-lactamase enzymes responsible for bacterial resistance to antibiotics are among the most important health threats to the human population today. Understanding the increasingly vast structural motifs responsible for the catalytic mechanism of β-lactamases will help improve the future design of new generation antibiotics and mechanism-based inhibitors of these enzymes. Here we report the construction of a large murine single chain fragment variable (scFv) phage display library of size 2.7 × 10⁹ with extended diversity by combining different mouse models. We have used two molecularly different inhibitors of the R-TEM β-lactamase as targets for selection of catalytic antibodies with β-lactamase activity. This novel methodology has led to the isolation of five antibody fragments, which are all capable of hydrolyzing the β-lactam ring. Structural modeling of the selected scFv has revealed the presence of different motifs in each of the antibody fragments potentially responsible for their catalytic activity. Our results confirm (a) the validity of using our two target inhibitors for the in vitro selection of catalytic antibodies endowed with β-lactamase activity, and (b) the plasticity of the β-lactamase active site responsible for the wide resistance of these enzymes to clinically available inhibitors and antibiotics.

Biochemical and Structural Analysis of a Novel Esterase from Caulobacter crescentus related to Penicillin-Binding Protein (PBP)

Considering that the prevalence of antibiotic-resistant pathogenic bacteria is largely increasing, a thorough understanding of penicillin-binding proteins (PBPs) is of great importance and crucial significance because this enzyme family is a main target of β-lactam-based antibiotics. In this work, combining biochemical and structural analysis, we present new findings that provide novel insights into PBPs. Here, a novel PBP homologue (CcEstA) from Caulobacter crescentus CB15 was characterized using native-PAGE, mass spectrometry, gel filtration, CD spectroscopy, fluorescence, reaction kinetics, and enzyme assays toward various substrates including nitrocefin. Furthermore, the crystal structure of CcEstA was determined at a 1.9 Å resolution. Structural analyses showed that CcEstA has two domains: a large α/β domain and a small α-helix domain. A nucleophilic serine (Ser⁶⁸) residue is located in a hydrophobic groove between the two domains along with other catalytic residues (Lys⁷¹ and Try¹⁵⁷). Two large flexible loops (UL and LL) of CcEstA are proposed to be involved in the binding of incoming substrates. In conclusion, CcEstA could be described as a paralog of the group that contains PBPs and β-lactamases. Therefore, this study could provide new structural and functional insights into the understanding this protein family.

High adaptability of the omega loop underlies the substrate-spectrum-extension evolution of a class A β-lactamase, PenL

The omega loop in β-lactamases plays a pivotal role in substrate recognition and catalysis, and some mutations in this loop affect the adaptability of the enzymes to new antibiotics. Various mutations, including substitutions, deletions, and intragenic duplications resulting in tandem repeats (TRs), have been associated with β-lactamase substrate spectrum extension. TRs are unique among the mutations as they cause severe structural perturbations in the enzymes. We explored the process by which TRs are accommodated in order to test the adaptability of the omega loop. Structures of the mutant enzymes showed that the extra amino acid residues in the omega loop were freed outward from the enzyme, thereby maintaining the overall enzyme integrity. This structural adjustment was accompanied by disruptions of the internal α-helix and hydrogen bonds that originally maintained the conformation of the omega loop and the active site. Consequently, the mutant enzymes had a relaxed binding cavity, allowing for access of new substrates, which regrouped upon substrate binding in an induced-fit manner for subsequent hydrolytic reactions. Together, the data demonstrate that the design of the binding cavity, including the omega loop with its enormous adaptive capacity, is the foundation of the continuous evolution of β-lactamases against new drugs.

Why reinvent the wheel? Building new proteins based on ready-made parts

We protein engineers are ambivalent about evolution: on the one hand, evolution inspires us with myriad examples of biomolecular binders, sensors, and catalysts; on the other hand, these examples are seldom well-adapted to the engineering tasks we have in mind. Protein engineers have therefore modified natural proteins by point substitutions and fragment exchanges in an effort to generate new functions. A counterpoint to such design efforts, which is being pursued now with greater success, is to completely eschew the starting materials provided by nature and to design new protein functions from scratch by using de novo molecular modeling and design. While important progress has been made in both directions, some areas of protein design are still beyond reach. To this end, we advocate a synthesis of these two strategies: by using design calculations to both recombine and optimize fragments from natural proteins, we can build stable and as of yet un-sampled structures, thereby granting access to an expanded repertoire of conformations and desired functions. We propose that future methods that combine phylogenetic analysis, structure and sequence bioinformatics, and atomistic modeling may well succeed where any one of these approaches has failed on its own.

The Evolution of New Catalytic Mechanisms for Xenobiotic Hydrolysis in Bacterial Metalloenzymes

An increasing number of bacterial metalloenzymes have been shown to catalyse the breakdown of xenobiotics in the environment, while others exhibit a variety of promiscuous xenobiotic-degrading activities. Several different evolutionary processes have allowed these enzymes to gain or enhance xenobiotic-degrading activity. In this review, we have surveyed the range of xenobiotic-degrading metalloenzymes, and discuss the molecular and catalytic basis for the development of new activities. We also highlight how our increased understanding of the natural evolution of xenobiotic-degrading metalloenzymes can be been applied to laboratory enzyme design.

A designed supramolecular protein assembly with in vivo enzymatic activity

The generation of new enzymatic activities has mainly relied on repurposing the interiors of preexisting protein folds because of the challenge in designing functional, three-dimensional protein structures from first principles. Here we report an artificial metallo-β-lactamase, constructed via the self-assembly of a structurally and functionally unrelated, monomeric redox protein into a tetrameric assembly that possesses catalytic zinc sites in its interfaces. The designed metallo-β-lactamase is functional in the Escherichia coli periplasm and enables the bacteria to survive treatment with ampicillin. In vivo screening of libraries has yielded a variant that displays a catalytic proficiency [(k(cat)/K(m))/k(uncat)] for ampicillin hydrolysis of 2.3 × 10(6) and features the emergence of a highly mobile loop near the active site, a key component of natural β-lactamases to enable substrate interactions.

Structural basis of substrate selectivity of E. coli prolidase

Prolidases, metalloproteases that catalyze the cleavage of Xaa-Pro dipeptides, are conserved enzymes found in prokaryotes and eukaryotes. In humans, prolidase is crucial for the recycling of collagen. To further characterize the essential elements of this enzyme, we utilized the Escherichia coli prolidase, PepQ, which shares striking similarity with eukaryotic prolidases. Through structural and bioinformatic insights, we have extended previous characterizations of the prolidase active site, uncovering a key component for substrate specificity. Here we report the structure of E. coli PepQ, solved at 2.0 Å resolution. The structure shows an antiparallel, dimeric protein, with each subunit containing N-terminal and C-terminal domains. The C-terminal domain is formed by the pita-bread fold typical for this family of metalloproteases, with two Mg(II) ions coordinated by five amino-acid ligands. Comparison of the E. coli PepQ structure and sequence with homologous structures and sequences from a diversity of organisms reveals distinctions between prolidases from Gram-positive eubacteria and archaea, and those from Gram-negative eubacteria, including the presence of loop regions in the E. coli protein that are conserved in eukaryotes. One such loop contains a completely conserved arginine near the catalytic site. This conserved arginine is predicted by docking simulations to interact with the C-terminus of the substrate dipeptide. Kinetic analysis using both a charge-neutralized substrate and a charge-reversed variant of PepQ support this conclusion, and allow for the designation of a new role for this key region of the enzyme active site.

Highly diverse protein library based on the ubiquitous (β/α)₈ enzyme fold yields well-structured proteins through in vitro folding selection

Proper protein folding is a prerequisite for protein stability and enzymatic activity. Although directed evolution can be a powerful tool to investigate enzymatic function and to isolate novel activities, well-designed libraries of folded proteins are essential. In vitro selection methods are particularly capable of searching for enzymatic activities in libraries of trillions of protein variants, yet high-quality libraries of well-folded enzymes with such high diversity are lacking. We describe the construction and detailed characterization of a folding-enriched protein library based on the ubiquitous (β/α)₈ barrel fold, which is found in five of the six enzyme classes. We introduced seven randomized loops on the catalytic face of the monomeric, thermostable (β/α)₈ barrel of glycerophosphodiester phosphodiesterase (GDPD) from Thermotoga maritima. We employed in vitro folding selection based on protease digestion to enrich intermediate libraries containing three to four randomized loops for folded variants, and then combined them to assemble the final library (10¹⁴ DNA sequences). The resulting library was analyzed by using the in vitro protease assay and an in vivo GFP-folding assay; it contains ∼10¹² soluble monomeric protein variants. We isolated six library members and demonstrated that these proteins are soluble, monomeric and show (β/α)₈-barrel fold-like secondary and tertiary structure. The quality of the folding-enriched library improved up to 50-fold compared to a control library that was assembled without the folding selection. To the best of our knowledge, this work is the first example of combining the ultra-high throughput mRNA display method with selection for folding. The resulting (β/α)₈ barrel libraries provide a valuable starting point to study the unique catalytic capabilities of the (β/α)₈ fold, and to isolate novel enzymes.

Reconstructing a missing link in the evolution of a recently diverged phosphotriesterase by active-site loop remodeling

Only decades after the introduction of organophosphate pesticides, bacterial phosphotriesterases (PTEs) have evolved to catalyze their degradation with remarkable efficiency. Their closest known relatives, lactonases, with promiscuous phosphotriasterase activity, dubbed PTE-like lactonases (PLLs), share only 30% sequence identity and also differ in the configuration of their active-site loops. PTE was therefore presumed to have evolved from a yet unknown PLL whose primary activity was the hydrolysis of quorum sensing homoserine lactones (HSLs) (Afriat et al. (2006) Biochemistry 45, 13677-13686). However, how PTEs diverged from this presumed PLL remains a mystery. In this study we investigated loop remodeling as a means of reconstructing a homoserine lactonase ancestor that relates to PTE by few mutational steps. Although, in nature, loop remodeling is a common mechanism of divergence of enzymatic functions, reproducing this process in the laboratory is a challenge. Structural and phylogenetic analyses enabled us to remodel one of PTE's active-site loops into a PLL-like configuration. A deletion in loop 7, combined with an adjacent, highly epistatic, point mutation led to the emergence of an HSLase activity that is undetectable in PTE (k(cat)/K(M) values of up to 2 × 10(4)). The appearance of the HSLase activity was accompanied by only a minor decrease in PTE's paraoxonase activity. This specificity change demonstrates the potential role of bifunctional intermediates in the divergence of new enzymatic functions and highlights the critical contribution of loop remodeling to the rapid divergence of new enzyme functions.

Establishing a toolkit for precursor-directed polyketide biosynthesis: exploring substrate promiscuities of acid-CoA ligases

Polyketides are chemically diverse and medicinally important biochemicals that are biosynthesized from acyl-CoA precursors by polyketide synthases. One of the limitations to combinatorial biosynthesis of polyketides has been the lack of a toolkit that describes the means of delivering novel acyl-CoA precursors necessary for polyketide biosynthesis. Using five acid-CoA ligases obtained from various plants and microorganisms, we biosynthesized an initial library of 79 acyl-CoA thioesters by screening each of the acid-CoA ligases against a library of 123 carboxylic acids. The library of acyl-CoA thioesters includes derivatives of cinnamyl-CoA, 3-phenylpropanoyl-CoA, benzoyl-CoA, phenylacetyl-CoA, malonyl-CoA, saturated and unsaturated aliphatic CoA thioesters, and bicyclic aromatic CoA thioesters. In our search for the biosynthetic routes of novel acyl-CoA precursors, we discovered two previously unreported malonyl-CoA derivatives (3-thiophenemalonyl-CoA and phenylmalonyl-CoA) that cannot be produced by canonical malonyl-CoA synthetases. This report highlights the utility and importance of determining substrate promiscuities beyond conventional substrate pools and describes novel enzymatic routes for the establishment of precursor-directed combinatorial polyketide biosynthesis.

Exploring the evolution of novel enzyme functions within structurally defined protein superfamilies

In order to understand the evolution of enzyme reactions and to gain an overview of biological catalysis we have combined sequence and structural data to generate phylogenetic trees in an analysis of 276 structurally defined enzyme superfamilies, and used these to study how enzyme functions have evolved. We describe in detail the analysis of two superfamilies to illustrate different paradigms of enzyme evolution. Gathering together data from all the superfamilies supports and develops the observation that they have all evolved to act on a diverse set of substrates, whilst the evolution of new chemistry is much less common. Despite that, by bringing together so much data, we can provide a comprehensive overview of the most common and rare types of changes in function. Our analysis demonstrates on a larger scale than previously studied, that modifications in overall chemistry still occur, with all possible changes at the primary level of the Enzyme Commission (E.C.) classification observed to a greater or lesser extent. The phylogenetic trees map out the evolutionary route taken within a superfamily, as well as all the possible changes within a superfamily. This has been used to generate a matrix of observed exchanges from one enzyme function to another, revealing the scale and nature of enzyme evolution and that some types of exchanges between and within E.C. classes are more prevalent than others. Surprisingly a large proportion (71%) of all known enzyme functions are performed by this relatively small set of 276 superfamilies. This reinforces the hypothesis that relatively few ancient enzymatic domain superfamilies were progenitors for most of the chemistry required for life.

Enzymes, as biological catalysts, are crucial to life. Understanding how enzymes have evolved to perform the wide variety of reactions found across all kingdoms of life is fundamental to a broad range of biological studies, especially those leading to new therapeutics. To unravel the evolution of novel enzyme function requires combining information on protein structure, sequence, phylogeny and chemistry (in terms of interacting small molecules and reaction mechanisms). We have developed a protocol for integrating this wide range of data, which we have applied to a relatively large number of families comprising some very diverse relatives. This has permitted us to present an initial overview of the evolution of novel enzyme functions, in which we observe that some changes in function between relatives are more common than others, with most of the functionality observed in nature confined to relatively few families. Moreover, we are able to identify the evolutionary route taken within a superfamily to change the enzyme function from one reaction to another. This information may help in predicting the function of an enzyme that has yet to be experimentally characterised as well as in designing new enzymes for industrial and medical purposes.

Evolution to carbapenem-hydrolyzing activity in noncarbapenemase class D β-lactamase OXA-10 by rational protein design

Class D β-lactamases with carbapenemase activity are emerging as carbapenem-resistance determinants in gram-negative bacterial pathogens, mostly Acinetobacter baumannii and Klebsiella pneumoniae. Carbapenemase activity is an unusual feature among class D β-lactamases, and the structural elements responsible for this activity remain unclear. Based on structural and molecular dynamics data, we previously hypothesized a potential role of the residues located in the short-loop connecting strands β5 and β6 (the β5-β6 loop) in conferring the carbapenemase activity of the OXA-48 enzyme. In this work, the narrow-spectrum OXA-10 class D β-lactamase, which is unable to hydrolyze carbapenems, was used as a model to investigate the possibility of evolving carbapenemase activity by replacement of the β5-β6 loop with those present in three different lineages of class D carbapenemases (OXA-23, OXA-24, and OXA-48). Biological assays and kinetic measurements showed that all three OXA-10-derived hybrids acquired significant carbapenemase activity. Structural analysis of the OXA-10loop24 and OXA-10loop48 hybrids revealed no significant changes in the molecular fold of the enzyme, except for the orientation of the substituted β5-β6 loops, which was reminiscent of that found in their parental enzymes. These results demonstrate the crucial role of the β5-β6 loop in the carbapenemase activity of class D β-lactamases, and provide previously unexplored insights into the mechanism by which these enzymes can evolve carbapenemase activity.

Directed evolution of a thermostable quorum-quenching lactonase from the amidohydrolase superfamily

A thermostable quorum-quenching lactonase from Geobacillus kaustophilus HTA426 (GI: 56420041) was used as an initial template for in vitro directed evolution experiments. This enzyme belongs to the phosphotriesterase-like lactonase (PLL) group of enzymes within the amidohydrolase superfamily that hydrolyze N-acylhomoserine lactones (AHLs) that are involved in virulence pathways of quorum-sensing pathogenic bacteria. Here we have determined the N-butyryl-L-homoserine lactone-liganded structure of the catalytically inactive D266N mutant of this enzyme to a resolution of 1.6 Å. Using a tunable, bioluminescence-based quorum-quenching molecular circuit, the catalytic efficiency was enhanced, and the AHL substrate range increased through two point mutations on the loops at the C-terminal ends of the third and seventh β-strands. This E101N/R230I mutant had an increased value of k(cat)/K(m) of 72-fold toward 3-oxo-N-dodecanoyl-L-homoserine lactone. The evolved mutant also exhibited lactonase activity toward N-butyryl-L-homoserine lactone, an AHL that was previously not hydrolyzed by the wild-type enzyme. Both the purified wild-type and mutant enzymes contain a mixture of zinc and iron and are colored purple and brown, respectively, at high concentrations. The origin of this coloration is suggested to be because of a charge transfer complex involving the β-cation and Tyr-99 within the enzyme active site. Modulation of the charge transfer complex alters the lactonase activity of the mutant enzymes and is reflected in enzyme coloration changes. We attribute the observed enhancement in catalytic reactivity of the evolved enzyme to favorable modulations of the active site architecture toward productive geometries required for chemical catalysis.

The origin and evolution of lignin biosynthesis

Lignin, a phenolic polymer derived mainly from hydroxycinnamyl alcohols, is ubiquitously present in tracheophytes. The development of lignin biosynthesis has been considered to be one of the key factors that allowed land plants to flourish in terrestrial ecosystems. Lignin provides structural rigidity for tracheophytes to stand upright, and strengthens the cell wall of their water-conducting tracheary elements to withstand the negative pressure generated during transpiration. In this review, we discuss a number of aspects regarding the origin and evolution of lignin biosynthesis during land plant evolution, including the establishment of its monomer biosynthetic scaffold, potential precursors to the lignin polymer, as well as the emergence of the polymerization machinery and regulatory system. The accumulated knowledge on the topic, as summarized here, provides us with an evolutionary view on how this complex metabolic system emerged and developed.

A portable hot spot recognition loop transfers sequence preferences from APOBEC family members to activation-induced cytidine deaminase

Enzymes of the AID/APOBEC family, characterized by the targeted deamination of cytosine to generate uracil within DNA, mediate numerous critical immune responses. One family member, activation-induced cytidine deaminase (AID), selectively introduces uracil into antibody variable and switch regions, promoting antibody diversity through somatic hypermutation or class switching. Other family members, including APOBEC3F and APOBEC3G, play an important role in retroviral defense by acting on viral reverse transcripts. These enzymes are distinguished from one another by targeting cytosine within different DNA sequence contexts; however, the reason for these differences is not known. Here, we report the identification of a recognition loop of 9-11 amino acids that contributes significantly to the distinct sequence motifs of individual family members. When this recognition loop is grafted from the donor APOBEC3F or 3G proteins into the acceptor scaffold of AID, the mutational signature of AID changes toward that of the donor proteins. These loop-graft mutants of AID provide useful tools for dissecting the biological impact of DNA sequence preferences upon generation of antibody diversity, and the results have implications for the evolution and specialization of the AID/APOBEC family.

Directed evolution of a quorum-quenching lactonase from Mycobacterium avium subsp. paratuberculosis K-10 in the amidohydrolase superfamily

The PLL(PTE-like lactonase)-group of enzymes within the amidohydrolase superfamily hydrolyze N-acyl-homoserine lactones (AHLs) that are involved in bacterial quorum-sensing pathways. These enzymes possess the (beta/alpha)(8)-barrel fold and serve as attractive templates for in vitro evolution and engineering of quorum-quenching biological molecules that can serve as antivirulence therapeutic agents. Using a quorum-quenching lactonase from Mycobacterium avium subsp. paratuberculosis K-10 (GI: 41409766) as the initial template for in vitro evolution experiments, we enhanced the catalytic efficiency and increased the substrate range of the wild-type enzyme through a single point mutation on the loop at the C-terminal end of the eighth beta-strand. This N266Y mutant had an increased value of k(cat)/K(M) of 30- and 32-fold toward 3-oxo-N-octanoyl-l-homoserine lactone and N-hexanoyl-l-homoserine lactone, respectively; the evolved mutant also exhibited lactonase activity toward 3-oxo-N-hexanoyl-l-homoserine lactone and N-butyryl-l-homoserine lactone, AHLs that were previously not hydrolyzed by the wild-type enzyme. This article reinforces the evolutionary potential of the (beta/alpha)(8)-barrel fold and highlights the possibility of using quorum-quenching lactonases in the amidohydrolase superfamily as templates for engineering biomolecules of therapeutic use.

Enzyme (re)design: lessons from natural evolution and computation

The (re)design of enzymes to catalyze 'new' reactions is a topic of considerable practical and intellectual interest. Directed evolution (random mutagenesis followed by screening/selection) has been used widely to identify novel biocatalysts. However, 'rational' approaches using either natural divergent evolution or computational predictions based on chemical principles have been less successful. This review summarizes recent progress in evolution-based and computation-based (re)design.

Protein design through systematic catalytic loop exchange in the (beta/alpha)8 fold

Protein engineering by directed evolution has proven effective in achieving various functional modifications, but the well-established protocols for the introduction of variability, typically limited to random point mutations, seriously restrict the scope of the approach. In an attempt to overcome this limitation, we sought to explore variant libraries with richer diversity at regions recognized as functionally important through an exchange of natural components, thus combining design with combinatorial diversity. With this approach, we expected to maintain interactions important for protein stability while directing the introduction of variability to areas important for catalysis. Our strategy consisted in loop exchange over a (beta/alpha)(8) fold. Phosphoribosylanthranilate isomerase was chosen as scaffold, and we investigated its tolerance to loop exchange by fusing variant libraries to the chloramphenicol acetyl transferase coding gene as an in vivo folding reporter. We replaced loops 2, 4, and 6 of phosphoribosylanthranilate isomerase with loops of varied types and sizes from enzymes sharing the same fold. To allow for a better structural fit, saturation mutagenesis was adopted at two amino acid positions preceding the exchanged loop. Our results showed that 30% to 90% of the generated mutants in the different libraries were folded. Some variants were selected for further characterization after removal of chloramphenicol acetyl transferase gene, and their stability was studied by circular dichroism and fluorescence spectroscopy. The sequences of 545 clones show that the introduction of variability at "hinges" connecting the loops with the scaffold exhibited a noticeable effect on the appearance of folded proteins. Also, we observed that each position accepted foreign loops of different sizes and sequences. We believe our work provides the basis of a general method of exchanging variably sized loops within the (beta/alpha)(8) fold, affording a novel starting point for the screening of novel activities as well as modest diversions from an original activity.

How protein stability and new functions trade off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (ΔΔG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (ΔΔG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied ΔΔG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to ΔΔG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average ΔΔG = +0.9 kcal/mol), and are almost as destabilizing as the “average” mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently “silent” mutations in regions of the protein that are irrelevant to its function.

To perform its function, a protein must fold into a complex, three-dimensional structure that is maintained by a network of interactions between its amino acid residues. Evolution of a new protein function will be driven by mutation of amino acids in key positions (new-function mutations). Such mutation can also hamper interactions that ensure the stability of a protein's fold—sometimes to a degree that renders the protein non-functional. Indeed, previous studies have noted that the evolution of new enzymatic functions is accompanied by significant losses in protein stability, suggesting a “tradeoff” between acquisition of new enzymatic functions and stability. But since most mutations are destabilizing, we sought to compare new-function mutations with other types of mutations. We performed a comprehensive analysis of the type, location, and stability effects of mutations that have conferred new enzymatic functions in laboratory evolution experiments. We found that stability changes (ΔΔG) of new-function mutations are similar to those of all other mutations, but are weaker than those of mutations that characterize neutral evolutionary changes (mutations that accumulate with no change of structure and function). Our analysis also revealed the important role of neutral (i.e., “non-functional”) mutations in compensating for the destabilizing effects of the “new-function” mutations.

Enzyme optimization: moving from blind evolution to statistical exploration of sequence-function space

Directed evolution is a powerful tool for the creation of commercially useful enzymes, particularly those approaches that are based on in vitro recombination methods, such as DNA shuffling. Although these types of search algorithms are extraordinarily efficient compared with purely random methods, they do not explicitly represent or interrogate the genotype-phenotype relationship and are essentially blind in nature. Recently, however, researchers have begun to apply multivariate statistical techniques to model protein sequence-function relationships and guide the evolutionary process by rapidly identifying beneficial diversity for recombination. In conjunction with state-of-the-art library generation methods, the statistical approach to sequence optimization is now being used routinely to create enzymes efficiently for industrial applications.