Cytochrome P450 monooxygenases play a crucial role in the biosynthesis of many natural products and in the human metabolism of numerous pharmaceuticals. This has inspired synthetic organic and medicinal chemists to exploit them as catalysts in regio- and stereoselective CH-activating oxidation of structurally simple and complex organic compounds such as steroids. However, levels of regio- and stereoselectivity as well as activity are not routinely high enough for real applications. Protein engineering using rational design or directed evolution has helped in many respects, but simultaneous engineering of multiple catalytic traits such as activity, regioselectivity and stereoselectivity, while overcoming tradeoffs and diminishing returns, remains a challenge. Here we show that the exploitation of information derived from mutability landscapes and molecular dynamics simulations for rationally designing iterative saturation mutagenesis constitutes a viable directed evolution strategy. This combined approach is illustrated by the evolution of P450BM3 mutants which enable nearly perfect regio- and diastereoselective hydroxylation of five different steroids specifically at the C16-position with unusually high activity, while avoiding activity-selectivity tradeoffs as well as keeping the screening effort relatively low. The C16 alcohols are of practical interest as components of biologically active glucocorticoids.

Steroidal C7β alcohols and their respective esters have shown significant promise as neuroprotective and anti‐inflammatory agents to treat chronic neuronal damage like stroke, brain trauma, and cerebral ischemia. Since C7 is spatially far away from any functional groups that could direct C−H activation, these transformations are not readily accessible using modern synthetic organic techniques. Reported here are P450‐BM3 mutants that catalyze the oxidative hydroxylation of six different steroids with pronounced C7 regioselectivities and β stereoselectivities, as well as high activities. These challenging transformations were achieved by a focused mutagenesis strategy and application of a novel technology for protein library construction based on DNA assembly and USER (Uracil‐Specific Excision Reagent) cloning. Upscaling reactions enabled the purification of the respective steroidal alcohols in moderate to excellent yields. The high‐resolution X‐ray structure and molecular dynamics simulations of the best mutant unveil the origin of regio‐ and stereoselectivity.

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Steroids are the most widely marketed products by the pharmaceutical industry after antibiotics. Steroid hydroxylation is one of the most important functionalizations because their derivatives enable a higher biological activity compared to their less polar non-hydroxylated analogs. Bacterial cytochrome P450s constitute promising biocatalysts for steroid hydroxylation due to their high expression level in common workhorses like Escherichia coli. However, they often suffer from wrong or insufficient regio- and/or stereoselectivity, low activity, narrow substrate range as well as insufficient thermostability, which hampers their industrial application. Fortunately, these problems can be generally solved by protein engineering based on directed evolution and rational design. In this work, an overview of recent developments on the engineering of bacterial cytochrome P450s for steroid hydroxylation is presented.