Development of a Novel In Silico Classification Model to Assess Reactive Metabolite Formation in the Cysteine Trapping Assay and Investigation of Important Substructures.

Predicting whether a compound can cause drug-induced liver injury (DILI) is difficult due to the complexity of drug mechanism. The cysteine trapping assay is a method for detecting reactive metabolites that bind to microsomes covalently. However, it is cumbersome to use 35S isotope-labeled cysteine for this assay. Therefore, we constructed an in silico classification model for predicting a positive/negative outcome in the cysteine trapping assay. We collected 475 compounds (436 in-house compounds and 39 publicly available drugs) based on experimental data performed in this study, and the composition of the results showed 248 positives and 227 negatives. Using a Message Passing Neural Network (MPNN) and Random Forest (RF) with extended connectivity fingerprint (ECFP) 4, we built machine learning models to predict the covalent binding risk of compounds. In the time-split dataset, AUC-ROC of MPNN and RF were 0.625 and 0.559 in the hold-out test, restrictively. This result suggests that the MPNN model has a higher predictivity than RF in the time-split dataset. Hence, we conclude that the in silico MPNN classification model for the cysteine trapping assay has a better predictive power. Furthermore, most of the substructures that contributed positively to the cysteine trapping assay were consistent with previous results.

Development of Novel Methods for QSAR Modeling by Machine Learning Repeatedly: A Case Study on Drug Distribution to Each Tissue.

Artificial intelligence is expected to help identify excellent candidates in drug discovery. However, we face a lack of data, as it is time-consuming and expensive to acquire raw data perfectly for many compounds. Hence, we tried to develop a novel quantitative structure-activity relationship (QSAR) method to predict a parameter more precisely from an incomplete data set via optimizing data handling by making use of predicted explanatory variables. As a case study we focused on the tissue-to-plasma partition coefficient (Kp), which is an important parameter for understanding drug distribution in tissues and building the physiologically based pharmacokinetic model and is a representative of small and sparse data sets. In this study, we predicted the Kp values of 119 compounds in nine tissues (adipose, brain, gut, heart, kidney, liver, lung, muscle, and skin), although some of these were not available. To fill the missing values in Kp for each tissue, first we predicted those Kp values by the nonmissing data set using a random forest (RF) model with in vitro parameters (log P, fu, Drug Class, and fi) like a classical prediction by a QSAR model. Next, to predict the tissue-specific Kp values in a test data set, we constructed a second RF model with not only in vitro parameters but also the Kp values of other tissues (i.e., other than target tissues) predicted by the first RF model as explanatory variables. Furthermore, we tested all possible combinations of explanatory variables and selected the model with the highest predictability from the test data set as the final model. The evaluation of Kp prediction accuracy based on the root-mean-square error and R2 value revealed that the proposed models outperformed other machine learning methods such as the conventional RF and message-passing neural networks. Significant improvements were observed in the Kp values of adipose tissue, brain, kidney, liver, and skin. These improvements indicated that the Kp information on other tissues can be used to predict the same for a specific tissue. Additionally, we found a novel relationship between each tissue by evaluating all combinations of explanatory variables. In conclusion, we developed a novel RF model to predict Kp values. We hope that this method will be applied to various problems in the field of experimental biology which often contains missing values in the near future.

On the difficulty of validating molecular generative models realistically: a case study on public and proprietary data.

While a multitude of deep generative models have recently emerged there exists no best practice for their practically relevant validation. On the one hand, novel de novo-generated molecules cannot be refuted by retrospective validation (so that this type of validation is biased); but on the other hand prospective validation is expensive and then often biased by the human selection process. In this case study, we frame retrospective validation as the ability to mimic human drug design, by answering the following question: Can a generative model trained on early-stage project compounds generate middle/late-stage compounds de novo? To this end, we used experimental data that contains the elapsed time of a synthetic expansion following hit identification from five public (where the time series was pre-processed to better reflect realistic synthetic expansions) and six in-house project datasets, and used REINVENT as a widely adopted RNN-based generative model. After splitting the dataset and training REINVENT on early-stage compounds, we found that rediscovery of middle/late-stage compounds was much higher in public projects (at 1.60%, 0.64%, and 0.21% of the top 100, 500, and 5000 scored generated compounds) than in in-house projects (where the values were 0.00%, 0.03%, and 0.04%, respectively). Similarly, average single nearest neighbour similarity between early- and middle/late-stage compounds in public projects was higher between active compounds than inactive compounds; however, for in-house projects the converse was true, which makes rediscovery (if so desired) more difficult. We hence show that the generative model recovers very few middle/late-stage compounds from real-world drug discovery projects, highlighting the fundamental difference between purely algorithmic design and drug discovery as a real-world process. Evaluating de novo compound design approaches appears, based on the current study, difficult or even impossible to do retrospectively.Scientific Contribution This contribution hence illustrates aspects of evaluating the performance of generative models in a real-world setting which have not been extensively described previously and which hopefully contribute to their further future development.

Development of a 2D-QSAR Model for Tissue-to-Plasma Partition Coefficient Value with High Accuracy Using Machine Learning Method, Minimum Required Experimental Values, and Physicochemical Descriptors.

BACKGROUND: The demand for physiologically based pharmacokinetic (PBPK) model is increasing currently. New drug application (NDA) of many compounds is submitted with PBPK models for efficient drug development. Tissue-to-plasma partition coefficient (Kp) is a key parameter for the PBPK model to describe differential equations. However, it is difficult to obtain the Kp value experimentally because the measurement of drug concentration in the tissue is much harder than that in plasma. OBJECTIVE: Instead of experiments, many researchers have sought in silico methods. Today, most of the models for Kp prediction are using in vitro and in vivo parameters as explanatory variables. We thought of physicochemical descriptors that could improve the predictability. Therefore, we aimed to develop the two-dimensional quantitative structure-activity relationship (2D-QSAR) model for Kp using physicochemical descriptors instead of in vivo experimental data as explanatory variables. METHODS: We compared our model with the conventional models using 20-fold cross-validation according to the published method (Yun et al. J Pharmacokinet Pharmacodyn 41:1-14, 2014). We used random forest algorithm, which is known to be one of the best predictors for the 2D-QSAR model. Finally, we combined minimum in vitro experimental values and physiochemical descriptors. Thus, the prediction method for Kp value using a few in vitro parameters and physicochemical descriptors was developed; this is a multimodal model. RESULTS: Its accuracy was found to be superior to that of the conventional models. Results of this research suggest that multimodality is useful for the 2D-QSAR model [RMSE and % of two-fold error: 0.66 and 42.2% (Berezohkovsky), 0.52 and 52.2% (Rodgers), 0.65 and 34.6% (Schmitt), 0.44 and 61.1% (published model), 0.41 and 62.1% (traditional model), 0.39 and 64.5% (multimodal model)]. CONCLUSION: We could develop a 2D-QSAR model for Kp value with the highest accuracy using a few in vitro experimental data and physicochemical descriptors.

Prediction of Compound Plasma Concentration-Time Profiles in Mice Using Random Forest.

Pharmacokinetic (PK) parameters such as clearance (CL) and volume of distribution (Vd) have been the subject of previous in silico predictive models. However, having information of the concentration over time profile explicitly can provide additional value like time above MIC or AUC, etc., to understand both the efficacy and safety-related aspects of a compound. In this work, we developed machine learning models for plasma concentration-time profiles after both i.v. and p.o. dosing for a series of 17 in-house projects. For explanatory variables, MACCS Keys chemical descriptors as well as in silico and experimental in vitro PK parameters were used. The predictive accuracy of random forest (RF), message passing neural network, 2-compartment models using estimated CL and Vdss, and an average model (as a control experiment) was investigated using 5-fold cross-validation (5-fold CV) and leave-one-project-out validation (LOPO-V). The predictive accuracy of RF in 5-fold CV for i.v. and p.o. plasma concentration-time profiles was the best among the models studied, with an RMSE for i.v. dosing at 0.08, 1, and 8 h of 0.245, 0.474, and 0.462, respectively, and an RMSE for p.o. dosing at 0.25, 1, and 8 h of 0.500, 0.612, and 0.509, respectively. Furthermore, by investigating the importance of the in vitro PK parameters using the Gini index, we observed that the general prior knowledge in ADME research was reflected well in the respective feature importance of in vitro parameters such as predicted human Vd (hVd) for the initial distribution, mouse intrinsic CL and unbound fraction of mouse plasma for the elimination process, and Caco2 permeability for the absorption process. Also, this model is the first model that can predict twin peaks in the concentration-time profile much better than a baseline compartment model. Because of its combination of sufficient accuracy and speed of prediction, we found the model to be fit-for-purpose for practical lead optimization.

Significance of Basal Membrane Permeability of Epithelial Cells in Predicting Intestinal Drug Absorption.

Drug absorption from the gastrointestinal tract is often restricted by efflux transport by P-glycoprotein (P-gp) and metabolism by CYP3A4. Both localize in the epithelial cells, and thus, their activities are directly affected by the intracellular drug concentration, which should be regulated by the ratio of permeability between apical (A) and basal (B) membranes. In this study, using Caco-2 cells with forced expression of CYP3A4, we assessed the transcellular permeation of A-to-B and B-to-A directions and the efflux from the preloaded cells to both sides of 12 representative P-gp or CYP3A4 substrate drugs and obtained the parameters for permeabilities, transport, metabolism, and unbound fraction in the enterocytes (fent) using simultaneous and dynamic model analysis. The membrane permeability ratios for B to A (RBA) and fent varied by 8.8-fold and by more than 3000-fold, respectively, among the drugs. The RBA values for digoxin, repaglinide, fexofenadine, and atorvastatin were greater than 1.0 (3.44, 2.39, 2.27, and 1.90, respectively) in the presence of a P-gp inhibitor, thus suggesting the potential involvement of transporters in the B membrane. The Michaelis constant for quinidine for P-gp transport was 0.077 µM for the intracellular unbound concentration. These parameters were used to predict overall intestinal availability (FAFG) by applying an intestinal pharmacokinetic model, advanced translocation model (ATOM), in which permeability of A and B membranes accounted separately. The model predicted changes in the absorption location for P-gp substrates according to its inhibition, and FAFG values of 10 of 12 drugs, including quinidine at varying doses, were explained appropriately. SIGNIFICANCE STATEMENT: Pharmacokinetics has improved predictability by identifying the molecular entities of metabolism and transport and by using mathematical models to appropriately describe drug concentrations at the locations where they act. However, analyses of intestinal absorption so far have not been able to accurately consider the concentrations in the epithelial cells where P-glycoprotein and CYP3A4 exert effects. In this study, the limitation was removed by measuring the apical and basal membrane permeability separately and then analyzing these values using new appropriate models.

Stabilization and quantitative measurement of nicotinamide adenine dinucleotide in human whole blood using dried blood spot sampling.

Nicotinamide adenine dinucleotide (NAD+) is a coenzyme essential for energy production. Recently, associations between NAD+ and aging-related diseases have been reported, and NAD+ precursors that increase NAD+ concentration in the body have been acknowledged as anti-aging supplements. However, there have been only a few studies on the link between aging or aging-related diseases and human blood NAD+ concentration because NAD+ and its precursors are unstable in blood and difficult to measure. Therefore, we aimed to construct a quantitative NAD+ measurement method that is simpler than the existing methods. The calibration standards of NAD+ showed good linearity (0.9936 to 0.9990) in the range of 0.25 to 200 µM, and the lower limit of quantification was 0.5 to 2 µM. We found that QIAcard FTA DMPK-B maintained NAD+ stability of 85% or more for at least 2 weeks at 4 °C and 1 week at room temperature using the dried blood spot method. Additionally, NAD+ stability in the blood extraction solution was more than 90% for 2 months. To our knowledge, there has been no report on a quantitative NAD+ measurement method in human whole blood that can be performed with as little as 5 µL of blood and can be easily implemented at both medical clinics and private homes. Our simple and convenient method has the potential to become the gold standard for NAD+ measurement in blood. It is expected to contribute to the acceleration of research on the correlation between aging or aging-related diseases and NAD+ concentration in human blood.

Cerebral and extracerebral distribution of radioactivity associated with oxytocin in rabbits after intranasal administration: Comparison of TTA-121, a newly developed oxytocin formulation, with Syntocinon.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder associated with deficits in social interactions/communication. Despite the large number of ASD patients, there is no drug approved to treat its core symptoms. Recently, Syntocinon (oxytocin nasal spray) has been reported to have a therapeutic effect on ASD. However, the disadvantage of Syntocinon for ASD treatment is that 6 puffs/administration are required to achieve the effective pharmacological dose. Furthermore, there are no published reports evaluating the cerebral distribution profile of oxytocin after intranasal administration. TTA-121 is a newly developed intranasal oxytocin formulation with high bioavailability produced by optimizing the physicochemical properties. In this study, we prepared the same formula as Syntocinon as the control formulation (CF), and the cerebral and extracerebral distribution of oxytocin in rabbits after single intranasal administration of 3H-labeled oxytocin formulations-[3H]TTA-121 and [3H]CF were examined and compared. The area under the concentration-time curve to the time of the last quantifiable concentration (AUCt) in the whole brain was 3.6-fold higher in the [3H]TTA-121 group than in the [3H]CF group, indicating increased delivery of radioactivity to the brain by TTA-121 than by CF. Since the distribution profiles showed no notable differences in radioactivity between the olfactory bulb and trigeminal nerve, intranasally-administered oxytocin was probably transferred to the brain via both pathways. The results also showed an increase in radioactivity in the prefrontal area and the precuneus, which are probable sites of pharmacological action as shown in clinical studies using functional magnetic resonance imaging (fMRI), confirming that intranasally-administered oxytocin could reach these tissues.

Relationship between excretion clearance of rhodamine 123 and P-glycoprotein (Pgp) expression induced by representative Pgp inducers.

P-Glycoprotein (Pgp) locates in several tissues in the living body and acts as an efflux pump for many drugs. In this study, the usefulness of intravenous rhodamine 123 (Rho123) administration as a marker for detecting the inducing effect of Pgp by drugs was identified, and the relationship between excretion clearances of Rho123 via Pgp and its expression during treatment with the representative Pgp inducers rifampicin (RFP), dexamethasone (DEX) and St. John's Wort (SJW) were examined in rat liver, intestine and kidney. After pretreatment with RFP (10 mg/kg/d) for 4 d, DEX (50 mg/kg/d) for 4 d or SJW (15 mg/kg/d) for 7 d orally, the biliary excretion of Rho123 after intravenous administration (0.2 mg/kg) increased significantly by 40%, 55% and 14%, respectively, and the intestinal excretion increased significantly by 24%, 50% and 27%, respectively, as compared with the controls. In contrast, there were no notable changes in the urinary excretion of Rho123 among rats that received these inducers. Western blot analysis with a monoclonal antibody for Pgp (C219) showed that Pgp levels in the small intestine and liver in the inducer-treated rats increased markedly as compared with the controls. In addition, there was a significant correlation between the induction levels of Pgp in the liver or small intestine and their clearance ratios (r2=0.7583, p<0.05), but not in the kidney. These observations suggest that the excretion clearances of Rho123 from blood circulation to the small intestine or to the bile after its intravenous administration are useful indicators to assess the Pgp function in the presence of Pgp inducers.

In vivo effects of cyclosporin A and ketoconazole on the pharmacokinetics of representative substrates for P-glycoprotein and cytochrome P450 (CYP) 3A in rats.

In this study, the in vivo effects of cyclosporin A (CsA) and ketoconazole (KCZ), which are used as inhibitors of P-glycoprotein (Pgp) and cytochrome P450 (CYP) 3A, respectively, on the pharmacokinetics of rhodamine 123 (Rho123), nelfinavir (NFV) and erythromycin (EM) were evaluated in rats. The biliary excretion clearance (Clbile) of a known Pgp substrate, Rho123, after intravenous pretreatment with CsA or KCZ (0-20 mg/kg after i.v. administration) showed maximum reduction by 85.6 or 54.1%, respectively, suggesting that the inhibitory potency of KCZ is about half that of Pgp in the liver. Without pretreatment with CsA or KCZ, the clearance ratio of Clbile relative to the total body clearances of Rho123, NFV and EM was 10.5, 0.07 and 31.1%, respectively. After CsA pretreatment, these ratios decreased markedly in a manner dependent on the dose of CsA, while after CZ pretreatment the clearance ratios of NFV and EM increased significantly in a manner dependent on the dose of KCZ. However, in the liver, the contribution of Pgp to the changes in the pharmacokinetic parameters of Rho123, NFV and EM after intravenous administration was much less than that of CYP3A. The portal levels of Rho123 and EM but not NFV after intra-loop administration in the presence of 10 microM CsA in the jejunal loop increased significantly, while in the presence of 25 microM KCZ in the jejunal loop, the portal levels of those substrates showed no notable change as compared to the control levels. In conclusion, KCZ had dual potency to inhibit CYP3A and Pgp, and its inhibitory potency for Pgp was half that of CsA in the rat liver. In addition, metabolism via CYP3A contributed more significantly to the clearance of these substrates that did excretion via Pgp in the liver. In the small intestine, the contribution of Pgp is a more important factor in determining the oral bioavailability of EM than metabolism via CYP enzymes. The elimination of NFV is mainly dependent on liver metabolism via CYP3A, and the Pgp efflux mechanism in the liver and intestine did not contribute as importantly to the oral bioavailability of NFV under in vivo conditions, although NFV has been demonstrated to be a substrate of Pgp under in vitro conditions.

Effect of chronic administration of ritonavir on function of cytochrome P450 3A and P-glycoprotein in rats.

Ritonavir (RTV) is well known as an inhibitor of many drugs that are metabolized by cytochrome P450 (CYP) 3A or fluxed via P-glycoprotein (Pgp), although it is also reported that RTV is a potent inducer for them. In this study, to elucidate these contradictory phenomena, functional changes of CYP3A or Pgp during chronic administration of RTV were examined in rats. After pretreatment with RTV for indicated days (day 3-day 14), rats were used in the experiments. The area under the plasma drug concentration vs. time curve (AUC(0-infinity)) after oral administration of RTV (20 mg/kg) to these rats showed an RTV-treatment period-dependent decrease, and the mean AUC(0-infinity) of RTV in Day 14 rats decreased significantly by 57% as compared to the control. The AUC(0-infinity) after intravenous (i.v.) administration of RTV to Day 3 and Day 5 rats increased significantly by 28% and 22%, respectively, while there were no significant changes in the AUC(0-infinity) in Day 7 and Day 14 rats as compared to the control. As for i.v. administration of erythromycin (EM) or midazolam (MDZ) to RTV-treated rats, the AUC(0-infinity)in Day 3 and Day 5 rats increased significantly as compared to the control, while in Day 7 rats and rifampicin-treated rats, the AUC(0-infinity) of EM decreased significantly by 82% and 42%, respectively, as compared to the control. For MDZ, there were no significant changes in the AUC(0-infinity) in Day 7 or Day 14 rats. After i.v. administration of rhodamine123 (Rho123), the excretion clearances from blood circulation to the intestinal lumen and the biliary excretion clearances in Day 14 rats increased markedly by 2.2-fold and 2.6-fold as compared to the control. It has been confirmed that RTV is not only a potent inhibitor but also a potent inducer of CYP3A, and that RTV is a potent inducer of intestinal Pgp. This property of RTV is responsible for regulating the oral bioavailability of drugs that are mediated by CYP3A and Pgp.

Evaluation of factors to decrease plasma concentration of an HIV protease inhibitor, saquinavir in ethanol-treated rats.

Although alcohol consumption is a factor in which the bioavailability of saquinavir (SQV) are retarded, the cause for this phenomenon remains to be uncertain. In the presence study, we examined factors to decrease plasma concentration of SQV in ethanol-treated rats. The ethanol-treated rats were prepared by making them freely access to 15% ethanol solution for 14 d (Day 14 rats). The exsorption clearance of SQV from the blood circulation to the jejunal lumen in the Day 14 rats increased by 6-fold as compared to ethanol non-treated (NT) rats. In the presence of 25 microM ketoconazole (KCZ) or 10 microM cyclosporin A (CsA) in the jejunal lumen, the plasma concentration of SQV in the portal vein increased significantly, and this effect of 10 microM CsA was superior to that of 25 microM KCZ. The biliary excretion clearance of SQV in Day 14 rats also increased by 1.8-fold as compared to that in the NT rats. The metabolic clearance rate (V(max)/K(m)) of SQV in the intestinal microsomes from the Day 14 rats increased significantly, while in the liver microsomes the V(max)/K(m) did not change. The phase II metabolism processes in the Day 14 rats based on UDP-glucuronosyltransferases and gultathion-S-tnrasferase activities were activated, however, they were not likely to be one of factors to decrease the bioavailability of oral SQV, because CYP3A activity in the liver and intestine was not activated to such an extent and SQV itself was not conjugated. These observations suggest that a main possible factor to explain the reducing effect on the SQV oral bioavailability during ethanol consumption is an enhanced efflux of SQV at the intestine and liver, where it is suggested that functional enhancement or excessive expression of P-glycoprotein is caused by ethanol consumption.

Pharmacokinetic characterization of a human immunodeficiency virus protease inhibitor, saquinavir, during ethanol intake in rats.

Throughout therapeutic drug monitoring of human immunodeficiency virus (HIV) protease inhibitors in HIV-infected patients, it was found that plasma concentrations of saquinavir (SQV) were reduced in patients who had a habit of alcohol intake during double protease therapy with SQV and ritonavir (RTV). This study confirmed the pharmacokinetic profiles of SQV during ethanol intake in rats. After oral administration of SQV alone (20 mg/kg) in rats prepared by free access to 15% ethanol solution for 14 days (day 14 rats), the area under the concentration vs time curves (AUC) showed a significant decrease (p<0.01) in comparison with control rats from 0.78+/-0.10 to 0.38+/-0.03 microg h/ml. For intravenous administration of SQV alone (5 mg/kg) to day 14 rats, the total body clearance increased significantly by 1.4-fold (p<0.05), whereas for intracolonic administration of SQV alone, no significant differences in the values of pharmacokinetic parameters were found between control and day 14 rats. With RTV, which has the strongest inhibitory effect on the CYP3A enzyme of the current HIV protease inhibitors, the AUC values of SQV at RTV doses of 2 and 20 mg/kg in day 14 rats also decreased significantly (p<0.01) from 1.30+/-0.06 to 0.57+/-0.05 microg h/ml and from 17.63+/-1.66 to 4.18+/-0.94 microg h/ml, respectively, indicating that the degree of the decrease of AUC values after oral administration with RTV after ethanol intake was larger than the mono-therapy with SQV. This study showed that ethanol-intake decreases the bioavailability of SQV after oral administration alone or with RTV. These observations provide useful information for the treatment of HIV-infected patients when they receive a combination therapy with SQV and RTV, and arouse attention for the effects of alcohol intake.

Effect of indinavir on the intestinal exsorption of amprenavir, saquinavir and nelfinavir after intravenous administration in rats.

To elucidate drug interaction between human immunodeficiency virus (HIV) protease inhibitors (PIs), the effect of indinavir (IDV) on the intestinal exsorption of other HIV PIs, amprenavir (APV), saquinavir (SQV) and nelfinavir (NFV) was investigated in rats using an in situ single perfusion method. IDV in the intestinal perfusate inhibited the exsorption of rhodamine 123 (Rho123), a known P-glycoprotein (P-gp) substrate, from blood into intestinal lumen in a concentration-dependent manner, and the inhibitory potency of 10 micro M IDV in the perfusate was close to that of 10 micro M cyclosporin A (CsA) in the perfusate. Ten micro M of IDV in the intestinal perfusate also decreased significantly the exsorption clearance of Rho123 after intravenous administration. The IDV concentration in this system was not likely to cause hepatic interaction between HIV PIs, because the plasma IDV concentration was far below its inhibition constants for other HIV PIs in the liver microsomes. Thus, 10 micro M of IDV was chosen to investigate the effect of this inhibition on the exsorption of APV, SQV and NFV. IDV in the intestinal perfusate markedly increased the exsorbed amounts of SQV and NFV but not APV after intravenous administrations. Their exsorption clearances, however, showed only a slight increasing tendency or remained unchanged. These findings suggest that in addition to P-gp inhibition, other factors such as CYP3A inhibition might be important in the drug interaction of IDV with APV, SQV and NFV after intravenous administration in rat small intestine. The results obtained in this study will provide useful information to discuss the interactions among PIs when a double protease therapy is used for in HIV-infected patients.

Hepatic and intestinal contributions to pharmacokinetic interaction of indinavir with amprenavir, nelfinavir and saquinavir in rats.

To elucidate the aspects of pharmacokinetic interactions among HIV protease inhibitors (PIs), we investigated the effects of indinavir (IDV) on the hepatic and intestinal first-pass metabolism of other HIV PIs, amprenavir (APV), saquinavir (SQV) and nelfinavir (NFV), in rats. After oral co-administration with IDV, the area under the concentration versus time curves (AUC) of APV, SQV and NFV increased significantly by 1.6-, 9.5- and 2.3-fold, respectively, compared with mono-administration. After intravenous administration, the AUC of APV, SQV and NFV also increased in the presence of IDV by 1.4-, 1.2- and 1.5-fold, respectively. Mean concentrations of APV, SQV and NFV in the liver extracellular fluid, measured using a liver microdialysis method, were very low compared with their Michaelis constants regardless of co-administration of IDV, suggesting that APV, SQV and NFV metabolism follows linear kinetics in the liver. This finding also indicates that metabolism of PIs depended on the metabolic clearance rate in the liver microsomes. The oral bioavailability of SQV in the presence of IDV increased markedly by 8.5-fold, and that of APV and NFV also increased by 1.2- and 1.5-fold, respectively. On the basis of the well-stirred model, the hepatic availabilities of APV, SQV and NFV in the presence of IDV increased by 1.1-, 1.4- and 1.5-fold, and the intestinal availabilities increased by 1.1-, 6.2- and 1.1-fold, respectively. These results suggest that both hepatic and intestinal metabolism were essentially involved in the interactions between IDV and other HIV PIs, and the degree of those contributions varied with each combination of HIV PIs.

Sensitive and simultaneous determination of HIV protease inhibitors in rat biological samples by liquid chromatography-mass spectrometry.

A sensitive and simultaneous liquid chromatographic-mass spectrometric (LC/MS) method for the determination of current four HIV protease inhibitors (PIs), indinavir (IDV), saquinavir (SQV), nelfinavir (NFV) and amprenavir (APV) in rat plasma and liver dialysate by a microdialysis method was described. An isocratic LC/MS method in combination with atmospheric pressure chemical ionization was developed for the determination of these four PIs in biological samples in the same run. The analytes including an internal standard were extracted from 100 microL of plasma or 150 microL of liver dialysate samples by salting-out with 100 microL of ice-cold 2 M K(3)PO(4) followed by ether extraction. The separation of analytes was carried out on a reversed-phase semi-micro column using 50% of acetonitrile containing 1% acetic acid as mobile phase at a flow rate of 0.2mL/min(-1). The separation was completed within 5 min. Precision, recovery and limits of detection indicated that the method was suitable for the quantitative determination of these PIs in rat plasma or liver dialysate. This simple, sensitive and highly specific LC/MS method is suitable for pharmacokinetic studies and therapeutic drug monitoring in AIDS patients who receive double protease therapy.