Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography, docking, molecular dynamics, and free energy calculations

Abstract

Cytochrome P450 BM3 (CYP102A1) mutant M11 is able to metabolize a wide range of drugs and drug-like compounds. Among these, M11 was recently found to be able to catalyze formation of human metabolites of mefenamic acid and other nonsteroidal anti-inflammatory drugs (NSAIDs). Interestingly, single active-site mutations such as V87I were reported to invert regioselectivity in NSAID hydroxylation. In this work, we combine crystallography and molecular simulation to study the effect of single mutations on binding and regioselective metabolism of mefenamic acid by M11 mutants. The heme domain of the protein mutant M11 was expressed, purified, and crystallized, and its X-ray structure was used as template for modeling. A multistep approach was used that combines molecular docking, molecular dynamics (MD) simulation, and binding free-energy calculations to address protein flexibility. In this way, preferred binding modes that are consistent with oxidation at the experimentally observed sites of metabolism (SOMs) were identified. Whereas docking could not be used to retrospectively predict experimental trends in regioselectivity, we were able to rank binding modes in line with the preferred SOMs of mefenamic acid by M11 and its mutants by including protein flexibility and dynamics in free-energy computation. In addition, we could obtain structural insights into the change in regioselectivity of mefenamic acid hydroxylation due to single active-site mutations. Our findings confirm that use of MD and binding free-energy calculation is useful for studying biocatalysis in those cases in which enzyme binding is a critical event in determining the selective metabolism of a substrate.

Human Cytochrome P450 3A4 as a Biocatalyst: Effects of the Engineered Linker in Modulation of Coupling Efficiency in 3A4-BMR Chimeras

“…This has led to a wide range of protein engineering studies on this enzyme to widen its catalytic abilities (Tsotsou et al, 2012, 2013; Ryu et al, 2014; Di Nardo et al, 2015; Ren et al, 2015; Capoferri et al, 2016). …”

Role of Protein-Protein Interactions in Metabolism: Genetics, Structure, Function

“…This has led to a wide range of protein engineering studies on this enzyme to widen its catalytic abilities (Tsotsou et al, 2012(Tsotsou et al, , 2013Ryu et al, 2014;Di Nardo et al, 2015;Ren et al, 2015;Capoferri et al, 2016).…”

Advances in Computational Prediction of Regioselective and Isoform-Specific Drug Metabolism Catalyzed by CYP450s.

“…Differences in the regioselective metabolism of mefenamic acid by mutants of CYP450 BM3 were recently explained by Capoferri et al. using a combination of docking, MD and free energy calculations . Here docking alone could not explain regioselectivity patterns and inclusion of protein flexibility and accurate estimation of binding energies was found crucial.…”