1Faculty of Pharmacy, Silpakorn University, Nakhon Pathom, Thailand, 73000, 2Department of Medicine, Yuwaprasart Waithayopathum Child Psychiatric Hospital, Samut Prakarn, Thailand 10270, 3Department of Pharmacy, Yuwaprasart Waithayopathum Child Psychiatric Hospital, Samut Prakarn, Thailand 10270
Email: tewthanom_k@silpakorn.edu
Received: 04 Sep 2018 Revised and Accepted: 23 Nov 2018
ABSTRACT
Objective: The purpose of the study was to estimate the pharmacokinetic parameters and calculate the precision (%RMSE) of the predicted lithium concentration equation.
Methods: This research was studied from blood Lithium levels of Children who had visited Yuwaprasart Waithayopathum Child Psychiatric Hospital from 31 December 2009-1 January 2011 The accuracy and precision of the equation were evaluated by the mathematical principle. The assist package software (WIN-NONLIN) was used to create pharmacokinetic parameters.
Results: Twenty-nine patients were recruited. The characteristics (mean+SD) presented as the following; age = 15.79+2.64 y, weight = 69.75+22.28 kg and the daily dose = 858.62+274.53 mg. The trough Lithium concentrations (mean+SD) = 0.56+19 mg/l. The kinetic parameter (Mean+SD) presented as the following; t1/2 =7.23+3.38 hr., ke =0.12+0.07 hr-1, Vd = 48.83+15.60 L, AUC0-12 =11.01+5.61 mg x hr/l. The population pharmacokinetic parameters; ke, and AUC (0-∞) = 0.151 hr-1 and 169.81, mg x hr./l respectively. The modified equation from Yukawa for lithium clearance calculation and prediction lithium concentrations were CL (mL/min) = [36.5+(0.242x BW (kg)-7.79]/Scr (mg/dL) and lithium concentration (mg/l) = 4 x Dose (mg)/CL (mL/min)x73.89, respectively with 6.39 %RMSE (test for 20 patients data).
Conclusion: In conclusion, the modified Yukawa equation may be used for the prediction of lithium concentration with 6.39% RMSE.
Keywords: Pharmacokinetics parameters, Lithium, Children, Accuracy
© 2019 The Authors. Published by Innovare Academic Sciences Pvt Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/)
DOI: http://dx.doi.org/10.22159/ijpps.2019v11i1.29517
Lithium is the drug of choice in treating patients with mood disorders, especially in children [1-3]. However, the difficulty for drug use is adverse drug reactions [4-7] that needed to monitor serum drug concentrations [8-13]. Yuwaprasart Waithayopathum Child Psychiatric Hospital which one of the hospital that responsible for caring child psychiatric disorders patient, for example; autistic child also use lithium. While lithium prescribing is most common in this group, it has a problem of initiative lithium dosing design and monitoring. The application of the individual pharmacokinetic model for dosing adjustment by estimation pharmacokinetic parameters in individualized patients will be useful for lithium therapy monitoring.
Therefore, the aim of the study is to estimate the pharmacokinetic parameters and calculate the precision (% RMSE) of the predicted lithium concentration equation in child psychiatric patients.
Population
The information of serum lithium concentration in children less than 18 y who used lithium treatment in Yuwaprasart Waithayopathum Child Psychiatric Hospital was collected since 1 January 2009 to 31 December 2011.
Inclusion criteria
The information of the following patients was used.
Patients who aged less than 18 y old.
Patients who used lithium continuously more than 1 w.
Patients who had information of renal and hepatic laboratory examination values.
Patients who had good compliance of drug therapy.
Exclusion criteria
Patients who had in-completed information such as laboratory values, drugs co-administration therapy.
Research design
The research design in this study was retrospective. The information was collected from children who received Lithium treatment from Yuwaprasart Waithayopathum Child Psychiatric Hospital during1 January 2009 to 31 December 2011. The accuracy and precision of the equation were evaluated was evaluated based on the mathematical principle, and WIN-NONLIN program (with authorized code PHXWNL092) was used for calculation of the pharmacokinetic parameters of the population. This research was approved by the ethical committee of the hospital. (Reg. no. 0 07).
Data analysis
The information was analyzed by Excel program version 2013. The data presented in percentage, mean and standard deviation. The percent of Root mean square errors (% RMSE) was used for evaluating the precision of the equation, and Student-paired t-test statistic was used for comparison of predicted and observed serum lithium concentrations with the level of significance equal to 0.05.
Patient characteristic
Twenty-nine patients information were collected as present in table 1. From the results, the mean age of participants was 16 y with the average weight of 70 kg. The mean lithium dose was 900 mg/day. The serum lithium concentration (mean+SD) was 0.56+0.19 mg/l. All patients had normal serum creatinine (mean+SD) was 0.60+0.19 mg/dl.
The number of patients who had serum lithium within therapeutic range (0.6-1.2 mg/l) was 13 (44.82%). While there were 16 patients (55.18%) had serum lithium outside therapeutic range.
According to use the win-nonlin program to estimate population pharmacokinetics parameters and individual pharmacokinetics parameters via mathematical formulations, the individual patients’ values were presented in table 2 with population value was shown below the table. In addition, lithium clearance and predicted serum lithium concentrations of patients were presented in table 3. With the modified Yugawa equation with adjusted by measurement unit for calculation of lithium clearance as seen below.
Then, it can calculate serum lithium concentration by the following equation.
*73.89 is the molecular weight of lithium carbonate
The lithium clearance and predicted serum lithium concentration from 20 patients who had completed data calculated from the above equation illustrated in table 4. Then the determination precision of estimation steady-state serum lithium concentration by comparison of observed value is available in table 5. The graph between observed serum lithium concentration and predicted serum lithium concentration demonstrated in fig. 1.
The distribution of serum lithium concentration by lithium dose was presented in fig. 2 and 3, respectively.
The distribution of calculated pharmacokinetic parameters compared with mean were explained by the graph in fig. 4-8.
Table 1: Patients’ characteristics
No | Age (y) | Weight (kg) | lithium dose (mg) | Serum lithium concentration (mg/l) | Serum creatinine (mg/dl) |
1 | 15 | 74.0 | 900 | 0.49 | NA |
2 | 15 | 51.0 | 1200 | 0.76 | 0.5 |
3 | 14 | 68.0 | 900 | 0.7 | 0.7 |
4 | 18 | 44.0 | 600 | 0.59 | 0.4 |
5 | 17 | 118.0 | 1200 | 0.49 | 0.7 |
6 | 17 | 65.0 | 900 | 0.32 | NA |
7 | 16 | 70.0 | 900 | 0.54 | 0.5 |
8 | 24 | 92.4 | 900 | 0.55 | 0.5 |
9 | 15 | 51.0 | 1200 | 0.76 | 0.5 |
10 | 14 | 49.0 | 900 | 0.6 | NA |
11 | 15 | 56.0 | 900 | 0.8 | 0.5 |
12 | 14 | 71.0 | 1200 | 0.73 | NA |
13 | 14 | 54.0 | 300 | 0.5 | NA |
14 | 13 | 57.5 | 300 | 0.58 | 0.6 |
15 | 14 | 65.0 | 600 | 0.1 | 0.49 |
16 | 17 | 93.4 | 600 | 0.64 | 0.8 |
17 | 15 | 58.0 | 300 | 0.16 | 0.6 |
18 | 8 | 38.0 | 600 | 0.39 | 0.3 |
19 | 15 | 50.0 | 900 | 0.8 | NA |
20 | 18 | 78.0 | 1200 | 0.78 | 0.8 |
21 | 19 | 125.0 | 600 | 0.29 | NA |
22 | 16 | 65.0 | 900 | 0.65 | 0.5 |
23 | 19 | 62.0 | 900 | 0.63 | 0.2 |
24 | 15 | 62.5 | 900 | 0.42 | NA |
25 | 16 | 118.0 | 1200 | 0.49 | 0.7 |
26 | 15 | 75.0 | 900 | 0.77 | NA |
27 | 17 | 88.0 | 900 | 0.75 | 0.8 |
28 | 17 | 77.0 | 1200 | 0.47 | 1 |
29 | 16 | 47.0 | 900 | 0.49 | 0.4 |
Mean | 15.7 | 69.8 | 859 | 0.56 | 0.6 |
SD | 2.6 | 22.28 | 274 | 0.19 | 0.19 |
Fig. 1: Graph of observed steady-state serum lithium concentration (Css observed, mg/l) and predicted serum lithium from the modified Yukawa equation. (C predicted)
Table 2: The estimation of pharmacokinetic parameters
No. | t1/2 (h) | ke (h) | Vd (L) | AUC (0-12) ( mg h/l) |
1 | NA | NA | 51.80 | NA |
2 | 5.02 | 0.14 | 35.70 | 14.62 |
3 | 8.52 | 0.08 | 47.60 | 13.95 |
4 | 3.62 | 0.19 | 30.80 | 6.10 |
5 | 11.66 | 0.06 | 82.60 | 14.67 |
6 | NA | NA | 45.50 | NA |
7 | 6.20 | 0.11 | 49.00 | 9.86 |
8 | 7.31 | 0.09 | 64.68 | 8.81 |
9 | 5.02 | 0.14 | 35.70 | 14.62 |
10 | NA | NA | 34.30 | NA |
11 | 5.36 | 0.13 | 39.20 | 10.65 |
12 | NA | NA | 49.70 | NA |
13 | NA | NA | 37.80 | NA |
14 | 6.54 | 0.11 | 40.25 | 4.22 |
15 | 5.79 | 0.12 | 45.50 | 6.62 |
16 | 11.77 | 0.06 | 65.38 | 9.35 |
17 | 6.58 | 0.11 | 40.60 | 4.21 |
18 | 2.43 | 0.29 | 26.60 | 4.75 |
19 | NA | NA | 35.00 | NA |
20 | 10.60 | 0.07 | 54.60 | 20.17 |
21 | NA | NA | 87.50 | NA |
22 | 5.91 | 0.12 | 45.50 | 10.13 |
23 | 2.29 | 0.30 | 43.40 | 4.12 |
24 | NA | NA | 43.75 | NA |
25 | 11.66 | 0.06 | 82.60 | 14.67 |
26 | NA | NA | 52.50 | NA |
27 | 11.38 | 0.06 | 61.60 | 14.40 |
28 | 13.15 | 0.05 | 53.90 | 25.35 |
29 | 3.79 | 0.18 | 32.90 | 8.98 |
Mean | 7.23 | 0.12 | 48.83 | 11.01 |
SD | 3.38 | 0.07 | 15.60 | 5.61 |
NA= not available, the population pharmacokinetic parameters, the population pharmacokinetic parameters as shown below, ke (h-1) = 0.151429 AUC (0-∞) (mghr/l) = 169.8113
Table 3: Lithium clearance and predicted serum lithium concentration
No. | lithium clearance (ml/min) | Predicted lithium concentration (mg/l) |
2 | 82.10 | 0.79 |
3 | 64.52 | 0.76 |
4 | 98.40 | 0.33 |
5 | 81.81 | 0.79 |
7 | 91.30 | 0.53 |
8 | 102.14 | 0.48 |
9 | 82.10 | 0.79 |
11 | 84.52 | 0.58 |
14 | 71.04 | 0.23 |
15 | 90.69 | 0.36 |
16 | 64.14 | 0.51 |
17 | 71.24 | 0.23 |
18 | 126.35 | 0.26 |
20 | 59.48 | 1.09 |
22 | 88.88 | 0.55 |
23 | 218.57 | 0.22 |
25 | 81.81 | 0.79 |
27 | 62.51 | 0.78 |
28 | 47.34 | 1.37 |
29 | 100.21 | 0.49 |
Mean | 88.46 | 0.60 |
SD | 35.51 | 0.30 |
Fig. 2: The distribution of serum lithium concentration (mg/l) versus dose
Table 4: Observed serum lithium concentration and predicted serum lithium concentration
No. | Observed serum lithium concentration (ml/min) | Predicted lithium concentration (mg/l) | Difference |
2 | 0.76 | 0.79 | 0.03 |
3 | 0.70 | 0.76 | 0.06 |
4 | 0.59 | 0.33 | -0.26 |
5 | 0.49 | 0.79 | 0.30 |
7 | 0.54 | 0.53 | -0.01 |
8 | 0.55 | 0.48 | -0.07 |
9 | 0.76 | 0.79 | 0.03 |
11 | 0.80 | 0.58 | -0.22 |
14 | 0.58 | 0.23 | -0.35 |
15 | 0.10 | 0.36 | 0.26 |
16 | 0.64 | 0.51 | -0.13 |
17 | 0.16 | 0.23 | 0.07 |
18 | 0.39 | 0.26 | -0.13 |
20 | 0.78 | 1.09 | 0.31 |
22 | 0.65 | 0.55 | -0.10 |
23 | 0.63 | 0.22 | -0.41 |
25 | 0.49 | 0.79 | 0.30 |
27 | 0.75 | 0.78 | 0.03 |
28 | 0.47 | 1.37 | 0.90 |
29 | 0.49 | 0.49 | 0.00 |
%RMSE | 0.0639 (6.39) | ||
t-test (p-value) | 0.72 |
Fig. 3: The graph of serum lithium concentrations (mg/l) versus dose (mg)
Fig. 4: The distribution of the volume distribution of lithium (L)
Fig. 5: The distribution of the elimination rate constant (ke) of lithium (h-1)
Fig. 6: The distribution half-life of lithium (h)
Fig. 7: The distribution of lithium clearance (ml/min)
Fig. 8: The distribution of area under the concentration-time curve at 0-12 h (AUC) of lithium (mg h/l)
The results from this study demonstrated that the modified Yukawa equation can use to predict serum lithium concentration with 93.61% accuracy (error 6.39%). It is not statistically significant differences (p=0.72) between observed and predicted serum lithium concentration. However, the variations among this group also found. The reason is because the various factors of patients, especially weight. Comparison with the previous report of the population clearance value found that this study had more lithium clearance than the previous report (10-40 ml/min) [13]. Therefore, using lithium dose beyond population clearance may result in treatment failure. Motoki et al. [14] found that Yukawa equation produced 90.44% accuracy (error 9.56%) in prediction lithium clearance. However, when the modified Yukawa equation in Thai subject study was tested in Suwan P. study. The accuracy of this equation was 92.05% (error 7.95%). It is quite similar to this study which found 93.6% accuracy of modified Yukawa equation (error 6.39%). Therefore, the modified Yukawa equation in this study may appropriate for lithium dose calculation in children population.
Other mathematical models were reviewed [15-32]. The model that study in Thai population only established in Suwan P. study [26] and the results found that the modified Yukawa could be used for lithium clearance estimation in an Adult Thai population with good precision. Therefore, Yukawa equation was chosen to modify in this study for children population. Cornelia et al. [33] applied the pharmacokinetics and pharmacodynamics knowledge for find appropriate lithium dosage adjustment in children. However, the research design was retrospective data collection, and there was a limitation in the application of pharmacodynamics parameter to the present study because of incomplete data. Since the pharmacodynamics parameter of Cornelia et al. study was Yong Mania Rating Scale and the hospital setting in this study did not use for routine work. Another discussion point is in Cornelia et al. study, it was focus only in children with bipolar disorders which is a difference from this study. Therefore, the results from Cornelia et al. may not apply in the wider group of children with other psychological disorder as this study. The further study is needed.
According to there has been limited use of Lithium in psychiatric patients, especially in children in a healthcare setting. Therefore, leading this study had some limitation in term of "sample size". In addition, in-completed data through retrospective design were available. Moreover, it takes time and lots of expenditure if the prospective randomized design will be conducted.
The modified Yukawa equation with use to calculate Lithium clearance as the following CL (ml/min) = [36.5+(0.242x BW (kg)-7.79]/Scr (mg/dl).
This modified clearance equation can use to calculate serum lithium concentration as the following equation with % %RMSE = 6.39.
The acknowledgment is sent to the Silpakorn University research institute for funding support.
The thankful appreciation brought to the healthcare officer at Yuwaprasart Waithayopathum Child Psychiatric Hospital, pharmacy students for kindly assistance about the convenient in searching medical records, collecting data.
Assistant Prof. Dr. Karunrat Tewthanom contributed to set the design activity, evaluation, interpretation and writing the manuscript for submission.
Dr. Rinsook Ongarjsakulman contributed to a screening of medical recorded, provide information on ordering serum lithium concentration monitoring.
Miss Theerarat Tankhum, RPh. contribute to coordinate with other healthcare professions to authorize using medical record, submit a proposal to the hospital ethical committee.
The author declares no conflict of interest
Scahill L, Farkas L, Hamrin V. Lithium in children and adolescents. J Child Adolesc Psychiatr Nurs 2001;14:89-93.
Lopez Larson M, Frazier JA. Empirical evidence for the use of lithium and anticonvulsants in children with psychiatric disorders. Harv Rev Psychiatry 2006;14:285-304.
Masi G, Milone A, Manfredi A, Pari C, Paziente A, Millepiedi S. Effectiveness of lithium in children and adolescents with conduct disorder: a retrospective naturalistic study. CNS Drugs 2009;23:59-69.
Micromedex, Lithium. Thomson Health care; 2012.
Bramness JG, Grøholt B, Engeland A, Furu K. The use of lithium, valproate or lamotrigine for psychiatric conditions in children and adolescents in Norway 2004-2007-a prescription database study. J Affect Disord 2009;117:208-11.
Harvard Medical School. Lithium-induced kidney problems. Serious problems are unusual, but monitoring is key. Harv Ment Health Lett 2009;26:6-7.
Grandjean EM, JM Aubry. Lithium: updated human knowledge using an evidence-based approach. Part II: clinical pharmacology and therapeutic monitoring. CNS Drugs 2009; 23:331-49.
Sharma S, S Joshi, RK Chadda. Therapeutic drug monitoring of lithium in patients with bipolar affective disorder: experiences from a tertiary care hospital in India. Am J Ther 2009;16:393-7.
Reischer H, CR Pfeffer. Lithium pharmacokinetics. J Am Acad Child Adolesc Psychiatry 1996;35:130-1.
Arancibia A. Pharmacokinetics of lithium in healthy volunteers after exposure to high altitude. Int J Clin Pharmacol Ther 2003;41:200-6.
Chen KP, WW Shen, ML Lu. Implication of serum concentration monitoring in patients with lithium intoxication. Psychiatry Clin Neurosci 2004;58:25-9.
Lewis AM. Monitoring patients on lithium in the primary care setting. Nurs Times 2004;100:35.
Ratanajamit C, Soorapan S, Doang-Ngern T, Waenwaisart W, Suwanchavalit L, Suwansiri S, et al. Appropriateness of therapeutic drug monitoring for lithium. J Med Assoc Thai 2006;89:1954-60.
Motoki T. Calculation of lithium clearance for clinical use based on renal excretion in Japanese patients. Int J Clin Pharmacol Pharmacother 2016;1:107.
Abou Auda HS, Al Yamani MJ, Abou Shaaban RR, Khoshhal SI. A new accurate method for predicting lithium clearance and daily dosage requirements in adult psychiatric patients. Bipolar Disord 2008;10:369-76.
ElDesoky ES, Kumar V, Alorainy MS, Hamdi MM, Derendorf H. Estimation of lithium clearance from routine clinical data in Egyptian bipolar patients. A population pharmacokinetic approach. Int J Clin Pharmacol Ther 2008;46:617-26.
Findling RL, Landersdorfer CB, Kafantaris V, Pavuluri M, McNamara NK, McClellan J, et al. First-dose pharmacokinetics of lithium carbonate in children and adolescents. J Clin Psychopharmacol 2010;30:404-10.
Huang HC, Chang YL, Lan TH, Chiu HJ, Liu WM, Lee TJ. Prediction of optimal litium dose for Taiwanese psychiatric patients. J Can Pet Technol 2008;33:115-21.
Higuchi S, Fukuoka R, Aoyama T, Horioka M. Predicting serum lithium concentration using bayesian method: a comparison with other methods. J Pharmacobiodyn 1988;11:158-74.
Jermain DM, ML Crismon, ES Martin. 3rd. Population pharmacokinetics of lithium. Clin Pharm 1991;10:376-81.
Kerbusch T, Mathot RA, Otten HM, Meesters EW, Van Kan HJ, Schellens JH, et al. Bayesian pharmacokinetics of lithium after an acute self-intoxication and subsequent hemodialysis: a case report. Pharmacol Toxicol 2002;90:243-5.
Pumpaisarnchai W, Kittirattanapaiboon P, Sanichwannakul K. A simple mathematical model for prediction lithium maintenance dose requirement for individual Thai patents. [Research Report], Changmai; Suan prung psychiatric hospital; 1996.
Schou M, K Thomsen, P Vestergaard. The renal lithium clearance and its correlations with other biological variables: observations in a large group of physically healthy persons. Clin Nephrol 1986;25:207-11.
Smith EG. Lithium dosing schedule and renal insufficiency. J Clin Psychiatry 2005;66:396-7.
Sproule BA, Bazoon M, Shulman KI, Turksen IB, Naranjo CA. Fuzzy logic pharmacokinetic modeling: application to lithium concentration prediction. Clin Pharmacol Ther 1997;62:29-40.
Suwan P. Prediction of steady-state serum lithium concentration by population model, [Thesis] in the pharmacy department. Faculty of Pharmacy, Silpakorn University: Meaung Nakhon Pathom, Thailand; 2009.
Taright N, Mentre F, Mallet A, Jouvent R. Nonparametric estimation of population characteristics of the kinetics of lithium from observational and experimental data: individualization of chronic dosing regimen using a new Bayesian approach. Ther Drug Monit 1994;16:258-69.
Tyrer SP, Grof P, Kalvar M, Shopsin B. Estimation of lithium dose requirement by lithium clearance, serum lithium and saliva lithium following a loading dose of lithium carbonate. Neuropsychobiology 1981;7:152-8.
Nichols TA, Drayton SJ, Borckardt J, Taber DJ. Lithium dosing equations: are they accurate? Ann Pharmacother 2014;48:596–600.
Williams PJ, JL Browne, RA Patel. Bayesian forecasting of serum lithium concentrations. Comparison with traditional methods. Clin Pharmacokinet 1989;17:45-52.
Wing YK1, Chan E, Chan K, Lee S, Shek CC. Lithium pharmacokinetics in Chinese manic-depressive patients. J Clin Psychopharmacol 1997;17:179-84.
Yukawa E, Nomiyama N, Higuchi S, Aoyama T. Lithium population pharmacokinetics from routine clinical data: the role of patient characteristics for estimating dosing regimens. Ther Drug Monit 1993;15:75-82.
Cornelia B. Lithium in pediatric patients with bipolar disorder: implications for selection of dosage regimens via population pharmacokinetics/pharmacodynamics. Clin Pharmacokinet 2017; 56:77-90.