Search

Purpose: To describe the pharmacokinetic properties of the [18F]fluoro-polyethylene glycol(PEG)-folate radiotracer in PET/CT imaging of patients with advanced stage epithelial ovarian cancer (EOC). Procedures: In five patients with advanced EOC (FIGO stage IIIB/IIIC, Fédération Internationale de Gynécologie et d'Obstétrique), a 90-min dynamic PET acquisition of the pelvis was performed directly after i.v. administration of 185 MBq [18F]fluoro-PEG6-folate. Arterial blood samples collected at nineteen timepoints were used to determine the plasma input function. A static volume of interest (VOI) for included tumor lesions was drawn manually on the PET images. Modelling was performed using PMOD software. Three different models (a 1-tissue compartment model (1T2k) and two 2-tissue compartment models, irreversible (2T3k) and reversible (2T4k)) were compared in goodness of fit with the time activity curves by means of the Akaike information criterion. Results: The pharmacokinetic analysis in the pelvic area has proven to be much more challenging than expected. Only four out of 22 tumor lesions in five patients were considered suitable to perform modelling on. The remaining tumor lesions were inapt due to either low tracer uptake, small size, proximity to other [18F]fluoro-PEG6-folate -avid structures and/or displacement by abdominal organ motion in the dynamic scan. Data from the four analyzed tumor lesions suggest that the irreversible 2T3k may best describe the pharmacokinetics. All 22 lesions were immunohistochemically stained positive for the folate receptor alpha (FRα) after resection. Conclusion: Performing pharmacokinetic analysis in the abdominal pelvic region is very challenging. This brief article describes the challenges and pitfalls in pharmacokinetic analysis of a tracer with high physiological accumulation in the intestines, in case of lesions of limited size in the abdominal pelvic area. [ABSTRACT FROM AUTHOR]

Aciclovir is considered the first‐line treatment against Herpes simplex virus (HSV) infections in new‐borns and infants. As renal excretion is the major route of elimination, in renally‐impaired patients, aciclovir doses are adjusted according to the degree of impairment. However, limited attention has been given to the implications of immature renal function or dysfunction due to the viral disease itself. The aim of this investigation was to characterize the pharmacokinetics of aciclovir taking into account maturation and disease processes in the neonatal population. Pharmacokinetic data obtained from 2 previously published clinical trials (n = 28) were analyzed using a nonlinear mixed effects modeling approach. Post‐menstrual age (PMA) and creatinine clearance (CLCR) were assessed as descriptors of maturation and renal function. Simulation scenarios were also implemented to illustrate the use of pharmacokinetic data to extrapolate efficacy from adults. Aciclovir pharmacokinetics was described by a one‐compartment model with first‐order elimination. Body weight and diagnosis (systemic infection) were statistically significant covariates on the volume of distribution, whereas body weight, CLCR and PMA had a significant effect on clearance. Median clearance varied from 0.2 to 1.0 L/h in subjects with PMA <34 or ≥34 weeks, respectively. Population estimate for volume of distribution was 1.93 L with systemic infection increasing this value by almost 3‐fold (2.67 times higher). A suitable model parameterization was identified, which discriminates the effects of developmental growth, maturation, and organ function. Exposure to aciclovir was found to increase with decreasing PMA and renal function (CLCR), suggesting different dosing requirement for pre‐term neonates. [ABSTRACT FROM AUTHOR]

PurposeModel-constrained reconstruction with Fourier-based undersampling (MoReFUn) is introduced to accelerate the acquisition of dynamic MRI using hyperpolarized [1-13 C]-pyruvate.MethodsThe MoReFUn method resolves spatial aliasing using constraints introduced by a pharmacokinetic model that describes the signal evolution of both pyruvate and lactate. Acceleration was evaluated on three single-channel data sets: a numerical digital phantom that is used to validate the accuracy of reconstruction and model parameter restoration under various SNR and undersampling ratios, prospectively and retrospectively sampled data of an in vitro dynamic multispectral phantom, and retrospectively undersampled imaging data from a prostate cancer patient to test the fidelity of reconstructed metabolite time series.ResultsAll three data sets showed successful reconstruction using MoReFUn. In simulation and retrospective phantom data, the restored time series of pyruvate and lactate maintained the image details, and the mean square residual error of the accelerated reconstruction increased only slightly (

Introduction: Cannabis is widely used for recreational and medical purposes, but its therapeutic efficacy remains unresolved for many applications as data from retrospective studies show dramatic discrepancy. We hypothesized that false self-reporting of cannabis use and lack of differentiation of heavy users from light or occasional users contribute to the conflicting outcomes. Objective: The goal of this study was to develop an objective biomarker of cannabis use and test how application of such biomarker impacts clinical study outcomes and dose-response measures. Methods and Analysis: Population pharmacokinetic (PK) models of (-)-trans-Δ9-tetrahydrocannabinol (THC) and its metabolites 11-hydroxy-Δ9-tetrahydrocannabinol (11-OH-THC) and 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (11-COOH-THC) were developed based on published studies reporting cannabinoid disposition in individual subjects following intravenous administration or smoking of cannabis. Plasma 11-COOH-THC concentration distributions in different cannabis user groups smoking cannabis were generated via Monte Carlo simulations, and plasma concentration cutoff values of 11-COOH-THC were developed to differentiate light and heavy daily cannabis users in clinical studies. The developed cutoff value was then applied to a retrospective study that assessed the impact of cannabis use on T cell activation in subjects with HIV who self-reported as either nonuser or daily user of cannabis. Results: The developed population PK models established plasma 11-COOH-THC concentration of 73.1 μg/L as a cutoff value to identify heavy daily users, with a positive predictive value of 80% in a mixed population of equal proportions of once daily and three times a day users. The stratification allowed detection of changes in T cell activation in heavy users which was not detected based on self-reporting or detectability of plasma cannabinoids. A proof-of-concept power analysis demonstrated that implementation of such cutoff value greatly increases study power and sensitivity to detect pharmacological effects of cannabis use. Conclusions: This study shows that the use of plasma 11-COOH-THC concentration cutoff value as an objective measure to classify cannabis use in target populations is critical for study sensitivity and specificity and provides much needed clarity for addressing dose-response relationships and therapeutic effects of cannabis.

Drug dose response curves are ubiquitous in cancer biology, but these curves are often used to measure differential response in first-order effects: the effectiveness of increasing the cumulative dose delivered. In contrast, second-order effects (the variance of drug dose) are often ignored. Knowledge of second-order effects may improve the design of chemotherapy scheduling protocols, leading to improvements in tumor response without changing the total dose delivered. By considering treatment schedules with identical cumulative dose delivered, we characterize differential treatment outcomes resulting from high variance schedules (e.g. high dose, low dose) and low variance schedules (constant dose). We extend a previous framework used to quantify second-order effects, known as antifragility theory, to investigate the role of drug pharmacokinetics. Using a simple one-compartment model, we find that high variance schedules are effective for a wide range of cumulative dose values. Next, using a mouse-parameterized two-compartment model of 5-fluorouracil, we show that schedule viability depends on initial tumor volume. Finally, we illustrate the trade-off between tumor response and lean mass preservation. Mathematical modeling indicates that high variance dose schedules provide a potential path forward in mitigating the risk of chemotherapy-associated cachexia by preserving lean mass without sacrificing tumor response. [ABSTRACT FROM AUTHOR]

Background and purpose: In our previous studies, drug nanocrystals were directly prepared by solution crystallization, possessing uniform particle size and morphology suitable for intravenous (IV) injection. These nanocrystals accumulated in a small percentage of their injected dose in tumor-bearing mice but showed similar anti-tumor effectiveness and much-reduced side effects compared with current commercial solubilized and encapsulated delivery systems. Experimental approach: In this study, we aimed to delineate possible controlling factors for the pharmacokinetics (PK) and biodistribution behaviors of paclitaxel (PTX) nanocrystals tested in mice by applying physiologically based pharmacokinetics (PBPK) modeling, coupled with pharmacodynamics (PD) simulation, to the data. Key Results: Our results show that clearance of the drug plays a significant, if not the most important, role in determining tissue distribution, including tumor accumulation of PTX nanocrystals. Surface treatment of drug nanocrystals with polymeric surfactants also appeared to affect PK profiles and PD outcomes. Importantly, when scaled to model human parameters, our PK/PD simulations suggest that drug distribution in humans, as opposed to animal models, was significantly influenced by tissue partitioning rather than drug clearance. This finding could facilitate the design and development of future drug delivery systems. Conclusion: Drug nanocrystals deposited in tissues, including tumors, could therefore act as depots, releasing the drug back into the circulation, possibly contributing to extended treatment, as well as any detrimental effects.

Abstract Aciclovir is considered the first‐line treatment against Herpes simplex virus (HSV) infections in new‐borns and infants. As renal excretion is the major route of elimination, in renally‐impaired patients, aciclovir doses are adjusted according to the degree of impairment. However, limited attention has been given to the implications of immature renal function or dysfunction due to the viral disease itself. The aim of this investigation was to characterize the pharmacokinetics of aciclovir taking into account maturation and disease processes in the neonatal population. Pharmacokinetic data obtained from 2 previously published clinical trials (n = 28) were analyzed using a nonlinear mixed effects modeling approach. Post‐menstrual age (PMA) and creatinine clearance (CLCR) were assessed as descriptors of maturation and renal function. Simulation scenarios were also implemented to illustrate the use of pharmacokinetic data to extrapolate efficacy from adults. Aciclovir pharmacokinetics was described by a one‐compartment model with first‐order elimination. Body weight and diagnosis (systemic infection) were statistically significant covariates on the volume of distribution, whereas body weight, CLCR and PMA had a significant effect on clearance. Median clearance varied from 0.2 to 1.0 L/h in subjects with PMA

Introduction/Justification: Magnetic nanoparticles (MNPs) have been explored as a new potential theranostic agent, and several studies have devoted significant efforts to employ MNPs in biomedicine, including their applications as imaging contrast agents in alternating current biosusceptometry (ACB). In this field, the novelty of multichannel alternating current biosusceptometry (MC-ACB) devices allow the generation of real-time magnetic images of processes of biodistribution of MNPs in essays with animal models, including experiments focused on liver diseases. Objectives: This study refers to the in vivo biodistribution of MNPs detected by the multichannel ACB system, aimed at how the pharmacokinetics of MNPs is affected in the case of hepatocellular carcinoma. This evaluation has already been presented in a previous paper in terms of experimental results, but not in terms of pharmacokinetic modeling. Materials and Methods: In order to quantitatively describe how this disease may alter the biodistribution of MNPs, two different groups of animals are addressed here: a control group of healthy animals (SAL) and a group of animals with hepatocellular carcinoma (DEN/TAA). The MC-ACB system was used to simultaneously record the transit of MNPs in the heart and their accumulation in the liver. Pharmacokinetic rates of change of MNPs are reported here by proceeding with population parameter estimation for two groups of animals: cancer (DEN/TAA) and control (SAL). They refer to the change of MNPs from heart to liver (k1), from liver to heart (k2), and as the irreversible uptake of MNPs by Kupffer cells within a liver subcompartment (k3). All animal experiments were previously approved and performed according to the protocol 7571041120 by the Ethics Committee on Animal Use of the State University of São Paulo (IBB/UNESP). Results: Notably, both k2 and k3 were found to differ between groups, but not k1, consistent with the fact that MNP liver pharmacokinetics are expected to be affected by the chemically induced cirrhosis-associated hepatocarcinogenesis of the DEN/TAA group. The fact that k2 and k3 are higher for the DEN/TAA group is related to earlier MNP liver saturation in cirrhosis due to higher blood volume, and cirrhosis-associated hepatocarcinogenesis is revealed to affect the biodistribution of magnetic nanoparticles. Conclusion: The modeling approach proposed in the research provides a powerful tool to quantitatively describe the biodistribution of magnetic nanoparticles in healthy and cirrhotic rats, allowing the estimation of pharmacokinetic rate parameters related to the distribution of magnetic nanoparticles. The introduced modeling robustly describes the concentration of nanoparticles in compartments over time, which may assist in the development of targeted drug delivery strategies. Finally, this study highlights the relevance of mathematical modeling for understanding complex phenomena such as nanoparticle interactions with living systems, particularly in the development of effective theranostic applications.

BackgroundA knowledge gap exists for dolutegravir (DTG) pharmacokinetics and safety during the first 4 weeks of life, preventing safe and effective DTG use in neonates.SettingPopulation pharmacokinetic modeling and simulation were used to assess newborn DTG dosing requirements during the first few days of life as a function of maternal DTG dosing history before delivery.MethodsDTG PK data were obtained from pregnant women and infants enrolled in the International Maternal Pediatric Adolescent AIDS Clinical Trials Network P1026S study. Maternal and neonate population pharmacokinetic models were separately developed. Monte Carlo simulations were performed to simulate neonatal concentrations after 2 doses of DTG after birth for infants born to mothers either receiving or not receiving DTG before delivery.ResultsIn DTG-naïve infants, a 5-mg DTG dose at birth with a second dose after 48 hours maintained median concentrations above the lower bound of the target range (0.77 μg/mL) and below the upper bound of the target range (7.34 μg/mL representing 2-fold above the adult Cmax value). In DTG-exposed infants, a 5-mg DTG dose at 24 hours after birth with a second dose after 48 hours maintained median concentrations within or nearly within the target range, even if the last maternal DTG dose was taken as soon as 6 hours or as long as 24 hours before delivery.ConclusionsNewborn DTG dosing requirements during the first few days of life depend on maternal DTG dosing history before delivery. These results may help the design of future clinical studies of DTG in the neonatal population.

Background: Current acceptable chemical exposure levels (e.g., tolerable daily intake) are mainly based on animal experiments, which are costly, time-consuming, considered non-ethical by many, and may poorly predict adverse outcomes in humans. Objective: To evaluate a method using human in vitro data and biological modeling to calculate an acceptable exposure level through a case study on 2,2′,4,4′-tetrabromodiphenyl ether (BDE-47) developmental neurotoxicity (DNT). Methods: We reviewed the literature on in vitro assays studying BDE-47-induced DNT. Using the most sensitive endpoint, we derived a point of departure using a mass-balance in vitro disposition model and benchmark dose modeling for a 5% response (BMC05) in cells. We subsequently used a pharmacokinetic model of gestation and lactation to estimate administered equivalent doses leading to four different metrics of child brain concentration (i.e., average prenatal, average postnatal, average overall, and maximum concentration) equal to the point of departure. The administered equivalent doses were translated into tolerable daily intakes using uncertainty factors. Finally, we calculated biomonitoring equivalents for maternal serum and compared them to published epidemiological studies of DNT. Results: We calculated a BMC05 of 164 μg/kg of cells for BDE-47 induced alteration of differentiation in neural progenitor cells. We estimated administered equivalent doses of 0.925–3.767 μg/kg/day in mothers, and tolerable daily intakes of 0.009–0.038 μg/kg/day (composite uncertainty factor: 100). The lowest derived biomonitoring equivalent was 19.75 ng/g lipids, which was consistent with reported median (0.9–23 ng/g lipids) and geometric mean (7.02–26.9 ng/g lipids) maternal serum concentrations from epidemiological studies. Conclusion: This case study supports using in vitro data and biological modeling as a viable alternative to animal testing to derive acceptable exposure levels.

Abstract Albendazole is an anthelmintic drug commonly used in parenchymal neurocysticercosis and cystic echinococcosis. The aim of this study was to explore whether disparities in the dissolution profiles of albendazole products lead to significant differences in pharmacokinetic parameters. Three generic products and the innovator were evaluated in vitro. Quality control tests were performed, and dissolution profiles were obtained according to the Mexican Pharmacopeia. Although all products passed the quality control tests, none of the generic products complied with the similarity factor (f 2). The product with the lowest f 2 value in respect to the reference was chosen for in vivo evaluation. The study was carried out in 12 healthy volunteers who received 400 mg of the generic or reference product according to a crossover design. No significant differences were found in Cmax and AUC for albendazole and its main metabolite, albendazole sulfoxide, between products. Two absorption peaks were observed in the pharmacokinetic profile, and a population (22%) with different absorption rates and delay time for the the second peak was found. Based on the results, due to the high variability in the absorption process the differences observed in vitro could not be observed in vivo.

Objective: Precisely localizing the seizure onset zone remains a challenging task in drug-resistant epilepsy (DRE) patients especially given its critical role in successful surgery and effective management. This study aimed to investigate the kinetic parameters of regional 18F-fluorodeoxyglucose (FDG) uptake in DRE patients, aiming to identify the kinetic parameters best enabling the identification of the epileptogenic region.Methods: Consecutive DRE patients with clinically mandated interictal 18F-FDG PET/CT were recruited from October 2019 to September 2020 for pre-surgical evaluation. Immediately after injecting 18F-FDG of 112–179 MBq, dynamic data were acquired for 90 min. The motion correction and resampling to the Montreal atlas was performed in order to generate a transformation matrix. 116 volume of interests (VOIs) and regional time-activity curves (TACs) were generated by employing the automated anatomical labeling (AAL) template using PMOD software. Kinetic parameters of FDG unidirectional blood-brain clearance (K1), efflux (k2), phosphorylation (k3), and net metabolic flux (Ki) were derived using irreversible 2-tissue-compartment model with an image-derived input function (IDIF). The kinetic parameters values obtained from all regions were ranked and compared with the presumed epileptogenic zone (EZ).Results: Eleven DRE patients (5 males, 6 females, mean age 35.1 ± 10.2 years) were analyzed. We found that the region with the lowest values of Ki provided correct lateralization in 7/7 (100%) of patient with temporal lobe epilepsy (TLE) and the region with the lowest Ki and k3 parameters showed concordance with the EZ in 100% and 71.4% of patients, respectively.Conclusion: The present parametric approach to the evaluation of FDG-PET may be more sensitive than semi-quantitative approaches for the detection of pathophysiology in the EZ of patients with medically unresponsive TLE in addition to the routine clinical investigations.

Background: Perfluorooctanoic acid (PFOA) is detected in the blood of virtually all biomonitoring study participants. Assessing health risks associated with blood PFOA levels is challenging because exposure guidance values (EGVs) are typically expressed in terms of external dose. Biomonitoring equivalents (BEs) consistent with EGVs could facilitate health-based interpretations. Objective: To i) derive BEs for serum/plasma PFOA corresponding to non-cancer EGVs of the U.S. Environmental Protection Agency (U.S. EPA), the Agency for Toxic Substances and Disease Registry (ATSDR) and Health Canada, and ii) compare with PFOA concentrations from national biomonitoring surveys. Methods: Starting from EGV points of departure, we employed pharmacokinetic data/models and uncertainty factors. Points of departure in pregnant rodents (U.S. EPA 2016, ATSDR) were converted into fetus and pup serum concentrations using an animal gestation/lactation pharmacokinetic model, and equivalent human fetus and child concentrations were converted into BEs in maternal serum using a human gestation/lactation model. The point of departure in adult rodents (Health Canada) was converted into a BE using experimental data. For epidemiology-based EGVs (U.S. EPA 2023, draft), BEs were directly based on epidemiological data or derived using a human gestation/lactation pharmacokinetic model. BEs were compared with Canadian/U.S. biomonitoring data. Results: Non-cancer BEs (ng/mL) were 684 (Health Canada, 2018) or ranged from 15 to 29 (U.S. EPA, 2016), 6–10 (ATSDR, 2021) and 0.2–0.8 (U.S. EPA, 2023, draft). Ninety-fifth percentiles of serum levels from the 2018–2019 Canadian Health Measures Survey (CHMS) and the 2017–2018 National Health and Nutrition Examination Survey (NHANES) were slightly below the BE for ATSDR, and geometric means were above the non-cancer BEs for the U.S. EPA (2023, draft). Conclusion: Non-cancer BEs spanned three orders of magnitude. The lowest BEs were for EGVs based on developmental endpoints in epidemiological studies. Concentrations in Canadian/U.S. national surveys were higher than or close to BEs for the most recent non-cancer EGVs.

Golimumab was recently evaluated in a phase 2a, randomized, double‐blind, placebo‐controlled, multicenter study (T1GER study) for safety and efficacy in children and young adults with newly diagnosed type 1 diabetes (T1D). Golimumab showed significant treatment effect where endogenous insulin production was preserved and clinical and metabolic parameters were improved. The objective of this report was to evaluate pharmacokinetic (PK) and pharmacodynamic data from the T1GER study by developing a population pharmacokinetic (PopPK) model and performing exposure‐response (ER) analyses. The PopPK model was developed using data from the T1D study and 2 other pediatric studies with golimumab in ulcerative colitis and in polyarticular juvenile idiopathic arthritis. A 1‐compartment model with first‐order absorption and elimination was applied to describe the concentration‐time profiles. Typical parameters normalized to the values in subjects with a standard weight of 70 kg were apparent clearance, 0.850 L/day; apparent volume of distribution, 16.0 L; and absorption rate constant, 1.01/day. From the ER analyses, no clear trends were observed for changes in both C‐peptide area under the concentration‐time curve and hemoglobin A1c levels for the relatively narrow exposure ranges following the body surface area–based dosing regimen used in this study. In conclusion, the developed PopPK model was able to adequately describe the observed PK of golimumab in patients with T1D. Although golimumab showed significant treatment effect, the ER analyses did not show clear trends within the active treatment group, which may indicate that the exposure from this T1D‐specific dosing regimen was at the plateau of the ER curve. [ABSTRACT FROM AUTHOR]

Purpose: Model‐constrained reconstruction with Fourier‐based undersampling (MoReFUn) is introduced to accelerate the acquisition of dynamic MRI using hyperpolarized [1‐13C]‐pyruvate. Methods: The MoReFUn method resolves spatial aliasing using constraints introduced by a pharmacokinetic model that describes the signal evolution of both pyruvate and lactate. Acceleration was evaluated on three single‐channel data sets: a numerical digital phantom that is used to validate the accuracy of reconstruction and model parameter restoration under various SNR and undersampling ratios, prospectively and retrospectively sampled data of an in vitro dynamic multispectral phantom, and retrospectively undersampled imaging data from a prostate cancer patient to test the fidelity of reconstructed metabolite time series. Results: All three data sets showed successful reconstruction using MoReFUn. In simulation and retrospective phantom data, the restored time series of pyruvate and lactate maintained the image details, and the mean square residual error of the accelerated reconstruction increased only slightly (< 10%) at a reduction factor up to 8. In prostate data, the quantitative estimation of the conversion‐rate constant of pyruvate to lactate was achieved with high accuracy of less than 10% error at a reduction factor of 2 compared with the conversion rate derived from unaccelerated data. Conclusion: The MoReFUn technique can be used as an effective and reliable imaging acceleration method for metabolic imaging using hyperpolarized [1‐13C]‐pyruvate. [ABSTRACT FROM AUTHOR]

Purpose: To propose a robust time and space invariant deep learning (DL) method to directly estimate the pharmacokinetic/tracer kinetic (PK/TK) parameters from undersampled dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) data. Methods: DCE-MRI consists of 4D (3D-spatial + temporal) data and has been utilized to estimate 3D (spatial) tracer kinetic maps.Existing DL architecture for this task needs retraining for variation in temporal and/or spatial dimensions. This work proposes a DL algorithm that is invariant to training and testing in both temporal and spatial dimensions. The proposed network was based on a 2.5-dimensional Unet architecture,where the encoder consists of a 3D convolutional layer and the decoder consists of a 2D convolutional layer. The proposed VTDCE-Net was evaluated for solving the ill-posed inverse problem of directly estimating TK parameters from undersampled k − t space data of breast cancer patients, and the results were systematically compared with a total variation (TV) regularization based direct parameter estimation scheme. In the breast dataset, the training was performed on patients with 32 time samples, and testing was carried out on patients with 26 and 32 time samples. Translation of the proposed VTDCE-Net for brain dataset to show the generalizability was also carried out. Undersampling rates (R) of 8×, 12×, and 20× were utilized with PSNR and SSIM as the figures of merit in this evaluation. TK parameter maps estimated from fully sampled data were utilized as ground truth. Results: Experiments carried out in this work demonstrate that the proposed VTDCE-Net outperforms the TV scheme on both breast and brain datasets across all undersampling rates. For Ktrans and Vp maps, the improvement over TV is as high as 2 and 5 dB, respectively, using the proposed VTDCE-Net. Conclusion: Temporal points invariant DL network that was proposed in this work to estimate the TK-parameters using DCE-MRI data has provided state-of-the-art performance compared to standard image reconstruction methods and is shown to work across all undersampling rates. [ABSTRACT FROM AUTHOR]

Purpose: A method is presented to select the optimal time points at which to measure DCE‐MRI signal intensities, leaving time in the MR exam for high‐spatial resolution image acquisition. Theory: Simplicial complexes are generated from the Kety‐Tofts model pharmacokinetic parameters Ktrans and ve. A geometric search selects optimal time points for accurate estimation of perfusion parameters. Methods: The DCE‐MRI data acquired in women with invasive breast cancer (N = 27) were used to retrospectively compare parameter maps fit to full and subsampled time courses. Simplicial complexes were generated for a fixed range of Kety‐Tofts model parameters and for the parameter ranges weighted by estimates from the fully sampled data. The largest‐area manifolds determined the optimal three time points for each case. Simulations were performed along with retrospectively subsampled data fits. The agreement was computed between the model parameters fit to three points and those fit to all points. Results: The optimal three‐point sample times were from the data‐informed simplicial complex analysis and determined to be 65, 204, and 393 s after arrival of the contrast agent to breast tissue. In the patient data, tumor‐median parameter values fit using all points and the three selected time points agreed with concordance correlation coefficients of 0.97 for Ktrans and 0.67 for ve. Conclusion: It is possible to accurately estimate pharmacokinetic parameters from three properly selected time points inserted into a clinical DCE‐MRI breast exam. This technique can provide guidance on when to capture images for quantitative data between high‐spatial‐resolution DCE‐MRI images. [ABSTRACT FROM AUTHOR]

Thanh D Dinh,1,2 Hung N Nguyen,1 Ba Hai Le,1 Thuy TT Nguyen,1 Huong TL Nguyen1 1Department of Clinical Pharmacy, Hanoi University of Pharmacy, Hanoi, Vietnam; 2Phu Tho General Hospital, Việt Trì, Phu Tho Province, VietnamCorrespondence: Huong TL Nguyen, Department of Clinical Pharmacy, Hanoi University of Pharmacy, 13, 15 – Lê Thánh Tông, Hoàn Kiếm, HàNội, 100000, Vietnam, Tel +84904308406, Email huongntl@hup.edu.vnPurpose: The purpose of this study was to characterize the population-based pharmacokinetic (POP-PK) profile of imipenem in Vietnamese adult patients and to assess the probability of target attainment (PTA) of the pharmacokinetic/pharmacodynamic (PK/PD) parameter to determine the optimal dose.Patients and Methods: A POP-PK model of imipenem was developed in patients with severe infection from a 1500-bed general hospital in Vietnam, using MONOLIX 2019. After the initial dose infusion, 6 blood samples per patient were collected to measure plasma imipenem levels. Eight covariates (eg, age, weight) were investigated to ascertain their influence on imipenem’s PK. Monte Carlo simulations were performed to determine the PTA for the time in which the total steady-state imipenem concentrations remained above the MIC (T>MIC) for 40% and 100% of the dosing interval.Results: The best fit to the PK data was a two-compartment model with inter-individual variability (IIV) in clearance (CL), central volume of distribution (Vc), intercompartmental clearance (Q), and peripheral volume of distribution (Vp). Only creatinine clearance was retained as a covariate on CL in the final model. The typical value of CL and Vc were estimated to be 4.79 L/h and 11.1 L, respectively. The between-subject variability in this population was noted to be high (> 38%, especially for IIV on Q at 110%). Prolonged or continuous infusion demonstrated efficacy (40% T>MIC) against bacteria with a MIC of 4mg/L. To achieve 100% T>MIC or bacteria with MIC> 4 mg/L, however, the number of doses must be increased while maintaining the same daily dose for the 3-hour prolonged infusion regimen.Conclusion: A population pharmacokinetic model of imipenem was developed for Vietnamese adult patients with severe illness. By using Monte Carlo simulation, the appropriate dose has been suggested based on the bacterial MIC value and the targeted PK/PD goal.Keywords: imipenem, pharmacokinetic modeling, PK/PD, Monte Carlo, carbapenem

Nanoparticle (NP) technology holds significant promise to mediate targeted drug delivery to specific organs in the body. Understanding the 3D biodistribution of NPs in heterogeneous environments such as the tumor tissue can provide crucial information on efficacy, safety and potential clinical outcomes. Here we present a novel end-to-end workflow, VIOLA, which makes use of tissue clearing methodology in conjunction with high resolution imaging and advanced 3D image processing to quantify the spatiotemporal 3D biodistribution of fluorescently labeled ACCURIN® NPs. Specifically, we investigate the spatiotemporal biodistribution of NPs in three different murine tumor models (CT26, EMT6, and KPC-GEM) of increasing complexity and translational relevance. We have developed new endpoints to characterize NP biodistribution at multiple length scales. Our observations reveal that the macroscale NP biodistribution is spatially heterogeneous and exhibits a gradient with relatively high accumulation at the tumor periphery that progressively decreases towards the tumor core in all the tumor models. Microscale analysis revealed that NP extravasation from blood vessels increases in a time dependent manner and plateaus at 72 h post injection. Volumetric analysis and pharmacokinetic modeling of NP biodistribution in the vicinity of the blood vessels revealed that the local NP density exhibits a distance dependent spatiotemporal biodistribution which provide insights into the dynamics of NP extravasation in the tumor tissue. Our data represents a comprehensive analysis of NP biodistribution at multiple length scales in different tumor models providing unique insights into their spatiotemporal dynamics. Specifically, our results show that NPs exhibit a dynamic equilibrium with macroscale heterogeneity combined with microscale homogeneity. [Display omitted] • Novel tissue clearing/image analysis pipeline for quantifying NP biodistribution • Comprehensive 3D biodistribution evaluation of PLA-PEG NPs in murine tumor models • Observed 3D heterogeneity in NP and tumor vasculature macroscale distribution • NP extravasation from blood vessels show distinct, time-dependent kinetics • Local NP density around blood vessels exhibit time-dependent equilibration [ABSTRACT FROM AUTHOR]

Although human urinary aniline and 2,6-dimethylaniline were unexpectedly detected in biomonitoring data, little is known about the daily intake doses of aniline and 2,6-dimethylaniline in the living environment or their relation to tolerable daily intake (TDI) values in humans. In the current study, to evaluate the daily oral intake of aniline and 2,6-dimethylaniline in humans, forward and reverse dosimetry was carried out using simplified in silico physiologically based pharmacokinetic (PBPK) modeling established using in vivo experimental pharmacokinetic data. These data were from humanized-liver mice after single oral doses of 100 mg/kg aniline (previously determined) and 116 mg/kg 2,6-dimethylanine (currently investigated). The in vivo elimination rates of 2,6-dimethylaniline from plasma in humanized-liver mice were generally slow compared with those of aniline. Faster in vitro metabolic elimination rates of aniline mediated by liver 9000 × g supernatant fractions from rats than those from humans may suggest the existence of higher first-pass effects in rats than in humanized-liver mice. In silico aniline and 2,6-dimethylaniline concentration curves in human urine after virtual oral administrations were estimated by human PBPK models created with data from humanized-liver mice. Reverse dosimetry analysis using human PBPK models estimated the daily intake of aniline, based on reported human urinary concentrations in biomonitoring data, to be roughly similar to the aniline TDI level. These results suggest that forward and reverse dosimetry using simplified human PBPK models founded on data from humanized-liver mice can be used to evaluate possible higher than expected exposures of aniline and 2,6-dimethylaniline in humans. [ABSTRACT FROM AUTHOR]

Developing long-acting products and formulations for infectious diseases is a nontrivial undertaking that is frequently classified as high risk and low reward by the pharmaceutical industry. The Long-Acting/Extended Release Antiretroviral Research Resource Program (LEAP) was founded in 2015 with the support of the National Institutes of Health to encourage, promote, and accelerate the development of such products. Assessment methodology for any new proposal brought to this group is part of a framework—the LEAP Process—that includes a landscape analysis of what is currently available in the public domain. This is followed by in silico modeling and simulation offered as a service to the relevant scientific community. A variety of preclinical and clinical outcome metrics are applied to each new agent as part of a continuous feedback loop to improve product characteristics. This allows us to catalog knowledge gaps and barriers that can be addressed by engaged stakeholders. Results are communicated in scientific articles, reviews, and position papers. This undertaking serves to de-risk discovery, development, and implementation by bridging the gaps between academic, regulatory, and industrial investigators, and by engaging those in the community who will be the eventual users of these medicines. The LEAP Process has supported formulations now approved for human immunodeficiency virus, as well as products in clinical and preclinical development for tuberculosis and hepatitis viruses B and C. [ABSTRACT FROM AUTHOR]

Background: I-123 meta-iodobenzylguanidine (MIBG) imaging has long been employed to noninvasively assess the integrity of human norepinephrine transporter-1 and, hence, myocardial sympathetic innervation. Positron-emitting F-18 meta-fluorobenzylguanidine (MFBG) has recently been developed for potentially superior quantitative characterization. We assessed the feasibility of MFBG imaging of myocardial sympathetic innervation. Methods: 16 patients were imaged with MFBG PET (30-minute dynamic imaging of chest, followed by 3 whole-body acquisitions between 30 minutes and 4-hour post-injection). Blood kinetics were assessed from multiple samples. Pharmacokinetic modeling with reversible 1- and 2-compartment models was performed. Kinetic rate constants were re-calculated from truncated datasets. All patients underwent concurrent MIBG SPECT. Results: MFBG myocardial uptake was rapid and sustained; the mean standardized uptake value (SUV (mean ± standard deviation)) was 5.1 ± 2.2 and 3.4 ± 1.9 at 1 hour and 3-4-hour post-injection, respectively. The mean K1 and distribution volume (VT) were 1.1 ± 0.6 mL/min/g and 34 ± 22 mL/cm3, respectively. Both were reproducible when re-calculated from truncated 1-hour datasets (Intraclass Correlation Coefficient of 0.99 and 0.91, respectively). Spearman's ϱ = 0.86 between MFBG SUV and VT and 0.80 between MFBG PET-derived VT and MIBG SPECT-derived heart-to-mediastinum activity concentration ratio. Conclusion: MFBG is a promising PET radiotracer for the assessment of myocardial sympathetic innervation. [ABSTRACT FROM AUTHOR]

ABSTRACT: Poultry production is linked with the use of veterinary medicinal products to manage diseases. Ionophore coccidiostats have been permitted for use as feed additives within the European Union (EU) for the prevention of coccidiosis in various species of poultry with except of laying hens. The presence of chemical residues in eggs is a matter of major concern for consumers’ health. Despite such prohibition of use in laying hens, they were identified as the most common non-target poultry species being frequently exposed to these class of coccidiostats. Many factors can influence the presence of residues in eggs. Carryover of these class of coccidiostat feed additives in the feed of laying hens has been identified as the main reason of their occurrence in commercial poultry eggs. The physicochemical properties of individual compounds, the physiology of the laying hen, and the biology of egg formation are believed to govern the residue transfer rate and its distribution between the egg white and yolk compartments. This paper reviews the causes of occurrence of residues of ionophore coccidiostats in eggs within the EU with special emphasis on their disposition kinetics in laying hens, and residue transfer into eggs. Additional effort was made to highlight future modeling perspectives on the potential application of pharmacokinetic modeling in predicting drug residue transfer and its concentration in eggs.

Monoclonal antibodies targeting the programmed cell death protein 1 (PD-1) remain the most prevalent cancer immunotherapy both as a monotherapy and in combination with additional therapies. Despite the extensive success of anti-PD-1 monoclonal antibodies in the clinic, the experimental relationship between binding affinity and functional potency for anti-PD-1 antibodies in vivo has not been reported. Anti-PD-1 antibodies with higher and lower affinity than nivolumab or pembrolizumab are entering the clinic and show varied preclinical efficacy. Here, we explore the role of broad-ranging affinity variation within a single lineage in a syngeneic immunocompetent mouse model. By developing a panel of murine anti-PD-1 antibodies with varying affinity (ranging from KD = 20 pM – 15 nM), we find that there is a threshold affinity required for maximum efficacy at a given dose in the treatment of the MC38 adenocarcinoma model with anti-PD-1 immunotherapy. Physiologically based pharmacokinetic modeling complements interpretation of the experimental results and highlights the direct relationship between dose, affinity, and PD-1 target saturation in the tumor.

Purpose: Contrast‐enhanced ultrasound (CEUS) by injection of microbubbles (MBs) has shown promise as a cost‐effective imaging modality for prostate cancer (PCa) detection. More recently, nanobubbles (NBs) have been proposed as novel ultrasound contrast agents. Unlike MBs, which are intravascular ultrasound contrast agents, the smaller diameter of NBs allows them to cross the vessel wall and target specific receptors on cancer cells such as the prostate‐specific membrane antigen (PSMA). It has been demonstrated that PSMA‐targeted NBs can bind to the receptors of PCa cells and show a prolonged retention effect in dual‐tumor mice models. However, the analysis of the prolonged retention effect has so far been limited to qualitative or semi‐quantitative approaches. Methods: This work introduces two pharmacokinetics models for quantitative analysis of time–intensity curves (TICs) obtained from the CEUS loops. The first model is based on describing the vascular input by the modified local density random walk (mLDRW) model and independently interprets TICs from each tumor lesion. Differently, the second model is based on the reference‐tissue model, previously proposed in the context of nuclear imaging, and describes the binding kinetics of an indicator in a target tissue by using a reference tissue where binding does not occur. Results: Our results show that four estimated parameters, β, β/λ$\beta /\lambda $, β+/β−${\beta }_ + /{\beta }_ - $, for the mLDRW‐input model, and γ for the reference‐based model, were significantly different (p‐value <0.05) between free NBs and PSMA‐NBs. These parameters estimated by the two models demonstrate different behaviors between PSMA‐targeted and free NBs. Conclusions: These promising results encourage further quantitative analysis of targeted NBs for improved cancer diagnostics and characterization. [ABSTRACT FROM AUTHOR]

Patients with tuberculosis (TB) coinfected with HIV are more likely to have low blood concentrations of the first‐line anti‐TB drugs (associated with poor outcomes). Therapeutic drug monitoring (TDM) is recommended for certain patient populations with TB at increased risk for a poor outcome. Our objective was to estimate the diagnostic accuracy of a 2‐hour TDM serum sample for the first‐line anti‐TB drugs among patients with HIV/TB and evaluate the information gained by an additional 6‐hour sample. We created a virtual (n = 1000) HIV/TB patient population and performed pharmacokinetic simulations using published population models for isoniazid, rifampin, pyrazinamide, and ethambutol. We performed receiver operating characteristic analysis to compare the diagnostic performance of a single 2‐hour serum sample with samples obtained at 2 and 6 hours after dosing. The sensitivity of a single 2‐hour serum concentration to identify patients with HIV/TB with adequate serum exposures was lowest for rifampin (54.9%; 95%CI, 50.79%‐59.41%) and highest for ethambutol (70.8%; 95%CI, 66.06%‐72.61%) for maximum concentration (Cmax) targets. Diagnostic accuracy of a single 2‐hour serum sample for the area under the concentration‐time curve (AUC) from time 0 to 24 hours target was highest for isoniazid (93%; 95%CI, 90.9%‐94.1%) and lowest for pyrazinamide (66.3%; 95%CI, 62.6%‐70.0%). In summary, the diagnostic performance of TDM for Cmax and AUC from time 0 to 24 hours targets demonstrated variability across the first‐line anti‐TB drugs. The addition of a 6‐hour serum sample led to the highest statistically significant improvement (P <.001) and highest increase in diagnostic accuracy (area under the receiver operating characteristic curve) for rifampin for Cmax and AUC. The other first‐line drugs had modest/negligible increases in diagnostic accuracy. [ABSTRACT FROM AUTHOR]

The placenta acts as a barrier, excluding noxious substances whilst actively transferring nutrients to the fetus, mediated by various transporters. This study quantified the expression of key placental transporters in term human placenta (n=5) and BeWo, BeWo b30, and JEG-3 placenta cell lines. Combining these results with pregnancy physiologically-based pharmacokinetic (PBPK) modeling, we demonstrate the utility of proteomic analysis for predicting placental drug disposition and fetal exposure. Using targeted proteomics with QconCAT standards, we found significant expression of P-gp, BCRP, MRP2, MRP4, and MRP6 in the human placenta (0.05 - 0.25 pmol/mg membrane protein) with only regional differences observed for P-gp. Unexpectedly, both P-gp and BCRP were below the limit of quantification in the regularly used BeWo cells, indicating that this cell line may not be suitable for the study of placental P-gp and BCRP-mediated transport. In cellular and vesicular overexpression systems, P-gp and BCRP were detectable as expected. Vesicle batches showed consistent P-gp expression correlating with functional activity (N-methyl-quinidine (NMQ) transport). However, BCRP activity (Estrone 3-sulfate (E1S) transport) did not consistently align with expression levels. Incorporating in vitro transporter kinetic data, along with placental transporter abundance, into a PBPK model enabled the evaluation of fetal exposure. Simulation with a hypothetical drug indicated that estimating fetal exposure relies on the intrinsic clearances of relevant transporters. To minimize interlaboratory discrepancies, expression data was generated using consistent proteomic methodologies in the same lab. Integration of this data in pregnancy-PBPK modeling offers a promising tool to investigate maternal, placental and fetal drug exposure. Significance Statement This study quantified the expression of key transporters in human placenta and various placental cell lines, revealing significant expression variations. By integrating these data with PBPK modeling, the study highlights the importance of transporter abundance data in understanding and predicting placental drug disposition., (© 2024 The Authors. This is an open access article under the terms of the Creative Commons Attribution CC BY License, which permits use, distribution, and reproduction in any medium, provided that the original work is properly cited.)

Abstract Background Caffeine (0.1 g) is used as a central nervous system stimulant and as a nontoxic phenotyping probe for cytochrome P450 1A2. However, an increasing number of suicide attempts by caffeine overdose have been recently reported. Case presentation A 25-year-old woman (body weight, 43 kg) who intentionally took an overdose of 5.9 g caffeine as a suicide attempt was emergently admitted to Kyoto Medical Center. The plasma concentrations of caffeine and its primary metabolite, N-demethylated paraxanthine, in the current case were 100 and 7.3 μg/mL, 81 and 9.9 μg/mL, 63 and 12 μg/mL, and 21 and 14 μg/mL, at 12, 20, 30, and 56 h after oral overdose, respectively. The observed apparent terminal elimination half-life of caffeine during days 1 and 2 of hospitalization was 27 h, which is several times longer than the reported normal value. This finding implied nonlinearity of caffeine pharmacokinetics over such a wide dose range, which could affect the accuracy of values simulated by a simplified physiologically based pharmacokinetic model founded on a normal dose of 100 mg. Low serum potassium levels (2.9 and 3.5 mM) on days 1 and 2 may have been caused by the caffeine overdose in the current case. Conclusions The patient underwent infusion with bicarbonate Ringer’s solution and potassium chloride and was discharged on the third day of hospitalization despite taking a potentially lethal dose of caffeine. The virtual plasma exposures of caffeine estimated using the current simplified PBPK model were higher than the measured values. The present results based on drug monitoring data and additional pharmacokinetic predictions could serve as a useful guide in cases of caffeine overdose.

Abstract Background Loxoprofen is a propionic acid derivative and is the most widely prescribed non-steroidal anti-inflammatory drug in Japan. Loxoprofen is generally considered to be relatively nontoxic. Case presentation A 33-year-old man (body weight, 55 kg) who intentionally took an overdose of 100 tablets of loxoprofen (6000 mg) as a suicide attempt was emergently admitted to Kyoto Medical Center. On arrival, the patient was suffering disorders of consciousness. His plasma concentrations of loxoprofen and its reduced trans-alcohol metabolite were 52 and 24 μg/mL, 3.7 and 2.3 μg/mL, 0.81 and 0.54 μg/mL, and 0.015 and 0.011 μg/mL, respectively, at 4, 26, 50, and 121 h after the oral overdose. The observed apparent terminal elimination half-life of loxoprofen during days 1 and 2 of hospitalization was in the range 6–12 h, which is several times longer than the reported normal value. This finding implied nonlinearity of loxoprofen pharmacokinetics over the current 100-fold dose range, which could affect the accuracy of values simulated by a simplified physiologically based pharmacokinetic (PBPK) model founded on data from a normal dose of 60 mg. The reasons for the delayed eliminations from plasma of loxoprofen and its trans-alcohol metabolite in this case are uncertain, but slight renal impairment (low eGFR values) developed on the second and third hospital days and could be a causal factor. Conclusions Because the patient’s level of consciousness had gradually improved, he was discharged on the fourth day of hospitalization. The virtual plasma exposures of loxoprofen and its reduced trans-alcohol metabolite estimated using the current simplified PBPK model were lower than the measured values in the overdose case. The present results based on drug monitoring data and pharmacokinetic predictions could serve as a useful guide in cases of loxoprofen overdose.

Alan Morrison,1 Melissa E Stauffer,1,2 Anna S Kaufman1 1ScribCo, Albrightsville, PA, USA; 2Scientific Editing Solutions, Walworth, WI, USACorrespondence: Alan MorrisonScribCo, 111 Gower Road, Albrightsville, PA, 18210, USATel +1 610 389-0843Email scribcopersonnel@gmail.comIntroduction: Patients rarely, if ever, take their medications exactly as prescribed. The extent to which missed doses interfere with a drug’s therapeutic effect remains unclear.Methods: After weekly oral dosing of methotrexate (MTX) for rheumatoid arthritis, its polyglutamate derivatives (MTXglu) accumulate in red blood cells, where they are markers for the drug’s therapeutic effectiveness. We used Medication Event Monitoring System data and pharmacokinetic modeling to analyze whether missing MTX doses causes the MTXglu level in red blood cells to fall below the range associated with the drug’s clinical effect.Results: For patients initiating oral MTX, the threshold for clinical effectiveness and the steady state level were reached in medians of 6 weeks and 22 weeks, respectively. For patients at steady state who discontinued MTX, the MTXglu level fell below the therapeutic threshold after a median of 3 weeks. After initiating MTX, single missed doses did not cause a loss of therapeutic effect in the median patient if they occurred after 10 weeks, while runs of ≥ 3 consecutive missed doses did cause the MTXglu level to fall below the therapeutic threshold.Conclusion: While there is considerable variation between patients, pharmacokinetic modeling indicates that instances of isolated single missed doses of MTX typically will not cause polyglutamated methotrexate levels in red blood cells to fall below the range associated with the therapeutic effect. Runs of ≥ 3 consecutive missed doses, however, are typically expected to result in a loss of the therapeutic effect.Keywords: therapeutic effect, adherence, pharmacokinetic modeling, antirheumatic agents, rheumatoid arthritis, drug therapy

Pharmacokinetic (PK) modeling is a useful method for investigating drug absorption, distribution, metabolism, and excretion. The most commonly used mathematical models in PK modeling are the compartment model and physiologically based pharmacokinetic (PBPK) model. Although the theoretical characteristics of each model are well known, there have been few comparative studies of the compatibility of the models. Therefore, we evaluated the compatibility of PBPK and compartment models using the lumping method with 20 model compounds. The PBPK model was theoretically reduced to the lumped model using the principle of grouping tissues and organs that show similar kinetic behaviors. The area under the concentration-time curve (AUC) based on the simulated concentration and PK parameters (drug clearance [CL], central volume of distribution [Vc], peripheral volume of distribution [Vp]) in each model were compared, assuming administration to humans. The AUC and PK parameters in the PBPK model were similar to those in the lumped model within the 2-fold range for 17 of 20 model compounds (85%). In addition, the relationship of the calculated Vd/fu (volume of distribution [Vd], drug-unbound fraction [fu]) and the accuracy of AUC between the lumped model and compartment model confirmed their compatibility. Accordingly, the compatibility between PBPK and compartment models was confirmed by the lumping method. This method can be applied depending on the requirement of compatibility between the two models. [ABSTRACT FROM AUTHOR]

Background: Patients with adrenal insufficiency are treated with oral hydrocortisone (HC) to compensate for the loss of endogenous cortisol production. Intrinsic imperfections of cortisol replacement strategies in mimicking normal cortisol secretion are the underlying cause of the increased morbidity and mortality of patients suffering from secondary adrenal insufficiency (SAI). To improve oral hydrocortisone substitution therapy, a better understanding of its pharmacokinetics (PK) is necessary. The previous PK model did not include protein binding. It is known that protein binding can impact hydrocortisone pharmacokinetics. The aim of this study is to describe HC pharmacokinetics including the protein-binding state using Edsim++ (Mediware, Prague) pharmacokinetic modeling software, paving the way for an in-silico tool suitable for drug delivery design. Methods: A total of 46 patients with SAI participated in a randomized double-blind crossover study Patients randomly received a low dose of HC (0.2–0.3 mg/kg body weight/day) for 10 weeks, followed by a high dose (0.4–0.6 mg/kg body weight/day) for another 10 weeks, or vice versa. Plasma samples were obtained and analyzed for free and total hydrocortisone. Single compartment population pharmacokinetic analysis was performed using an extended Werumeus-Buning model built in Edsim++. This model includes a mathematical approach for estimating free cortisol by Nguyen et al., taking the protein binding of HC to albumin and hydrocortisone-binding globulin (CBG, transcortin) into consideration, as well as different states of CBG which affect binding kinetics to HC. The goodness of fit for observed versus predicted values was calculated. Results and conclusions: Nguyen's formula for free cortisol estimation was successfully implemented in a pharmacokinetic model. The model shows high Spearman's correlation for observed versus predicted hydrocortisone concentrations. Significantly higher correlations (Spearman's r, 0.901 vs. 0.836) between total and free hydrocortisone AUC24 (area-under the curve over 24 h) are found when comparing new and old models. This new model was used to simulate the plasma concentration–time behavior of a more suitable hydrocortisone formulation. [ABSTRACT FROM AUTHOR]

Fexuprazan is a potassium-competitive acid blocker (P-CAB). The compounds in this newly developed drug family suppress intragastric acidity. As there are already other acid-suppressing drugs on the market, such as H2 antagonists and proton pump inhibitors (PPIs), it would be informative to compare the biological effects of fexuprazan against another approved drug with the same indication. The drug concentration predicted by the pharmacokinetic (PK) model could serve as an input function for a pharmacodynamic (PD) model. The apparent pharmacokinetics of fexuprazan could be described by a simpler model. However, a physiologically based pharmacokinetic (PBPK) model was developed in a previous study. A one-compartment model was also proposed in the present study. Both the newly suggested model and the previously validated PBPK model were used as input functions of the PD models. Our simulation revealed that the effects of fexuprazan could be effectively simulated by the proposed PK–PD models. A PK–PD model was also proposed for the oral administration of the PPI reference drug esomeprazole. A model-based analysis was then performed for intragastric pH using several dosing methods. The expected pH could be predicted for both drugs under several dosing regimens using the proposed PK–PD models. [ABSTRACT FROM AUTHOR]

Objective: Ultra-rapid insulin formulations control postprandial hyperglycemia; however, inadequate understanding of injection site absorption mechanisms is limiting further advancement. We used photoacoustic imaging to investigate the injection site dynamics of dye-labeled insulin lispro in the Humalog® and Lyumjev® formulations using the murine ear cutaneous model and correlated it with results from unlabeled insulin lispro in pig subcutaneous injection model. Methods: We employed dual-wavelength optical-resolution photoacoustic microscopy to study the absorption and diffusion of the near-infrared dye-labeled insulin lispro in the Humalog and Lyumjev formulations in mouse ears. We mathematically modeled the experimental data to calculate the absorption rate constants and diffusion coefficients. We studied the pharmacokinetics of the unlabeled insulin lispro in both the Humalog and Lyumjev formulations as well as a formulation lacking both the zinc and phenolic preservative in pigs. The association state of insulin lispro in each of the formulations was characterized using SV-AUC and NMR spectroscopy. Results: Through experiments using murine and swine models, we show that the hexamer dissociation rate of insulin lispro is not the absorption rate-limiting step. We demonstrated that the excipients in the Lyumjev formulation produce local tissue expansion and speed both insulin diffusion and microvascular absorption. We also show that the diffusion of insulin lispro at the injection site drives its initial absorption; however, the rate at which the insulin lispro crosses the blood vessels is its overall absorption rate-limiting step. Conclusions: This study provides insights into injection site dynamics of insulin lispro and the impact of formulation excipients. It also demonstrates photoacoustic microscopy as a promising tool for studying protein therapeutics. The results from this study address critical questions around the subcutaneous behavior of insulin lispro and the formulation excipients, which could be useful to make faster and better controlled insulin formulations in the future.

Abstract Background Exhaled propofol concentrations correlate with propofol concentrations in adult human blood and the brain tissue of rats, as well as with electroencephalography (EEG) based indices of anesthetic depth. The pharmacokinetics of propofol are however different in children compared to adults. The value of exhaled propofol measurements in pediatric anesthesia has not yet been investigated. Breathing system filters and breathing circuits can also interfere with the measurements. In this study, we investigated correlations between exhaled propofol (exP) concentrations and the Narkotrend Index (NI) as well as calculated propofol plasma concentrations. Methods A multi-capillary-column (MCC) combined with ion mobility spectrometry (IMS) was used to determine exP. Optimal positioning of breathing system filters (near-patient or patient-distant) and sample line (proximal or distal to filter) were investigated. Measurements were taken during induction (I), maintenance (M) and emergence (E) of children under total intravenous anesthesia (TIVA). Correlations between ExP concentrations and NI and predicted plasma propofol concentrations (using pediatric pharmacokinetic models Kataria and Paedfusor) were assessed using Pearson correlation and regression analysis. Results Near-patient positioning of breathing system filters led to continuously rising exP values when exP was measured proximal to the filters, and lower concentrations when exP was measured distal to the filters. The breathing system filters were therefore subsequently attached between the breathing system tubes and the inspiratory and expiratory limbs of the anesthetic machine. ExP concentrations significantly correlated with NI and propofol concentrations predicted by pharmacokinetic models during induction and maintenance of anesthesia. During emergence, exP significantly correlated with predicted propofol concentrations, but not with NI. Conclusion In this study, we demonstrated that exP correlates with calculated propofol concentrations and NI during induction and maintenance in pediatric patients. However, the correlations are highly variable and there are substantial obstacles: Without patient proximal placement of filters, the breathing circuit tubing must be changed after each patient, and furthermore, during ventilation, a considerable additional loss of heat and moisture can occur. Adhesion of propofol to plastic parts (endotracheal tube, breathing circle) may especially be problematic during emergence. Trial Registration The study was registered in the German registry of clinical studies (DRKS-ID: DRKS00015795 ).

In brain tumor, dynamic contrast‐enhanced magnetic resonance imaging (DCE‐MRI) spatiotemporally resolved high‐quality reconstruction which is required for quantitative analysis of some physiological characteristics of brain tissue. By exploiting some kind of sparsity priori, compressed sensing methods can achieve high spatiotemporal DCE‐MRI image reconstruction from undersampled k‐space data. Recently, as a kind of priori information about the contrast agent (CA) concentration dynamics, Pharmacokinetic (PK) models have been explored for undersampled DCE‐MRI reconstruction. This paper presents a novel dictionary learning‐based reconstruction method with Adaptive Pharmaco‐Kinetic Model Constraints (APKMC). In APKMC, the priori knowledge about CA dynamics is incorporated into a novel dictionary, which consists of PK model‐based atoms and adaptive atoms. The PK atoms are constructed based on Patlak model and K‐SVD dimension reduction algorithm, and the adaptive ones are used to resolve PK model inconsistencies. To solve APKMC, an optimization algorithm based on variable splitting and alternating iterative optimization is presented. The proposed method has been validated on three brain tumor DCE‐MRI data sets by comparing with two state‐of‐the‐art methods. As demonstrated by the quantitative and qualitative analysis of results, APKMC achieved substantially better quality in the reconstruction of brain DCE‐MRI images, as well as in the reconstruction of PK model parameter maps. [ABSTRACT FROM AUTHOR]

Abstract Background Combinations of antidepressant duloxetine (at doses of 40–60 mg/day) and other antipsychotics are frequently used in clinical treatment; however, several fatal and nonfatal cases of duloxetine overdose have been documented. We experienced a patient who had taken an overdose of duloxetine (780 mg) in combination with other drugs in a suicide attempt. Case presentation The patient was a 37-year-old man (body weight, 64 kg) with a history of gender identity disorder and depression. He intentionally took an overdose of duloxetine in combination with three other antipsychotic drugs (18 mg flunitrazepam, 850 mg quetiapine, and 1100 mg trazodone) and was emergently admitted to Kyoto Medical Center. The patient’s plasma concentration of duloxetine during ambulance transport was 57 ng/ml, and the level was still as high as 126 ng/mL at 32 h after administration. Duloxetine disappeared most slowly from plasma, in contrast to quetiapine, which was the fastest to clear among the four medicines determined in this patient. The observed concentrations of duloxetine in this overdose patient were generally within the 95% confidence intervals of the plasma concentration curves predicted using a physiologically based pharmacokinetic (PBPK) model. Conclusion Even if more than 1 h (the generally recommended period) has passed after administration of duloxetine in such overdose cases, gastric lavage and/or administration of activated charcoal may be effective in clinical practice up to 6 h because of the typically slow elimination behavior illustrated by the PBPK model. Pharmacokinetic profiles visualized using PBPK modeling can inform treatment decisions in cases of drug overdose for medicines such as duloxetine in emergency clinical practice.

The effect of transporters and enzymes on drug pharmacokinetics is increasingly evaluated using genetically modified animals that have these proteins either knocked-out or their human orthologues transgenically expressed. Analysis of pharmacokinetic data obtained in such experiments is typically performed using non-compartmental analysis (NCA), which has limitations such as not being able to identify the PK parameter that is affected by the genetic modification of the enzymes or transporters and the requirement of intense and homogeneous sampling of all subjects. Here we used a compartmental population pharmacokinetic modeling approach using PK data from a series of genetically modified mouse experiments with lorlatinib to extend the results and conclusions from previously reported NCA analyses. A compartmental population pharmacokinetic model was built and physiologically plausible covariates were evaluated for the different mouse strains. With the model, similar effects of the strains on the area under the concentration-time curve (AUC) from 0 to 8 hours were found as for the NCA. Additionally, the differences in AUC between the strains were explained by specific effects on clearance and bioavailability for the strain with human expressing CYP3A4. Finally, effects of multidrug efflux transporters ATP-binding cassette (ABC) sub-family B member 1 (ABCB1) and G member 2 (ABCG2) on brain efflux were quantified. Use of compartmental population PK modeling yielded additional insight into the role of drug-metabolizing enzymes and drug transporters in mouse experiments compared to the NCA. Furthermore, these models allowed analysis of heterogeneous pooled datasets and the sparse organ concentration data in contrast to classical NCA analyses. [ABSTRACT FROM AUTHOR]

Modeling of metabolite kinetics after oral administration of ketamine is of special interest because of the higher concentrations of active metabolites because of the hepatic first‐pass effect. This holds especially in view of the potential analgesic and antidepressant effects of 2R,6R‐ and 2S,6S‐hydroxynorketamine at low doses of ketamine. Therefore, a 9‐compartment model was developed to analyze the pharmacokinetics of ketamine enantiomers and their metabolites after racemic ketamine administered intravenously (5 mg) and as 4 doses (10, 20, 40, and 80 mg) of a prolonged‐release formulation (PR‐ketamine). Using a population approach, the serum concentration‐time data of the enantiomers of ketamine, norketamine, dehydronorketamine, and 2,6‐hydroxynorketamine obtained in 15 healthy volunteers could be adequately fitted. The estimated model parameters were used to simulate serum concentration‐time profiles; after multiple dosing of PR‐ketamine (2 daily doses of 20 mg), the steady‐state concentrations of R‐ and S‐ketamine were 1.4 and 1.3 ng/mL, respectively. The steady‐state concentration of 2R,6R‐hydroxynorketamine exceeded those of R‐norketamine (4‐fold), R‐dehydonorketamine (8‐fold), and R‐ketamine (46‐fold), whereas that of 2S,6S‐hydroxynorketamine exceeded that of S‐ketamine by 14‐fold. The model may be useful for identifying dosing regimens aiming at optimal plasma concentrations of 2,6‐hydroxynorketamines. [ABSTRACT FROM AUTHOR]

Purpose: Pharmacokinetic modeling can be applied to quantify the kinetics of fluorescently labeled compounds using longitudinal micro-computed tomography and fluorescence-mediated tomography (μCT-FMT). However, fluorescence blurring from neighboring organs or tissues and the vasculature within tissues impede the accuracy in the estimation of kinetic parameters. Contributions of elimination and retention activities of fluorescent probes inside the kidneys and liver can be hard to distinguish by a kinetic model. This study proposes a deconvolution approach using a mixing matrix to model fluorescence contributions to improve whole-body pharmacokinetic modeling. Procedures: In the kinetic model, a mixing matrix was applied to unmix the fluorescence blurring from neighboring tissues and blood vessels and unmix the fluorescence contributions of elimination and retention in the kidney and liver compartments. Accordingly, the kinetic parameters of the hepatobiliary and renal elimination routes and five major retention sites (the kidneys, liver, bone, spleen, and lung) were investigated in simulations and in an in vivo study. In the latter, the pharmacokinetics of four fluorescently labeled compounds (indocyanine green (ICG), HITC-iodide-microbubbles (MB), Cy7-nanogels (NG), and OsteoSense 750 EX (OS)) were evaluated in BALB/c nude mice. Results: In the simulations, the corrected modeling resulted in lower relative errors and stronger linear relationships (slopes close to 1) between the estimated and simulated parameters, compared to the uncorrected modeling. For the in vivo study, MB and NG showed significantly higher hepatic retention rates (P<0.05 and P<0.05, respectively), while OS had smaller renal and hepatic retention rates (P<0.01 and P<0.01, respectively). Additionally, the bone retention rate of OS was significantly higher (P<0.01). Conclusions: The mixing matrix correction improves pharmacokinetic modeling and thus enables a more accurate assessment of the biodistribution of fluorescently labeled pharmaceuticals by μCT-FMT. [ABSTRACT FROM AUTHOR]

Purpose: Crizotinib, a potent oral tyrosine kinase inhibitor, was evaluated in combination with dasatinib in a phase 1 trial (NCT01644773) in children with progressive or recurrent high-grade and diffuse intrinsic pontine gliomas (HGG and DIPG). This study aimed to characterize the pharmacokinetics of crizotinib in this population and identify significant covariates.Methods: Patients (N = 36, age range 2.9-21.3 years) were treated orally once or twice-daily with 100-215 mg/m2 crizotinib and 50-65 mg/m2 dasatinib. Pharmacokinetic studies were performed for crizotinib alone after the first dose and at steady state, and for the drug combination at steady state. Crizotinib plasma concentrations were measured using a validated LC-MS/MS method. Population modeling was performed (Monolix) and the impact of factors including patient demographics and co-medications were investigated on crizotinib pharmacokinetics.Results: Crizotinib concentrations were described with a linear two-compartment model and absorption lag time. Concomitant dasatinib and overweight/obese status significantly influenced crizotinib pharmacokinetics, resulting in clinically relevant impact (> 20%) on drug exposure. Crizotinib mean apparent clearance (CL/F) was 66.7 L/h/m2 after single-dose and decreased to 26.5 L/h/m2 at steady state when given alone, but not when combined with dasatinib (mean 60.8 L/h/m2). Overweight/obese patients exhibited lower crizotinib CL/F and apparent volume V1/F (mean 46.2 L/h/m2 and 73.3 L/m2) compared to other patients (mean 75.5 L/h/m2 and 119.3 L/m2, p < 0.001).Conclusion: A potential pharmacokinetic interaction was observed between crizotinib and dasatinib in children with HGG and DIPG. Further, crizotinib exposure was significantly higher in overweight/obese patients, who may require a dosing adjustment. [ABSTRACT FROM AUTHOR]

Abstract Background Circumventricular organs (CVOs) are small structures without a blood–brain barrier surrounding the brain ventricles that serve homeostasic functions and facilitate communication between the blood, cerebrospinal fluid and brain. Secretory CVOs release peptides and sensory CVOs regulate signal transmission. However, pathogens may enter the brain through the CVOs and trigger neuroinflammation and neurodegeneration. We investigated the feasibility of dynamic contrast-enhanced (DCE) MRI to assess the CVO permeability characteristics in vivo, and expected significant contrast uptake in these regions, due to blood–brain barrier absence. Methods Twenty healthy, middle-aged to older males underwent brain DCE MRI. Pharmacokinetic modeling was applied to contrast concentration time-courses of CVOs, and in reference to white and gray matter. We investigated whether a significant and positive transfer from blood to brain could be measured in the CVOs, and whether this differed between secretory and sensory CVOs or from normal-appearing brain matter. Results In both the secretory and sensory CVOs, the transfer constants were significantly positive, and all secretory CVOs had significantly higher transfer than each sensory CVO. The transfer constants in both the secretory and sensory CVOs were higher than in the white and gray matter. Conclusions Current measurements confirm the often-held assumption of highly permeable CVOs, of which the secretory types have the strongest blood-to-brain transfer. The current study suggests that DCE MRI could be a promising technique to further assess the function of the CVOs and how pathogens can potentially enter the brain via these structures. Trial registration: Netherlands Trial Register number: NL6358, date of registration: 2017-03-24

Abstract Background The anticoagulant edoxaban is used clinically at doses of 30–60 mg/day; however, we experienced a patient who had taken an overdose of edoxaban of 750 mg. We investigated the pharmacokinetics of edoxaban in this patient by using liquid chromatography–tandem spectrometry to estimate the follow-up period in emergency clinical practice with this medicine. Case presentation The patient was a 57-year-old woman (body weight, 69 kg) who had taken a single oral dose of 750 mg of edoxaban in a suicide attempt. She was emergently admitted to Kyoto Medical Center. The patient’s edoxaban plasma concentrations during ambulance transport (8 h after oral administration) were ~ 4900 ng/ml, and the concentration gradually decreased to ~ 10 ng/mL and to detectable but unmeasurable levels of ~ 1.0 ng/mL at 60 h and 100 h, respectively. The linear range of the relationship between the dose and plasma concentration was assumed to have been exceeded during the first 8 h; however, the measured elimination rate of edoxaban was similar to that visualized curves predicted by a simplified physiologically based pharmacokinetic model previously established. Conclusion Simplified physiologically based pharmacokinetic models for creating visualized curves have proven to be useful not only during drug discovery or chemical risk assessment but also in cases of medical poisoning. We used a physiologically based pharmacokinetic model previously established for edoxaban to predict the pharmacokinetics in the current case. It is hoped that the results of this study, which encompass drug monitoring data in the patient and visualized pharmacokinetic prediction, will serve as an index when setting the treatment and follow-up period in cases of drug overdose for medicines such as edoxaban in emergency clinical practice.

The single-injection start regimen for aripiprazole once-monthly 400 mg (AOM 400) in patients with schizophrenia requires a single intramuscular injection in the gluteal or deltoid site and 14 days of concurrent oral therapy. A simplified, single-day regimen of two injections at separate gluteal and/or deltoid injection sites, together with a single 20-mg dose of oral aripiprazole on the 1st day, was assessed. A previously developed population-pharmacokinetic (popPK) model for characterizing aripiprazole PK following oral administration and gluteal intramuscular depot injection was expanded to include deltoid injection. Simulations were conducted to assess PK profiles following various (including two-injection) start regimens. Postmarketing data on patients who received higher-than-recommended AOM doses were used to assess overall safety/tolerability. The two-injection start regimen with a single concurrent oral dose displayed a comparable PK profile to the single-injection start regimen with concurrent 14-day oral administration in simulations. The safety assessment indicated the two-injection start regimen was unlikely to be associated with safety concerns beyond those expected with a single-injection start regimen. These data support use of the two-injection start regimen in clinical practice to reduce reliance on daily oral administration and optimize the therapeutic benefits of AOM 400 in patients with schizophrenia. [ABSTRACT FROM AUTHOR]

Background: Caffeine (0.1 g) is used as a central nervous system stimulant and as a nontoxic phenotyping probe for cytochrome P450 1A2. However, an increasing number of suicide attempts by caffeine overdose have been recently reported. Case presentation: A 25-year-old woman (body weight, 43 kg) who intentionally took an overdose of 5.9 g caffeine as a suicide attempt was emergently admitted to Kyoto Medical Center. The plasma concentrations of caffeine and its primary metabolite, N-demethylated paraxanthine, in the current case were 100 and 7.3 μg/mL, 81 and 9.9 μg/mL, 63 and 12 μg/mL, and 21 and 14 μg/mL, at 12, 20, 30, and 56 h after oral overdose, respectively. The observed apparent terminal elimination half-life of caffeine during days 1 and 2 of hospitalization was 27 h, which is several times longer than the reported normal value. This finding implied nonlinearity of caffeine pharmacokinetics over such a wide dose range, which could affect the accuracy of values simulated by a simplified physiologically based pharmacokinetic model founded on a normal dose of 100 mg. Low serum potassium levels (2.9 and 3.5 mM) on days 1 and 2 may have been caused by the caffeine overdose in the current case. Conclusions: The patient underwent infusion with bicarbonate Ringer's solution and potassium chloride and was discharged on the third day of hospitalization despite taking a potentially lethal dose of caffeine. The virtual plasma exposures of caffeine estimated using the current simplified PBPK model were higher than the measured values. The present results based on drug monitoring data and additional pharmacokinetic predictions could serve as a useful guide in cases of caffeine overdose. [ABSTRACT FROM AUTHOR]