Next Article in Journal
Mesenchymal Stromal Cells: Potential Option for COVID-19 Treatment
Next Article in Special Issue
Mechanistically Coupled PK (MCPK) Model to Describe Enzyme Induction and Occupancy Dependent DDI of Dabrafenib Metabolism
Previous Article in Journal
Permeation Challenges of Drugs for Treatment of Neurological Tuberculosis and HIV and the Application of Magneto-Electric Nanoparticle Drug Delivery Systems
Previous Article in Special Issue
Mathematical Modeling of Remdesivir to Treat COVID-19: Can Dosing Be Optimized?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Physiologically Based Pharmacokinetic Model for Predicting Diazepam Pharmacokinetics after Intravenous, Oral, Intranasal, and Rectal Applications

1
Department of Pharmacy Practice, Faculty of Pharmacy, Bahauddin Zakariya University, Multan 60800, Pakistan
2
Department of Pharmacology, Faculty of Pharmacy, Bahauddin Zakariya University, Multan 60800, Pakistan
3
Section of Pharmaceutics, University College of Pharmacy, Allama Iqbal Campus, University of the Punjab, Lahore 54000, Pakistan
4
Institute for Advanced Biosciences (IAB), Centre National de la Recherche Scientifique, Unite Mixte de Recherche 5309, Institut National de la Sante et de la Recherche Medicale U1209, Grenoble Alpes University, 38700 La Tronche, France
5
Department of Pharmaceutics, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
6
Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Pharmaceutics 2021, 13(9), 1480; https://doi.org/10.3390/pharmaceutics13091480
Submission received: 18 August 2021 / Revised: 7 September 2021 / Accepted: 13 September 2021 / Published: 15 September 2021
(This article belongs to the Special Issue Model-Informed Drug Discovery and Development)

Abstract

:
Diazepam is one of the most prescribed anxiolytic and anticonvulsant that is administered through intravenous (IV), oral, intramuscular, intranasal, and rectal routes. To facilitate the clinical use of diazepam, there is a need to develop formulations that are convenient to administer in ambulatory settings. The present study aimed to develop and evaluate a physiologically based pharmacokinetic (PBPK) model for diazepam that is capable of predicting its pharmacokinetics (PK) after IV, oral, intranasal, and rectal applications using a whole-body population-based PBPK simulator, Simcyp®. The model evaluation was carried out using visual predictive checks, observed/predicted ratios (Robs/pred), and the average fold error (AFE) of PK parameters. The Diazepam PBPK model successfully predicted diazepam PK in an adult population after doses were administered through IV, oral, intranasal, and rectal routes, as the Robs/pred of all PK parameters were within a two-fold error range. The developed model can be used for the development and optimization of novel diazepam dosage forms, and it can be extended to simulate drug response in situations where no clinical data are available (healthy and disease).

Graphical Abstract

1. Introduction

The amnesic and anxiolytic properties of benzodiazepines make this pharmacological class clinically useful in several indications, including insomnia, anxiety, muscle relaxation, and epilepsy, but their use is also associated with several unwanted side effects [1,2]. Among the benzodiazepines, diazepam is the most commonly prescribed medium potency long-acting drug [3,4]. Diazepam is uniquely metabolized in the liver into a pharmacologically active metabolite, N-desmethyldiazepam, by demethylation and hydroxylation [2,5,6]. Diazepam shows wide interindividual variability in its metabolism, which results in marked differences in its systemic concentrations. Moreover, the genetic variability in the expressions of CYP2C19 and CYP3A4 can affect its therapeutic efficacy and may lead to adverse drug reactions [7,8,9,10].
Diazepam is commonly employed as an anxiolytic and anticonvulsant in intravenous (IV), oral, intramuscular, intranasal, and rectal dosage forms [3,4]. It is administered through different routes of drug administration for the management of acute and chronic diseases [11]. Parenteral diazepam is preferred in hospitals because of its rapid action (especially to suppress seizures), but this route of administration is not convenient in ambulatory settings and is also associated with serious side effects, i.e., hypotension, dysrhythmias, etc. [12,13]. The other routes of administration, i.e., oral, intranasal, and rectal are considered more valuable when the IV route becomes inconvenient. Diazepam’s high permeability leads to rapid drug absorption via the rectal route, and diazepam is the only approved drug by U.S Food and Drug Administration (FDA) for rectal administration for out-patient treatment of early status epilepticus specifically in pediatrics [14,15,16]. On the other hand, the intranasal passage is an alternative to the oral and rectal route for diazepam administration for treating acute seizure management [17,18]. For any treatment, the therapeutic goal is only achieved if the drug is effective through the specified route of administration without causing any harm. To facilitate the clinical use of diazepam, there is a need to develop other non-oral formulations for humans when the oral route of administration is not convenient such as transdermal applications. Although such formulations have been examined and assessed in animals, they have not yet been tested in humans [19,20,21]. The development of non-oral dosage forms in humans can be facilitated by utilizing Physiologically Based Pharmacokinetics (PBPK), which presents countless opportunities for improvement in drug development in healthy and diseased populations [22].
PBPK models incorporate drug-dependent parameters along with the population-dependent system parameters in the presence of intrinsic or extrinsic factors to estimate the pharmacokinetics (PK) of the drug [23]. The concept of PBPK was established in 1937 [24], almost a decade ago, but for the last few years, PBPK modeling has been used as a tool for drug development and discovery, as described in several literature reviews [25,26,27,28]. PBPK models not only predict the clinical PK of a chemical entity but also reveal variables that may influence drug development, i.e., drug–drug interactions (DDI) or disease-states [29,30]. By allowing for the incorporation of in vitro drug release data, the PBPK modeling platforms [31,32] can be used to estimate the absorption, distribution, metabolism, and excretion (ADME) of the drug, which can assist in designing and optimizing novel dosage forms [33].
There are few published reports for the PBPK models for diazepam, which are focused on predicting its PK and DDIs in both humans and animals [34,35,36]. One study reported a PBPK model for diazepam to evaluate its DDIs with opioids (oxycodone, buprenorphine, and fentanyl) [35]. In contrast, another study aimed to develop a diazepam PBPK model by using different modeling and simulation tools after IV application only [34]. The third reported diazepam PBPK model was focused on the prediction of drug parameters in rats [36]. Since diazepam is administered through different routes of administration (IV, oral, intranasal, and rectal), if a PBPK model that can predict its PK after application of different dosage forms is developed, it may have many clinical applications. Keeping this in mind, a PBPK diazepam model was developed using a systematic model-building approach [37,38] that was capable of predicting its ADME through different routes of drug administration. This study aims to develop a PBPK model in a healthy population for the prediction and evaluation of diazepam pharmacokinetics after the administration of the drug through different routes, i.e., IV, oral, intranasal, and rectal.

2. Methods

2.1. Modeling Software

Simcyp® version 19 release 1 (Simcyp Ltd., Sheffield, UK), a population-based simulator, was employed for diazepam PBPK model development.

2.2. Modeling Approach

The process of PBPK model development commenced with the literature search for the extraction of drug-specific and population-specific parameters. The obtained data, i.e., molecular weight, fraction absorbed (fa), blood to plasma ratio (B:P), etc. were incorporated into the program to simulate and evaluate the IV profiles in the healthy adults. After the IV model was successfully developed in the adult population, oral PK profiles were simulated by incorporating oral drug absorption parameters, i.e., human jejunum permeability (Peff), etc. After successful oral data assessment, the clinical PK profiles of other routes of administration for diazepam were evaluated, i.e., intranasal and rectal. The layout for the developed diazepam PBPK model is given in Figure 1.

2.3. Model Parameters

The parameterization of the model began with a thorough review of the published literature, including in vitro and in vivo data for diazepam. The physicochemical properties of diazepam including molecular weight, octanol–water partition coefficient (LogPo:w), and acid dissociation constant (pKa) are given in Table 1. For estimating drug absorption via oral, intranasal, and rectal routes, the advanced dissolution, absorption, and metabolism (ADAM) model was utilized. The Peff was predicted using the number of hydrogen bond donors (HBD) and polar surface area (PSA). The predicted Peff value for diazepam was 12.434 × 10−4 cm/s by incorporating PSA and HBD values of 32.67 Ao and 0, respectively [39]. For the prediction of diazepam distribution in an adult population, a whole-body full PBPK model was applied for estimating the tissue-to-plasma partition co-efficient and volume of distribution at steady-state (Vss) using “Method-3” (The Rodger and Rowland method + ion membrane permeability) within the Simcyp®. The Michaelis-Menten constant (Km) and maximum rate of reaction (Vmax) for CYP2B6, CYP2C19, CYP3A4, and CYP3A5 were used for predicting drug clearance. The final diazepam-specific input parameters for the developed PBPK model are given in Table 1.

2.4. Pharmacokinetic Data

A detailed online literature search was executed to categorize the PK profiles (plasma concentration vs. time profiles or data) of diazepam in the adult population for model development. The studies were selected based on the PK profiles, the route of administration (IV, oral, intranasal, and rectal), the dose administered, age, the female proportion, and the duration of the study. The observed data were extracted by scanning PK profiles through GetData Graph Digitizer V.2.26.0 (http://getdata-graph-digitizer.com, accessed on 12 September 2021) [45].
In conclusion, 8 clinical studies and 20 concentrations vs. time profiles were added for diazepam model development; out of these 20 profiles, 6 were based on IV administration (71 individuals), 5 were profiles of oral (88 individuals), 6 were of intranasal administration (66 individuals), and 3 were of rectal diazepam administration (38 individuals) were used. Among all of these PK profiles, about one-third (two IV, two oral, two intranasal, and one rectal) were employed for diazepam PBPK model development and the remaining two-third (four IV, three oral, four intranasal, and two rectal) were utilized for verification of the model. For model evaluation, all of the observed clinical PK data were used. Table 2 shows the demographic data utilized for diazepam model construction and evaluation in the adult populations.

2.5. Model Evaluation

For every PK profile, 100 individuals were selected as a virtual population to perform a simulation by incorporating the parameters, i.e., the administration route, age range, female proportion, dosing, and study duration, as mentioned in the above-published studies. The evaluation of the developed diazepam PBPK model was initially performed with visual predictive checks by overlaying mean observed PK profiles on predicted data, which include the mean, 5th percentile and 95th percentile, and maximum and minimum predicted values. A non-compartmental analysis (NCA) was performed for the comparison of observed and predicted values of PK parameters including clearance (CL), area under the plasma concentration vs. time curve (AUC0–t), and maximum plasma concentration (Cmax) values. Moreover, the fold-error (observed/predicted ratios) and average fold error (AFE) of AUC0–t, Cmax, and CL were calculated using Equations (1) and (2). A two-fold error range (within 0.5−2-fold range) was used for the evaluating ratios (Robs/pred) for PK parameters [53].
Mean observed/predicted ratio
ratio obs pred = observed   value   of   PK   parameter predicted   value   of   PK   parameter
Average fold error
AFE = 10 log fold   error N
Additionally, the values of bioavailability parameters such as Fa (fraction of oral dose absorbed from the intestinal lumen), Fg (fraction of drug that escaped both intestinal first-pass metabolism and transporter-secretion available at hepatic portal blood), Fh (fraction of drug escaping hepatic first-pass elimination) were predicted after simulating oral and intranasal profiles. Through these parameters, the oral and intranasal data can be evaluated using predicted bioavailability, which can be calculated using the following equation:
F = F a · F g · F h

3. Results

3.1. Intravenous Dose Administration in the Adult Population

The observed and predicted systemic diazepam concentration profiles can be seen in Figure 2. The visual predictive checks showed that the developed model successfully predicted diazepam PK after its IV administration within the dose range of 2–10 mg. The mean AUC0–t and Cmax Robs/pred were 0.94 (95% CI 0.75−1.13) and 0.95 (95% CI 0.82−1.08), respectively. All of the PK parameters were within the two-fold error range (Table 3 and Table 4, Figure 3).

3.2. Oral Dose Administration in the Adult Population

The observed and predicted concentration profiles after the oral administration of a 2–10 mg dose of diazepam are given in Figure 4. The observed clinical data was in between the maximum and minimum values, as perceived by visual predictive checks. The mean Robs/pred of Cmax and CL were within the range of the two-fold error, as shown by the values of 0.93 (95% CI 0.75–1.11) and 1.14 (95% CI 0.90−1.39), respectively (Table 3 and Table 4, and Figure 3).

3.3. Intranasal Dose Administered in the Adult Population

The developed model effectively predicts diazepam concentration–time profiles after administration of the intranasal dose of 2–10 mg in adult individuals, as shown in Figure 5. The visual predictive checks represented that these predictions were in accordance with the observed data. Additionally, the PK parameters depicted the values within a two-fold range (0.5–2.0); the mean Robs/pred of AUC0–t was 0.93 (95% CI 0.76–1.10) and the Cmax value was 0.91 (95% CI 0.72–1.10), which can be seen in Table 3 and Table 4, and Figure 3.

3.4. Rectal Dose Administered in the Adult Population

The predicted diazepam concentration–time profiles after rectal administration doses of 10–15 mg were in complete agreement with the observed PK data (Figure 6). These results were further confirmed by the mean Robs/pred of the PK parameters of diazepam being within the two-fold error range (Table 3 and Table 4, and Figure 3) as AUC0–t after the rectal application was 1.18 (95% CI 0.80–1.56). Additionally, the AFE values represented that the model effectively predicted PK parameters after rectal administration.

4. Discussion

In the presented work, a PBPK model was developed to predict the ADME of diazepam after IV, oral, intranasal, and rectal administration. The systematic model building was commenced by predicting IV clinical profiles in adults by incorporating all drug and population-specific parameters in Simcyp®. After successful evaluation of IV profiles, absorption parameters were incorporated for the prediction and evaluation of oral PK clinical data. Intranasal and rectal PK profiles were also evaluated in the same way as IV and oral data.
The developed diazepam PBPK model represented comparable drug disposition in adults after IV administration as the mean Robs/pred value of AUC0–t with 95% CI was 0.94 (95% CI 0.75–1.13). Although diazepam is used as a first-line agent in the emergency management of seizures, it is difficult to administer IV diazepam to stop seizures in an ambulatory setup. The suitable diazepam alternatives in the outpatient setting are non-IV formulations.
The presented model successfully simulated oral diazepam PK data, and the observed and predicted systemic concentrations were in close association with the reported data as the mean Robs/pred Cmax value with 95% CI was 0.93 (95% CI 0.75–1.11). Moreover, the predicted bioavailability after the oral application was in the range of 76–97%, which is comparable with the reported value of 94% [54]. The oral route of administration is not feasible in suppressing seizures due to its late onset of action and congestion problems. Due to this reason, clinicians prefer other routes of administration, i.e., rectal and intranasal, for diazepam administration to suppress seizures using rapid diazepam action.
Apart from the IV and oral routes, the intranasal route of administration was used as an alternative for drug administration from the early 1980s [55]. The intranasal route is commonly employed because of its accessibility, its rapid action without first-pass metabolism, its patient conformity, and its noninvasiveness [56,57,58]. Diazepam is considered effective for intranasal administration due to its physicochemical properties and provides extended action compared with other benzodiazepines [59]. This intranasal route for diazepam is more attractive for the user compared with some other formulations, i.e., rectal [11]. For intranasal PK profiles of diazepam, the abovementioned strategy of IV and oral formulations were implemented with additional parameters, i.e., fa and ka. The presented model includes multiple intranasal formulations, i.e., solution, suspension, and supersaturated solutions. The supersaturated solution of diazepam is expected to be rapidly absorbed across the nasal mucosa soon after administration, which leads to the shorter residence time of the drug at the absorption site along with nasal drainage, which contributes to a decrease amount of drug available for absorption over time [60,61]. Therefore, in the case of a supersaturated solution, compared with other intranasal formulations, ka changed to accommodate the increased drug absorption because of rapid permeation in a short time. The model effectively predicted the overall PK parameters as Robs/pred with a 95% CI of AUC0–t was calculated as 0.93 (95% CI 0.76–1.10). Thus, a solution with a drug given through the intranasal route is absorbed more rapidly than suspension [62], but the developed model showed no significant increase in drug disposition as the AUC0–t of the intranasal solution (1532 ng.h/mL) was slightly higher than that of the intranasal suspension (1523 ng.h/mL). Furthermore, the reported bioavailability after intranasal administration was 97% [48,63], and the model predictions were in the range of 81–97%.
The other non-oral route for diazepam administration that has rapid absorption and showed up to 80% to 90% bioavailability is the rectal route. The therapeutic concentration of diazepam after rectal administration was achieved within 5 to 10 min, similar to that of IV administration. However, the use of a rectal formulation has social limitations in terms of patient conformity [16,64,65]. In the recent past, a diazepam gel for rectal administration was the only drug available with an immediate action for suppressing seizures in ambulatory settings [66]. The presented study, successfully predicted the PK profiles of different rectal formulations, i.e., suppository, solution, and gel, that can be recognized with comparable values of PK parameters. The mean Robs/pred with a 95% CI of CL was calculated as 0.85 (95% CI 0.56–1.13). The predicted AUC0-8 of the rectal gel (2028.11 ng.h/mL) [49] was lower than that of the AUC0-24 of the suppositories (2666.01 ng.h/mL) [15]. Although IV diazepam is commonly used to suppress prolonged seizures, seizures need to be stopped immediately at home before hospitalization. Apart from the intranasal route of diazepam, its rectal dosage forms are commonly employed to control seizures, as the only FDA-approved diazepam rectal gel (Diastat®, Bausch Health NJ, USA) is intended for use in pediatrics [67]. The developed PBPK model may be used to optimize the other non-oral formulations (transdermal) in humans by utilizing the in vitro data (dissolution data) of the drug using the already reported strategies [68].
Ideally, the treatment option should be chosen based on disease type, rapid onset of action, easy administration, and extended duration with minimal adverse effects [11]. In terms of route of administration, diazepam has various pros and cons, where IV administration has a swift onset of action, but at the same time, it is not convenient in the ambulatory setting, and on the other hand, the oral route has slow absorption and administration issues, i.e., swallowing in seizure [16]. Intranasal and rectal routes of administration also have unwanted effects, i.e., irritation at the site of administration, but they are commonly employed for treating seizures because of their rapid action and minimal side effects [11]. The comparative evaluation of the intranasal and the rectal formulations suggests that intranasal dosage forms are associated with a high variability in drug exposure and they often require the optimization of spray devices [60].
The previously published PBPK models for diazepam were focused on predicting its PK and DDIs in both humans and animals [34,35,36]. One study has used two Bayesian software for the development of a PBPK model of diazepam after the administration of an IV infusion only. Both tools produced very good fits at the individual and population levels [34]. The other study reported the development of full and minimal PBPK models for the evaluation of DDIs between opioids (fentanyl, oxycodone, and buprenorphine) and benzodiazepines (alprazolam, diazepam, midazolam, and triazolam). Full PBPK models were applied to diazepam and all opioids, while minimal PBPK models were developed for all benzodiazepines except diazepam [35]. The third reported PBPK model accounted for the parameter variability and uncertainty in the presence of qualitative and semi-quantitative data. The two approaches for the incorporation of parameter variability and uncertainty used in the study included the fuzzy and Monte Carlo simulations [36]. However, to date, there is no report for a PBPK model for diazepam that has been developed and evaluated in the adult population after implementing a systematic model-building approach. This is the first report of a whole-body full PBPK model for diazepam that has successfully predicted its PK after IV, oral, intranasal, and rectal administration.

5. Conclusions

Our diazepam PBPK model successfully predicted diazepam PK in the adult population after doses administered through IV, oral, intranasal, and rectal routes. The non-oral routes for diazepam administration are preferred in an ambulatory setting, especially as the IV route in seizing patients is not convenient in an outpatient setting and as it may lead to adverse effects. The developed model can be used for the development and optimization of novel diazepam dosage forms, and it can be extended to simulate drug response in situations where no clinical data are available as in diseased states and special populations (pediatrics).

Author Contributions

Conceptualization, S.K., M.F.R., I.I., H.S. and F.A.; formal analysis, S.K., M.F.R., A.M., A.u.R., W.A. and T.A.; funding acquisition, Y.A.B.J. and F.A.; Investigation, S.K., M.F.R., A.M., A.u.R., T.A., Y.A.B.J. and F.A.; methodology, S.K., M.F.R., I.I., H.S. and W.A.; project administration, S.K., M.F.R., I.I., A.M., Y.A.B.J. and F.A.; validation, S.K., M.F.R., H.S., A.u.R., W.A. and T.A.; writing—original draft, S.K., M.F.R., A.M., A.u.R., W.A. and F.A.; writing—review and editing, S.K., M.F.R., I.I., H.S., T.A. and Y.A.B.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Distinguished Scientist Fellowship program at King Saud University, Riyadh, Saudi Arabia, through research supporting project number (RSP-2021/131).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Distinguished Scientist Fellowship program at King Saud University, Riyadh, Saudi Arabia, for funding this work through research supporting project number (RSP-2021/131).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Crestani, F.; Löw, K.; Keist, R.; Mandelli, M.-J.; Möhler, H.; Rudolph, U. Molecular targets for the myorelaxant action of diazepam. Mol. Pharmacol. 2001, 59, 442–445. [Google Scholar] [CrossRef]
  2. Griffin, C.E.; Kaye, A.M.; Bueno, F.R.; Kaye, A.D. Benzodiazepine pharmacology and central nervous system–mediated effects. Ochsner J. 2013, 13, 214–223. [Google Scholar] [PubMed]
  3. Mandelli, M.; Tognoni, G.; Garattini, S. Clinical pharmacokinetics of diazepam. Clin. Pharmacokinet. 1978, 3, 72–91. [Google Scholar] [CrossRef] [PubMed]
  4. Manchikanti, L.; Christo, P.J.; Trescot, A.; Falco, F.J.E. Pain Medicine & Interventional Pain Management: A Comprehensive Review: Clinical Aspects; ASIPP Publishing: Paducah, KY, USA, 2011. [Google Scholar]
  5. Andersson, T.; Miners, J.; Veronese, M.; Birkett, D. Diazepam metabolism by human liver microsomes is mediated by both S-mephenytoin hydroxylase and CYP3A isoforms. Br. J. Clin. Pharmacol. 1994, 38, 131–137. [Google Scholar] [CrossRef] [Green Version]
  6. Kaplan, S.; Jack, M.; Alexander, K.; Weinfeld, R. Pharmacokinetic profile of diazepam in man following single intravenous and oral and chronic oral administrations. J. Pharm. Sci. 1973, 62, 1789–1796. [Google Scholar] [CrossRef]
  7. Fukasawa, T.; Suzuki, A.; Otani, K. Effects of genetic polymorphism of cytochrome P450 enzymes on the pharmacokinetics of benzodiazepines. J. Clin. Pharm. Ther. 2007, 32, 333–341. [Google Scholar] [CrossRef]
  8. Qin, X.P.; Xie, H.G.; Wang, W.; He, N.; Huang, S.L.; Xu, Z.H.; Ou-Yang, D.S.; Wang, Y.J.; Zhou, H.H. Effect of the gene dosage of CYP2C19 on diazepam metabolism in Chinese subjects. Clin. Pharmacol. Ther. 1999, 66, 642–646. [Google Scholar] [CrossRef] [PubMed]
  9. Bertilsson, L.; Henthorn, T.K.; Sanz, E.; Tybring, G.; Säwe, J.; Villén, T. Importance of genetic factors in the regulation of diazepam metabolism: Relationship to S-mephenytoin, but not debrisoquin, hydroxylation phenotype. Clin. Pharmacol. Ther. 1989, 45, 348–355. [Google Scholar] [CrossRef] [PubMed]
  10. Skryabin, V.Y.; Zastrozhin, M.S.; Torrado, M.V.; Grishina, E.A.; Ryzhikova, K.A.; Shipitsyn, V.V.; Galaktionova, T.E.; Sorokin, A.S.; Bryun, E.A.; Sychev, D.A. How do CYP2C19* 2 and CYP2C19* 17 genetic polymorphisms affect the efficacy and safety of diazepam in patients with alcohol withdrawal syndrome? Drug Metab. Pers. Ther. 2020, 35, 84–92. [Google Scholar] [CrossRef] [PubMed]
  11. Boddu, S.H.; Kumari, S. A Short Review on the Intranasal Delivery of Diazepam for Treating Acute Repetitive Seizures. Pharmaceutics 2020, 12, 1167. [Google Scholar] [CrossRef] [PubMed]
  12. Rogawski, M.A.; Heller, A.H. Diazepam buccal film for the treatment of acute seizures. Epilepsy Behav. 2019, 101, 106537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Greenblatt, D.; Arendt, R.; Abernethy, D.R.; Giles, H.; Sellers, E.; Shader, R. In vitro quantitation of benzodiazepine lipophilicity: Relation to in vivo distribution. BJA Br. J. Anaesth. 1983, 55, 985–989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Anderson, G.D.; Saneto, R.P. Current oral and non-oral routes of antiepileptic drug delivery. Adv. Drug Deliv. Rev. 2012, 64, 911–918. [Google Scholar] [CrossRef]
  15. Moolenaar, F.; Bakker, S.; Visser, J.; Huizinga, T. Biopharmaceutics of rectal administration of drugs in man IX. Comparative biopharmaceutics of diazepam after single rectal, oral, intramuscular and intravenous administration in man. Int. J. Pharm. 1980, 5, 127–137. [Google Scholar] [CrossRef]
  16. Maglalang, P.D.; Rautiola, D.; Siegel, R.A.; Fine, J.M.; Hanson, L.R.; Coles, L.D.; Cloyd, J.C. Rescue therapies for seizure emergencies: New modes of administration. Epilepsia 2018, 59, 207–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Wermeling, D.P.H.; Miller, J.L.; Archer, S.M.; Manaligod, J.M.; Rudy, A.C. Bioavailability and pharmacokinetics of lorazepam after intranasal, intravenous, and intramuscular administration. J. Clin. Pharmacol. 2001, 41, 1225–1231. [Google Scholar] [CrossRef] [Green Version]
  18. Vyas, T.K.; Shahiwala, A.; Marathe, S.; Misra, A. Intranasal drug delivery for brain targeting. Curr. Drug Deliv. 2005, 2, 165–175. [Google Scholar] [CrossRef]
  19. Schwarz, J.S.; Weisspapir, M.R.; Friedman, D.I. Enhanced transdermal delivery of diazepam by submicron emulsion (SME) creams. Pharm. Res. 1995, 12, 687–692. [Google Scholar] [CrossRef]
  20. Mehmood, Y. Preparation of diazipam delayed release patch, for anxiolytic treatment. Int. J. 2014, 3, 19–23. [Google Scholar]
  21. Xu, Q.; Zan, Y.; Gao, C. Preparation of diazepam transdermal gel and its bioavailability. Zhejiang Da Xue Xue Bao Yi Xue Ban = J. Zhejiang Univ. Med Sci. 2012, 41, 441–444. [Google Scholar]
  22. Cottura, N.; Howarth, A.; Rajoli, R.K.; Siccardi, M. The Current Landscape of Novel Formulations and the Role of Mathematical Modeling in Their Development. J. Clin. Pharmacol. 2020, 60, S77–S97. [Google Scholar] [CrossRef] [PubMed]
  23. Jones, H.; Chen, Y.; Gibson, C.; Heimbach, T.; Parrott, N.; Peters, S.; Snoeys, J.; Upreti, V.; Zheng, M.; Hall, S. Physiologically based pharmacokinetic modeling in drug discovery and development: A pharmaceutical industry perspective. Clin. Pharmacol. Ther. 2015, 97, 247–262. [Google Scholar] [CrossRef]
  24. Teorell, T. Kinetics of distribution of substances administered to the body, I: The extravascular modes of administration. Arch. Int. De Pharmacodyn. Et De Ther. 1937, 57, 205–225. [Google Scholar]
  25. Poulin, P.; Theil, F.-P. Prediction of pharmacokinetics prior to in vivo studies. II. Generic physiologically based pharmacokinetic models of drug disposition. J. Pharm. Sci. 2002, 91, 1358–1370. [Google Scholar] [CrossRef]
  26. Theil, F.-P.; Guentert, T.W.; Haddad, S.; Poulin, P. Utility of physiologically based pharmacokinetic models to drug development and rational drug discovery candidate selection. Toxicol. Lett. 2003, 138, 29–49. [Google Scholar] [CrossRef]
  27. Edginton, A.N.; Theil, F.-P.; Schmitt, W.; Willmann, S. Whole body physiologically-based pharmacokinetic models: Their use in clinical drug development. Expert Opin. Drug Metab. Toxicol. 2008, 4, 1143–1152. [Google Scholar] [CrossRef] [PubMed]
  28. Rowland, M.; Peck, C.; Tucker, G. Physiologically-based pharmacokinetics in drug development and regulatory science. Annu. Rev. Pharmacol. Toxicol. 2011, 51, 45–73. [Google Scholar] [CrossRef] [PubMed]
  29. De Buck, S.S.; Mackie, C.E. Physiologically based approaches towards the prediction of pharmacokinetics: In vitro–in vivo extrapolation. Expert Opin. Drug Metab. Toxicol. 2007, 3, 865–878. [Google Scholar] [CrossRef]
  30. Rasool, M.F.; Läer, S. Development and evaluation of a physiologically based pharmacokinetic model to predict carvedilol-paroxetine metabolic drug–drug interaction in healthy adults and its extrapolation to virtual chronic heart failure patients for dose optimization. Expert Opin. Drug Metab. Toxicol. 2021, 17, 717–724. [Google Scholar] [CrossRef] [PubMed]
  31. Jamei, M.; Marciniak, S.; Feng, K.; Barnett, A.; Tucker, G.; Rostami-Hodjegan, A. The Simcyp® population-based ADME simulator. Expert Opin. Drug Metab. Toxicol. 2009, 5, 211–223. [Google Scholar] [CrossRef] [PubMed]
  32. Zhuang, X.; Lu, C. PBPK modeling and simulation in drug research and development. Acta Pharm. Sin. B 2016, 6, 430–440. [Google Scholar] [CrossRef] [Green Version]
  33. Jamei, M.; Turner, D.; Yang, J.; Neuhoff, S.; Polak, S.; Rostami-Hodjegan, A.; Tucker, G. Population-based mechanistic prediction of oral drug absorption. AAPS J. 2009, 11, 225–237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Tsiros, P.; Bois, F.Y.; Dokoumetzidis, A.; Tsiliki, G.; Sarimveis, H. Population pharmacokinetic reanalysis of a Diazepam PBPK model: A comparison of Stan and GNU MCSim. J. Pharmacokinet. Pharmacodyn. 2019, 46, 173–192. [Google Scholar] [CrossRef] [PubMed]
  35. Ji, B.; Liu, S.; Xue, Y.; He, X.; Man, V.H.; Xie, X.-Q.; Wang, J. Prediction of drug–drug interactions between opioids and overdosed benzodiazepines using physiologically based pharmacokinetic (PBPK) modeling and simulation. Drugs RD 2019, 19, 297–305. [Google Scholar] [CrossRef] [Green Version]
  36. Gueorguieva, I.I.; Nestorov, I.A.; Rowland, M. Fuzzy simulation of pharmacokinetic models: Case study of whole body physiologically based model of diazepam. J. Pharmacokinet. Pharmacodyn. 2004, 31, 185–213. [Google Scholar] [CrossRef] [PubMed]
  37. Rasool, M.F.; Khalid, S.; Majeed, A.; Saeed, H.; Imran, I.; Mohany, M.; Al-Rejaie, S.S.; Alqahtani, F. Development and Evaluation of Physiologically Based Pharmacokinetic Drug–Disease Models for Predicting Rifampicin Exposure in Tuberculosis and Cirrhosis Populations. Pharmaceutics 2019, 11, 578. [Google Scholar] [CrossRef] [Green Version]
  38. Rasool, M.F.; Khalid, R.; Imran, I.; Majeed, A.; Saeed, H.; Alasmari, F.; Alanazi, M.M.; Alqahtani, F. Investigating the role of altered systemic albumin concentration on the disposition of theophylline in adult and pediatric patients with asthma by using the physiologically based pharmacokinetic approach. Drug Metab. Dispos. 2020, 48, 570–579. [Google Scholar] [CrossRef]
  39. Drugbank (Diazepam Compound Summary). Available online: https://go.drugbank.com/drugs/DB00829 (accessed on 14 December 2020).
  40. Pubchem (Diazepam Compound Summary). Available online: https://pubchem.ncbi.nlm.nih.gov/compound/Diazepam#section=Computed-Properties (accessed on 6 January 2021).
  41. Ye, M.; Nagar, S.; Korzekwa, K. A physiologically based pharmacokinetic model to predict the pharmacokinetics of highly protein-bound drugs and the impact of errors in plasma protein binding. Biopharm. Drug Dispos. 2016, 37, 123–141. [Google Scholar] [CrossRef]
  42. Gonzalez, F.; Coughtrie, M.; Tukey, R. Goodman and Gilman’s The Pharmacological Basis of Therapeutics; The McGraw-Hill Companies, Inc.: New York, NY, USA, 2006; pp. 71–92. [Google Scholar]
  43. FDA (Valium Label). Available online: https://www.accessdata.fda.gov/drugsatfda_docs/label/2016/013263s094lbl.pdf (accessed on 15 January 2021).
  44. Onof, S.; Hatanaka, T.; Miyazawa, S.; Tsutsui, M.; Aoyama, T.; Gonzalez, F.; Satoh, T. Human liver microsomal diazepam metabolism using cDNA-expressed cytochrome P450s: Role of CYP2B6, 2C19 and the 3A subfamily. Xenobiotica 1996, 26, 1155–1166. [Google Scholar] [CrossRef]
  45. GDG Digitizer. GetData Graph Digitizer (2.26); GDG Digitizer: 2013. Available online: http://getdata-graph-digitizer.com, (accessed on 12 September 2021).
  46. Gizurarson, S.; Gudbrandsson, F.K.; Jonsson, H.; BECHGAARD, E. Intranasal administration of diazepam aiming at the treatment of acute seizures: Clinical trials in healthy volunteers. Biol. Pharm. Bull. 1999, 22, 425–427. [Google Scholar] [CrossRef] [Green Version]
  47. Ivaturi, V.D.; Riss, J.R.; Kriel, R.L.; Siegel, R.A.; Cloyd, J.C. Bioavailability and tolerability of intranasal diazepam in healthy adult volunteers. Epilepsy Res. 2009, 84, 120–126. [Google Scholar] [CrossRef]
  48. Agarwal, S.K.; Kriel, R.L.; Brundage, R.C.; Ivaturi, V.D.; Cloyd, J.C. A pilot study assessing the bioavailability and pharmacokinetics of diazepam after intranasal and intravenous administration in healthy volunteers. Epilepsy Res. 2013, 105, 362–367. [Google Scholar] [CrossRef]
  49. Cloyd, J.C.; Lalonde, R.L.; Beniak, T.E.; Novack, G.D. A single-blind, crossover comparison of the pharmacokinetics and cognitive effects of a new diazepam rectal gel with intravenous diazepam. Epilepsia 1998, 39, 520–526. [Google Scholar] [CrossRef]
  50. Lau, S.; Slattery, J.T. Absorption of diazepam and lorazepam following intranasal administration. Int. J. Pharm. 1989, 54, 171–174. [Google Scholar] [CrossRef]
  51. Friedman, H.; Greenblatt, D.J.; Peters, G.R.; Metzler, C.M.; Charlton, M.D.; Harmatz, J.S.; Antal, E.J.; Sanborn, E.C.; Francom, S.F. Pharmacokinetics and pharmacodynamics of oral diazepam: Effect of dose, plasma concentration, and time. Clin. Pharmacol. Ther. 1992, 52, 139–150. [Google Scholar] [CrossRef] [PubMed]
  52. Greenblatt, D.J.; Harmatz, J.S.; Friedman, H.; Locniskar, A.; Shader, R.I. A large-sample study of diazepam pharmacokinetics. Ther. Drug Monit. 1989, 11, 652–657. [Google Scholar] [CrossRef]
  53. Rasool, M.F.; Ali, S.; Khalid, S.; Khalid, R.; Majeed, A.; Imran, I.; Saeed, H.; Usman, M.; Ali, M.; Alali, A.S. Development and evaluation of physiologically based pharmacokinetic drug-disease models for predicting captopril pharmacokinetics in chronic diseases. Sci. Rep. 2021, 11, 8589. [Google Scholar] [CrossRef] [PubMed]
  54. Divoll, M.; Greenblatt, D.J.; Ochs, H.R.; Shader, R.I. Absolute bioavailability of oral and intramuscular diazepam: Effects of age and sex. Anesth. Analg. 1983, 62, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ali, J.; Ali, M.; Baboota, S.; Kaur Sahni, J.; Ramassamy, C.; Dao, L. Potential of nanoparticulate drug delivery systems by intranasal administration. Curr. Pharm. Des. 2010, 16, 1644–1653. [Google Scholar] [CrossRef]
  56. Choi, H.-G.; Jung, J.-H.; Ryu, J.-M.; Yoon, S.-J.; Oh, Y.-K.; Kim, C.-K. Development of in situ-gelling and mucoadhesive acetaminophen liquid suppository. Int. J. Pharm. 1998, 165, 33–44. [Google Scholar] [CrossRef]
  57. Ugwoke, M.I.; Verbeke, N.; Kinget, R. The biopharmaceutical aspects of nasal mucoadhesive drug delivery. J. Pharm. Pharmacol. 2001, 53, 3–22. [Google Scholar] [CrossRef]
  58. Arora, P.; Sharma, S.; Garg, S. Permeability issues in nasal drug delivery. Drug Discov. Today 2002, 7, 967–975. [Google Scholar] [CrossRef]
  59. Schrier, L.; Zuiker, R.; Merkus, F.W.; Klaassen, E.S.; Guan, Z.; Tuk, B.; van Gerven, J.M.; van der Geest, R.; Groeneveld, G.J. Pharmacokinetics and pharmacodynamics of a new highly concentrated intranasal midazolam formulation for conscious sedation. Br. J. Clin. Pharmacol. 2017, 83, 721–731. [Google Scholar] [CrossRef] [Green Version]
  60. Ivaturi, V.D. Intranasal and Rectal Diazepam for Rescue Therapy: Assessment of Pharmacokinetics and Tolerability; University of Minnesota: Minneapolis, MN, USA, 2010. [Google Scholar]
  61. Hou, H.; Siegel, R.A. Enhanced permeation of diazepam through artificial membranes from supersaturated solutions. J. Pharm. Sci. 2006, 95, 896–905. [Google Scholar] [CrossRef]
  62. FDA. Nasal Spray and Inhalation Solution, Suspension, and Spray Drug Products—Chemistry, Manufacturing and Controls Documentation; FDA: New Hampshire, MD, USA, 2002. [Google Scholar]
  63. Hogan, R.E.; Tarquinio, D.; Sperling, M.R.; Klein, P.; Miller, I.; Segal, E.B.; Rabinowicz, A.L.; Carrazana, E. Pharmacokinetics and safety of VALTOCO (NRL-1; diazepam nasal spray) in patients with epilepsy during seizure (ictal/peri-ictal) and nonseizure (interictal) conditions: A phase 1, open-label study. Epilepsia 2020, 61, 935–943. [Google Scholar] [CrossRef]
  64. Detyniecki, K.; Van Ess, P.J.; Sequeira, D.J.; Wheless, J.W.; Meng, T.C.; Pullman, W.E. Safety and efficacy of midazolam nasal spray in the outpatient treatment of patients with seizure clusters—A randomized, double-blind, placebo-controlled trial. Epilepsia 2019, 60, 1797–1808. [Google Scholar] [CrossRef] [PubMed]
  65. Anderson, M. Buccal midazolam for pediatric convulsive seizures: Efficacy, safety, and patient acceptability. Patient Prefer. Adherence 2013, 7, 27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Sarma, A.K.; Khandker, N.; Kurczewski, L.; Brophy, G.M. Medical management of epileptic seizures: Challenges and solutions. Neuropsychiatr. Dis. Treat. 2016, 12, 467. [Google Scholar] [PubMed] [Green Version]
  67. Pellock, J.M. Safety of Diastat®, a rectal gel formulation of diazepam for acute seizure treatment. Drug Saf. 2004, 27, 383–392. [Google Scholar] [CrossRef]
  68. Chetty, M.; Rose, R.H.; Abduljalil, K.; Patel, N.; Lu, G.; Cain, T.; Jamei, M.; Rostami-Hodjegan, A. Applications of linking PBPK and PD models to predict the impact of genotypic variability, formulation differences, differences in target binding capacity and target site drug concentrations on drug responses and variability. Front. Pharmacol. 2014, 5, 258. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Systematic diagram for the development of the PBPK model of diazepam. Human jejunum permeability (Peff), the volume of distribution at steady-state (Vss), pharmacokinetic (PK), and intravenous (IV). The figure was produced using “Servier Medical Art” (https://www.smart.servier.com, accessed on 12 September 2021).
Figure 1. Systematic diagram for the development of the PBPK model of diazepam. Human jejunum permeability (Peff), the volume of distribution at steady-state (Vss), pharmacokinetic (PK), and intravenous (IV). The figure was produced using “Servier Medical Art” (https://www.smart.servier.com, accessed on 12 September 2021).
Pharmaceutics 13 01480 g001
Figure 2. Observed and predicted systemic concentration profiles of diazepam after IV application. The solid line (―) represents the mean value; the dotted lines (….) indicate the 5th and 95th percentiles; the dash line (---) shows the maximum and minimum predicted values; and the red filled circles () represent the mean observed data along with the standard deviation, where available. The figure contains systemic diazepam concentration profiles after administering IV doses of (A) 2 mg [46], (B) 5 mg [47], (C) 5 mg [48], (D) 7.5 mg [49], (E) 10 mg [15], and (F) 10 mg [50].
Figure 2. Observed and predicted systemic concentration profiles of diazepam after IV application. The solid line (―) represents the mean value; the dotted lines (….) indicate the 5th and 95th percentiles; the dash line (---) shows the maximum and minimum predicted values; and the red filled circles () represent the mean observed data along with the standard deviation, where available. The figure contains systemic diazepam concentration profiles after administering IV doses of (A) 2 mg [46], (B) 5 mg [47], (C) 5 mg [48], (D) 7.5 mg [49], (E) 10 mg [15], and (F) 10 mg [50].
Pharmaceutics 13 01480 g002
Figure 3. Mean ratio observed/predicted for the pharmacokinetic parameters of diazepam along with their 95% confidence intervals after IV, oral, intranasal, and rectal application. (A) Area under the curve (AUC0-1), (B) maximum systemic concentration (Cmax), (C) clearance.
Figure 3. Mean ratio observed/predicted for the pharmacokinetic parameters of diazepam along with their 95% confidence intervals after IV, oral, intranasal, and rectal application. (A) Area under the curve (AUC0-1), (B) maximum systemic concentration (Cmax), (C) clearance.
Pharmaceutics 13 01480 g003
Figure 4. Observed and predicted systemic concentration profiles of diazepam after oral application. The solid line (―) represents the mean value, the dotted lines (….) indicate the 5th and 95th percentiles, the dash line (---) shows the maximum and minimum predicted values, and the red filled circles () represent the mean observed data. The figure contains systemic diazepam concentration profiles after administering oral doses of (A) 2 mg [51], (B) 5 mg [51], (C) 10 mg [15], (D) 10 mg [52], and (E) 10 mg [51].
Figure 4. Observed and predicted systemic concentration profiles of diazepam after oral application. The solid line (―) represents the mean value, the dotted lines (….) indicate the 5th and 95th percentiles, the dash line (---) shows the maximum and minimum predicted values, and the red filled circles () represent the mean observed data. The figure contains systemic diazepam concentration profiles after administering oral doses of (A) 2 mg [51], (B) 5 mg [51], (C) 10 mg [15], (D) 10 mg [52], and (E) 10 mg [51].
Pharmaceutics 13 01480 g004
Figure 5. Observed and predicted systemic concentration profiles of diazepam after intranasal application. The solid line (―) represents the mean value; the dotted lines (….) indicate the 5th and 95th percentiles; the dash line (---) shows the maximum and minimum predicted values; and the red filled circles () represent the mean observed data along with the standard deviation, where available. The figure contains systemic diazepam concentration profiles after administering intranasal doses of (A) 2 mg [46], (B) 5 mg [47], (C) 10 mg [50], (D) 10 mg [47], (E) 10 mg [48], and (F) 10 mg [48].
Figure 5. Observed and predicted systemic concentration profiles of diazepam after intranasal application. The solid line (―) represents the mean value; the dotted lines (….) indicate the 5th and 95th percentiles; the dash line (---) shows the maximum and minimum predicted values; and the red filled circles () represent the mean observed data along with the standard deviation, where available. The figure contains systemic diazepam concentration profiles after administering intranasal doses of (A) 2 mg [46], (B) 5 mg [47], (C) 10 mg [50], (D) 10 mg [47], (E) 10 mg [48], and (F) 10 mg [48].
Pharmaceutics 13 01480 g005
Figure 6. Observed and predicted systemic concentration profiles of diazepam after rectal application. The solid line (―) represents the mean value; the dotted lines (….) indicate the 5th and 95th percentiles; the dash line (---) shows the maximum and minimum predicted values; and the red filled circles () represent the mean observed data along with the standard deviation, where available. The figure contains systemic diazepam concentration profiles after administering rectal doses of (A) 10 mg [15] (B) 10 mg [15], and (C) 15 mg [49].
Figure 6. Observed and predicted systemic concentration profiles of diazepam after rectal application. The solid line (―) represents the mean value; the dotted lines (….) indicate the 5th and 95th percentiles; the dash line (---) shows the maximum and minimum predicted values; and the red filled circles () represent the mean observed data along with the standard deviation, where available. The figure contains systemic diazepam concentration profiles after administering rectal doses of (A) 10 mg [15] (B) 10 mg [15], and (C) 15 mg [49].
Pharmaceutics 13 01480 g006
Table 1. Drug-specific input parameters of diazepam.
Table 1. Drug-specific input parameters of diazepam.
ParametersValuesRef.
Physicochemical Properties
Molecular Weight (g/mol) 284.74 [40]
LogPo:w2.82[35]
pKa3.4[41]
Absorption
Model
Peff,man (cm/s)12.434 × 10−4 a,b[35]
PSA (Ao)32.67[39]
HBD0[39]
Distribution
Model Prediction MethodThe Rodger and Rowland Method + Ion Membrane Permeability
B/P ratio0.58[41]
fu0.03[3,42]
Vss (L/kg)0.66 a,c[35,43]
For Intranasal Administration
Lung fa0.7 d
Lung ka1.6 d,e
Elimination
fumic0.59 a
N-demethylationVmax (pmol/min/pmol)Km (μM)
CYP2B63.6113[44]
CYP2C192.332
CYP3A414.81828
CYP3A51.8293
3-hydroxylation
CYP2B60.1150[44]
CYP2C1920.2846
CYP3A4151.32235
CYP3A548.4316
a Value predicted by Simcyp, b reported value, c the reported value range was 0.59−1 L/kg, d adjusted by manual optimization, e adjusted to 3.2 while simulating supersaturated solution, blood-to-plasma ratio (B:P), fraction unbound (fu), Michaelis–Menten constant (Km), maximum rate of metabolic formation (Vmax), the volume of distribution at steady-state (Vss), human jejunum permeability (Peff), polar surface area (PSA), hydrogen bonding donor (HBD), acid dissociation constant (pKa), octanol–water partition coefficient (LogPo:w), fraction absorbed (fa), and absorption rate constant (ka).
Table 2. Population data for IV, oral, intranasal, and rectal administration of diazepam in the adult population.
Table 2. Population data for IV, oral, intranasal, and rectal administration of diazepam in the adult population.
No.PopulationNo. of SubjectsDose
(mg)
Dosage FormAge (Years)Weight
(kg)
Female ProportionRef.
IV
1Healthy92IV solution20−3058−800[46]
2Healthy85IV solutionMean 28.3-0.33[47]
3Healthy245IV solution18−45Mean 71.80.2[48]
4Healthy207.5IV solutionMean 28.8Mean 72.60[49]
5Healthy91IV solution18−2558−700.33[15]
6Healthy110IV solution26−3760−850.2[50]
Oral
7Healthy112-19−35-0.27[51]
8Healthy115-19−35-0.27[51]
9Healthy4810Tablets18−4459.1−950[52]
10Healthy1110-19−35-0.27[51]
11Healthy910Tablet18−2558−700.33[15]
Intranasal
12Healthy92Solution20−3058–800[46]
13Healthy85Supersaturated solutionMean 28.3-0.33[47]
14Healthy2410Solution18−45Mean 71.80.2[48]
15Healthy2410Suspension18−45Mean 71.80.2[48]
16Healthy810Supersaturated solutionMean 28.3-0.33[47]
17Healthy110Solution26−3760−850.2[50]
Rectal
18Healthy910Solution18−2558−700.33[15]
19Healthy910Suppository18−2558−700.33[15]
19Healthy2015GelMean 28.8Mean 173.90[49]
Table 3. Observed and predicted values of PK parameters after IV, oral, intranasal, and rectal routes of administration.
Table 3. Observed and predicted values of PK parameters after IV, oral, intranasal, and rectal routes of administration.
Dose (mg)Cmax (ng/mL)RatioAUC0–t (ng.h/mL)RatioCL (L/h)RatioRef.
ObservedPredicted ObservedPredicted ObservedPredicted
IV
12157.06154.391.01285.45256.601.117.007.790.89[46]
25332.28452.920.73455.11522.200.8710.989.571.14[47]
35492.1479.771.02851.741305.300.655.873.831.53[48]
47.5536.44605.140.8812441428.080.876.275.461.14[49]
510623586.111.063505.373117.921.122.853.20.89[15]
610642.95634.131.014491.904364.791.022.222.290.97[50]
Oral
1262.0763.970.97330406.550.816.064.911.23[51]
25125.24159.950.78784.191012.400.776.374.931.29[51]
310325311.871.0434872965.751.172.863.370.80[15]
410255.53322.620.801540.282026.760.756.494.931.31[51]
510352318.991.101445.961493.590.966.916.661.03[52]
Intranasal
1246.2539.661.16161.73159.941.0112.3612.50.98[46]
25102.66127.340.80309.82352.620.8716.1314.181.13[47]
310218.15222.080.981513.781532.920.986.606.521.01[48]
410172.99222.480.7710841523.590.709.226.561.41[48]
510180.51254.690.70588.64710.800.8216.9814.061.20[47]
610236.82222.701.063558.443033.841.172.813.290.85[50]
Rectal
110297224.171.323566.352699.951.322.803.700.75[15]
210254229.401.103258.122666.011.223.063.750.81[15]
315384.86350.741.092070.022028.111.027.247.390.97[49]
Table 4. Robs/pred and AFE values of PK parameters in an adult population after IV, oral, intranasal, and rectal administration of diazepam.
Table 4. Robs/pred and AFE values of PK parameters in an adult population after IV, oral, intranasal, and rectal administration of diazepam.
ParametersRobs/predAFE
IV
AUC0–t0.940.93
CL1.101.08
Cmax0.960.95
Oral
AUC0–t0.900.89
CL1.151.13
Cmax0.940.93
Intranasal
AUC0–t0.930.92
CL1.101.09
Cmax0.920.90
Rectal
AUC0–t1.191.18
CL0.850.85
Cmax1.181.17
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khalid, S.; Rasool, M.F.; Imran, I.; Majeed, A.; Saeed, H.; Rehman, A.u.; Ashraf, W.; Ahmad, T.; Bin Jardan, Y.A.; Alqahtani, F. A Physiologically Based Pharmacokinetic Model for Predicting Diazepam Pharmacokinetics after Intravenous, Oral, Intranasal, and Rectal Applications. Pharmaceutics 2021, 13, 1480. https://doi.org/10.3390/pharmaceutics13091480

AMA Style

Khalid S, Rasool MF, Imran I, Majeed A, Saeed H, Rehman Au, Ashraf W, Ahmad T, Bin Jardan YA, Alqahtani F. A Physiologically Based Pharmacokinetic Model for Predicting Diazepam Pharmacokinetics after Intravenous, Oral, Intranasal, and Rectal Applications. Pharmaceutics. 2021; 13(9):1480. https://doi.org/10.3390/pharmaceutics13091480

Chicago/Turabian Style

Khalid, Sundus, Muhammad Fawad Rasool, Imran Imran, Abdul Majeed, Hamid Saeed, Anees ur Rehman, Waseem Ashraf, Tanveer Ahmad, Yousef A. Bin Jardan, and Faleh Alqahtani. 2021. "A Physiologically Based Pharmacokinetic Model for Predicting Diazepam Pharmacokinetics after Intravenous, Oral, Intranasal, and Rectal Applications" Pharmaceutics 13, no. 9: 1480. https://doi.org/10.3390/pharmaceutics13091480

APA Style

Khalid, S., Rasool, M. F., Imran, I., Majeed, A., Saeed, H., Rehman, A. u., Ashraf, W., Ahmad, T., Bin Jardan, Y. A., & Alqahtani, F. (2021). A Physiologically Based Pharmacokinetic Model for Predicting Diazepam Pharmacokinetics after Intravenous, Oral, Intranasal, and Rectal Applications. Pharmaceutics, 13(9), 1480. https://doi.org/10.3390/pharmaceutics13091480

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop