ABSTRACT Large enzyme families catalyze metabolic diversification by virtue of their ability to use diverse chemical scaffolds. How enzyme families attain such functional diversity is not clear. Here, we addressed this question using BAHD acyltransferases as a model, and identified the routes by which duplication, promiscuity and sequence changes influenced BAHD diversification. This fast-evolving family expanded drastically during land plant evolution from 1-5 copies in algae to ∼100 copies in diploid angiosperm genomes. In vitro characterization of fourteen BAHDs against a substrate panel and compilation of &gt;160 published activities revealed the wide prevalence of promiscuity among BAHDs. Using phylogenetic analysis, we predicted the substrate classes that the ancestral enzymes were likely capable of using prior to land plant origins. While the anthocyanin acylation activity was fixed in BAHDs later near the origin of angiosperms, in vitro testing of BAHDs from non-seed plant lineages suggested that the ability to acylate anthocyanins likely existed promiscuously millions of years prior to its fixation. Motif enrichment analysis in anthocyanin-acylating BAHDs identified two motifs fixed in the largest anthocyanin acylating clade. Molecular dynamic simulations and enzyme kinetics revealed the important role of an active site tryptophan, whose bulkiness, hydrophobicity and aromaticity are critical for anthocyanin acylation. Our results thus describe the molecular processes in robust, evolvable enzymes that drive emergence of functional diversity in enzyme families. One sentence summary Using a combination of phylogenetics, biochemistry and protein structure analysis, we investigated how the BAHD acyltransferase family evolved to use a structurally diverse array of substrates.

Multiple routes of functional diversification of the plant BAHD acyltransferase family revealed by comparative biochemical and genomic analyses