Handling Delayed or Missed Dose of Antiseizure Medications: A Model-Informed Individual Remedial Dosing.

BACKGROUND AND OBJECTIVES: Delayed or missed antiseizure medications (ASMs) doses are common during long-term or lifelong antiepilepsy treatment. This study aims to explore optimal individualized remedial dosing regimens for delayed or missed doses of 11 commonly used ASMs. METHODS: To explore remedial dosing regimens, Monte Carlo simulation was used based on previously identified and published population pharmacokinetic models. Six remedial strategies for delayed or missed doses were investigated. The deviation time outside the individual therapeutic range was used to evaluate each remedial regimen. The influences of patients' demographics, concomitant medication, and scheduled dosing intervals on remedial regimens were assessed. RxODE and Shiny in R were used to perform Monte Carlo simulation and recommend individual remedial regimens. RESULTS: The recommended remedial regimens were highly correlated with delayed time, scheduled dosing interval, and half-life of the ASM. Moreover, the optimal remedial regimens for pediatric and adult patients were different. The renal function, along with concomitant medication that affects the clearance of the ASM, may also influence the remedial regimens. A web-based dashboard was developed to provide individualized remedial regimens for the delayed or missed dose, and a user-defined module with all parameters that could be defined flexibly by the user was also built. DISCUSSION: Monte Carlo simulation based on population pharmacokinetic models may provide a rational approach to propose remedial regimens for delayed or missed doses of ASMs in pediatric and adult patients with epilepsy.

The positive role of trade openness in industrial green total factor productivity-provincial evidence from China.

The sustainable development of the economy is a key issue of global concern. Green total factor productivity (GTFP) combining economic growth with resources and the environment can evaluate the quality of economic development comprehensively and accurately. In this paper, super slack-based measure (Super-SBM) model and Malmquist-Luenberger (ML) index were used to calculate GTFP. The trend of industrial GTFP in China's 30 provinces from 2006-2015 was analyzed. Furtherly, a dynamic panel model was used to discuss the impact of trade openness on GTFP. The results showed that (1) the growth rate of GTFP rose from 2007 to 2011 and declined significantly from 2011 to 2015, and GTFP only achieved positive growth in 2011; (2) the growth rate of GTFP in the eastern region was higher than that in the central and western regions; (3) the trend of technical progress change (MLTECH) index was highly consistent with that of ML index. That was, technical progress played a major role in the variation of GTFP; (4) trade openness could significantly improve China's GTFP. Every 1% increase in trade openness could increase GTFP by 0.097% on average. It is advisable to implement differentiated economic development and environmental policies in different regions. Meanwhile, relevant measures can be taken to promote import and export trade, such as encouraging companies to increase investment in green technology research and development, optimizing the trade environment.

Assessment and Feedback Control of Paving Quality of Earth-Rock Dam Based on OODA Loop.

Paving thickness and evenness are two key factors that affect the paving operation quality of earth-rock dams. However, in the recent study, both of the key factors characterising the paving quality were measured using finite point random sampling, which resulted in subjectivity in the detection and a lag in the feedback control. At the same time, the on-site control of the paving operation quality based on experience results in a poor and unreliable paving quality. To address the above issues, in this study, a novel assessment and feedback control framework for the paving operation quality based on the observe-orient-decide-act (OODA) loop is presented. First, in the observation module, a cellular automaton is used to convert the location of the bulldozer obtained by monitoring devices into the paving thickness of the levelling layer. Second, in the orient module, the learning automaton is used to update the state of the corresponding and surrounding cells. Third, in the decision module, an overall path planning method is developed to realise feedback control of the paving thickness and evenness. Finally, in the act module, the paving thickness and evenness of the entire work unit are calculated and compared to their control thresholds to determine whether to proceed with the next OODA loop. The experiments show that the proposed method can maintain the paving thickness less than the designed standard value and effectively prevent the occurrence of ultra-thick or ultra-thin phenomena. Furthermore, the paving evenness is improved by 21.5% as compared to that obtained with the conventional paving quality control method. The framework of the paving quality assessment and feedback control proposed in this paper has extensive popularisation and application value for the same paving construction scene.

How to handle the delayed or missed dose of rivaroxaban in patients with non-valvular atrial fibrillation: model-informed remedial dosing.

BACKGROUND: Rivaroxaban is an oral anticoagulant widely used for stroke prevention in patients with non-valvular atrial fibrillation (NVAF). During long-term anticoagulant therapy, delayed or missed doses are common. This study aimed to explore appropriate remedial dosing regimens for non-adherent rivaroxaban-treated patients. METHODS: Monte Carlo simulation based on a previously established rivaroxaban population pharmacokinetic/pharmacodynamic (PK/PD) model for patients with NVAF was employed to design remedial dosing regimens. The proposed regimens were compared with remedial strategies in the European Heart Rhythm Association (EHRA) guide by assessing deviation time in terms of drug concentration, factor Xa activity, and prothrombin time. RESULTS: The proposed remedial dosing regimens were dependent on delay duration. The missed dose should be taken immediately when the delay does not exceed 6 h; a half dose is advisable when the delay is between 6 and 20 h. A missed dose should be skipped if less than 4 h remains before the next dose. The proposed regimens resulted in shorter deviation time than that of the EHRA guide. CONCLUSION: PK/PD modeling and simulation provide valid evidence on the remedial dosing regimen of rivaroxaban, which could help to minimize the risk of bleeding and thromboembolism.

Establishment of a population pharmacokinetics model of vancomycin in 94 infants with septicemia and its application in individualized therapy.

BACKGROUND: We aim to develop a population pharmacokinetics (PopPK) model of vancomycin for the treatment of septicemia in infants younger than one year. Factors influence of the PK was investigated to optimize vancomycin dosing regimen. METHODS: The nonlinear mixed effects modelling software (NONMEM) was used to develop the PopPK model of vancomycin. The stability and predictive ability of the final model were assessed by using normalized prediction distribution errors (NPDE) and bootstrap methods. The final model was subjected to Monte Carlo simulation in order to determine the optimal dose. RESULTS: A total of 205 trough and peak concentrations in 94 infants (0-1 year of age) with septicemia were analyzed. The interindividual variability of the PK parameter was described by the exponential model. Residual error was better described by the proportional model than the mixed proportional and addition models. Serum creatinine concentration and body weight are the major factors that affect the PK parameters of vancomycin. The clearance was shown to be higher when ceftriaxone was co-treated. More than two model evaluation methods showed better stability than the base model, with superior predictive performance, which can develop individualized dosing regimens for clinical reference. Through prediction of final model, the trough concentration was more likely < 5 mg/L when a routine dose of 10 mg/kg is administered every 6 h to 3-9-month-old infants. Therefore, the dose should be increased in the treatment of infant septicemia. CONCLUSIONS: The stable and effective PopPK model of vancomycin in Chinese infants with septicemia was established. This model has satisfactory predictive ability for clinically individualized dosing regimens in this vulnerable population.

Population pharmacokinetic analysis and dosing guidelines for tacrolimus co-administration with Wuzhi capsule in Chinese renal transplant recipients.

WHAT IS KNOWN AND OBJECTIVES: Tacrolimus (TAC) is a first-line immunosuppressant which is used to prevent transplant rejection after solid organ transplantation (SOT). However, it has a narrow therapeutic index and high individual variability in pharmacokinetics (PK) and pharmacogenomics (PG). It has been reported that the metabolism of TAC can be affected by genetic factors, leading to different rates of metabolism in different subjects. Wuzhi Capsule (WZC) is a commonly used TAC-sparing agent in Chinese SOT to reduce TAC dosing due to its inhibitory effect on TAC metabolism by enzymes of the CYP3A subfamily. The aims of this study were to assess the effect of TAC+WZC co-administration and genetic polymorphism on the pharmacokinetics of TAC, by using a population pharmacokinetic (PPK) model. A dosing guideline for individualized TAC dosing is proposed based on the PPK study. METHODS: The medical records of 165 adult patients with kidney transplant and their 824 TAC concentrations from two kidney transplantation centres were reviewed. The genotypes of four single-nucleotide polymorphisms (SNPs) in CYP3A5*3 and ABCB1 (rs1128503, rs2032582 and rs1045642) were tested by MASSARRAY. A PPK model was constructed by nonlinear mixed effect model (NONMEM® , Version 7.3). Finally, Monte Carlo simulations were employed to design initial dosing regimens based on the final model. RESULTS AND DISCUSSION: The one-compartmental PPK model with first-order absorption and elimination of TAC was established in kidney transplant recipients (KTRs). CYP3A5*3 had significant impact on the PPK model. The haematocrit (HCT), postoperative time (POD) and CYP3A5*3 genotypes had a significant influence on TAC clearance when combined with WZC. The model was expressed as 23.4 × (HCT/0.3)-0.729  × 0.837 (combination with WZC) × e-0.0875(POD/12.6) ×1.18 (CYP3A5 expressors). For patients carrying the CYP3A5*3/*3 allele and with 30% HCT, the required TAC dose to achieve target trough concentrations of 10-15 ng/ml was 4 mg twice daily (q12h). For patients with the CYP3A5*3/*3 allele, the required dose was 3 mg TAC q12h when combined with WZC, and for patients with the CYP3A5*1/*1 or *1/*3 allele, the required dose was 4 mg of TAC q12h when co-administered with WZC. WHAT IS NEW AND CONCLUSION: Wuzhi Capsule co-administration and CYP3A5 variants affect the PK of TAC Dosing guidelines are made based on the PPK model to allow individualized administration of TAC, especially when co-administered with WZC.

Determination of vancomycin exposure target and individualised dosing recommendations for neonates: model-informed precision dosing.

INTRODUCTION: Few studies incorporating population pharmacokinetic/pharmacodynamic (Pop-PK/PD) modelling have been conducted to quantify the exposure target of vancomycin in neonates. A retrospective observational cohort study was undertaken in neonates to determine this target and dosing recommendations (chictr.org.cn, ChiCTR1900027919). METHODS: A Pop-PK model was developed to estimate PK parameters. Causalities between acute kidney injury (AKI) occurrence and vancomycin use were verified using Naranjo criteria. Thresholds of vancomycin exposure in predicting AKI or efficacy were identified via classification and regression tree analysis. Associations between exposure thresholds and clinical outcomes, including AKI and efficacy, were analysed by logistic regression. Dosing recommendations were designed using Monte Carlo simulations based on the optimised exposure target. RESULTS: Pop-PK modelling included 182 neonates with 411 observations. On covariate analysis, neonatal physiological maturation, renal function and concomitant use of vasoactive agents (VAS) significantly affected vancomycin PK. Seven cases of vancomycin-induced AKI were detected. Area under the concentration-time curve from 0-24 hours (AUC0-24) ≥ 485 mg•h/L was an independent risk factor for AKI after adjusting for VAS co-administration. The clinical efficacy of vancomycin was analysed in 42 patients with blood culture-proven staphylococcal sepsis. AUC0-24 to minimum inhibitory concentration (AUC0-24/MIC) ≥ 234 was the only significant predictor of clinical effectiveness. Monte Carlo simulations indicated that regimens in Neonatal Formulary 7 and Red Book (2018) were unsuitable for all neonates. CONCLUSION: An AUC0-24 of 240-480 (assuming MIC = 1 mg/L) is a recommended exposure target of vancomycin in neonates. Model-informed dosing regimens are valuable in clinical practice.

Population pharmacokinetics of unbound valproic acid in pediatric epilepsy patients in China: a protein binding model.

PURPOSE: The purpose of this study was to establish a protein binding model of unbound valproic acid (VPA) based on Chinese pediatric patients with epilepsy and provide a reference for clinical medication. METHODS: A total of 313 patients were included and both their total and unbound VPA concentrations (375 pairs of concentrations) were measured. NONMEM software was used for population pharmacokinetic modeling. The stepwise method was used to screen the potential covariates. Goodness-of-fit plot, bootstrap, and visual predictive check were used for model evaluation. In addition, dose recommendations for typical patients aged 0 to 16 years were proposed by Monte Carlo simulations. RESULTS: A one-compartment model of first-order absorption and first-order elimination was used to describe the pharmacokinetic characteristics of unbound VPA, and the linear non-saturable binding equation was introduced to describe the protein binding. Body weight, age-based maturation, and co-medicated with lamotrigine could affect the CL/F of unbound and bound VPA. Model evaluation showed satisfactory robustness of the final model. The dosing regimens for children aged 0 to 16 years were proposed based on the final established model. CONCLUSION: We developed a population pharmacokinetic model of unbound and bound VPA that took account of protein binding. The VPA dosing regimen in pediatric patients with epilepsy needs to be optimized by the body weight, age, and co-medications.

Population pharmacokinetic and optimization of polymyxin B dosing in adult patients with various renal functions.

AIMS: Current FDA-approved label recommends that the dosage of polymyxin B should be adjusted according to renal function. However, the correlation between polymyxin B pharmacokinetics (PK) and creatinine clearance (CrCL) is poor. This study aimed to develop a population PK model of polymyxin B in adult patients with various renal functions and to identify a dosing strategy. METHODS: A retrospective PK study was performed in 32 adult patients with various renal function. Nonlinear mixed effects modelling was applied to build a population PK model of polymyxin B followed by Monte Carlo simulations which designed polymyxin B dosing regimens across various renal function. RESULTS: Polymyxin B PK analyses included 112 polymyxin B concentrations at steady state from 32 adult patients, in which 71.9% of them were critically ill. In the final PK model, CrCL was the significant covariate on CL (typical value 1.59 L/h; between-subject variability 13%). The mean (SD) individual empirical Bayesian estimate of CL was 1.75 (0.43) L/h. In addition, a new dosing strategy combining the PK/pharmacodynamic (PD) targets and Monte Carlo simulation indicated that the reduction of polymyxin B dose in patients with renal insufficiency improved the probability of achieving optimal exposure. For severe infections caused by organisms with minimum inhibitory concentration (MIC) ≥ 2 mg/L, a high daily dose of polymyxin B might be possible for bacterial eradication, but the risk of nephrotoxicity is increased. CONCLUSIONS: Renal function plays a significant role in polymyxin B PK, and the dose of polymyxin B should be adjusted according to CrCL in patients with renal insufficiency.

Population pharmacokinetics of the anti-PD-1 antibody camrelizumab in patients with multiple tumor types and model-informed dosing strategy.

Camrelizumab, a programmed cell death 1 (PD-1) inhibitor, has been approved for the treatment of patients with relapsed or refractory classical Hodgkin lymphoma, nasopharyngeal cancer and non-small cell lung cancer. The aim of this study was to perform a population pharmacokinetic (PK) analysis of camrelizumab to quantify the impact of patient characteristics and to investigate the appropriateness of a flat dose in the dosing regimen. A total of 3092 camrelizumab concentrations from 133 patients in four clinical trials with advanced melanoma, relapsed or refractory classical Hodgkin lymphoma and other solid tumor types were analyzed using nonlinear mixed effects modeling. The PKs of camrelizumab were properly described using a two-compartment model with parallel linear and nonlinear clearance. Then, covariate model building was conducted using stepwise forward addition and backward elimination. The results showed that baseline albumin had significant effects on linear clearance, while actual body weight affected intercompartmental clearance. However, their impacts were limited, and no dose adjustments were required. The final model was further evaluated by goodness-of-fit plots, bootstrap procedures, and visual predictive checks and showed satisfactory model performance. Moreover, dosing regimens of 200 mg every 2 weeks and 3 mg/kg every 2 weeks provided similar exposure distributions by model-based Monte Carlo simulation. The population analyses demonstrated that patient characteristics have no clinically meaningful impact on the PKs of camrelizumab and present evidence for no advantage of either the flat dose or weight-based dose regimen for most patients with advanced solid tumors.

Evidence-based Guideline for Therapeutic Drug Monitoring of Vancomycin: 2020 Update by the Division of Therapeutic Drug Monitoring, Chinese Pharmacological Society.

BACKGROUND: Clinical practice guidelines or recommendations often require timely and regular updating as new evidence emerges, because this can alter the risk-benefit trade-off. The scientific process of developing and updating guidelines accompanied by adequate implementation can improve outcomes. To promote better management of patients receiving vancomycin therapy, we updated the guideline for the therapeutic drug monitoring (TDM) of vancomycin published in 2015. METHODS: Our updated recommendations complied with standards for developing trustworthy guidelines, including timeliness and rigor of the updating process, as well as the use of the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) approach. We also followed the methodology handbook published by the National Institute for Health and Clinical Excellence and the Spanish National Health System. RESULTS: We partially updated the 2015 guideline. Apart from adults, the updated guideline also focuses on pediatric patients and neonates requiring intravenous vancomycin therapy. The guideline recommendations involve a broadened range of patients requiring TDM, modified index of TDM (both 24-hour area under the curve and trough concentration), addition regarding the necessity and timing of repeated TDM, and initial dose for specific subpopulations. Overall, 1 recommendation was deleted and 3 recommendations were modified. Eleven new recommendations were added, and no recommendation was made for 2 clinical questions. CONCLUSIONS: We updated an evidence-based guideline regarding the TDM of vancomycin using a rigorous and multidisciplinary approach. The updated guideline provides more comprehensive recommendations to inform rational and optimized vancomycin use and is thus of greater applicability.

Remedial dosing recommendations for delayed or missed doses of valproic acid in patients with epilepsy based on Monte Carlo simulations.

OBJECTIVE: Delayed or missed doses are unavoidable in the pharmacotherapy of epilepsy and significantly compromise the efficacy of antiepileptic drug treatment. An inappropriate remedial regimen can cause seizure relapse or serious adverse events. This study investigated the effect of delayed or missed doses on the pharmacokinetics (PK) of valproic acid (VPA) in patients with epilepsy and established remedial dosing recommendations for nonadherent patients. METHODS: Monte Carlo simulations are based on all previous population pharmacokinetic models for pediatric, adult and elderly patients with epilepsy. The following four remedial strategies were investigated for each delayed dose: A) A partial dose or a regular dose is taken immediately; a regular dose is taken at the next scheduled time. B) The delayed dose was administered immediately, followed by a partial dose at the next scheduled time. C) The delayed dose and a partial dose are taken; the next scheduled time is skipped, and the regular regimen is resumed. D) Double doses are taken when missed one dose or two doses, and the regular regimen at the subsequent scheduled time is resumed. RESULTS: The recommended remedial dose was related to the delay duration and daily dose. Remedial dosing strategies A and B were almost equivalent, whereas Strategy C was recommended when the delayed dose was close to the next scheduled dose. Strategy D was only suggested for delayed two doses. CONCLUSION: Simulations provide quantitative insight into the remedial regimens for nonadherent patients, and clinicians should select the optimal regimen for each patient based on the individual's status.

Remedial dosing recommendations for delayed or missed doses of lamotrigine in pediatric patients with epilepsy using Monte Carlo simulations.

OBJECTIVE: This study investigated the effect of delayed or missed doses on the pharmacokinetics (PK) of lamotrigine (LTG) in children with epilepsy and established remedial dosing recommendations for nonadherent patients. METHODS: The Monte Carlo simulation based on a published LTG population PK model was used to assess the effect of different scenarios of nonadherence and the subsequently administered remedial regimens. The following three remedial approaches were investigated for each delayed dose: A) A partial dose was administered immediately, and the regular dose was administered at the next scheduled time. B) The delayed dose was administered immediately, followed by a partial dose at the next scheduled time. C) The delayed and partial doses were coadministered immediately, the next scheduled dose was skipped, and the regular dosing was resumed at the subsequent scheduled time. The most appropriate remedial regimen was that with the shortest deviation time from the individual therapeutic window. RESULTS: The effect of nonadherence on PK was dependent on the delay duration and daily dose, and the recommended remedial dose was related to the delay duration and concomitant antiepileptic drugs. Remedial dosing strategies A and B were almost equivalent, whereas C showed a larger PK deviation time. If one dose was missed, double doses were not recommended for the next scheduled time. CONCLUSIONS: Simulations provide quantitative insight into the remedial regimens for nonadherent patients, and clinicians should select the optimal regimen based on the status of patients.

[Effects of multiple-trough sampling design and algorithm on the estimation of population and individual pharmacokinetic parameters].

The purpose of this study is to investigate the effects of multiple-trough sampling design and nonlinear mixed effect modeling (NONMEM) algorithm on the estimation of population and individual pharmacokinetic parameters. Oxcarbazepine and tacrolimus were used as one-compartment and two-compartment model drugs, respectively. Seven sampling designs were investigated using various number of trough concentrations per individual ranging from 1-4. Monte Carlo simulations were performed to produce state-steady trough concentrations. One-compartment model was used to fit simulated data from oxcarbazepine and tacrolimus. The accuracy and precision of the estimated parameters were evaluated using the median prediction error (PE), the median absolute PE and boxplot. The results indicated that trough concentrations could yield reliable estimates of apparent clearance (CL/F). For oxcarbazepine, as the number of trough concentrations per subject increased, the accuracy and precision of CL/F, between-subject variability (BSV) of CL/F and residual variability (RUV) tended to be improved. For tacrolimus, however, although no improvement were observed in the accuracy of CL/F and BSV of CL/F, the PE distribution ranges were significantly narrowed and the RUV estimates were less bias and imprecise. In terms of algorithm, Monte Carlo importance sampling (IMP) and IMP assisted by mode a posteriori estimation (IMPMAP) were consistently better than other methods. Additionally, the sampling design had no significant effects on the individual parameter estimates, which were only depended on the interaction between BSV and RUV in various algorithms. Decreased in BSV and RUV levels can improve the accuracy and precision of the estimation for both population and individual pharmacokinetic parameter estimates.

The effect of poor compliance on the pharmacokinetics of carbamazepine and its epoxide metabolite using Monte Carlo simulation.

AIM: To study the effects of delayed and missed doses (poor compliance) on the pharmacokinetics of carbamazepine (CBZ) and its main active metabolite carbamazepine-10,11-epoxide (CBZE) in Chinese epilepsy patients using Monte Carlo simulation. METHODS: CBZ and CBZE time-concentration profiles in various scenarios were generated based on a population pharmacokinetic study in Chinese epilepsy patients using Monte Carlo simulation. The scenarios included patients given multiple doses of CBZ that ranged from 100 to 300 mg three times daily or from 200 to 300 mg every 12 h. The therapeutic range of CBZ and CBZE for each scenario was estimated to assess the effect of delayed or missed doses and to design corresponding rescue regimens. Moreover, the impact of body weight, absorption rate and co-therapy with other antiepileptic drugs (phenytoin, phenobarbital and valproic acid) on the dosage recommendation was investigated in the event of poor compliance. RESULTS: The risk for a sub-therapeutic range of CBZ and CBZE was increased in a dose-dependent manner in both two and three times daily regimens when delayed or missed doses occurred. The effects of poor compliance was less prominent on the lower daily doses compared with those on the higher daily doses. The dose recommendations, in the event of poor compliance, were time related and dose dependent. Patient body weight, absorption rate and co-therapy with phenytoin, phenobarbital and valproic acid had no significant impact on the dose recommendation. CONCLUSION: Patients with epilepsy should take the delayed doses as soon as they remember, and partial missed doses may need to be taken near or at the next scheduled time.

[Estimation of individual pharmacokinetic parameters using maximum a posteriori Bayesian method with D-optimal sampling strategy].

This study was aimed to develop a maximum a posteriori Bayesian (MAPB) estimation method to estimate individual pharmacokinetic parameters based on D-optimal sampling strategy. Meanwhile, the performance of MAPB was compared with the multiple linear regression (MLR) method in terms of accuracy and precision. Pharmacokinetic study of pioglitazone was employed as the example case. The population pharmacokinetics was characterized by nonlinear mixed effects model (NONMEM). The sparse sampling strategy (1-4 points) was identified by D-optimal algorithm using WinPOPT software. The simulated data generated by Monte Carlo method were used to access the performance of MAPB and MLR. As the number of samples per subject decreased, the accuracy and precision of MAPB method tended to get worse. The estimation for CL and Vby MAPB using D-optimal two-point design had less bias with low inter-individual variability, and had more bias and imprecision with high residue variability. The estimation of AUC by MAPB using D-optimal 2 points design had similar accuracy and precision to MLR. However, MAPB estimation was better than MLR while adjusting the sampling time to one hour. Overall, the MAPB method had similar predictive performance as MLR, but MAPB could provide more pharmacokinetic information with higher sampling flexibility.

Insight on fractal assessment strategies for tin dioxide thin films.

Tin oxide is a unique material of widespread technological applications, particularly in the field of environmental functional materials. New strategies of fractal assessment for tin dioxide thin films formed at different substrate temperatures are of fundamental importance in the development of microdevices, such as gas sensors for the detection of environmental pollutants. Here, tin dioxide thin films with interesting fractal features were successfully prepared by pulsed laser deposition techniques under different substrate temperatures. Fractal method has been first applied to the evaluation of this material. The measurements of carbon monoxide gas sensitivity confirmed that the gas sensing behavior is sensitively dependent on fractal dimensions, fractal densities, and average sizes of the fractal clusters. The random tunneling junction network mechanism was proposed to provide a rational explanation for this gas sensing behavior. The formation process of tin dioxide nanocrystals and fractal clusters could be reasonably described by a novel model.

[Development and identifiability analysis of parent-metabolite pharmacokinetic model for risperidone and its main active metabolite 9-hydroxyrisperidone].

To develop a parent-metabolite pharmacokinetic model for risperidone (RIP) and its major active metabolite (9-hydroxyrisperidone) and investigate their pharmacokinetics characteristics in healthy male volunteers, twenty-two healthy volunteers were orally given a single dose of 2 mg RIP. Plasma samples were collected in the period of 96 hours and concentrations of RIP and 9-hydroxyrisperidone were measured by a validated HPLC/MS method. CYP2D6 phenotypes were identified by the T1/2 of RIP and 9-hydroxyrisperidone according to the literature. Model structure identifiability analysis was performed by the similarity transformation approach to investigate whether the unknown parameters of the proposed model could be estimated from the designed experiment. Pharmacokinetics parameters were estimated using weighted least squares method, and the final pharmacokinetics model were tested and evaluated by Monte Carlo simulation. Eighteen volunteers were phenotyped as extensive metabolizers (EM) and four volunteers were identified as intermediate metabolizers (IM). The final model included central and peripheral compartment for both parent (RIP) and metabolite (9-hydroxyrisperidone) respectively. Model structure identifiability analysis indicated that the proposed model was local identifiable. However, if the ratio of RIP converted to 9-hydroxyrisperidone was assumed to be 32% in EM, and 22% in IM, the model could be globally identifiable. The predicted time-concentration curve and AUC(0-t), C(max), T(max) of RIP and 9-hydroxyrisperidone estimated by the established model were in agreement with the observations and noncompartment analysis. Rate constant of RIP conversion to 9-hydroxyrisperidone was (0.12 +/- 0.08) h(-1) and (0.014 +/- 0.007) h(-1) for EM and IM, respectively. Elimination rate constants of RIP were (0.25 +/- 0.18) and (0.05 +/- 0.23) h(-1) for EM and IM, respectively. Model validation result showed that all parameters derived from the concentration data fitted well with the theoretical value, with mean prediction error of most PK parameter within +/- 15%. The established model well defined the disposition of RIP and 9-hydroxyrisperidone simultaneously and showed large inter-individual pharmacokinetics variation in different CYP2D6 phenotype. The model also provide a useful approach to characterize pharmacokinetics of other parent-metabolite drugs.

[Bioequivalence assessment of pioglitazone hydrochloride oral preparation by limited sampling strategy].

AIM: To develop limited sampling strategy (LSS) for estimation of C(max) and AUC(0-t) and assessing the bioequivalence of two pioglitazone hydrochloride (PGT) preparations. METHODS: Healthy subjects (n = 20), enrolled in a bioequivalence study, were received 30 mg PGT po of reference or test formulation. The plasma concentration of PGT was determined by the validated HPLC method. A multiple linear regression analysis of the Cmax and AUC(0-t) against the PGT concentration for the reference formulation was carried out to develop LSS models to estimate these parameters. The models were internally validated by the Jackknife method and externally validated using simulated sets generated by Monte Carlo method. The best model was employed to assess bioequivalence of the two PGT formulations. RESULTS: The linear relationship between pharmacokinetics parameters and single concentration point was poor. Several models for these parameters estimation met the predefined criteria (r2 > 0.9). The Jackknife validation procedure revealed that LSS models based on two sampling times (C1, C2.5 and C1.5, C2.5 for C(max); C1.5, C9 and C2.5, C9 for AUC(0-t) predict accurately. Mean prediction errors (MPE) were less than 3%, and mean absolute prediction error (MAE) were less than 9%. The prediction error (PE) beyond 20% was less than 5% of total samples. Model external validation by Monte Carlo simulated data indicated that the most informative sampling combinations were C1.5, C2.5 for C(max), and C1.5, C9 for AUC(0-t), respectively. MPE and MAE of the proposed models were less than 5% , and 9% respectively. The PE beyond 20% was less than 5% of the total. Bioequivalence assessment of the two PGT formulations, based on the best LSS models, provided results similar to those obtained using all the observed concentration-time data points, and indicated that the two PGT formulations were bioequivalent. CONCLUSION: The LSS method for bioequivalence assessment of PGT formulations was established and proved to be applicable and accurate. Thus, it could be considered appropriate for PGT bioequivalence study with inexpensive cost of sampling acquisition and analysis. Key words: pioglitazone hydrochloride; limited sampling strategy; Monte Carlo simulation; bioequivalence