Insights into the evolution of enzyme substrate promiscuity after the discovery of (βα)₈ isomerase evolutionary intermediates from a diverse metagenome

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

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

Taking Advantage of Promiscuity of Cold-Active Enzymes

Cold-active enzymes increase their catalytic efficiency at low-temperature, introducing structural flexibility at or near the active sites. Inevitably, this feat seems to be accompanied by lower thermal stability. These characteristics have made cold-active enzymes into attractive targets for the industrial applications, since they could reduce the energy cost in the reaction, attenuate side-reactions, and simply be inactivated. In addition, the increased structural flexibility could result in broad substrate specificity for various non-native substrates, which is called substrate promiscuity. In this perspective, we deal with a less addressed aspect of cold-active enzymes, substrate promiscuity, which has enormous potential for semi-synthesis or enzymatic modification of fine chemicals and drugs. Further structural and directed-evolutional studies on substrate promiscuity of cold-active enzymes will provide a new workhorse in white biotechnology.

How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation

Promiscuity is the coincidental ability of an enzyme to catalyze its native reaction and additional reactions that are not biological functions in the same active site. Promiscuity plays a central role in enzyme evolution and is thus a useful property for protein and metabolic engineering. This review examines enzyme evolution holistically, beginning with evaluating biochemical support for four enzyme evolution models. As expected, there is strong biochemical support for the subfunctionalization and innovation-amplification-divergence models, in which promiscuity is a central feature. In many cases, however, enzyme evolution is more complex than the models indicate, suggesting much is yet to be learned about selective pressures on enzyme function. A complete understanding of enzyme evolution must also explain the ability of metabolic networks to integrate new enzyme activities. Hidden within metabolic networks are underground metabolic pathways constructed from promiscuous activities. We discuss efforts to determine the diversity and pervasiveness of underground metabolism. Remarkably, several studies have discovered that some metabolic defects can be repaired via multiple underground routes. In prokaryotes, metabolic innovation is driven by connecting enzymes acquired by horizontal gene transfer (HGT) into the metabolic network. Thus, we end the review by discussing how the combination of promiscuity and HGT contribute to evolution of metabolism in prokaryotes. Future studies investigating the contribution of promiscuity to enzyme and metabolic evolution will need to integrate deeper probes into the influence of evolution on protein biophysics, enzymology, and metabolism with more complex and realistic evolutionary models. ENZYMES: lactate dehydrogenase (EC 1.1.1.27), malate dehydrogenase (EC 1.1.1.37), OSBS (EC 4.2.1.113), HisA (EC 5.3.1.16), TrpF, PriA (EC 5.3.1.24), R-mandelonitrile lyase (EC 4.1.2.10), Maleylacetate reductase (EC 1.3.1.32).

CiPerGenesis, A Mutagenesis Approach that Produces Small Libraries of Circularly Permuted Proteins Randomly Opened at a Focused Region: Testing on the Green Fluorescent Protein

Circularly permuted proteins (cpPs) represent a novel type of mutant proteins with original termini that are covalently linked through a peptide connector and opened at any other place of the polypeptide backbone to create new ends. cpPs are finding wide applications in biotechnology because their properties may be quite different from those of the parental protein. However, the actual challenge for the creation of successful cpPs is to identify those peptide bonds that can be broken to create new termini and ensure functional and well-folded cpPs. Herein, we describe CiPerGenesis, a combinatorial mutagenesis approach that uses two oligonucleotide libraries to amplify a circularized gene by PCR, starting and ending from a focused target region. This approach creates small libraries of circularly permuted genes that are easily cloned in the correct direction and frame using two different restriction sites encoded in the oligonucleotides. Once expressed, the protein libraries exhibit a unique sequence diversity, comprising cpPs that exhibit ordinary breakpoints between adjacent amino acids localized at the target region as well as cpPs with new termini containing user-defined truncations and repeats of some amino acids. CiPerGenesis was tested at the lid region G134-H148 of green fluorescent protein (GFP), revealing that the most fluorescent variants were those starting at Leu141 and ending at amino acids Tyr145, Tyr143, Glu142, Leu141, Lys140, and H139. Purification and biochemical characterization of some variants suggested a differential expression, solubility and maturation extent of the mutant proteins as the likely cause for the variability in fluorescence intensity observed in colonies.

Plant-based chimeric HPV-virus-like particles bearing amyloid-β epitopes elicit antibodies able to recognize amyloid plaques in APP-tg mouse and Alzheimer's disease brains

The main amyloid-beta (Aβ) variants detected in the human brain are full-length Aβ1-40 and Aβ1-42 peptides; however, a significant proportion of AD brain Aβ consists also of N-terminal truncated/modified species. The majority of the previous immunotherapeutic strategies targeted the N-terminal immunodominant epitope of the full-length Aβ; however, most of the pathological N-truncated forms of Aβ lack this critical B cell epitope. Recently, virus-like particles (VLPs), self-assembled structures with highly ordered repetitive patterns on their surface and capable of inducing robust immune responses, were applied as a promising platform for various antigen expressions. In this study, we expressed in plants two chimeric HPV16 L1 capsid proteins obtained by introduction of the β-amyloid 11-28 epitope (Aβ 11-28) into the h4 helix or into the coil regions of the L1 protein. The Aβ 11-28 epitope was chosen because it is present in the full-length Aβ 1-42 as well as in the truncated/modified amyloid peptide species. After expression, we assembled the chimerical L1/Aβ 11-28 into a VLP in which the Aβ 11-28 epitope is exposed at very high density (360 times) on the surface of the VLP. The chimeric VLPs elicited in mice Aβ-specific antibodies binding to β-amyloid plaques in APP-tg mouse and AD brains. Our study is the first to demonstrate a successful production in plants and immunogenic properties in mice of chimeric HPV16 L1 VLPs bearing Aβ epitope that may be of potential relevance for the development of multivalent vaccines for a multifactorial disease such as AD.

Evolution of substrate specificity in a retained enzyme driven by gene loss

The connection between gene loss and the functional adaptation of retained proteins is still poorly understood. We apply phylogenomics and metabolic modeling to detect bacterial species that are evolving by gene loss, with the finding that Actinomycetaceae genomes from human cavities are undergoing sizable reductions, including loss of L-histidine and L-tryptophan biosynthesis. We observe that the dual-substrate phosphoribosyl isomerase A or priA gene, at which these pathways converge, appears to coevolve with the occurrence of trp and his genes. Characterization of a dozen PriA homologs shows that these enzymes adapt from bifunctionality in the largest genomes, to a monofunctional, yet not necessarily specialized, inefficient form in genomes undergoing reduction. These functional changes are accomplished via mutations, which result from relaxation of purifying selection, in residues structurally mapped after sequence and X-ray structural analyses. Our results show how gene loss can drive the evolution of substrate specificity from retained enzymes.

DOI:http://dx.doi.org/10.7554/eLife.22679.001

Gene-specific amplicons from metagenomes as an alternative to directed evolution for enzyme screening: a case study using phenylacetaldehyde reductases

Screening gene‐specific amplicons from metagenomes (S‐GAM) is a highly promising technique for the isolation of genes encoding enzymes for biochemical and industrial applications. From metagenomes, we isolated phenylacetaldehyde reductase (par) genes, which code for an enzyme that catalyzes the production of various Prelog's chiral alcohols. Nearly full‐length par genes were amplified by PCR from metagenomic DNA, the products of which were fused with engineered par sequences at both terminal regions of the expression vector to ensure proper expression and then used to construct Escherichia coli plasmid libraries. Sequence‐ and activity‐based screening of these libraries identified different homologous par genes, Hpar‐001 to ‐036, which shared more than 97% amino acid sequence identity with PAR. Comparative characterization of these active homologs revealed a wide variety of enzymatic properties including activity, substrate specificity, and thermal stability. Moreover, amino acid substitutions in these genes coincided with those of Sar268 and Har1 genes, which were independently engineered by error‐prone PCR to exhibit increased activity in the presence of concentrated 2‐propanol. The comparative data from both approaches suggest that sequence information from homologs isolated from metagenomes is quite useful for enzyme engineering. Furthermore, by examining the GAM‐based sequence dataset derived from soil metagenomes, we easily found amino acid substitutions that increase the thermal stability of PAR/PAR homologs. Thus, GAM‐based approaches can provide not only useful homologous enzymes but also an alternative to directed evolution methodologies.

Co-occurrence of analogous enzymes determines evolution of a novel (βα)8-isomerase sub-family after non-conserved mutations in flexible loop

We investigate the evolution of co-occurring analogous enzymes involved in L-tryptophan and L-histidine biosynthesis in Actinobacteria Phylogenetic analysis of trpF homologues, a missing gene in certain clades of this lineage whose absence is complemented by a dual-substrate HisA homologue, termed PriA, found that they fall into three categories: (i) trpF-1, an L-tryptophan biosynthetic gene horizontally acquired by certain Corynebacterium species; (ii) trpF-2, a paralogue known to be involved in synthesizing a pyrrolopyrrole moiety and (iii) trpF-3, a variable non-conserved orthologue of trpF-1 We previously investigated the effect of trpF-1 upon the evolution of PriA substrate specificity, but nothing is known about the relationship between trpF-3 and priA After in vitro steady-state enzyme kinetics we found that trpF-3 encodes a phosphoribosyl anthranilate isomerase. However, mutation of this gene in Streptomyces sviceus did not lead to auxothrophy, as expected from the biosynthetic role of trpF-1 Biochemical characterization of a dozen co-occurring TrpF-2 or TrpF-3, with PriA homologues, explained the prototrophic phenotype, and unveiled an enzyme activity trade-off between TrpF and PriA. X-ray structural analysis suggests that the function of these PriA homologues is mediated by non-conserved mutations in the flexible L5 loop, which may be responsible for different substrate affinities. Thus, the PriA homologues that co-occur with TrpF-3 represent a novel enzyme family, termed PriB, which evolved in response to PRA isomerase activity. The characterization of co-occurring enzymes provides insights into the influence of functional redundancy on the evolution of enzyme function, which could be useful for enzyme functional annotation.

Increasing Metagenomic Resolution of Microbiome Interactions Through Functional Phylogenomics and Bacterial Sub-Communities

The genomic composition of the microbiome and its relationship with the environment is an exciting open question in biology. Metagenomics is a useful tool in the discovery of previously unknown taxa, but its use to understand the functional and ecological capacities of the microbiome is limited until taxonomy and function are understood in the context of the community. We suggest that this can be achieved using a combined functional phylogenomics and co-culture-based experimental strategy that can increase our capacity to measure sub-community interactions. Functional phylogenomics can identify and partition the genome such that hidden gene functions and gene clusters with unique evolutionary signals are revealed. We can test these phylogenomic predictions using an experimental model based on sub-community populations that represent a subset of the diversity directly obtained from environmental samples. These populations increase the detection of mechanisms that drive functional forces in the assembly of the microbiome, in particular the role of metabolites from key taxa in community interactions. Our combined approach leverages the potential of metagenomics to address biological questions from ecological systems.

Long-Term Persistence of Bi-functionality Contributes to the Robustness of Microbial Life through Exaptation

Modern enzymes are highly optimized biocatalysts that process their substrates with extreme efficiency. Many enzymes catalyze more than one reaction; however, the persistence of such ambiguities, their consequences and evolutionary causes are largely unknown. As a paradigmatic case, we study the history of bi-functionality for a time span of approximately two billion years for the sugar isomerase HisA from histidine biosynthesis. To look back in time, we computationally reconstructed and experimentally characterized three HisA predecessors. We show that these ancient enzymes catalyze not only the HisA reaction but also the isomerization of a similar substrate, which is commonly processed by the isomerase TrpF in tryptophan biosynthesis. Moreover, we found that three modern-day HisA enzymes from Proteobacteria and Thermotogae also possess low TrpF activity. We conclude that this bi-functionality was conserved for at least two billion years, most likely without any evolutionary pressure. Although not actively selected for, this trait can become advantageous in the case of a gene loss. Such exaptation is exemplified by the Actinobacteria that have lost the trpF gene but possess the bi-functional HisA homolog PriA, which adopts the roles of both HisA and TrpF. Our findings demonstrate that bi-functionality can perpetuate in the absence of selection for very long time-spans.

The term exaptation describes the process by which a trait that is initially just a by-product of another function may become important in a later evolutionary phase. For example, feathers served to insulate dinosaurs before helping birds fly. On the level of enzymes, bi-functionality can contribute to microbial evolution through exaptation. However, bi-functional enzymes may cause metabolic conflicts, if they are involved in different metabolic pathways. By characterizing properties of modern and computationally reconstructed ancestral variants of the sugar isomerase HisA, we demonstrate that it has been a bi-functional enzyme for the last two billion years. Most likely, bi-functionality persisted because the remaining TrpF activity is not harmful or its elimination would concurrently compromise HisA activity. Moreover, this substrate ambiguity is advantageous, as it allows compensating a gene loss as exemplified by the Actinobacteria. These microbes have lost the isomerase TrpF but possess the bi-functional HisA homolog PriA, which takes over the roles of both HisA and TrpF. Our results argue to view bi-functionality not as an evolutionary disadvantage but rather as a contribution to the evolvability of novel functions via exaptation.

Two-step Ligand Binding in a (βα)8 Barrel Enzyme: SUBSTRATE-BOUND STRUCTURES SHED NEW LIGHT ON THE CATALYTIC CYCLE OF HisA

HisA is a (βα)8 barrel enzyme that catalyzes the Amadori rearrangement of N'-[(5'-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (ProFAR) to N'-((5'-phosphoribulosyl) formimino)-5-aminoimidazole-4-carboxamide-ribonucleotide (PRFAR) in the histidine biosynthesis pathway, and it is a paradigm for the study of enzyme evolution. Still, its exact catalytic mechanism has remained unclear. Here, we present crystal structures of wild type Salmonella enterica HisA (SeHisA) in its apo-state and of mutants D7N and D7N/D176A in complex with two different conformations of the labile substrate ProFAR, which was structurally visualized for the first time. Site-directed mutagenesis and kinetics demonstrated that Asp-7 acts as the catalytic base, and Asp-176 acts as the catalytic acid. The SeHisA structures with ProFAR display two different states of the long loops on the catalytic face of the structure and demonstrate that initial binding of ProFAR to the active site is independent of loop interactions. When the long loops enclose the substrate, ProFAR adopts an extended conformation where its non-reacting half is in a product-like conformation. This change is associated with shifts in a hydrogen bond network including His-47, Asp-129, Thr-171, and Ser-202, all shown to be functionally important. The closed conformation structure is highly similar to the bifunctional HisA homologue PriA in complex with PRFAR, thus proving that structure and mechanism are conserved between HisA and PriA. This study clarifies the mechanistic cycle of HisA and provides a striking example of how an enzyme and its substrate can undergo coordinated conformational changes before catalysis.