The Role of Gene Elongation in the Evolution of Histidine Biosynthetic Genes

Targeting imidazole-glycerol phosphate dehydratase in plants: novel approach for structural and functional studies, and inhibitor blueprinting

The histidine biosynthetic pathway (HBP) is targeted for herbicide design with preliminary success only regarding imidazole-glycerol phosphate dehydratase (IGPD, EC 4.2.1.19), or HISN5, as referred to in plants. HISN5 catalyzes the sixth step of the HBP, in which imidazole-glycerol phosphate (IGP) is dehydrated to imidazole-acetol phosphate. In this work, we present high-resolution cryoEM and crystal structures of Medicago truncatula HISN5 (MtHISN5) in complexes with an inactive IGP diastereoisomer and with various other ligands. MtHISN5 can serve as a new model for plant HISN5 structural studies, as it enables resolving protein-ligand interactions at high (2.2 Å) resolution using cryoEM. We identified ligand-binding hotspots and characterized the features of plant HISN5 enzymes in the context of the HISN5-targeted inhibitor design. Virtual screening performed against millions of small molecules not only revealed candidate molecules but also identified linkers for fragments that were experimentally confirmed to bind. Based on experimental and computational approaches, this study provides guidelines for designing symmetric HISN5 inhibitors that can reach two neighboring active sites. Finally, we conducted analyses of sequence similarity networks revealing that plant HISN5 enzymes derive from cyanobacteria. We also adopted a new approach to measure MtHISN5 enzymatic activity using isothermal titration calorimetry and enzymatically synthesized IGP.

Improving enzyme functional annotation by integrating in vitro and in silico approaches: The example of histidinol phosphate phosphatases

Advances in sequencing technologies have led to a rapid growth of public protein sequence databases, whereby the fraction of proteins with experimentally verified function continuously decreases. This problem is currently addressed by automated functional annotations with computational tools, which however lack the accuracy of experimental approaches and are susceptible to error propagation. Here, we present an approach that combines the efficiency of functional annotation by in silico methods with the rigor of enzyme characterization in vitro. First, a thorough experimental analysis of a representative enzyme of a group of homologues is performed which includes a focused alanine scan of the active site to determine a fingerprint of function-determining residues. In a second step, this fingerprint is used in combination with a sequence similarity network to identify putative isofunctional enzymes among the homologues. Using this approach in a proof-of-principle study, homologues of the histidinol phosphate phosphatase (HolPase) from Pseudomonas aeruginosa, many of which were annotated as phosphoserine phosphatases, were predicted to be HolPases. This functional annotation of the homologues was verified by in vitro testing of several representatives and an analysis of the occurrence of annotated HolPases in the corresponding phylogenetic groups. Moreover, the application of the same approach to the homologues of the HolPase from the archaeon Nitrosopumilus maritimus, which is not related to the HolPase from P. aeruginosa and was newly discovered in the course of this work, led to the annotation of the putative HolPase from various archaeal species.

Mimicking the last step of gene elongation: hints from the bacterial hisF gene

Gene elongation consists in an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. Many present-day proteins show internal repeats of amino acid sequences, generated by gene elongation events; however, gene elongation is still a poorly studied evolutionary molecular mechanism. The most documented case is that of the histidine biosynthetic genes hisA and hisF, which derive from the gene elongation of an ancestral gene half the size of the extant ones. The aim of this work was to experimentally simulate the possible last step of the gene elongation event occurred during hisF gene evolution under selective pressure conditions. Azospirillum brasilense hisF gene, carrying a single nucleotide mutation that generates a stop codon between the two halves of the gene, was used to transform the histidine-auxotrophic Escherichia coli strain FB182 (hisF892). The transformed strain was subjected to selective pressure (i.e., low concentration/absence of histidine in the growth medium) and the obtained mutants were characterized. The restoration of prototrophy was strongly dependent on the time of incubation and on the strength of the selective pressure. The mutations involved the introduced stop codon with a single base substitution and none of the mutants restored the wild-type codon. Possible correlations between the different mutations and i) E. coli codon usage, ii) three-dimensional structures of the mutated HisF proteins, and iii) growth ability of the mutants were investigated. On the contrary, when the experiment was repeated by mutating a more conserved codon, only a synonymous substitution was obtained. Thus, experiments performed in this study allowed to mimic a possible gene elongation event occurred during the evolution of hisF gene, evidencing the ability of bacterial cells to modify their genome in short times under selective conditions.

The Operon as a Conundrum of Gene Dynamics and Biochemical Constraints: What We Have Learned from Histidine Biosynthesis

Operons represent one of the leading strategies of gene organization in prokaryotes, having a crucial influence on the regulation of gene expression and on bacterial chromosome organization. However, there is no consensus yet on why, how, and when operons are formed and conserved, and many different theories have been proposed. Histidine biosynthesis is a highly studied metabolic pathway, and many of the models suggested to explain operons origin and evolution can be applied to the histidine pathway, making this route an attractive model for the study of operon evolution. Indeed, the organization of his genes in operons can be due to a progressive clustering of biosynthetic genes during evolution, coupled with a horizontal transfer of these gene clusters. The necessity of physical interactions among the His enzymes could also have had a role in favoring gene closeness, of particular importance in extreme environmental conditions. In addition, the presence in this pathway of paralogous genes, heterodimeric enzymes and complex regulatory networks also support other operon evolution hypotheses. It is possible that histidine biosynthesis, and in general all bacterial operons, may result from a mixture of several models, being shaped by different forces and mechanisms during evolution.

Application of Cloning-Free Genome Engineering to Escherichia coli

The propagation of foreign DNA in Escherichia coli is central to molecular biology. Recent advances have dramatically expanded the ability to engineer (bacterial) cells; however, most of these techniques remain time-consuming. The aim of the present work was to explore the possibility to use the cloning-free genome editing (CFGE) approach, proposed by Döhlemann and coworkers (2016), for E. coli genetics, and to deepen the knowledge about the homologous recombination mechanism. The E. coli auxotrophic mutant strains FB182 (hisF892) and FB181 (hisI903) were transformed with the circularized wild-type E. coli (i) hisF gene and hisF gene fragments of decreasing length, and (ii) hisIE gene, respectively. His⁺ clones were selected based on their ability to grow in the absence of histidine, and their hisF/hisIE gene sequences were characterized. CFGE method allowed the recombination of wild-type his genes (or fragments of them) within the mutated chromosomal copy, with a different recombination frequency based on the fragment length, and the generation of clones with a variable number of in tandem his genes copies. Data obtained pave the way to further evolutionary studies concerning the homologous recombination mechanism and the fate of in tandem duplicated genes.

Experimental and computational analysis of the ancestry of an evolutionary young enzyme from histidine biosynthesis

The conservation of fold and chemistry of the enzymes associated with histidine biosynthesis suggests that this pathway evolved prior to the diversification of Bacteria, Archaea, and Eukaryotes. The only exception is the histidinol phosphate phosphatase (HolPase). So far, non-homologous HolPases that possess distinct folds and belong to three different protein superfamilies have been identified in various phylogenetic clades. However, their evolution has remained unknown to date. Here, we analyzed the evolutionary history of the HolPase from γ-Proteobacteria (HisB-N). It has been argued that HisB-N and its closest homologue d-glycero-d-manno-heptose-1,7-bisphosphate 7-phosphatase (GmhB) have emerged from the same promiscuous ancestral phosphatase. GmhB variants catalyze the hydrolysis of the anomeric d-glycero-d-manno-heptose-1,7-bisphosphate (αHBP or βHBP) with a strong preference for one anomer (αGmhB or βGmhB). We found that HisB-N from Escherichia coli shows promiscuous activity for βHBP but not αHBP, while βGmhB from Crassaminicella sp. shows promiscuous activity for HolP. Accordingly, a combined phylogenetic tree of αGmhBs, βGmhBs, and HisB-N sequences revealed that HisB-Ns form a compact subcluster derived from βGmhBs. Ancestral sequence reconstruction and in vitro analysis revealed a promiscuous HolPase activity in the resurrected enzymes prior to functional divergence of the successors. The following increase in catalytic efficiency of the HolP turnover is reflected in the shape and electrostatics of the active site predicted by AlphaFold. An analysis of the phylogenetic tree led to a revised evolutionary model that proposes the horizontal gene transfer of a promiscuous βGmhB from δ- to γ-Proteobacteria where it evolved to the modern HisB-N.

The Azurin Coding Gene: Origin and Phylogenetic Distribution

Azurin is a bacterial-derived cupredoxin, which is mainly involved in electron transport reactions. Interest in azurin protein has risen in recent years due to its anticancer activity and its possible applications in anticancer therapies. Nevertheless, the attention of the scientific community only focused on the azurin protein found in Pseudomonas aeruginosa (Proteobacteria, Gammaproteobacteria). In this work, we performed the first comprehensive screening of all the bacterial genomes available in online repositories to assess azurin distribution in the three domains of life. The Azurin coding gene was not detected in the domains Archaea and Eucarya, whereas it was detected in phyla other than Proteobacteria, such as Bacteroidetes, Verrucomicrobia and Chloroflexi, and a phylogenetic analysis of the retrieved sequences was performed. Observed patchy distribution and phylogenetic data suggest that once it appeared in the bacterial domain, the azurin coding gene was lost in several bacterial phyla and/or anciently horizontally transferred between different phyla, even though a vertical inheritance appeared to be the major force driving the transmission of this gene. Interestingly, a shared conserved domain has been found among azurin members of all the investigated phyla. This domain is already known in P. aeruginosa as p28 domain and its importance for azurin anticancer activity has been widely explored. These findings may open a new and intriguing perspective in deciphering the azurin anticancer mechanisms and to develop new tools for treating cancer diseases.

The Histidine Biosynthetic Genes in the Superphylum Bacteroidota-Rhodothermota-Balneolota-Chlorobiota: Insights into the Evolution of Gene Structure and Organization

One of the most studied metabolic routes is the biosynthesis of histidine, especially in enterobacteria where a single compact operon composed of eight adjacent genes encodes the complete set of biosynthetic enzymes. It is still not clear how his genes were organized in the genome of the last universal common ancestor community. The aim of this work was to analyze the structure, organization, phylogenetic distribution, and degree of horizontal gene transfer (HGT) of his genes in the Bacteroidota-Rhodothermota-Balneolota-Chlorobiota superphylum, a group of phylogenetically close bacteria with different surviving strategies. The analysis of the large variety of his gene structures and organizations revealed different scenarios with genes organized in more or less compact—heterogeneous or homogeneous—operons, in suboperons, or in regulons. The organization of his genes in the extant members of the superphylum suggests that in the common ancestor of this group, genes were scattered throughout the chromosome and that different forces have driven the assembly of his genes in compact operons. Gene fusion events and/or paralog formation, HGT of single genes or entire operons between strains of the same or different taxonomic groups, and other molecular rearrangements shaped the his gene structure in this superphylum.

Structural and mechanistic insights into the bifunctional HISN2 enzyme catalyzing the second and third steps of histidine biosynthesis in plants

The second and third steps of the histidine biosynthetic pathway (HBP) in plants are catalyzed by a bifunctional enzyme–HISN2. The enzyme consists of two distinct domains, active respectively as a phosphoribosyl-AMP cyclohydrolase (PRA-CH) and phosphoribosyl-ATP pyrophosphatase (PRA-PH). The domains are analogous to single-domain enzymes encoded by bacterial hisI and hisE genes, respectively. The calculated sequence similarity networks between HISN2 analogs from prokaryotes and eukaryotes suggest that the plant enzymes are closest relatives of those in the class of Deltaproteobacteria. In this work, we obtained crystal structures of HISN2 enzyme from Medicago truncatula (MtHISN2) and described its architecture and interactions with AMP. The AMP molecule bound to the PRA-PH domain shows positioning of the N1-phosphoribosyl relevant to catalysis. AMP bound to the PRA-CH domain mimics a part of the substrate, giving insights into the reaction mechanism. The latter interaction also arises as a possible second-tier regulatory mechanism of the HBP flux, as indicated by inhibition assays and isothermal titration calorimetry.

Exploring the role of the histidine biosynthetic hisF gene in cellular metabolism and in the evolution of (ancestral) genes: from LUCA to the extant (micro)organisms

Histidine biosynthesis is an ancestral pathway that was assembled before the appearance of the Last Universal Common Ancestor; afterwards, it remained unaltered in all the extant histidine-synthesizing (micro)organisms. It is a metabolic cross-road interconnecting histidine biosynthesis to nitrogen metabolism and the de novo synthesis of purines. This interconnection is due to the reaction catalyzed by the products of hisH and hisF genes. The latter gene is an excellent model to study which trajectories have been followed by primordial cells to build the first metabolic routes, since its evolution is the result of different molecular rearrangement events, i.e. gene duplication, gene fusion, gene elongation, and domain shuffling. Additionally, this review summarizes data concerning the involvement of hisF and its product in other different cellular processes, revealing that HisF very likely plays a role also in cell division control and involvement in virulence and nodule development in different bacteria. From the metabolic viewpoint, these results suggest that HisF plays a central role in cellular metabolism, highlighting the interconnections of different metabolic pathways.