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

Distal Mutations in the β-Clamp of DNA Polymerase III* Disrupt DNA Orientation and Affect Exonuclease Activity

DNA polymerases are responsible for the replication and repair of DNA found in all DNA-based organisms. DNA Polymerase III is the main replicative polymerase of E. coli and is composed of over 10 proteins. A subset of these proteins (Pol III*) includes the polymerase (α), exonuclease (ϵ), clamp (β), and accessory protein (θ). Mutations of residues in, or around the active site of the catalytic subunits (α and ϵ), can have a significant impact on catalysis. However, the effects of distal mutations in noncatalytic subunits on the activity of catalytic subunits are less well-characterized. Here, we investigate the effects of two Pol III* variants, β-L82E/L82'E and β-L82D/L82'D, on the proofreading reaction catalyzed by ϵ. MD simulations reveal major changes in the dynamics of Pol III*, which extend throughout the complex. These changes are mostly induced by a shift in the position of the DNA substrate inside the β-clamp, although no major structural changes are observed in the protein complex. Quantum mechanics/molecular mechanics (QM/MM) calculations indicate that the β-L82D/L82'D variant has reduced catalytic proficiency due to highly endoergic reaction energies resulting from structural changes in the active site and differences in the electric field at the active site arising from the protein and substrate. Conversely, the β-L82E/L82'E variant is predicted to maintain proofreading activity, exhibiting a similar reaction barrier for nucleotide excision compared with the WT system. However, significant differences in the reaction mechanism are obtained due to the changes induced by the mutations on the β-clamp.

The sphingolipid anteome: implications for evolution of the sphingolipid metabolic pathway

Modern cell membranes contain a bewildering complexity of lipids, among them sphingolipids (SLs). Advances in mass spectrometry have led to the realization that the number and combinatorial complexity of lipids, including SLs, is much greater than previously appreciated. SLs are generated de novo by four enzymes, namely serine palmitoyltransferase, 3-ketodihydrosphingosine reductase, ceramide synthase and dihydroceramide Δ4-desaturase 1. Some of these enzymes depend on the availability of specific substrates and cofactors, which are themselves supplied by other complex metabolic pathways. The evolution of these four enzymes is poorly understood and likely depends on the co-evolution of the metabolic pathways that supply the other essential reaction components. Here, we introduce the concept of the 'anteome', from the Latin ante ('before') to describe the network of metabolic ('omic') pathways that must have converged in order for these pathways to co-evolve and permit SL synthesis. We also suggest that current origin of life and evolutionary models lack appropriate experimental support to explain the appearance of this complex metabolic pathway and its anteome.

Functional divergence of oligoadenylate synthetase 1 (OAS1) proteins in Tetrapods

OASs play critical roles in immune response against virus infection by polymerizing ATP into 2-5As, which initiate the classical OAS/RNase L pathway and induce degradation of viral RNA. OAS members are functionally diverged in four known innate immune pathways (OAS/RNase L, OASL/IRF7, OASL/RIG-I, and OASL/cGAS), but how they functionally diverged is unclear. Here, we focus on evolutionary patterns and explore the link between evolutionary processes and functional divergence of Tetrapod OAS1. We show that Palaeognathae and Primate OAS1 genes are conserved in genomic and protein structures but differ in function. The former (i.e., ostrich) efficiently synthesized long 2-5A and activated RNase L, while the latter (i.e., human) synthesized short 2-5A and did not activate RNase L. We predicted and verified that two in-frame indels and one positively selected site in the active site pocket contributed to the functional divergence of Palaeognathae and Primate OAS1. Moreover, we discovered and validated that an in-frame indel in the C-terminus of Palaeognathae OAS1 affected the binding affinity of dsRNA and enzymatic activity, and contributed to the functional divergence of Palaeognathae OAS1 proteins. Our findings unravel the molecular mechanism for functional divergence and give insights into the emergence of novel functions in Tetrapod OAS1.

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.

Horizontal-Acquisition of a Promiscuous Peptidoglycan-Recycling Enzyme Enables Aphids To Influence Symbiont Cell Wall Metabolism

During evolution, enzymes can undergo shifts in preferred substrates or in catalytic activities. An intriguing question is how enzyme function changes following horizontal gene transfer, especially for bacterial genes that have moved to animal genomes. Some insects have acquired genes that encode enzymes for the biosynthesis of bacterial cell wall components and that appear to function to support or control their obligate endosymbiotic bacteria. In aphids, the bacterial endosymbiont Buchnera aphidicola provides essential amino acids for aphid hosts but lacks most genes for remodeling of the bacterial cell wall. The aphid genome has acquired seven genes with putative functions in cell wall metabolism that are primarily expressed in the aphid cells harboring Buchnera. In analyses of aphid homogenates, we detected peptidoglycan (PGN) muropeptides indicative of the reactions of PGN hydrolases encoded by horizontally acquired aphid genes but not by Buchnera genes. We produced one such host enzyme, ApLdcA, and characterized its activity with both cell wall derived and synthetic PGN. Both ApLdcA and the homologous enzyme in Escherichia coli, which functions as an l,d-carboxypeptidase in the cytoplasmic PGN recycling pathway, exhibit turnover of PGN substrates containing stem pentapeptides and cross-linkages via l,d-endopeptidase activity, consistent with a potential role in cell wall remodeling. Our results suggest that ApLdcA derives its functions from the promiscuous activities of an ancestral LdcA enzyme, whose acquisition by the aphid genome may have enabled hosts to influence Buchnera cell wall metabolism as a means to control symbiont growth and division. IMPORTANCE Most enzymes are capable of performing biologically irrelevant side reactions. During evolution, promiscuous enzyme activities may acquire new biological roles, especially after horizontal gene transfer to new organisms. Pea aphids harbor obligate bacterial symbionts called Buchnera and encode horizontally acquired bacterial genes with putative roles in cell wall metabolism. Though Buchnera lacks cell wall endopeptidase genes, we found evidence of endopeptidase activity among peptidoglycan muropeptides purified from aphids. We characterized a multifunctional, aphid-encoded enzyme, ApLdcA, which displays l,d-endopeptidase activities considered promiscuous for the Escherichia coli homolog, for which these activities do not contribute to its native role in peptidoglycan recycling. These results exemplify the roles of enzyme promiscuity and horizontal gene transfer in enzyme evolution and demonstrate how aphids influence symbiont cell wall metabolism.

Rhodococcus comparative genomics reveals a phylogenomic-dependent non-ribosomal peptide synthetase distribution: insights into biosynthetic gene cluster connection to an orphan metabolite

Natural products (NPs) are synthesized by biosynthetic gene clusters (BGCs), whose genes are involved in producing one or a family of chemically related metabolites. Advances in comparative genomics have been favourable for exploiting huge amounts of data and discovering previously unknown BGCs. Nonetheless, studying distribution patterns of novel BGCs and elucidating the biosynthesis of orphan metabolites remains a challenge. To fill this knowledge gap, our study developed a pipeline for high-quality comparative genomics for the actinomycete genus Rhodococcus , which is metabolically versatile, yet understudied in terms of NPs, leading to a total of 110 genomes, 1891 BGCs and 717 non-ribosomal peptide synthetases (NRPSs). Phylogenomic inferences showed four major clades retrieved from strains of several ecological habitats. BiG-SCAPE sequence similarity BGC networking revealed 44 unidentified gene cluster families (GCFs) for NRPS, which presented a phylogenomic-dependent evolution pattern, supporting the hypothesis of vertical gene transfer. As a proof of concept, we analysed in-depth one of our marine strains, Rhodococcus sp. H-CA8f, which revealed a unique BGC distribution within its phylogenomic clade, involved in producing a chloramphenicol-related compound. While this BGC is part of the most abundant and widely distributed NRPS GCF, corason analysis unveiled major differences regarding its genetic context, co-occurrence patterns and modularity. This BGC is composed of three sections, two well-conserved right/left arms flanking a very variable middle section, composed of nrps genes. The presence of two non-canonical domains in H-CA8f’s BGC may contribute to adding chemical diversity to this family of NPs. Liquid chromatography-high resolution MS and dereplication efforts retrieved a set of related orphan metabolites, the corynecins, which to our knowledge are reported here for the first time in Rhodococcus . Overall, our data provide insights to connect BGC uniqueness with orphan metabolites, by revealing key comparative genomic features supported by models of BGC distribution along phylogeny.

In silico functional and evolutionary analyses of rubber oxygenases (RoxA and RoxB)

The study presents an in silico identification of poly (cis-1,4-isoprene) cleaving enzymes, viz., RoxA and RoxB in bacteria, followed by their functional and evolutionary exploration using comparative genomics. The orthologs of these proteins were found to be restricted to Gram-negative beta-, gamma-, and delta-proteobacteria. Toward the evolutionary propagation, the RoxA and RoxB genes were predicted to have evolved via a common interclass route of horizontal gene transfer in the phylum Proteobacteria (delta → gamma → beta). Besides, recombination, mutation, and gene conversion were also detected in both the genes leading to their diversification. Further, the differential selective pressure is predicted to be operating on entire RoxA and RoxB genes such that the former is diversifying further, whereas the latter is evolving to reduce its genetic diversity. However, the structurally and functionally important sites/residues of these genes were found to be preventing changes implying their evolutionary conservation. Further, the phylogenetic analysis demonstrated a sharp split between the RoxA and RoxB orthologs and indicated the emergence of their variant as another type of putative rubber oxygenase (RoxC) in the class Gammaproteobacteria. A detailed in silico analysis of the signature motifs and residues of Rox sequences exhibited important differences as well as similarities among the RoxA, RoxB, and putative RoxC sequences. Although RoxC appears to be a hybrid of RoxA and RoxB, the signature motifs and residues of RoxC are more similar to RoxB.

3D domain swapping in the TIM barrel of the α subunit of Streptococcus pneumoniae tryptophan synthase

Tryptophan synthase catalyzes the last two steps of tryptophan biosynthesis in plants, fungi and bacteria. It consists of two protein chains, designated α and β, encoded by trpA and trpB genes, that function as an αββα complex. Structural and functional features of tryptophan synthase have been extensively studied, explaining the roles of individual residues in the two active sites in catalysis and allosteric regulation. TrpA serves as a model for protein-folding studies. In 1969, Jackson and Yanofsky observed that the typically monomeric TrpA forms a small population of dimers. Dimerization was postulated to take place through an exchange of structural elements of the monomeric chains, a phenomenon later termed 3D domain swapping. The structural details of the TrpA dimer have remained unknown. Here, the crystal structure of the Streptococcus pneumoniae TrpA homodimer is reported, demonstrating 3D domain swapping in a TIM-barrel fold for the first time. The N-terminal domain comprising the H0-S1-H1-S2 elements is exchanged, while the hinge region corresponds to loop L2 linking strand S2 to helix H2'. The structural elements S2 and L2 carry the catalytic residues Glu52 and Asp63. As the S2 element is part of the swapped domain, the architecture of the catalytic apparatus in the dimer is recreated from two protein chains. The homodimer interface overlaps with the α-β interface of the tryptophan synthase αββα heterotetramer, suggesting that the 3D domain-swapped dimer cannot form a complex with the β subunit. In the crystal, the dimers assemble into a decamer comprising two pentameric rings.

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).

Evolutionary dynamics of natural product biosynthesis in bacteria

Covering: 2008 up to 2019The forces of biochemical adaptive evolution operate at the level of genes, manifesting in complex phenotypes and the global biodiversity of proteins and metabolites. While evolutionary histories have been deciphered for some other complex traits, the origins of natural product biosynthesis largely remain a mystery. This fundamental knowledge gap is surprising given the many decades of research probing the genetic, chemical, and biophysical mechanisms of bacterial natural product biosynthesis. Recently, evolutionary thinking has begun to permeate this otherwise mechanistically dominated field. Natural products are now sometimes referred to as 'specialized' rather than 'secondary' metabolites, reinforcing the importance of their biological and ecological functions. Here, we review known evolutionary mechanisms underlying the overwhelming chemical diversity of bacterial secondary metabolism, focusing on enzyme promiscuity and the evolution of enzymatic domains that enable metabolic traits. We discuss the mechanisms that drive the assembly of natural product biosynthetic gene clusters and propose formal definitions for 'specialized' and 'secondary' metabolism. We further explore how biosynthetic gene clusters evolve to synthesize related molecular species, and in turn how the biological and ecological roles that emerge from metabolic diversity are acted on by selection. Finally, we reconcile chemical, functional, and genetic data into an evolutionary model, the dynamic chemical matrix evolutionary hypothesis, in which the relationships between chemical distance, biomolecular activity, and relative fitness shape adaptive landscapes.

Mutations that improve efficiency of a weak-link enzyme are rare compared to adaptive mutations elsewhere in the genome

New enzymes often evolve by gene amplification and divergence. Previous experimental studies have followed the evolutionary trajectory of an amplified gene, but have not considered mutations elsewhere in the genome when fitness is limited by an evolving gene. We have evolved a strain of Escherichia coli in which a secondary promiscuous activity has been recruited to serve an essential function. The gene encoding the ‘weak-link’ enzyme amplified in all eight populations, but mutations improving the newly needed activity occurred in only one. Most adaptive mutations occurred elsewhere in the genome. Some mutations increase expression of the enzyme upstream of the weak-link enzyme, pushing material through the dysfunctional metabolic pathway. Others enhance production of a co-substrate for a downstream enzyme, thereby pulling material through the pathway. Most of these latter mutations are detrimental in wild-type E. coli, and thus would require reversion or compensation once a sufficient new activity has evolved.

Functional Profiling and Crystal Structures of Isothiocyanate Hydrolases Found in Gut-Associated and Plant-Pathogenic Bacteria

Isothiocyanates (ITCs) are produced by cruciferous plants to protect them against herbivores and infection by microbes. These compounds are of particular interest due to their antimicrobial and anticarcinogenic properties. The breakdown of ITCs in nature is catalyzed by isothiocyanate hydrolases (ITCases), a novel family within the metallo-β-lactamase (MBL)-fold superfamily of proteins. saxA genes that code for ITCases are particularly widespread in insect- and plant-associated bacteria. Enzymatic characterization of seven phylogenetically related but distinct ITCases revealed similar activities on six selected ITCs, suggesting that phylogenetic diversity does not determine the substrate specificity of ITCases. X-ray crystallography studies of two ITCases sharing 42% amino acid sequence identity revealed a highly conserved tertiary structure. Notable features of ITCases include a hydrophobic active site with two Zn²⁺ ions coordinating water/hydroxide and a flexible cap that is implicated in substrate recognition and covers the active site. This report reveals the function and structure of the previously uncharacterized family of isothiocyanate hydrolases within the otherwise relatively well-studied superfamily of metallo-β-lactamases.IMPORTANCE This study explores a newly discovered protein in the β-lactamase superfamily, namely, SaxA, or isothiocyanate hydrolase. Isothiocyanates are defensive compounds found in many cabbage-related crop plants and are currently being investigated for their antimicrobial and anticarcinogenic properties. We show that isothiocyanate hydrolases are responsible for the breakdown of several of these plant defensive chemicals in vitro and suggest their potential for mitigating the beneficial effects of isothiocyanates in crop protection and cancer prevention.

Chlamydia trachomatis dapF Encodes a Bifunctional Enzyme Capable of Both d-Glutamate Racemase and Diaminopimelate Epimerase Activities

Peptidoglycan is a sugar/amino acid polymer unique to bacteria and essential for division and cell shape maintenance. The d-amino acids that make up its cross-linked stem peptides are not abundant in nature and must be synthesized by bacteria de novo d-Glutamate is present at the second position of the pentapeptide stem and is strictly conserved in all bacterial species. In Gram-negative bacteria, d-glutamate is generated via the racemization of l-glutamate by glutamate racemase (MurI). Chlamydia trachomatis is the leading cause of infectious blindness and sexually transmitted bacterial infections worldwide. While its genome encodes a majority of the enzymes involved in peptidoglycan synthesis, no murI homologue has ever been annotated. Recent studies have revealed the presence of peptidoglycan in C. trachomatis and confirmed that its pentapeptide includes d-glutamate. In this study, we show that C. trachomatis synthesizes d-glutamate by utilizing a novel, bifunctional homologue of diaminopimelate epimerase (DapF). DapF catalyzes the final step in the synthesis of meso-diaminopimelate, another amino acid unique to peptidoglycan. Genetic complementation of an Escherichia coli murI mutant demonstrated that Chlamydia DapF can generate d-glutamate. Biochemical analysis showed robust activity, but unlike canonical glutamate racemases, activity was dependent on the cofactor pyridoxal phosphate. Genetic complementation, enzymatic characterization, and bioinformatic analyses indicate that chlamydial DapF shares characteristics with other promiscuous/primordial enzymes, presenting a potential mechanism for d-glutamate synthesis not only in Chlamydia but also numerous other genera within the Planctomycetes-Verrucomicrobiae-Chlamydiae superphylum that lack recognized glutamate racemases.IMPORTANCE Here we describe one of the last remaining "missing" steps in peptidoglycan synthesis in pathogenic Chlamydia species, the synthesis of d-glutamate. We have determined that the diaminopimelate epimerase (DapF) encoded by Chlamydia trachomatis is capable of carrying out both the epimerization of DAP and the pyridoxal phosphate-dependent racemization of glutamate. Enzyme promiscuity is thought to be the hallmark of early microbial life on this planet, and there is currently an active debate as to whether "moonlighting enzymes" represent primordial evolutionary relics or are a product of more recent reductionist evolutionary pressures. Given the large number of Chlamydia species (as well as members of the Planctomycetes-Verrucomicrobiae-Chlamydiae superphylum) that possess DapF but lack homologues of MurI, it is likely that DapF is a primordial isomerase that functions as both racemase and epimerase in these organisms, suggesting that specialized d-glutamate racemase enzymes never evolved in these microbes.