Molecular evolution is moving from statistical descriptions of adaptive molecular changes toward predicting the fitness effects of mutations. Here, we characterize the fitness landscape of the six amino acids controlling coenzyme use in isopropylmalate dehydrogenase (IMDH). Although all natural IMDHs use nicotinamide adenine dinucleotide (NAD) as a coenzyme, they can be engineered to use nicotinamide adenine dinucleotide phosphate (NADP) instead. Intermediates between these two phenotypic extremes show that each amino acid contributes additively to enzyme function, with epistatic contributions confined to fitness. The genotype-phenotype-fitness map shows that NAD use is a global optimum.

The functional effects of most amino acid replacements accumulated during molecular evolution are unknown, because most are not observed naturally and the possible combinations are too numerous. We created 168 single mutations in wild-type Escherichia coli isopropymalate dehydrogenase (IMDH) that match the differences found in wild-type Pseudomonas aeruginosa IMDH. 104 mutant enzymes performed similarly to E. coli wild-type IMDH, one was functionally enhanced, and 63 were functionally compromised. The transition from E. coli IMDH, or an ancestral form, to the functional wild-type P. aeruginosa IMDH requires extensive epistasis to ameliorate the combined effects of the deleterious mutations. This result stands in marked contrast with a basic assumption of molecular phylogenetics, that sites in sequences evolve independently of each other. Residues that affect function are scattered haphazardly throughout the IMDH structure. We screened for compensatory mutations at three sites, all of which lie near the active site and all of which are among the least active mutants. No compensatory mutations were found at two sites indicating that a single site may engage in compound epistatic interactions. One complete and three partial compensatory mutations of the third site are remote and lie in a different domain. This demonstrates that epistatic interactions can occur between distant (>20Å) sites. Phylogenetic analysis shows that incompatible mutations were fixed in different lineages.

Catalytic promiscuity is a useful, but accidental, enzyme property, so finding catalytically promiscuous enzymes in nature is inefficient. Some ancestral enzymes were branch points in the evolution of new enzymes and are hypothesized to have been promiscuous. To test the hypothesis that ancestral enzymes were more promiscuous than their modern descendants, we reconstructed ancestral enzymes at four branch points in the divergence hydroxynitrile lyases (HNL’s) from esterases ~100 million years ago. Both enzyme types are α/β-hydrolase-fold enzymes and have the same catalytic triad, but differ in reaction type and mechanism. Esterases catalyze hydrolysis via an acyl enzyme intermediate, while lyases catalyze an elimination without an intermediate. Screening ancestral enzymes and their modern descendants with six esterase substrates and six lyase substrates found higher catalytic promiscuity among the ancestral enzymes (P <0.01). Ancestral esterases were more likely to catalyze a lyase reaction than modern esterases and the ancestral HNL was more likely to catalyze ester hydrolysis than modern HNL’s. One ancestral enzyme (HNL1) along the path from esterase to hydroxynitrile lyases was especially promiscuous and catalyzed both hydrolysis and lyase reactions with many substrates. A broader screen tested mechanistically related reactions that were not selected for by evolution: decarboxylation, Michael addition, γ-lactam hydrolysis and 1,5-diketone hydrolysis. The ancestral enzymes were more promiscuous than their modern descendants (P = 0.04). Thus, these reconstructed ancestral enzyme are catalytically promiscuous, but HNL1 is especially so.