Directed evolution of (βα)8-barrel enzymes

Study of peripheral domains in structure-function of isocitrate lyase (ICL) from Pseudomonas aeruginosa

The capacity of Pseudomonas aeruginosa to assimilate nutrients is essential for niche colonization and contributes to its pathogenicity. Isocitrate lyase (ICL), the first enzyme of the glyoxylate cycle, redirects isocitrate from the tricarboxylic acid cycle to render glyoxylate and succinate. P. aeruginosa ICL (PaICL) is regarded as a virulence factor due to its role in carbon assimilation during infection. The AceA/ICL protein family shares the catalytic domain I, triosephosphate isomerase barrel (TIM-barrel). The carboxyl terminus of domain I is essential for Escherichia coli ICL (EcICL) of subfamily 1. PaICL, which belongs to subfamily 3, has domain II inserted at the periphery of domain I, which is believed to participate in enzyme oligomerization. In addition, PaICL has the α13-loop-α14 (extended motif), which protrudes from the enzyme core, being of unknown function. This study investigates the role of domain II, the extended motif, and the carboxyl-terminus (C-ICL) and amino-terminus (N-ICL) regions in the function of the PaICL enzyme, also as their involvement in the virulence of P. aeruginosa PAO1. Deletion of domain II and the extended motif results in enzyme inactivation and structural instability of the enzyme. The His6-tag fusion at the C-ICL protein produced a less efficient enzyme than fusion at the N-ICL, but without affecting the acetate assimilation or virulence. The PaICL homotetrameric structure of the enzyme was more stable in the N-His6-ICL than in the C-His6-ICL, suggesting that the C-terminus is critical for the ICL quaternary conformation. The ICL-mutant A39 complemented with the recombinant proteins N-His6-ICL or C-His6-ICL were more virulent than the WT PAO1 strain. The findings indicate that the domain II and the extended motif are essential for the ICL structure/function, and the C-terminus is involved in its quaternary structure conformation, confirming that in P. aeruginosa, the ICL is essential for acetate assimilation and virulence.

Active site center redesign increases protein stability preserving catalysis in thioredoxin

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

An allosteric pathway explains beneficial fitness in yeast for long-range mutations in an essential TIM barrel enzyme

Protein evolution proceeds by a complex response of organismal fitness to mutations that can simultaneously affect protein stability, structure, and enzymatic activity. To probe the relationship between genotype and phenotype, we chose a fundamental paradigm for protein evolution, folding, and design, the (βα)₈ TIM barrel fold. Here, we demonstrate the role of long-range allosteric interactions in the adaptation of an essential hyperthermophilic TIM barrel enzyme to mesophilic conditions in a yeast host. Beneficial fitness effects observed with single and double mutations of the canonical βα-hairpin clamps and the α-helical shell distal to the active site revealed an underlying energy network between opposite faces of the cylindrical β-barrel. We experimentally determined the fitness of multiple mutants in the energetic phase plane, contrasting the energy barrier of the chemical reaction and the folding free energy of the protein. For the system studied, the reaction energy barrier was the primary determinant of organism fitness. Our observations of long-range epistatic interactions uncovered an allosteric pathway in an ancient and ubiquitous enzyme that may provide a novel way of designing proteins with a desired activity and stability profile.

Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 β-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the β-barrel core is highly rigid, while the α-helices that flank the β-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.

De novo design of a four-fold symmetric TIM-barrel protein with atomic-level accuracy

Despite efforts for over 25 years, de novo protein design has not succeeded in achieving the TIM-barrel fold. Here we describe the computational design of four-fold symmetrical (β/α)8 barrels guided by geometrical and chemical principles. Experimental characterization of 33 designs revealed the importance of side chain-backbone hydrogen bonds for defining the strand register between repeat units. The X-ray crystal structure of a designed thermostable 184-residue protein is nearly identical to that of the designed TIM-barrel model. PSI-BLAST searches do not identify sequence similarities to known TIM-barrel proteins, and sensitive profile-profile searches indicate that the design sequence is distant from other naturally occurring TIM-barrel superfamilies, suggesting that Nature has sampled only a subset of the sequence space available to the TIM-barrel fold. The ability to design TIM barrels de novo opens new possibilities for custom-made enzymes.

Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently

The amino acid sequence of a protein affects both its structure and its function. Thus, the ability to modify the sequence, and hence the structure and activity, of individual proteins in a systematic way, opens up many opportunities, both scientifically and (as we focus on here) for exploitation in biocatalysis. Modern methods of synthetic biology, whereby increasingly large sequences of DNA can be synthesised de novo, allow an unprecedented ability to engineer proteins with novel functions. However, the number of possible proteins is far too large to test individually, so we need means for navigating the 'search space' of possible protein sequences efficiently and reliably in order to find desirable activities and other properties. Enzymologists distinguish binding (Kd) and catalytic (kcat) steps. In a similar way, judicious strategies have blended design (for binding, specificity and active site modelling) with the more empirical methods of classical directed evolution (DE) for improving kcat (where natural evolution rarely seeks the highest values), especially with regard to residues distant from the active site and where the functional linkages underpinning enzyme dynamics are both unknown and hard to predict. Epistasis (where the 'best' amino acid at one site depends on that or those at others) is a notable feature of directed evolution. The aim of this review is to highlight some of the approaches that are being developed to allow us to use directed evolution to improve enzyme properties, often dramatically. We note that directed evolution differs in a number of ways from natural evolution, including in particular the available mechanisms and the likely selection pressures. Thus, we stress the opportunities afforded by techniques that enable one to map sequence to (structure and) activity in silico, as an effective means of modelling and exploring protein landscapes. Because known landscapes may be assessed and reasoned about as a whole, simultaneously, this offers opportunities for protein improvement not readily available to natural evolution on rapid timescales. Intelligent landscape navigation, informed by sequence-activity relationships and coupled to the emerging methods of synthetic biology, offers scope for the development of novel biocatalysts that are both highly active and robust.

Three amino acid changes contribute markedly to the thermostability of β-glucosidase BglC from Thermobifida fusca

Thermostability of β-glucosidase was enhanced by family shuffling, site saturation mutagenesis, and site-directed mutagenesis. Family shuffling was carried out based on β-glucosidase BglC from Thermobifida fusca and β-glucosidase BglB from Paebibacillus polymxyxa with the help of synthetic primers. High-throughput screening revealed mutants with higher thermostability than both parental enzymes. The most thermostable mutant VM2 containing three key amino acid changes in L444Y/G447S/A433V had a 144-fold increase in half-life of inactivation as compared to the parental enzyme BglC. The mutant VM2 showed 28% and 94% increase in k(cat) towards p-nitrophenyl-β-D-glucopyranoside (pNPG) and cellobiose, respectively. The mutant with enhanced stability would facilitate the recycle of β-glucosidase in the bioconversion of cellulosic biomass.

From amino acid to glucosinolate biosynthesis: protein sequence changes in the evolution of methylthioalkylmalate synthase in Arabidopsis

Methylthioalkylmalate synthase (MAM) catalyzes the committed step in the side chain elongation of Met, yielding important precursors for glucosinolate biosynthesis in Arabidopsis thaliana and other Brassicaceae species. MAM is believed to have evolved from isopropylmalate synthase (IPMS), an enzyme involved in Leu biosynthesis, based on phylogenetic analyses and an overlap of catalytic abilities. Here, we investigated the changes in protein structure that have occurred during the recruitment of IPMS from amino acid to glucosinolate metabolism. The major sequence difference between IPMS and MAM is the absence of 120 amino acids at the C-terminal end of MAM that constitute a regulatory domain for Leu-mediated feedback inhibition. Truncation of this domain in Arabidopsis IPMS2 results in loss of Leu feedback inhibition and quaternary structure, two features common to MAM enzymes, plus an 8.4-fold increase in the k(cat)/K(m) for a MAM substrate. Additional exchange of two amino acids in the active site resulted in a MAM-like enzyme that had little residual IPMS activity. Hence, combination of the loss of the regulatory domain and a few additional amino acid exchanges can explain the evolution of MAM from IPMS during its recruitment from primary to secondary metabolism.

Triosephosphate isomerase: a highly evolved biocatalyst

Triosephosphate isomerase (TIM) is a perfectly evolved enzyme which very fast interconverts dihydroxyacetone phosphate and D: -glyceraldehyde-3-phosphate. Its catalytic site is at the dimer interface, but the four catalytic residues, Asn11, Lys13, His95 and Glu167, are from the same subunit. Glu167 is the catalytic base. An important feature of the TIM active site is the concerted closure of loop-6 and loop-7 on ligand binding, shielding the catalytic site from bulk solvent. The buried active site stabilises the enediolate intermediate. The catalytic residue Glu167 is at the beginning of loop-6. On closure of loop-6, the Glu167 carboxylate moiety moves approximately 2 Å to the substrate. The dynamic properties of the Glu167 side chain in the enzyme substrate complex are a key feature of the proton shuttling mechanism. Two proton shuttling mechanisms, the classical and the criss-cross mechanism, are responsible for the interconversion of the substrates of this enolising enzyme.

A Paenibacillus sp. dextranase mutant pool with improved thermostability and activity

Random mutagenesis was used to create a library of chimeric dextranase (dex1) genes. A plate-screening protocol was developed with improved thermostability as a selection criterion. The mutant library was screened for active dextranase variants by observing clearing zones on dextran-blue agar plates at 50 degrees C after exposure to 68 degrees C for 2 h, a temperature regime at which wild-type activity was abolished. A number of potentially improved variants were identified by this strategy, five of which were further characterised. DNA sequencing revealed ten nucleotide substitutions, ranging from one to four per variant. Thermal inactivation studies showed reduced (2.9-fold) thermostability for one variant and similar thermostability for a second variant, but confirmed improved thermostability for three mutants with 2.3- (28.9 min) to 6.9-fold (86.6 min) increases in half-lives at 62 degrees C compared to that of the wild-type enzyme (12.6 min). Using a 10-min assay, apparent temperature optima of the variants were similar to that of the wild type (T (opt) 60 degrees C). However, one of these variants had increased enzyme activity. Therefore, the first-generation dextranase mutant pool obtained in this study has sufficient molecular diversity for further improvements in both thermostability and activity through recombination (gene shuffling).