AN ANALYSIS OF THE REGIOSELECTIVITY OF AROMATIC HYDROXYLATION AND N-OXYGENATION BY CYTOCHROME P450 ENZYMES

Enzymatic approaches to site-selective oxidation of quinoline and derivatives

Enzyme-mediated oxidation has been a green and efficient strategy for preparation of derivative chemicals from quinoline and its structural analogues. Herein, we report the progress made to date in enzymatic methods to oxidation of the pyridine moieties of quinoline and its structural analogues 1,2,3,4-tetrahydroquinoline, isoquinoline and 1,2,3,4-tetrahydroisoquinoline, including whole cell- and isolated enzyme-based transformations. In addition, methods to tune the site selectivity of the course of enzymatic transformation are also addressed, in particular the protein engineering approaches.

Toluene Dioxygenase-Catalyzed cis-Dihydroxylation of Quinolines: A Molecular Docking Study and Chemoenzymatic Synthesis of Quinoline Arene Oxides

Molecular docking studies of quinoline and 2-chloroquinoline substrates at the active site of toluene dioxygenase (TDO), were conducted using Autodock Vina, to identify novel edge-to-face interactions and to rationalize the observed stereoselective cis-dihydroxylation of carbocyclic rings and formation of isolable cis-dihydrodiol metabolites. These in silico docking results of quinoline and pyridine substrates, with TDO, also provided support for the postulated cis-dihydroxylation of electron-deficient pyridyl rings, to give transient cis-dihydrodiol intermediates and the derived hydroxyquinolines. 2-Chloroquinoline cis-dihydrodiol metabolites were used as precursors in the chemoenzymatic synthesis of enantiopure arene oxide and arene dioxide derivatives of quinoline, in the context of its possible mammalian metabolism and carcinogenicity.

Decarboxylative C-H Alkylation of Heteroarene N-Oxides by Visible Light/Copper Catalysis

This paper reports a highly site-selective alkylation of heteroarene N-oxides using hypervalent iodine(III) carboxylates to serve as an alkylating agent in the presence of a cheap copper catalyst under visible light conditions. This mild method proceeds at room temperature in an air atmosphere and can withstand various heteroarene N-oxides as well as various primary, secondary, and tertiary alkyl carboxylic acids. It also provides a practical method for enabling the rapid conversion of commercially available raw materials into medically relevant "drug-like" molecules.

Site-Directed Mutagenesis at the Molybdenum Pterin Cofactor Site of the Human Aldehyde Oxidase: Interrogating the Kinetic Differences Between Human and Cynomolgus Monkey

The estimation of the drug clearance by aldehyde oxidase (AO) has been complicated because of this enzyme's atypical kinetics and species and substrate specificity. Since human AO (hAO) and cynomolgus monkey AO (mAO) have a 95.1% sequence identity, cynomolgus monkeys may be the best species for estimating AO clearance in humans. Here, O⁶-benzylguanine (O6BG) and dantrolene were used under anaerobic conditions, as oxidative and reductive substrates of AO, respectively, to compare and contrast the kinetics of these two species through numerical modeling. Whereas dantrolene reduction followed the same linear kinetics in both species, the oxidation rate of O6BG was also linear in mAO and did not follow the already established biphasic kinetics of hAO. In an attempt to determine why hAO and mAO are kinetically distinct, we have altered the hAO V811 and F885 amino acids at the oxidation site adjacent to the molybdenum pterin cofactor to the corresponding alanine and leucine in mAO, respectively. Although some shift to a more monkey-like kinetics was observed for the V811A mutant, five more mutations around the AO cofactors still need to be investigated for this purpose. In comparing the oxidative and reductive rates of metabolism under anaerobic conditions, we have come to the conclusion that despite having similar rates of reduction (4-fold difference), the oxidation rate in mAO is more than 50-fold slower than hAO. This finding implies that the presence of nonlinearity in AO kinetics is dependent upon the degree of imbalance between the rates of oxidation and reduction in this enzyme. SIGNIFICANCE STATEMENT: Although they have as much as 95.1% sequence identity, human and cynomolgus monkey aldehyde oxidase are kinetically distinct. Therefore, monkeys may not be good estimators of drug clearance in humans.

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.

Aldehyde oxidase and its role as a drug metabolizing enzyme

Aldehyde oxidase (AO) is a cytosolic enzyme that belongs to the family of structurally related molybdoflavoproteins like xanthine oxidase (XO). The enzyme is characterized by broad substrate specificity and marked species differences. It catalyzes the oxidation of aromatic and aliphatic aldehydes and various heteroaromatic rings as well as reduction of several functional groups. The references to AO and its role in metabolism date back to the 1950s, but the importance of this enzyme in the metabolism of drugs has emerged in the past fifteen years. Several reviews on the role of AO in drug metabolism have been published in the past decade indicative of the growing interest in the enzyme and its influence in drug metabolism. Here, we present a comprehensive monograph of AO as a drug metabolizing enzyme with emphasis on marketed drugs as well as other xenobiotics, as substrates and inhibitors. Although the number of drugs that are primarily metabolized by AO are few, the impact of AO on drug development has been extensive. We also discuss the effect of AO on the systemic exposure and clearance these clinical candidates. The review provides a comprehensive analysis of drug discovery compounds involving AO with the focus on developmental candidates that were reported in the past five years with regards to pharmacokinetics and toxicity. While there is only one known report of AO-mediated clinically relevant drug-drug interaction (DDI), a detailed description of inhibitors and inducers of AO known to date has been presented here and the potential risks associated with DDI. The increasing recognition of the importance of AO has led to significant progress in predicting the site of AO-mediated metabolism using computational methods. Additionally, marked species difference in expression of AO makes it is difficult to predict human clearance with high confidence. The progress made towards developing in vivo, in vitro and in silico approaches for predicting AO metabolism and estimating human clearance of compounds that are metabolized by AO have also been discussed.

Degradation of pristine and oxidized single wall carbon nanotubes by CYP3A4

Carbon nanotubes (CNTs) are a class of carbon based nanomaterials which have attracted substantial attention in recent years as they exhibit outstanding physical, mechanical and optical properties. In the last decade many studies have emerged of the underlying mechanisms behind CNT toxicity including malignant transformation, the formation of granulomas, inflammatory responses, oxidative stress, DNA damage and mutation. In the present investigation, we studied the biodegradation of single-walled carbon nanotubes (SWCNTs) by Cytochrome P450 enzymes (CYP3A4) through using Raman spectroscopy. CYP3A4 is known isozyme accountable for metabolizing various endogenous and exogenous xenobiotics. CYP3A4 is expressed dominantly in the liver and other organs including the lungs. Our results suggest that CYP3A4 has a higher affinity for p-SWNTs compared to c-SWNTs. HEK293 cellular viability was not compromised when incubated with SWNT. However, CYP3A4 transfected HEK293 cell line showed no digestion of c-SWNTs after incubation for 96 h. Cellular uptake of c-SWNTs was observed by electron microscopy and localization of c-SWNTs was confirmed in endosomal vesicles and in the cytoplasm. This is the first study CYP3A4 degrading both p-SWNTs and c-SWNTs in an in vitro setup. Interestingly, our results show that CYP3A4 is more proficient in degrading p-SWNTs than c-SWNTs. We also employed computational modeling and docking assessments to develop a further understanding of the molecular interaction mechanism.

Differences and comparisons of the properties and reactivities of iron(III)-hydroperoxo complexes with saturated coordination sphere

Heme and nonheme monoxygenases and dioxygenases catalyze important oxygen atom transfer reactions to substrates in the body. It is now well established that the cytochrome P450 enzymes react through the formation of a high-valent iron(IV)-oxo heme cation radical. Its precursor in the catalytic cycle, the iron(III)-hydroperoxo complex, was tested for catalytic activity and found to be a sluggish oxidant of hydroxylation, epoxidation and sulfoxidation reactions. In a recent twist of events, evidence has emerged of several nonheme iron(III)-hydroperoxo complexes that appear to react with substrates via oxygen atom transfer processes. Although it was not clear from these studies whether the iron(III)-hydroperoxo reacted directly with substrates or that an initial O-O bond cleavage preceded the reaction. Clearly, the catalytic activity of heme and nonheme iron(III)-hydroperoxo complexes is substantially different, but the origins of this are still poorly understood and warrant a detailed analysis. In this work, an extensive computational analysis of aromatic hydroxylation by biomimetic nonheme and heme iron systems is presented, starting from an iron(III)-hydroperoxo complex with pentadentate ligand system (L5(2)). Direct C-O bond formation by an iron(III)-hydroperoxo complex is investigated, as well as the initial heterolytic and homolytic bond cleavage of the hydroperoxo group. The calculations show that [(L5(2))Fe(III)(OOH)](2+) should be able to initiate an aromatic hydroxylation process, although a low-energy homolytic cleavage pathway is only slightly higher in energy. A detailed valence bond and thermochemical analysis rationalizes the differences in chemical reactivity of heme and nonheme iron(III)-hydroperoxo and show that the main reason for this particular nonheme complex to be reactive comes from the fact that they homolytically split the O-O bond, whereas a heterolytic O-O bond breaking in heme iron(III)-hydroperoxo is found.

Understanding a substrate's product regioselectivity in a family of enzymes: a case study of acetaminophen binding in cytochrome P450s

Product regioselectivity as influenced by molecular recognition is a key aspect of enzyme catalysis. We applied large-scale two-dimensional (2D) umbrella sampling (USP) simulations to characterize acetaminophen (APAP) binding in the active sites of the family of Cytochrome P450 (CYP) enzymes as a case study to show the different regioselectivity exhibited by a single substrate in comparative enzymes. Our results successfully explain the experimentally observed product regioselectivity for all five human CYPs included in this study, demonstrating that binding events play an important role in determining regioselectivity. In CYP2C9 and CYP3A4, weak interactions in an overall large active site cavity result in a fairly small binding free energy difference between APAP reactive binding states, consistent with experimental results that show little preference for resulting metabolites. In contrast, in CYP1A2 and CYP2E1, APAP is strongly restrained by a compact binding pocket, leading to a preferred binding conformation. The calculated binding equilibrium of APAP within the compact active site of CYP2A6 is able to predict the experimentally documented product ratios and is also applied to explain APAP regioselectivity in CYP1A2 and CYP2C9. APAP regioselectivity seems to be related to the selectivity for one binding conformation over another binding conformation as dictated by the size and shape of the active site. Additionally, unlike docking and molecular dynamics (MD), our free energy calculations successfully reproduced a unique APAP pose in CYP3A4 that had been reported experimentally, suggesting this approach is well suited to find the realistic binding pose and the lowest-energy starting structure for studying the chemical reaction step in the future.

Evaluation of in vitro metabolic systems for common drugs of abuse. 1. Cocaine

This study examined the efficacy of four common in vitro assay systems in producing metabolic profiles consistent with in vivo data for drugs of abuse. Cocaine (COC) was selected for this study because of its complex biotransformation pathways, diverse metabolic processes and because extensive Phase I and Phase II metabolomic examination of COC has not yet been reported by means of in vitro assay. COC metabolism was assessed with a series of common in vitro assay systems (human liver microsomes, cytosol and human liver S9 fraction and horseradish peroxidase) using liquid chromatography-tandem mass spectrometry with multiple reaction monitoring. Qualitative and quantitative differences in analyte production were noted among the various active Phase I and Phase II metabolic systems. Assay incubation time was found to be a determining factor in metabolic profile, specifically with primary versus secondary metabolite formation. Regioselective arene hydroxylation of COC was conclusively documented in human hepatic metabolic models, while peroxidase-based assay systems displayed less selectivity in oxidative aryl biotransformation. Results demonstrate the applicability of in vitro systems in studying COC metabolite production and the impact of assay selection and variation in method parameters on metabolite profiles for this important drug of abuse.

Construction of a CYP2E1-template system for prediction of the metabolism on both site and preference order

We have constructed an in silico system for the prediction of CYP2E1-mediated reaction using a two-dimensional template derived from substrate structures. Although CYP2E1 prefers small-size molecules for the substrates, the enzyme mediates oxidations of large-size molecules, such as benzo[a]pyrene. Overlays of these substrates, to assemble their sites of oxidation into a specific area, suggested a range of regions frequently occupied. The region, having a benzo[a]pyrene-like shape, was thus used as a CYP2E1 template. In this system, atoms in substrates, except for hydrogen atoms, were placed on corners of honeycomb structures of the template after having expanded the structures. Using published data for the metabolism on more than 80 substrates of CYP2E1, the core template was further refined to verify the adjacent area and to define the relative contribution of template positions for the catalysis. The positions on the template were classified into four different point (0-3) groups, depending on relative usage. In addition, we set independent points (-5 to 3) for specific positions to incorporate three-dimensional or functional information. Total scores from both position-occupancy and -function points were calculated for all the orientations of possible conformers of test substrates, and the scores were found to predict the relative abundance (i.e., order) as well as the regioselectivity of human CYP2E1 reactions with high fidelities.

Comparative study of the affinity and metabolism of type I and type II binding quinoline carboxamide analogues by cytochrome P450 3A4

Compounds that coordinate to the heme-iron of cytochrome P450 (CYP) enzymes are assumed to increase metabolic stability. However, recently we observed that the type II binding quinoline carboxamide (QCA) compounds were metabolically less stable. To test if the higher intrinsic clearance of type II binding compounds relative to type I binding compounds is general for other metabolic transformations, we synthesized a library of QCA compounds that could undergo N-dealkylation, O-dealkylation, benzylic hydroxylation, and aromatic hydroxylation. The results demonstrated that type II binding QCA analogues were metabolically less stable (2- to 12-fold) at subsaturating concentration compared to type I binding counterparts for all the transformations. When the rates of different metabolic transformations between type I and type II binding compounds were compared, they were found to be in the order of N-demethylation > benzylic hydroxylation> O-demethylation > aromatic hydroxylation. Finally, for the QCA analogues with aza-heteroaromatic rings, we did not detect metabolism in aza-aromatic rings (pyridine, pyrazine, pyrimidine), indicating that electronegativity of the nitrogen can change regioselectivity in CYP metabolism.

The effects of nitrogen-heme-iron coordination on substrate affinities for cytochrome P450 2E1

A descriptor based computational model was developed for cytochrome P450 2E1 (CYP2E1) based on inhibition constants determined for inhibition of chlorzoxazone, or 4-nitrophenol, metabolism. An empirical descriptor for type II binding was developed and tested for a series of CYP2E1 inhibitors. Inhibition constants where measured for 51 different compounds. A fast 2-dimensional predictive model was developed based on 40 compounds, and tested on 8 compounds of diverse structure. The trained model (n=40) had an r(2) value of 0.76 and an RMSE of 0.48. The correlation between the predicted and actual pK(i) values of the test set of compounds not included in the model gives an r(2) value of 0.78. The features that described binding include heme coordination (type II binding), molecular volume, octanol/water partition coefficient, solvent accessible surface area, and the sum of the atomic polarizabilities. The heme coordination parameter assigns an integer between 0 and 6 depending on structure, and is a new descriptor, based on simple quantum chemical calculations with correction for steric effects. The type II binding parameter was found to be important in obtaining a good correlation between predicted and experimental inhibition constants increasing the r(2) value from 0.38 to 0.77.

Prediction of activation energies for aromatic oxidation by cytochrome P450

We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 for a diverse set of 17 substrates using state-of-the-art density functional theory (B3LYP) with large basis sets. The activation energies vary from 60 to 87 kJ/mol. We then test if these results can be reproduced by computationally less demanding methods. The best methods (a B3LYP calculation of the activation energy of a methoxy-radical model or a partial least-squares model of the semiempirical AM1 bond dissociation energies and spin densities of the tetrahedral intermediate for both a hydroxyl-cation and a hydroxyl-radical model) give correlations with r(2) of 0.8 and mean absolute deviations of 3 kJ/mol. Finally, we apply these simpler methods on several sets of reactions for which experimental data are available and show that we can predict the reactive sites by combining calculations of the activation energies with the solvent-accessible surface area of each site.

EaMEAD: Activation energy prediction of cytochrome P450 mediated metabolism with effective atomic descriptors

In an effort to improve drug design and predictions for pharmacokinetics (PK), an empirical model was developed to predict the activation energies (Ea) of cytochrome P450 (CYP450) mediated metabolism. The model, EaMEAD (Activation energy of Metabolism reactions with Effective Atomic Descriptors), predicts the Ea of four major metabolic reactions of the CYP450 enzyme: aliphatic hydroxylation, N-dealkylation, O-dealkylation, and aromatic hydroxylation. To build and validate the empirical model, the E(a) values of the substrates with diverse chemical structures (394 metabolic sites for aliphatic hydroxylation, 27 metabolic sites for N-dealkylation, 9 metabolic sites for O-dealkylation, and 85 metabolic sites for aromatic hydroxylation) were calculated by AM1 molecular orbital (MO). Empirical equations, Quantitative Structure Activity Relationship (QSAR) models, were derived using effective atomic charge, effective atomic polarizability, and bond dipole moments of the substrates as descriptors. EaMEAD is shown to accurately predict Ea with a correlation coefficient (R) of 0.94 and root-mean-square error (RMSE, unit is kcal/mol) of 0.70 for aliphatic hydroxylation, N-dealkylation, and O-dealkylation, and R of 0.83 and RMSE of 0.80 for aromatic hydroxylation, respectively. Physical origin and the role of the effective atomic descriptors of the models are presented in detail. With this model, the Ea of the metabolism can be rapidly predicted without any experimental parameters or time-consuming QM calculation. Regioselectivity prediction with our model is presented in the case of CYP3A4 metabolism. The reliability and ease of use of this model will greatly facilitate early stage PK predictions and rational drug design. Moreover, the model can be applied to develop the Ea prediction model of various types of chemical reactions.

Purification and mechanism of human aldehyde oxidase expressed in Escherichia coli

Human aldehyde oxidase 1 (AOX1) has been subcloned into a vector suitable for expression in Escherichia coli, and the protein has been expressed. The resulting protein is active, with sulfur being incorporated in the molybdopterin cofactor. Expression levels are modest, but 1 liter of cells supplies enough protein for both biochemical and kinetic characterization. Partial purification is achieved by nickel affinity chromatography through the addition of six histidines to the amino-terminal end of the protein. Kinetic analysis, including kinetic isotope effects and comparison with xanthine oxidase, reveal similar mechanisms, with some subtle differences. This expression system will allow for the interrogation of human aldehyde oxidase structure/function relationships by site-directed mutagenesis and provide protein for characterizing the role of AOX1 in drug metabolism.

The biochemistry of drug metabolism--an introduction: Part 2. Redox reactions and their enzymes

This review continues a general presentation of the metabolism of drugs and other xenobiotics started in a recent issue of Chemistry & Biodiversity. This Part 2 presents the numerous oxidoreductases involved, their nomenclature, relevant biochemical properties, catalytic mechanisms, and the very diverse reactions they catalyze. Many medicinally, environmentally, and toxicologically relevant examples are presented and discussed. Cytochromes P450 occupy a majority of the pages of Part 2, but a large number of relevant oxidoreductases are also considered, e.g., flavin-containing monooxygenases, amine oxidases, molybdenum hydroxylases, peroxidases, and the innumerable dehydrogenases/reductases.

Visible spectra of type II cytochrome P450-drug complexes: evidence that "incomplete" heme coordination is common

The visible spectrum of a ligand-bound cytochrome P450 is often used to determine the nature of the interaction between the ligand and the P450. One particularly characteristic form of spectra arises from the coordination of nitrogen-containing ligands to the P450 heme iron. These type II ligands tend to be inhibitors because they stabilize the low reduction potential P450 and prevent oxygen binding to the heme. Yet, several type II ligands containing aniline, imidazole, and triazole moieties are also known to be substrates of P450, although P450 binding spectra are not often scrutinized to make this distinction. Therefore, the three nitrogenous ligands aniline, imidazole, and triazole were used as binding spectra standards with purified human CYP3A4 and CYP2C9, because their small size should not present any steric limitations in their accessing the heme prosthetic group. Next, the spectra of P450 with drugs containing the three nitrogenous groups were collected for comparison. The absolute spectra demonstrated that the red-shift of the low-spin Soret band is mostly dependent on the electronic properties of the nitrogen ligand since they tended to match their respective standards, aniline, imidazole, and triazole. On the other hand, difference spectra seemed to be more sensitive to the steric properties of the ligand because they facilitated comparison of the spectral amplitudes achieved with the drugs versus those with the standard nitrogen ligands. Therefore, difference spectra may help reveal "weak" coordination to the heme that results from suboptimal orientation or ligand binding to more remote locations within the P450 active sites.

Prediction of activation energies for hydrogen abstraction by cytochrome p450

We have estimated the activation energy for hydrogen abstraction by compound I in cytochrome P450 for a diverse set of 24 small organic substrates using state-of-the-art density functional theory (B3LYP). We then show that these results can be reproduced by computationally less demanding methods, for example, by using small organic mimics of compound I with both B3LYP and the semiempirical AM1 method (mean absolute error of 3-4 kJ/mol) or by calculating the bond dissociation energy, without relaxation of the radical (B3LYP) or estimated from three-point fit to a Morse potential (AM1; errors of 4 and 5 kJ/mol, respectively). We can assign activation energies of 74, 61, 53, 47, and 30 kJ/mol to primary carbons, secondary/tertiary carbons, carbons with adjacent sp(2) or aromatic groups, ethers/thioethers, and amines, respectively, which gives a very simple and predictive model. Finally, some of the less demanding methods are applied to study the CYP3A4 metabolism of progesterone and dextromethorphan.

KINETIC ISOTOPE EFFECTS IMPLICATE A SINGLE OXIDANT FOR CYTOCHROME P450-MEDIATEDO-DEALKYLATION,N-OXYGENATION, AND AROMATIC HYDROXYLATION OF 6-METHOXYQUINOLINE

One major point of controversy in the area of cytochrome P450 (P450)-mediated oxidation reactions is the nature of the active-oxygen species. A number of hypotheses have been advanced which implicate a second oxidant besides the iron-oxo species designated as compound I (Cpd 1). This oxygen is thought to be either an iron-hydroperoxy species (Cpd 0) or a second spin-state of Cpd 1. Very little information is available on what fraction of P450 oxidations is mediated by the two different oxidants. Herein, we report results on three cytochrome P450-mediated reactions: O-dealkylation, N-oxygenation, and aromatic hydroxylation, which occur by three distinct chemical mechanisms. We have used kinetic isotope effects to test for branching from O-demethylation to N-oxygenation and aromatic hydroxylation, using 6-methoxyquinoline and 2H3-6-methoxyquinoline as substrates for P4501A2. Identical large inverse isotope effects on Vmax/Km are obtained for the formation of both the N-oxide and the phenol. This indicates that all three reactions occur through the same enzyme-substrate complex and, thus, through a single iron-oxygen species. The nature of the iron-oxygen species is less certain but is more likely to be iron-oxo Cpd 1, given the energetics of these reactions.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Regiospecificity of human cytochrome P450 1A1-mediated oxidations: the role of steric effects

Cytochrome P450 1A1 oxidizes a diverse range of substrates, including the procarcinogenic xenobiotic benzo[a]pyrene (B[a]P) and endogenous fatty acid precursors of prostaglandins, such as arachidonic acid (AA) and eicosapentaenoic acid (EA). We have investigated the extent to which enzyme-substrate interactions govern regio- and stereoselectivity of oxidation of these compounds by using docking and molecular dynamics (MD) simulations to examine the likelihood of substrate oxidation at various sites. Due to structural differences between the substrates analyzed, B[a]P and its diols (planar, rigid), and the fatty acids AA and EA (long, flexible), different docking strategies were required. B[a]P, B[a]P-7,8-diols, (+) 7S,8S- and (-) 7R,8R-diols, were docked into the active site of a homology model of P450 1A1 using an automated routine, Affinity (Accelrys, San Diego, CA). AA and EA, on the other hand, required a series of restrained MD simulations to obtain a variety of productive binding modes. All complexes were evaluated by MD-based in silico site scoring to predict product profiles based on certain geometric criteria, such as angle and distance of a given substrate atom from the ferryl oxygen. For all substrates studied, the in vitro profiles were generally reflected by the in silico scores, which suggests that steric factors play a key role in determining regiospecificity in P450 1A1-mediated oxidations. We have also shown that molecular dynamics simulations may be very useful in determination of product profiles for structurally diverse substrates of P450 enzymes.

A New Statistical Approach to Predicting Aromatic Hydroxylation Sites. Comparison with Model-Based Approaches

A new approach is described that is able to predict the most probable metabolic sites on the basis of a statistical analysis of various metabolic transformations reported in the literature. The approach is applied to the prediction of aromatic hydroxylation sites for diverse sets of substrates. Training is performed using the aromatic hydroxylation reactions from the Metabolism database (Accelrys). Validation is carried out on heterogeneous sets of aromatic compounds reported in the Metabolite database (MDL). The average accuracy of prediction of experimentally observed hydroxylation sites estimated for 1552 substrates from Metabolite is 84.5%. The proposed approach is compared with two electronic models for P450 mediated aromatic hydroxylation: the oxenoid model using the atomic oxygen and the model using the methoxy radical as a model for the heme active oxygen species. For benzene derivatives, the proposed method is inferior to the oxenoid model and as accurate as the methoxy-radical model. For hetero- and polycyclic compounds, the oxenoid model is not applicable, and the statistical method is the most accurate. Broad applicability and high speed of calculations provide the basis for using the proposed statistical approach for high-throughput metabolism prediction in the early stages of drug discovery.