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

BACKGROUND

Current sequence-based approaches to identify enzyme functional shifts, such as enzyme promiscuity, have proven to be highly dependent on a priori functional knowledge, hampering our ability to reconstruct evolutionary history behind these mechanisms. Hidden Markov Model (HMM) profiles, broadly used to classify enzyme families, can be useful to distinguish between closely related enzyme families with different specificities. The (βα)8-isomerase HisA/PriA enzyme family, involved in L-histidine (HisA, mono-substrate) biosynthesis in most bacteria and plants, but also in L-tryptophan (HisA/TrpF or PriA, dual-substrate) biosynthesis in most Actinobacteria, has been used as model system to explore evolutionary hypotheses and therefore has a considerable amount of evolutionary, functional and structural knowledge available. We searched for functional evolutionary intermediates between the HisA and PriA enzyme families in order to understand the functional divergence between these families.

RESULTS

We constructed a HMM profile that correctly classifies sequences of unknown function into the HisA and PriA enzyme sub-families. Using this HMM profile, we mined a large metagenome to identify plausible evolutionary intermediate sequences between HisA and PriA. These sequences were used to perform phylogenetic reconstructions and to identify functionally conserved amino acids. Biochemical characterization of one selected enzyme (CAM1) with a mutation within the functionally essential N-terminus phosphate-binding site, namely, an alanine instead of a glycine in HisA or a serine in PriA, showed that this evolutionary intermediate has dual-substrate specificity. Moreover, site-directed mutagenesis of this alanine residue, either backwards into a glycine or forward into a serine, revealed the robustness of this enzyme. None of these mutations, presumably upon functionally essential amino acids, significantly abolished its enzyme activities. A truncated version of this enzyme (CAM2) predicted to adopt a (βα)6-fold, and thus entirely lacking a C-terminus phosphate-binding site, was identified and shown to have HisA activity.

CONCLUSION

As expected, reconstruction of the evolution of PriA from HisA with HMM profiles suggest that functional shifts involve mutations in evolutionarily intermediate enzymes of otherwise functionally essential residues or motifs. These results are in agreement with a link between promiscuous enzymes and intragenic epistasis. HMM provides a convenient approach for gaining insights into these evolutionary processes.

Alternative splice variants in TIM barrel proteins from human genome correlate with the structural and evolutionary modularity of this versatile protein fold

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)₈ barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)₂, (βα)₄, (βα)₅, (βα)₆, and (βα)₇ subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.

Protein design through systematic catalytic loop exchange in the (beta/alpha)8 fold

Protein engineering by directed evolution has proven effective in achieving various functional modifications, but the well-established protocols for the introduction of variability, typically limited to random point mutations, seriously restrict the scope of the approach. In an attempt to overcome this limitation, we sought to explore variant libraries with richer diversity at regions recognized as functionally important through an exchange of natural components, thus combining design with combinatorial diversity. With this approach, we expected to maintain interactions important for protein stability while directing the introduction of variability to areas important for catalysis. Our strategy consisted in loop exchange over a (beta/alpha)(8) fold. Phosphoribosylanthranilate isomerase was chosen as scaffold, and we investigated its tolerance to loop exchange by fusing variant libraries to the chloramphenicol acetyl transferase coding gene as an in vivo folding reporter. We replaced loops 2, 4, and 6 of phosphoribosylanthranilate isomerase with loops of varied types and sizes from enzymes sharing the same fold. To allow for a better structural fit, saturation mutagenesis was adopted at two amino acid positions preceding the exchanged loop. Our results showed that 30% to 90% of the generated mutants in the different libraries were folded. Some variants were selected for further characterization after removal of chloramphenicol acetyl transferase gene, and their stability was studied by circular dichroism and fluorescence spectroscopy. The sequences of 545 clones show that the introduction of variability at "hinges" connecting the loops with the scaffold exhibited a noticeable effect on the appearance of folded proteins. Also, we observed that each position accepted foreign loops of different sizes and sequences. We believe our work provides the basis of a general method of exchanging variably sized loops within the (beta/alpha)(8) fold, affording a novel starting point for the screening of novel activities as well as modest diversions from an original activity.

Generation of variability by in vivo recombination of halves of a (beta/alpha)8 barrel protein

Similar to what has been achieved with nucleic acids, directed evolution of proteins would be greatly facilitated by the availability of large libraries and efficient selection methods. So far, host cell transformation efficiency has been a bottleneck, practically limiting libraries to sizes less than 10(9). One way to circumvent this problem has been implemented with antibody systems, where contribution to the binding site is provided by two different polypeptides (light and heavy chains). The central concept is the construction of binary systems in which the gene from the two chains are separated by a cre-lox recombinase recognition site, packaged in a phage, and subsequently introduced, by multiple infection, into a recombinase expressing cell [Sblattero D, Bradbury A. Nat Biotechnol 2000;18(1):75-80]. Here, we describe the development of a system which applies the same concept to a single-domain enzyme, the cytoplasmic (beta/alpha)8 barrel protein phosphoribosyl anthranilate isomerase (PRAI) from E. coli. For that purpose, we identified the site at which a loop containing the recognition sequence for cre-lox recombinase could be inserted yielding a functional enzyme. We evaluated the effect of this insertion on the capability of the engineered gene to complement a trp F-E. coli strain and the efficiency of the system to recover the original sequence from an abundance of non-functional mutant genes.

Convergent evolution of hydroxylation mechanisms in the fungal kingdom: molybdenum cofactor-independent hydroxylation of xanthine via α-ketoglutarate-dependent dioxygenases

The xanthine oxidases and dehydrogenases are among the most conserved enzymes in all living kingdoms. They contain the molybdopterin cofactor Moco. We show here that in the fungi, in addition to xanthine dehydrogenase, a completely different enzyme is able to catalyse the oxidation of xanthine to uric acid. In Aspergillus nidulans this enzyme is coded by the xanA gene. We have cloned the xanA gene and determined its sequence. A deletion of the gene has the same phenotype as the previously known xanA1 miss-sense mutation. Homologues of xanA exist only in the fungal kingdom. We have inactivated the cognate gene of Schizosaccharomyces pombe and this results in strongly impaired xanthine utilization as a nitrogen source. We have shown that the Neurospora crassa homologue is functionally equivalent to xanA. The enzyme coded by xanA is an alpha-ketoglutarate- and Fe(II)-dependent dioxygenase which shares a number of properties with other enzymes of this group. This work shows that only in the fungal kingdom, an alternative mechanism of xanthine oxidation, not involving Moco, has evolved using the dioxygenase scaffold.