Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 21 October 2025

Sec. Intensive Care Medicine and Anesthesiology

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1670994

Increased propofol consumption with later anesthesia start times in sedated gastrointestinal endoscopy: insights from regression and machine learning models

Qiong Lan&#x;Qiong Lan1Zhuonan Sun&#x;Zhuonan Sun1Tian Wang&#x;Tian Wang1Zhuya HuangZhuya Huang2Dengyang HanDengyang Han1Taotao LiuTaotao Liu1Hua ZhangHua Zhang3Ye WangYe Wang4Rui ZhangRui Zhang1Binlong LiBinlong Li1Ning YangNing Yang1Yinyin QuYinyin Qu1Huili Liu
Huili Liu1*Mao Xu
Mao Xu1*
  • 1Department of Anesthesiology, Peking University Third Hospital, Beijing, China
  • 2School of Electronic Engineering, Beijing University of Posts and Telecommunications, Beijing, China
  • 3Research Center of Clinical Epidemiology, Peking University Third Hospital, Beijing, China
  • 4Department of Gastroenterology, Peking University Third Hospital, Beijing, China

Background: Chronopharmacology is an important but underexplored aspect of propofol administration. Despite the implementation of propofol administration models, none have yet incorporated temporal variables. This study aims to investigate the impact of temporal variations on propofol administration during sedated gastrointestinal (GI) endoscopy. Moreover, we aim to develop regression models to predict manually-controlled propofol administration that integrate temporal variables.

Methods: This prospective single-center cohort study enrolled patients undergoing sedated GI endoscopy. For analysis, patients were categorized into 4 groups based on the anesthesia start time: Group 1 (8:00–10:00), Group 2 (10:00–12:00), Group 3 (13:00–15:00), and Group 4 (15:00–17:00). Perioperative characteristics and propofol doses were compared across groups. Correlation analysis was conducted to evaluate the relationship between propofol dose and the anesthesia start time. Subsequently, linear regression models were developed for manually-controlled propofol administration.

Results: A total of 146 cases were included in the statistical analysis. Significant differences were found for all parameters related to propofol dose across the 4 different groups, including induction dose, maintenance dose, total dose and these doses per kilogram per hour. Furthermore, there were positive correlations between the anesthesia start time and all parameters. In the linear regression models, the induction dose equation incorporated the anesthesia start time, age and weight as variables. The model of the maintenance dose per kilogram per hour included the anesthesia start time, duration and weight as variables.

Conclusion: The results suggest that propofol dose increases with later anesthesia start times. Therefore, further clinical administration of propofol should incorporate a heightened consideration of temporal factors.

Trial registration: This prospective study has been registered in the Chinese Clinical Trial Registry (Registration date: December 3, registry number ChiCTR2400093328).

Introduction

Over the past decade, there has been a substantial increase in the global population undergoing sedated examinations and procedures, including but not limited to sedated endoscopy, hysteroscopy, and oocyte retrieval. In China, about 7,000,000 patients undergo sedated gastrointestinal endoscopy GI annually (1). The advantages of sedated procedures, such as rapid recovery, enhanced comfort, and early discharge have significantly improved patient satisfaction and compliance (2, 3). Correspondingly, the precise administration of sedatives and anesthetics is crucial for optimizing procedural efficiency and patient turnover (4, 5).

Propofol is the most extensively applied anesthetic in sedated examinations and procedures (6). Various factors influence the dose of propofol, including age (7), gender (7), body weight (7, 8), circadian rhythm (9, 10), insomnia (11), anxiety (12), smoking (13), and others. Recently, the impact of circadian rhythm on propofol administration has attracted increasing attention (14).

Circadian rhythms are endogenous, entrainable oscillations of physiological structure and function that exhibit a periodicity of approximately 24 h. Studies have demonstrated that circadian rhythm influences propofol administration although these studies have limitations. For instance, in patients undergoing laparoscopic surgeries, the Narcotrend index, mean arterial pressure (MAP) and heart rate (HR) during procedures between 08:00 and 18:00 were observed to be higher than those during procedures between 22:00 and 5:00 (10). Another study demonstrated that patients undergoing elective laparoscopic abdominal surgeries between 8:00 and 12:00 required a higher dose of propofol compared to those undergoing surgeries between 18:00 and 22:00 (15). These studies compared surgeries conducted during daytime hours with those performed at night, thereby limiting their applicability. Another study illustrated that gynecological outpatient short-duration surgeries (abortion and curettage) performed between 8:30 and 11:30 required a lower dose of propofol compared to surgeries conducted between 14:00 and 17:00 (16). All these qualitative studies confirmed the temporal variation in propofol administration. However, because they lack quantitative computation, they were unable to offer practical guidance for clinical practice in propofol administration.

Therefore, this study aims to investigate the variations in propofol dose among patients undergoing sedated GI endoscopy during daytime hours and to construct a feasible model of propofol administration that includes circadian rhythm considerations. The anticipated outcomes are expected to facilitate faster patient turnover and recovery, as well as reduce adverse reactions among outpatients.

Materials and methods

A prospective, single-center cohort study was conducted at the Endoscopic Ambulatory Surgery Center of Peking University Third Hospital from December 2024 to March 2025. The primary aim of this investigation was to determine whether temporal variations exist in the dose of propofol administered during sedated GI endoscopy. A secondary objective was to examine the variables influencing propofol dose in patients undergoing sedated GI endoscopy during daytime hours and to develop viable models for manually-controlled propofol administration that incorporate circadian rhythm considerations.

Ethical statement

This study was conducted in compliance with the Helsinki Declaration and with the approval of the Ethics Committee at Peking University Third Hospital (Grant number M2024956). Furthermore, this prospective study has been registered in the Chinese Clinical Trial Registry (registry number ChiCTR2400093328).

Study population

The study sought to recruit adult patients (18–65 years) with an American Society of Anesthesiologists (ASA) classification of I or II. ASA I patients were otherwise healthy except for the surgery, and ASA II patients had with mild systemic disease and good tolerance. Patients classified as ASA III and IV, who have severe systemic disease and substantive functional limitations or a constant threat to life, respectively, were not included. These patients underwent diagnostic esophagogastroduodenoscopy and colonoscopy under general anesthesia without intubation.

Informed consent was obtained from all patients prior to their inclusion in the study. Each participant was provided with a thorough explanation of the study’s objectives and nature and received a signed copy of the consent document.

The exclusion criteria were as follows: documented allergic reaction to propofol or opioids, history of opioids or propofol misuse, administration of general anesthesia within the preceding 30 days, the necessity for therapeutic procedures during GI endoscopy, anesthesia duration exceeding 60 min, incomplete colonoscopy due to inadequate bowel preparation, and unwillingness to participate in the research.

Data collection

Pre-operative baseline collection

A total of 162 patients scheduled for GI endoscopy were recruited. Patients were required to undergo bowel preparation and were prescribed a laxative solution before GI endoscopy. The bowel preparation protocol applied in our hospital involved a divided administration method: half of the solution was taken the night before GI endoscopy, and the second half was taken early in the morning on the day of the procedure, between 2:00 to 7:00 a.m. The administration time of the second half of the solution depended on the specific situation of each patient.

Before GI endoscopy, patients received a preoperative interview to collect baseline data. This data included demographic information such as age, gender, height, and weight, as well as a medical history of chronic cardiovascular and cerebrovascular diseases, including hypertension, diabetes, cerebrovascular disease, hyperlipidemia, and coronary heart disease. Additionally, the interview gathered a history of other factors history that may affect propofol dosage, including depression, anxiety, insomnia, recent smoking, and alcoholism. GI endoscopy were conducted consecutively from 8:00 to 12:00 in the morning and 13:00 to 17:00 in the afternoon, except for a one-hour lunch break from 12:00 to 13:00.

Perioperative management

The administration of general anesthesia was conducted by two anesthesiologists (Q, L and Z, S), in accordance with institutional protocols and current guidelines (2, 3, 6). The intravenous catheterization was placed in the preparation room prior to GI endoscopy, followed by the infusion of Ringer’s lactate solution. Subsequently, prior to gastroscopy, dyclonine hydrochloride (10 mL: 0.1 g) was administered orally for mucosal anesthesia and lubrication.

Upon entry into the operating room and throughout the entire procedure, patients underwent continuous monitoring of HR, heart rhythm, oxygen saturation (SpO2), non-invasive blood pressure (NIBP), and the bispectral index (BIS). These parameters were recorded at 5 time points: before induction (T0), after the gastroscopy traversed the esophageal orifice (T1), after the gastroscopy passed through the cardia (T2), at the commencement of colonoscopy (T3) and before departure from the operating room (T4). Patients received oxygen via a nasal cannula at a flow rate of 5 L/min while positioned in the left lateral decubitus position.

The administration of propofol for GI endoscopy involved two phases: an initial bolus dose for induction, followed by a continuous infusion facilitated by a syringe pump. To determine the accurate dose, propofol administration adhered to a strategy of starting with a relatively low dose within the recommended range and gradually increasing it to meet clinical requirements (7, 17). Specifically, for induction, propofol was titrated to the onset of hypnotic effect with an infusion rate of 800 mL/h, alongside fentanyl at a dosage of 0.001 mg/kg (17). Following this, for maintenance, propofol was continuously infused at a rate of 4–5 mg/kg/h to maintain a BIS between 40 and 65 (18, 19). Additional doses of propofol (0.5–1 mg/kg) based on the extent was administered each time the patient exhibited a BIS ≥65 and/or physical reflexes, such as a gag reflex, coughing or body movement. If two consecutive supplementary doses of propofol were required, the maintenance rate of propofol was increased by 0.5 mg/kg/h (20). In cases when three administrations of propofol were insufficient to suppress the aforementioned reflexes, an additional dose of fentanyl was administered, and the patient was excluded from the data analysis to eliminate the confounding effect of varying fentanyl doses on propofol administration. The infusion rate of propofol for maintenance would be reduced by 0.5–1 mg/kg/h if BIS < 40.

All endoscopic procedures were carried out by experienced endoscopists who possess a minimum of 5 years of expertise in performing colonoscopies (21). Propofol infusion was discontinued immediately after the colonoscope reached the ileum and preparations for withdrawal began. Following removal of the endoscope from the anus, the attending clinician gently tapped patients on the shoulder and verbally prompted them to awaken. The awakening time was defined as the interval between the cessation of propofol infusion and the attainment of a modified Aldrete score >8. Once a modified Aldrete score of >8 was achieved, patients were transferred to the post-anesthesia care unit for further recovery.

The adverse reactions were defined as SpO2 below 90%, HR below 50 beats per minute (bpm), systolic blood pressure (SBP) below 90 mmHg or BIS below 40. When SpO2 dropped below 90%, we performed interventions such as the jaw-thrust maneuver, inserted an oropharyngeal or nasopharyngeal airway, or assisted ventilation via a facial mask according to the procedural requirements and the patient’s condition. In cases where the HR dropped below 50 bpm, atropine (0.5 mg) was administered (20) when SBP fell below 90 mmHg, ephedrine 5–10 mg was given and 200 mL Ringer’s lactate was also administered (20). The maintenance infusion rate of propofol would be reduced by 0.5–1 mg/kg/h if adverse reactions recurred after the initial event and intervention.

Sample-size estimation

The study was designed to compare induction and maintenance doses across 4 distinct periods: Group 1 (8:00–10:00), Group 2 (10:00–12:00), Group 3 (13:00–15:00), and Group 4 (15:00–17:00), categorized by anesthesia start time. Due to clinical constraints, however, the distribution of participants among different groups was uneven, with a higher number of patients allocated to Group 1 and Group 3. In our pilot study, the maintenance doses per kilogram per hour were 6.35 ± 2.09 mg/kg/h and 7.51 ± 2.00 mg/kg/h (n = 10) for Groups 1 and 3, respectively. To achieve statistical significance with a two-sided α value of 0.05 and 85% power, while accounting for a 10% drop-out rate, we determined that a sample size of 55 participants was necessary in each of Groups 1 and 3.

Statistical analysis

The statistical analysis was performed using SPSS 25.0 software. The Kolmogorov–Smirnov test was used to assess the normality of the data distribution. Categorical variables were presented as numbers and percentages, while continuous variables were reported as mean ± standard deviation (SD) or median ± interquartile range (IQR). A one-way analysis of variance (ANOVA) was used to compare the 4 groups differentiated by varying initial times of anesthesia induction, followed by the Least Significant Difference (LSD) test for post-hoc pairwise comparisons. For categorical variables, the χ2 test or the Fisher’s exact test was applied. Statistical significance was set at a two-sided p-value < 0.05.

The correlations between the anesthesia start time and all 6 parameters of propofol dosage were analyzed using linear regression. The association between the number of patients requiring supplemental administration of propofol and the anesthesia start time was examined using a two-way ordered chi-square test with Kendell’s tau-b.

To develop linear regression models for manually-controlled propofol administration, the dataset was divided into constructing and testing subsets using a random 80%:20% split. The constructing subset was used to develop the predictive model, while the testing subset was used to evaluate the model’s performance in predicting propofol administration. For the regression analysis, univariate linear regression was initially employed to identify potential risk factors for the induction dose and maintenance doses per kilogram per hour. All candidate variables with a p-value <0.2 were subsequently included in multiple linear regression analysis. These candidate variables were evaluated for collinearity and, if collinearity was present, the variables were selected according to clinical significance. Following this, multivariable logistic regressions were performed to further analyze the data.

Machine learning regression

The dataset was also divided into training and testing subsets using a random 80%:20% split. Two separate neural networks were constructed for predicting the induction dose and the maintenance dose per kilogram per hour, respectively.

For the output induction dose, a neural network comprised of 3 layers was constructed: the first layer contained 64 neurons, the second layer 32 neurons, and the last layer produced a single output value. The Rectified Linear Unit (ReLU) was used as the activation function, and a 20% dropout layer was incorporated after the first two layers to mitigate overfitting. The Adam optimizer was used for model training, with a learning rate of 0.001 and L2 regularization (weight decay = 1e-5).

To predict the maintenance dose per kilogram per hour, a neural network comprised of 4 layers was constructed. The first layer consisted of 128 neurons; the second layer, 64 neurons; the third layer, 32 neurons; and the final layer produced a single output value. The ReLU activation function was applied after each fully connected layer to introduce nonlinearity. Additionally, a dropout rate of 0.01 was implemented after the first and second layers to reduce overfitting. The RMSprop optimizer was used for parameter optimization, with a learning rate of 0.005, and L2 regularization (weight decay = 1e-4).

Both models were trained over a total of 300 epochs, using the mean squared error (MSE) as the loss function. The training loss was recorded at each epoch. An early stopping strategy was applied based on the validation set MSE. If the validation loss did not decrease for several consecutive epochs, training was terminated early, and the model parameters with the best validation performance were saved.

Comparison of the predictive performance metrics among regression models

The predictive performances of the multiple linear and machine-learning models were compared using mean absolute error (MAE), root mean square error (RMSE), and the coefficient of determination (R2). The optimal model is characterized by the highest R2 and lowest errors (22).

Results

Demographic characteristics and perioperative conditions

A total of 162 patients were initially screened for the study. Ultimately, 146 patients were eligible for data analysis (Figure 1). To determine whether variations exist in the dosage of propofol during sedated GI endoscopy at different time intervals, the data were categorized into four groups for analysis, based on the anesthesia start time: Group 1 (8:00–10:00), Group 2 (10:00–12:00), Group 3 (13:00–15:00), and Group 4 (15:00–17:00). The demographic and baseline characteristics of the patients are shown in Table 1. No significant differences in baseline characteristics were observed across the different groups.

Figure 1
Flowchart depicting patient selection for a study on painless gastrointestinal procedures using propofol. It starts with 162 adult patients (ages eighteen to sixty-five). Two were excluded for not signing informed consent, resulting in 160 patients enrolled. Fourteen more were excluded due to various reasons: extra fentanyl use (one), transfer to therapeutic procedures (four), inadequate bowel preparation (three), discontinued medication due to adverse events (three), and records with missing data (three). Finally, 146 patients were eligible for data analysis.

Figure 1. Study flow chart.

Table 1
www.frontiersin.org

Table 1. Demographic and baseline characteristics.

Hemodynamic parameters (HR, SpO2, BIS, SBP, and MAP) and adverse reactions are detailed in Supplementary Material 1 and Table 2, respectively. No significant differences in hemodynamic parameters (Supplementary Material 1) or adverse reactions (Table 2) were observed across groups except for the baseline SBP measured prior to the anesthesia. These results indicate a comparable depth and efficacy of anesthesia among the groups.

Table 2
www.frontiersin.org

Table 2. Hemodynamic adverse reactions during anesthesia.

Propofol dosage increased with a later initial time of anesthesia induction

To determine the temporal variation in propofol dosage, we analyzed the duration and dosage of propofol administration, as well as the awakening time after the cessation of propofol infusion (Table 3). The results indicated no significant differences in the duration of propofol infusion or the awakening time after cessation of propofol infusion (Table 3). Together with the absence of significant differences in hemodynamic parameters and adverse reactions, as previously noted, these results suggest comparable GI endoscopy and anesthesia efficacy across different time groups.

Table 3
www.frontiersin.org

Table 3. Intraoperative condition of propofol administration.

Next, we compared the dose of propofol across different time groups. Significant differences were observed among the groups for all 6 parameters of propofol dosage: total dose, propofol dose per kilogram per minute, induction dose, induction dose per kilogram, maintenance dose, and maintenance dose per kilogram per hour (Table 3). Pairwise comparisons revealed that more pronounced differences were found between groups with larger intervals in anesthesia start time while smaller differences were observed between groups with adjacent anesthesia start time. For example, all 6 parameters were significantly higher in Group 3 compared to Group 1 (Figures 2AF), and 3 parameters were significantly higher in Group 4 compared to Group 1 (Figures 2A,C,D). When comparing Group 2 with Group 3 (Figures 2A,C,D) or Group 4 (Figures 2A,C,D), significant differences were observed in 3 parameters. One parameter demonstrated significant differences when comparing Group 1 with Group 2 (Figure 2F) or Group 3 with Group 4 (Figure 2C). These results indicate an increasing trend in propofol dosage with later anesthesia start times.

Figure 2
Bar charts A to F compare propofol administration metrics across groups G1 to G4, corresponding to different time slots: 8:00-10:00, 10:00-12:00, 13:00-15:00, and 15:00-17:00. Significant differences are indicated with p-values. Chart A shows total propofol infusion amounts; B: infusion rates per kilogram per hour; C, induction doses; D, induction doses per kilogram; E, maintenance doses; and F, maintenance dose per kilogram per hour. Error bars indicate variation.

Figure 2. Parameters of propofol dose trends to increase with later anesthesia start times. The parameters include: (A) total propofol infusion dose; (B) propofol infusion rate; (C) induction dose of propofol; (D) induction dose of propofol per kilogram; (E) maintenance dose; (F) maintenance dose per kilogram per hour.

To confirm the correlation between the anesthesia start time and propofol dose, correlation analysis was conducted. The results demonstrated a positive correlation between the anesthesia start time and all six parameters of propofol dose: total propofol dose (r2 = 0.06, p = 0.003) (Figure 3A), propofol dose per kilogram per hour (r2 = 0.06, p = 0.004) (Figure 3B), induction dose (r2 = 0.05, p = 0.007) (Figure 3C), induction dose per kilogram (r2 = 0.05, p = 0.007) (Figure 3D), maintenance dose (r2 = 0.03, p = 0.047) (Figure 3E), and the maintenance dose of propofol per kilogram per hour (r2 = 0.04, p = 0.02) (Figure 3F).

Figure 3
Six scatter plots labeled A to F demonstrate relationships between anesthesia start time and various measures of propofol dosage. Each plot displays a positive trend with R squared values, indicating statistically significant positive correlations (p-values below 0.05). The panels represent: A) Total propofol dose, B) Propofol infusion rate per kilogram per hour, C) Induction dose of propofol, D) Induction dose per kilogram, E) Maintenance dose, and F) Maintenance dose per kilogram per hour.

Figure 3. A positive correlation between propofol dosage and the anesthesia start time. The data discontinuity observed between 12:00–13:00 corresponded to the scheduled lunch break. (A) Total propofol dose; (B) propofol infusion rate; (C) induction dose of propofol; (D) induction dose of propofol per kilogram; (E) maintenance dose; (F) maintenance dose per kilogram per hour.

The percentage and frequency of patients requiring supplementary dosing increased with later anesthesia start time

Our propofol administration strategy, initiating with a relatively low dose and gradually increasing it, facilitated accurate propofol delivery and enable the investigation of propofol dose, although it was accompanied by increased BIS and/or physical responses. Among 146 patients, 95 (65.1%) required additional propofol administration. Specifically, 55 patients (57.9%) required one additional administration of propofol, 27 patients (28.4%) required two, and 13 patients (13.7%) required three. Group 1 had the highest percentage of patients who did not require extra dosing and the lowest percentage of patients who required 3 times of extra dosing. Later anesthesia start times, the percentage and frequency of patients requiring extra dosing of propofol increased (p = 0.038) (Figure 4). The mean dose of supplementary administration did not demonstrate statistical significance among the 4 groups (Supplementary Material 2).

Figure 4
Four pie charts depict the distribution of patients requiring 0 to 3 supplemental doses across the four time groups. A clear trend is observed: the proportion of patients requiring no supplemental doses decreases, while the proportion requiring more frequent dosing (2 or 3 times) increases as the time of day progresses from Group 1(08:00-10:00) to Group 4 (15:00-17:00).

Figure 4. The percentage and frequency of patients requiring extra propofol doses increased with a later anesthesia start time. Values of 0, 1, 2, and 3 represent the number of additional propofol doses required. Panels (A–D) represent Groups 1, 2, 3, and 4, respectively. Panel (E) presents the number and percentage of patients requiring additional doses among the four groups.

The multiple linear regression models exhibited a positive correlation between the start time and propofol dose

Previous studies have developed models to predict manually-controlled propofol administration in GI endoscopy (23). However, to our knowledge, no prior study has examined the impact of circadian rhythm on these regression models (24). Our results confirmed the influence of the anesthesia start time on propofol administration. We subsequently developed 2 sets of regression models for manually-controlled propofol administration that incorporated this variable. Based on the practice of manually-controlled propofol administration, the models predicted two dependent variables: the induction dose and the maintenance dose per kilogram per hour.

First, linear regression models were developed using SPSS. For the induction dose, a univariable linear regression analysis was performed (Supplementary Material 3), followed by an assessment of collinearity (Supplementary Material 4). The analysis identified 5 impact factors for inclusion in the multiple linear regression model: age, gender, weight, the anesthesia start time and the presence of chronic diseases. Consequently, the multiple linear regression was represented by the following equation: Induction dose per kilogram = 53.113–0.447 × Age+ 0.656 × Weight + 2.003 × The anesthesia start time (R2 = 0.34, p < 0.001) (Table 4 and Figure 5).

Table 4
www.frontiersin.org

Table 4. Multiple linear regression for the induction dose as the dependent variable.

Figure 5
Scatter plot showing the relationship between induction dose (mg) and three variables: weight in kilograms (blue dots, upward trend), age in years (red squares, downward trend), and anesthesia start time (black triangles, flat trend). The coefficient of determination is R² = 0.34, with a p-value less than 0.001.

Figure 5. Graphical representation of the linear regression model of the induction dose. The anesthesia start time, originally recorded as hh:mm, was converted into decimal hours using the following formula: decimal time = hh + (mm/60).

For the maintenance dose per kilogram per hour, 6 impact factors were incorporated into the multiple linear regression model based on the univariable linear regression (Supplementary Material 5) and the collinearity analysis (Supplementary Material 6). These factors included gender, weight, other relevant factors affecting propofol dose, the anesthesia start time, duration and the induction dose per kilogram. Consequently, the multiple linear regression model was represented by the following equation: Maintenance dose per kilogram per hour = 11.455–0.074 × Weight −0.076 × Duration + 0.172 × The anesthesia start time (R2 = 0.21, p < 0.001) (Table 5 and Figure 6).

Table 5
www.frontiersin.org

Table 5. Multiple linear regression for the dependent variable of maintenance dose per kilogram per hour.

Figure 6
Scatter plot depicting the relationship between maintenance dose per kilogram per hour and three variables: weight in kilograms (blue), duration in minutes (red), and anesthesia start time in converted clock time (black). The coefficient of determination is R² = 0.21, with a p value less than 0.001.

Figure 6. Graphical representation of the linear regression model of the maintenance dose per kilogram per hour. The anesthesia start time, originally recorded as hh:mm, was converted into decimal hours using the formula: decimal time = hh + (mm/60).

Machine learning regression models

Machine learning regression models were also developed. Extensive model experiments were performed (Supplementary Material 7). Given the volume of data and analysis requirements, the selected models demonstrated optimal suitability for this study by exhibiting lower errors and higher R2 values. For predicting the induction dose, the model’s input features included: age, gender, height, weight, ASA classification, presence of any chronic diseases, other relevant factors and the time of anesthesia induction. For predicting the maintenance dose per kilogram per hour, the input features of the predictive model included: age, gender, BMI, ASA classification, presence of any chronic diseases, other relevant factors, duration and the time of anesthesia induction. The statistical parameters of the models (RMSE, MAE, R2) were presented in Table 5. The feature importance and the training loss curves were illustrated in Supplementary Material 8.

Comparison of the linear regression models and machine learning regression models

The linear regression models were compared with the machine learning regression models. For predicting induction dose, the linear regression model demonstrated enhanced predictive performance compared to the machine learning model, characterized by lower RMSE, MAE and higher R2. Conversely, for predicting the maintenance dose per kilogram per hour, the machine learning model demonstrated superior performance compared to the linear regression model based on RMSE, MAE, and R2 parameters (Table 6).

Table 6
www.frontiersin.org

Table 6. Comparison of the linear regression models and machine learning regression models.

Discussion

Chronopharmacology represents a crucial yet insufficiently explored dimension of propofol administration. Previous studies have investigated the influence of circadian rhythm on propofol administration; however, these studies had limitations, such as small sample sizes and lack of clear clinical practice (10, 15, 16). This cohort study prospectively collected data from 162 cases during sedated GI endoscopy using propofol, confirmed the temporal variation in propofol administration, and developed regression models for manually-controlled propofol administration in these contexts.

Previous researches have indicated that under target-controlled infusion (TCI) of propofol, BIS was significantly lower in the nocturnal group (22:00–2:00) compared to the diurnal group (8:00–12:00) (25). Furthermore, when BIS values were maintained at comparable levels, a greater dose of propofol was required during daytime administration than at nighttime 10. Nonetheless, these studies offer limited clinical guidance, as the majority of surgeries are conducted during daytime hours. There is a paucity of research examining temporal variations in propofol administration during the day. One study indicated that compared with short-duration gynecological procedures conducted during 8:30–11:30, propofol dose was higher in procedures undergoing during 14:00–17:00 (16). Our study demonstrates that, within a comparable population and consistent GI endoscopy stimulation across different anesthesia start time, the dose of propofol increased with later anesthesia start times, consistent with previous research.

Chronopharmacology can be divided into two primary areas: chronopharmacokinetics and chronopharmacodynamics. In the context of pharmacokinetics, propofol undergoes hepatic metabolism (17). The hepatic blood flow (26) and cytochrome P450 monooxygenase (5, 27), the enzyme responsible for the metabolism of numerous drugs including propofol, exhibit circadian fluctuations. Studies have shown that P450 monooxygenase activity increased at night and decreased during the day in a rat model (28). These results suggest that the pharmacokinetics of propofol exhibit a circadian rhythm. Within the realm of pharmacodynamics, the central nervous system depression induced by propofol is mediated through gamma--aminobutyric acid (GABA) (17, 29) on GABAergic neurons. Studies have shown that the morphology and function of neurons are subject to circadian rhythm. For instance, variations in the dendritic structure and spine density of neurons in the rat frontal cortex have been observed to follow circadian patterns (30, 31). Another study showed that GABAergic activity peaked at nocturnal hours and reached a trough at diurnal hours in the cerebral cortex of Syrian hamsters (32). These studies support the impact of circadian rhythm on neuronal activity targeted by propofol.

To our knowledge, no existing models of propofol administration, whether manually-controlled or TCI, have been developed that incorporate circadian rhythm factors. In this study, linear regression models and machine learning algorithms were developed to facilitate manually-controlled propofol administration. The linear regression model for induction dose included age, weight and the anesthesia start time as influencing factors, with age and weight exerting a more substantial impact on propofol dosage than the anesthesia start time. However, age was not included in the model of the maintenance dose per kilogram per hour. In that model, the duration exerted a relatively smaller influence compared to the anesthesia start time. These models align with the clinical practice and experience.

According to expert consensus (7, 33) and the authors’ experience, anesthesiologists consider age as a crucial factor influencing propofol induction. Consequently, they tend to select a lower induction dose per kilogram within the recommended range (1.5–2.5 mg/kg) for elderly patients. However, age appears to be a less critical consideration when determining the maintenance dose per kilogram per hour. Regarding the maintenance dose, it is widely accepted that, due to the redistribution of propofol, the maintenance dose decreases as the duration of administration extends (17). These findings were all consistent with our models. However, the influence of circadian rhythm on propofol dosing has not been thoroughly investigated before. Our models indicate that the timing of anesthesia induction exerts a more significant impact than the duration of administration. These results underscore the need to consider circadian rhythm in determining propofol dosing.

However, this study has certain limitations. First, due to clinical constraints, our data collection was restricted to specific time periods. Previous animal studies have indicated that anesthetic pharmacodynamics may adhere to a Cosine pattern rather than a linear model when data are collected over 24 h. For instance, one rat study demonstrated that the longest pharmacological effect (the loss of righting reflex) of propofol injection occurred shortly afternoon during the rest phase and the shortest effect shortly after midnight during the active phase (9). In our study, the maintenance dose of propofol per kilogram per hour during 15:00–17:00 was lower than that during 13:00–15:00, which seems to be consistent with the former study. To determine whether propofol administration follows a cosine pattern, further studies are required.

Additionally, the models were developed from data from patients undergoing sedated GI endoscopy. The applicability of these models to anesthetic management and surgeries exceeding 1 h in duration requires further investigation. Nonetheless, the present models provide a quantitative framework for propofol dosing in sedated GI endoscopy and highlight the potential influence of circadian rhythm on drug requirements.

Another limitation may involve sleep deprivation induced by the bowel preparation protocol, which required patients to ingest the second solution early in the morning, potentially disrupting sleep. Although animal studies suggest that sleep deprivation can alter propofol consumption (34, 35), this has not been confirmed in human clinical studies. Future controlled trials are needed to isolate this potential confounding factor.

Finally, due to the limited sample size in Groups 2 and 4, some parameter estimates may be underpowered. That said, the primary conclusions of this study are supported by other inter-group comparisons with sufficient statistical power (all >0.8), indicating that the overall findings are robust.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

The studies involving humans were approved by Ethics Committee of Scientific Research, Peking University Third Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

QL: Formal analysis, Methodology, Writing – original draft. ZS: Formal analysis, Methodology, Writing – original draft. TW: Data curation, Project administration, Writing – review & editing. ZH: Formal analysis, Software, Writing – original draft. DH: Writing – review & editing. TL: Writing – review & editing. HZ: Methodology, Writing – review & editing. YW: Methodology, Writing – review & editing. RZ: Writing – review & editing. BL: Writing – review & editing. NY: Writing – review & editing. YQ: Funding acquisition, Writing – review & editing. HL: Supervision, Writing – review & editing. MX: Funding acquisition, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by Wu Jieping Medical Foundation (320.6750.2024-05-17); National Natural Science Foundation of China (grant number: 82201337); and the Peking University Third Hospital Returned Scholar Startup Fund (Y85477-01).

Acknowledgments

The authors thank all staff of the Endoscopy Center in Peking University Third Hospital for support of the study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1670994/full#supplementary-material

References

1. Zhou, S, Zhu, Z, Dai, W, Qi, S, Tian, W, Zhang, Y, et al. National survey on sedation for gastrointestinal endoscopy in 2758 Chinese hospitals. Br J Anaesth. (2021) 127:56–64. doi: 10.1016/j.bja.2021.01.028

PubMed Abstract | Crossref Full Text | Google Scholar

2. Sidhu, R, Turnbull, D, Haboubi, H, Leeds, JS, Healey, C, Hebbar, S, et al. British Society of Gastroenterology guidelines on sedation in gastrointestinal endoscopy. Gut. (2024) 73:1–27. doi: 10.1136/gutjnl-2023-330396

PubMed Abstract | Crossref Full Text | Google Scholar

3. Hara, T, Ozawa, A, Shibutani, K, Tsujino, K, Miyauchi, Y, Kawano, T, et al. Practical guide for safe sedation. J Anesth. (2023) 37:340–56. doi: 10.1007/s00540-023-03177-5

PubMed Abstract | Crossref Full Text | Google Scholar

4. Cortínez, LI, De la Fuente, N, Eleveld, DJ, Oliveros, A, Crovari, F, Sepulveda, P, et al. Performance of propofol target-controlled infusion models in the obese: pharmacokinetic and pharmacodynamic analysis. Anesth Analg. (2014) 119:302–10. doi: 10.1213/ane.0000000000000317

Crossref Full Text | Google Scholar

5. Nam, S, Yoo, S, Park, SK, Kim, Y, and Kim, J-T. Relationship between preinduction electroencephalogram patterns and propofol sensitivity in adult patients. J Clin Monit Comput. (2024) 38:1069–77. doi: 10.1007/s10877-024-01149-y

PubMed Abstract | Crossref Full Text | Google Scholar

6. Gotoda, T, Akamatsu, T, Abe, S, Shimatani, M, Nakai, Y, Hatta, W, et al. Guidelines for sedation in gastroenterological endoscopy (second edition). Dig Endosc. (2021) 33:21–53. doi: 10.1111/den.13882

PubMed Abstract | Crossref Full Text | Google Scholar

7. Nishizawa, T, and Suzuki, H. Propofol for gastrointestinal endoscopy. United European Gastroenterol J. (2018) 6:801–5. doi: 10.1177/2050640618767594

PubMed Abstract | Crossref Full Text | Google Scholar

8. Kanaya, A, Sato, T, Fuse, N, Yamaguchi, H, Mano, N, and Yamauchi, M. Impact of clinical factors and UGT1A9 and CYP2B6 genotype on inter-individual differences in propofol pharmacokinetics. J Anesth. (2018) 32:236–43. doi: 10.1007/s00540-018-2470-3

PubMed Abstract | Crossref Full Text | Google Scholar

9. Challet, E, Gourmelen, S, Pevet, P, Oberling, P, and Pain, L. Reciprocal relationships between general (Propofol) anesthesia and circadian time in rats. Neuropsychopharmacology. (2007) 32:728–35. doi: 10.1038/sj.npp.1301081

PubMed Abstract | Crossref Full Text | Google Scholar

10. Shen, JH, Ye, M, Chen, Q, Chen, Y, Zhao, HL, Khan, A, et al. Effects of circadian rhythm on Narcotrend index and target-controlled infusion concentration of propofol anesthesia. BMC Anesthesiol. (2021) 21:215. doi: 10.1186/s12871-021-01445-z

PubMed Abstract | Crossref Full Text | Google Scholar

11. Brown, EN, Lydic, R, and Schiff, ND. General anesthesia, sleep, and coma. N Engl J Med. (2010) 363:2638–50. doi: 10.1056/NEJMra0808281

PubMed Abstract | Crossref Full Text | Google Scholar

12. Uysal, A, Altıparmak, B, Korkmaz Toker, M, Uysal, Aİ, Dede, G, Sezgin, Ç, et al. The effect of preoperative anxiety level on mean platelet volume and propofol consumption. BMC Anesthesiol. (2020) 20:34. doi: 10.1186/s12871-020-0955-8

PubMed Abstract | Crossref Full Text | Google Scholar

13. Öztürk, E, Aydoğan, MS, Karaaslan, K, Doğan, Z, and Topuz, U. Does smoking increase the anesthetic requirement? Turk J Med Sci. (2019) 49:1271–6. doi: 10.3906/sag-1602-57

PubMed Abstract | Crossref Full Text | Google Scholar

14. Vullo, PA, Real Navacerrada, MÍI, and Navarro Suay, R. Hemodynamic impact of increasing time between fentanyl and propofol administration during anesthesia induction: a randomised, clinical trial. Braz J Anesthesiol. (2024) 74:744230. doi: 10.1016/j.bjane.2021.07.009

PubMed Abstract | Crossref Full Text | Google Scholar

15. Song, B, Li, Y, Teng, X, Li, X, Yang, Y, and Zhu, J. Comparison of morning and evening operation under general anesthesia on intraoperative anesthetic requirement, postoperative sleep quality, and Pain: A randomized controlled trial. Nat Sci Sleep. (2020) 12:467–75. doi: 10.2147/nss.s257896

PubMed Abstract | Crossref Full Text | Google Scholar

16. Chen, X. Comparison of propofol dosage in morning and at afternoon undergoing gynecological clinical short surgery. J Clin Anesthesiol. (2016) 32:654–6. Available at: https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CJFD&dbname=CJFDLAST2016&filename=LCMZ201607008

Google Scholar

17. Sahinovic, MM, Struys, M, and Absalom, AR. Clinical pharmacokinetics and pharmacodynamics of propofol. Clin Pharmacokinet. (2018) 57:1539–58. doi: 10.1007/s40262-018-0672-3

PubMed Abstract | Crossref Full Text | Google Scholar

18. Johansen, JW, and Sebel, PS. Development and clinical application of electroencephalographic bispectrum monitoring. Anesthesiology. (2000) 93:1336–44. doi: 10.1097/00000542-200011000-00029

PubMed Abstract | Crossref Full Text | Google Scholar

19. Gu, A, Cao, W, Zhang, J, Qin, Z, Chen, B, Huang, P, et al. Effect of genetic polymorphism of CYP2B6 * 6 on pharmacokinetics and pharmacodynamics of single dose propofol in Chinese healthy subjects Chin J Clin Pharmacol (2018) 34:1052–5. doi: 10.13699/j.cnki.1001-6821.2018.09.013

Crossref Full Text | Google Scholar

20. Steur, RJ, Perez, RS, and De Lange, JJ. Dosage scheme for propofol in children under 3 years of age. Paediatr Anaesth. (2004) 14:462–7. doi: 10.1111/j.1460-9592.2004.01238.x

PubMed Abstract | Crossref Full Text | Google Scholar

21. Qiu, D, Wang, XM, Yang, JJ, Chen, S, Yue, CB, Hashimoto, K, et al. Effect of intraoperative Esketamine infusion on postoperative sleep disturbance after gynecological laparoscopy: A randomized clinical trial. JAMA Netw Open. (2022) 5:e2244514. doi: 10.1001/jamanetworkopen.2022.44514

PubMed Abstract | Crossref Full Text | Google Scholar

22. Alboukadel, K. Machine learning essentials: Practical guide in R STHDA (2017). Available at: http://www.sthda.com/en/wiki/machine-learning-essentials-practical-guide-in-r

Google Scholar

23. Rocha, C, Mendonça, T, Silva, ME, and Gambús, P. Individualizing propofol dosage: a multivariate linear model approach. J Clin Monit Comput. (2014) 28:525–36. doi: 10.1007/s10877-013-9510-1

PubMed Abstract | Crossref Full Text | Google Scholar

24. Yun, WJ, Shin, M, Jung, S, Ko, JG, Lee, HC, and Kim, J. Deep reinforcement learning-based propofol infusion control for anesthesia: A feasibility study with a 3000-subject dataset. Comput Biol Med. (2023) 156:106739. doi: 10.1016/j.compbiomed.2023.106739

PubMed Abstract | Crossref Full Text | Google Scholar

25. Zhu, GQ, Yan, L, and You, KZ. Influence of circadian rhythm on depth of propofol anesthesia. Chin JMAP. (2013) 30:1239–42. doi: 10.13748/j.cnki.issn1007-7693.2013.11.017

Crossref Full Text | Google Scholar

26. Lemmer, B, and Nold, G. Circadian changes in estimated hepatic blood flow in healthy subjects. Br J Clin Pharmacol. (1991) 32:627–9. doi: 10.1111/j.1365-2125.1991.tb03964.x

PubMed Abstract | Crossref Full Text | Google Scholar

27. Košir, R, Španinger, K, and Rozman, D. Circadian events in human diseases and in cytochrome P450-related drug metabolism and therapy. IUBMB Life. (2013) 65:487–96. doi: 10.1002/iub.1160

PubMed Abstract | Crossref Full Text | Google Scholar

28. Hirao, J, Arakawa, S, Watanabe, K, Ito, K, and Furukawa, T. Effects of restricted feeding on daily fluctuations of hepatic functions including p450 monooxygenase activities in rats. J Biol Chem. (2006) 281:3165–71. doi: 10.1074/jbc.M511194200

PubMed Abstract | Crossref Full Text | Google Scholar

29. Lingamaneni, R, Birch, ML, and Hemmings, HC Jr. Widespread inhibition of sodium channel-dependent glutamate release from isolated nerve terminals by isoflurane and propofol. Anesthesiology. (2001) 95:1460–6. doi: 10.1097/00000542-200112000-00027

PubMed Abstract | Crossref Full Text | Google Scholar

30. Perez-Cruz, C, Simon, M, Flügge, G, Fuchs, E, and Czéh, B. Diurnal rhythm and stress regulate dendritic architecture and spine density of pyramidal neurons in the rat infralimbic cortex. Behav Brain Res. (2009) 205:406–13. doi: 10.1016/j.bbr.2009.07.021

PubMed Abstract | Crossref Full Text | Google Scholar

31. Qi, XR, Kamphuis, W, and Shan, L. Astrocyte changes in the prefrontal cortex from aged non-suicidal depressed patients. Front Cell Neurosci. (2019) 13:503. doi: 10.3389/fncel.2019.00503

PubMed Abstract | Crossref Full Text | Google Scholar

32. Kanterewicz, BI, Rosenstein, RE, Golombek, DA, Yannielli, PC, and Cardinali, DP. Daily variations in GABA receptor function in Syrian hamster cerebral cortex. Neurosci Lett. (1995) 200:211–3. doi: 10.1016/0304-3940(95)12112-h

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ho, AM, and Mizubuti, GB. Co-induction with a vasopressor "chaser" to mitigate propofol-induced hypotension when intubating critically ill/frail patients-A questionable practice. J Crit Care. (2019) 54:256–60. doi: 10.1016/j.jcrc.2019.09.015

PubMed Abstract | Crossref Full Text | Google Scholar

34. Yan, Z, Ha, L, Chen, H, Xiao, YF, Chen, M, Wu, B, et al. Sleep deprivation alters hepatic UGT1A9 and propofol metabolism in mice. Biochem Pharmacol. (2025) 232:116713. doi: 10.1016/j.bcp.2024.116713

PubMed Abstract | Crossref Full Text | Google Scholar

35. Tung, A, Szafran, MJ, Bluhm, B, and Mendelson, WB. Sleep deprivation potentiates the onset and duration of loss of righting reflex induced by propofol and isoflurane. Anesthesiology. (2002) 97:906–11. doi: 10.1097/00000542-200210000-00024

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: propofol, temporal variations, sedated gastrointestinal endoscopy, manually-controlled administration, machine learning model

Citation: Lan Q, Sun Z, Wang T, Huang Z, Han D, Liu T, Zhang H, Wang Y, Zhang R, Li B, Yang N, Qu Y, Liu H and Xu M (2025) Increased propofol consumption with later anesthesia start times in sedated gastrointestinal endoscopy: insights from regression and machine learning models. Front. Med. 12:1670994. doi: 10.3389/fmed.2025.1670994

Received: 22 July 2025; Accepted: 06 October 2025;
Published: 21 October 2025.

Edited by:

Katherine N. Theken, University of Pennsylvania, United States

Reviewed by:

Ronan Lordan, University of Pennsylvania, United States
Laurent Bairy, CHU UCL Namur Anesthesiologie, Belgium

Copyright © 2025 Lan, Sun, Wang, Huang, Han, Liu, Zhang, Wang, Zhang, Li, Yang, Qu, Liu and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mao Xu, YW5hZUAxNjMuY29t; Huili Liu, bGl1aHVpbGlAYmptdS5lZHUuY24=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.