Human GST P1-1 Redesigned for Enhanced Catalytic Activity with the Anticancer Prodrug Telcyta and Improved Thermostability.

Protein engineering can be used to tailor enzymes for medical purposes, including antibody-directed enzyme prodrug therapy (ADEPT), which can act as a tumor-targeted alternative to conventional chemotherapy for cancer. In ADEPT, the antibody serves as a vector, delivering a drug-activating enzyme selectively to the tumor site. Glutathione transferases (GSTs) are a family of naturally occurring detoxication enzymes, and the finding that some of them are overexpressed in tumors has been exploited to develop GST-activated prodrugs. The prodrug Telcyta is activated by GST P1-1, which is the GST most commonly elevated in cancer cells, implying that tumors overexpressing GST P1-1 should be particularly vulnerable to Telcyta. Promising antitumor activity has been noted in clinical trials, but the wildtype enzyme has modest activity with Telcyta, and further functional improvement would enhance its usefulness for ADEPT. We utilized protein engineering to construct human GST P1-1 gene variants in the search for enzymes with enhanced activity with Telcyta. The variant Y109H displayed a 2.9-fold higher enzyme activity compared to the wild-type GST P1-1. However, increased catalytic potency was accompanied by decreased thermal stability of the Y109H enzyme, losing 99% of its activity in 8 min at 50 °C. Thermal stability was restored by four additional mutations simultaneously introduced without loss of the enhanced activity with Telcyta. The mutation Q85R was identified as an important contributor to the regained thermostability. These results represent a first step towards a functional ADEPT application for Telcyta.

"Site and Mutation"-Specific Predictions Enable Minimal Directed Evolution Libraries.

Directed evolution experiments designed to improve the activity of a biocatalyst have increased in sophistication from the early days of completely random mutagenesis. Sequence-based and structure-based methods have been developed to identify "hotspot" positions that when randomized provide a higher frequency of beneficial mutations that improve activity. These focused mutagenesis methods reduce library sizes and therefore reduce screening burden, accelerating the rate of finding improved enzymes. Looking for further acceleration in finding improved enzymes, we investigated whether two existing methods, one sequence-based (Protein GPS) and one structure-based (using Bioluminate and MOE), were sufficiently predictive to provide not just the hotspot position, but also the amino acid substitution that improved activity at that position. By limiting the libraries to variants that contained only specific amino acid substitutions, library sizes were kept to less than 100 variants. For an initial round of ATA-117 R-selective transaminase evolution, we found that the methods used produced libraries where 9% and 18% of the amino acid substitutions chosen were amino acids that improved reaction performance in lysates. The ability to create combinations of mutations as part of the initial design was confounded by the relatively large number of predicted mutations that were inactivating (30% and 45% for the sequence-based and structure-based methods, respectively). Despite this, combining several mutations identified within a given method produced variant lysates 7- and 9-fold more active than the wild-type lysate, highlighting the capability of mutations chosen this way to generate large advances in activity in addition to the reductions in screening.

Exploring sequence-function space of a poplar glutathione transferase using designed information-rich gene variants.

Exploring the vicinity around a locus of a protein in sequence space may identify homologs with enhanced properties, which could become valuable in biotechnical and other applications. A rational approach to this pursuit is the use of 'infologs', i.e. synthetic sequences with specific substitutions capturing maximal sequence information derived from the evolutionary history of the protein family. Ninety-five such infolog genes of poplar glutathione transferase were synthesized and expressed in Escherichia coli, and the catalytic activities of the proteins determined with alternative substrates. Sequence-activity relationships derived from the infologs were used to design a second set of 47 infologs in which 90% of the members exceeded wild-type properties. Two mutants, C2 (V55I/E95D/D108E/A160V) and G5 (F13L/C70A/G122E), were further functionally characterized. The activities of the infologs with the alternative substrates 1-chloro-2,4-dinitrobenzene and phenethyl isothiocyanate, subject to different chemistries, were positively correlated, indicating that the examined mutations were affecting the overall catalytic competence without major shift in substrate discrimination. By contrast, the enhanced protein expressivity observed in many of the mutants were not similarly correlated with the activities. In conclusion, small libraries of well-defined infologs can be used to systematically explore sequence space to optimize proteins in multidimensional functional space.

Mapping of amino acid substitutions conferring herbicide resistance in wheat glutathione transferase.

We have used design of experiments (DOE) and systematic variance to efficiently explore glutathione transferase substrate specificities caused by amino acid substitutions. Amino acid substitutions selected using phylogenetic analysis were synthetically combined using a DOE design to create an information-rich set of gene variants, termed infologs. We used machine learning to identify and quantify protein sequence-function relationships against 14 different substrates. The resulting models were quantitative and predictive, serving as a guide for engineering of glutathione transferase activity toward a diverse set of herbicides. Predictive quantitative models like those presented here have broad applicability for bioengineering.

Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Improved biocatalysts from a synthetic circular permutation library of the flavin-dependent oxidoreductase old yellow enzyme.

Members of the old yellow enzyme (OYE) family are widely used, effective biocatalysts for the stereoselective trans-hydrogenation of activated alkenes. To further expand their substrate scope and improve catalytic performance, we have applied a protein engineering strategy called circular permutation (CP) to enhance the function of OYE1 from Saccharomyces pastorianus. CP can influence a biocatalyst's function by altering protein backbone flexibility and active site accessibility, both critical performance features because the catalytic cycle for OYE1 is thought to involve rate-limiting conformational changes. To explore the impact of CP throughout the OYE1 protein sequence, we implemented a highly efficient approach for cell-free cpOYE library preparation by combining whole-gene synthesis with in vitro transcription/translation. The versatility of such an ex vivo system was further demonstrated by the rapid and reliable functional evaluation of library members under variable environmental conditions with three reference substrates ketoisophorone, cinnamaldehyde, and (S)-carvone. Library analysis identified over 70 functional OYE1 variants with several biocatalysts exhibiting over an order of magnitude improved catalytic activity. Although catalytic gains of individual cpOYE library members vary by substrate, the locations of new protein termini in functional variants for all tested substates fall within the same four distinct loop/lid regions near the active site. Our findings demonstrate the importance of these structural elements in enzyme function and support the hypothesis of conformational flexibility as a limiting factor for catalysis in wild type OYE.

Redesigning and characterizing the substrate specificity and activity of Vibrio fluvialis aminotransferase for the synthesis of imagabalin.

Several protein engineering approaches were combined to optimize the selectivity and activity of Vibrio fluvialis aminotransferase (Vfat) for the synthesis of (3S,5R)-ethyl 3-amino-5-methyloctanoate; a key intermediate in the synthesis of imagabalin, an advanced candidate for the treatment of generalized anxiety disorder. Starting from wild-type Vfat, which had extremely low activity catalyzing the desired reaction, we engineered an improved enzyme with a 60-fold increase in initial reaction velocity for transamination of (R)-ethyl 5-methyl 3-oxooctanoate to (3S,5R)-ethyl 3-amino-5-methyloctanoate. To achieve this, <450 variants were screened, which allowed accurate assessment of enzyme performance using a low-throughput ultra performance liquid chromatography assay. During the course of this work, crystal structures of Vfat wild type and an improved variant (Vfat variant r414) were solved and they are reported here for the first time. This work also provides insight into the critical residues for substrate specificity for the transamination of (R)-ethyl 5-methyl 3-oxooctanoate and structurally related ß-ketoesters.

Patterning of antibodies using flexographic printing.

Antibodies were patterned onto flexible plastic films using the flexographic printing process. An ink formulation was developed using high molecular weight polyvinyl alcohol in carbonate-bicarbonate buffer. In order to aid both antibody adhesion and the quality of definition in the printed features, a nitrocellulose coating was developed that was capable of being discretely patterned, thus increasing the signal-to-noise ratio of an antibody array. Printing antibody features such as dots, squares, text, and fine lines were reproduced effectively. Furthermore, this process could be easily adapted for printing of other biological materials, including, but not limited to, enzymes, DNA, proteins, aptamers, and cells.

Engineering genes for predictable protein expression.

The DNA sequence used to encode a polypeptide can have dramatic effects on its expression. Lack of readily available tools has until recently inhibited meaningful experimental investigation of this phenomenon. Advances in synthetic biology and the application of modern engineering approaches now provide the tools for systematic analysis of the sequence variables affecting heterologous expression of recombinant proteins. We here discuss how these new tools are being applied and how they circumvent the constraints of previous approaches, highlighting some of the surprising and promising results emerging from the developing field of gene engineering.

Analysis of the ketosynthase-chain length factor heterodimer from the fredericamycin polyketide synthase.

The pentadecaketide fredericamycin has the longest carbon chain backbone among polycyclic aromatic polyketide antibiotics whose biosynthetic genes have been sequenced. This backbone is synthesized by the bimodular fdm polyketide synthase (PKS). Here, we demonstrate that the bimodular fdm PKS as well as its elongation module alone synthesize undecaketides and dodecaketides. Thus, unlike other homologs, the fdm ketosynthase-chain length factor (KS-CLF) heterodimer does not exclusively control the backbone length of its natural product. Using sequence- and structure-based approaches, 48 CLF multiple mutants were engineered and analyzed. Unexpectedly, the I134F mutant was unable to turn over but could initiate and partially elongate the polyketide chain. This unprecedented mutant suggests that the KS-CLF heterodimer harbors an as yet uncharacterized chain termination mechanism. Together, our findings reveal fundamental mechanistic differences between the fdm PKS and its well-studied homologs.

Reconstructed evolutionary adaptive paths give polymerases accepting reversible terminators for sequencing and SNP detection.

Any system, natural or human-made, is better understood if we analyze both its history and its structure. Here we combine structural analyses with a "Reconstructed Evolutionary Adaptive Path" (REAP) analysis that used the evolutionary and functional history of DNA polymerases to replace amino acids to enable polymerases to accept a new class of triphosphate substrates, those having their 3'-OH ends blocked as a 3(')-ONH(2) group (dNTP-ONH(2)). Analogous to widely used 2',3'-dideoxynucleoside triphosphates (ddNTPs), dNTP-ONH(2)s terminate primer extension. Unlike ddNTPs, however, primer extension can be resumed by cleaving an O-N bond to restore an -OH group to the 3'-end of the primer. REAP combined with crystallographic analyses identified 35 sites where replacements might improve the ability of Taq to accept dNTP-ONH(2)s. A library of 93 Taq variants, each having replacements at three or four of these sites, held eight variants having improved ability to accept dNTP-ONH(2) substrates. Two of these (A597T, L616A, F667Y, E745H, and E520G, K540I, L616A) performed notably well. The second variant incorporated both dNTP-ONH(2)sand ddNTPs faithfully and efficiently, supporting extension-cleavage-extension cycles applicable in parallel sequencing and in SNP detection through competition between reversible and irreversible terminators. Dissecting these results showed that one replacement (L616A), not previously identified, allows Taq to incorporate both reversible and irreversible terminators. Modeling showed how L616A might open space behind Phe-667, allowing it to move to accommodate the larger 3'-substituent. This work provides polymerases for DNA analyses and shows how evolutionary analyses help explore relationships between structure and function in proteins.

Design parameters to control synthetic gene expression in Escherichia coli.

BACKGROUND: Production of proteins as therapeutic agents, research reagents and molecular tools frequently depends on expression in heterologous hosts. Synthetic genes are increasingly used for protein production because sequence information is easier to obtain than the corresponding physical DNA. Protein-coding sequences are commonly re-designed to enhance expression, but there are no experimentally supported design principles. PRINCIPAL FINDINGS: To identify sequence features that affect protein expression we synthesized and expressed in E. coli two sets of 40 genes encoding two commercially valuable proteins, a DNA polymerase and a single chain antibody. Genes differing only in synonymous codon usage expressed protein at levels ranging from undetectable to 30% of cellular protein. Using partial least squares regression we tested the correlation of protein production levels with parameters that have been reported to affect expression. We found that the amount of protein produced in E. coli was strongly dependent on the codons used to encode a subset of amino acids. Favorable codons were predominantly those read by tRNAs that are most highly charged during amino acid starvation, not codons that are most abundant in highly expressed E. coli proteins. Finally we confirmed the validity of our models by designing, synthesizing and testing new genes using codon biases predicted to perform well. CONCLUSION: The systematic analysis of gene design parameters shown in this study has allowed us to identify codon usage within a gene as a critical determinant of achievable protein expression levels in E. coli. We propose a biochemical basis for this, as well as design algorithms to ensure high protein production from synthetic genes. Replication of this methodology should allow similar design algorithms to be empirically derived for any expression system.

Engineering the Salmonella type III secretion system to export spider silk monomers.

The type III secretion system (T3SS) exports proteins from the cytoplasm, through both the inner and outer membranes, to the external environment. Here, a system is constructed to harness the T3SS encoded within Salmonella Pathogeneity Island 1 to export proteins of biotechnological interest. The system is composed of an operon containing the target protein fused to an N-terminal secretion tag and its cognate chaperone. Transcription is controlled by a genetic circuit that only turns on when the cell is actively secreting protein. The system is refined using a small human protein (DH domain) and demonstrated by exporting three silk monomers (ADF-1, -2, and -3), representative of different types of spider silk. Synthetic genes encoding silk monomers were designed to enhance genetic stability and codon usage, constructed by automated DNA synthesis, and cloned into the secretion control system. Secretion rates up to 1.8 mg l(-1) h(-1) are demonstrated with up to 14% of expressed protein secreted. This work introduces new parts to control protein secretion in Gram-negative bacteria, which will be broadly applicable to problems in biotechnology.

SCHEMA recombination of a fungal cellulase uncovers a single mutation that contributes markedly to stability.

A quantitative linear model accurately (R(2) = 0.88) describes the thermostabilities of 54 characterized members of a family of fungal cellobiohydrolase class II (CBH II) cellulase chimeras made by SCHEMA recombination of three fungal enzymes, demonstrating that the contributions of SCHEMA sequence blocks to stability are predominantly additive. Thirty-one of 31 predicted thermostable CBH II chimeras have thermal inactivation temperatures higher than the most thermostable parent CBH II, from Humicola insolens, and the model predicts that hundreds more CBH II chimeras share this superior thermostability. Eight of eight thermostable chimeras assayed hydrolyze the solid cellulosic substrate Avicel at temperatures at least 5 degrees C above the most stable parent, and seven of these showed superior activity in 16-h Avicel hydrolysis assays. The sequence-stability model identified a single block of sequence that adds 8.5 degrees C to chimera thermostability. Mutating individual residues in this block identified the C313S substitution as responsible for the entire thermostabilizing effect. Introducing this mutation into the two recombination parent CBH IIs not featuring it (Hypocrea jecorina and H. insolens) decreased inactivation, increased maximum Avicel hydrolysis temperature, and improved long time hydrolysis performance. This mutation also stabilized and improved Avicel hydrolysis by Phanerochaete chrysosporium CBH II, which is only 55-56% identical to recombination parent CBH IIs. Furthermore, the C313S mutation increased total H. jecorina CBH II activity secreted by the Saccharomyces cerevisiae expression host more than 10-fold. Our results show that SCHEMA structure-guided recombination enables quantitative prediction of cellulase chimera thermostability and efficient identification of stabilizing mutations.

A family of thermostable fungal cellulases created by structure-guided recombination.

SCHEMA structure-guided recombination of 3 fungal class II cellobiohydrolases (CBH II cellulases) has yielded a collection of highly thermostable CBH II chimeras. Twenty-three of 48 genes sampled from the 6,561 possible chimeric sequences were secreted by the Saccharomyces cerevisiae heterologous host in catalytically active form. Five of these chimeras have half-lives of thermal inactivation at 63 degrees C that are greater than the most stable parent, CBH II enzyme from the thermophilic fungus Humicola insolens, which suggests that this chimera collection contains hundreds of highly stable cellulases. Twenty-five new sequences were designed based on mathematical modeling of the thermostabilities for the first set of chimeras. Ten of these sequences were expressed in active form; all 10 retained more activity than H. insolens CBH II after incubation at 63 degrees C. The total of 15 validated thermostable CBH II enzymes have high sequence diversity, differing from their closest natural homologs at up to 63 amino acid positions. Selected purified thermostable chimeras hydrolyzed phosphoric acid swollen cellulose at temperatures 7 to 15 degrees C higher than the parent enzymes. These chimeras also hydrolyzed as much or more cellulose than the parent CBH II enzymes in long-time cellulose hydrolysis assays and had pH/activity profiles as broad, or broader than, the parent enzymes. Generating this group of diverse, thermostable fungal CBH II chimeras is the first step in building an inventory of stable cellulases from which optimized enzyme mixtures for biomass conversion can be formulated.

Protein engineering of improved prolyl endopeptidases for celiac sprue therapy.

Due to their unique ability to cleave immunotoxic gluten peptides endoproteolytically, prolyl endopeptidases (PEPs) are attractive oral therapeutic candidates for protecting celiac sprue patients from the toxic effects of dietary gluten. Enhancing the activity and stability of PEPs under gastric conditions (low pH, high pepsin concentration) is a challenge for protein engineers. Using a combination of sequence- and structure-based approaches together with machine learning algorithms, we have identified improved variants of the Sphingomonas capsulata PEP, a target of clinical relevance. Through two rounds of iterative mutagenesis and analysis, variants with as much as 20% enhanced specific activity at pH 4.5 and 200-fold greater resistance to pepsin were identified. Our results vividly reinforce the concept that conservative changes in proteins, especially in hydrophobic residues within tightly packed regions, can profoundly influence protein structure and function in ways that are difficult to predict entirely from first principles and must therefore be optimized through iterative design and analytical cycles. Incubation with whole wheat bread under simulated gastric conditions also suggests that some variants have pharmacologically significant improvements in gluten detoxification activity.

Palaeotemperature trend for Precambrian life inferred from resurrected proteins.

Biosignatures and structures in the geological record indicate that microbial life has inhabited Earth for the past 3.5 billion years or so. Research in the physical sciences has been able to generate statements about the ancient environment that hosted this life. These include the chemical compositions and temperatures of the early ocean and atmosphere. Only recently have the natural sciences been able to provide experimental results describing the environments of ancient life. Our previous work with resurrected proteins indicated that ancient life lived in a hot environment. Here we expand the timescale of resurrected proteins to provide a palaeotemperature trend of the environments that hosted life from 3.5 to 0.5 billion years ago. The thermostability of more than 25 phylogenetically dispersed ancestral elongation factors suggest that the environment supporting ancient life cooled progressively by 30 degrees C during that period. Here we show that our results are robust to potential statistical bias associated with the posterior distribution of inferred character states, phylogenetic ambiguity, and uncertainties in the amino-acid equilibrium frequencies used by evolutionary models. Our results are further supported by a nearly identical cooling trend for the ancient ocean as inferred from the deposition of oxygen isotopes. The convergence of results from natural and physical sciences suggest that ancient life has continually adapted to changes in environmental temperatures throughout its evolutionary history.

Directed evolution of gene-shuffled IFN-alpha molecules with activity profiles tailored for treatment of chronic viral diseases.

Type I IFNs are unusually pleiotropic cytokines that bind to a single heterodimeric receptor and have potent antiviral, antiproliferative, and immune modulatory activities. The diverse effects of the type I IFNs are of differential therapeutic importance; in cancer therapy, an enhanced antiproliferative effect may be beneficial, whereas in the therapy of viral infections (such as hepatitis B and hepatitis C), the antiproliferative effects lead to dose limiting bone marrow suppression. Studies have shown that various members of the natural IFN-alpha family and engineered variants, such as IFN-con1, vary in the ratios between various IFN-mediated cellular activities. We used DNA shuffling to explore and confirm the hypothesis that one could simultaneously increase the antiviral and Th1-inducing activity and decrease the antiproliferative activity. We report IFN-alpha hybrids wherein the ratio of antiviral:antiproliferative and Th1-inducing: antiproliferative potencies are markedly increased with respsect to IFN-con1 (75- and 80-fold, respectively). A four-residue motif that overlaps with the IFNAR1 binding site and is derived by cross breeding with a pseudogene contributes significantly to this phenotype. These IFN-alphas have an activity profile that may result in an improved therapeutic index and, consequently, better clinical efficacy for the treatment of chronic viral diseases such as hepatitis B virus, human papilloma virus, HIV, or chronic hepatitis C.

Engineering proteinase K using machine learning and synthetic genes.

BACKGROUND: Altering a protein's function by changing its sequence allows natural proteins to be converted into useful molecular tools. Current protein engineering methods are limited by a lack of high throughput physical or computational tests that can accurately predict protein activity under conditions relevant to its final application. Here we describe a new synthetic biology approach to protein engineering that avoids these limitations by combining high throughput gene synthesis with machine learning-based design algorithms. RESULTS: We selected 24 amino acid substitutions to make in proteinase K from alignments of homologous sequences. We then designed and synthesized 59 specific proteinase K variants containing different combinations of the selected substitutions. The 59 variants were tested for their ability to hydrolyze a tetrapeptide substrate after the enzyme was first heated to 68 degrees C for 5 minutes. Sequence and activity data was analyzed using machine learning algorithms. This analysis was used to design a new set of variants predicted to have increased activity over the training set, that were then synthesized and tested. By performing two cycles of machine learning analysis and variant design we obtained 20-fold improved proteinase K variants while only testing a total of 95 variant enzymes. CONCLUSION: The number of protein variants that must be tested to obtain significant functional improvements determines the type of tests that can be performed. Protein engineers wishing to modify the property of a protein to shrink tumours or catalyze chemical reactions under industrial conditions have until now been forced to accept high throughput surrogate screens to measure protein properties that they hope will correlate with the functionalities that they intend to modify. By reducing the number of variants that must be tested to fewer than 100, machine learning algorithms make it possible to use more complex and expensive tests so that only protein properties that are directly relevant to the desired application need to be measured. Protein design algorithms that only require the testing of a small number of variants represent a significant step towards a generic, resource-optimized protein engineering process.

Gene Designer: a synthetic biology tool for constructing artificial DNA segments.

BACKGROUND: Direct synthesis of genes is rapidly becoming the most efficient way to make functional genetic constructs and enables applications such as codon optimization, RNAi resistant genes and protein engineering. Here we introduce a software tool that drastically facilitates the design of synthetic genes. RESULTS: Gene Designer is a stand-alone software for fast and easy design of synthetic DNA segments. Users can easily add, edit and combine genetic elements such as promoters, open reading frames and tags through an intuitive drag-and-drop graphic interface and a hierarchical DNA/Protein object map. Using advanced optimization algorithms, open reading frames within the DNA construct can readily be codon optimized for protein expression in any host organism. Gene Designer also includes features such as a real-time sliding calculator of oligonucleotide annealing temperatures, sequencing primer generator, tools for avoidance or inclusion of restriction sites, and options to maximize or minimize sequence identity to a reference. CONCLUSION: Gene Designer is an expandable Synthetic Biology workbench suitable for molecular biologists interested in the de novo creation of genetic constructs.