Prediction of Drug Terminal Half-Life and Terminal Volume of Distribution After Intravenous Dosing Based on Drug Clearance, Steady-State Volume of Distribution, and Physiological Parameters of the Body

Interspecies Brain PBPK Modeling Platform to Predict Passive Transport through the Blood-Brain Barrier and Assess Target Site Disposition

The high failure rate of central nervous system (CNS) drugs is partly associated with an insufficient understanding of target site exposure. Blood-brain barrier (BBB) permeability evaluation tools are needed to explore drugs' ability to access the CNS. An outstanding aspect of physiologically based pharmacokinetic (PBPK) models is the integration of knowledge on drug-specific and system-specific characteristics, allowing the identification of the relevant factors involved in target site distribution. We aimed to qualify a PBPK platform model to be used as a tool to predict CNS concentrations when significant transporter activity is absent and human data are sparse or unavailable. Data from the literature on the plasma and CNS of rats and humans regarding acetaminophen, oxycodone, lacosamide, ibuprofen, and levetiracetam were collected. Human BBB permeability values were extrapolated from rats using inter-species differences in BBB surface area. The percentage of predicted AUC and Cmax within the 1.25-fold criterion was 85% and 100% for rats and humans, respectively, with an overall GMFE of <1.25 in all cases. This work demonstrated the successful application of the PBPK platform for predicting human CNS concentrations of drugs passively crossing the BBB. Future applications include the selection of promising CNS drug candidates and the evaluation of new posologies for existing drugs.

Theoretical Examination Seeking Tangible Physical Meanings of Slopes and Intercepts of Plasma Concentration-Time Relationships in Minimal Physiologically Based Pharmacokinetic Models

In minimal physiologically based pharmacokinetic (mPBPK) models, physiological (e.g., cardiac output) and anatomical (e.g., blood/tissue volumes) variables are utilized in the domain of differential equations (DEs) for mechanistic understanding of the plasma concentration-time relationships [Formula: see text]. Although fundamental biopharmaceutical variables in terms of distribution (e.g., [Formula: see text] and [Formula: see text]) and elimination kinetics (e.g., [Formula: see text]) in mPBPK provide greater insights in comparison to classical compartment models, an absence of kinetic elucidation of slopes and intercepts in light of such DE model parameters hinders more intuitive appreciation of [Formula: see text]. Therefore, this study seeks the tangible physical meanings of slopes and intercepts of the plasma concentration-time relationships in one- and two-tissue mPBPK models (i.e., m2CM and m3CM), with respect to time parameters that are readily understandable in PK analyses, i.e., the mean residence ([Formula: see text]) and transit ([Formula: see text]) times. Utilizing the explicit equations (EEs) for the slopes, intercepts, and areas of each exponential phase in the m2CM and m3CM, we theoretically and numerically examined the limiting/boundary conditions of such kinetic properties, based on the ratio of the longest tissue [Formula: see text] to the [Formula: see text] in the body (i.e., [Formula: see text]) that is useful for dissecting complex PBPK systems. The kinetic contribution of the area of each exponential phase to the total drug exposure was assessed to identify the elimination phase between the terminal and non-terminal phases of the [Formula: see text] in the m2CM and m3CM. This assessment provides improved understanding of the complexities inherent in all PBPK profiles and models.

Determinants of Biological Half-Lives and Terminal Slopes in Physiologically Based Pharmacokinetic Systems: Assessment of Limiting Conditions

In pharmacokinetic (PK) analyses, the biological half-life T_1/2 is usually determined in the terminal phase after drug administration, which is readily calculated from the relationship T_1/2 = ln2/λ_z where λ_z is the terminal-phase slope obtainable from non-compartmental analysis (NCA). Since kinetic understanding of λ_z has been limited to the theory of a one-compartment model, this study seeks kinetic determinants of λ_z in more complex plasma concentration-time profiles. We utilized physiologically based pharmacokinetic (PBPK) systems that are consistent with the assumptions of NCA (e.g., linear PK and elimination occurring from plasma) to interrelate λ_z and disposition kinetic parameters of PBPK models. In a mammillary form of PBPK models, the two boundary conditions of λ_z are the inverses of the mean residence time in the body (1/MRT_B = CL/V_SS) and the mean transit time through the kinetically largest tissue (1/MTT_max = Q_Tf_dR_b/V_TK_p). Importantly, the limiting conditions of λ_z between 1/MRT_B and 1/MTT_max are dependent on a simple product MRT_Bλ_z (P_det) and a simple ratio MTT_max/MRT_B (K_det), leading to introduction of the unitless product-ratio plot for determination of the limiting condition of λ_z in linear PK. We found that the MRT_Bλ_z value of 0.5 serves as a practical threshold determining whether λ_z is more closely associated with 1/MRT_B or 1/MTT_max. The current theory was applied for assessment of the terminal slope λ_z for observed PK data of various compounds in man and rat.

Determination of the Number of Tissue Groups of Kinetically Distinct Transit Time in Whole-Body Physiologically Based Pharmacokinetic (PBPK) Models II: Practical Application of Tissue Lumping Theories for Pharmacokinetics of Various Compounds

In our companion paper, we described the theoretical basis for tissue lumping in whole-body physiologically based pharmacokinetic (WB-PBPK) models and found that K_det, a coefficient for determining the number of tissue groups of distinct transit time in WB-PBPK models, was related to the fractional change in the terminal slope (FCT) when tissues were progressively lumped from the longest transit time to shorter ones. This study was conducted to identify the practical threshold of K_det by applying the lumping theory to plasma/blood concentration-time relationships of 113 model compounds collected from the literature. We found that drugs having K_det < 0.3 were associated with FCT < 0.1 even when all peripheral tissues were lumped, resulting in comparable plasma concentration-time profiles between one-tissue minimal PBPK (mPBPK) and WB-PBPK models. For drugs with K_det ≥ 1, WB-PBPK profiles appeared similar with two-tissue mPBPK models by applying the rule of FCT < 0.1 for lumping slowly equilibrating tissues. The two-tissue mPBPK model also appeared appropriate in terms of concentration-time profiles for drugs with 0.3 ≤ K_det < 1, although, some compounds (15.9% of the total cases), but not all, in this range showed a slight (maximum of 18.9% of the total AUC) deviation from WB-PBPK models, indicating that the two-tissue model, with caution, could still be used for those cases. Comparison of kinetic parameters between traditional (model-fitting) and current (theoretical calculation) mPBPK analyses revealed their significant correlations. Collectively, these observations suggest that the number of tissue groups could be determined based on the K_det/FCT criteria, and plasma concentration-time profiles from WB-PBPK could be calculated using equations significantly less complex.

Aqueous humour concentrations after topical apPlication of combinEd levofloxacin-dexamethasone eye dRops and of its single components: a randoMised, assEssor-blinded, parallel-group study in patients undergoing cataract surgery: the iPERME study

Purpose

To evaluate the penetration of levofloxacin and dexamethasone sodium phosphate into the aqueous humour (AH) after administration in combination and as single molecules. Evaluation of the penetration of those agents in the site of action and their pharmacodynamic potential activity in view of the intended clinical use after cataract surgery.

Methods

Randomised, assessor-blinded, parallel-group. Patients scheduled for cataract surgery were assigned in a 1:1:1 ratio to: levofloxacin + dexamethasone sodium phosphate (L-DSP), Levofloxacin (L) or Dexamethasone sodium phosphate (DSP) eye drops. Either test or reference drugs were instilled in the cul-de-sac twice, 90 and 60 min before paracentesis.

Results

A total of 125 patients completed the study. Fraction of dose absorbed in the anterior chamber was 3.8–4.2 · 10⁻⁴ for levofloxacin and 0.3–0.4 · 10⁻⁴ for dexamethasone, respectively. No notable differences in concentration of levofloxacin were found between L-DSP arm (1.970 nmol/ml) and L arm (2.151 nmol/ml). The concentrations of levofloxacin were well above the MICs for the most frequent Gram-positive and Gram-negative eye pathogens. Dexamethasone concentrations were slightly lower in L-DSP arm (0.030 nmol/ml) than in DSP arm (0.042 nmol/ml), but still in the pharmacodynamically active range in the site of action. The difference was not clinically relevant. DSP was not detected in any HA sample, suggesting its full hydrolysis to free dexamethasone.

Conclusion

Our results confirm that no interaction is evident on the corneal penetration of levofloxacin and dexamethasone which reach pharmacologically active concentrations when instilled as fixed combination eye drops to patients undergoing cataract surgery.

Trial registration

ClinicalTrials.gov Identifier: NCT03740659

Electronic supplementary material

The online version of this article (10.1007/s00228-020-02863-7) contains supplementary material, which is available to authorized users.

Systems Pharmacology for Investigation of the Mechanisms of Action of Traditional Chinese Medicine in Drug Discovery

As a traditional medical intervention in Asia and a complementary and alternative medicine in western countries, traditional Chinese medicine (TCM) has attracted global attention in the life science field. TCM provides extensive natural resources for medicinal compounds, and these resources are generally regarded as effective and safe for use in drug discovery. However, owing to the complexity of compounds and their related multiple targets of TCM, it remains difficult to dissect the mechanisms of action of herbal medicines at a holistic level. To solve the issue, in the review, we proposed a novel approach of systems pharmacology to identify the bioactive compounds, predict their related targets, and illustrate the molecular mechanisms of action of TCM. With a predominant focus on the mechanisms of actions of TCM, we also highlighted the application of the systems pharmacology approach for the prediction of drug combination and dynamic analysis, the synergistic effects of TCMs, formula dissection, and theory analysis. In summary, the systems pharmacology method contributes to understand the complex interactions among biological systems, drugs, and complex diseases from a network perspective. Consequently, systems pharmacology provides a novel approach to promote drug discovery in a precise manner and a systems level, thus facilitating the modernization of TCM.

Clearance kinetics of the VGF-derived neuropeptide TLQP-21

TLQP-21 is a multifunctional neuropeptide and a promising new medicinal target for cardiometabolic and neurological diseases. However, to date its clearance kinetics and plasma stability have not been studied. The presence of four arginine residues led us to hypothesize that its half-life is relatively short. Conversely, its biological activities led us to hypothesize that the peptide is still taken up by adipose tissues effectively. [¹²⁵I]TLQP-21 was i.v. administered in rats followed by chasing the plasma radioactivity and assessing tissue uptake. Plasma stability was measured using LC-MS. In vivo lipolysis was assessed by the palmitate rate of appearance.

RESULTS

A small single i.v. dose of [¹²⁵I]TLQP-21 had a terminal half-life of 110 min with a terminal clearance rate constant, k_t, of 0.0063/min, and an initial half-life of 0.97 min with an initial clearance rate constant, k_i, of 0.71/min. The total net uptake by adipose tissue accounts for 4.4% of the entire dose equivalent while the liver, pancreas and adrenal gland showed higher uptake. Uptake by the brain was negligible, suggesting that i.v.-injected peptide does not cross the blood-brain-barrier. TLQP-21 sustained isoproterenol-stimulated lipolysis in vivo. Finally, TLQP-21 was rapidly degraded producing several N-terminal and central sequence fragments after 10 and 60 min in plasma in vitro. This study investigated the clearance and stability of TLQP-21 peptide for the first time. While its pro-lipolytic effect supports and extends previous findings, its short half-life and sequential cleavage in the plasma suggest strategies for chemical modifications in order to enhance its stability and therapeutic efficacy.

Determining potential adverse effects in marine fish exposed to pharmaceuticals and personal care products with the fish plasma model and whole-body tissue concentrations

The Fish Plasma Model (FPM) was applied to water exposure and tissue concentrations in fish collected from two wastewater treatment plant impacted estuarine sites. In this study we compared predicted fish plasma concentrations to Cmax values for humans, which represents the maximum plasma concentration for the minimum therapeutic dose. The results of this study show that predictions of plasma concentrations for a variety of pharmaceutical and personal care products (PPCPs) from effluent concentrations resulted in 37 compounds (54%) exceeding the response ratio (RR = Fish [Plasma]/1%Cmax_total) of 1 compared to 3 compounds (14%) detected with values generated with estuarine receiving water concentrations. When plasma concentrations were modeled from observed whole-body tissue residues, 16 compounds out of 24 detected for Chinook (67%) and 7 of 14 (50%) for sculpin resulted in an RR_tissue value greater than 1, which highlights the importance of this dose metric over that using estuarine water. Because the tissue residue approach resulted in a high percentage of compounds with calculated response ratios exceeding a value of unity, we believe this is a more accurate representation for exposure in the field. Predicting plasma concentrations from tissue residues improves our ability to assess the potential for adverse effects in fish because exposure from all sources is captured. Tissue residues are also more likely to represent steady-state conditions compared to those from water exposure because of the inherent reduction in variability usually observed for field data and the time course for bioaccumulation. We also examined the RR in a toxic unit approach to highlight the importance of considering multiple compounds exhibiting a similar mechanism of action.

The Use of Gene Ontology Term and KEGG Pathway Enrichment for Analysis of Drug Half-Life

A drug's biological half-life is defined as the time required for the human body to metabolize or eliminate 50% of the initial drug dosage. Correctly measuring the half-life of a given drug is helpful for the safe and accurate usage of the drug. In this study, we investigated which gene ontology (GO) terms and biological pathways were highly related to the determination of drug half-life. The investigated drugs, with known half-lives, were analyzed based on their enrichment scores for associated GO terms and KEGG pathways. These scores indicate which GO terms or KEGG pathways the drug targets. The feature selection method, minimum redundancy maximum relevance, was used to analyze these GO terms and KEGG pathways and to identify important GO terms and pathways, such as sodium-independent organic anion transmembrane transporter activity (GO:0015347), monoamine transmembrane transporter activity (GO:0008504), negative regulation of synaptic transmission (GO:0050805), neuroactive ligand-receptor interaction (hsa04080), serotonergic synapse (hsa04726), and linoleic acid metabolism (hsa00591), among others. This analysis confirmed our results and may show evidence for a new method in studying drug half-lives and building effective computational methods for the prediction of drug half-lives.

Effect of Chinese herbs on CYP3A4 activity and expression in vitro

ETHNOPHARMACOLOGICAL RELEVANCE

Traditional Chinese Medicine (TCM) has become more popular among cancer patients in the Western world, who often use Chinese herbs as adjuvant therapy to reduce the adverse effects of conventional chemotherapy. However, pharmacokinetic (PK) interactions between Chinese herbs and anticancer drugs can occur and have dramatic consequences for these patients. Currently, only a few possible PK interactions between Chinese herbs and conventional Western drugs have been documented.

AIM OF THE STUDY

Since the drug-metabolizing enzyme cytochrome P450 3A4 (CYP3A4) contributes to most of the PK interactions with (anticancer) drugs, the effect of four Chinese herbs (Oldenlandia diffusa, Codonopsis tangshen, Rehmannia glutinosa and Astragalus propinquus) on the activity and expression of CYP3A4 was investigated in vitro.

MATERIALS AND METHODS

Ethanol and water-ethanol extracts of the four Chinese herbs were prepared from raw material. CYP3A4 inhibition was assessed by the use of Supersomes™ in a fluorescence assay. Furthermore, CYP3A4 induction was evaluated in a human pregnane X receptor (hPXR)-mediated CYP3A4 reporter gene assay and a quantitative real time PCR assay, both in human colon adenocarcinoma-derived LS180 cells (LS180).

RESULTS

Extracts of Oldenlandia diffusa, Codonopsis tangshen, Rehmannia glutinosa and Astragalus propinquus inhibited CYP3A4 in human CYP3A4 Supersomes™ (IC50 values: 17-83 µg/mL). Oldenlandia diffusa and Rehmannia glutinosa significantly induced PXR-mediated CYP3A4 (p<0.001). Oldenlandia diffusa also significantly induced CYP3A4 mRNA levels (p<0.001 at 250 µg/mL).

CONCLUSIONS

Concomitant use of Oldenlandia diffusa and Rehmannia glutinosa could result in induction of CYP3A4, leading to a reduced efficacy of drugs that are CYP3A4 substrates and have a narrow therapeutic window. Because of the possible enhanced toxicity caused by CYP3A4 inhibition, clinical effects of CYP3A4 inhibition by Astragalus propinquus and Codonopsis tangshen must also be taken into account. In conclusion, herb-drug interactions between Chinese herbs and various CYP3A4 substrates can occur. Further research to investigate the clinical relevance of the interactions caused by Oldenlandia diffusa, Codonopsis tangshen, Rehmannia glutinosa and Astragalus propinquus is required.

On the accuracy of a one-compartment approach for determination of drug terminal half-life

The drug terminal half-life (t(1/2)) is commonly predicted by a simplified one-compartment approach (t(1/2) = ln 2V(ss)/CL), where V(ss) and CL are the steady-state volume of distribution and the total body clearance of drug, respectively. The analysis of the accuracy of this approach is provided. It turns out that most often a simplified one-compartment calculation underestimates t(1/2) by no more than 25% for human, 26% for dog, 20% for monkey, 19% for rat, and 23% for mouse. Thus, the application of a one-compartment calculation of t(1/2) is well justifiable, except for the rare cases of very high drug clearance (CL/(rQ) ≳ 0.5), where r is the equilibrium blood-plasma concentration ratio, and Q is the cardiac output.