A Modified Arrhenius Approach to Thermodynamically Study Regioselectivity in Cytochrome P450-Catalyzed Substrate Conversion.

The regio- (and stereo-)selectivity and specific activity of cytochrome P450s are determined by the accessibility of potential sites of metabolism (SOMs) of the bound substrate relative to the heme, and the activation barrier of the regioselective oxidation reaction(s). The accessibility of potential SOMs depends on the relative binding free energy (ΔΔGbind ) of the catalytically active substrate-binding poses, and the probability of the substrate to adopt a transition-state geometry. An established experimental method to measure activation energies of enzymatic reactions is the analysis of reaction rate constants at different temperatures and the construction of Arrhenius plots. This is a challenge for multistep P450-catalyzed processes that involve redox partners. We introduce a modified Arrhenius approach to overcome the limitations in studying P450 selectivity, which can be applied in multiproduct enzyme catalysis. Our approach gives combined information on relative activation energies, ΔΔGbind values, and collision entropies, yielding direct insight into the basis of selectivity in substrate conversion.

Structural analysis of Cytochrome P450 BM3 mutant M11 in complex with dithiothreitol.

The bacterial Cytochrome P450 (CYP) BM3 (CYP102A1) is one of the most active CYP isoforms. BM3 mutants can serve as a model for human drug-metabolizing CYPs and/or as biocatalyst for selective formation of drug metabolites. Hence, molecular and computational biologists have in the last two decades shown strong interest in the discovery and design of novel BM3 variants with optimized activity and selectivity for substrate conversion. This led e.g. to the discovery of mutant M11 that is able to metabolize a variety of drugs and drug-like compounds with relatively high activity. In order to further improve our understanding of CYP binding and reactions, we performed a co-crystallization study of mutant M11 and report here the three-dimensional structure M11 in complex with dithiothreitol (DTT) at a resolution of 2.16 Å. The structure shows that DTT can coordinate to the Fe atom in the heme group. UV/Vis spectroscopy and molecular dynamics simulation studies underline this finding and as first structure of the CYP BM3 mutant M11 in complex with a ligand, it offers a basis for structure-based design of novel mutants.

Linear Interaction Energy Based Prediction of Cytochrome P450 1A2 Binding Affinities with Reliability Estimation.

Prediction of human Cytochrome P450 (CYP) binding affinities of small ligands, i.e., substrates and inhibitors, represents an important task for predicting drug-drug interactions. A quantitative assessment of the ligand binding affinity towards different CYPs can provide an estimate of inhibitory activity or an indication of isoforms prone to interact with the substrate of inhibitors. However, the accuracy of global quantitative models for CYP substrate binding or inhibition based on traditional molecular descriptors can be limited, because of the lack of information on the structure and flexibility of the catalytic site of CYPs. Here we describe the application of a method that combines protein-ligand docking, Molecular Dynamics (MD) simulations and Linear Interaction Energy (LIE) theory, to allow for quantitative CYP affinity prediction. Using this combined approach, a LIE model for human CYP 1A2 was developed and evaluated, based on a structurally diverse dataset for which the estimated experimental uncertainty was 3.3 kJ mol-1. For the computed CYP 1A2 binding affinities, the model showed a root mean square error (RMSE) of 4.1 kJ mol-1 and a standard error in prediction (SDEP) in cross-validation of 4.3 kJ mol-1. A novel approach that includes information on both structural ligand description and protein-ligand interaction was developed for estimating the reliability of predictions, and was able to identify compounds from an external test set with a SDEP for the predicted affinities of 4.6 kJ mol-1 (corresponding to 0.8 pKi units).

Application of engineered cytochrome P450 mutants as biocatalysts for the synthesis of benzylic and aromatic metabolites of fenamic acid NSAIDs.

Cytochrome P450 BM3 mutants are promising biocatalysts for the production of drug metabolites. In the present study, the ability of cytochrome P450 BM3 mutants to produce oxidative metabolites of structurally related NSAIDs meclofenamic acid, mefenamic acid and tolfenamic acid was investigated. A library of engineered P450 BM3 mutants was screened with meclofenamic acid (1) to identify catalytically active and selective mutants. Three mono-hydroxylated metabolites were identified for 1. The hydroxylated products were confirmed by NMR analysis to be 3'-OH-methyl-meclofenamic acid (1a), 5-OH-meclofenamic acid (1b) and 4'-OH-meclofenamic acid (1c) which are human relevant metabolites. P450 BM3 variants containing V87I and V87F mutation showed high selectivity for benzylic and aromatic hydroxylation of 1 respectively. The applicability of these mutants to selectively hydroxylate structurally similar drugs such as mefenamic acid (2) and tolfenamic acid (3) was also investigated. The tested variants showed high total turnover numbers in the order of 4000-6000 and can be used as biocatalysts for preparative scale synthesis. Both 1 and 2 could undergo benzylic and aromatic hydroxylation by the P450 BM3 mutants, whereas 3 was hydroxylated only on aromatic rings. The P450 BM3 variant M11 V87F hydroxylated the aromatic ring at 4' position of all three drugs tested with high regioselectivity. Reference metabolites produced by P450 BM3 mutants allowed the characterisation of activity and regioselectivity of metabolism of all three NSAIDs by thirteen recombinant human P450s. In conclusion, engineered P450 BM3 mutants that are capable of benzylic or aromatic hydroxylation of fenamic acid containing NSAIDs, with high selectivity and turnover numbers have been identified. This shows their potential use as a greener alternative for the generation of drug metabolites.

Identification of overlapping but differential binding sites for the high-affinity CXCR3 antagonists NBI-74330 and VUF11211.

CXC chemokine receptor CXCR3 and/or its main three ligands CXCL9, CXCL10, and CXCL11 are highly upregulated in a variety of diseases. As such, considerable efforts have been made to develop small-molecule receptor CXCR3 antagonists, yielding distinct chemical classes of antagonists blocking binding and/or function of CXCR3 chemokines. Although it is suggested that these compounds bind in an allosteric fashion, thus far no evidence has been provided regarding the molecular details of their interaction with CXCR3. Using site-directed mutagenesis complemented with in silico homology modeling, we report the binding modes of two high-affinity CXCR3 antagonists of distinct chemotypes: VUF11211 [(S)-5-chloro-6-(4-(1-(4-chlorobenzyl)piperidin-4-yl)-3-ethylpiperazin-1-yl)-N-ethylnicotinamide] (piperazinyl-piperidine) with a rigid elongated structure containing two basic groups and NBI-74330 [(R)-N-(1-(3-(4-ethoxyphenyl)-4-oxo-3,4-dihydropyrido[2,3-d]pyrimidin-2-yl)ethyl)-2-(4-fluoro-3-(trifluoromethyl)phenyl)-N-(pyridin-3-ylmethyl)acetamide] (8-azaquinazolinone) without any basic group. Here we show that NBI-74330 is anchored in the transmembrane minor pocket lined by helices 2 (W2.60, D2.63), 3 (F3.32), and 7 (S7.39, Y7.43), whereas VUF11211 extends from the minor pocket into the major pocket of the transmembrane domains, located between residues in helices 1 (Y1.39), 2 (W2.60), 3 (F3.32), 4 (D4.60), 6 (Y6.51), and 7 (S7.39, Y7.43). Mutation of these residues did not affect CXCL11 binding significantly, confirming the allosteric nature of the interaction of these small molecules with CXCR3. Moreover, the model derived from our in silico-guided studies fits well with the already published structure-activity relationship data on these ligands. Altogether, in this study, we show overlapping, yet different binding sites for two high-affinity CXCR3 antagonists, which offer new opportunities for the structure-based design of allosteric modulators for CXCR3.