Free energy calculations elucidate substrate binding, gating mechanism, and tolerance-promoting mutations in herbicide target 4-hydroxyphenylpyruvate dioxygenase.

4-Hydroxyphenylpyruvate dioxygenase (HPPD) catalyzes the second reaction in the tyrosine catabolism and is linked to the production of cofactors plastoquinone and tocopherol in plants. This important biological role has put HPPD in the focus of current herbicide design efforts including the development of herbicide-tolerant mutants. However, the molecular mechanisms of substrate binding and herbicide tolerance have yet to be elucidated. In this work, we performed molecular dynamics simulations and free energy calculations to characterize active site gating by the C-terminal helix H11 in HPPD. We compared gating equilibria in Arabidopsis thaliana (At) and Zea mays (Zm) wild-type proteins retrieving the experimentally observed preferred orientations from the simulations. We investigated the influence of substrate and product binding on the open-closed transition and discovered a ligand-mediated conformational switch in H11 that mediates rapid substrate access followed by active site closing and efficient product release through H11 opening. We further studied H11 gating in At mutant HPPD, and found large differences with correlation to experimentally measured herbicide tolerance. The computational findings were then used to design a new At mutant HPPD protein that showed increased tolerance to six commercially available HPPD inhibitors in biochemical in vitro experiments. Our results underline the importance of protein flexibility and conformational transitions in substrate recognition and enzyme inhibition by herbicides.