Redesigning the substrate specificity of an enzyme by cumulative effects of the mutations of non-active site residues

Investigating role of positively selected genes and mutation sites of ERG11 in drug resistance of Candida albicans

The steep increase in acquired drug resistance in Candida isolates has posed a great challenge in the clinical management of candidiasis globally. Information of genes and codon sites that are positively selected during evolution can provide insights into the mechanisms driving antifungal resistance in Candida. This study aimed to create a manually curated list of genes of Candida spp. reported to be associated with antifungal resistance in literature, and further investigate the structure-function implications of positively selected genes and mutation sites. Sequence analysis of antifungal drug resistance associated gene sequences from various species and strains of Candida revealed that ERG11 and MRR1 of C. albicans were positively selected during evolution. Four sites in ERG11 and two sites in MRR1 of C. albicans were positively selected and associated with drug resistance. These four sites (132, 405, 450, and 464) of ERG11 are predictive markers for azole resistance and have evolved over time. A well-characterized crystal structure of sterol-14-α-demethylase (CYP51) encoded by ERG11 is available in PDB. Therefore, the stability of CYP51 in complex with fluconazole was evaluated using MD simulations and molecular docking studies for two mutations (Y132F and Y132H) reported to be associated with azole resistance in literature. These mutations induced high flexibility in functional motifs of CYP51. It was also observed that residues such as I304, G308, and I379 of CYP51 play a critical role in fluconazole binding affinity. The insights gained from this study can further guide drug design strategies addressing antimicrobial resistance.

Substrate specificity of a branch of aromatic dioxygenases determined by three distinct motifs

The inversion of substrate size specificity is an evolutionary roadblock for proteins. The Duf4243 dioxygenases GedK and BTG13 are known to catalyze the aromatic cleavage of bulky tricyclic hydroquinone. In this study, we discover a Duf4243 dioxygenase PaD that favors small monocyclic hydroquinones from the penicillic-acid biosynthetic pathway. Sequence alignments between PaD and GedK and BTG13 suggest PaD has three additional motifs, namely motifs 1-3, distributed at different positions in the protein sequence. X-ray crystal structures of PaD with the substrate at high resolution show motifs 1-3 determine three loops (loops 1-3). Most intriguing, loops 1-3 stack together at the top of the pocket, creating a lid-like tertiary structure with a narrow channel and a clearly constricted opening. This drastically changes the substrate specificity by determining the entry and binding of much smaller substrates. Further genome mining suggests Duf4243 dioxygenases with motifs 1-3 belong to an evolutionary branch that is extensively involved in the biosynthesis of natural products and has the ability to degrade diverse monocyclic hydroquinone pollutants. This study showcases how natural enzymes alter the substrate specificity fundamentally by incorporating new small motifs, with a fixed overall scaffold-architecture. It will also offer a theoretical foundation for the engineering of substrate specificity in enzymes and act as a guide for the identification of aromatic dioxygenases with distinct substrate specificities.

Structural insights into the substrate recognition of serine palmitoyltransferase from Sphingobacterium multivorum

Serine palmitoyltransferase (SPT) is a key enzyme of sphingolipid biosynthesis, which catalyzes the pyridoxal-5'-phosphate-dependent decarboxylative condensation reaction of L-serine (L-Ser) and palmitoyl-CoA (PalCoA) to form 3-ketodihydrosphingosine called long chain base (LCB). SPT is also able to metabolize L-alanine (L-Ala) and glycine (Gly), albeit with much lower efficiency. Human SPT is a membrane-bound large protein complex containing SPTLC1/SPTLC2 heterodimer as the core subunits, and it is known that mutations of the SPTLC1/SPTLC2 genes increase the formation of deoxy-type of LCBs derived from L-Ala and Gly to cause some neurodegenerative diseases. In order to study the substrate recognition of SPT, we examined the reactivity of Sphingobacterium multivorum SPT on various amino acids in the presence of PalCoA. The S. multivorum SPT could convert not only L-Ala and Gly but also L-homoserine, in addition to L-Ser, into the corresponding LCBs. Furthermore, we obtained high-quality crystals of the ligand-free form and the binary complexes with a series of amino acids, including a nonproductive amino acid, L-threonine, and determined the structures at 1.40-1.55 Å resolutions. The S. multivorum SPT accommodated various amino acid substrates through subtle rearrangements of the active-site amino acid residues and water molecules. It was also suggested that non-active-site residues mutated in the human SPT genes might indirectly influence the substrate specificity by affecting the hydrogen-bonding networks involving the bound substrate, water molecules, and amino acid residues in the active site of this enzyme. Collectively, our results highlight SPT structural features affecting substrate specificity for this stage of sphingolipid biosynthesis.

Logistic Regression-Guided Identification of Cofactor Specificity-Contributing Residues in Enzyme with Sequence Datasets Partitioned by Catalytic Properties

Changing the substrate/cofactor specificity of an enzyme requires multiple mutations at spatially adjacent positions around the substrate pocket. However, this is challenging when solely based on crystal structure information because enzymes undergo dynamic conformational changes during the reaction process. Herein, we proposed a method for estimating the contribution of each amino acid residue to substrate specificity by deploying a phylogenetic analysis with logistic regression. Since this method can estimate the candidate amino acids for mutation by ranking, it is readable and can be used in protein engineering. We demonstrated our concept using redox cofactor conversion of the Escherichia coli malic enzyme as a model, which still lacks crystal structure elucidation. The use of logistic regression with amino acid sequences classified by cofactor specificity showed that the NADP⁺-dependent malic enzyme completely switched cofactor specificity to NAD⁺ dependence without the need for a practical screening step. The model showed that surrounding residues made a greater contribution to cofactor specificity than those in the interior of the substrate pocket. These residues might be difficult to identify from crystal structure observations. We show that a highly accurate and inferential machine learning model was obtained using amino acid sequences of structurally homologous and functionally distinct enzymes as input data.

Implications of divergence of methionine adenosyltransferase in archaea

Methionine adenosyltransferase (MAT) catalyzes the biosynthesis of S-adenosyl methionine from l-methionine and ATP. MAT enzymes are ancient, believed to share a common ancestor, and are highly conserved in all three domains of life. However, the sequences of archaeal MATs show considerable divergence compared with their bacterial and eukaryotic counterparts. Furthermore, the structural significance and functional significance of this sequence divergence are not well understood. In the present study, we employed structural analysis and ancestral sequence reconstruction to investigate archaeal MAT divergence. We observed that the dimer interface containing the active site (which is usually well conserved) diverged considerably between the bacterial/eukaryotic MATs and archaeal MAT. A detailed investigation of the available structures supports the sequence analysis outcome: The protein domains and subdomains of bacterial and eukaryotic MAT are more similar than those of archaea. Finally, we resurrected archaeal MAT ancestors. Interestingly, archaeal MAT ancestors show substrate specificity, which is lost during evolution. This observation supports the hypothesis of a common MAT ancestor for the three domains of life. In conclusion, we have demonstrated that archaeal MAT is an ideal system for studying an enzyme family that evolved differently in one domain compared with others while maintaining the same catalytic activity.

Novel Biorecognition Elements against Pathogens in the Design of State-of-the-Art Diagnostics

Infectious agents, especially bacteria and viruses, account for a vast number of hospitalisations and mortality worldwide. Providing effective and timely diagnostics for the multiplicity of infectious diseases is challenging. Conventional diagnostic solutions, although technologically advanced, are highly complex and often inaccessible in resource-limited settings. An alternative strategy involves convenient rapid diagnostics which can be easily administered at the point-of-care (POC) and at low cost without sacrificing reliability. Biosensors and other rapid POC diagnostic tools which require biorecognition elements to precisely identify the causative pathogen are being developed. The effectiveness of these devices is highly dependent on their biorecognition capabilities. Naturally occurring biorecognition elements include antibodies, bacteriophages and enzymes. Recently, modified molecules such as DNAzymes, peptide nucleic acids and molecules which suffer a selective screening like aptamers and peptides are gaining interest for their biorecognition capabilities and other advantages over purely natural ones, such as robustness and lower production costs. Antimicrobials with a broad-spectrum activity against pathogens, such as antibiotics, are also used in dual diagnostic and therapeutic strategies. Other successful pathogen identification strategies use chemical ligands, molecularly imprinted polymers and Clustered Regularly Interspaced Short Palindromic Repeats-associated nuclease. Herein, the latest developments regarding biorecognition elements and strategies to use them in the design of new biosensors for pathogens detection are reviewed.

Redesign of (R)-Omega-Transaminase and Its Application for Synthesizing Amino Acids with Bulky Side Chain

ω-Transaminase (ω-TA) is an attractive biocatalyst for stereospecific preparation of amino acids and derivatives, but low catalytic efficiency and unfavorable substrate specificity hamper their industrial application. In this work, to obtain applicable (R)-ω-TA responsible for amination of α-keto acids substrates, the reactivities of eight previously synthesized ω-TAs toward pyruvate using (R)-α-methylbenzylamine ((R)-α-MBA) as amine donor were investigated, and Gibberella zeae TA (GzTA) with the highest (R)-TA activity and stereoselectivity was selected as starting scaffold for engineering. Site-directed mutagenesis around enzymatic active pocket and access tunnel identified three positive mutation sites, S214A, F113L, and V60A. Kinetic analysis synchronously with molecular docking revealed that these mutations afforded desirable alleviation of steric hindrance for pyruvate and α-MBA. Furthermore, the constructed single-, double-, and triple-mutant exhibited varying degrees of improved specificities toward bulkier α-keto acids. Using 2-oxo-2-phenylacetic acid (1d) as substrate, the conversion rate of triple-mutant F113L/V60A/S214A increased by 3.8-fold relative to that of wide-type GzTA. This study provided a practical engineering strategy for improving catalytic efficiency and substrate specificity of (R)-ω-TA. The obtained experience shed light on creating more industrial ω-TAs mutants that can accommodate structurally diverse substrates.

Biomolecular QM/MM Simulations: What Are Some of the "Burning Issues"?

QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.

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.

Structural Determinants for Substrate Selectivity in Guanine Deaminase Enzymes of the Amidohydrolase Superfamily

Guanine deaminase is a metabolic enzyme, found in all forms of life, which catalyzes the conversion of guanine to xanthine. Despite the availability of several crystal structures, the molecular determinants of substrate orientation and mechanism remain to be elucidated for the amidohydrolase family of guanine deaminase enzymes. Here, we report the crystal structures of Escherichia coli and Saccharomyces cerevisiae guanine deaminase enzymes (EcGuaD and Gud1, respectively), both members of the amidohydrolase superfamily. EcGuaD and Gud1 retain the overall TIM barrel tertiary structure conserved among amidohydrolase enzymes. Both proteins also possess a single zinc cation with trigonal bipyrimidal coordination geometry within their active sites. We also determined a liganded structure of Gud1 bound to the product, xanthine. Analysis of this structure, along with kinetic data of native and site-directed mutants of EcGuaD, identifies several key residues that are responsible for substrate recognition and catalysis. In addition, after a small library of compounds had been screened, two guanine derivatives, 8-azaguanine and 1-methylguanine, were identified as EcGuaD substrates. Interestingly, both EcGuaD and Gud1 also exhibit secondary ammeline deaminase activity. Overall, this work details key structural features of substrate recognition and catalysis of the amidohydrolase family of guanine deaminase enzymes in support of our long-term goal to engineer these enzymes with altered activity and substrate specificity.

Guide to Selecting a Biorecognition Element for Biosensors

Biosensors are powerful diagnostic tools defined as having a biorecognition element for analyte specificity and a transducer for a quantifiable signal. There are a variety of different biorecognition elements, each with unique characteristics. Understanding the advantages and disadvantages of each biorecognition element and their influence on overall biosensor performance is crucial in the planning stages to promote the success of novel biosensor development. Therefore, this review will focus on selecting the optimal biorecognition element in the preliminary design phase for novel biosensors. Included is a review of the typical characteristics and binding mechanisms of various biorecognition elements, and how they relate to biosensor performance characteristics, specifically sensitivity, selectivity, reproducibility, and reusability. The goal is to point toward language needed to improve the design and development of biosensors toward clinical success.

Structural complexity and functional diversity of plant NADPH oxidases

Plant NADPH oxidases also known as respiratory burst oxidase homologs (Rbohs) are a family of membrane-bound enzymes that play diverse roles in the defense response and morphogenetic processes via regulated generation of reactive oxygen species. Rbohs are associated with a variety of functions, although the reason for this is not clear. To evaluate using bioinformatics, the possible mechanisms for the observed functional diversity within the plant kingdom, 127 Rboh protein sequences representing 26 plant species were analyzed. Multiple clusters were identified with gene duplications that were both dicot as well as monocot-specific. The N-terminal sequences were observed to be highly variable. The conserved cysteine (equivalent of Cys890) in C-terminal of AtRbohD suggested that the redox-based modification like S-nitrosylation may regulate the activity of other Rbohs. Three-dimensional models corresponding to the N-terminal domain for Rbohs from Arabidopsis thaliana and Oryza sativa were constructed and molecular dynamics studies were carried out to study the role of Ca2+ in the folding of Rboh proteins. Certain mutations indicated possibly affect the structure and function of the plant NADPH oxidases, thereby providing the rationale for further experimental validation.

Stereoselective Nanozyme Based on Ceria Nanoparticles Engineered with Amino Acids

Stereoselectivity towards substrates is one of the most important characteristics of enzymes. Amino acids, as cofactors of many enzymes, play important roles in stereochemistry. Herein, chiral nanozymes were constructed by grafting a series of d- or l-amino acids onto the surfaces of ceria (cerium oxide) nanoparticles. We selected the most commonly used drug for combating Parkinson's disease, that is, 3,4-dihydroxyphenylalanine (DOPA) enantiomers, as examples for chiral catalysis. Through detailed kinetic studies of cerium oxide nanoparticles (CeNPs) modified with different eight amino acids, we found that phenylalanine-modified CeNP was optimal for the DOPA oxidation reaction and showed excellent stereoselectivity towards its enantiomers. l-Phenylalanine-modified CeNPs showed higher catalytic ability for oxidation of d-DOPA, while d-phenylalanine-modified CeNPs were more effective towards l-DOPA. Taken together, the results indicated that stereoselective nanozyme can be constructed by grafting nanoparticles with chiral molecules. This work may inspire better design of chiral nanozymes.

Enzyme engineering: reaching the maximal catalytic efficiency peak

The practical need for highly efficient enzymes presents new challenges in enzyme engineering, in particular, the need to improve catalytic turnover (k_cat) or efficiency (k_cat/K_M) by several orders of magnitude. However, optimizing catalysis demands navigation through complex and rugged fitness landscapes, with optimization trajectories often leading to strong diminishing returns and dead-ends. When no further improvements are observed in library screens or selections, it remains unclear whether the maximal catalytic efficiency of the enzyme (the catalytic 'fitness peak') has been reached; or perhaps, an alternative combination of mutations exists that could yield additional improvements. Here, we discuss fundamental aspects of the process of catalytic optimization, and offer practical solutions with respect to overcoming optimization plateaus.

Single-mutation fitness landscapes for an enzyme on multiple substrates reveal specificity is globally encoded

Our lack of total understanding of the intricacies of how enzymes behave has constrained our ability to robustly engineer substrate specificity. Furthermore, the mechanisms of natural evolution leading to improved or novel substrate specificities are not wholly defined. Here we generate near-comprehensive single-mutation fitness landscapes comprising >96.3% of all possible single nonsynonymous mutations for hydrolysis activity of an amidase expressed in E. coli with three different substrates. For all three selections, we find that the distribution of beneficial mutations can be described as exponential, supporting a current hypothesis for adaptive molecular evolution. Beneficial mutations in one selection have essentially no correlation with fitness for other selections and are dispersed throughout the protein sequence and structure. Our results further demonstrate the dependence of local fitness landscapes on substrate identity and provide an example of globally distributed sequence-specificity determinants for an enzyme.

Systematically understanding the sequence determinants to substrate specificity for enzymes has implications in areas from evolutionary biology to biocatalysis. Here, Whitehead and colleagues generate and analyse near-comprehensive single-mutation fitness landscapes for an amidase with three different substrates.

Directed evolution to improve the catalytic efficiency of urate oxidase from Bacillus subtilis

Urate oxidase is a key enzyme in purine metabolism and catalyzes the oxidation of uric acid to allantoin. It is used to treat hyperuricemia and gout, and also in a diagnostic kit. In this study, error-prone polymerase chain reaction and staggered extension process was used to generate a mutant urate oxidase with improved enzyme activity from Bacillus subtilis. After several rounds of mutagenesis and screening, two mutants 6E9 and 8E279 were obtained which exhibited 2.99 and 3.43 times higher catalytic efficiency, respectively. They also exhibited lower optimal reaction temperature and higher thermo-stability. D44V, Q268R and K285Q were identified as the three most beneficial amino acid substitutions introduced by site-directed mutagenesis. D44V/Q268R, which was obtained through random combination of the three mutants, displayed the highest catalytic activity. The K_m,k_cat/K_m and enzyme activity of D44V/Q268R increased by 68%, 83% and 129% respectively, compared with that of wild-type urate oxidase. Structural modeling indicated that mutations far from the active site can have significant effects on activity. For many of them, the underlying mechanisms are still difficult to explain from the static structural model. We also compared the effects of the same set of single point mutations on the wild type and on the final mutant. The results indicate strong effects of epistasis, which may imply that the mutations affect catalysis through influences on protein dynamics besides equilibrium structures.

Quantum Chemical Study of Dual-Substrate Recognition in ω-Transaminase

ω-Transaminases are attractive biocatalysts for the production of chiral amines. These enzymes usually have a broad substrate range. Their substrates include hydrophobic amines as well as amino acids, a feature referred to as dual-substrate recognition. In the present study, the reaction mechanism for the half-transamination of l-alanine to pyruvate in (S)-selective Chromobacterium violaceum ω-transaminase is investigated using density functional theory calculations. The role of a flexible arginine residue, Arg416, in the dual-substrate recognition is investigated by employing two active-site models, one including this residue and one lacking it. The results of this study are compared to those of the mechanism of the conversion of (S)-1-phenylethylamine to acetophenone. The calculations suggest that the deaminations of amino acids and hydrophobic amines follow essentially the same mechanism, but the energetics of the reactions differ significantly. It is shown that the amine is kinetically favored in the half-transamination of l-alanine/pyruvate, whereas the ketone is kinetically favored in the half-transamination of (S)-1-phenylethylamine/acetophenone. The calculations further support the proposal that the arginine residue facilitates the dual-substrate recognition by functioning as an arginine switch, where the side chain is positioned inside or outside of the active site depending on the substrate. Arg416 participates in the binding of l-alanine by forming a salt bridge to the carboxylate moiety, whereas the conversion of (S)-1-phenylethylamine is feasible in the absence of Arg416, which here represents the case in which the side chain of Arg416 is positioned outside of the active site.

In vitro Engineering of Novel Bioactivity in the Natural Enzymes

Enzymes catalyze various biochemical functions with high efficiency and specificity. In vitro design of the enzyme leads to novel bioactivity in this natural biomolecule that give answers of some vital questions like crucial residues in binding with substrate, molecular evolution, cofactor specificity etc. Enzyme engineering technology involves directed evolution, rational designing, semi-rational designing, and structure-based designing using chemical modifications. Similarly, combined computational and in vitro evolution approaches together help in artificial designing of novel bioactivity in the natural enzyme. DNA shuffling, error prone PCR and staggered extension process are used to artificially redesign active site of enzyme, which can alter its efficiency and specificity. Modifications of the enzyme can lead to the discovery of new path of molecular evolution, designing of efficient enzymes, locating active sites and crucial residues, shift in substrate, and cofactor specificity. The methods and thermodynamics of in vitro designing of the enzyme are also discussed. Similarly, engineered thermophilic and psychrophilic enzymes attain substrate specificity and activity of mesophilic enzymes that may also be beneficial for industry and therapeutics.

Conformational Tinkering Drives Evolution of a Promiscuous Activity through Indirect Mutational Effects

How remote mutations can lead to changes in enzyme function at a molecular level is a central question in evolutionary biochemistry and biophysics. Here, we combine laboratory evolution with biochemical, structural, genetic, and computational analysis to dissect the molecular basis for the functional optimization of phosphotriesterase activity in a bacterial lactonase (AiiA) from the metallo-β-lactamase (MBL) superfamily. We show that a 1000-fold increase in phosphotriesterase activity is caused by a more favorable catalytic binding position of the paraoxon substrate in the evolved enzyme that resulted from conformational tinkering of the active site through peripheral mutations. A nonmutated active site residue, Phe68, was displaced by ∼3 Å through the indirect effects of two second-shell trajectory mutations, allowing molecular interactions between the residue and paraoxon. Comparative mutational scanning, i.e., examining the effects of alanine mutagenesis on different genetic backgrounds, revealed significant changes in the functional roles of Phe68 and other nonmutated active site residues caused by the indirect effects of trajectory mutations. Our work provides a quantitative measurement of the impact of second-shell mutations on the catalytic contributions of nonmutated residues and unveils the underlying intramolecular network of strong epistatic mutational relationships between active site residues and more remote residues. Defining these long-range conformational and functional epistatic relationships has allowed us to better understand the subtle, but cumulatively significant, role of second- and third-shell mutations in evolution.

Integrating error-prone PCR and DNA shuffling as an effective molecular evolution strategy for the production of α-ketoglutaric acid byl-amino acid deaminase

l-Amino acid deaminases (LAADs; EC 1.4.3.2) belong to a family of amino acid dehydrogenases that catalyze the formation of α-keto acids froml-amino acids.

Prediction of distal residue participation in enzyme catalysis

A scoring method for the prediction of catalytically important residues in enzyme structures is presented and used to examine the participation of distal residues in enzyme catalysis. Scores are based on the Partial Order Optimum Likelihood (POOL) machine learning method, using computed electrostatic properties, surface geometric features, and information obtained from the phylogenetic tree as input features. Predictions of distal residue participation in catalysis are compared with experimental kinetics data from the literature on variants of the featured enzymes; some additional kinetics measurements are reported for variants of Pseudomonas putida nitrile hydratase (ppNH) and for Escherichia coli alkaline phosphatase (AP). The multilayer active sites of P. putida nitrile hydratase and of human phosphoglucose isomerase are predicted by the POOL log ZP scores, as is the single-layer active site of P. putida ketosteroid isomerase. The log ZP score cutoff utilized here results in over-prediction of distal residue involvement in E. coli alkaline phosphatase. While fewer experimental data points are available for P. putida mandelate racemase and for human carbonic anhydrase II, the POOL log ZP scores properly predict the previously reported participation of distal residues.

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.

Directed evolution of the substrate specificity of dialkylglycine decarboxylase

Dialkylglycine decarboxylase (DGD) is an unusual pyridoxal phosphate dependent enzyme that catalyzes decarboxylation in the first and transamination in the second half-reaction of its ping-pong catalytic cycle. Directed evolution was employed to alter the substrate specificity of DGD from 2-aminoisobutyrate (AIB) to 1-aminocyclohexane-1-carboxylate (AC6C). Four rounds of directed evolution led to the identification of several mutants, with clones in the final rounds containing five persistent mutations. The best clones show ~2.5-fold decrease in KM and ~2-fold increase in kcat, giving a modest ~5-fold increase in catalytic efficiency for AC6C. Additional rounds of directed evolution did not improve catalytic activity toward AC6C. Only one (S306F) of the five persistent mutations is close to the active site. S306F was observed in all 33 clones except one, and the mutation is shown to stabilize the enzyme toward denaturation. The other four persistent mutations are near the surface of the enzyme. The S306F mutation and the distal mutations all have significant effects on the kinetic parameters for AIB and AC6C. Molecular dynamics simulations suggest that the mutations alter the conformational landscape of the enzyme, favoring a more open active site conformation that facilitates the reactivity of the larger substrate. We speculate that the small increases in kcat/KM for AC6C are due to two constraints. The first is the mechanistic requirement for catalyzing oxidative decarboxylation via a concerted decarboxylation/proton transfer transition state. The second is that DGD must catalyze transamination at the same active site in the second half-reaction of the ping-pong catalytic cycle.

Reduced amino acid specificity of mammalian tyrosyl-tRNA synthetase is associated with elevated mistranslation of Tyr codons

Quality control operates at different steps in translation to limit errors to approximately one mistranslated codon per 10,000 codons during mRNA-directed protein synthesis. Recent studies have suggested that error rates may actually vary considerably during translation under different growth conditions. Here we examined the misincorporation of Phe at Tyr codons during synthesis of a recombinant antibody produced in tyrosine-limited Chinese hamster ovary (CHO) cells. Tyr to Phe replacements were previously found to occur throughout the antibody at a rate of up to 0.7% irrespective of the identity or context of the Tyr codon translated. Despite this comparatively high mistranslation rate, no significant change in cellular viability was observed. Monitoring of Phe and Tyr levels revealed that changes in error rates correlated with changes in amino acid pools, suggesting that mischarging of tRNA(Tyr) with noncognate Phe by tyrosyl-tRNA synthetase was responsible for mistranslation. Steady-state kinetic analyses of CHO cytoplasmic tyrosyl-tRNA synthetase revealed a 25-fold lower specificity for Tyr over Phe as compared with previously characterized bacterial enzymes, consistent with the observed increase in translation error rates during tyrosine limitation. Functional comparisons of mammalian and bacterial tyrosyl-tRNA synthetase revealed key differences at residues responsible for amino acid recognition, highlighting differences in evolutionary constraints for translation quality control.

Experimental and computational mutagenesis to investigate the positioning of a general base within an enzyme active site

The positioning of catalytic groups within proteins plays an important role in enzyme catalysis, and here we investigate the positioning of the general base in the enzyme ketosteroid isomerase (KSI). The oxygen atoms of Asp38, the general base in KSI, were previously shown to be involved in anion–aromatic interactions with two neighboring Phe residues. Here we ask whether those interactions are sufficient, within the overall protein architecture, to position Asp38 for catalysis or whether the side chains that pack against Asp38 and/or the residues of the structured loop that is capped by Asp38 are necessary to achieve optimal positioning for catalysis. To test positioning, we mutated each of the aforementioned residues, alone and in combinations, in a background with the native Asp general base and in a D38E mutant background, as Glu at position 38 was previously shown to be mispositioned for general base catalysis. These double-mutant cycles reveal positioning effects as large as 10³-fold, indicating that structural features in addition to the overall protein architecture and the Phe residues neighboring the carboxylate oxygen atoms play roles in positioning. X-ray crystallography and molecular dynamics simulations suggest that the functional effects arise from both restricting dynamic fluctuations and disfavoring potential mispositioned states. Whereas it may have been anticipated that multiple interactions would be necessary for optimal general base positioning, the energetic contributions from positioning and the nonadditive nature of these interactions are not revealed by structural inspection and require functional dissection. Recognizing the extent, type, and energetic interconnectivity of interactions that contribute to positioning catalytic groups has implications for enzyme evolution and may help reveal the nature and extent of interactions required to design enzymes that rival those found in biology.

Biocatalytic and structural properties of a highly engineered halohydrin dehalogenase

Two highly engineered halohydrin dehalogenase variants were characterized in terms of their performance in dehalogenation and epoxide cyanolysis reactions. Both enzyme variants outperformed the wild-type enzyme in the cyanolysis of ethyl (S)-3,4-epoxybutyrate, a conversion yielding ethyl (R)-4-cyano-3-hydroxybutyrate, an important chiral building block for statin synthesis. One of the enzyme variants, HheC2360, displayed catalytic rates for this cyanolysis reaction enhanced up to tenfold. Furthermore, the enantioselectivity of this variant was the opposite of that of the wild-type enzyme, both for dehalogenation and for cyanolysis reactions. The 37-fold mutant HheC2360 showed an increase in thermal stability of 8 °C relative to the wild-type enzyme. Crystal structures of this enzyme were elucidated with chloride and ethyl (S)-3,4-epoxybutyrate or with ethyl (R)-4-cyano-3-hydroxybutyrate bound in the active site. The observed increase in temperature stability was explained in terms of a substantial increase in buried surface area relative to the wild-type HheC, together with enhanced interfacial interactions between the subunits that form the tetramer. The structures also revealed that the substrate binding pocket was modified both by substitutions and by backbone movements in loops surrounding the active site. The observed changes in the mutant structures are partly governed by coupled mutations, some of which are necessary to remove steric clashes or to allow backbone movements to occur. The importance of interactions between substitutions suggests that efficient directed evolution strategies should allow for compensating and synergistic mutations during library design.

Evidence for niche partitioning revealed by the distribution of sulfur oxidation genes collected from areas of a terrestrial sulfidic spring with differing geochemical conditions

The diversity and phylogenetic significance of bacterial genes in the environment has been well studied, but comparatively little attention has been devoted to understanding the functional significance of different variations of the same metabolic gene that occur in the same environment. We analyzed the geographic distribution of 16S rRNA pyrosequences and soxB genes along a geochemical gradient in a terrestrial sulfidic spring to identify how different taxonomic variations of the soxB gene were naturally distributed within the spring outflow channel and to identify possible evidence for altered SoxB enzyme function in nature. Distinct compositional differences between bacteria that utilize their SoxB enzyme in the Paracoccus sulfide oxidation pathway (e.g., Bradyrhizobium, Paracoccus, and Rhodovulum) and bacteria that utilize their SoxB enzyme in the branched pathway (e.g., Chlorobium, Thiothrix, Thiobacillus, Halothiobacillus, and Thiomonas) were identified. Different variations of the soxB genes were present at different locations within the spring outflow channel in a manner that significantly corresponded to geochemical conditions. The distribution of the different soxB gene sequence variations suggests that the enzymes encoded by these genes are functionally different and could be optimized to specific geochemical conditions that define niche space for bacteria capable of oxidizing reduced sulfur compounds.

Structure and mechanism of a cysteine sulfinate desulfinase engineered on the aspartate aminotransferase scaffold

The joint substitution of three active-site residues in Escherichia coli (L)-aspartate aminotransferase increases the ratio of l-cysteine sulfinate desulfinase to transaminase activity 10(5)-fold. This change in reaction specificity results from combining a tyrosine-shift double mutation (Y214Q/R280Y) with a non-conservative substitution of a substrate-binding residue (I33Q). Tyr214 hydrogen bonds with O3 of the cofactor and is close to Arg374 which binds the α-carboxylate group of the substrate; Arg280 interacts with the distal carboxylate group of the substrate; and Ile33 is part of the hydrophobic patch near the entrance to the active site, presumably participating in the domain closure essential for the transamination reaction. In the triple-mutant enzyme, k(cat)' for desulfination of l-cysteine sulfinate increased to 0.5s(-1) (from 0.05s(-1) in wild-type enzyme), whereas k(cat)' for transamination of the same substrate was reduced from 510s(-1) to 0.05s(-1). Similarly, k(cat)' for β-decarboxylation of l-aspartate increased from<0.0001s(-1) to 0.07s(-1), whereas k(cat)' for transamination was reduced from 530s(-1) to 0.13s(-1). l-Aspartate aminotransferase had thus been converted into an l-cysteine sulfinate desulfinase that catalyzes transamination and l-aspartate β-decarboxylation as side reactions. The X-ray structures of the engineered l-cysteine sulfinate desulfinase in its pyridoxal-5'-phosphate and pyridoxamine-5'-phosphate form or liganded with a covalent coenzyme-substrate adduct identified the subtle structural changes that suffice for generating desulfinase activity and concomitantly abolishing transaminase activity toward dicarboxylic amino acids. Apparently, the triple mutation impairs the domain closure thus favoring reprotonation of alternative acceptor sites in coenzyme-substrate intermediates by bulk water.

Structure-function analysis from the outside in: long-range tertiary contacts in RNA exhibit distinct catalytic roles

The conserved catalytic core of the Tetrahymena group I ribozyme is encircled by peripheral elements. We have conducted a detailed structure-function study of the five long-range tertiary contacts that fasten these distal elements together. Mutational ablation of each of the tertiary contacts destabilizes the folded ribozyme, indicating a role of the peripheral elements in overall stability. Once folded, three of the five tertiary contact mutants exhibit defects in overall catalysis that range from 20- to 100-fold. These and the subsequent results indicate that the structural ring of peripheral elements does not act as a unitary element; rather, individual connections have distinct roles as further revealed by kinetic and thermodynamic dissection of the individual reaction steps. Ablation of P14 or the metal ion core/metal ion core receptor (MC/MCR) destabilizes docking of the substrate-containing P1 helix into tertiary interactions with the ribozyme's conserved core. In contrast, ablation of the L9/P5 contact weakens binding of the guanosine nucleophile by slowing its association, without affecting P1 docking. The P13 and tetraloop/tetraloop receptor (TL/TLR) mutations had little functional effect and small, local structural changes, as revealed by hydroxyl radical footprinting, whereas the P14, MC/MCR, and L9/P5 mutants show structural changes distal from the mutation site. These changes extended into regions of the catalytic core involved in docking or guanosine binding. Thus, distinct allosteric pathways couple the long-range tertiary contacts to functional sites within the conserved core. This modular functional specialization may represent a fundamental strategy in RNA structure-function interrelationships.

Tightening of active site interactions en route to the transition state revealed by single-atom substitution in the guanosine-binding site of the Tetrahymena group I ribozyme

Protein enzymes establish intricate networks of interactions to bind and position substrates and catalytic groups within active sites, enabling stabilization of the chemical transition state. Crystal structures of several RNA enzymes also suggest extensive interaction networks, despite RNA's structural limitations, but there is little information on the functional and the energetic properties of these inferred networks. We used double mutant cycles and presteady-state kinetic analyses to probe the putative interaction between the exocyclic amino group of the guanosine nucleophile and the N7 atom of residue G264 of the Tetrahymena group I ribozyme. As expected, the results supported the presence of this interaction, but remarkably, the energetic penalty for introducing a CH group at the 7-position of residue G264 accumulates as the reaction proceeds toward the chemical transition state to a total of 6.2 kcal/mol. Functional tests of neighboring interactions revealed that the presence of the CH group compromises multiple contacts within the interaction network that encompass the reactive elements, apparently forcing the nucleophile to bind and attack from an altered, suboptimal orientation. The energetic consequences of this indirect disruption of neighboring interactions as the reaction proceeds demonstrate that linkage between binding interactions and catalysis hinges critically on the precise structural integrity of a network of interacting groups.

Enhancement of the latent 3-isopropylmalate dehydrogenase activity of promiscuous homoisocitrate dehydrogenase by directed evolution

HICDH (homoisocitrate dehydrogenase), which is involved in lysine biosynthesis through α-aminoadipate, is a paralogue of IPMDH [3-IPM (3-isopropylmalate) dehydrogenase], which is involved in leucine biosynthesis. TtHICDH (Thermus thermophilus HICDH) can recognize isocitrate, as well as homoisocitrate, as the substrate, and also shows IPMDH activity, although at a considerably decreased rate. In the present study, the promiscuous TtHICDH was evolved into an enzyme showing distinct IPMDH activity by directed evolution using a DNA-shuffling technique. Through five repeats of DNA shuffling/screening, variants that allowed Escherichia coli C600 (leuB−) to grow on a minimal medium in 2 days were obtained. One of the variants LR5–1, with eight amino acid replacements, was found to possess a 65-fold increased kcat/Km value for 3-IPM, compared with TtHICDH. Introduction of a single back-replacement H15Y change caused a further increase in the kcat/Km value and a partial recovery of the decreased thermotolerance of LR5–1. Site-directed mutagenesis revealed that most of the amino acid replacements found in LR5–1 effectively increased IPMDH activity; replacements around the substrate-binding site contributed to the improved recognition for 3-IPM, and other replacements at sites away from the substrate-binding site enhanced the turnover number for the IPMDH reaction. The crystal structure of LR5–1 was determined at 2.4 Å resolution and revealed that helix α4 was displaced in a manner suitable for recognition of the hydrophobic γ-moiety of 3-IPM. On the basis of the crystal structure, possible reasons for enhancement of the turnover number are discussed.

Structural insights into substrate specificity in variants of N-acetylneuraminic Acid lyase produced by directed evolution

The substrate specificity of Escherichia coli N-acetylneuraminic acid lyase was previously switched from the natural condensation of pyruvate with N-acetylmannosamine, yielding N-acetylneuraminic acid, to the aldol condensation generating N-alkylcarboxamide analogues of N-acetylneuraminic acid. This was achieved by a single mutation of Glu192 to Asn. In order to analyze the structural changes involved and to more fully understand the basis of this switch in specificity, we have isolated all 20 variants of the enzyme at position 192 and determined the activities with a range of substrates. We have also determined five high-resolution crystal structures: the structures of wild-type E. coli N-acetylneuraminic acid lyase in the presence and in the absence of pyruvate, the structures of the E192N variant in the presence and in the absence of pyruvate, and the structure of the E192N variant in the presence of pyruvate and a competitive inhibitor (2R,3R)-2,3,4-trihydroxy-N,N-dipropylbutanamide. All structures were solved in space group P2₁ at resolutions ranging from 1.65 Å to 2.2 Å. A comparison of these structures, in combination with the specificity profiles of the variants, reveals subtle differences that explain the details of the specificity changes. This work demonstrates the subtleties of enzyme–substrate interactions and the importance of determining the structures of enzymes produced by directed evolution, where the specificity determinants may change from one substrate to another.

Evolution of enzymatic activities of testis-specific short-chain dehydrogenase/reductase in Drosophila

The testis-specific gene Jingwei (jgw) is a newly evolved short-chain dehydrogenase/reductase in Drosophila. Preliminary substrate screening indicated that JGW prefers long-chain primary alcohols as substrates, including several exotic alcohols such as farnesol and geraniol. Using steady-state kinetics analyses and molecular docking, we not only confirmed JGW's substrate specificity, but also demonstrated that the new enzymatic activities of JGW evolved extensively after exon-shuffling from a preexisting enzyme. Analysis of JGW orthologs in sister species shows that subsequent evolutionary changes following the birth of JGW altered substrate specificities and enzyme stabilities. Our results lend support to a general mechanism for the evolution of a new enzyme, in which catalytic chemistry evolves first followed by diversification of substrate utilization.

NMR and bioinformatics discovery of exosites that tune metalloelastase specificity for solubilized elastin and collagen triple helices

The catalytic domain of metalloelastase (matrix metalloproteinase-12 or MMP-12) is unique among MMPs in exerting high proteolytic activity upon fibrils that resist hydrolysis, especially elastin from lungs afflicted with chronic obstructive pulmonary disease or arteries with aneurysms. How does the MMP-12 catalytic domain achieve this specificity? NMR interface mapping suggests that α-elastin species cover the primed subsites, a strip across the β-sheet from β-strand IV to the II-III loop, and a broad bowl from helix A to helix C. The many contacts may account for the comparatively high affinity, as well as embedding of MMP-12 in damaged elastin fibrils in vivo. We developed a strategy called BINDSIght, for bioinformatics and NMR discovery of specificity of interactions, to evaluate MMP-12 specificity without a structure of a complex. BINDSIght integration of the interface mapping with other ambiguous information from sequences guided choice mutations in binding regions nearer the active site. Single substitutions at each of ten locations impair specific activity toward solubilized elastin. Five of them impair release of peptides from intact elastin fibrils. Eight lesions also impair specific activity toward triple helices from collagen IV or V. Eight sites map to the "primed" side in the III-IV, V-B, and S1' specificity loops. Two map to the "unprimed" side in the IV-V and B-C loops. The ten key residues circumscribe the catalytic cleft, form an exosite, and are distinctive features available for targeting by new diagnostics or therapeutics.

High-throughput metabolic engineering: advances in small-molecule screening and selection

Metabolic engineering for the overproduction of high-value small molecules is dependent upon techniques in directed evolution to improve production titers. The majority of small molecules targeted for overproduction are inconspicuous and cannot be readily obtained by screening. We provide a review on the development of high-throughput colorimetric, fluorescent, and growth-coupled screening techniques, enabling inconspicuous small-molecule detection. We first outline constraints on throughput imposed during the standard directed evolution workflow (library construction, transformation, and screening) and establish a screening and selection ladder on the basis of small-molecule assay throughput and sensitivity. An in-depth analysis of demonstrated screening and selection approaches for small-molecule detection is provided. Particular focus is placed on in vivo biosensor-based detection methods that reduce or eliminate in vitro assay manipulations and increase throughput. We conclude by providing our prospectus for the future, focusing on transcription factor-based detection systems as a natural microbial mode of small-molecule detection.

Induced allostery in the directed evolution of an enantioselective Baeyer-Villiger monooxygenase

The molecular basis of allosteric effects, known to be caused by an effector docking to an enzyme at a site distal from the binding pocket, has been studied recently by applying directed evolution. Here, we utilize laboratory evolution in a different way, namely to induce allostery by introducing appropriate distal mutations that cause domain movements with concomitant reshaping of the binding pocket in the absence of an effector. To test this concept, the thermostable Baeyer-Villiger monooxygenase, phenylacetone monooxygenase (PAMO), was chosen as the enzyme to be employed in asymmetric Baeyer-Villiger reactions of substrates that are not accepted by the wild type. By using the known X-ray structure of PAMO, a decision was made regarding an appropriate site at which saturation mutagenesis is most likely to generate mutants capable of inducing allostery without any effector compound being present. After screening only 400 transformants, a double mutant was discovered that catalyzes the asymmetric oxidative kinetic resolution of a set of structurally different 2-substituted cyclohexanone derivatives as well as the desymmetrization of three different 4-substituted cyclohexanones, all with high enantioselectivity. Molecular dynamics (MD) simulations and covariance maps unveiled the origin of increased substrate scope as being due to allostery. Large domain movements occur that expose and reshape the binding pocket. This type of focused library production, aimed at inducing significant allosteric effects, is a viable alternative to traditional approaches to "designed" directed evolution that address the binding site directly.

Directed evolution and structural characterization of a simvastatin synthase

Enzymes from natural product biosynthetic pathways are attractive candidates for creating tailored biocatalysts to produce semisynthetic pharmaceutical compounds. LovD is an acyltransferase that converts the inactive monacolin J acid (MJA) into the cholesterol-lowering lovastatin. LovD can also synthesize the blockbuster drug simvastatin using MJA and a synthetic alpha-dimethylbutyryl thioester, albeit with suboptimal properties as a biocatalyst. Here we used directed evolution to improve the properties of LovD toward semisynthesis of simvastatin. Mutants with improved catalytic efficiency, solubility, and thermal stability were obtained, with the best mutant displaying an approximately 11-fold increase in an Escherichia coli-based biocatalytic platform. To understand the structural basis of LovD enzymology, seven X-ray crystal structures were determined, including the parent LovD, an improved mutant G5, and G5 cocrystallized with ligands. Comparisons between the structures reveal that beneficial mutations stabilize the structure of G5 in a more compact conformation that is favorable for catalysis.

Enzyme engineering for enantioselectivity: from trial-and-error to rational design?

The availability of tailored enzymes is crucial for the implementation of biocatalysis in organic chemistry. Enantioselectivity is one key parameter defining the usefulness of an enzyme and, therefore, the competitiveness of the corresponding industrial process. Hence, identification of enzymes with high enantioselectivity in the desired transformation is important. Currently, this is achieved by screening collections and libraries comprising natural or man-made diversity for the desired trait. Recently, a variety of improved methods have been developed to generate and screen this diversity more efficiently. Here, we present and discuss the most important advances in both library generation and screening. We also evaluate future trends, such as moving from random evolution to more rational.

In silico modelling of directed evolution: Implications for experimental design and stepwise evolution

We model the process of directed evolution (DE) in silico using genetic algorithms. Making use of the NK fitness landscape model, we analyse the effects of mutation rate, crossover and selection pressure on the performance of DE. A range of values of K, the epistatic interaction of the landscape, are considered, and high- and low-throughput modes of evolution are compared. Our findings suggest that for runs of or around ten generations' duration-as is typical in DE-there is little difference between the way in which DE needs to be configured in the high- and low-throughput regimes, nor across different degrees of landscape epistasis. In all cases, a high selection pressure (but not an extreme one) combined with a moderately high mutation rate works best, while crossover provides some benefit but only on the less rugged landscapes. These genetic algorithms were also compared with a "model-based approach" from the literature, which uses sequential fixing of the problem parameters based on fitting a linear model. Overall, we find that purely evolutionary techniques fare better than do model-based approaches across all but the smoothest landscapes.