Directed evolution of new catalytic activity using the alpha/beta-barrel scaffold

Enzyme modification using mutation site prediction method for enhancing the regioselectivity of substrate reaction sites

Enzymes with low regioselectivity of substrate reaction sites may produce multiple products from a single substrate. When a target product is produced industrially using these enzymes, the production of non-target products (byproducts) causes adverse effects such as increased processing costs for purification and the amount of raw material. Thus it is required the development of modified enzymes to reduce the amount of byproducts’ production. In this paper, we report a method called mutation site prediction for enhancing the regioselectivity of substrate reaction sites (MSPER). MSPER takes conformational data for docking poses of an enzyme and a substrate as input and automatically generates a ranked list of mutation sites to destabilize docking poses for byproducts while maintaining those for target products in silico. We applied MSPER to the enzyme cytochrome P450 CYP102A1 (BM3) and the two substrates to enhance the regioselectivity for four target products with different reaction sites. The 13 of the total 14 top-ranked mutation sites predicted by MSPER for the four target products succeeded in selectively enhancing the regioselectivity up to 6.4-fold. The results indicate that MSPER can distinguish differences of substrate structures and the reaction sites, and can accurately predict mutation sites to enhance regioselectivity without selection by directed evolution screening.

Biomimetic Molecular Clamp Nanopores for Simultaneous Quantifications of NAD+ and NADH

NADH/NAD⁺ is pivotal to fundamental biochemistry research and molecular diagnosis, but recognition and detection for them are a big challenge at the single-molecule level. Inspired by the biological system, here, we designed and synthesized a biomimetic NAD⁺/NADH molecular clamp (MC), octakis-(6-amino-6-deoxy)-γ-cyclomaltooctaose, and harbored in the engineered α-HL(M113R)₇ nanopore, forming a novel single-molecule biosensor. The single-molecule measurement possesses high selectivity and a high signal-to-noise ratio, allowing to simultaneously recognize and detect for sensing NADH/NAD⁺ and their transformations.

Innovations in CAZyme gene diversity and its modification for biorefinery applications

For sustainable growth, concept of biorefineries as recourse to the "fossil derived" energy source is important. Here, the Carbohydrate Active enZymes (CAZymes) play decisive role in generation of biofuels and related sugar-based products utilizing lignocellulose as a carbon source. Given their industrial significance, extensive studies on the evolution of CAZymes have been carried out. Various bacterial and fungal organisms have been scrutinized for the development of CAZymes, where advance techniques for strain enhancement such as CRISPR and analysis of specific expression systems have been deployed. Specific Omic-based techniques along with protein engineering have been adopted to unearth novel CAZymes and improve applicability of existing enzymes. In-Silico computational research and functional annotation of new CAZymes to synergy experiments are being carried out to devise cocktails of enzymes for use in biorefineries. Thus, with the establishment of these technologies, increased diversity of CAZymes with broad span of functions and applications is seen.

Underlying Principles of Natural Selection in Network Evolution: Systems Biology Approach

Systems biology is a rapidly expanding field that integrates diverse areas of science such as physics, engineering, computer science, mathematics, and biology toward the goal of elucidating the underlying principles of hierarchical metabolic and regulatory systems in the cell, and ultimately leading to predictive understanding of cellular response to perturbations. Because post-genomics research is taking place throughout the tree of life, comparative approaches offer a way for combining data from many organisms to shed light on the evolution and function of biological networks from the gene to the organismal level. Therefore, systems biology can build on decades of theoretical work in evolutionary biology, and at the same time evolutionary biology can use the systems biology approach to go in new uncharted directions. In this study, we present a review of how the post-genomics era is adopting comparative approaches and dynamic system methods to understand the underlying design principles of network evolution and to shape the nascent field of evolutionary systems biology. Finally, the application of evolutionary systems biology to robust biological network designs is also discussed from the synthetic biology perspective.

On the Adaptive Design Rules of Biochemical Networks in Evolution

Biochemical networks are the backbones of physiological systems of organisms. Therefore, a biochemical network should be sufficiently robust (not sensitive) to tolerate genetic mutations and environmental changes in the evolutionary process. In this study, based on the robustness and sensitivity criteria of biochemical networks, the adaptive design rules are developed for natural selection in the evolutionary process. This will provide insights into the robust adaptive mechanism of biochemical networks in the evolutionary process.

We find that if a mutated biochemical network satisfies the robustness and sensitivity criteria of natural selection, there is a high probability for the biochemical network to prevail during natural selection in the evolutionary process. Since there are various mutated biochemical networks that can satisfy these criteria but have some differences in phenotype, the biochemical networks increase their diversities in the evolutionary process. The robustness of a biochemical network enables co-option so that new phenotypes can be generated in evolution. The proposed robust adaptive design rules of natural selection gain much insight into the evolutionary mechanism and provide a systematic robust biochemical circuit design method of biochemical networks for biotechnological and therapeutic purposes in the future.

Achievements and Challenges in Computational Protein Design

Computational protein design (CPD), a yet evolving field, includes computer-aided engineering for partial or full de novo designs of proteins of interest. Designs are defined by a requested structure, function, or working environment. This chapter describes the birth and maturation of the field by presenting 101 CPD examples in a chronological order emphasizing achievements and pending challenges. Integrating these aspects presents the plethora of CPD approaches with the hope of providing a "CPD 101". These reflect on the broader structural bioinformatics and computational biophysics field and include: (1) integration of knowledge-based and energy-based methods, (2) hierarchical designated approach towards local, regional, and global motifs and the integration of high- and low-resolution design schemes that fit each such region, (3) systematic differential approaches towards different protein regions, (4) identification of key hot-spot residues and the relative effect of remote regions, (5) assessment of shape-complementarity, electrostatics and solvation effects, (6) integration of thermal plasticity and functional dynamics, (7) negative design, (8) systematic integration of experimental approaches, (9) objective cross-assessment of methods, and (10) successful ranking of potential designs. Future challenges also include dissemination of CPD software to the general use of life-sciences researchers and the emphasis of success within an in vivo milieu. CPD increases our understanding of protein structure and function and the relationships between the two along with the application of such know-how for the benefit of mankind. Applied aspects range from biological drugs, via healthier and tastier food products to nanotechnology and environmentally friendly enzymes replacing toxic chemicals utilized in the industry.

Computational enzyme design approaches with significant biological outcomes: progress and challenges

Enzymes are powerful biocatalysts, however, so far there is still a large gap between the number of enzyme-based practical applications and that of naturally occurring enzymes. Multiple experimental approaches have been applied to generate nearly all possible mutations of target enzymes, allowing the identification of desirable variants with improved properties to meet the practical needs. Meanwhile, an increasing number of computational methods have been developed to assist in the modification of enzymes during the past few decades. With the development of bioinformatic algorithms, computational approaches are now able to provide more precise guidance for enzyme engineering and make it more efficient and less laborious. In this review, we summarize the recent advances of method development with significant biological outcomes to provide important insights into successful computational protein designs. We also discuss the limitations and challenges of existing methods and the future directions that should improve them.

Metabolic engineering for production of biorenewable fuels and chemicals: contributions of synthetic biology

Production of fuels and chemicals through microbial fermentation of plant material is a desirable alternative to petrochemical-based production. Fermentative production of biorenewable fuels and chemicals requires the engineering of biocatalysts that can quickly and efficiently convert sugars to target products at a cost that is competitive with existing petrochemical-based processes. It is also important that biocatalysts be robust to extreme fermentation conditions, biomass-derived inhibitors, and their target products. Traditional metabolic engineering has made great advances in this area, but synthetic biology has contributed and will continue to contribute to this field, particularly with next-generation biofuels. This work reviews the use of metabolic engineering and synthetic biology in biocatalyst engineering for biorenewable fuels and chemicals production, such as ethanol, butanol, acetate, lactate, succinate, alanine, and xylitol. We also examine the existing challenges in this area and discuss strategies for improving biocatalyst tolerance to chemical inhibitors.

Betaalpha-hairpin clamps brace betaalphabeta modules and can make substantive contributions to the stability of TIM barrel proteins

Non-local hydrogen bonding interactions between main chain amide hydrogen atoms and polar side chain acceptors that bracket consecutive βα or αβ elements of secondary structure in αTS from E. coli, a TIM barrel protein, have previously been found to contribute 4–6 kcal mol⁻¹ to the stability of the native conformation. Experimental analysis of similar βα-hairpin clamps in a homologous pair of TIM barrel proteins of low sequence identity, IGPS from S. solfataricus and E. coli, reveals that this dramatic enhancement of stability is not unique to αTS. A survey of 71 TIM barrel proteins demonstrates a 4-fold symmetry for the placement of βα-hairpin clamps, bracing the fundamental βαβ building block and defining its register in the (βα)₈ motif. The preferred sequences and locations of βα-hairpin clamps will enhance structure prediction algorithms and provide a strategy for engineering stability in TIM barrel proteins.

Directed evolution of a lysosomal enzyme with enhanced activity at neutral pH by mammalian cell-surface display

Human beta-glucuronidase, due to low intrinsic immunogenicity in humans, is an attractive enzyme for tumor-specific prodrug activation, but its utility is hindered by low activity at physiological pH. Here we describe the development of a high-throughput screening procedure for enzymatic activity based on the stable retention of fluorescent reaction product in mammalian cells expressing properly folded glycoproteins on their surface. We utilized this procedure on error-prone PCR and saturation mutagenesis libraries to isolate beta-glucuronidase tetramers that were up to 60-fold more active (k(cat)/K(m)) at pH 7.0 and were up to an order of magnitude more effective at catalyzing the conversion of two structurally disparate glucuronide prodrugs to anticancer agents. The screening procedure described here can facilitate investigation of eukaryotic enzymes requiring posttranslational modifications for biological activity.

On the robust circuit design schemes of biochemical networks: steady-state approach

Based on the steady-state analyses of the synergism and saturation system (S-system) model, a robust control method is proposed for biochemical networks via feedback and feedforward biochemical circuits. Two robust biochemical circuit design schemes are developed. One scheme is to improve the system's structural stability so as to tolerate larger kinetic parameter variations, whereas the other is to compensate for the kinetic parameter variations to eliminate their effects. In addition, a multi-objective biochemical circuit design scheme is introduced for both the robust design against kinetic parameter variations and a desired sensitivity design to eliminate the effect of external disturbance simultaneously. The proposed robust circuit design schemes will provide a systematic method with potential applications in synthetic circuit design for biotechnological purpose and drug design purpose. Recent advances in both metabolic and genetic engineering have made the robust biochemical circuit control approach feasible through the design and implementation of synthetic biological networks amenable to mathematical modeling and quantitative analysis. Finally, several examples including the robust circuit design of the tricarboxylic acid cycle are used in silico to illustrate the design procedure and to confirm the performance of the proposed design method.

Multiobjective optimization in bioinformatics and computational biology

This paper reviews the application of multiobjective optimization in the fields of bioinformatics and computational biology. A survey of existing work, organized by application area, forms the main body of the review, following an introduction to the key concepts in multiobjective optimization. An original contribution of the review is the identification of five distinct "contexts," giving rise to multiple objectives: These are used to explain the reasons behind the use of multiobjective optimization in each application area and also to point the way to potential future uses of the technique.

Directed evolution: an approach to engineer enzymes

Directed evolution is being used increasingly in industrial and academic laboratories to modify and improve commercially important enzymes. Laboratory evolution is thought to make its biggest contribution in explorations of non-natural functions, by allowing us to distinguish the properties nurtured by evolution. In this review we report the significant advances achieved with respect to the methods of biocatalyst improvement and some critical properties and applications of the modified enzymes. The application of directed evolution has been elaborately demonstrated for protein solubility, stability and catalytic efficiency. Modification of certain enzymes for their application in enantioselective catalysis has also been elucidated. By providing a simple and reliable route to enzyme improvement, directed evolution has emerged as a key technology for enzyme engineering and biocatalysis.

Oestrogen-receptors (ER) are likely to be promiscuous: wider role for oestrogens and mimics

The anti-breast cancer drug tamoxifen that binds to ER is metabolised in human liver by CYP2D6 isoenzyme, whilst the metabolism of 17beta-oestradiol (by hydroxylation) is by phase I biotransformation in the liver to 2-hydroxyoestradiol and to 4-hydroxyoestradiol respectively by two isoenzymes of this mixed function oxidase CYP cytochromes P450 (EC 1.14.14.1); CYP1A2 and by CYP1B1. Nevertheless, it appears that the receptor (AhR) itself causes the expression of oestrogen-regulated target genes (studied by binding of dioxin). This is the result of an unknown signalling mechanism at the genome that is triggered directly by this receptor by binding promiscuously to ER (alpha or beta) sites. This has been observed even in the absence of oestrogens or mimics therefore in genome-binding investigations of target tissues such as uterus: oestrogen-receptor (ER) is likely to be promiscuous therefore. Furthermore, AhR (polycyclic aromatic hydrocarbon receptor), when activated by the binding of aromatic hydrocarbons (Ah) forms a complex with the aryl hydrocarbon nuclear-translocator chaperone protein (Arnt). It is this binding to xenobiotic response elements in DNA that initiates expression of the appropriate oestrogen-regulated target-genes in the uterus and other target tissues (including mammary, ovaries, and brain). The likely promiscuity of oestrogen receptors is proposed to be the cause of numerous side effects when oestrogen is involved in therapy, these can be manifest in hormone replacement therapy (HRT) and in the incorporation of synthetic oestrogens in the wide varieties of oral contraceptives now available.

Oxygen-induced reperfusion-injury is caused by ROS: Amelioration is possible by recombinant-DNA antioxidant enzymes and mimics in selected tissues

Billions of years of photolytic cleavage of the water molecule has led to a build-up at sea-level of a "toxic" oxygen concentration in the atmosphere (to 20%). Unfortunately, this dioxygen is likely to be converted in the mitochondria (cell organelles) during cellular respiration to generate reactive oxygen species (ROS); including free-radicals such as superoxide anion (.O2') and hydroxyl radical (.OH) (peroxide O2") a ROS is not a free radical because its electrons are paired). Development of recombinant-DNA improved isoenzyme forms (or mimics) of antioxidant enzymes such as superoxide dismutases is predicted in this hypothesis to be utilised to ameliorate reperfusion-injury (and other oxygen-induced molecular pathology). Introduction into the human genome of the genes for expression of antioxidant enzymes, to order in particular tissues, is imminent. Furthermore, O2 itself will be recognised universally as a harmful gas that can subject the cell to oxidative stress; because it produces ROS, such as the superoxide anion when it acts as the terminal electron acceptor in cellular respiration in the production of water from hydrogen ions in mitochondria. In conclusion, therefore, oxygen-induced injury in humans can no longer be accepted during medical techniques such as reperfusion procedures, because this is associated with reperfusion-injury that can be the cause of several serious medical conditions arising from biomolecular pathology. Some of this is caused by isoforms of cytochromes P450 (CYP; EC 1.14.14.1) such as 3A4, 2D6 and 2C19. These can generate ROS in the liver at low substrate concentrations by futile recycling of oxygen.

In vitro selection and characterization of a stable subdomain of phosphoribosylanthranilate isomerase

The (beta(alpha))8-barrel is the most common enzyme fold and it is capable of catalyzing an enormous diversity of reactions. It follows that this scaffold should be an ideal starting point for engineering novel enzymes by directed evolution. However, experiments to date have utilized in vivo screens or selections and the compatibility of (beta(alpha))8-barrels with in vitro selection methods remains largely untested. We have investigated plasmid display as a suitable in vitro format by engineering a variant of phosphoribosylanthranilate isomerase (PRAI) that carried the FLAG epitope in active-site-forming loop 6. Trial enrichments for binding of mAb M2 (a mAb to FLAG) demonstrated that FLAG-PRAI could be identified from a 10(6)-fold excess of a FLAG-negative competitor in three rounds of in vitro selection. These results suggest PRAI as a useful scaffold for epitope and peptide grafting experiments. Further, we constructed a FLAG-PRAI loop library of approximately 10(7) clones, in which the epitope residues most critical for binding mAb M2 were randomized. Four rounds of selection for antibody binding identified and enriched for a variant in which a single nucleotide insertion produced a truncated (beta(alpha))8-barrel consisting of (beta(alpha))1-5beta6. Biophysical characterization of this clone, trPRAI, demonstrated that it was selected because of a 21-fold increase in mAb M2 affinity compared with full-length FLAG-PRAI. Remarkably, this truncated barrel was found to be soluble, structured, thermostable and monomeric, implying that it represents a genuine subdomain of PRAI and providing further evidence that such subdomains have played an important role in the evolution of the (beta(alpha))8-barrel fold.

Modeling and Analysis of Protein Design under Resource Constraints

The potency, or fitness, of a protein-based drug can be enhanced by changing the sequence of its underlying protein. We present a novel stochastic model for the sequence-fitness relation, and estimate its four parameters from industrial data. Using this model, we formulate and analyze two variants of the protein design problem. In the single-period design problem, the designer needs to decide under capacity constraints which set of sequences to screen in order to maximize the expected fitness of the best sequence in the set. In the more general two-period design problem, the designer can afford two screening rounds and needs to allocate resources optimally across the two periods to maximize the same objective function. Analytical and simulation results allow us to assess the utility of the proposed design strategies for various parameter regimes.

p53 Protein is activated by Pin1: And also by Cu-SOD prion-like enzyme

Peptidyl-proline isomerase (Pin1) is able to trigger some conformationally important change in the p53 protein: there is notable protection by p53 (tumour suppressor protein) of human cells that prevents their entry into the carcinogenesis-committed routeway. Pin1 controls the ready (low energy change) equilibrium between the cis and trans distinctive folding configurations differentially at a proline residue: this amino acid residue in proteins is unique in bending sharply its peptide chain (to 90 degrees change): in the cis rather than trans orientation with respect to the peptide bond to residue X "upstream" linked as XCONHR. Moreover p53 protein can arrest a cell cycle progression (or trigger apoptosis) by acting as a transcription factor to nuclear DNA acting at p53 nuclear responsive element controlling a larger number of genes that produce proteins that stop cell growth or stimulate apoptosis, in stressed cells. Oxidative stress by reactive oxygen species (ROS) is carcinogenic but also stops cell growth and triggers apoptosis, Cu-SOD removes ROS (see figure). Could superoxide dismutase (Cu-SOD), therefore, provide the DNA-damage direct second route (first route is binding of Pin1) in DNA-damaged cells to p53 activation? The p53 protein that prevents carcinogenesis is activated by Pin1. In addition, this p53 tumour suppressor protein is activated by Cu-SOD.

Therapeutic-antagonists of oestrogens can be produced for cancer and other therapies using cytochromes P450 (CYP)

Oestrogens such as 17beta-oestradiol initiates nuclear-gene transcription in gender-specified tissues such as the ovaries and mammaries; and unfortunately too in cancer cells derived from target tissues. Consequently, there has been the development of novel agents for particular cancer therapies that are antagonists of oestrogens for oestrogen-receptor (ER) binding and of drugs with ER-specific interference RNA (RNAi) abilities. Therapeutic-antagonists of oestrogens will be re-designed and biosynthesised and deployed to circumvent the gene DNA-transcription abilities of oestrogens and mimics: and their metabolites in oestrogen-target tissues (see above). Furthermore, opportunities will emerge for adjunct-chemotherapy of particular tissue cancers: and in the prevention of recurrence outcomes. Cytochromes P450 can play an important part in these developments especially for the production of novel metabolites of oestrogens as therapeutic-antagonists of oestrogen-stimulated cancers.

Relation of Enzyme Activity to Local/Global Stability of Murine Adenosine Deaminase: 19F NMR Studies

Adenosine deaminase (ADA, EC 3.5.4.4) is a ubiquitous (beta/alpha)8-barrel enzyme crucial for purine metabolism and normal immune competence. In this study, it was observed that loss of enzyme activity of murine ADA (mADA) precedes the global secondary and tertiary structure transition when the protein is exposed to denaturant. The structural mechanism for this phenomenon was probed using site-specific 19F NMR spectroscopy in combination with [6-19F]tryptophan labeling and inhibitor binding. There are four tryptophan residues in mADA and all are located more than 12 A from the catalytic site. The 19F NMR spectra of [6-19F]Trp-labelled mADA show that the urea-induced chemical shift change of 19F resonance of W161, one of the four tryptophan 19F nuclei, correlates with the loss of enzyme activity. The urea-induced chemical shift change of another 19F resonance of W117 correlates with the change of the apparent rate constant for the binding of transition-state analogue inhibitor deoxycoformycin to the enzyme. On the other hand, the chemical environment of the local region around W264 does not change significantly, as a consequence of perturbation by low concentrations of urea or substrate analog. The results indicate that different regions of mADA have different local stability, which controls the activity and stability of the enzyme. The results provide new insights into the relationship between the function of a protein and its conformational flexibility as well as its global stability. This study illustrates the advantage of 19F NMR spectroscopy in probing site-related conformational change information in ligand binding, enzymatic activity and protein folding.

Aminocyclopentadienyl Ruthenium Complexes as Racemization Catalysts for Dynamic Kinetic Resolution of Secondary Alcohols at Ambient Temperature

Aminocyclopentadienyl ruthenium complexes, which can be used as room-temperature racemization catalysts with lipases in the dynamic kinetic resolution (DKR) of secondary alcohols, were synthesized from cyclopenta-2,4-dienimines, Ru(3)(CO)(12), and CHCl(3): [2,3,4,5-Ph(4)(eta(5)-C(4)CNHR)]Ru(CO)(2)Cl (4: R = i-Pr; 5: R = n-Pr; 6: R = t-Bu), [2,5-Me(2)-3,4-Ph(2)(eta(5)-C(4)CNHR)]Ru(CO)(2)Cl (7: R = i-Pr; 8: R = Ph), and [2,3,4,5-Ph(4)(eta(5)-C(4)CNHAr)]Ru(CO)(2)Cl (9: Ar = p-NO(2)C(6)H(4); 10: Ar = p-ClC(6)H(4); 11: Ar = Ph; 12: Ar = p-OMeC(6)H(4); 13: Ar = p-NMe(2)C(6)H(4)). The tests in the racemization of (S)-4-phenyl-2-butanol showed that 7 is the most active catalyst, although the difference decreased in the DKR. Complex 4 was used in the DKR of various alcohols; at room temperature, not only simple alcohols but also functionalized ones such as allylic alcohols, alkynyl alcohols, diols, hydroxyl esters, and chlorohydrins were successfully transformed to chiral acetates. In mechanistic studies for the catalytic racemization, ruthenium hydride 14 appeared to be a key species. It was the major organometallic species in the racemization of (S)-1-phenylethanol with 4 and potassium tert-butoxide. In a separate experiment, (S)-1-phenylethanol was racemized catalytically by 14 in the presence of acetophenone.

Novel cytochromes P450 applications arising from the directed-evolution of recombinant micro-organisms

AIMS

Directed (forced) evolution of cytochromes P450 (overall 2700 CYP isoforms in non-recombinant biota) is a method that has been investigated in yeasts (and other micro-organisms) by aerobically growing brewers' yeast Saccharomyces cerevisiae in very high glucose (20%) media.

METHODS AND RESULTS

Mitochondrial repression subverts cytochrome oxidase biosynthesis into manifest cytochromes P450 accumulation in brewers' yeast. A similar phenomenon is observed with the acridine-induced petit mutant. Cytochromes P450 EC 1.14.14.1 (and mimics) display a range of redox iron-mediated bioconversions in food processing, with mixed function oxidase (O2:mono-oxygenase) intervention results. Unfortunately these enzymes generate reactive oxygen species (ROS) through redox electron recycling, whilst isoform CYP 1A1 can activate precarcinogens such as benzo(a)pyrene to the ultimate (proximate) carcinogen that binds to nuclear DNA.

CONCLUSIONS

In conclusion, another 5000 CYP isoforms, for example, might be identified in micro-organisms and many more made to order through recombinant DNA technology and utilized both in vitro and in vivo for aimed bioconversions in industry and in the environment, as part of the impact of greener-approach supporting strategies to minimize global pollution.

ELISA: structure-function inferences based on statistically significant and evolutionarily inspired observations

The problem of functional annotation based on homology modeling is primary to current bioinformatics research. Researchers have noted regularities in sequence, structure and even chromosome organization that allow valid functional cross-annotation. However, these methods provide a lot of false negatives due to limited specificity inherent in the system. We want to create an evolutionarily inspired organization of data that would approach the issue of structure-function correlation from a new, probabilistic perspective. Such organization has possible applications in phylogeny, modeling of functional evolution and structural determination. ELISA (Evolutionary Lineage Inferred from Structural Analysis, http://romi.bu.edu/elisa) is an online database that combines functional annotation with structure and sequence homology modeling to place proteins into sequence-structure-function "neighborhoods". The atomic unit of the database is a set of sequences and structural templates that those sequences encode. A graph that is built from the structural comparison of these templates is called PDUG (protein domain universe graph). We introduce a method of functional inference through a probabilistic calculation done on an arbitrary set of PDUG nodes. Further, all PDUG structures are mapped onto all fully sequenced proteomes allowing an easy interface for evolutionary analysis and research into comparative proteomics. ELISA is the first database with applicability to evolutionary structural genomics explicitly in mind.

AVAILABILITY

The database is available at http://romi.bu.edu/elisa.

Ancient origin of the tryptophan operon and the dynamics of evolutionary change

The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting features that can be distinguished. As additional genomes are thoroughly analyzed, an increasingly refined resolution of the sequential evolutionary steps is clearly possible. These comparisons suggest that present-day trp operons that possess finely tuned regulatory features are under strong positive selection and are able to resist the disruptive evolutionary events that may be experienced by simpler, poorly regulated operons.

Diversity, taxonomy and evolution of medium-chain dehydrogenase/reductase superfamily

A comprehensive, structural and functional, in silico analysis of the medium-chain dehydrogenase/reductase (MDR) superfamily, including 583 proteins, was carried out by use of extensive database mining and the blastp program in an iterative manner to identify all known members of the superfamily. Based on phylogenetic, sequence, and functional similarities, the protein members of the MDR superfamily were classified into three different taxonomic categories: (a) subfamilies, consisting of a closed group containing a set of ideally orthologous proteins that perform the same function; (b) families, each comprising a cluster of monophyletic subfamilies that possess significant sequence identity among them and might share or not common substrates or mechanisms of reaction; and (c) macrofamilies, each comprising a cluster of monophyletic protein families with protein members from the three domains of life, which includes at least one subfamily member that displays activity related to a very ancient metabolic pathway. In this context, a superfamily is a group of homologous protein families (and/or macrofamilies) with monophyletic origin that shares at least a barely detectable sequence similarity, but showing the same 3D fold. The MDR superfamily encloses three macrofamilies, with eight families and 49 subfamilies. These subfamilies exhibit great functional diversity including noncatalytic members with different subcellular, phylogenetic, and species distributions. This results from constant enzymogenesis and proteinogenesis within each kingdom, and highlights the huge plasticity that MDR superfamily members possess. Thus, through evolution a great number of taxa-specific new functions were acquired by MDRs. The generation of new functions fulfilled by proteins, can be considered as the essence of protein evolution. The mechanisms of protein evolution inside MDR are not constrained to conserve substrate specificity and/or chemistry of catalysis. In consequence, MDR functional diversity is more complex than sequence diversity. MDR is a very ancient protein superfamily that existed in the last universal common ancestor. It had at least two (and probably three) different ancestral activities related to formaldehyde metabolism and alcoholic fermentation. Eukaryotic members of this superfamily are more related to bacterial than to archaeal members; horizontal gene transfer among the domains of life appears to be a rare event in modern organisms.

Structure-based rational design of streptavidin mutants with pseudo-catalytic activity

Introduction of enzymatic activity into proteins or other types of polymers by rational design is a major objective in the life sciences. To date, relatively low levels of enzymatic activity could be introduced into antibodies by using transition-state analogues of haptens. In the present study, we identify the structural elements that contribute to the observed hydrolytic activity in egg white avidin, which promote the cleavage of active biotin esters (notably biotinyl p-nitrophenyl ester). The latter elements were then incorporated into bacterial streptavidin via genetic engineering. The streptavidin molecule was thus converted from a protector to an enhancer of hydrolysis of biotin esters. The conversion was accomplished by the combined replacement of a "lid-like loop" (L3,4) and a leucine-to-arginine point mutation in streptavidin. Interestingly, neither of these elements play a direct role in the hydrolytic reaction. The latter features were thus shown to be responsible for enhanced substrate hydrolysis. This work indicates that structural and non-catalytic elements of a protein can be modified to promote the induced fit of a substrate for subsequent interaction with either a catalytic residue or water molecules. This approach complements the conventional design of active sites that involves direct modifications of catalytic residues.

De novo backbone and sequence design of an idealized alpha/beta-barrel protein: evidence of stable tertiary structure

We have designed, synthesized, and characterized a 216 amino acid residue sequence encoding a putative idealized alpha/beta-barrel protein. The design was elaborated in two steps. First, the idealized backbone was defined with geometric parameters representing our target fold: a central eight parallel-stranded beta-sheet surrounded by eight parallel alpha-helices, connected together with short structural turns on both sides of the barrel. An automated sequence selection algorithm, based on the dead-end elimination theorem, was used to find the optimal amino acid sequence fitting the target structure. A synthetic gene coding for the designed sequence was constructed and the recombinant artificial protein was expressed in bacteria, purified and characterized. Far-UV CD spectra with prominent bands at 222nm and 208nm revealed the presence of alpha-helix secondary structures (50%) in fairly good agreement with the model. A pronounced absorption band in the near-UV CD region, arising from immobilized aromatic side-chains, showed that the artificial protein is folded in solution. Chemical unfolding monitored by tryptophan fluorescence revealed a conformational stability (DeltaG(H2O)) of 35kJ/mol. Thermal unfolding monitored by near-UV CD revealed a cooperative transition with an apparent T(m) of 65 degrees C. Moreover, the artificial protein did not exhibit any affinity for the hydrophobic fluorescent probe 1-anilinonaphthalene-8-sulfonic acid (ANS), providing additional evidence that the artificial barrel is not in the molten globule state, contrary to previously designed artificial alpha/beta-barrels. Finally, 1H NMR spectra of the folded and unfolded proteins provided evidence for specific interactions in the folded protein. Taken together, the results indicate that the de novo designed alpha/beta-barrel protein adopts a stable three-dimensional structure in solution. These encouraging results show that de novo design of an idealized protein structure of more than 200 amino acid residues is now possible, from construction of a particular backbone conformation to determination of an amino acid sequence with an automated sequence selection algorithm.

Structure-based mutagenesis approaches toward expanding the substrate specificity of D-2-deoxyribose-5-phosphate aldolase

2-Deoxyribose-5-phosphate aldolase (DERA, EC 4.1.2.4) catalyzes the reversible aldol reaction between acetaldehyde and D-glyceraldehyde-3-phosphate to generate D-2-deoxyribose-5-phosphate. It is unique among the aldolases as it catalyzes the reversible asymmetric aldol addition reaction of two aldehydes. In order to expand the substrate scope and stereoselectivity of DERA, structure-based substrate design as well as site-specific mutation has been investigated. Using the 1.05 A crystal structure of DERA in complex with its natural substrate as a guide, five site-directed mutants were designed in order to improve its activity with the unnatural nonphosphorylated substrate, D-2-deoxyribose. Of these, the S238D variant exhibited a 2.5-fold improvement over the wild-type enzyme in the retroaldol reaction of 2-deoxyribose. Interestingly, this S238D mutant enzyme was shown to accept 3-azidopropinaldehyde as a substrate in a sequential asymmetric aldol reaction to form a deoxy-azidoethyl pyranose, which is a precursor to the corresponding lactone and the cholesterol-lowering agent Lipitor. This azidoaldehyde is not a substrate for the wild-type enzyme. Another structure-based design of new nonphosphorylated substrates was focused on the aldol reaction with inversion in enantioselectivity using the wild type or the S238D variant as the catalyst and 2-methyl-substituted aldehydes as substrates. An example was demonstrated in the asymmetric synthesis of a deoxypyranose as a new effective synthon for the total synthesis of epothilones. In addition, to facilitate the discovery of new enzymatic reactions, the engineered E. coli strain SELECT (Deltaace, adhC, DE3) was developed to be used in the future for selection of DERA variants with novel nonphosphorylated acceptor specificity.

Enhancement of transglycosylation activity by construction of chimeras between mesophilic and thermophilic beta-glucosidase

The family 3 beta-glucosidase from Thermotoga maritima is a highly thermostable enzyme (85 degrees C) that displays transglycosylation activity. In contrast, the beta-glucosidase from Cellvibrio gilvus is mesophilic (35 degrees C) and displays no such transglycosylation activity. Both enzymes consist of two domains, an N-terminal and a C-terminal domain, and the amino acid identities between the two enzymes in these domains are 32.4 and 36.4%, respectively. In an attempt to identify the molecular basis underpinning the display of transglycosylation activity and the requirements for thermal stability, eight chimeric genes were constructed by shuffling the two parental beta-glucosidase genes at four selected borders, two in the N-terminal domain and two in the C-terminal domain. Of the eight chimeric genes constructed, only two chimeric enzymes (Tm578/606Cg and Tm638/666Cg) gave catalytically active forms and these were the ones shuffled in the C-terminal domain. For these active chimeric enzymes, 80% (Tm578/606Cg) and 88% (Tm638/666Cg) of their amino acid sequences originated from T. maritima. With regard to their thermal profiles, the two active chimeric enzymes, Tm578/606Cg and Tm638/666Cg, displayed profiles intermediate to those of the two parental enzymes as they were optimally active at 65 and 70 degrees C, respectively. These two chimeric enzymes were optimally active at pH 4.1 and 3.9, which is closer to that observed for the T. maritima enzyme (pH 3.2-3.5) than that for the C. gilvus enzyme (pH 6.2-6.5). Kinetic parameters for the chimeric enzymes were investigated with five different substrates including pNP-beta-D-glucopyranoside. The kinetic parameters obtained for the chimeric enzymes were closer to those of the T. maritima enzyme than to those of the C. gilvus enzyme. Transglycosylation activity was observed for both chimeric enzymes and the activity of the Tm578/606Cg chimera was at a level twice that observed with the T. maritima enzyme. This study is an effective demonstration of the usefulness of chimeric enzymes in altering the characteristics of an enzyme.

In vitro selection for catalytic activity with ribosome display

We report what is, to our knowledge, the first in vitro selection for catalytic activity based on catalytic turnover by using ribosome display, a method which does not involve living cells at any step. RTEM-beta-lactamase was functionally displayed on ribosomes as a complex with its encoding mRNA. We designed and synthesized a mechanism-based inhibitor of beta-lactamase, biotinylated ampicillin sulfone, appropriate for selection of catalytic activity of the ribosome-displayed beta-lactamase. This derivative of ampicillin inactivated beta-lactamase in a specific and irreversible manner. Under appropriate selection conditions, active RTEM-beta-lactamase was enriched relative to an inactive point mutant over 100-fold per ribosome display selection cycle. Selection for binding, carried out with beta-lactamase inhibitory protein (BLIP), gave results similar to selection with the suicide inhibitor, indicating that ribosome display is similarly efficient in catalytic activity and affinity selections. In the future, the capacity to select directly for enzymatic activity using an entirely in vitro process may allow for a significant increase in the explorable sequence space relative to existing strategies.

Industrial enzyme applications

The effective catalytic properties of enzymes have already promoted their introduction into several industrial products and processes. Recent developments in biotechnology, particularly in areas such as protein engineering and directed evolution, have provided important tools for the efficient development of new enzymes. This has resulted in the development of enzymes with improved properties for established technical applications and in the production of new enzymes tailor-made for entirely new areas of application where enzymes have not previously been used.

Plasticity of enzyme active sites

The expectation is that any similarity in reaction chemistry shared by enzyme homologues is mediated by common functional groups conserved through evolution. However, detailed enzyme studies have revealed the flexibility of many active sites, in that different functional groups, unconserved with respect to position in the primary sequence, mediate the same mechanistic role. Nevertheless, the catalytic atoms might be spatially equivalent. More rarely, the active sites have completely different locations in the protein scaffold. This variability could result from: (1) the hopping of functional groups from one position to another to optimize catalysis; (2) the independent specialization of a low-activity primordial enzyme in different phylogenetic lineages; (3) functional convergence after evolutionary divergence; or (4) circular permutation events.

In vitro evolution of recognition specificity mediated by SH3 domains reveals target recognition rules

We have designed a repertoire of 10(7) different SH3 domains by grafting the residues that are represented in the binding surfaces of natural SH3 domains onto the scaffold of the human Abl-SH3 domain. This phage-displayed library was screened by affinity selection for SH3 domains that bind to the synthetic peptides, APTYPPPLPP and LSSRPLPTLPSP, which are peptide ligands for the human Abl or Src SH3 domains, respectively. By characterizing the isolates, we have observed that as few as two or three amino acid substitutions lead to dramatic changes in recognition specificity. We propose that the ability to shift recognition specificity with a small number of amino acid replacements is an important evolutionary characteristic of protein binding modules. Furthermore, we have used the information obtained by these in vitro evolution experiments to generate a scoring matrix that evaluates the probability that any SH3 domain binds to the peptide ligands for the Abl and Src SH3 domains. A table of predictions for the 28 SH3 domains of baker's yeast is presented.

Directed evolution of enzymes and pathways for industrial biocatalysis

Directed evolution has become a powerful tool for developing enzyme and whole cell based biocatalysts. Significant recent advances include the creation of novel enzyme functions and the development of several new efficient directed evolution methods. The combination of directed evolution and rational design promises to accelerate the development of biocatalysts for applications in the pharmaceutical, chemical and food industries.

De novo design of biocatalysts

The challenging field of de novo enzyme design is beginning to produce exciting results. The application of powerful computational methods to functional protein design has recently succeeded at engineering target activities. In addition, efforts in directed evolution continue to expand the transformations that can be accomplished by existing enzymes. The engineering of completely novel catalytic activity requires traversing inactive sequence space in a fitness landscape, a feat that is better suited to computational design. Optimizing activity, which can include subtle alterations in backbone conformation and protein motion, is better suited to directed evolution, which is highly effective at scaling fitness landscapes towards maxima. Improved rational design efforts coupled with directed evolution should dramatically improve the scope of de novo enzyme design.

Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies

The protein sequence and structure databases are now sufficiently representative that strategies nature uses to evolve new catalytic functions can be identified. Groups of divergently related enzymes whose members catalyze different reactions but share a common partial reaction, intermediate, or transition state (mechanistically diverse superfamilies) have been discovered, including the enolase, amidohydrolase, thiyl radical, crotonase, vicinal-oxygen-chelate, and Fe-dependent oxidase superfamilies. Other groups of divergently related enzymes whose members catalyze different overall reactions that do not share a common mechanistic strategy (functionally distinct suprafamilies) have also been identified: (a) functionally distinct suprafamilies whose members catalyze successive transformations in the tryptophan and histidine biosynthetic pathways and (b) functionally distinct suprafamilies whose members catalyze different reactions in different metabolic pathways. An understanding of the structural bases for the catalytic diversity observed in super- and suprafamilies may provide the basis for discovering the functions of proteins and enzymes in new genomes as well as provide guidance for in vitro evolution/engineering of new enzymes.

The structure of 3-methylaspartase from Clostridium tetanomorphum functions via the common enolase chemical step

Methylaspartate ammonia-lyase (3-methylaspartase, MAL; EC ) catalyzes the reversible anti elimination of ammonia from L-threo-(2S,3S)-3-methylaspartic acid to give mesaconic acid. This reaction lies on the main catabolic pathway for glutamate in Clostridium tetanomorphum. MAL requires monovalent and divalent cation cofactors for full catalytic activity. The enzyme has attracted interest because of its potential use as a biocatalyst. The structure of C. tetanomorphum MAL has been solved to 1.9-A resolution by the single-wavelength anomalous diffraction method. A divalent metal ion complex of the protein has also been determined. MAL is a homodimer with each monomer consisting of two domains. One is an alpha/beta-barrel, and the other smaller domain is mainly beta-strands. The smaller domain partially occludes the C terminus of the barrel and forms a large cleft. The structure identifies MAL as belonging to the enolase superfamily of enzymes. The metal ion site is located in a large cleft between the domains. Potential active site residues have been identified based on a combination of their proximity to a metal ion site, molecular modeling, and sequence homology. In common with all members of the enolase superfamily, the carboxylic acid of the substrate is co-ordinated by the metal ions, and a proton adjacent to a carboxylic acid group of the substrate is abstracted by a base. In MAL, it appears that Lys(331) removes the alpha-proton of methylaspartic acid. This motif is the defining mechanistic characteristic of the enolase superfamily of which all have a common fold. The degree of structural conservation is remarkable given only four residues are absolutely conserved.

Stabilization of a (betaalpha)8-barrel protein by an engineered disulfide bridge

The aim of this study was to increase the stability of the thermolabile (betaalpha)8-barrel enzyme indoleglycerol phosphate synthase from Escherichia coli by the introduction of disulfide bridges. For the design of such variants, we selected two out of 12 candidates, in which newly introduced cysteines potentially form optimal disulfide bonds. These variants avoid short-range connections, substitutions near catalytic residues, and crosslinks between the new and the three parental cysteines. The variant linking residues 3 and 189 fastens the N-terminus to the (betaalpha)8-barrel. The rate of thermal inactivation at 50 degrees C of this variant with a closed disulfide bridge is 65-fold slower than that of the reference dithiol form, but only 13-fold slower than that of the parental protein. The near-ultraviolet CD spectrum, the reactivity of parental buried cysteines with Ellman's reagent as well as the decreased turnover number indicate that the protein structure is rigidified. To confirm these data, we have solved the X-ray structure to 2.1-A resolution. The second variant was designed to crosslink the terminal modules betaalpha1 and betaalpha8. However, not even the dithiol form acquired the native fold, possibly because one of the targeted residues is solvent-inaccessible in the parental protein.

Recombinatoric exploration of novel folded structures: a heteropolymer-based model of protein evolutionary landscapes

The role of recombination in evolution is compared with that of point mutations (substitutions) in the context of a simple, polymer physics-based model mapping between sequence (genotype) and conformational (phenotype) spaces. Crossovers and point mutations of lattice chains with a hydrophobic polar code are investigated. Sequences encoding for a single ground-state conformation are considered viable and used as model proteins. Point mutations lead to diffusive walks on the evolutionary landscape, whereas crossovers can "tunnel" through barriers of diminished fitness. The degree to which crossovers allow for more efficient sequence and structural exploration depends on the relative rates of point mutations versus that of crossovers and the dispersion in fitness that characterizes the ruggedness of the evolutionary landscape. The probability that a crossover between a pair of viable sequences results in viable sequences is an order of magnitude higher than random, implying that a sequence's overall propensity to encode uniquely is embodied partially in local signals. Consistent with this observation, certain hydrophobicity patterns are significantly more favored than others among fragments (i.e., subsequences) of sequences that encode uniquely, and examples reminiscent of autonomous folding units in real proteins are found. The number of structures explored by both crossovers and point mutations is always substantially larger than that via point mutations alone, but the corresponding numbers of sequences explored can be comparable when the evolutionary landscape is rugged. Efficient structural exploration requires intermediate nonextreme ratios between point-mutation and crossover rates.

A functional role for a flexible loop containing Glu182 in the class II fructose-1,6-bisphosphate aldolase from Escherichia coli

Class II fructose 1,6-bisphosphate aldolases (FBP-aldolases) catalyse the zinc-dependent, reversible aldol condensation of dihydroxyacetone phosphate (DHAP) and glyceraldehyde 3-phosphate (G3P) to form fructose 1,6-bisphosphate (FBP). Analysis of the structure of the enzyme from Escherichia coli in complex with a transition state analogue (phosphoglycolohydroxamate, PGH) suggested that substrate binding caused a conformational change in the beta5-alpha7 loop of the enzyme and that this caused the relocation of two glutamate residues (Glu181 and Glu182) into the proximity of the active site. Site-directed mutagenesis of these two glutamate residues (E181A and E182A) along with another active site glutamate (Glu174) was carried out and the mutant enzymes characterised using steady-state kinetics. Mutation of Glu174 (E174A) resulted in an enzyme which was severely crippled in catalysis, in agreement with its position as a zinc ligand in the enzyme's structure. The E181A mutant showed the same properties as the wild-type enzyme indicating that the residue played no major role in substrate binding or enzyme catalysis. In contrast, mutation of Glu182 (E182A) demonstrated that Glu182 is important in the catalytic cycle of the enzyme. Furthermore, the measurement of deuterium kinetic isotope effects using [1(S)-(2)H]DHAP showed that, for the wild-type enzyme, proton abstraction was not the rate determining step, whereas in the case of the E182A mutant this step had become rate limiting, providing evidence for the role of Glu182 in abstraction of the C1 proton from DHAP in the condensation direction of the reaction. Glu182 lies in a loop of polypeptide which contains four glycine residues (Gly176, Gly179, Gly180 and Gly184) and a quadruple mutant (where each glycine was converted to alanine) showed that flexibility of this loop was important for the correct functioning of the enzyme, probably to change the microenvironment of Glu182 in order to perturb its pK(a) to a value suitable for its role in proton abstraction. These results highlight the need for further studies of the dynamics of the enzyme in order to fully understand the complexities of loop closure and catalysis in this enzyme.