Role of Positively Charged Residues Lys267, Lys270, and Arg411 of Cytochrome P450scc (Cyp11A1) in Interaction with Adrenodoxin

Cytochrome P450scc and adrenodoxin are redox proteins of the electron transfer chain of the inner mitochondrial membrane steroid hydroxylases. In the present work site-directed mutagenesis of the charged residues of cytochrome P450scc and adrenodoxin, which might be involved in interaction, was used to study the nature of electrostatic contacts between the hemeprotein and the ferredoxin. The target residues for mutagenesis were selected based on the theoretical model of cytochrome P450scc-adrenodoxin complex and previously reported chemical modification studies of cytochrome P450scc. In the present work, to clarify the molecular mechanism of hemeprotein interaction with ferredoxin, we constructed cytochrome P450scc Lys267, Lys270, and Arg411 mutants and Glu47 mutant of adrenodoxin and analyzed their possible role in electrostatic interaction and the role of these residues in the functional activity of the proteins. Charge neutralization at positions Lys267 or Lys270 of cytochrome P450scc causes no significant effect on the physicochemical and functional properties of cytochrome P450scc. However, cytochrome P450scc mutant Arg411Gln was found to exhibit decreased binding affinity to adrenodoxin and lower activity in the cholesterol side chain cleavage reaction. Studies of the functional properties of Glu47Gln and Glu47Arg adrenodoxin mutants indicate that a negatively charged residue in the loop covering the Fe2S2 cluster, being important for maintenance of the correct architecture of these structural elements of ferredoxin, is not directly involved in electrostatic interaction with cytochrome P450scc. Moreover, our results indicate the presence of at least two different binding (contact) sites on the proximal surface of cytochrome P450scc with different electrostatic input to interaction with adrenodoxin. In the binary complex, the positively charged sites of the proximal surface of cytochrome P450scc well correspond to the two negatively charged sites of adrenodoxin: the "interaction" domain site and the "core" domain site.     

Immobilization and activity assay of cytochrome P450 on patterned lipid membranes

We report on a methodology for immobilizing cytochrome P450 on the surface of micropatterned lipid bilayer membranes and measuring the enzymatic activity. The patterned bilayer comprised a matrix of polymeric lipid bilayers and embedded fluid lipid bilayers. The polymeric lipid bilayer domains act as a barrier to confine fluid lipid bilayers in defined areas and as a framework to stabilize embedded membranes. The fluid bilayer domains, on the other hand, can contain lipid compositions that facilitate the fusion between lipid membranes, and are intended to be used as the binding agent of microsomes containing rat CYP1A1. By optimizing the membrane compositions of the fluid bilayers, we could selectively immobilize microsomal membranes on these domains. The enzymatic activity was significantly higher on lipid bilayer substrates compared with direct adsorption on glass. Furthermore, competitive assay experiment between two fluorogenic substrates demonstrated the feasibility of bioassays based on immobilized P450s.     

Homology modeling of the three membrane proteins of the dhurrin metabolon: catalytic sites, membrane surface association and protein-protein interactions

Formation of metabolons (macromolecular enzyme complexes) facilitates the channelling of substrates in biosynthetic pathways. Metabolon formation is a dynamic process in which transient structures mediated by weak protein-protein interactions are formed. In Sorghum, the cyanogenic glucoside dhurrin is derived from l-tyrosine in a pathway involving the two cytochromes P450 (CYPs) CYP79A1 and CYP71E1, a glucosyltransferase (UGT85B1), and the redox partner NADPH-dependent cytochrome P450 reductase (CPR). Experimental evidence suggests that the enzymes of this pathway form a metabolon. Homology modeling of the three membrane bound proteins was carried out using the Sybyl software and available relevant crystal structures. Residues involved in tight positioning of the substrates and intermediates in the active sites of CYP79A1 and CYP71E1 were identified. In both CYPs, hydrophobic surface domains close to the N-terminal trans-membrane anchor and between the F′ and G helices were identified as involved in membrane anchoring. The proximal surface of both CYPs showed positively charged patches complementary to a negatively charged bulge on CPR carrying the FMN domain. A patch of surface exposed, positively charged amino acid residues positioned on the opposite face of the membrane anchor was identified in CYP71E1 and might be involved in binding UGT85B1 via a hypervariable negatively charged loop in this protein.     

Interactions of mammalian cytochrome P450, NADPH-cytochrome P450 reductase, and cytochrome b(5) enzymes

An immobilized system was developed to detect interactions of human cytochromes P450 (P450) with the accessory proteins NADPH-P450 reductase and cytochrome b(5) (b(5)) using an enzyme-linked affinity approach. Purified enzymes were first bound to wells of a polystyrene plate, and biotinylated partner enzymes were added and bound. A streptavidin-peroxidase complex was added, and protein-protein binding was monitored by measuring peroxidase activity of the bound biotinylated proteins. In a model study, we examined protein-protein interactions of Pseudomonas putida putidaredoxin (Pdx) and putidaredoxin reductase (PdR). A linear relationship (r(2)=0.96) was observed for binding of PdR-biotin to immobilized Pdx compared with binding of Pdx-biotin to immobilized PdR (the estimated K(d) value for the Pdx.PdR complex was 0.054muM). Human P450 2A6 interacted strongly with NADPH-P450 reductase; the K(d) values (with the reductase) ranged between 0.005 and 0.1muM for P450s 2C19, 2D6, and 3A4. Relatively weak interaction was found between holo-b(5) or apo-b(5) (devoid of heme) with NADPH-P450 reductase. Among the rat, rabbit, and human P450 1A2 enzymes, the rat enzyme showed the tightest interaction with b(5), although no increases in 7-ethoxyresorufin O-deethylation activities were observed with any of the P450 1A2 enzymes. Human P450s 2A6, 2D6, 2E1, and 3A4 interacted well with b(5), with P450 3A4 yielding the lowest K(d) values followed by P450s 2A6 and 2D6. No appreciable increases in interaction between human P450s with b(5) or NADPH-P450 reductase were observed when typical substrates for the P450s were included. We also found that NADPH-P450 reductase did not cause changes in the P450.substrate K(d) values estimated from substrate-induced UV-visible spectral changes with rabbit P450 1A2 or human P450 2A6, 2D6, or 3A4. Collectively, the results show direct and tight interactions between P450 enzymes and the accessory proteins NADPH-P450 reductase and b(5), with different affinities, and that ligand binding to mammalian P450s did not lead to increased interaction between P450s and the reductase.     

Comparative study of topogenesis of cytochrome P450scc (CYP11A1) and its hybrids with adrenodoxin expressed in <Emphasis Type="Italic">Escherichia coli</Emphasis> cells

Hybrid proteins consisting of the mature form of cytochrome P450scc (mP) and adrenodoxin (Ad), attached to either the NH2- or COOH-terminus (Ad-mP and mP-Ad, respectively), were expressed in E. coli. Spectral and catalytic properties of P450scc were studied using the membrane fraction of E. coli cells. It has been shown that the Ad amino acid sequence attached to the termini of the P450scc-domain neither affects the insertion of a hybrid protein into the cytoplasmic membrane nor influences its heme binding ability. The results suggest that Ad attached to the NH2-terminus does not markedly affect the folding of the P450scc-domain, but cholesterol hydroxylase/lyase activity of the Ad-mP hybrid was found to be much lower than that of the native P450scc enzyme. The modification of the COOH-terminus does not alter the specific P450scc activity, but results in a dramatic increase in the amount of hybrid protein with incorrectly folded P450scc domain.     

Involvement of cytochrome P450 3A in the metabolism and covalent binding of 14C-monocrotaline in rat liver microsomes

The metabolism and covalent binding of ¹⁴C-monocrotaline in Sprague–Dawley (SD) rat liver microsomes was investigated using the inducers dexamethasone, clotrimazole, pregnenolone-16α-carbonitrile, and phenobarbital. Monocrotaline is a pyrrolizidine alkaloid (PA) that causes a syndrome in rats that is a model for human primary pulmonary hypertension. It has been documented that bioactivation of PAs (dehydrogenation to reactive pyrroles) in the liver by cytochromes P450 is required for their toxicity. Covalent binding of these reactive pyrroles to tissue macromolecules has been hypothesized to correspond to PA toxicosis. We correlated metabolism and total microsomal covalent binding of ¹⁴C-monocrotaline with cytochrome P450 3A using the aforementioned inducers, troleandomycin (a cytochrome P450 3A inhibitor), erythromycin N-demethylase assay of cytochrome P450 3A activity, and Western blots employing anti-rat cytochrome P450 3A antibodies. In addition, autoradiography of membranes electroblotted from SDS-PAGE demonstrated the formation of radiolabeled adducts with specific protein(s). The most intensely radiolabeled protein bands have an apparent molecular weight of ∼52 kDa, which was similar to the molecular weight detected by anti-rat cytochrome P450 3A antibodies in the Western blots. No radiolabeled proteins were detected in microsomes pretreated with troleandomycin. © 1998 John Wiley & Sons, Inc. J Biochem Toxicol 12: 157–166, 1998     

cytochrome p450 superfamily in Arabidopsis thaliana: isolation of cDNAs,Differential Expression,and RFLP mapping of multiple cytochromes P450

We have isolated multiple cDNAs encoding cytochromes P450 (P450s) from Arabidopsis thaliana employing a PCR strategy. Degenerate oligonucleotide primers were designed from amino acid sequences conserved between two plant P450s, CYP71A1 and CYP73A2, including the heme-binding site and the proline-rich motif found in the N-terminal region, and 11 putative P450 fragments were amplified from first-strand cDNA from 7-day-old Arabidopsis as a template. With these PCR fragments as hybridization probes, 13 full-length and 3 partial cDNAs encoding different P450s have been isolated from an Arabidopsis cDNA library. These P450s have been assigned to either one of the established subfamilies: CYP71B, CYP73A, and CYP83A; or novel subfamilies: CYP76C, CYP83B, and CYP91A. The primary protein structures predicted from the cDNA sequences revealed that the regions around both the heme-binding site and the proline-rich motif were highly conserved among all these P450s. The N-terminal structures of the predicted P450 proteins suggested that these Arabidopsis P450s were located at the endoplasmic reticulum membrane. The loci of four P450 genes were determined by RFLP mapping. One of the clones, CYP71B2, was located at a position very close to the ga4 and gai mutations. RNA blot analysis showed expression patterns unique to each of the P450s in terms of tissue specificity and responsiveness to wounding and light/dark cycle, implicating involvement of these P450s in diverse metabolic processes.     

Site-Directed Mutagenesis of Cytochrome b5 for Studies of Its Interaction with Cytochrome P450

We have shown earlier that microsomal cytochrome b ₅ can form a specific complex with mitochondrial cytochrome P450 (cytochrome P450scc). The formation of the complex between these two heme proteins was proved spectrophotometrically, by affinity chromatography on immobilized cytochrome b ₅, and by measuring the cholesterol side-chain cleavage activity of cytochrome P450scc in a reconstituted system in the presence of cytochrome b ₅. To further study the mechanism of interaction of these heme proteins and evaluate the role of negatively charged amino acid residues Glu42, Glu48, and Asp65 of cytochrome b ₅, which are located at the site responsible for interaction with electron transfer partners, we used sitedirected mutagenesis to replace residues Glu42 and Glu48 with lysine and residue Asp65 with alanine. The resulting mutant forms of cytochrome b ₅ were expressed in E. coli, and full-length and truncated forms (shortened from the C-terminal sequence due to cleavage of 40 amino acid residues) of these cytochrome b ₅ mutants were purified. Addition of the truncated forms of cytochrome b ₅ (which do not contain the hydrophobic C-terminal sequence responsible for interaction with the membrane) to the reconstituted system containing cytochrome P450scc caused practically no stimulation of catalytic activity, indicating an important role of the hydrophobic fragment of cytochrome b ₅ in its interaction with cytochrome P450scc. However, full-length cytochrome b ₅ and the full-length Glu48Lys and Asp65Ala mutant forms of cytochrome b ₅ stimulated the cholesterol side-chain cleavage reaction catalyzed by cytochrome P450scc by 100%, suggesting that residues Glu48 and Asp65 of cytochrome b ₅ are not directly involved in its interaction with cytochrome P450scc. The replacement of Glu42 for lysine, however, made the Glu42Lys mutant form of cytochrome b ₅ about 40% less effective in stimulation of the cholesterol side-chain cleavage activity of cytochrome P450scc, indicating that residue Glu42 of cytochrome b ₅ is involved in electrostatic interactions with cytochrome P450scc. Residues Glu42 and Glu48 of cytochrome b ₅ appear to participate in electrostatic interaction with microsomal type cytochrome P450.     

Site-Directed Mutagenesis of Cytochrome P450scc. II. Effect of Replacement of the Arg425 and Arg426 Residues on the Structural and Functional Properties of the Cytochrome P450scc

Cytochrome P450-dependent monooxygenases, in spite of their wide distribution, can be simply divided into a few groups differing in the location of the electron transfer chain and their composition. The two main groups of cytochrome P450-dependent monooxygenases are the mitochondrial and microsomal enzymes. While in two-component microsomal cytochrome P450-dependent monooxygenases electrons are supplied to cytochrome P450 by a flavoprotein (NADPH-cytochrome P450 reductase), in three-component mitochondrial monooxygenases the electrons are supplied to cytochrome P450 by a low molecular weight protein (ferredoxin). The interaction of cytochrome P450 with NADPH-cytochrome P450 reductase and ferredoxin is the subject of intensive studies. Using chemical modification, chemical cross-linking, and sitedirected mutagenesis, we identified surface exposed positively charged residues of cytochrome P450scc which might be important for interaction with adrenodoxin. Theoretical analysis of the distribution of surface electrostatic potential in cytochrome P450 indicates that in contrast to microsomal monooxygenases, cytochromes P450 of mitochondrial type, and cholesterol side-chain cleavage cytochrome P450 (P450scc) in part, carry on the proximal surface an evidently positively charged site that is formed by residues Arg425 and Arg426. In the present work, to estimate the functional role of Arg425 and Arg426 of cytochrome P450scc, we used site-directed mutagenesis to replace these residues with glutamine. The results indicate that residues Arg425 and Arg426 are involved in the formation of a heme-binding center and electrostatic interaction of cytochrome P450scc with its physiological electron-transfer partner, adrenodoxin.     

Structural characterization of human cholesterol 7α-hydroxylase

Hepatic conversion to bile acids is a major elimination route for cholesterol in mammals. CYP7A1 catalyzes the first and rate-limiting step in classic bile acid biosynthesis, converting cholesterol to 7α-hydroxycholesterol. To identify the structural determinants that govern the stereospecific hydroxylation of cholesterol, we solved the crystal structure of CYP7A1 in the ligand-free state. The structure-based mutation T104L in the B′ helix, corresponding to the nonpolar residue of CYP7B1, was used to obtain crystals of complexes with cholest-4-en-3-one and with cholesterol oxidation product 7-ketocholesterol (7KCh). The structures reveal a motif of residues that promote cholest-4-en-3-one binding parallel to the heme, thus positioning the C7 atom for hydroxylation. Additional regions of the binding cavity (most distant from the access channel) are involved to accommodate the elongated conformation of the aliphatic side chain. Structural complex with 7KCh shows an active site rigidity and provides an explanation for its inhibitory effect. Based on our previously published data, we proposed a model of cholesterol abstraction from the membrane by CYP7A1 for metabolism. CYP7A1 structural data provide a molecular basis for understanding of the diversity of 7α-hydroxylases, on the one hand, and cholesterol-metabolizing enzymes adapted for their specific activity, on the other hand.     

Substrate-modulated Cytochrome P450 17A1 and Cytochrome b5 Interactions Revealed by NMR

The membrane heme protein cytochrome b₅ (b₅) can enhance, inhibit, or have no effect on cytochrome P450 (P450) catalysis, depending on the specific P450, substrate, and reaction conditions, but the structural basis remains unclear. Here the interactions between the soluble domain of microsomal b₅ and the catalytic domain of the bifunctional steroidogenic cytochrome P450 17A1 (CYP17A1) were investigated. CYP17A1 performs both steroid hydroxylation, which is unaffected by b₅, and an androgen-forming lyase reaction that is facilitated 10-fold by b₅. NMR chemical shift mapping of b₅ titrations with CYP17A1 indicates that the interaction occurs in an intermediate exchange regime and identifies charged surface residues involved in the protein/protein interface. The role of these residues is confirmed by disruption of the complex upon mutagenesis of either the anionic b₅ residues (Glu-48 or Glu-49) or the corresponding cationic CYP17A1 residues (Arg-347, Arg-358, or Arg-449). Cytochrome b₅ binding to CYP17A1 is also mutually exclusive with binding of NADPH-cytochrome P450 reductase. To probe the differential effects of b₅ on the two CYP17A1-mediated reactions and, thus, communication between the superficial b₅ binding site and the buried CYP17A1 active site, CYP17A1/b₅ complex formation was characterized with either hydroxylase or lyase substrates bound to CYP17A1. Significantly, the CYP17A1/b₅ interaction is stronger when the hydroxylase substrate pregnenolone is present in the CYP17A1 active site than when the lyase substrate 17α-hydroxypregnenolone is in the active site. These findings form the basis for a clearer understanding of this important interaction by directly measuring the reversible binding of the two proteins, providing evidence of communication between the CYP17A1 active site and the superficial proximal b₅ binding site.     

Tomato cytochrome P450 CYP734A7 functions in brassinosteroid catabolism

Several cytochrome P450 monooxygenases (P450s) catalyze essential oxidative reactions in brassinosteroid (BR) biosynthesis as well as in BR catabolism; however, only limited information exists on the P450s involved in the BR catabolic pathway. Here, we report the characterization of two P450 mRNAs, CYP734A7 and CYP734A8, from Lycopersicon esculentum. These P450s show high homology with Arabidopsis CYP734A1/BAS1 (formerly CYP72B1), which inactivates BRs via C-26 hydroxylation. Transgenic tobacco plants that constitutively overexpressed CYP734A7 showed an extreme dwarf phenotype similar to BR deficiency. Quantitative gas chromatography-mass spectrometry analysis of endogenous BRs in the transgenic plants showed that the levels of castasterone and 6-deoxocastasterone significantly decreased in comparison with those in wild-type plants. By measuring the Type I substrate-binding spectra using recombinant CYP734A7, the dissociation constants for castasterone, brassinolide, and 6-deoxocastasterone were determined to be 6.7, 12, and 12 microM, respectively. In an in vitro assay, CYP734A7 was confirmed to metabolize castasterone to 26-hydroxycastasterone. In addition, 28-norcastasterone and brassinolide were converted to the hydroxylated products. The expression of CYP734A7 and CYP734A8 genes in tomato seedlings was upregulated by exogenous application of bioactive BRs. These results indicated that CYP734A7 is a C-26 hydroxylase of BRs and is likely involved in BR catabolism in tomato. The presence of the CYP734A subfamily in various plant species suggests that oxidative inactivation of BRs by these proteins is a widespread phenomenon in plants.     

Molecular definitions of fatty acid hydroxylases in Arabidopsis thaliana

Towards defining the function of Arabidopsis thaliana fatty acid hydroxylases, five members of the CYP86A subfamily have been heterologously expressed in baculovirus-infected Sf9 cells and tested for their ability to bind a range of fatty acids including unsubstituted (lauric acid (C12:0) and oleic acid (C18:1)) and oxygenated (9,10-epoxystearic acid and 9,10-dihydroxystearic acid). Comparison between these five P450s at constant P450 content over a range of concentrations for individual fatty acids indicates that binding of different fatty acids to CYP86A2 always results in a higher proportion of high spin state heme than binding titrations conducted with CYP86A1 or CYP86A4. In comparison to these three, CYP86A7 and CYP86A8 produce extremely low proportions of high spin state heme even with the most effectively bound fatty acids. In addition to their previously demonstrated lauric acid hydroxylase activities, all CYP86A proteins are capable of hydroxylating oleic acid but not oxygenated 9,10-epoxystearic acid. Homology models have been built for these five enzymes that metabolize unsubstituted fatty acids and sometimes bind oxygenated fatty acids. Comparison of the substrate binding modes and predicted substrate access channels indicate that all use channel pw2a consistent with the crystal structures and models of other fatty acid-metabolizing P450s in bacteria and mammals. Among these P450s, those that bind internally oxygenated fatty acids contain polar residues in their substrate binding cavity that help stabilize these charged/polar groups within their largely hydrophobic catalytic site.     

Site-Directed Mutagenesis of Cytochrome P450scc (CYP11A1). Effect of Lysine Residue Substitution on Its Structural and Functional Properties

Our previous chemical modification and cross-linking studies identified some positively charged amino acid residues of cytochrome P450scc that may be important for its interaction with adrenodoxin and for its functional activity. The present study was undertaken to further evaluate the role of these residues in the interaction of cytochrome P450scc with adrenodoxin using site-directed mutagenesis. Six cytochrome P450scc mutants containing replacements of the surface-exposed positively charged residues (Lys103Gln, Lys110Gln, Lys145Gln, Lys394Gln, Lys403Gln, and Lys405Gln) were expressed in E. coli cells, purified as a substrate-bound high-spin form, and characterized as compared to the wild-type protein. The replacement of the surface Lys residues does not dramatically change the protein folding or the heme pocket environment as judged from limited proteolysis and spectral studies of the cytochrome P450 mutants. The replacement of Lys in the N-terminal sequence of P450scc does not dramatically affect the activity of the heme protein. However, mutant Lys405Gln revealed rather dramatic loss of cholesterol side-chain cleavage activity, efficiency of enzymatic reduction in a reconstituted system, and apparent dissociation constant for adrenodoxin binding. The present results, together with previous findings, suggest that the changes in functional activity of mutant Lys405Gln may reflect the direct participation of this amino acid residue in the electrostatic interaction of cytochrome P450scc with its physiological partner, adrenodoxin.     

Functional reconstitution of monomeric CYP3A4 with multiple cytochrome P450 reductase molecules in Nanodiscs

Traditional reconstitution of membrane cytochromes P450 monooxygenase system requires efficient solubilization of both P450 heme enzymes and redox partner NADPH dependent reductase, CPR, either in mixed micellar solution or by incorporation in liposomes. Here we describe a simple alternative approach to assembly of soluble complexes of monomeric human hepatic cytochrome P450 CYP3A4 with CPR by co-incorporation into nanoscale POPC bilayer Nanodiscs. Stable and fully functional complexes with different CPR:CYP3A4 stoichiometric ratios are formed within several minutes after addition of the full-length CPR to the solution of CYP3A4 preassembled into POPC Nanodiscs at 37 °C. We find that the steady state rates of NADPH oxidation and testosterone hydroxylation strongly depend on CPR:CYP3A4 ratio and reach maximum at tenfold molar access of CPR. The binding of CPR to CYP3A4 in Nanodiscs is tight, such that complexes with different stoichiometry can be separated by size-exclusion chromatography. Reconstitution systems based on the co-incorporation of CPR into preformed Nanodiscs with different human cytochromes P450 are suitable for high-throughput screening of substrates and inhibitors and for drug-drug interaction studies.     

Establishment of a yeast system that stably expresses human cytochrome P450 reductase: application for the study of drug metabolism of cytochrome P450s in vitro

Cytochrome P450s (CYPs) hold a balance in studying pharmacokinetics, toxico-kinetics, drug metabolism, and drug-drug interactions, which require association with cytochrome P450 reductase (CPR) to achieve optimal activity. A novel system of Saccharomyces cerevisiae useful for expression studies of mammalian microsomal CYPs was established. Human CPR (hCPR) was co-expressed with human CYP3A4 (hCYP3A4) in this system, and two expression plasmids pTpLC and pYeplac195-3A4 containing the cDNA of hCPR and hCYP3A4 were constructed, respectively. The two plasmids were applied first and controlled by phosphoglycerate kinase (PGK) promoter. S. cerevisiae BWG1-7alpha transformed with the expression plasmids produced the respective proteins in the expected molecular sizes reactive with both anti-hCYP3A4 immunoglobulin (Ig) and anti-hCPR Ig. The activity of hCPR in yeast BWG-CPR was 443.2 nmol reduced cytochrome c/min/mg, which was about three times the CPR activity of the microsome prepared from the parental yeast. The protein amount of hCYP3A4 in BWG-CPR/3A4 was 35.53 pmol/mg, and the 6beta-hydroxylation testosterone formation activity of hCYP3A4 expressed was 7.5 nmol/min/nmol CYP, 30 times higher than the activity of hCYP3A4 expressed in the parental yeast, and almost two times the activity of hCYP3A4 from homologous human liver microsome. Meanwhile, BWG-CPR/3A4 retained 100 generations under nonselective culture conditions, indicating this yeast was a mitotically stable transformant. BWG-CPR was further tested daily by the PCR amplification of hCPR of yeast genome, Western blot analysis, and the activity assay of hCPR of yeast microsome. This special expression host for CYPs was validated to be stable and efficient for the expression of CYPs, applying as an effective selection model for the drug metabolism in vitro.