Evolutionary fitness landscapes of certain antibiotic target enzymes have been comprehensively mapped showing strong high order epistasis between mutations, but understanding these effects at the biochemical and molecular levels remained open. Here, we carried out an extensive experimental and computational study to quantitatively understand the evolutionary dynamics of Escherichia coli dihydrofolate reductase (DHFR) enzyme in the presence of trimethoprim induced selection. Biochemical and structural characterization of resistance-conferring mutations targeting a total of ten residues spanning the substrate binding pocket of DHFR revealed distinct resistance mechanisms. Next, we experimentally measured biochemical parameters (Km, Ki, and kcat) for a mutant library carrying all possible combinations of six resistance-conferring DHFR mutations and quantified epistatic interactions between them. We found that the epistasis between DHFR mutations is high-order for catalytic power of DHFR (kcat and Km), but less prevalent for trimethoprim affinity (Ki). Taken together our data provide a concrete illustration of how epistatic coupling at the level of biochemical parameters can give rise to complex fitness landscapes, and suggest new strategies for developing mutant specific inhibitors.

High-order epistasis in catalytic power of dihydrofolate reductase gives rise to a rugged fitness landscape in the presence of trimethoprim selection