Skip to main content

Predicting quetiapine dose in patients with depression using machine learning techniques based on real-world evidence

Abstract

Background

Being one of the most widespread, pervasive, and troublesome illnesses in the world, depression causes dysfunction in various spheres of individual and social life. Regrettably, despite obtaining evidence-based antidepressant medication, up to 70% of people are going to continue to experience troublesome symptoms. Quetiapine, as one of the most commonly prescribed antipsychotic medication worldwide, has been reported as an effective augmentation strategy to antidepressants. The right quetiapine dose and personalized quetiapine treatment are frequently challenging for clinicians. This study aimed to identify important influencing variables for quetiapine dose by maximizing the use of data from real world, and develop a predictive model of quetiapine dose through machine learning techniques to support selections for treatment regimens.

Methods

The study comprised 308 depressed patients who were medicated with quetiapine and hospitalized in the First Hospital of Hebei Medical University, from November 1, 2019, to August 31, 2022. To identify the important variables influencing the dose of quetiapine, a univariate analysis was applied. The prediction abilities of nine machine learning models (XGBoost, LightGBM, RF, GBDT, SVM, LR, ANN, DT) were compared. Algorithm with the optimal model performance was chosen to develop the prediction model.

Results

Four predictors were selected from 38 variables by the univariate analysis (p < 0.05), including quetiapine TDM value, age, mean corpuscular hemoglobin concentration, and total bile acid. Ultimately, the XGBoost algorithm was used to create a prediction model for quetiapine dose that had the greatest predictive performance (accuracy = 0.69) out of nine models. In the testing cohort (62 cases), a total of 43 cases were correctly predicted of the quetiapine dose regimen. In dose subgroup analysis, AUROC for patients with daily dose of 100 mg, 200 mg, 300 mg and 400 mg were 0.99, 0.75, 0.93 and 0.86, respectively.

Conclusions

In this work, machine learning techniques are used for the first time to estimate the dose of quetiapine for patients with depression, which is valuable for the clinical drug recommendations.

Background

Depression is a severe affective mental disorder that is accompanied by a lack of pleasure, and the impairment of cognition, behavior and autonomic nerve function, which causes dysfunction in various spheres of individual and social life, severely limits psychosocial functioning, and diminishes quality of life [1, 2]. Being one of the most widespread, pervasive, and troublesome illnesses in the world [3,4,5], depression can affect individuals of any age. By 2020, depression is anticipated to overtake heart disease as the second-leading cause of disability or early death, according to estimates from the World Health Organization (WHO) [6]. As a common and disabling mental disorder [7], it is a serious global public health issue that not only results in personal misery for those affected but also places a large economic burden on both the patients and the entire society [8, 9]. When it comes to medicinal therapy for depressive disorders, the American Psychiatric Association recommends selective serotonin reuptake inhibitors (SSRI, such as sertraline) and serotonin–norepinephrine reuptake inhibitors (SNRI, such as duloxetine), as well as noradrenergic and specific serotonergic antidepressants (NaSSA, such as mirtazapine) [10, 11]. Regrettably, despite obtaining evidence-based antidepressant medication, up to 70% of people are going to continue to experience troublesome symptoms [12, 13].

According to the Canadian Network for Mood and Anxiety Treatments (CANMAT) guidelines and American Psychiatric Association Practice guidelines, atypical antipsychotics (AA), specifically the use of quetiapine has been reported as an effective augmentation strategy to antidepressants. Quetiapine is an atypical antipsychotic agent, which was first introduced in the pharmaceutical market in 1997 [14]. In 2010, the European Medicine Agency (EMA) approved the extended-release formulation of the drug, quetiapine XR, as an add-on to antidepressants when monotherapy gives suboptimal response [15]. Studies have shown that quetiapine (mean dose, 156.74 ± 97.6 mg/day) showed significant benefits for both response and remission rates compared to placebo [16, 17]. Despite its high effectiveness, its optimal use is limited by widely variant individual factors, including height, weight, age, medical history, and the CYP3A4 and CYP2D6 enzymes and so on [18]. Before achieving the quetiapine maintenance dose, these influencing factors make it challenging to reach the narrow therapeutic window, which is monitored by the therapeutic drug monitoring from AGNP and a sub- or supra-therapeutic recommended therapeutic reference range (200–750 ng/ml). This may render treatment ineffective or increase the risk of sedation, hypotension, dry mouth, constipation, and tachycardia. Therefore, it is critical to help clinicians select the appropriate quetiapine dose and individualized quetiapine treatment using prediction models.

Recently, there has been a trend toward using machine learning and deep learning methods to create customized medications based on research from real-world situations [19]. With the help of large-scale complex algorithms and datasets, machine learning and deep learning algorithms, a branch of artificial intelligence, are able to predict clinical outcomes with high accuracy [20, 21]. When predicting from a variety of variables, they can assess data-driven estimation and derive nonlinear variable linkages [20, 21]. Several studies have utilized machine learning and deep learning approaches to improve the model depiction of the complex link between individual characteristics and drug dose, such as a vancomycin treatment prediction system using Extreme Gradient Boosting (XGBoost) [22], and a brand-new warfarin maintaining dose prediction system using Light Gradient Boosting Machine (LightGBM) [23].

Herein, our goal was to build a prediction model of quetiapine adjusted dose in a stationary state using algorithms based on machine learning and deep learning to support clinical prescription decisions. We did this by maximizing the use of real-world data to find significant influencing variables for quetiapine dose.

Methods

Patients and data

We included 474 patients with depression, who were treated with quetiapine and hospitalized in the First Hospital of Hebei Medical University, from November 1, 2019, to August 31, 2022.

The inclusion criteria included the following: (1) patients who were diagnosed with depression and (2) patients who took quetiapine orally for a long time (at least for 3 days) at the same dose, and the blood concentration reached steady state at the time of blood collection. The exclusion criteria were as following: (1) patients older than 60 years were deleted; (2) patients with missing information (e.g., patient ID, medication record, etc.) were deleted; (3) samples that contained quetiapine at levels below the lower limit of quantification of 20 ng·ml−1 were eliminated; and (4) patients who were diagnosed with organic mental disorders or took psychoactive drug substance were deleted. The International Classification of diseases-10 (ICD-10) was used for diagnosis, and the supervising doctor made a diagnosis of depression. According to the Chinese Guidelines for the Diagnosis and Treatment of Mental Disorders 2020 Edition, antidepressants should be used as single as possible for patients with depression. When changing medicine is ineffective, combination therapy may be considered. The combination of two antidepressants with different mechanisms of action can be used, and other combinations include the combination of second-generation antipsychotics and lithium [24]. The First Hospital of Hebei Medical University is a tertiary hospital in local, and most of the patients who admitted in our hospital had poor effect after single antidepressant treatment and changing medicines in primary hospitals. Therefore, depending on the patient’s condition and guidelines, antidepressants combined with second generation antipsychotics (such as quetiapine) were commonly prescribed. Herein, the primary purpose of using quetiapine is the synergistic treatment of depression. All data were gathered from clinical paper records and computerized medical records held by the hospital for patients. Eventually, 308 eligible individuals were enrolled in this study. Figure 1 provides an illustration of the sample selection workflow.

Fig. 1
figure 1

Workflow of sample selection

Data collection and processing

Figure 2 provides an illustration of data collecting and processing. First, based on the database's available data, we collected 47 clinical variables, including quetiapine administration information (e.g., daily dose and concentration), demographic information (e.g., age, gender, weight, height), comorbidities (e.g., hypertension, diabetes, hyperlipidemia), combination medication (e.g., CYP3A4 enzyme inhibitors, and CYP2D6/CYP3A4 competitive substrates) and laboratory parameters (e.g., regular blood test, liver function, and renal function). Considering the missing rates or extremely unbalanced variables, we preprocessed the obtained data, and the variables' missing values were filled with the mean.

Fig. 2
figure 2

Process for establishing models and analyzing data

Variable selection and model establishment

As depicted in Fig. 2, univariate analysis was used to screen the variables after data collection from all relevant samples. Ultimately, variables which had p < 0.05 were selected. Based on the final dataset, the whole dataset was randomly divided into training cohort and testing cohort at the ratio of 80%: 20%. The data of the training cohort is used to train the model, and the test cohort is used to verify the final effect of the model. In this study, 246 subjects were in the training cohort and 62 subjects were in the testing cohort. Following that, key variables with p < 0.05 were chosen, and the daily dose of quetiapine was defined as the target variable. We created and evaluated nine different machine learning and deep learning models to compare the prediction abilities, including XGBoost, LightGBM, Random Forest (RF), Gradient Boosting (GBDT), Artificial Neural Network (ANN), Lasso Regression (LR), Support Vector Machine (SVM), TabNet and Decision Tree (DT). Assessment indicators were used for model evaluation, including precision, recall, F1-score, accuracy, sensitivity, and specificity. At the same time, we evaluated the effectiveness (AUROC) of quetiapine at various doses (100 mg/d, 200 mg/d, 300 mg/d, and 400 mg/d). Among these evaluation indicators, precision denotes the proportion of false positives [25]. Recall/sensitivity measures false negatives against true positives [25]. The F1-score is the harmonic average of the precision and recall [25]. Specificity measures false positives against true negatives [25]. The area under the ROC curve, or AUROC, is a comprehensive measurement that reflects the sensitivity and specificity of continuous variables. Accuracy is the proportion of correct predictions over the output results [25]. By contrasting the models' overall average accuracy, we may assess how well these models perform in terms of classification. After that, the confusion matrix, a special table used to view a classification model's performance, was then used to evaluate the prediction findings [26]. To avoid model overfitting and reduce bias, we used grid search combined with tenfold cross validation for hyperparameter tuning. Parameters of all nine models are displayed in Additional file 2: Table S1.

Statistical analysis

IBM SPSS version 26.0 was used for statistical research. (IBM Corporation, Armonk, New York, USA). In the comparison between training cohort and testing cohort, Mann–Whitney U test (non-normal distribution) and independent t test (normal distribution) were used to analyze the various continuous factors. Categorical data were analyzed by the Chi-squared test (n ≥ 5) or Fisher's exact test (n < 5). Statistical significance was set at p value < 0.05. Windows Python 3.9.12 was used to create each machine learning model.

Results

Baseline information

Table 1 displays the distribution of features across the complete dataset. This research included 308 depressed patients in total, 131 of whom were men and 177 of whom were women. Median (interquartile range, IQR) was used to characterize continuous variables, and frequency (percentage, %) was used to describe categorical variables. Patients’ average age was 19.00 (IQR 15.00–36.25) years. The median height and weight were 166.00 (IQR 160.00–172.00) cm and 65.00 (IQR 55.00–76.00) kg. Based on their daily quetiapine dose, the patients were separated into several groups, with 57 (18.51%) receiving a dose of 100 mg, 108 (35.06) receiving a dose of 200 mg, 74 (24.03%) receiving a dose of 300 mg, and 69 (22.03%) receiving a dose of 400 mg. The median value of the serum levels in the dataset was 213.06 (IQR 124.23–370.24) ng·ml−1. Comorbidities including hypertension, diabetes, and hyperlipidemia occupied 8.44%, 10.00%, and 7.14%, respectively. Combination medicine usage rates for CYP3A4 enzyme inhibitors were 0.32%, CYP3A4 competitive substrates were 5.19%, and CYP2D6 competitive substrates were 9.74%.

Table 1 Description of the study samples

Variable analysis

Considering extremely unbalanced variables, including hypertension, diabetes, hyperlipidemia, CYP3A4 enzyme inhibitors/inducers/competitive substrates, CYP2D6 enzyme inhibitors/competitive substrates, and variables with a missing rate greater than 50%, including weight, height may influence the predicted results of quetiapine, we preprocessed the obtained data before determining the significant associations between univariates. This led to a total of 38 candidate predictors, and finally four variables were selected which had p < 0.05, including quetiapine TDM value, age, mean corpuscular hemoglobin concentration (MCHC), and total bile acid (TBA), described in Table 2.

Table 2 Outcomes of the univariate analysis

Model establishment and validation

We developed and validated prediction models based on the selected features using nine algorithms (including XGBoost, LightGBM, RF, GBDT, SVM, LR, ANN, TabNet, and DT). Table 3 displays the performance of these models in testing cohort. The metrics of the XGBoost model outperformed those of other models and achieved the best overall performance, with precision = 0.91 ± 0.07, recall = 0.68 ± 0.1, F1 score = 0.78 ± 0.09, AUROC = 0.93 ± 0.04, sensitivity = 0.68 ± 0.1, and specificity = 0.98 ± 0.01 for predicting the daily dose of 100 mg quetiapine; precision = 0.67 ± 0.05, recall = 0.76 ± 0.09, F1 score = 0.71 ± 0.03, AUROC = 0.77 ± 0.04, sensitivity = 0.76 ± 0.09, and specificity = 0.78 ± 0.05 for predicting the daily dose of 200 mg quetiapine; precision = 0.67 ± 0.17, recall = 0.47 ± 0.09, F1 score = 0.54 ± 0.1, AUROC = 0.77 ± 0.08, sensitivity = 0.47 ± 0.09, and specificity = 0.93 ± 0.04 for predicting the daily dose of 300 mg quetiapine; precision = 0.64 ± 0.1, recall = 0.79 ± 0.08, F1 score = 0.7 ± 0.07, AUROC = 0.86 ± 0.06, sensitivity = 0.79 ± 0.08, and specificity = 0.86 ± 0.05 for predicting the daily dose of 400 mg quetiapine, and accuracy = 0.69 ± 0.03 for the entire XGBoost model. As a result, XGBoost was chosen to forecast the daily dose of quetiapine.

Table 3 Nine different algorithms' model performance metrics

On this basis, XGBoost calculated and ranked the importance scores of four selected variables, as shown in Table 4. Among them, the most important feature in the prediction model was discovered to be the quetiapine TDM value (importance score = 0.41 ± 0.02), followed by AGE (importance score = 0.23 ± 0.01), MCHC (importance score = 0.19 ± 0.01) and TBA (importance score = 0.18 ± 0.01).

Table 4 Importance score ranking of variables by XGBoost

Then, we evaluated the performance of XGBoost model with 4 variables (quetiapine TDM value, AGE, MCHC, and TBA) using a testing cohort of 62 patients. Figure 3 shows the AUROC values for XGBoost under different groups according to the daily dose of quetiapine. Typically, an AUROC has a value between 0.5 and 1.0, and the larger AUROC indicates the greater model classification effect. Based on different dose intervals, the patients were separated into four subgroups: those with a daily dose of 100 mg (11 cases), 200 mg (23 cases), 300 mg (14 cases), and 400 mg (14 cases). In different subgroups according to the quetiapine daily dose of 100 mg, 200 mg, 300 mg and 400 mg, AUROC were 0.99, 0.75, 0.93, and 0.86, respectively.

Fig. 3
figure 3

ROC curve at different doses. Class 0 indicates patients with daily dose of 100 mg, Class 1 indicates patients with daily dose of 200 mg, Class 2 indicates patients with daily dose of 300 mg, and Class 3 indicates patients with daily dose of 400 mg

Figure 4 summarizes the model’s performance in the testing cohort (62 cases) through confusion matrix. The model accurately predicted the dose regimen of 100 mg, 200 mg, 300 mg, and 400 mg quetiapine for 9, 15, 9, and 10 individuals, respectively. The evaluation indicators of four subgroups in the XGBoost model were calculated. The model can predict the dose regimen of 100 mg quetiapin with 100% precision and 82% recall rate; the dose regimen of 200 mg with 75% precision and 65% recall rate; the dose regimen of 300 mg with 69% precision and 64% recall rate; and the dose regimen of 400 mg with 50% precision and 71% recall rate, respectively. The results showed that the predicted quetiapine dose metrics agreed well with those from the clinically delivered plans for these patients.

Fig. 4
figure 4

Confusion matrix in the CatBoost model

Discussion

One of the most popular and efficient ways to treat depression in a therapeutic context is with antidepressants, which also has the ability to successfully slow the onset of disease in depressed patients. However, research indicates that only about half of major depressive disorder (MDD) patients receive antidepressants that work well for them, and only about a third of them experience remission [27]. The use of AAs as first-line medicines, notably quetiapine, is recommended by various current pharmacological augmentation guidelines for treating depression [17, 28, 29]. The right quetiapine dose and personalized quetiapine treatment are frequently challenging for clinicians.

To better estimate quetiapine dose during depression treatment and to find valid and accurate predictors, we compared the prediction abilities of quetiapine dose by applying nine machine learning and deep learning techniques for patients with depression. Ultimately, the XGBoost algorithm with the best performance (accuracy = 0.69) among nine models was selected to build the prediction model. Afterward, it can be observed that a number of 43 instances of quetiapine dose were properly predicted in the testing cohort. The overall accuracy of the model was 0.69. The moderate accuracy demonstrates that the effect of accurately predicting quetiapine dose is acceptable, and our findings may offer clinicians recommendations for prompt drug regimen adjustments. In addition, we performed dose subgroup analyses to show individual predictive performance across dose levels and to help refine model performance with continued recruitment of data for a given range of daily doses.

Calculations of the area under the concentration–time curve (AUC), for instance, provide the basis of classic pharmacokinetic studies. However, if the data are insufficient or cannot support a pharmacokinetic modeling technique, the model is erroneous [30]. Recently, it has been noticed that there is growing interest in novel statistical techniques, such as population pharmacokinetic (popPK) analysis. Nonlinear mixed-effects modeling (NONMEM) is the most popular method for this type of pharmacokinetic data analysis [31, 32]. Nevertheless, the PPK model is relatively inflexible to apply because of the explicit mathematical model used, and adding or removing a parameter may be challenging [33]. Machine learning, in contrast, is renowned for its self-organizational and learning skills, which let computers learn from “experience” without being explicitly taught [34, 35]. It is a form of artificial intelligence that enables systems to examine a wide range of data gathered from electronic health records (EHRs) and automatically learn from them using cutting-edge statistical and probabilistic techniques to make more precise predictions by building clever and efficient predictive models [36]. Recent years have seen a significant increase in study interest in the use of machine learning for clinical drug therapies, which leads to an increasingly significant impact on the development of personalized dosing, particularly in the choice of drug dose [37]. A few studies on the use of machine learning to forecast drug doses or blood concentrations have been reported [38,39,40,41,42,43].

In this study, we innovatively used machine learning and deep learning techniques to predict quetiapine dose based on real-world data. Machine learning models can be updated by automatically extracting EHR data and continuously monitoring physiological data, and are effective approaches to modeling real-world data. The commonly used PPK models have some limitations, such as difficulty in modeling, less consideration of influencing factors, and low accuracy. Herein, multi-level data mining was conducted by machine learning to screen out a variety of real-world influencing factors, to construct a more practical and accurate dose prediction model. Therefore, the combination of machine learning and dose prediction can help to improve the level of precision medicine in clinical.

We considered multiple algorithms for model establishment. DT is simple and easy to understand, but there is a risk of overfitting. RF uses bagging sampling, random attribute selection and model ensemble to address excessive risk decision tree learning. On the basis of RF, GBDT combined with Boosting establishes the connections between trees, making the forest an ordered collective decision-making system. XGBoost goes a step further than GBDT by adding regular terms to the objective function at each iteration to reduce the risk of overfitting, and it can integrate multiple decision trees to achieve the goal of regression or classification [44]. For models such as ANN and XGBoost, they perform quite well on large-scale datasets. However, good prediction results can also be obtained on small data sets by adjusting hyperparameters to avoid overfitting. Each algorithm has its advantages and disadvantages, the performance of different algorithms depends on the characteristics of the dataset, and the final selection of the algorithm is based on the computational results. Herein, we used grid search combined with tenfold cross validation to find the optimal hyperparameters and avoid overfitting to obtain the optimal model.

The significant predictor for predicting quetiapine dose is the quetiapine concentration. Several studies on psychotic disorders have identified that dose affects quetiapine concentration. According to a review, quetiapine had linear pharmacokinetics in the studied dose range, and had predictable pharmacokinetics [45]. Albantakis et al. have also quantified the relationship between daily dose and serum concentration in children and adolescents with psychotic and mood disorders. Between the daily dose and quetiapine serum levels (from trough samples) in the entire sample, they discovered a statistically significant, positive, but flimsy linear connection [46]. Among the crucial parameters we chose for our study's prediction model, the concentration was the most prominent influencing variable, and it was positively associated with quetiapine dose, which was in line with earlier research.

The effect of age on the metabolism of second-generation antipsychotics has been described in a few prior investigations. One study revealed that dose-adjusted concentrations of quetiapine increased by an average of 13% per decade from the age of 20 [47], while another found that the average concentrations were 67% higher in patients over the age of 70 compared to those between the ages of 18 and 69 [48]. Another study found that patients aged 65 and above had 50% higher plasma concentrations than younger patients [49]. For children and adolescents (10–17 years of age), at steady state, the pharmacokinetics of the parent compound were similar to adults. However, when adjusted for dose and weight, AUC and Cmax of the parent compound were 41% and 39% lower, respectively, in children and adolescents than in adults [50, 51]. In our study, patients older than 60 years were excluded because of the small number of senior patients that model can only learn little information. The ability of the elderly to metabolize and excrete drugs may be reduced, which may lead to the accumulation of drugs in the body, and liver and kidney function may also be affected. As a result, older people tend to require smaller doses of drugs. In this study, age is one of the most important feature in the final prediction model. In the following study, we will include more patients older than 60 years in the model to verify its generalizability.

In addition, some previous studies have indicated that low MCHC increases the likelihood of developing pathological disorders, such as poor functional status, dementia, and cognitive decline as well as morbidity and death [50,51,52,53]. Poor functional status, such as decreased ability to carry oxygen, may lead to changes in the pharmacokinetics of quetiapine and thus affect the dose of quetiapine. Meanwhile, because it is a measure determined from the haemoglobin concentration (HGB) divided by mean cellular volume (MCV) and red blood cell count (RBC), the MCHC is a good indicator to detect anaemia [54]. Depending on the demographic data investigated, anemia, a condition marked by a deficiency in hemoglobin in the blood, affects an estimated percentage from 2.9% to 60.1% of older persons [55]. Many illnesses, including malnutrition, obesity, cancer, chronic renal disease, are linked to anemic people, which may lead to changes in the pharmacokinetics of quetiapine and thus affect its dose.

Furthermore, hepatic metabolism accounts for the majority of quetiapine elimination, and less than 1% of the amount taken orally after a single administration was excreted unaltered, showing quetiapine is rapidly metabolized [56, 57]. According to studies, people with liver disease (n = 8) had a 30% lower mean oral clearance of quetiapine than patients with normal liver function. Two of 8 patients with hepatic impairment experienced a threefold increase in AUC and Cmax compared with healthy patients [56, 57]. TBA is closely related to liver function and abnormally high value suggests poor liver health. Abnormal TBA levels indicate that patients may have impaired liver function, which may inhibit metabolism of quetiapine in the liver, resulting in high quetiapine concentration and dose adjustment may be needed. In one word, quetiapine TDM value, age, MCHC, and TBA, show important associations with quetiapine dose, which could be used as the predictors in the individualized medication model of quetiapine, to help clinicians choose the reasonable regimen.

In different dose groups, according to Additional file 1: Figure S1, the blood concentration points of some patients with a dose of 200 mg are extreme outliers, and there is a crossover with the upper quartile concentration points of patients with a dose of 400 mg. Also, there is a crossover between the upper quartile concentration points of patients with a dose of 200 mg and the lower quartile concentration points of patients with a dose of 300 mg. All the situations of crossover may affect the clinician’s regimen choice and the prediction outcome in 200 mg group. Therefore, the AUROC for 200 mg group is lower than other dose groups.

Our model has some notable flaws. First, due to the availability of data, such as extremely uneven distribution, lots of missing values and so on, some variables were excluded. A future goal is to improve the model when a great deal of samples may be used to thoroughly study the factors. Second, our model has not been sufficiently tested on additional data sets. By using the model on a larger pooled data set, future studies could delve deeper into these problems. The identification of more potent predictors and the improvement of prediction accuracy are likely to result from the input of additional data. Third, due to the constraints of the test conditions, several pertinent patient characteristics (such as CYP450 polymorphisms) were excluded. Last, in this study, some underlying confounding factors were not analyzed, such as the using duration of quetiapine, prior use of antipsychotics, mood stabilizers and antidepressants before admission, drug combination of benzodiazepines, anxiolytics, and lithium, and complex clinical situations including severity of illness and multiple complications [58]. There is a drawback of real-world study that there exist some unknown confounders from real clinical settings. In future study, we expect to apply propensity score matching and stratified analysis for reducing confounding bias.

According to our knowledge, this research is the initial to use XGBoost algorithm for estimating the dose of quetiapine for patients with depression. Our study could identify important influencing variables for quetiapine dose by maximizing the use of real-world data to support quetiapine dose adjustments for each patient. In clinical applications, we expect to develop a web tool for drug dose calculation that can automatically generate recommended quetiapine doses by entering the values of key variables (such as quetiapine TDM value, age, MCHC, and TBA) based on electronic medical records, blood tests and TDM, providing clinical decision support to improve therapeutic response and reduce patient’s burden.

Conclusion

In this work, machine learning techniques are used for the first time to estimate the dose of quetiapine for patients with depression, which is important and valuable for the clinical drug recommendations. Our model was designed as a real-time assisting clinical decision support tool to balance the effect of quetiapine dose on both treatment efficacy and toxicity outcomes, and to maximize the benefit of treatment for each patient. Therefore, our study fills the gap in this research field.

Data availability

All data generated and analyzed during this study are included in this published article.

Abbreviations

XGBoost:

Extreme Gradient Boosting

LightGBM:

Light Gradient Boosting Machine

RF:

Random Forest

GBDT:

Gradient Boosting

SVM:

Support vector machine

LR:

Lasso Regression

ANN:

Artificial neural network

TabNet:

Attentive Interpretable Tabular Learning

DT:

Decision tree

TDM:

Therapeutic drug monitoring

AUROC:

Area under the receiver operating characteristic

WHO:

World Health Organization

SSRI:

Serotonin reuptake inhibitors

SNRI:

Serotonin–norepinephrine reuptake inhibitors

NaSSA:

Noradrenergic and specific serotonergic antidepressants

CANMAT:

Canadian network for mood and anxiety treatments

AA:

Atypical antipsychotics

EMA:

European Medicine Agency

XR:

Extended-release

CYP:

CytochromeP450

IQR:

Interquartile range

AFU:

A-L-fucosidase

HBDH:

Alpha-hydroxybutyrate dehydrogenase

GGT:

γ Glutamyl transpeptidase

ALT:

Alanine aminotransferase

NEUR:

Percentage of neutrophils

LDH:

Lactic dehydrogenase

MONOR:

Percentage of monocytes

BASOR:

Percentage of basophilic granulocyte

EOSR:

Percentage of eosinophils granulocyte

AST:

Aspartate aminotransferase

ALT:

Alanine aminotransferase

BUN:

Blood urea nitrogen

UA:

Uric acid

MCV:

Mean corpuscular volume

MPV:

Mean platelet volume

MCHC:

Mean corpuscular hemoglobin concentration

MCH:

Mean corpuscular hemoglobin

TBA:

Total bile acid

TBIL:

Total bilirubin

TP:

Total protein

LYP:

Percentage of lymphocytes

GLB:

Globulin

WBC:

White blood cells

AIB:

Albumin

DBIL:

Direct bilirubin

RDW:

Red cell distribution width

HCT:

Hematocrit

RBC:

Red blood cell count

Cr:

Creatinine

CK:

Creatine kinase

ChE:

Cholinesterase

ADA:

Adenosine deaminase

Hb:

Hemoglobin

IBIL:

Indirect Bilirubin

MCHC:

Mean corpuscular hemoglobin concentration

TBA:

Total bile acid (TBA)

NONMEM:

Nonlinear mixed-effects modeling (NONMEM)

AUC:

Concentration–time curve

popPK:

Population pharmacokinetic

EHRs:

Electronic health records

HGB:

Haemoglobin concentration

MCV:

Mean cellular volume

RBC:

Red blood cell count

References

  1. McCarron RM, Shapiro B, Rawles J, Luo J. Depression. Ann Intern Med. 2021;174:ITC65–80.

    Article  PubMed  Google Scholar 

  2. Malhi GS, Mann JJ. Depression. Lancet. 2018;392(10161):2299–312.

    Article  PubMed  Google Scholar 

  3. Greenberg PE, Fournier A, Sisitsky T, Pike CT, Kessler RC. The economic burden of adults with major depressive disorder in the United States (2005 and 2010). J Clin Psychiatry. 2015;76(2):155–62.

    Article  PubMed  Google Scholar 

  4. Holvast F, Massoudi B, Voshaar RC, Verhaak PFM. Non-pharmacological treatment for depressed older patients in primary care: a systematic review and meta-analysis. PLoS ONE. 2017;12(9):e0184666.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Davydow DS, Fenger-Grøn M, Ribe A, Pedersen H, Prior A, Vedsted P, et al. Depression and risk of hospitalisations and rehospitalisations for ambulatory care-sensitive conditions in Denmark: a population-based cohort study. BMJ Open. 2015;5:e009878.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Park SC, Oh HS, Oh DH, Jung SA, Na KS, Lee HY, et al. Evidence-Based, non-pharmacological treatment guideline for depression in Korea. J Korean Med Sci. 2014;29:12–22.

    Article  PubMed  CAS  Google Scholar 

  7. Hasin DS, Sarvet AL, Meyers JL, Saha TD, Ruan WJ, Stohl M, Grant BF. Epidemiology of adult DSM-5 major depressive disorder and its specifiers in the United States. JAMA Psychiat. 2018;75:336–46.

    Article  Google Scholar 

  8. Kleine-Budde K, Müller R, Kawohl W, Bramesfeld A, Moock J, Rössler W. The cost of depression—A cost analysis from a large database. J Affect Disord. 2013;147:137–43.

    Article  PubMed  Google Scholar 

  9. Okumura Y, Higuchi T. Cost of depression among adults in Japan. Prim Care Companion CNS Disord. 2011;13(3):26159.

    Google Scholar 

  10. Zhang L, Chen Y, Yue L, Liu Q, Montgomery W, Zhi L, et al. Medication use patterns, health care resource utilization, and economic burden for patients with major depressive disorder in Beijing, People’s Republic of China. Neuropsychiatr Dis Trea. 2016;20(12):941–9.

    Google Scholar 

  11. Han C, Wang SM, Lee SJ, Patkar AA, Masand PS, Pae CU. Second-generation antipsychotics in the treatment of major depressive disorder: current evidence. Expert Rev Neurother. 2013;13(7):851–70.

    Article  PubMed  CAS  Google Scholar 

  12. Valenstein M. Keeping our eyes on STAR*D. AJP. 2006;163:1484–6.

    Article  Google Scholar 

  13. Wiles N, Taylor A, Turner N, Barnes M, Campbell J, Lewis G, Morrison J, Peters TJ, Thomas L, Turner K, et al. Management of treatment-resistant depression in primary care: a mixed-methods study. Br J Gen Pr. 2018;68:e673–81.

    Article  Google Scholar 

  14. Rege S, Sura S, Aparasu RR. Atypical antipsychotic prescribing in elderly patients with depression. Res Social Adm Pharm. 2018;14:645–52.

    Article  PubMed  Google Scholar 

  15. EMA, Questions and answers on Seroquel XR and associated names (50, 150, 200, 300 and 400 mg prolonged-release tablets containing quetiapine), EMA Website. 2010. https://www.ema.europa.eu/en/documents/referral/questions-answers-seroquel-xr-associated-names-50-150-200-300-400-mg-prolonged-release-tablets_en.pdf. Accessed 27 Sep 2019.

  16. Cleare A, Pariante CM, Young AH, et al. Evidence-based guidelines for treating depressive disorders with antidepressants: a revision of the 2008 British Association for Psychopharmacology guidelines. J Psychopharmacol. 2015;29:459–525.

    Article  PubMed  CAS  Google Scholar 

  17. Kennedy SH, Lam RW, McIntyre RS, et al. Canadian network for mood and anxiety treatments (CANMAT) 2016 clinical guidelines for the management of adults with major depressive disorder: section 3 Pharmacological treatments. Can J Psychiatry. 2016;61:540–60.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Hiemke C, Bergemann N, Clement HW, Conca A, Deckert J, Domschke K, Eckermann G, Egberts K, Gerlach M, Greiner C, Gründer G, Haen E, Havemann-Reinecke U, Hefner G, Helmer R, Janssen G, Jaquenoud E, Laux G, Messer T, Mössner R, Müller MJ, Paulzen M, Pfuhlmann B, Riederer P, Saria A, Schoppek B, Schoretsanitis G, Schwarz M, Gracia MS, Stegmann B, Steimer W, Stingl JC, Uhr M, Ulrich S, Unterecker S, Waschgler R, Zernig G, Zurek G, Baumann P. Consensus guidelines for therapeutic drug monitoring in neuropsychopharmacology: update 2017. Pharmacopsychiatry. 2018;51(1–02):e1. https://0-doi-org.brum.beds.ac.uk/10.1055/s-0037-1600991. (Epub 2018 Feb 1. Erratum for: Pharmacopsychiatry. 2018 Jan;51(1–02):9–62).

    Article  PubMed  CAS  Google Scholar 

  19. Avanzo M, Wei L, Stancanello J, Vallieres M, Rao A, Morin O, Mattonen SA, El Naga I. Machine and deep learning methods for radiomics. Med Phys. 2020;47(5):e185–202.

    Article  PubMed  Google Scholar 

  20. Gautier T, Ziegler LB, Gerber MS, Campos-Náñez E, Patek SD. Artificial intelligence and diabetes technology: a review. Metabolism. 2021;124:154872.

    Article  PubMed  CAS  Google Scholar 

  21. Rani P, Kotwal S, Manhas J, Sharma V, Sharma S. Machine learning and deep learning based computational approaches in automatic microorganisms image recognition: methodologies, challenges, and developments. Arch Comput Methods. 2021;29:1–37.

    Google Scholar 

  22. Huang X, Yu Z, Wei X, Shi J, Wang Y, Wang Z, Chen J, Bu S, Li L, Gao F, Zhang J, Xu A. Prediction of vancomycin dose on high-dimensional data using machine learning techniques. Expert Rev Clin Pharmacol. 2021;14(6):761–71.

    Article  PubMed  CAS  Google Scholar 

  23. Liu Y, Chen J, You Y, Xu A, Li P, Wang Y, Sun J, Yu Z, Gao F, Zhang J. An ensemble learning based framework to estimate warfarin maintenance dose with cross-over variables exploration on incomplete data set. Comput Biol Med. 2021;131:104242.

    Article  PubMed  CAS  Google Scholar 

  24. National Health Commission of the People’s Republic of China. Code for Diagnosis and Treatment of Mental Disorders 2020 [M]. The National Health Commission of the People’s Republic of China. 2020;5:164.

  25. Chen H, Ma Y, Hong N, Wang H, Su L, Liu C, He J, Jiang H, Long Y, Zhu W. Early warning of citric acid overdose and timely adjustment of regional citrate anticoagulation based on machine learning methods. BMC Med Inform Decis Mak. 2021;21(Suppl 2):126.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Powers DMW. Evaluation: from precision, recall and F-measure to ROC, informedness, markedness and correlation. arXiv. Preprint posted online. Accessed 11 Oct 2020.

  27. Trivedi MH, Rush AJ, Wisniewski SR, Nierenberg AA, Warden D, Ritz L, Fava M. Evaluation of outcomes with citalopram for depression using measurement-based care in STAR*D: Implications for clinical practice. Am J Psychiatry. 2006;163(1):28–40.

    Article  PubMed  Google Scholar 

  28. National Collaborating Centre for Mental Health (UK). Depression: The Treatment and Management of Depression in Adults (Updated Edition). Leicester (UK): British Psychological Society. 2010

  29. Seshadri A, Wermers ML, Habermann TJ, et al. Long-term efficacy and tolerability of adjunctive aripiprazole for major depressive disorder: systematic review and meta-analysis. Prim Care Companion CNS Disord. 2021;23(4):34898.

    Article  Google Scholar 

  30. You W, Widmer N, De Micheli G. Example-based support vector machine for drug concentration analysis. In: Conf Proc IEEE Eng Med Biol Soc. 2011, 153–157.

  31. Ludden TM. Population pharmacokinetics. J Clin Pharmacol. 1988;28:1059–63.

    Article  PubMed  CAS  Google Scholar 

  32. Johansson ÅM, Ueckert S, Plan EL, Hooker AC, Karlsson MO. Evaluation of bias, precision, robustness and runtime for estimation methods in NONMEM 7. J Pharmacokinet Pharmacodyn. 2014;41:223–38.

    Article  PubMed  CAS  Google Scholar 

  33. Sibieude E, Khandelwal A, Girard P, Hesthaven JS, Terranova N. Population pharmacokinetic model selection assisted by machine learning. J Pharmacokinet Pharmacodyn. 2021. https://0-doi-org.brum.beds.ac.uk/10.1007/s10928-021-09793-6.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Huang X, Yu Z, Bu S, Lin Z, Hao X, He W, et al. An ensemble model for prediction of vancomycin trough concentrations in pediatric patients. Drug Des Devel Ther. 2021;15:1549–59.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Poynton MR, Choi BM, Kim YM, Park IS, Noh GJ, Hong SO, et al. Machine learning methods applied to pharmacokinetic modelling of remifentanil in healthy volunteers: a multi-method comparison. J Int Med Res. 2009;37:1680–91.

    Article  PubMed  CAS  Google Scholar 

  36. Shatte A, Hutchinson DM, Teague SJ. Machine learning in mental health: a scoping review of methods and applications. Psychol Med. 2019;49:1426–48.

    Article  PubMed  Google Scholar 

  37. Meng HY, Jin WL, Yan CK, Yang H. The application of machine learning techniques in clinical drug therapy. Curr Comput Aided Drug Des. 2019;15:111–9.

    Article  PubMed  CAS  Google Scholar 

  38. Jovanović M, et al. Application of counter-propagation artificial neural networks in prediction of topiramate concentration in patients with epilepsy. J Pharm Pharm Sci. 2015;18:856–62. https://0-doi-org.brum.beds.ac.uk/10.18433/j33031.

    Article  PubMed  Google Scholar 

  39. Tang J, et al. Application of machine-learning models to predict tacrolimus stable dose in renal transplant recipients. Sci Rep. 2017;7:42192.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Liu R, Li X, Zhang W, Zhou HH. Comparison of nine statistical model based warfarin pharmacogenetic dosing algorithms using the racially diverse international warfarin pharmacogenetic consortium cohort database. PLoS ONE. 2015;10:e0135784.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ma Z, Wang P, Gao Z, Wang R, Khalighi K. Ensemble of machine learning algorithms using the stacked generalization approach to estimate the warfarin dose. PLoS ONE. 2018;13:e0205872.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Roche-Lima A, et al. Machine learning algorithm for predicting warfarin dose in Caribbean hispanics using pharmacogenetic data. Front Pharmacol. 2020;10:1550.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Chen SS, et al. Optimizing levothyroxine dose adjustment after thyroidectomy with a decision tree. J Surg Res. 2019;244:102–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Chen T, Guestrin C. In: Proceedings of the 22nd ACM SIGKDD international conference on knowledge discovery and data mining. XGBoost: a scalable tree boosting system. San Francisco, CA: ACM, 2016. p. 785–94.

  45. Bui K, Earley W, Nyberg S. Pharmacokinetic profile of the extended-release formulation of quetiapine fumarate (quetiapine XR): clinical implications. Curr Med Res Opin. 2013;29(7):813–25.

    Article  PubMed  CAS  Google Scholar 

  46. Albantakis L, Egberts K, Burger R, Kulpok C, Mehler-Wex C, Taurines R, Unterecker S, Wewetzer C, Romanos M, Gerlach M. Relationship between daily dose, serum concentration, and clinical response to quetiapine in children and adolescents with psychotic and mood disorders. Pharmacopsychiatry. 2017;50(6):248–55.

    Article  PubMed  CAS  Google Scholar 

  47. Aichhorn W, Marksteiner J, Walch T, Zernig G, Saria A, Kemmler G. Influence of age, gender, body weight and valproate comedication on quetiapine plasma concentrations. Int Clin Psychopharmacol. 2006;21:81–5.

    Article  PubMed  Google Scholar 

  48. Castberg I, Skogvoll E, Spigset O. Quetiapine and drug interactions: evidence from a routine therapeutic drug monitoring service. J Clin Psychiatry. 2007;68:1540–5.

    Article  PubMed  CAS  Google Scholar 

  49. Bakken GV, Rudberg I, Molden E, Refsum H, Hermann M. Pharmacokinetic variability of quetiapine and the active metabolite N-desalkylquetiapine in psychiatric patients. Ther Drug Monit. 2011;33:222–6.

    Article  PubMed  CAS  Google Scholar 

  50. den Elzen WPJ, Willems JM, Westendorp RGJ, et al. Effect of anemia and comorbidity on functional satus and mortality in old age: results from the Leiden 85-plus Study. CMAJ. 2009;181:151–7.

    Article  Google Scholar 

  51. Riva E, Tettamanti M, Mosconi P, et al. Association of mild anemia with hospitalization and mortality in the elderly: the health and anemia population-based study. Haematologica. 2009;94:22–8.

    Article  PubMed  Google Scholar 

  52. Patel KV. Epidemiology of anemia in older adults. Semin Hematol. 2008;45:210–7.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Peters R, Burch L, Warner J, et al. Haemoglobin, anemia, dementia and cognitive decline in the elderly, a systematic review. BMC Geriatr. 2008;8:18.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Katsogiannou E, Athanasiou L, Christodoulopoulos G, Plozopoulou Z. Diagnostic approach of anemia in ruminants. J Hell Vet Med Soc. 2018;69:1033–46.

    Article  Google Scholar 

  55. Beghe C, Wilson A, Ershler WB. Prevalence and outcomes of anemia in geriatrics: a systematic review of the literature. Am J Med. 2004;116(Suppl 7A):3S-10S.

    Article  PubMed  Google Scholar 

  56. Wilmington DE. Product Information: SEROQUEL(R) oral tablets, quetiapine fumarate oral tablets. AstraZeneca Pharmaceuticals LP (per FDA). 2013.

  57. Wilmington DE. Product Information: SEROQUEL XR(R) oral extended-release tablets, quetiapine fumarate oral extended-release tablets. AstraZeneca Pharmaceuticals LP (per FDA). 2013.

  58. Zhao Y, Wen SW, Li M, et al. Dose-response association of acute-phase quetiapine treatment with risk of new-onset hypothyroidism in schizophrenia patients. Br J Clin Pharmacol. 2021;87(12):4823–30.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported by Hebei Provincial Department of science and technology in China (Grant Number 22377782D) and Medical science research project of Hebei Health Commission (Grant Number 20221434).

Author information

Authors and Affiliations

Authors

Contributions

XP and JY planned and designed the study. YH, JZ, LY analyzed data. CZ and ZY organized and coordinated the work. YH drafted the manuscript. YH and XH revised the manuscript.

Corresponding authors

Correspondence to Fei Gao or Chunhua Zhou.

Ethics declarations

Ethics approval and consent to participate

The Hospital Ethics Committee of Hebei Medical University's First Hospital (20220403) approved the study protocol, which was performed in compliance with the Helsinki Declaration. The ethics committee waived the requirement of written informed consent for the participants because study data have been fully deidentified, and confidential information of patients has been deleted, in accordance with the CIOMS/WHO International Ethical Guidelines for Health-related Research Involving Humans (2016).

Consent for publication

Not applicable.

Competing interests

None.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Figure S1. Boxplot of different doses.

Additional file 2.

Table S1. Parameters of nine models.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hao, Y., Zhang, J., Yu, J. et al. Predicting quetiapine dose in patients with depression using machine learning techniques based on real-world evidence. Ann Gen Psychiatry 23, 5 (2024). https://0-doi-org.brum.beds.ac.uk/10.1186/s12991-023-00483-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12991-023-00483-w

Keywords