BACKGROUND:The individual character of pharmacokinetics is of great importance in the risk assessment of new drug leads in pharmacological research. Amongst others, it is severely influenced by the properties and inter-individual variability of the enzymes and transporters of the drug detoxification system of the liver. Predicting individual drug biotransformation capacity requires quantitative and detailed models.RESULTS:In this contribution we present the de novo deterministic modeling of atorvastatin biotransformation based on comprehensive published knowledge on involved metabolic and transport pathways as well as physicochemical properties. The model was evaluated on primary human hepatocytes and parameter identifiability analysis was performed under multiple experimental constraints. Dynamic simulations of atorvastatin biotransformation considering the inter-individual variability of the two major involved enzymes CYP3A4 and UGT1A3 based on quantitative protein expression data in a large human liver bank (n = 150) highlighted the variability in the individual biotransformation profiles and therefore also points to the individuality of pharmacokinetics.CONCLUSIONS:A dynamic model for the biotransformation of atorvastatin has been developed using quantitative metabolite measurements in primary human hepatocytes. The model comprises kinetics for transport processes and metabolic enzymes as well as population liver expression data allowing us to assess the impact of inter-individual variability of concentrations of key proteins. Application of computational tools for parameter sensitivity analysis enabled us to considerably improve the validity of the model and to create a consistent framework for precise computer-aided simulations in toxicology.

Atorvastatin δ-lactone, a major, pharmacologically inactive metabolite, has been associated with toxicity. In a previous study we showed that polymorphisms of UGT1A3 influence atorvastatin δ-lactone formation. Here we investigated the reverse reaction, atorvastatin δ-lactone hydrolysis, in a human liver bank. Screening of microarray data revealed paraoxonases PON1 and PON3 among 17 candidate esterases. Microsomal δ-lactone hydrolysis was significantly correlated to PON1 and PON3 protein (r(s) = 0.60; r(s) = 0.62, respectively; P < 0.0001). PON1 and PON3 were strongly correlated to each other (r(s) = 0.60) but PON1 was shown to be more extensively glycosylated than PON3. In addition a novel splice-variant of PON3 was identified. Genotyping of 40 polymorphisms within the PON-locus identified PON1 promoter polymorphisms (-108T > C, -832G > A, -1741G > A) and a tightly linked group of PON3 polymorphisms (-4984A > G, -4105G > A, -1091A > G, -746C > T, and F21F) to be associated with changes in atorvastatin δ-lactone hydrolysis and expression of PON1 but not PON3. However, carriers of the common PON1 polymorphisms L55M or Q192R showed no difference in δ-lactone hydrolysis or PON expression. Haplotype analysis revealed decreased δ-lactone hydrolysis in carriers of the most common haplotype *1 compared to carriers of haplotypes *2, *3, *4, and *7. Analysis of non-genetic factors showed association of hepatocellular and cholangiocellular carcinoma with decreased PON1 and PON3 expression, respectively. Increased C-reactive protein and γ-glutamyl transferase levels were associated with decreased protein expression of both enzymes, and increased bilirubin levels, cholestasis, and presurgical exposure to omeprazole or pantoprazole were related to decreased PON3 protein. In conclusion, PON-locus polymorphisms affect PON1 expression whereas non-genetic factors have an effect on PON1 and PON3 expression. This may influence response to therapy or adverse events in statin treatment.

Human hepatocytes in primary culture are a very useful model to directly assess induction of gene expression by xenobiotics. We developed a cytochrome P450 (P450) activity cocktail assay using model substrates for the seven important P450s 1A2 (phenacetin), 2B6 (bupropion), 2C8 (amodiaquine), 2C9 (tolbutamide), 2C19 (S-mephenytoin), 2D6 (propafenone), and 3A4 (atorvastatin). Metabolite formation was determined by liquid chromatography-tandem mass spectrometry in hepatocyte culture supernatants. Atorvastatin has not been previously assessed as a CYP3A probe drug. We demonstrate highly selective atorvastatin ortho-hydroxylation by CYP3A4 by recombinant P450s. In human liver microsomes ortho-hydroxyatorvastatin formation was highly correlated with CYP3A4 protein content (r(s) = 0.78, p < 0.0001, n = 150). We profiled induction of these P450 activities in primary human hepatocytes after treatment with 30 microM atorvastatin, lovastatin, pravastatin, rosuvastatin, and simvastatin for 24 to 72 h. Except for pravastatin, all statins induced P450 activities to various degrees, approximately in the order atorvastatin > simvastatin > lovastatin > rosuvastatin. Inducibility of P450s followed the order CYP2C8 > CYP3A4 > CYP2C9 > CYP2B6 > CYP2C19 approximately CYP2D6 > CYP1A2. The strongest induction was observed for amodiaquine N-desalkylation, which was induced approximately 20-fold. Quantitative reverse transcription-polymerase chain reaction confirmed corresponding changes on the mRNA level with even more dramatic induction up to almost 100-fold. These data suggest a broader inducing effect of statins on cytochrome P450s and possibly other absorption, distribution, metabolism, and excretion genes than previously known, thus further emphasizing their drug-drug interaction potential. Our cocktail assay should be helpful for economical use of human hepatocytes in the assessment of P450 induction by drugs and drug candidates.

In mammals, Cytochrome P450 (CYP) enzymes are bound to membranes of the endoplasmic reticulum and mitochondria, where they are responsible for the oxidative metabolism of many xenobiotics as well as organic endogenous compounds. In humans, 57 isoforms were identified which are classified based on sequence homology. In the present work, we demonstrate the performance of a mass spectrometry-based strategy to simultaneously detect and differentiate distinct human Cytochrome P450 (CYP) isoforms including the highly similar CYP3A4, CYP3A5, CYP3A7, as well as CYP2C8, CYP2C9, CYP2C18, CYP2C19, and CYP4F2, CYP4F3, CYP4F11, CYP4F12. Compared to commonly used immunodetection methods, mass spectrometry overcomes limitations such as low antibody specificity and offers high multiplexing possibilities. Furthermore, CYP phosphorylation, which may affect various biochemical and enzymatic properties of these enzymes, is still poorly analyzed, especially in human tissues. Using titanium dioxide resin combined with tandem mass spectrometry for phosphopeptide enrichment and sequencing, we discovered eight human P450 phosphorylation sites, seven of which were novel. The data from surgical human liver samples establish that the isoforms CYP1A2, CYP2A6, CYP2B6, CYP2E1, CYP2C8, CYP2D6, CYP3A4, CYP3A7, and CYP8B1 are phosphorylated in vivo. These results will aid in further investigation of the functional significance of protein phosphorylation for this important group of enzymes.