樋坂 章博

As Parkinson's disease (PD) progresses, there are multiple biomarker changes, and sex and genetic variants may influence the rate of progression. Data-driven, long-term disease progression model analysis may provide precise knowledge of the relationships between these risk factors and progression and would allow for the selection of appropriate diagnosis and treatment according to disease progression. To construct a long-term disease progression model of PD based on multiple biomarkers and evaluate the effects of sex and leucine-rich repeat kinase 2 (LRRK2) mutations, a technique derived from the nonlinear mixed-effects model (Statistical Restoration of Fragmented Time course [SReFT]) was applied to datasets of patients provided by the Parkinson's Progression Markers Initiative. Four biomarkers, including the Unified PD Rating Scale, were used, and a covariate analysis was performed to investigate the effects of sex and LRRK2-related mutations. A model of disease progression over ~30 years was successfully developed using patient data with a median of 6 years. Covariate analysis suggested that female sex and LRRK2 G2019S mutations were associated with 21.6% and 25.4% significantly slower progression, respectively. LRRK2 rs76904798 mutation also tended to delay disease progression by 10.4% but the difference was not significant. In conclusion, a long-term PD progression model was successfully constructed using SReFT from relatively short-term individual patient observations and depicted nonlinear changes in relevant biomarkers and their covariates, including sex and genetic variants.

In vitro methods are widely used in modern toxicological testing; however, the data cannot be directly employed for risk assessment. In vivo toxicity of chemicals can be predicted from in vitro data using physiologically based toxicokinetic (PBTK) modelling-facilitated reverse dosimetry (PBTK-RD). In this study, a minimal-PBTK model was constructed to predict the in-vivo kinetic profile of fenarimol (FNL) in rats and humans. The model was verified by comparing the observed and predicted pharmacokinetics of FNL for rats (calibrator) and further applied to humans. Using the PBTK-RD approach, the reported in vitro developmental toxicity data for FNL was translated to in vivo dose-response data to predict the assay equivalent oral dose in rats and humans. The predicted assay equivalent rat oral dose (36.46 mg/kg) was comparable to the literature reported in vivo BMD10 value (22.8 mg/kg). The model was also employed to derive the chemical-specific adjustment factor (CSAF) for interspecies toxicokinetics variability of FNL. Further, Monte Carlo simulations were performed to predict the population variability in the plasma concentration of FNL and to derive CSAF for intersubject human kinetic differences. The comparison of CSAF values for interspecies and intersubject toxicokinetic variability with their respective default values revealed that the applied uncertainty factors were adequately protective.

BACKGROUND: The aim of this study was to identify significant factors affecting the effectiveness of exercise training using information of the HF-ACTION (Heart Failure: A Controlled Trial Investigating Outcomes of Exercise Training) study. METHODS: Background factors influencing the effect of exercise training were comprehensively surveyed for 2,130 patients by multivariable Cox regression analysis with the stepwise variable selection, and only significant factors were selected that were statistically distinguished from dummy noise factors using the Boruta method. RESULTS: The analysis suggested that the use of beta-blockers, pulse pressure, hemoglobin level, electrocardiography findings, body mass index, and history of stroke at baseline potentially influenced the exercise effect on all-cause death (AD). Therefore, a hypothetical score to estimate the effect of exercise training was constructed based on the analysis. The analysis suggested that the score is useful in identifying patients for whom exercise training may be significantly effective in reducing all-caused death and hospitalization (ADH) as well as AD. Such a subpopulation accounted for approximately 40% of the overall study population. On the other hand, in approximately 45% of patients, the effect of exercise was unclear on either AD or ADH. In the remaining 15% of patients, it was estimated that the effect of exercise might be unclear for ADH and potentially rather increase AD. CONCLUSIONS: This study is the first analysis to comprehensively evaluate the effects of various factors on the outcome of exercise training in chronic heart failure, underscoring the need to carefully consider the patient's background before recommending exercise training. However, it should be noted that exercise training can improve many outcomes in a wide variety of diseases. Therefore, given the limitations involved in post-hoc analyses of a single clinical trial, the characteristics of patients to whom the results of this analysis can be applied need attention, and also further research is necessary on the relationship between the degree of exercise and the outcomes. A new clinical trial would be needed to confirm the factors detected and the appropriateness of the score.

In this study, the ethnic ratios (ERs) of oral clearance between Japanese and Western populations were subjected to model-based meta-analysis (MBMA) for 81 drugs evaluated in 673 clinical studies. The drugs were classified into eight groups according to the clearance mechanism, and the ER for each group was inferred together with interindividual variability (IIV), interstudy variability (ISV), and inter-drug variability within a group (IDV) using the Markov chain Monte Carlo (MCMC) method. The ER, IIV, ISV, and IDV were dependent on the clearance mechanism, and, except for particular groups such as drugs metabolized by polymorphic enzymes or their clearance mechanism is not confirmative, the ethnic difference was found to be generally small. The IIV was well-matched across ethnicities, and the ISV was approximately half of the IIV as the coefficient of variation. To adequately assess ethnic differences in oral clearance without false detections, phase I studies should be designed with full consideration of the mechanism of clearance. This study suggests that the methodology of classifying drugs based on the mechanism that causes ethnic differences and performing MBMA with statistical techniques such as MCMC analysis is helpful for a rational understanding of ethnic differences and for strategic drug development.

Background With the increased use of immune checkpoint inhibitors (ICIs), side effects and toxicity are a great concern. Anaphylaxis has been identified as a potential adverse event induced by ICIs. Anaphylaxis is a life-threatening medical emergency. However, the mechanisms and factors that can potentially influence the incidence and severity of anaphylaxis in patients with cancer remain unclear. Methods Healthy, murine colon 26, CT26, breast 4T1, EMT6, and renal RENCA tumor-bearing mice were treated with an anti-PD-L1 antibody (clone 10F.9G2). Symptoms of anaphylaxis were evaluated along with body temperature and mortality. The amounts of antidrug antibody and platelet-activating factor (PAF) in the blood were quantified via ELISA and liquid chromatography-mass spectrometry (LC-MS/MS). Immune cells were analyzed and isolated using a flow cytometer and magnetic-activated cell sorting, respectively. Results Repeated administration of the anti-PD-L1 antibody 10F.9G2 to tumor-bearing mice caused fatal anaphylaxis, depending on the type of tumor model. After administration, antidrug immunoglobulin G (IgG), but not IgE antibodies, were produced, and PAF was released as a chemical mediator during anaphylaxis, indicating that anaphylaxis was caused by an IgG-dependent pathway. Anaphylaxis induced by 10F.9G2 was treated with a PAF receptor antagonist. We identified that neutrophils and macrophages were PAF-producing effector cells during anaphylaxis, and the tumor-bearing models with increased numbers of neutrophils and macrophages showed lethal anaphylaxis after treatment with 10F.9G2. Depletion of both neutrophils and macrophages using clodronate liposomes prevented anaphylaxis in tumor-bearing mice. Conclusions Thus, increased numbers of neutrophils and macrophages associated with cancer progression may be risk factors for anaphylaxis. These findings may provide useful insights into the mechanism of anaphylaxis following the administration of immune checkpoint inhibitors in human subjects.

OBJECTIVE: This study aimed to examine the magnitude of age-related change in hepatic clearance by integrating the data of multiple drugs and to compare this with renal clearance, considering associations with age-related changes in organ weight and blood flow. METHODS: The results of multiple population pharmacokinetic analyses that detected age-related clearance changes in hepatically eliminated drugs were collected. The relationship between hepatic clearance of the unbound drug and age was then analyzed using the nonlinear least-squares method, adjusting for interdrug differences. The obtained change in hepatic clearance was compared with age-related changes in liver weight and hepatic blood flow in Japanese and Westerners. For comparison, the changes in renal clearance were analyzed similarly. RESULTS: In total, 18 drugs were analyzed. The hepatic unbound clearance decreased by 32% at age 80 years and by 40% at age 90 years, compared with age 40 years, suggesting that it decreased by 0.80% per year with aging. The rate of the decrease was consistent with decreases in hepatic weight per person or blood flow per person, regardless of ethnicity and sex. Since age-related change in body weight varied somewhat by sex or ethnicity, hepatic weight per body weight was less consistent to account for age-related change in hepatic clearance. As for an index of renal clearance, the changes in inulin clearance with age were similar to those in renal blood flow, with a decrease of 0.97% per year from the age of 40 years. CONCLUSIONS: Hepatic clearance consistently decreased by 0.80% per year from the age of 40 years, with aging for multiple drugs analyzed in this study. Changes in organ weight and blood flow are considered to be primarily responsible for the age-related changes in hepatic and renal clearance.

<title>Abstract</title>Despite progress in the use of hyperthermia in clinical practice, the thermosensitivity of cancer cells is poorly understood. In a previous study, we found that sensitivity to hyperthermia varied between ovarian and uterine cancer cell lines. Upon hyperthermia, glycolytic enzymes decreased in hyperthermia-resistant SKOV3 cells. However, the mechanisms of glycolysis inhibition and their relationship with thermoresistance remain to be explored. In this study, metabolomic analysis indicated the downregulation of glycolytic metabolites in SKOV3 cells after hyperthermia. Proteomic and pathway analyses predicted that the ubiquitin pathway was explicitly activated in resistant SKOV3 cells, compared with hyperthermia-sensitive A2780 cells, and STUB1, a ubiquitin ligase, potentially targeted PKM, a glycolytic rate-limiting enzyme. PKM is degraded via ubiquitination upon hyperthermia. Although glycolysis is inactivated by hyperthermia, ATP production is maintained. We observed that oxygen consumption and mitochondrial membrane potential were activated in SKOV3 cells but suppressed in A2780 cells. The activation of mitochondria could compensate for the loss of ATP production due to the suppression of glycolysis by hyperthermia. Although the physiological significance has not yet been elucidated, our results demonstrated that metabolomic adaptation from the Warburg effect to mitochondrial oxidative phosphorylation could contribute to thermoresistance in ovarian and uterine cancer cells.

The second messenger 2'3'-cyclic-GMP-AMP (cGAMP) is thought to be transmitted from brain carcinomas to astrocytes via gap junctions, which functions to promote metastasis in the brain parenchyma. In the current study, we established a method to introduce cGAMP into astrocytes, which simulates the state of astrocytes that have been invaded by cGAMP around tumors. Astrocytes incorporating cGAMP were analyzed by metabolomics, which demonstrated that cGAMP increased glutamate production and astrocyte secretion. The same trend was observed for γ-aminobutyric acid (GABA). Conversely, glutamine production and secretion were decreased by cGAMP treatment. Due to the fundamental role of astrocytes in regulation of the glutamine-glutamate cycle, such metabolic changes may represent a potential mechanism and therapeutic target for alteration of the central nervous system (CNS) environment and the malignant transformation of brain carcinomas.

The role of astrocytes in the periphery of metastatic brain tumors is unclear. Since astrocytes regulate central nervous metabolism, we hypothesized that changes in astrocytes induced by contact with cancer cells would appear in the metabolome of both cells and contribute to malignant transformation. Coculture of astrocytes with breast cancer cell supernatants altered glutamate (Glu)-centered arginine-proline metabolism. Similarly, the metabolome of cancer cells was also altered by astrocyte culture supernatants, and the changes were further amplified in astrocytes exposed to Glu. Inhibition of Glu uptake in astrocytes reduces the variability in cancer cells. Principal component analysis of the cancer cells revealed that all these changes were in the first principal component (PC1) axis, where the responsible metabolites were involved in the metabolism of the arginine-proline, pyrimidine, and pentose phosphate pathways. The contribution of these changes to the tumor microenvironment needs to be further pursued.

To characterize and compare various medicines for chronic heart failure (CHF), changes in circulatory physiological parameter during pharmacotherapy were investigated by a model-based meta-analysis (MBMA) of circulatory physiology. The clinical data from 61 studies mostly in patients with heart failure with reduced ejection fraction (HFrEF), reporting changes in heart rate, blood pressure, or ventricular volumes after treatment with carvedilol, metoprolol, bisoprolol, bucindolol, enalapril, aliskiren, or felodipine, were analyzed. Seven cardiac and vasculature function indices were estimated without invasive measurements using models based on appropriate assumptions, and their correlations with the mortality were assessed. Estimated myocardial oxygen consumption, a cardiac load index, correlated excellently with the mortality at 3, 6, and 12 months after treatment initiation, and it explained differences in mortality across the different medications. The analysis based on the present models were reasonably consistent with the hypothesis that the treatment of HFrEF with various medications is due to effectively reducing the cardiac load. Assessment of circulatory physiological parameters by using MBMA would be insightful for quantitative understanding of CHF treatment.

Precise prediction of drug absorption is key to the success of new drug development and efficacious pharmacotherapy. In this study, we developed a new absorption model, the advanced translocation model (ATOM), by extending our previous model, the translocation model. ATOM reproduces the translocation of a substance in the intestinal lumen using a partial differential equation with variable dispersion and convection terms to describe natural flow and micro-mixing within the intestine, under not only fasted but also fed conditions. In comparison with ATOM, it was suggested that a conventional absorption model, advanced compartmental absorption and transit model, tends to underestimate micro-mixing in the upper intestine, and it is difficult to adequately describe movements under the fasted and fed conditions. ATOM explains the observed nonlinear absorption of midazolam successfully, with a minimal number of scaling factors. Furthermore, ATOM considers the apical and basolateral membrane permeabilities of enterocytes separately and assumes compartmentation of the lamina propria, including blood vessels, to consider intestinal blood flow appropriately. ATOM estimates changes in the intestinal availability caused by drug interaction associated with inhibition of CYP3A and P-gp in the intestine. Additionally, ATOM can estimate the drug absorption in the fed state considering delayed intestinal drug flow. Therefore, ATOM is a useful tool for the analysis of local pharmacokinetics in the gastrointestinal tract, especially for the estimation of nonlinear drug absorption that may involve various interactions with intestinal contents or other drugs. Significance Statement The advanced translocation model (ATOM) was newly developed that precisely explains various movements of intestinal contents including the fasted and fed conditions which cannot be adequately described by the current physiological pharmacokinetic models.

Despite the promising anticancer effects of immune checkpoint inhibitors, their low objective response rate remains to be resolved; thus, combination therapies have been investigated. We investigated the combination of an anti-programmed cell death 1 (aPD-1) monoclonal antibody with the knockdown of vascular endothelial factor receptor 2 (VEGFR2) on tumor endothelial cells to overcome resistance to immune checkpoint inhibitors and improve the objective response rate. The successful delivery of small interfering RNA to tumor endothelial cells was achieved by RGD peptide-modified lipid nanoparticles composed of a novel, pH-sensitive, and biodegradable ssPalmO-Phe. RGD-modified lipid nanoparticles efficiently induced the knockdown of VEGFR2 in tumor endothelial cells (TECs), which induced vascular normalization. The combination of a PD-1 monoclonal antibody with Vegfr2 knockdown enhanced CD8+ T cell infiltration into tumors and successfully suppressed tumor growth and improved response rate compared with monotherapy. Our combination approach provides a promising strategy to improve therapeutic outcomes in immune checkpoint inhibitor-resistant cancers.

In pharmacokinetic (PK) analysis, conventional models are described by ordinary differential equations (ODE) that are generally solved in their Laplace transformed forms. The solution in the Laplace transformed forms is inverse Laplace transformed to derive an analytical solution. However, inverse Laplace transform is often mathematically difficult. Consequently, numerical inverse Laplace transform methods have been developed. In this study, we focus on extending the modeling functions of Nonlinear Mixed Effect Model (NONMEM), a standard software for PK and population pharmacokinetic (PPK) analyses, by adding the Fast Inversion of Laplace Transform (FILT) method, one of the representative numerical inverse Laplace transform methods. We implemented PREDFILT, a specialized PRED subroutine, which functions as an internal model unit in NONMEM to enable versatile FILT analysis with second-order precision. The calculation results of the compartment models and a dispersion model are in good agreement with the ordinary analytical solutions and theoretical values. Therefore, PREDFILT ensures enhanced flexibility in PK or PPK analyses under NONMEM environments.

The aim of this study was to elucidate the lifelong disease progression of chronic obstructive pulmonary disease (COPD) with biomarker changes and identify their influencing factors, by utilizing a new analysis method, Statistical Restoration of Fragmented Time-course (SReFT). Individual patient data (n = 1025) participating in the Study to Understand Mortality and MorbidITy (SUMMIT, NCT01313676), which was collected within the observational period of 4 years, were analyzed. The SReFT analysis suggested that scores of St. George's Respiratory Questionnaire and COPD assessment test, representative scores of the health-related quality of life (HRQOL) questionnaire, increased consistently for 30 years of disease progression, which was not detected by conventional analysis with a linear mixed effect model. It was estimated by the SReFT analysis that normalized forced expiratory volume in one second for age, sex, and body size (%FEV1) reduced for the initial 10 years from the onset of the disease but thereafter remained constant. The analysis of HRQOL scores and lung functions suggested that smoking cessation slowed COPD progression by approximately half and that exacerbation accelerated it considerably. In conclusion, this retrospective study utilizing SReFT elucidated the progression of COPD over 30 years and associated quantitative changes in the HRQOL scores and lung functions.

Global differences in changes in the numbers of population-adjusted daily test-positive cases (NPDP) and deaths (NPDD) by COVID-19 were analyzed for 49 countries, including developed and developing countries. The changes as a proportion of national population were compared, adjusting by the beginning of test-positive cases increase (BPI) or deaths increase (BDI). Remarkable regional differences of more than 100-fold in NPDP and NPDD were observed. The trajectories of NPDD after BDI increased exponentially within 20 days in most countries. Machine learning analysis suggested that NPDD on 30 days after BDI was the highest in developed Western countries (1180 persons per hundred million), followed by countries in the Middle East (128), Latin America (97), and Asia (7). Furthermore, in Western countries with positive rates of the PCR test of less than 7.0%, the increase in NPDP was slowing-down two weeks after BPI, and subsequent NPDD was only 15% compared with those with higher positive rates, which suggested that the situation of testing might have affected the velocity of COVID-19 spread. The causes behind remarkable differences between regions possibly include genetic factors of inhabitants because distributions of the race and of the observed infection increasing rates were in good agreement globally.

Ceramide kinase (CerK) phosphorylates ceramide to ceramide-1-phosphate (C1P), a bioactive sphingolipid. Since the mechanisms responsible for regulating the proliferation and migration/metastasis of cancer cells by the CerK/C1P pathway remain unclear, we conducted the present study. The knockdown of CerK in A549 lung and MCF-7 breast cancer cells (shCerK cells) increased the formation of lamellipodia, which are membrane protrusions coupled with cell migration. Mouse embryonic fibroblasts prepared from CerK-null mice also showed an enhanced formation of lamellipodia. The overexpression of CerK inhibited lamellipodium formation in A549 cells. The knockdown of CerK increased the number of cells having lamellipodia with Rac1 and the levels of active Rac1-GTP form, whereas the overexpression of CerK decreased them. CerK was located in lamellipodia after the epidermal growth factor treatment, indicating that CerK functioned there to inhibit Rac1. The migration of A549 cells was negatively regulated by CerK. An intravenous injection of A549-shCerK cells into nude mice resulted in markedly stronger metastatic responses in the lungs than an injection of control cells. The in vitro growth of A549 cells and in vivo expansion after the injection into mouse flanks were not affected by the CerK knockdown. These results suggest that the activation of CerK/C1P pathway has inhibitory roles on lamellipodium formation, migration, and metastasis of A549 lung cancer cells.

BACKGROUND: Recently, antiprogrammed cell death protein 1 (aPD-1) and antiprogrammed death-ligand 1 (aPD-L1) monoclonal antibodies (mAbs) have been approved. Even though aPD-1 and aPD-L1 mAbs target the same PD-1/PD-L1 axis, it is still unclear whether both mAbs exert equivalent pharmacological activity in patients who are sensitive to PD-1/PD-L1 blockade therapy, as there is no direct comparison of their pharmacokinetics (PK) and antitumor effects. Therefore, we evaluated the differences between both mAbs in PK and therapeutic effects in PD-1/PD-L1 blockade-sensitive mouse models. METHODS: Herein, murine breast MM48 and colon MC38 xenografts were used to analyze the pharmacological activity of aPD-1 and aPD-L1 mAbs. The PK of the mAbs in the tumor-bearing mice was investigated at low and high doses using two radioisotopes (Indium-111 and Iodine-125) to evaluate the accumulation and degradation of the mAbs. RESULTS: aPD-1 mAb showed antitumor effect in a dose-dependent manner, indicating that the tumor model was sensitive to PD-1/PD-L1 blockade therapy, whereas aPD-L1 mAb failed to suppress tumor growth. The PK study showed that aPD-L1 mAb was accumulated largely in normal organs such as the spleen, liver, and kidney, resulting in low blood concentration and low distributions to tumors at a low dose, even though the tumors expressed PD-L1. Sufficient accumulation of aPD-L1 mAb in tumors was achieved by administration at a high dose owing to the saturation of target-mediated binding in healthy organs. However, degradation of aPD-L1 mAb in tumors was greater than that of aPD-1 mAb, which resulted in poor outcome presumably due to less inhibition of PD-L1 by aPD-L1 mAb than that of PD-1 by aPD-1 mAb. CONCLUSION: According to the PK studies, aPD-1 mAb showed linear PK, whereas aPD-L1 mAb showed non-linear PK between low and high doses. Collectively, the poor PK characteristics of aPD-L1 mAb caused lower antitumor activity than of aPD-1 mAb. These results clearly indicated that aPD-L1 mAb required higher doses than aPD-1 mAb in clinical setting. Thus, targeting of PD-1 would be more advantageous than PD-L1 in terms of PK.

Extract of Cyclolepis genistoides D. Don (vernacular name Palo azul; Palo) are traditionally consumed in the Republic of Paraguay in South America for the treatment of diabetes and kidney disease, and is sold in Japan as dietary supplement. This study aimed to elucidate the mechanism of anti-diabetes activity of Palo, especially focused on insulin resistance. Palo promoted adipocytes differentiation and regulated adipokine profiles in 3T3-L1 adipocytes by modulation of PPARγ, a major regulator of adipose differentiation. Human adipocyte showed almost similar profile with 3T3-L1 against Palo treatment. Furthermore, Palo treatment (250 or 1000 mg/kg) was performed with C57BL/6J mice for 14 weeks, being fed high-fat-diet (HFD60) simultaneously. Palo 250 mg/kg exhibited a tendency to decrease subcutaneous adipose volume along with increase of PPARγ and its target, adiponectin mRNA expression. In addition, as the other insulin targeted cell, effect on muscle differentiation was examined. Palo increased differentiation of C2C12 mouse muscle myoblasts by increase of IGF-1, myogenin, and myosine heavy chain (MHC) as well as 5'-AMP-activated protein kinase (AMPK) activation. Palo subsequently promoted myotube formation under differentiation condition. From the above, it was clarified that Palo acts variously on the differentiation and maturation of both adipocytes and muscle cells, and from the viewpoint of the regulatory mechanism for adipocytes, PPARγ-inducing action was shown to be a mechanism that acts across species.

Aims: To use a mouse model of imiquimod-induced psoriasis to investigate the relationship between pruritus and mast cells, nerve growth factor (NGF) and endogenous pruritogenic peptides, which are highly expressed in the skin of psoriasis patients. Main methods: We developed a mouse model of imiquimod-induced psoriasis and measured the frequency and duration of the model animals' self-scratching behavior using the SCLABA®-Real real-time scratch counting system. We then harvested the ears and subjected them to toluidine blue staining and real-time PCR. Key findings: Topical application of imiquimod increased the Psoriasis Area and Severity Index score as well as the frequency and duration of self-scratching. Regarding internal factors, increases in mast cells number and mRNA expression of NGF and endogenous pruritogenic peptide precursor were confirmed. Significance: Self-scratching behavior is accompanied by increased number of mast cells and expression of NGF and endogenous pruritogenic peptides in our imiquimod-induced psoriasis model. The expression of these factors was consistent with the features in patients with pruritic psoriasis, suggesting that our model reflects at least some of the precipitating factors of pruritus found in humans.

Cancer cachexia is a systemic wasting syndrome characterized by anorexia and loss of body weight. The xanthine oxidase (XO) inhibitor febuxostat is one of the promising candidates for cancer cachexia treatment. However, cachexic symptoms were not alleviated by oral administration of febuxostat in our cancer cachexia model. Metabolomic analysis with brains of our cachexic model showed that purine metabolism was activated and XO activity was increased, and thus suggested that febuxostat would not reach the brain. Accordingly, targeting XO in the brain, which controls appetite, may be an effective strategy for treatment of cancer cachexia.

Although sunitinib is the first-line drug for progressive renal cell carcinoma (RCC), most patients experience its tolerance. One possible way of overcoming drug resistance is combination therapy. Epigenetic modifier is one of the candidate drug group. A recent evidence suggests that cell metabolism is regulated by epigenetic mechanisms. Epigenetic abnormalities lead to changes in metabolism and may contribute to drug resistance and progression of RCC. Consequently, we investigated whether trichostatin A (TSA), a potent histone-deacetylase (HDAC) inhibitor, alters sunitinib-induced cytotoxicity and metabolism in RCC cells at epigenetic regulatory concentrations. Combined metabolome and transcriptome analysis suggested that TSA impacts on energy productive metabolic pathways, such as those involving TCA cycle and nucleotide metabolism especially for increase of hyperphosphorylated form. Combination of sunitinib and TSA increased cell death with PARP cleavage, an early marker of mitochondrial apoptosis, whereas receptor tyrosine kinase signaling, which is the target of sunitinib, was not altered by TSA. Finally, the established sunitinib resistant-RCC cell (786-O Res) was also exposed to sunitinib and TSA combination, resulting in significant growth inhibition. In summary, it was suggested that TSA reduces sunitinib resistance by triggering intracellular metabolome shifts regarding energy metabolism, that is the first recognized mechanism as an HDAC inhibitor.

Clinical observations of patients with chronic diseases are often restricted in terms of duration. Therefore, obtaining a quantitative and comprehensive understanding of the chronology of chronic diseases is challenging, because of the inability to precisely estimate the patient's disease stage at the time point of observation. We developed a novel method to reconstitute long-term disease progression from temporally fragmented data by extending the nonlinear mixed-effects model to incorporate the estimation of "disease time" of each subject. Application of this method to sporadic Alzheimer's disease successfully depicted disease progression over 20 years. The covariate analysis revealed earlier onset of amyloid-β accumulation in male and female apolipoprotein E ε4 homozygotes, whereas disease progression was remarkably slower in female ε3 homozygotes compared with female ε4 carriers and males. Simulation of a clinical trial suggests patient recruitment using the information of precise disease time of each patient will decrease the sample size required for clinical trials.

The noninferiority of direct oral factor Xa (FXa) inhibitors (rivaroxaban, apixaban, and edoxaban) in treatment of atrial fibrillation were demonstrated compared with warfarin by several large clinical trials; however, subsequent meta-analyses reported a higher risk of major bleeding with rivaroxaban than with the other FXa inhibitors. In the present study, we first estimated the changes of prothrombin time (PT) in 5 randomized trials based on reported population pharmacokinetic and pharmacodynamic models and then carried out a model-based meta-analysis to obtain models describing the relationship between PT changes and the event rates of ischemic stroke/systemic embolism (SE) and of major bleeding. By using the models, we simulated the optimal therapeutic doses for each FXa inhibitor. It was suggested that dose reduction of rivaroxaban from the current 20 mg/d to 10 mg/d would decrease patient deaths from major bleeding (hazard ratio [HR], 0.69; 95% confidence interval [CI], 0.64-0.74) with little increase in those for ischemic stroke/SE (HR, 1.11; 95% CI, 1.07-1.20). The overall decrease in the mortality caused by both events was estimated as 5.81 per 10 000 patient-years (95% CI, 3.92-8.16), with an HR of 0.87 (95% CI, 0.83-0.91). For apixaban and edoxaban, no distinct change in the overall mortality was simulated by dose modification. This study suggested that the current dose of rivaroxaban might be excessive and would need to be reduced to decrease the excess risk of major bleeding.

Bosentan is a substrate of hepatic uptake transporter organic anion-transporting polypeptides (OATPs), and undergoes extensive hepatic metabolism by cytochrome P450 (P450), namely, CYP3A4 and CYP2C9. Several clinical investigations have reported a nonlinear relationship between bosentan doses and its systemic exposure, which likely involves the saturation of OATP-mediated uptake, P450-mediated metabolism, or both in the liver. Yet, the underlying causes for the nonlinear bosentan pharmacokinetics are not fully delineated. To address this, we performed physiologically based pharmacokinetic (PBPK) modeling analyses for bosentan after its intravenous administration at different doses. As a bottom-up approach, PBPK modeling analyses were performed using in vitro kinetic parameters, other relevant parameters, and scaling factors. As top-down approaches, three different types of PBPK models that incorporate the saturation of hepatic uptake, metabolism, or both were compared. The prediction from the bottom-up approach (models 1 and 2) yielded blood bosentan concentration-time profiles and their systemic clearance values that were not in good agreement with the clinically observed data. From top-down approaches (models 3, 4, 5-1, and 5-2), the prediction accuracy was best only with the incorporation of the saturable hepatic uptake for bosentan. Taken together, the PBPK models for bosentan were successfully established, and the comparison of different PBPK models identified the saturation of the hepatic uptake process as a major contributing factor for the nonlinear pharmacokinetics of bosentan.

PURPOSE: To investigate the effect of excipients on latanoprost penetration into the aqueous humor with clinically available 6 products with different solutions mainly in the types and concentrations of preservatives. METHODS: In 363 patients with cataracts, we instilled 1 latanoprost drop in 1 eye before surgery. The drop was randomly selected by brand name product (A) and 5 generic products (B-F) composed with different excipients. B contains similar excipients to A. C and D contain lower concentrations of benzalkonium chloride than A. E includes sodium benzoate, and F contains no preservatives. At 0.5-1, 3, and 6 h after instillation, samples of aqueous humor were collected to determine the latanoprost free acid by mass spectrometry. The time course of intraocular concentration and the areas under the aqueous humor latanoprost free acid concentration-time curves (AUCs) were calculated. RESULTS: At 0.5-1 h, the aqueous humor concentration of latanoprost free acid was 8.5 ± 1.0 ng/mL for C, which was significantly higher (P < 0.01) than that of A (3.4 ± 0.5 ng/mL). At 3 and 6 h, however, no significant difference was noted in the concentrations of latanoprost free acid between the brand name and generic products. For each of the generic products, the peak free acid concentration was above the known threshold concentration for biological activity. At 6 h postdose, however, the levels of latanoprost free acid were below the threshold for Products C, E, and F. Comparisons of AUC0-6h and AUC0-24h values showed that these parameters were the greatest with A, and E was significantly inferior to A (P < 0.05). CONCLUSIONS: Currently available latanoprost solutions with different preservatives showed sufficient intraocular concentration to activate the FP receptor, but different pharmacokinetic profiles of absorption or elimination.

Acquired resistance to sunitinib is a challenge in the treatment of renal cell carcinoma (RCC). The dysregulation of cellular metabolism is prevalent during resistance acquisition. It is known that in sunitinib-resistant RCC 786-O (786-O Res) cells sunitinib is mainly sequestered in the intracellular lysosomes. However, the relevance between sunitinib resistance and cellular metabolism has not been examined. In this study, we examined the metabolic changes in 786-O Res by using capillary electrophoresis-time of flight mass spectrometry. The cell line 786-O Res was established via persistent treatment with sunitinib, where increase in intracellular sunitinib, and sizes of lysosomes and nuclei were enhanced as compared with those in the parental 786-O (786-O Par) cells. Metabolic analyses revealed that out of the 110 metabolites examined, 13 were up-regulated and 4 were down-regulated in the 786-O Res cells. The glycolysis, tricarboxylic acid cycle and pentose phosphate pathway (PPP) were identified as being altered in the sunitinib-resistant cells, which resulted in the enhanced metabolisms of energy, nucleic acids, and glutathione redox cycle. As sunitinib was sequestered in the enlarged lysosomes in 786-O Res, the enriched energy metabolism might contribute to the maintenance of luminal pH in lysosomes via the H+ ATPase. The changes in the PPP could contribute to nuclei enlargement through up-regulation of nucleic acid biosynthesis and protect 786-O Res from cytotoxicity induced by sunitinib through up-regulation of reduced glutathione. Though the direct link between sunitinib resistance and metabolic alternation remains to be elucidated, this metabolomics study provides fundamental insights into acquisition of sunitinib resistance.

In this study, impact of a polymorphism of CYP2C19 on drug-drug interaction (DDI) was examined for etizolam. The effect of itraconazole (a strong CYP3A inhibitor) on the pharmacokinetics of etizolam (a substrate of CYP2C19 and CYP3A) was assessed in both extensive metabolizers (EMs) and poor metabolizers (PMs) of CYP2C19. Sixteen participants (8 EMs and 8 PMs) received a single oral dose of etizolam (0.25 mg) on day 1. The participants ingested itraconazole (200 mg twice a day) on days 2-5. On day 5, participants received an oral dose of etizolam (0.25 mg) again. Before coadministration of itraconazole (day 1), the area under the time-plasma concentration curve from time zero to infinity (AUC∞ ) of etizolam was higher in PMs than in EMs (2.65-fold, P < .01). Coadministration of itraconazole increased the AUC∞ of etizolam 1.66-fold and 2.34-fold in EMs and PMs, respectively (day 5). Consequently, AUC∞ was 6.18-fold higher in PMs with itraconazole than that in EMs without itraconazole. The increase by itraconazole was larger in PMs (P < .01). In heterozygous EMs (hEMs), AUC∞ was simulated to be 2.56-fold higher with itraconazole than that in EMs without itraconazole. We found that in vitro measurements of fraction metabolized (fm ) using the liver microsome prepared from PM donors would be helpful to predict polymorphism-dependent DDIs. These results suggest that the PMs and hEMs of a polymorphic CYP would be at higher risk of DDIs relative to EMs for drugs metabolized by both polymorphic and nonpolymorphic CYPs such as etizolam.

In this study, intestinal drug-drug interactions (DDIs) for substrate drugs of P-glycoprotein were simulated extensively using the extended QGut model and translocation model to explore the determinants of DDI. The results of analyses using both models suggested that permeability and active efflux clearance were the major factors that influenced the fraction absorbed (FA). The results of simulation for 100 virtual drugs in which parameters were generated considering the actual values of commercially available drugs suggested that the ratio of the pH-corrected passive permeability to the intrinsic efflux clearance (Pu/CLeff) relative to that of digoxin would be a useful and quantitative index of P-glycoprotein (P-gp)-mediated DDI risk at lower doses. At higher doses, such as 100 mg, the risk of P-gp-mediated DDI would be significantly reduced because of saturation of P-gp efflux. The simulation suggested that although drugs with lower permeability were more susceptible, even drugs with higher permeability than metoprolol, a representative highly permeable drug, such as BCS class 1 and 2, may experience DDIs owing to P-gp inhibition. Overall, this study demonstrated the usefulness of mathematical intestinal models when only limited observational data are available.

The constituent protein of gap junctions, connexin (Cx), interacts with various proteins via its C-terminus region, including kinases, cell-adhesion proteins, and a pro-apoptotic protein, Bax. This molecular interaction may affect expression and functioning of the interacting proteins and modulate the cellular physiology. In our previous work, Cx43 was found to interact directly with Bax and in the presence of sunitinib, lead to the Bax-mediated apoptosis in mesothelioma cells. In this study, we investigated the mechanism of how Cx43 promotes Bax-mediated apoptosis using the same cell line. Treatment with sunitinib increased the expression of the active conformation of the Bax protein, which was predominantly localized at the mitochondria, only in Cx43-transfected cells. Bax oligomerization and decrease in the mitochondrial membrane potential were also observed. The involvement of c-Jun N-terminal kinase (JNK) in the interaction of Cx43 and Bax was further examined. Treatment with sunitinib increased the expression of phosphorylated (active) form of JNK only in the Cx43-transfected cells. Phosphorylated JNK and active Bax were co-localized, and the co-localization was suppressed by the knockdown of Cx43. Moreover, JNK inhibition clearly suppressed Bax activation. In conclusion, we identified a novel Cx43-JNK-Bax axis regulating the process of apoptosis for the first time.

Aprepitant is a known inducer of CYP2C9, the main warfarin-metabolizing enzyme. Consequently, co-administration of these two drugs may result in reduction of the anticoagulation activity of warfarin. However, the nature and degree of time-dependent changes in prothrombin time international normalized ratio (PT-INR) after aprepitant and warfarin co-treatment in patients receiving anticancer chemotherapy has not been elucidated. We retrospectively examined the changes in warfarin dose, PT-INR, and warfarin sensitivity index (WSI; average of PT-INR value/average of daily warfarin dose) during four weeks, i.e., one week before and three weeks after aprepitant administration. The mean and standard deviation values of WSI for one week before and one, two, and three weeks after the beginning of aprepitant administration were 0.51±0.22 (1.00, n=34), 0.74±0.30 (1.53±0.59, n=30), 0.38±0.15 (0.82±0.22, n=28), and 0.46±0.29 (0.87±0.23, n=24), respectively. Values in parentheses represent relative changes versus WSI of one week before and number of subjects. Although the mean value of WSI significantly increased one week after aprepitant administration compared to that at one week before the administration, it in turn significantly decreased two weeks after compared to one week before (paired t-test, p<0.05 after Bonferoni correction). In patients taking warfarin, PT-INR should be carefully monitored for at least two weeks after the beginning of aprepitant administration because it may fluctuate with both aprepitant and chemotherapy during this period.

PURPOSE: Drug transfer into milk is of concern due to the unnecessary exposure of infants to drugs. Proposed prediction methods for such transfer assume only passive drug diffusion across the mammary epithelium. This study reorganized data from the literature to assess the contribution of carrier-mediated transport to drug transfer into milk, and to improve the predictability thereof. METHODS: Milk-to-plasma drug concentration ratios (M/Ps) in humans were exhaustively collected from the literature and converted into observed unbound concentration ratios (M/Punbound,obs). The ratios were also predicted based on passive diffusion across the mammary epithelium (M/Punbound,pred). An in vitro transport assay was performed for selected drugs in breast cancer resistance protein (BCRP)-expressing cell monolayers. RESULTS: M/Punbound,obs and M/Punbound,pred values were compared for 166 drugs. M/Punbound,obs values were 1.5 times or more higher than M/Punbound,pred values for as many as 13 out of 16 known BCRP substrates, reconfirming BCRP as the predominant transporter contributing to secretory transfer of drugs into milk. Predictability of M/P values for selected BCRP substrates and non-substrates was improved by considering in vitro-evaluated BCRP-mediated transport relative to passive diffusion alone. CONCLUSIONS: The current analysis improved the predictability of drug transfer into milk, particularly for BCRP substrates, based on an exhaustive data overhaul followed by focused in vitro transport experimentation.

This study aimed to construct a new local pharmacokinetic model of gastrointestinal absorption, the translocation model (TLM), using an anatomically relevant, minimally segmented structure to explain linear and nonlinear intestinal absorption, metabolism, and transport. The TLM was based on the concept of a single absorption site that flexibly moves, expands, and shrinks along with the length of the gastrointestinal tract after the intake of an oral dose. The structure of the small intestine is continuous, and various time- and location-dependent issues are freely incorporated in the analysis. Since the model has only one absorption site, understanding and modification of factors affecting absorption are simple. The absorption site is composed of four compartments: solid drug in the lumen, solution drug in the lumen, concentration in the enterocytes, and concentration in the lamina propria. The lamina propria includes the blood capillaries. Blood flow in the absorption site of the lamina propria appropriately accounts for the absorption. In the TLM, the permeability of the apical membrane and that of the basolateral membrane are distinct. By considering plicate, villi, and microvilli expansions of the surface area, the apparent permeability measured in Caco-2 experiments was converted to the effective permeability in vivo. The intestinal availability, bioavailability, and dose product of intestinal availability and absorption rate relationship of the model drugs were well explained using the TLM. The TLM would be a useful tool for the consideration of local pharmacokinetics in the gastrointestinal tract in various situations.

5-[(2-Chloro-6-fluorophenyl)acetylamino]-3-(4-fluorophenyl)-4-(4-pyrimidinyl)isoxazole (AKP-001) is a potent p38 mitogen-activated protein kinase inhibitor that is being developed to specifically target the intestines for the treatment of inflammatory bowel disease. According to the ante-drug concept, AKP-001 was designed to be metabolized to inactive forms via the first-pass metabolism to avoid undesirable systemic exposure. The purpose of this study is to investigate the pharmacokinetic characteristics of AKP-001 and its metabolites (M1 and M2) in rats, utilizing a simple physiologically based pharmacokinetic (PBPK) model. In vitro metabolic activity of AKP-001 in the S9 fraction of rat liver was examined, and plasma concentration-time profiles were developed following intravenous and/or oral administration of AKP-001 and its metabolites. AKP-001 was primarily metabolized to M1; however, M2 was not detected in liver S9 fractions. In accordance with this observation in vitro, M2 was detected in plasma after oral dosing of AKP-001 with a lag time of 1.5 hours, but not after intravenous dosing. To analyze pharmacokinetics in rats in vivo, a simple PBPK model was developed by simultaneous fitting of the plasma concentrations after treatment with AKP-001 and its metabolites. The observed plasma concentration-time profiles of AKP-001 and metabolites were described by the model adequately. Intestinal and systemic exposures of AKP-001 were simulated using the model to assess the relationship between pharmacokinetics and efficacy/safety. Model analysis suggested that oral bioavailability of intestine-targeting ante-drugs should be low to avoid systemic side effects. The pharmacokinetic properties of AKP-001 meet this criterion owing to extensive first-pass metabolism.

CASE: We describe two cases in which treatment with aprepitant persistently altered antithrombotic control in patients receiving warfarin. A 60-year-old man received 5 weekly cycles of chemotherapy. Aprepitant was administered as a 3-day regimen from the second cycle of chemotherapy. In each of the chemotherapy cycles that included aprepitant, the therapeutic international normalized ratio (INR) decreased markedly to <1.6, and slowly recovered over several weeks. A 47-year-old woman was treated with 4 weekly cycles of chemotherapy. Aprepitant was administered as a 3-day regimen. On day 8 of the first cycle of chemotherapy, the patient's INR fell markedly to 1.1. Although warfarin dosage was steadily increased over the four subsequent cycles of chemotherapy, therapeutic target range was not recovered. INR gradually returned to normal during the 2 months after the final cycle of chemotherapy. CONCLUSION: To our knowledge, this is the first case report to document the effects of aprepitant in cancer patients receiving anticoagulation therapy.

In this study, we developed the drug-drug interaction (DDI) method as a new assessment technique of intestinal availability (F(G), the fraction of drug transferred from the intestinal enterocytes into the liver, escaping from intestinal metabolism) based on the clearance theory. This method evaluates F(G) from changes caused by DDIs in the area under the blood concentration-time curve and in the elimination half-life of victim drugs. Application of the DDI method to data from the literature revealed that the mean and S.D. of F(G) values for 20 substrate drugs of CYP3A was 0.56 ± 0.29, whereas that for 8 substrate drugs of CYP2C9, CYP2C19, and CYP2D6 was 0.86 ± 0.11. These results were consistent with the fact that intestinal metabolism is mediated predominantly by CYP3A. The DDI method showed reasonable correlations with the conventional i.v./p.o. method and the grape fruit juice (GFJ) method (coefficients of determination of 0.41 and 0.81, respectively). The i.v./p.o. method was more susceptible to fluctuations in the hepatic blood flow rate compared with the DDI and GFJ methods. The DDI method evaluates F(G) separating from the absorption ratio (F(A)) although it requires approximation of F(A). Since preciseness of approximation of F(A) does not greatly affect the evaluation of F(G) by the DDI method, we proposed a reasonable approximation method of F(A) for the evaluation of F(G) in the DDI method. The DDI method would be applicable to a broad range of situations in which various DDI data are utilizable.

Drug transfer into milk is a general concern during lactation. So far, breast cancer resistance protein (Bcrp) is the only transporter known to be involved in this process, whereas participation of other transporters remains unclear. We investigated the importance of organic cation transporter (Oct) in drug transfer into milk in mice. The mammary glands of lactating versus nonlactating FVB strain mice revealed elevated mRNA levels of Oct1 and Bcrp, whereas Oct2 and Oct3 mRNA levels were decreased. Specific uptake of cimetidine, acyclovir, metformin, and terbutaline was observed in human embryonic kidney 293 cells transfected with murine Oct1 or Oct2. The milk-to-plasma concentration ratio (M/P) values of cimetidine and acyclovir were significantly decreased in Bcrp knockout and Oct1/2 double-knockout (DKO) mice compared with control FVB mice, whereas the M/P values of terbutaline and metformin were significantly decreased in Oct1/2 DKO mice alone. These are the first to suggest that Oct1 might be involved in secretory transfer of substrate drugs into milk.

OBJECTIVES: Thrombocytopenia is sometimes observed during linezolid therapy. Here, we aimed to investigate the factors affecting linezolid-induced thrombocytopenia. METHODS: A prospective observational study was performed between October 2009 and February 2011; 30 patients were included. Plasma linezolid trough concentrations were measured on days 3, 7 and 14 after initial drug administration. Platelet counts and haemoglobin levels were also monitored. RESULTS: Thrombocytopenia occurred in 17 patients (56.7%). Median linezolid trough concentrations on day 3 were significantly higher in patients with renal impairment (creatinine clearance <60 mL/min) than in patients without renal impairment (14.7 versus 4.8 mg/L; P < 0.0001). Median linezolid trough concentrations on day 3 in patients who developed thrombocytopenia were also significantly higher than those in patients who did not (13.4 versus 4.3 mg/L, P < 0.0001). Development of thrombocytopenia occurred significantly more frequently in patients with linezolid trough concentration >7.5 mg/L (OR, 90.0; P < 0.0001) and renal impairment (OR, 39.0; P = 0.0002). The Kaplan-Meier plot showed that the median time from the initiation of therapy to development of thrombocytopenia was 11 days. CONCLUSIONS: Patients with renal impairment are more likely to have a high plasma linezolid concentration. In addition, a high plasma linezolid concentration and renal impairment significantly affected the development of linezolid-induced thrombocytopenia. Further studies are required to evaluate whether therapeutic drug monitoring-guided dosage adjustment of linezolid decreases the adverse effects while maintaining treatment efficacy in patients with renal dysfunction.

PURPOSE: Drug transfer into milk is a general concern during lactation. Because data are limited in human subjects, particularly for new drugs, experimental animal models of lactational drug transfer are critical. This study analyzed drug transfer into milk in a mouse model, as well as the contribution of similar and dissimilar host factors. METHODS: Milk/plasma drug concentration ratios (M/P) in humans were obtained from the literature, while those in mice were determined experimentally after intraperitoneal implantation of osmotic pumps containing drugs of interest. Unbound drug fractions in plasma and milk were determined in vitro for both species. RESULTS: M/P values were determined for 27 drugs in mice and compared with those in human. These values were increased in mice for 21 drugs; the geometric mean ratio of M/P between mice and humans was 2.03 (95% CI, 1.42-2.89) for all 27 drugs. These results were reasonably explained by the relatively high protein and lipid content in mouse milk. Moreover, species-specific asymmetrical transport systems were suggested for 9 drugs. CONCLUSIONS: In addition to species-specific differences in milk protein and lipid content, variances in asymmetrical drug transport across the mammary epithelium may yield discordant M/P values in humans and mice.

The plasma concentration of repaglinide is reported to increase greatly when given after repeated oral administration of itraconazole and gemfibrozil. The present study analyzed this interaction based on a physiologically based pharmacokinetic (PBPK) model incorporating inhibition of the hepatic uptake transporter and metabolic enzymes involved in repaglinide disposition. Firstly, the plasma concentration profiles of inhibitors (itraconazole, gemfibrozil, and gemfibrozil glucuronide) were reproduced by a PBPK model to obtain their pharmacokinetic parameters. The plasma concentration profiles of repaglinide were then analyzed by a PBPK model, together with those of the inhibitors, assuming a competitive inhibition of CYP3A4 by itraconazole, mechanism-based inhibition of CYP2C8 by gemfibrozil glucuronide, and inhibition of organic anion transporting polypeptide (OATP) 1B1 by gemfibrozil and its glucuronide. The plasma concentration profiles of repaglinide were well reproduced by the PBPK model based on the above assumptions, and the optimized values for the inhibition constants (0.0676 nM for itraconazole against CYP3A4; 14.2 μM for gemfibrozil against OATP1B1; and 5.48 μM for gemfibrozil glucuronide against OATP1B1) and the fraction of repaglinide metabolized by CYP2C8 (0.801) were consistent with the reported values. The validity of the obtained parameters was further confirmed by sensitivity analyses and by reproducing the repaglinide concentration increase produced by concomitant gemfibrozil administration at various timings/doses. The present findings suggested that the reported concentration increase of repaglinide, suggestive of synergistic effects of the coadministered inhibitors, can be quantitatively explained by the simultaneous inhibition of the multiple clearance pathways of repaglinide.

Induction of cytochrome P450 3A4 (CYP3A4) expression is often implicated in clinically relevant drug-drug interactions (DDI), as metabolism catalyzed by this enzyme is the dominant route of elimination for many drugs. Although several DDI models have been proposed, none have comprehensively considered the effects of enzyme transcription/translation dynamics on induction-based DDI. Rifampicin is a well-known CYP3A4 inducer, and is commonly used as a positive control for evaluating the CYP3A4 induction potential of test compounds. Herein, we report the compilation of in vitro induction data for CYP3A4 by rifampicin in human hepatocytes, and the transcription/translation model developed for this enzyme using an extended least squares method that can account for inherent inter-individual variability. We also developed physiologically based pharmacokinetic (PBPK) models for the CYP3A4 inducer and CYP3A4 substrates. Finally, we demonstrated that rifampicin-induced DDI can be predicted with reasonable accuracy, and that a static model can be used to simulate DDI once the blood concentration of the inducer reaches a steady state following repeated dosing. This dynamic PBPK-based DDI model was implemented on a new multi-hierarchical physiology simulation platform named PhysioDesigner.

In clinical settings, decrease of renal function is frequently observed in patients treated with vancomycin (VCM). In this study, mutable covariates models (MCMs) were constructed to analyze the pharmacokinetics (PK) of VCM with Bayesian estimation, considering time-dependent decreases in creatinine clearance in 23 patients with decreasing renal function. The predicted mean percentage error (MPE) of VCM concentrations analyzed with a conventional fixed covariates model (FCM) was -19.1%, whereas the MPE was improved to 2.5% by applying MCM. Furthermore, a probable lag time between fluctuations in VCM clearance (CL(VCM) ) and serum creatinine (S(cr) ) was analyzed by MCM, MCM(Lag1d) , and MCM(Lag2d) , which considered lag times of 0, 1, and 2 days, respectively. Compared with FCM, all MCMs improved fitness with the significantly decreased root mean square percentage error (RMSPE) and MPE. However, RMSPE and MPE analyzed with MCM were not significantly different from those with MCM(Lag1d) and MCM(Lag2d) , indicating that lag times between alterations in CL(VCM) and S(cr) were obscure in these patients. Collectively, these results suggest that PK parameters of VCM were more accurately calculated by MCMs than by conventional FCM, and that VCM dosages calculated by FCM would be overestimated by approximately 20% in patients with decreasing renal function.

The aim of the study was to quantitatively predict the clearance of three antibiotics, amikacin, vancomycin, and teicoplanin, during continuous hemodiafiltration (CHDF) and to propose their optimal dosage in patients receiving CHDF. For this goal, in vitro CHDF experiments with a polyacrylonitrile (PAN) membrane were first performed using these antibiotics, and then the clearances were compared with in vivo CHDF situations determined in 16 critically ill patients. The in vitro CHDF clearances were described as the product of the outflow rate of a drain (Q(outflow)) and the drug unbound fraction in artificial plasma, indicating that drug adsorption to the PAN membrane has minor effect on drug clearance in our settings. The observed in vivo clearances also agreed very well with the predicted values, with a product of Q(outflow) and plasma unbound fraction, when residual creatinine clearance (CL(CR)) was taken into account (within a range of 0.67- to 1.5-fold for 15 of 16 patients). Based on these results, a nomogram of the optimized dosages of amikacin, vancomycin, and teicoplanin was proposed, and it was evident that Q(outflow) and residual CL(CR) are major determinants of the dosage and dosing interval for these antibiotics. Although the applicability needs to be confirmed with another type of membrane or higher Q(outflow), our nomogram can help determine the dosage setting in critically ill patients receiving CHDF.

Medication use during lactation is a matter of concern due to unnecessary exposure of infants to drugs. Although some studies have predicted the extent of drug transfer into milk from physicochemical parameters, drug concentration-time profiles in milk have not been predicted or even analyzed yet. In the present study, a drug transfer model was constructed by defining secretion and reuptake clearances (CL(sec) and CL(re), respectively) between milk and plasma based on unbound drug concentrations. Through the use of this model, drug concentration-time profiles were analyzed in human milk and plasma based on data collected from the literature. CL(sec) and CL(re) values were obtained successfully for 49 drugs. Because the CL(sec) and CL(re) values were in general similar for each drug, transport across the mammary epithelia was mediated by passive diffusion in most cases. This study demonstrated that the logarithmically transformed values of CL(sec) and CL(re) can be predicted from physicochemical parameters with adjusted R(2) values of 0.705 and 0.472, respectively. Moreover, 66.7 and 77.8% of predicted CL(sec) and CL(re) values were within 3-fold error ranges of the observed values for 45 and 27 drugs, respectively. Finally, time profiles of drug concentrations in milk were simulated from physicochemical parameters. The milk-to-plasma area under the concentration-time curve ratios also were predicted successfully within 3-fold error ranges of the observed values for 71.9% of the drugs analyzed. The method described herein therefore may be useful in predicting drug concentration-time profiles in human milk for newly developed drugs.