Comparison of high pressure and nanoscale confinement effects on crystallization of the molecular glass-forming liquid, dimethyl phthalate.

High pressure and nanoscopic confinement are two different strategies commonly employed to modify the physicochemical properties of various materials. Both strategies act mostly by changing the molecular packing. In this work, we performed a comparative study on the effect of compression and confined geometry on crystallization of a molecular liquid. Dielectric spectroscopy was employed to investigate the crystallization of the van der Waals liquid, dimethyl phthalate, in nanoporous alumina of different pore sizes as well as on increased pressure (up to 200 MPa). The analysis of the crystallization kinetics under varying thermodynamic conditions revealed that both strategies affect the crystallization behavior of the sample in very distinct ways. Compression shifts the maximum crystallization rate towards a higher temperature and broadens it. As a result, it is more challenging to avoid crystallization upon cooling the liquid at high pressure. In contrast, when the same material is incorporated into nanopores, crystallization significantly slows down and the maximum rate shifts towards a lower temperature with decreasing pore size. Finally, we show that crystallization in nanoporous alumina is accompanied by pre-crystallization effects upon which a shift of the α-relaxation peak is observed. An equilibration process prior to the initiation of crystallization was detected for the confined material both above and below the glass transition temperature of the interfacial layer, while not in the bulk.

Studying tautomerism in an important pharmaceutical glibenclamide confined in the thin nanometric layers.

The uniform thin films with variable thicknesses (d = 49, 120, 220 nm) of active pharmaceutical ingredient (API) glibenclamide (GCM) was spin-coated and investigated using broadband dielectric, grazing incident FTIR spectroscopies, atomic force microscopy, and ellipsometry. Data analysis revealed that nanoconfined systems consist of a mixture of amide and imidic acid forms of this pharmaceutical, wherein the ratios of both tautomeric forms in the thin films were different with respect to the molten supercooled bulk system. Moreover, changes in the populations of glibenclamide tautomers, i.e. higher amide to imides ratio in the spatially restricted API with respect to the bulk sample, had a strong impact on the character of the proton transfer reaction. In this context, the kinetic curves constructed on the base of infrared data for the bulk system follow the sigmoidal shape, characteristic for the autocatalytic reaction, while results obtained for the confined samples provide exponential character and indicate first-order transformation. This allows hypothesizing that the autocatalytic nature of the tautomerism in the bulk sample is most likely related to the formation of the amide tautomers which further catalyze the progress of imide-amide transformation. Our results are the first studies showing that the change in the thickness of the film may affect the properties and isomerization kinetics in a pharmaceutical systems. Finally, our data open a new perspective for developing new drug delivery systems.

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.

Regulation of proteins in the cholesterol side-chain cleavage system in JEG-3 and Y-1 cells.

The conversion of cholesterol to pregnenolone, the rate-limiting step in steroid hormone synthesis, occurs on mitochondrial cytochrome P450scc, which catalyzes this reaction by receiving electrons from NADPH via a flavoprotein [adrenodoxin reductase (AdRed)] and an iron sulfur protein [adrenodoxin (Adx)]. The behavior of the genes and mRNAs encoding these proteins has been studied in several systems, but little is known about the behavior of the human proteins. Using cloned cDNAs for human P450scc and AdRed, we constructed bacterial expression vectors to make milligram quantities of the corresponding proteins. These, plus purified human Adx similarly prepared by Dr. L. Vickery, were injected into rabbits to raise antiserum to each of the proteins. Each antiserum was highly specific and did not cross-react with other mitochondrial proteins detectable by Western blotting. Human JEG-3 choriocarcinoma cells and mouse Y-1 adrenocortical carcinoma cells were then incubated for 0-24 h with 1 mM 8-bromo-cAMP (8Br-cAMP) or 30 nM phorbol 12-myristate 13-acetate (PMA; phorbol ester) plus 1 microM A23187 (calcium ionophore) to activate the protein kinase-A and -C pathways, respectively. In JEG-3 cells, 8Br-cAMP increased and PMA/A23187 slightly decreased the abundance of P450scc and Adx, but neither treatment had a detectable effect on AdRed. The production of pregnenolone by these cells increased 3-fold in response to 8Br-cAMP and fell to one third in response to PMA/A23187. In Y-1 cells, 8Br-cAMP increased the abundance of all three proteins, while PMA/A23187 decreased the abundance of P450scc and Adx. The production of pregnenolone by these cells increased 9-fold in response to 8Br-cAMP and was unaffected by TPA/A23187. These studies show that the three proteins of the cholesterol side-chain cleavage system behave in response to 8Br-cAMP and PMA/A23187 as predicted from the study of their genes and mRNAs, indicating that the chronic regulation of steroidogenesis in these cell systems is regulated principally at the level of mRNA abundance.

Construction and function of fusion enzymes of the human cytochrome P450scc system.

Type I cytochrome P450 enzyme systems are found in mitochondria and consist of three components, a flavoprotein (adrenodoxin reductase, AdRed), an iron-sulfur protein (adrenodoxin, Adx), and the cytochrome P450; Type II P450 enzymes in the endoplasmic reticulum consist of only two components, P450 reductase and the P450. Genetically engineered fusion proteins of Type II cytochromes P450 (such as steroid 17 alpha- and 21-hydroxylases) produce enzymes with increased activity. To test the consequences of constructing fusions of Type I enzymes, we built fusion proteins based on the cholesterol side-chain cleavage enzyme, P450scc. We constructed expression vectors for three fusion proteins: NH2-P450scc-AdRed-COOH, P450-AdRed-Adx, and P450scc-Adx-AdRed. The various components were assembled from cassette-like cDNA fragments modified and amplified by polymerase chain reaction (PCR), subcloned into a specially tailored vector, and linked by DNA segments encoding hydrophilic linker peptides. The final vectors were transfected into COS-1 cells, incubated with 22R-hydroxycholesterol, and assayed by the secretion of pregnenolone into the culture medium. Triple transfection of three individual vectors expressing P450scc, AdRed, and Adx yielded more pregnenolone than did transfection with P450scc alone. The P450scc-AdRed and P450scc-Adx-AdRed fusion proteins produced levels of pregnenolone similar to the control triple transfection. However, the P450scc-AdRed-Adx fusion produced substantially more pregnenolone, having an apparent Vmax of 9.1 ng of pregnenolone produced per milliliter of medium per 24 hr, compared to a Vmax of 1.7 ng/ml per day for the triple transfection.(ABSTRACT TRUNCATED AT 250 WORDS)

Use of homology modeling in conjunction with site-directed mutagenesis for analysis of structure-function relationships of mammalian cytochromes P450.

In recent years, homology modeling has become an important tool to study cytochrome P450 function, especially in conjunction with experimental approaches such as site-directed mutagenesis. Molecular models of mammalian P450s can be constructed based on crystal structures of four bacterial enzymes, P450cam, P450 BM-3, P450terp and P450eryF, using molecular replacement or consensus methods. In a model built by molecular replacement, the coordinates are copied from those of a given template protein, while consensus methods utilize more then one protein as a template and are based on distance geometry calculations. The models can be used to identify or confirm key residues, evaluate enzyme-substrate interactions and explain changes in protein stability and/or regio- and stereospecificity of substrate oxidation upon residue substitution by site-directed mutagenesis. P450 models have also been utilized to analyze binding of inhibitors or activators, as well as alterations in inhibition and activation due to residue replacement.

Application of 3-dimensional homology modeling of cytochrome P450 2B1 for interpretation of site-directed mutagenesis results.

Three-dimensional structures of cytochrome P450 2B1 were modeled based on the crystallographic structure of P450cam. The effect of the alignment, loop choice, and minimization with or without water was assessed. Although final models were similar in overall structure, the identity of active site residues depended upon the alignment. An example is Phe-206, which may or may not form part of the active site. The choice of the loop conformation had a lesser effect, while including water in the final minimization step was essential for preserving the shape and size of the active site. The best model (model 2) was in good agreement with the data from site-directed mutagenesis studies, and correctly predicted the effect of substitutions at 9 out of 10 amino acid positions. Thus, residues important for P450 2B1 activity, such as Ile-114, Phe-206, Ile-290, Thr-302, Val-363, and Gly-478, constitute part of the active site and are able to interact with the substrate androstenedione through hydrophobic interactions. On the other hand, Ser-303, Ser-360 and Lys-473 are far from the active site and/or cannot interact with the substrate, in agreement with experimental data. The model indicates other residues likely to be important for enzyme function, such as Tyr-111, Leu-209, Ile-477, and Ile-480, which can be tested experimentally. The substrate may assume numerous binding orientations consistent with observed patterns of hydroxylation at C15 and C16. The replacement in the model of certain amino acid residues to mimic residue substitutions from site-directed mutagenesis studies and docking of the substrate into the modified active site allowed a plausible explanation for alterations in regio- and stereospecificities of some mutants of P450 2B1, such as Gly-478-->Ala or Val-363-->Ala.

First results of the search for neutrinoless double-beta decay with the NEMO 3 detector.

The NEMO 3 detector, which has been operating in the Fréjus underground laboratory since February 2003, is devoted to the search for neutrinoless double-beta decay (beta beta 0v). The half-lives of the two neutrino double-beta decay (beta beta 2v) have been measured for 100Mo and 82Se. After 389 effective days of data collection from February 2003 until September 2004 (phase I), no evidence for neutrinoless double-beta decay was found from approximately 7 kg of 100Mo and approximately 1 kg of 82Se. The corresponding limits are T1/2(beta beta0v) > 4.6 x 10(23) yr for 100Mo and T1/2(beta beta 0v) > 1.0 x 10(23) yr for 82Se (90% C.L.). Depending on the nuclear matrix element calculation, the limits for the effective Majorana neutrino mass are < 0.7-2.8 e/v for 100Mo and < 1.7-4.9 eV for 82Se.

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 beta-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.

Molecular basis of P450 inhibition and activation: implications for drug development and drug therapy.

Three-dimensional homology models of cytochromes P450 (P450) 2B1 and P450 3A4 have been utilized along with site-directed mutagenesis to elucidate the molecular determinants of substrate specificity. Most of the key residues identified in 2B enzymes fall within five substrate recognition sites (SRSs) and have counterparts in bacterial P450 residues that regulate substrate binding or access. Docking of inhibitors into 2B models has provided a plausible explanation for changes in susceptibility to mechanism-based inactivation that accompany particular amino acid side-chain replacements. These studies provide a basis for predicting drug interactions due to P450 inhibition and for rational inhibitor design. In addition, the location of P450 3A4 residues capable of influencing homotropic stimulation by substrates and heterotropic stimulation by flavonoids has been identified. Steroid hydroxylation by the wild-type enzyme exhibits sigmoidal kinetics, indicative of positive cooperativity. Based on the 3A4 model and single-site mutants, a double mutant in SRS-2 has been constructed that exhibits normal Michaelis-Menten kinetics. Results of modeling and mutagenesis studies suggest that the substrate and effector bind at adjacent sites within a single large cavity in P450 3A4. A thorough understanding of the location and structural requirements of the substrate-binding and effector sites in cytochrome P450 3A4 should prove valuable in rationalizing and predicting interactions among the multitude of drugs and other compounds that bind to the enzyme.

Phosphate exchange between high-energy phosphate compounds in resting crustacean muscle.

Single giant muscle fibres from Megabalamus psittacus was microinjected with ATP-32P or internally perfused with a solution containing 32P-PArg in an attempt to demonstrate, in the absence of inhibitiors, ATP utilization during mechanical activity with negative results. The rate constant of transfer of 32P from microinjected 32P-PArg to endogenous ATP was studied at two different temperatures, obtaining a Q10 of 1.08 +/- .012. It is concluded that the rate limiting step in the process of transference of phosphate in the resting muscle from balanus is the intracellular diffusion of PArg and not the rate at which the equivalent of Lohman's reaction in invertebrates proceeds.

Molecular modeling of mammalian cytochromes P450: application to study enzyme function.

Cytochromes P450 are important heme-containing enzymes that catalyze the oxidation of a vast array of endogenous and exogenous compounds, including drugs and carcinogens. One of the more successful approaches to study P450 function involves molecular modeling. Because none of the mammalian P450s have been crystallized, a number of homology models have been constructed based on the structures of known bacterial P450s. Molecular models, generated using molecular replacement or distance geometry methods, can be used to dock substrates and/or inhibitors in the active site to explain various aspects of enzyme function. The majority of modeling research has dealt with enzyme-substrate interactions in the active site. The analysis of these interactions has helped us to better understand the mechanism of P450 catalysis and provided the structural basis for the regio- and stereospecificity of substrate oxidation as well as susceptibility to inhibition or inactivation. The models have been utilized to identify and/or confirm key residues and to rationally interpret experimental data. The alteration in activity in a mutant P450 can be related to changes in enzyme-substrate/inhibitor interactions, such as the removal or appearance of van der Waals overlaps or changes in compound mobility. Homology models can also help to analyze P450-redox partner interactions and identify critical determinants of protein stability. We can expect further development of molecular modeling methods and their increasing contribution into research on P450 function as an integral part of a combined theoretical-experimental approach.

Site-directed mutagenesis as a tool for molecular modeling of cytochrome P450 2B1.

Prompted by our previous homology model of cytochrome P450 2B1 based on the 3-D structure of P450cam [Szklarz, G. D., Ornstein, R. L., & Halpert, J. R. (1994) J. Biomol. Struct. Dyn. 12, 61-78], we constructed 11 new site-directed mutants at positions 100, 111, 205, 209, 291, 477, and 480 and expressed the enzymes in Escherichia coli. The mutations at positions 209, 477, and 480 affected androstenedione and progesterone hydroxylation as predicted by the model. For example, the Ile-477-->Ala and Ile-480-->Ala mutants retained < or = 5% activity with androstenedione and progesterone but were active with benzphetamine, whereas the Leu-209-->Ala mutant catalyzed 21-hydroxylation of progesterone. Mutations at the other positions, i.e., 100, 111, 205, and 291, did not change enzyme activity, contrary to predictions. Therefore, an improved molecular model of cytochrome P450 2B1 was constructed. An alignment of the P450 2B1 sequence with P450 BM-3, P450cam, and P450terp was optimized using data from site-directed mutagenesis at 27 positions in various cytochromes P450 2B and docking of androstenedione into the active site of the known crystal structures. Because all three structures were found to be suitable templates for P450 2B1, the new model was formulated on the basis of the crystallographic coordinates of the three proteins using a consensus strategy, a modeling method based on distance geometry calculations. The new model provides a means to explain alterations in regio- and stereospecificity of steroid hydroxylation upon residue substitution at key amino acid positions, including positions 114, 206, 209, 290, 302, 363, 367, 477, 478, and 480 in P450 2B1.

Interconversion of the androstenedione hydroxylase specificities of cytochromes P450 2B4 and 2B5 upon simultaneous site-directed mutagenesis of four key substrate recognition residues.

Based on recent studies of single reciprocal mutants of cytochrome P450 2B4 and the highly related P450 2B5 at positions 114, 294, 363, and 367 [G. D. Szklarz, Y. Q. He, K. M. Kedzie, J. R. Halpert, and V. L. Burnett (1996) Arch. Biochem. Biophys. 327,308-318], a number of multiple mutants were constructed, expressed in Escherichia coli, and assayed with androstenedione, progesterone, and benzyloxyresorufin. Simultaneous substitutions of Ile-114, Ser-294, Ile-363, and Val-367 in cytochrome P450 2B4 with Phe, Thr, Val, and Ala, respectively from 2B5, resulted in a marked increase in androstenedione 15alpha- and 16alpha-hydroxylation compared with the wild-type enzyme and yielded a metabolite profile indistinguishable from that of cytochrome P450 2B5. Likewise, the reciprocal P450 2B5 quadruple mutant exhibited the specificity for 16beta-hydroxylation characteristic of the 2B4 wild type. The two reciprocal quadruple mutants of P450 2B4 and 2B5 also displayed benzyloxyresorufin dealkylase activities similar to those of the wild-type P450 2B5 and 2B4, respectively. However, the progesterone metabolite profile of P450 2B5 was not identical to that of the 2B4 quadruple mutant or of a quintuple mutant in which residue 370 was also mutated to the 2B5 residue. Therefore, the 17beta-acetyl group on progesterone as opposed to the oxo group on androstenedione may lead to interaction with additional amino acid residues.

The mitochondrial environment is required for activity of the cholesterol side-chain cleavage enzyme, cytochrome P450scc.

Steroidogenesis is initiated by the conversion of cholesterol to pregnenolone by mitochondrial cytochrome P450scc [cholesterol, reduced-adrenal-ferredoxin:oxygen oxidoreductase (side-chain-cleaving); EC 1.14.15.6]. Several subsequent steroidal conversions occur in the endoplasmic reticulum (ER), but the last step in the production of glucocorticoids and mineralocorticoids again occurs in the mitochondria. Although cellular compartmentalization of steroidogenic enzymes appears to be a feature of all steroidogenic pathways, some reports indicate that cholesterol can be converted to pregnenolone outside the mitochondria. To investigate whether P450scc can function outside the mitochondria, we constructed vectors producing P450scc and various fusion enzymes of P450scc with electron-transport proteins and directed their expression to either the ER or the mitochondria. Whether targeted to mitochondria or to the ER, plasmid vectors encoding P450scc and fusion proteins of P450scc with either mitochondrial or microsomal electron-transport proteins produced immunodetectable protein. When expressed in mitochondria, all of these constructions converted 22-hydroxycholesterol to pregnenolone, but when expressed in the ER none of them produced pregnenolone. These results show that P450scc can function only in the mitochondria. Furthermore, it appears to be the mitochondrial environment that is required, rather than the specific mitochondrial electron-transport intermediates.

Identification of three key residues in substrate recognition site 5 of human cytochrome P450 3A4 by cassette and site-directed mutagenesis.

Cassette mutagenesis and site-directed mutagenesis were used to investigate the importance of individual amino acid residues at positions 364-377 of cytochrome P450 3A4 in determining steroid hydroxylation or stimulation by alpha-naphthoflavone. The mutants were expressed in an Escherichia coli system, and solubilized membranes were prepared. All mutants except R365G and R365K exhibited anti-3A immunoreactivity on Western blotting, although R372S and R375K were not detected as the Fe2+-CO complex. Replacement of Arg-372 by Lys yielded a typical P450 spectrum. The results indicate that the highly conserved Arg residues at positions 365 and 375 may play a role in stabilizing the tertiary structure or in heme binding. Catalytic activities of 12 mutants were examined using progesterone and testosterone as substrates, and residues 369, 370, and 373 were found to play an important role in determining substrate specificity. Although the three mutants hydroxylated progesterone and testosterone primarily at the 6beta-position like the wild-type, replacement of Ile-369 by Val suppressed progesterone 16alpha-hydroxylase activity, whereas substitution of Ala-370 with Val enhanced progesterone 16alpha-hydroxylation. Interestingly, substitution of Leu-373 with His resulted in production of a new metabolite from both steroids. Moreover, the mutants at positions 369 and 373 were more and less responsive, respectively, than the wild-type to alpha-naphthoflavone stimulation. Alterations in activities or expression of several mutants were interpreted using a three-dimensional model of P450 3A4. The results suggest that analogy with mammalian family 2 and bacterial cytochromes P450 can be used to predict P450 3A residues that contribute to regiospecific steroid hydroxylation.

Structural determinants of progesterone hydroxylation by cytochrome P450 2B5: the role of nonsubstrate recognition site residues.

The highly related rabbit cytochromes P450 2B4 and 2B5 differ in only 12 amino acid positions, but only 2B5 has activity toward progesterone. Previously, simultaneous site-directed mutagenesis of four key substrate recognition site (SRS) residues (114, 294, 363, and 367) was shown to result in interconversion of the androstenedione hydroxylase specificities of cytochrome P450 2B4 and 2B5. However, the progesterone metabolite profiles of the 2B4 quadruple mutant or of a quintuple mutant in which residue 370 was also mutated to the 2B5 residue were not identical to that of P450 2B5. Therefore, single mutants of P450 2B5 at the remaining seven positions were constructed, expressed in Escherichia coli, and studied with progesterone as the substrate. The single mutants at positions 120 and 221, which are outside any known SRS, exhibited a significant alteration in progesterone hydroxylation. Based on these results, Ile-114, Arg-120, Ser-221, Ser-294, Ile-363, and Val-367 in cytochrome P450 2B4 were replaced simultaneously with Phe, His, Pro, Thr, Val, and Ala, respectively, from 2B5. This yielded a mutant with a very similar progesterone metabolite profile to that of 2B5, although the total activity was lower. An additional substitution at residue 370 produced a multiple mutant P450 2B4 I114F-R120H-S221P-S294T-I363V-V367A- T370M with very similar or identical substrate specificity, regio- and stereospecificity and kinetic properties to that of P450 2B5 wild type.

Secobarbital-mediated inactivation of cytochrome P450 2B1 and its active site mutants. Partitioning between heme and protein alkylation and epoxidation.

Secobarbital (SB) is a relatively selective mechanism-based inactivator of cytochrome P450 2B1, that partitions between epoxidation and heme and protein modification during its enzyme inactivation. The SB-2B1 heme adduct formed in situ in a functionally reconstituted system has been spectrally documented and structurally characterized as N-(5-(2-hydroxypropyl)-5-(1-methylbutyl)barbituric acid)protoporphyrin IX. The SB-protein modification has been localized to 2B1 peptide 277-323 corresponding to the active site helix I of cytochrome P450 101. The targeting of heme and this active site peptide suggests that the 2B1 active site topology could influence the course of its inactivation. To explore this possibility, the individual SB epoxidation, heme and protein modification, and corresponding molar partition ratios of the wild type and seven structural 2B1 mutants, site-directed at specific substrate recognition sites, and known to influence 2B1 catalysis were examined after Escherichia coli expression. These studies reveal that Thr-302 is critical for SB-mediated heme N-alkylation, whereas Val-367 is a critical determinant of 2B1 protein modification, and Val-363 is important for SB epoxidation. SB docking into a refined 2B1 homology model coupled with molecular dynamics analyses provide a logical rationale for these findings.

Probing the active site of cytochrome P450 2B1: metabolism of 7-alkoxycoumarins by the wild type and five site-directed mutants.

A series of 7-alkoxycoumarins (chain length of 1-7 carbon atoms) was utilized as active site probes of purified Escherichia coli-expressed cytochrome P450 2B1 wild type and five site-directed mutants (I114V, F206L, V363A, V363L, and G478S). The production of 7-hydroxycoumarin, the O-dealkylation product, by the wild-type enzyme exhibited a rank order of C2 > C4 > C3 > C1 > C5 > C6 = C7. The pattern observed for the P450 I114V mutant was similar to that of the wild-type enzyme, whereas with F206L and G478S mutants, the rate of O-dealkylation was low with all the compounds. In contrast, with V363A, the highest rate of product formation was observed with 7-butoxycoumarin. The V363L mutant preferentially catalyzed the O-dealkylation of 7-methoxy- and 7-ethoxycoumarin, and a further increase in the length of the alkyl chain led to a marked decrease in product formation. The stoichiometry of 7-butoxycoumarin oxidation by V363L suggested that products other than 7-hydroxycoumarin were also formed. HPLC and GC-EIMS analyses revealed that P450 2B1 V363L produced 7-(3-hydroxybutoxy)coumarin and 7-(4-hydroxybutoxy)coumarin as major oxidation products, while the V363A mutant mainly catalyzed the O-dealkylation of 7-butoxycoumarin. Docking of alkoxycoumarins into the active site of a P450 2B1 homology model confirmed the importance of the studied residues in substrate dealkylation and explained the formation of novel 7-butoxycoumarin products by the V363L mutant.