Molecular Dynamics Explorations of Active Site Structure in Designed and Evolved Enzymes

This Account describes the use of molecular dynamics (MD) simulations to reveal how mutations alter the structure and organization of enzyme active sites. As proposed by Pauling about 70 years ago and elaborated by many others since then, biocatalysis is efficient when functional groups in the active site of an enzyme are in optimal positions for transition state stabilization. Changes in mechanism and covalent interactions are often critical parts of enzyme catalysis. We describe our explorations of the dynamical preorganization of active sites using MD, studying the fluctuations between active and inactive conformations normally concealed to static crystallography. MD shows how the various arrangements of active site residues influence the free energy of the transition state and relates the populations of the catalytic conformational ensemble to the enzyme activity. This Account is organized around three case studies from our laboratory. We first describe the importance of dynamics in evaluating a series of computationally designed and experimentally evolved enzymes for the Kemp elimination, a popular subject in the enzyme design field. We find that the dynamics of the active site is influenced not only by the original sequence design and subsequent mutations but also by the nature of the ligand present in the active site. In the second example, we show how microsecond MD has been used to uncover the role of remote mutations in the active site dynamics and catalysis of a transesterase, LovD. This enzyme was evolved by Tang at UCLA and Codexis, Inc., and is a useful commercial catalyst for the production of the drug simvastatin. X-ray analysis of inactive and active mutants did not reveal differences in the active sites, but relatively long time scale MD in solution showed that the active site of the wild-type enzyme preorganizes only upon binding of the acyl carrier protein (ACP) that delivers the natural acyl group to the active site. In the absence of bound ACP, a noncatalytic arrangement of the catalytic triad is dominant. Unnatural truncated substrates are inactive because of the lack of protein-protein interactions provided by the ACP. Directed evolution is able to gradually restore the catalytic organization of the active site by motion of the protein backbone that alters the active site geometry. In the third case, we demonstrate the key role of MD in combination with crystallography to identify the origins of substrate-dependent stereoselectivities in a number of Codexis-engineered ketoreductases, one of which is used commercially for the production of the antibiotic sulopenem. Here, mutations alter the shape of the active site as well as the accessibility of water to different regions of it. Each of these examples reveals something different about how mutations can influence enzyme activity and shows that directed evolution, like natural evolution, can increase catalytic activity in a variety of remarkable and often subtle ways.