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

Structural, kinetic, and evolutionary peculiarities of HISN3, a plant 5'-ProFAR isomerase

Histidine biosynthesis is essential for the growth and development of plants, where it occurs within chloroplasts. The eleven reactions are catalyzed by eight enzymes, known as HISN1-8, each acting sequentially. Here, we present the crystal structures of a 5'-ProFAR isomerase (HISN3) from the model legume Medicago truncatula bound to its enzymatically synthesized substrate (ProFAR) and product (PrFAR). The active site of MtHISN3 contains a sodium cation that participates in ligand recognition, a feature not observed in bacterial and fungal structures of homologous enzymes. The steady-state kinetics of wild-type MtHISN3 revealed a slightly higher turnover rate compared to its bacterial homologs. Plant HISN3 sequences contain an unusually elongated Lys60-Ser91 fragment, while deletion of the 74-80 region resulted in a 30-fold loss in catalytic efficiency compared to the wild-type. Molecular dynamics simulations suggested that the fragment facilitates product release, thereby contributing to a higher kcat. Moreover, conservation analyses suggested a non-cyanobacterial origin for plant HISN3 enzymes, which is another instance of a non-cyanobacterial enzyme in the plant histidine biosynthetic pathway. Finally, a virtual screening campaign yielded five molecules, with the energy gains ranging between -13.6 and -13.1 kcal/mol, which provide new scaffolds for the future development of herbicides.

Conformational Modulation of a Mobile Loop Controls Catalysis in the (βα)8-Barrel Enzyme of Histidine Biosynthesis HisF

The overall significance of loop motions for enzymatic activity is generally accepted. However, it has largely remained unclear whether and how such motions can control different steps of catalysis. We have studied this problem on the example of the mobile active site β1α1-loop (loop1) of the (βα)8-barrel enzyme HisF, which is the cyclase subunit of imidazole glycerol phosphate synthase. Loop1 variants containing single mutations of conserved amino acids showed drastically reduced rates for the turnover of the substrates N'-[(5'-phosphoribulosyl) formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) and ammonia to the products imidazole glycerol phosphate (ImGP) and 5-aminoimidazole-4-carboxamide-ribotide (AICAR). A comprehensive mechanistic analysis including stopped-flow kinetics, X-ray crystallography, NMR spectroscopy, and molecular dynamics simulations detected three conformations of loop1 (open, detached, closed) whose populations differed between wild-type HisF and functionally affected loop1 variants. Transient stopped-flow kinetic experiments demonstrated that wt-HisF binds PrFAR by an induced-fit mechanism whereas catalytically impaired loop1 variants bind PrFAR by a simple two-state mechanism. Our findings suggest that PrFAR-induced formation of the closed conformation of loop1 brings active site residues in a productive orientation for chemical turnover, which we show to be the rate-limiting step of HisF catalysis. After the cyclase reaction, the closed loop conformation is destabilized, which favors the formation of detached and open conformations and hence facilitates the release of the products ImGP and AICAR. Our data demonstrate how different conformations of active site loops contribute to different catalytic steps, a finding that is presumably of broad relevance for the reaction mechanisms of (βα)8-barrel enzymes and beyond.

Structural Outlier Detection and Zernike-Canterakis Moments for Molecular Surface Meshes-Fast Implementation in Python

Object retrieval systems measure the degree of similarity of the shape of 3D models. They search for the elements of the 3D model databases that resemble the query model. In structural bioinformatics, the query model is a protein tertiary/quaternary structure and the objective is to find similarly shaped molecules in the Protein Data Bank. With the ever-growing size of the PDB, a direct atomic coordinate comparison with all its members is impractical. To overcome this problem, the shape of the molecules can be encoded by fixed-length feature vectors. The distance of a protein to the entire PDB can be measured in this low-dimensional domain in linear time. The state-of-the-art approaches utilize Zernike-Canterakis moments for the shape encoding and supply the retrieval process with geometric data of the input structures. The BioZernike descriptors are a standard utility of the PDB since 2020. However, when trying to calculate the ZC moments locally, the issue of the deficiency of libraries readily available for use in custom programs (i.e., without relying on external binaries) is encountered, in particular programs written in Python. Here, a fast and well-documented Python implementation of the Pozo-Koehl algorithm is presented. In contrast to the more popular algorithm by Novotni and Klein, which is based on the voxelized volume, the PK algorithm produces ZC moments directly from the triangular surface meshes of 3D models. In particular, it can accept the molecular surfaces of proteins as its input. In the presented PK-Zernike library, owing to Numba's just-in-time compilation, a mesh with 50,000 facets is processed by a single thread in a second at the moment order 20. Since this is the first time the PK algorithm is used in structural bioinformatics, it is employed in a novel, simple, but efficient protein structure retrieval pipeline. The elimination of the outlying chain fragments via a fast PCA-based subroutine improves the discrimination ability, allowing for this pipeline to achieve an 0.961 area under the ROC curve in the BioZernike validation suite (0.997 for the assemblies). The correlation between the results of the proposed approach and of the 3D Surfer program attains values up to 0.99.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

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.

Simulations reveal the key role of Arg15 in the promiscuous activity in the HisA enzyme

The HisA enzyme catalyzes the first step of histidine biosynthesis via the Amadori rearrangement of the substrate ProFAR. Since it possesses the most conserved and ancient TIM-barrel fold, it provides an ideal framework for bioengineering of a new function from ancestral enzymes. In the present study, first, the catalytic mechanism of HisA biosynthesis was elucidated using hybrid Quantum Mechanical/Molecular Mechanical calculations, and thereafter, key residues contributing towards the promiscuity for TrpF activity were revealed using several MD simulations of a wild type enzyme and its variant with the native (ProFAR) and promiscuous (PRA) substrates. Our study reveals that the two loops (βα)₁ and (βα)₅ on the catalytic site of the HisA enzyme have incredible adaptability for the native and promiscuous substrates. The conformational interplay between these two loops is substrate driven and precise bioengineering targeting these loops is key to the emergence of new functions. Furthermore, the study reveals a key role of the Arg 15 residue which is close to the catalytic center of the enzyme in the bifunctionality of the HisA enzyme by increasing the loop flexibility. Therefore, our study provides crucial information for future bioengineering work to use the HisA enzyme as a scaffold for new enzymatic activity.

Evolution of new enzymes by gene duplication and divergence

Thousands of new metabolic and regulatory enzymes have evolved by gene duplication and divergence since the dawn of life. New enzyme activities often originate from promiscuous secondary activities that have become important for fitness due to a change in the environment or a mutation. Mutations that make a promiscuous activity physiologically relevant can occur in the gene encoding the promiscuous enzyme itself, but can also occur elsewhere, resulting in increased expression of the enzyme or decreased competition between the native and novel substrates for the active site. If a newly useful activity is inefficient, gene duplication/amplification will set the stage for divergence of a new enzyme. Even a few mutations can increase the efficiency of a new activity by orders of magnitude. As efficiency increases, amplified gene arrays will shrink to provide two alleles, one encoding the original enzyme and one encoding the new enzyme. Ultimately, genomic rearrangements eliminate co-amplified genes and move newly evolved paralogs to a distant region of the genome.

Mutational Pathways and Trade-Offs Between HisA and TrpF Functions: Implications for Evolution via Gene Duplication and Divergence

When a new activity evolves by changes in a pre-existing enzyme this is likely to reduce the original activity, generating a functional trade-off. The properties of this trade-off will affect the continued evolution of both functions. If the trade-off is strong, gene duplication and subsequent divergence would be favored whereas if the trade-off is weak a bi-functional enzyme could evolve that performs both functions. We previously showed that when a bi-functional HisA enzyme was evolved under selection for both HisA and TrpF functions, evolution mainly proceeded via duplication-divergence and specialization, implying that the trade-off is strong between these two functions. Here, we examined this hypothesis by identifying the mutational pathways (i.e., the mutational landscape) in the Salmonella enterica HisA enzyme that conferred a TrpF-like activity, and examining the trade-offs between the original and new activity. For the HisA enzyme there are many different paths toward the new TrpF function, each with its own unique trade-off. A total of 16 single mutations resulted in HisA enzyme variants that acquired TrpF activity and only three of them maintained HisA activity. Twelve mutants were evolved further toward increased TrpF activity and during evolution toward improved TrpF activity the original HisA activity was completely lost in all lineages. We propose that, aside from various relevant ecological factors, two main genetic factors influence whether evolution of a new function proceeds via duplication – divergence (specialization) or by evolution of a generalist: (i) the relative mutation supply of the two pathways and (ii) the shape of the trade-off curve between the native and new function.

Identification and Analysis of Natural Building Blocks for Evolution-Guided Fragment-Based Protein Design

Natural evolution has generated an impressively diverse protein universe via duplication and recombination from a set of protein fragments that served as building blocks. The application of these concepts to the design of new proteins using subdomain-sized fragments from different folds has proven to be experimentally successful. To better understand how evolution has shaped our protein universe, we performed an all-against-all comparison of protein domains representing all naturally existing folds and identified conserved homologous protein fragments. Overall, we found more than 1000 protein fragments of various lengths among different folds through similarity network analysis. These fragments are present in very different protein environments and represent versatile building blocks for protein design. These data are available in our web server called F(old P)uzzle (fuzzle.uni-bayreuth.de), which allows to individually filter the dataset and create customized networks for folds of interest. We believe that our results serve as an invaluable resource for structural and evolutionary biologists and as raw material for the design of custom-made proteins.

Structural analysis of Shigella flexneri bi-functional enzyme HisIE in histidine biosynthesis

Histidine biosynthesis, which is absent in animals, was shown to be highly conserved among gram-negative bacteria, thus making it an attractive target for antibiotic design. There are many fusion forms of enzymes in the histidine biosynthetic pathway and people still have limited knowledge about their domain organizations and catalytic mechanisms, due to the lack of structural information. Here we report the first crystal structure of Shigella flexneri bi-functional enzyme HisIE (SfHisIE) that functions in the 2nd and 3rd steps in the histidine biosynthetic pathway. This structure shows that HisIE exists as dimers with two loops (fusion loop) connecting the individual dimer of HisE and HisI in its N-terminus and C-terminus respectively. Our mutagenesis study shows mutations in this fusion loop are lethal for bacteria indicating the advantage of gene fusion in Histidine biosynthesis. Structural analysis revealed several highly conserved residues in the putative ligand binding grooves of HisE and HisI, showing an evolutionarily conserved catalytic mechanism shared among gram negative-bacteria.

Experimental Determination and Prediction of the Fitness Effects of Random Point Mutations in the Biosynthetic Enzyme HisA

The distribution of fitness effects of mutations is a factor of fundamental importance in evolutionary biology. We determined the distribution of fitness effects of 510 mutants that each carried between 1 and 10 mutations (synonymous and nonsynonymous) in the hisA gene, encoding an essential enzyme in the l-histidine biosynthesis pathway of Salmonella enterica. For the full set of mutants, the distribution was bimodal with many apparently neutral mutations and many lethal mutations. For a subset of 81 single, nonsynonymous mutants most mutations appeared neutral at high expression levels, whereas at low expression levels only a few mutations were neutral. Furthermore, we examined how the magnitude of the observed fitness effects was correlated to several measures of biophysical properties and phylogenetic conservation.We conclude that for HisA: (i) The effect of mutations can be masked by high expression levels, such that mutations that are deleterious to the function of the protein can still be neutral with regard to organism fitness if the protein is expressed at a sufficiently high level; (ii) the shape of the fitness distribution is dependent on the extent to which the protein is rate-limiting for growth; (iii) negative epistatic interactions, on an average, amplified the combined effect of nonsynonymous mutations; and (iv) no single sequence-based predictor could confidently predict the fitness effects of mutations in HisA, but a combination of multiple predictors could predict the effect with a SD of 0.04 resulting in 80% of the mutations predicted within 12% of their observed selection coefficients.

Structural and functional innovations in the real-time evolution of new (βα)8 barrel enzymes

New genes can arise by duplication and divergence, but there is a fundamental gap in our understanding of the relationship between these genes, the evolving proteins they encode, and the fitness of the organism. Here we used crystallography, NMR dynamics, kinetics, and mass spectrometry to explain the molecular innovations that arose during a previous real-time evolution experiment. In that experiment, the (βα)₈ barrel enzyme HisA was under selection for two functions (HisA and TrpF), resulting in duplication and divergence of the hisA gene to encode TrpF specialists, HisA specialists, and bifunctional generalists. We found that selection affects enzyme structure and dynamics, and thus substrate preference, simultaneously and sequentially. Bifunctionality is associated with two distinct sets of loop conformations, each essential for one function. We observed two mechanisms for functional specialization: structural stabilization of each loop conformation and substrate-specific adaptation of the active site. Intracellular enzyme performance, calculated as the product of catalytic efficiency and relative expression level, was not linearly related to fitness. Instead, we observed thresholds for each activity above which further improvements in catalytic efficiency had little if any effect on growth rate. Overall, we have shown how beneficial substitutions selected during real-time evolution can lead to manifold changes in enzyme function and bacterial fitness. This work emphasizes the speed at which adaptive evolution can yield enzymes with sufficiently high activities such that they no longer limit the growth of their host organism, and confirms the (βα)₈ barrel as an inherently evolvable protein scaffold.

The Amadori Rearrangement for Carbohydrate Conjugation: Scope and Limitations

The Amadori rearrangement was investigated for the synthesis of C-glycosyl-type neoglycoconjugates. Various amines including diamines, amino-functionalized glycosides, lysine derivatives, and peptides were conjugated with two different heptoses to generate non-natural C-glycosyl-type glycoconjugates of the d-gluco and d-manno series. With these studies, the scope and limitations of the Amadori rearrangement as a conjugation method have been exemplified with respect to the carbohydrate substrate, as well as the amino components.

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.

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.