Kinetic study of cytochrome P450 3A4 activity on warfarin by capillary electrophoresis with fluorescence detection

The use of capillary electrophoresis (CE) for the determination of cytochrome P450 3A4 (CYP3A4) activity with R-warfarin as a substrate was investigated. CYP3A4 activity was determined by the quantitation of the product, 10-hydroxywarfarin, based on separation by CE. The separation conditions were as follows: capillary, 80.5 cm (75 microm i.d., 60 cm effective length); 50 mM sodium phosphate buffer (pH 6.5); 23 kV (90 microA) applied voltage; fluorescence detection, excitation wavelength, 310 nm, emission wavelength, 418 nm; capillary temperature, 37 degrees C. With the developed CYP3A4 activity assay and the Lineweaver-Burk equation, the Michaelis-Menten parameters Km and Vmax for formation of 10-hydroxywarfarin from R-warfarin in the presence of CYP3A4 were calculated to be 166 +/- 12 microM and 713 +/- 14 pmol/min/nmol (or 91.4 pmol/min/mg) CYP3A4, respectively.     

Molecular modeling of cytochrome P450 3A4

The three-dimensional structure of human cytochrome P450 3A4 was modeled based on crystallographic coordinates of four bacterial P450s: P450 BM-3, P450cam, P450terp, and P450eryF. The P450 3A4 sequence was aligned to those of the known proteins using a structure-based alignment of P450 BM-3, P450cam, P450terp, and P450eryF. The coordinates of the model were then calculated using a consensus strategy, and the final structure was optimized in the presence of water. The P450 3A4 model resembles P450 BM-3 the most, but the B helix is similar to that of P450eryF, which leads to an enlarged active site when compared with P450 BM-3, P450cam, and P450terp. The 3A4 residues equivalent to known substrate contact residues of the bacterial proteins and key residues of rat P450 2B1 are located in the active site or the substrate access channel. Docking of progesterone into the P450 3A4 model demonstrated that the substrate bound in a 6-orientation can interact with a number of active site residues, such as 114, 119, 301, 304, 305, 309, 370, 373, and 479, through hydrophobic interactions. The active site of the enzyme can also accommodate erythromycin, which, in addition to the residues listed for progesterone, also contacts residues 101, 104, 105, 214, 215, 217, 218, 374, and 478. The majority of 3A4 residues which interact with progesterone and/or erythromycin possess their equivalents in key residues of P450 2B enzymes, except for residues 297, 480 and 482, which do not contact either substrate in P450 3A4. The results from docking of progesterone and erythromycin into the enzyme model make it possible to pinpoint residues which may be important for 3A4 function and to target them for site-directed mutagenesis.     

Regioselectivity preference of testosterone hydroxylation by cytochrome P450 3A4

Theoretical studies are presented into the experimentally observed regioselectivity difference of testosterone hydroxylation by cytochrome P450 3A4 at the 1β, 2β, 6β, and 15β positions. Such regioselectivity is investigated by density functional theory calculations on a model system. The barrier heights of hydrogen abstraction, which are corrected by zero-point vibrational energies, are computed to be about 10.1, 13.6, 14.4, and 16.2 kcal/mol for the 6β-, 2β-, 15β-, and 1β-positions, respectively. The calculated barriers suggest the regioselectivity preference of 6β ? 2β > 15β > 1β, which is in good agreement with experimental findings.     

Chiral capillary electrophoretic analysis of verapamil metabolism by cytochrome P450 3A4

Cytochrome P450 (CYP), which is one of the most important enzymes in human liver, is responsible for a large portion of the first-pass metabolism of drugs. Many studies have focused on the determination of CYP activity by substrate assays. Most of them used liquid chromatography (LC) as analytical technique, while only a few studies used capillary electrophoresis (CE) for the separation and quantitation of reaction components. In this study, the feasibility of using CE in an in vitro metabolism study with CYP was tested. Verapamil was chosen as the substrate for CYP 3A4 isozyme (Supersome). A chiral capillary electrophoretic method was developed and validated for the simultaneous determination of R,S-verapamil (VER) and their major metabolites, R,S-norverapamil (NOR). A method for CYP 3A4 activity assay was proposed with VER as a probe. At the same time, the enantioselective metabolism of VER was studied. Michaelis-Menten constants of R- and S-VER were determined. S-VER was metabolised faster and more extensively than R-VER, with K(m)=167+/-23 microM, V(max)=3,418+/-234 pmol/min/mg for S-VER, and K(m)=168+/-35 microM, V(max)=2,502+/-275 pmol/min/mg for R-VER.     

QSAR modeling of in vitro inhibition of cytochrome P450 3A4

We report the QSAR modeling of cytochrome P450 3A4 (CYP3A4) enzyme inhibition using four large data sets of in vitro data. These data sets consist of marketed drugs and drug-like compounds all tested in four assays measuring the inhibition of the metabolism of four different substrates by the CYP3A4 enzyme. The four probe substrates are benzyloxycoumarin, testosterone, benzyloxyresorufin, and midazolam. We first show that using state-of-the-art QSAR modeling approaches applied to only one of these four data sets does not lead to predictive models that would be useful for in silico filtering of chemical libraries. We then present the development and the testing of a multiple pharmacophore hypothesis (MPH) that is formulated as a conceptual extension of the traditional QSAR approach to modeling the promiscuous binding of a large variety of drugs to CYP3A4. In the simplest form, the MPH approach takes advantage of the multiple substrate data sets and identifies the binding of test compounds as either proximal or distal relative to that of a given substrate. Application of the approach to the in silico filtering of test compounds for potential inhibitors of CYP3A4 is also presented. In addition to an improvement in the QSAR modeling for the inhibition of CYP3A4, the results from this modeling approach provide structural insights into the drug-enzyme interactions. The existence of multiple inhibition data sets in the BioPrint database motivates the original development of the concept of a multiple pharmacophore hypothesis and provides a unique opportunity for formulating alternative strategies of QSAR modeling of the inhibition of the in vitro metabolism of CYP3A4.     

Structural and dynamical basis of broad substrate specificity, catalytic mechanism, and inhibition of cytochrome P450 3A4

Cytochrome P450 (CYP) 3A4 is responsible for the oxidative degradation of more than 50% of clinically used drugs. By means of molecular dynamics simulations with the newly developed force field parameters for the heme-thiolate group and its dioxygen adduct, we examine the differences in structural and dynamic properties between CYP3A4 in the resting form and its complexes with the substrate progesterone and the inhibitor metyrapone. The results indicate that the broad substrate specificity of CYP3A4 stems from the malleability of a loop (residues 211-218) that resides in the vicinity of the channel connecting the active site and bulk solvent. However, the high-amplitude motion of the flexible loop is found to be damped out upon binding of the inhibitor or the substrate in the active site. In the resting form of CYP3A4, a structural water molecule is bound to the sixth coordination position of the heme iron, stabilizing the octahedral coordination geometry. In addition to the direct coordination of metyrapone to the heme iron, the hydrogen bond interaction between the inhibitor carbonyl group and the side chain of Ser119 also contributes significantly to stabilizing the CYP3A4-metyrapone complex. On the other hand, progesterone is stabilized in the active site by the formation of two hydrogen bonds with Ser119 and Arg106, as well as by the van der Waals interactions with the heme and hydrophobic residues. The structural and dynamic features of the CYP3A4-progesterone complex indicate that the oxidative degradation of progesterone occurs through hydroxylation at the C16 position by the reactive oxygen coordinated to the heme iron.     

Structural determinants of steroids for cytochrome P450 3A4-mediated metabolism

Cytochrome P450 3A4 (CYP3A4), a major member of cytochrome P450 isoenzymes, metabolizes the majority of steroids in 6β-position. For the purpose of determining requisite structural features of a series of structurally related steroids for CYP3A4-mediated metabolism, three-dimensional pharmacophore modeling as well as electrotopological state map were conducted for 15 steroids. Though prior studies speculated that the chemical reactivity of the allylic 6β-position might have a greater influence on CYP3A4 selective 6-hydroxylation than steric constraints in the enzyme, our results reveal that for CYP3A4 steroidal substrates, it is not the chemical reactivity of atoms at 6β-site, but the pharmacophoric features, i.e. the two hydrophobic rings together with two H-bond donors, that act as the key factors responsible for determining the CYP3A4 selective 6-hydroxylation of steroids.     

Resolution of rabbit liver cytochrome P450 3b and 3c (P450 IIC3 and IIIA4) from microsomal membrane by two-dimensional gel electrophoresis

Two-dimensional gel electrophoresis was applied to the separation of microsomal proteins and P450 isozymes. The experimental conditions of the first-dimensional isoelectric focusing were optimized. The best results were obtained under the following conditions: (i) a mixture of two detergents (Tergitol NP 10 and CHAPS) in the focusing gel and sample buffer, (ii) a carrier ampholyte gradient between pH 6.5 and 9.0, (iii) sample loading at the anodic end of the gel, and (iv) low protein loading in the sample (below 50 micrograms). Four P450 forms, including P450 3b and 3c, could be resolved under these conditions.     

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.     

Biotransformation of the sesquiterpene (+)-valencene by cytochrome P450cam and P450BM-3

The sesquiterpenoids are a large class of naturally occurring compounds with biological functions and desirable properties. Oxidation of the sesquiterpene (+)-valencene by wild type and mutants of P450cam from Pseudomonas putida, and of P450BM-3 from Bacillus megaterium, have been investigated as a potential route to (+)-nootkatone, a fine fragrance. Wild type P450cam did not oxidise (+)-valencene but the mutants showed activities up to 9.8 nmol (nmol P450)(-1) min(-1), with (+)-trans-nootkatol and (+)-nootkatone constituting >85% of the products. Wild type P450BM-3 and mutants had higher activities (up to 43 min(-1)) than P450cam but were much less selective. Of the many products, cis- and trans-(+)-nootkatol, (+)-nootkatone, cis-(+)-valencene-1,10-epoxide, trans-(+)-nootkaton-9-ol, and (+)-nootkatone-13S,14-epoxide were isolated from whole-cell reactions and characterised. The selectivity patterns suggest that (+)-valencene has one binding orientation in P450cam but multiple orientations in P450BM-3.     

Defining CYP3A4 structural responses to substrate binding. Raman spectroscopic studies of a nanodisc-incorporated mammalian cytochrome P450

Resonance Raman (RR) spectroscopy is used to help define active site structural responses of nanodisc-incorporated CYP3A4 to the binding of three substrates: bromocriptine (BC), erythromycin (ERY), and testosterone (TST). We demonstrate that nanodisc-incorporated assemblies reveal much more well-defined active site RR spectroscopic responses as compared to those normally obtained with the conventional, detergent-stabilized, sampling strategies. While ERY and BC are known to bind to CYP3A4 with a 1:1 stoichiometry, only the BC induces a substantial conversion from low- to high-spin state, as clearly manifested in the RR spectra acquired herein. The third substrate, TST, displays significant homotropic interactions within CYP3A4, the active site binding up to 3 molecules of this substrate, with the functional properties varying in response to binding of individual substrate molecules. While such behavior seemingly suggests the possibility that each substrate binding event induces functionally important heme structural changes, up to this time spectroscopic evidence for such structural changes has not been available. The current RR spectroscopic studies show clearly that accommodation of different size substrates, and different loading of TST, do not significantly affect the structure of the substrate-bound ferric heme. However, it is here demonstrated that the nature and number of bound substrates do have an extraordinary influence on the conformation of bound exogenous ligands, such as CO or dioxygen and its reduced forms, implying an effective mechanism whereby substrate structure can impact reactivity of intermediates so as to influence function, as reflected in the diverse reactivity of this drug metabolizing cytochrome.     

Observation of ligand binding to cytochrome P450 BM-3 by means of solid-state NMR spectroscopy

A broad understanding of the binding modes of ligands and inhibitors to cytochrome P450 is vital for the development of new drugs. We investigated ligand binding in a site-specific fashion on cytochrome P450 BM-3 from Bacillus megaterium, a 119 kDa paramagnetic enzyme, using solid-state magic angle spinning nuclear magnetic resonance methods. Selective labeling and longitudinal relaxation effects were utilized to identify the peaks in a site-specific fashion and to provide evidence for binding. Well-resolved one-dimensional and two-dimensional NMR spectra of cytochrome P450 BM-3 reveal shifts upon binding of its substrate, N-palmitoylglycine. These data are consistent with the crystallographic result that a biochemically important amino acid residue, Phe87, moves upon ligation. This experimental scheme provides a tool for probing ligand binding for complex systems.     

Pyrene.pyrene complexes at the active site of cytochrome P450 3A4: evidence for a multiple substrate binding site

Cytochrome P450 monooxygenases (CYPs) metabolize nearly all drugs and toxins. Recently, it has become clear that CYPs exhibit both homotropic and heterotropic allosteric kinetics for many substrates. However, the mechanism of cooperative kinetics has not been established for any specific human CYP/substrate combination. Suggested mechanisms include binding of multiple substrates within distinct, static, subsites of a single large active site or binding of multiple substrates within a single fluid active site. CYP3A4 hydroxylates pyrene with positive cooperativity. Therefore, experiments were designed to exploit the fluorescence properties of pyrene, which diagnostically distinguish between pyrene.pyrene complexes versus spatially separated pyrene substrates. Pyrene complexes (excimers) yield an emission spectrum clearly distinct from pyrene monomers. In lipid-free aqueous/glycerol solutions of CYP3A4, addition of pyrene affords a concentration-dependent low-spin to high-spin conversion of the CYP3A4 heme prosthetic group, indicating occupancy of the active site by pyrene. Under the same conditions, in the presence of CYP3A4 but not other heme proteins, the excimer/monomer ratio (E/M) of pyrene was decreased in emission spectra, compared to pyrene alone. However, excitation spectra indicate a CYP3A4-dependent increase in the wavelength shift for the excimer excitation spectrum versus the monomer excitation spectrum, as well as changes in the excimer excitation peak shape and vibronic structure. These changes are reversed by the CYP3A4 substrate testosterone. Together, the results demonstrate that pyrene.pyrene ground-state complexes occupy the CYP3A4 active site, and they provide the first spectroscopic evidence for substrate complexes within a single fluid active site. Functional implications include the possibility that turnover rate, regioselectivity, and stereoselectivity of the reaction are determined by the substrate.substrate complex rather than individual substrates.     

P450 versus P420: correlation between cyclic voltammetry and visible absorption spectroscopy of the immobilized heme domain of cytochrome P450 BM3

Cyclic voltabsorptometry is used for the first time to distinguish and characterize electrochemically the active (P450) and inactive (P420) forms of cytochromes P450 immobilized on an electrode during voltammetry experiments. This was achieved by using the heme domain (BMP) of the bacterial cytochrome P450 BM3 from Bacillus megaterium (CYP102A1) immobilized on mesopouros tin-oxide (SnO2) electrodes. We demonstrate that the formation of either the P450 form or the P420 one can be obtained by modifying the mesoporous electrode surface with polycations with different properties such as polyethylenimmine (PEI) and polydiallyldimethylammonium chloride (PDDA). Potential step spectroelectrochemistry allowed measurement of reduction potentials of the active P450 form. Values of -0.39+/-0.01 V and -0.58+/-0.01 V (both versus Ag/AgCl) were calculated for the active P450 form immobilized on the BMP/PDDA-SnO2 and BMP/PEI-SnO2 electrodes, respectively. The cyclic voltabsorptometric experiments showed how, when both the active and inactive forms are present on the PEI film, the inactive P420 species tends to dominate the cyclic voltammetric signal.     

Aromatic hydroxylation by cytochrome P450: model calculations of mechanism and substituent effects

The mechanism and selectivity of aromatic hydroxylation by cytochrome P450 enzymes is explored using new B3LYP density functional theory computations. The calculations, using a realistic porphyrin model system, show that rate-determining addition of compound I to an aromatic carbon atom proceeds via a transition state with partial radical and cationic character. Reactivity is shown to depend strongly on ring substituents, with both electron-withdrawing and -donating groups strongly decreasing the addition barrier in the para position, and it is shown that the calculated barrier heights can be reproduced by a new dual-parameter equation based on radical and cationic Hammett sigma parameters.     

Determination of midazolam and its metabolite as a probe for cytochrome P450 3A4 phenotype by liquid chromatography-mass spectrometry

This study demonstrated the analysis of midazolam and its metabolites by liquid chromatography-mass spectrometry (LC-MS) with a sonic spray ionization (SSI) interface. The analytical column was a YMC-Pak Pro C18 (50 mm x 2.0 mm i.d.) using 10 mM ammonium acetate (pH 4.8)-methanol (1:1) at a flow rate of 0.2 ml min(-1). The drift voltage was 100 V. The sampling aperture was heated at 110 degrees C and the shield temperature was 230 degrees C. The lower limits for the detection of midazolam and 1'-hydroxymidazolam were 26.3 and 112.76 pg injected, respectively. The calibration curves for midazolam and 1'-hydroxymidazolam were linear in the range of 0.1-5 microg ml(-1). Within-day relative standard deviations was less than 7%. The method was applied to the determination of midazolam in monkey plasma, and the analysis of midazolam and its metabolites in an in vitro study with recombinant cytochrome P450 (CYP) 3A4. This method is sufficiently sensitive and useful to elucidate the kinetics of midazolam metabolite formation. We also investigated the effect of propofol on the metabolism of midazolam using recombinant CYP3A4. Propofol competitively inhibited the metabolism of midazolam to 1'-hydroxymidazolam by CYP3A4.     

New combined model for the prediction of regioselectivity in cytochrome P450/3A4 mediated metabolism

Cytochrome P450 3A4 metabolizes nearly 50% of the drugs currently in clinical use with a broad range of substrate specificity. Early prediction of metabolites of xenobiotic compounds is crucial for cost efficient drug discovery and development. We developed a new combined model, MLite, for the prediction of regioselectivity in the cytochrome P450 3A4 mediated metabolism. In the model, the ensemble catalyticphore- based docking method was implemented for the accessibility prediction, and the activation energy estimation method of Korzekwa et al. was used for the reactivity prediction. Four major metabolic reactions, aliphatic hydroxylation, N-dealkylation, O-dealkylation, and aromatic hydroxylation reaction, were included and the reaction data, metabolite information, were collected for 72 well-known substrates. The 47 drug molecules were used as the training set, and the 25 well-known substrates were used as the test set for the ensemble catalyticphore-based docking method. MLite predicted correctly about 76% of the first two sites in the ranking list of the test set. This predictability is comparable with that of another combined model, MetaSite, and the recently published QSAR model proposed by Sheridan et al. MLite also offers information about binding configurations of the substrate-enzyme complex. This may be useful in drug modification by the structure-based drug design.