Background

Propofol causes significant cardiovascular depression and a slowing of neurophysiological activity. However, literature on its effect on the heart rate remains mixed, and it is not known whether cortical slow waves are related to cardiac activity in propofol anesthesia.

Methods

The authors performed a secondary analysis of electrocardiographic and electroencephalographic data collected as part of a previously published study where n = 16 healthy volunteers underwent a slow infusion of propofol up to an estimated effect-site concentration of 4 µg/ml. Heart rate, heart rate variability, and individual slow electroencephalographic waves were extracted for each subject. Timing between slow-wave start and the preceding R-wave was tested against a uniform random surrogate. Heart rate data were further examined as a post hoc analysis in n = 96 members of an American Society of Anesthesiologists Physical Status II/III older clinical population collected as part of the AlphaMax trial.

Results

The slow propofol infusion increased the heart rate in a dose-dependent manner (mean ± SD, increase of +4.2 ± 1.5 beats/min/[μg ml−1]; P < 0.001). The effect was smaller but still significant in the older clinical population. In healthy volunteers, propofol decreased the electrocardiogram R-wave amplitude (median [25th to 75th percentile], decrease of –83 [–245 to –28] μV; P < 0.001). Heart rate variability showed a loss of high-frequency parasympathetic activity. Individual cortical slow waves were coupled to the heartbeat. Heartbeat incidence peaked about 450 ms before slow-wave onset, and mean slow-wave frequency correlated with mean heart rate.

Conclusions

The authors observed a robust increase in heart rate with increasing propofol concentrations in healthy volunteers and patients. This was likely due to decreased parasympathetic cardioinhibition. Similar to non-rapid eye movement sleep, cortical slow waves are coupled to the cardiac rhythm, perhaps due to a common brainstem generator.

Editor’s Perspective
What We Already Know about This Topic
  • Propofol is a commonly used intravenously administered anesthetic. It can cause significant cardiovascular depression and resulting arterial hypotension.

  • However, it is unclear from previous studies whether propofol administration causes bradycardia.

What This Article Tells Us That Is New
  • In this observational study of healthy volunteers, propofol administration resulted in increased mean heart rate that was proportional to propofol concentration. This increased heart rate was more pronounced in younger versus older subjects.

  • Frontal cortical slow waves preferentially occurred coupled to the heart rhythm. Additional research is needed to elucidate the mechanism and role of observed cortical-cardiac coupling.

Propofol is the most widely used intravenous anesthetic hypnotic drug due to its favorable kinetics, low adverse effects incidence, and smooth induction profile.1  It also causes significant cardiovascular depression, manifesting mainly as arterial hypotension.2  However, despite decades of use, the effect of propofol on heart rate (HR) remains controversial. In clinical settings, propofol administration has been reported to carry a risk of bradycardia,1,3  and several texts state that propofol decreases the HR as an accepted fact.4,5  Others, however, find propofol to have no effect on the HR,2,6  and much of the literature, especially in laboratory settings, appears to show significant increases in HR.7–12  Clinical research is complicated by common coadministration of opioids with some bradycardic effects. On theoretical grounds, propofol’s HR effect may be due to modulation of γ-aminobutyric acid (GABA)–mediated neurotransmission to cardiac parasympathetic neurons in the brainstem.11 

Propofol also affects autonomic nervous activity. Because the autonomic nervous system regulates the HR, this can be indexed by HR variability, the beat-to-beat variation in HR (distinct from the mean HR). “High-frequency” (approximately 0.5- to 7-s) fluctuations in HR variability are dominated by parasympathetically mediated respiratory sinus arrythmia and are a measure of vagal tone. Propofol decreases short-term HR variability, in part through this lowered parasympathetic tone.6  Lower-frequency fluctuations (seconds to 10s of seconds) are driven by other, largely sympathetic, factors.

In the brain, propofol causes neuronal hyperpolarization by prolonging GABA-activated chloride channel opening. At the network level, this causes the cortex to switch between high-activity up states and silent down states. This switching can be observed as slow (approximately 1-Hz) waves on the electroencephalogram (EEG).13  As propofol dose increases, power in the slow-wave band (typically 0.5 to 1.5 Hz) saturates, and the thalamocortical system becomes isolated from environmental stimuli.8  By disrupting cortical information processing, slow waves may have a causal role in sustaining unconsciousness.14  Similar slow waves are observed in non–rapid eye movement (REM) sleep,15  where they have been linked to changes in autonomic activity including individual heartbeats.16–18 

In this study, we first performed an advanced secondary analysis of EEG and electrocardiographic data collected in n = 16 healthy volunteers undergoing a slow propofol infusion up to 4 μg/ml estimated effect-site concentration. As effects of propofol may depend on induction speed,19  our ultraslow infusion provides a unique perspective without influences of concomitant medication. We hypothesized that propofol would increase HR and decrease parasympathetic effects, as indexed by high-frequency HR variability. Our secondary aim was to explore a possible link between electrocardiographic activity and the frontal cortical slow waves seen in the EEG. We hypothesized that, like non-REM sleep, slow waves would preferentially occur time-locked to individual heart beats. Finally, to see if our HR findings held clinical validity, we also performed a post hoc analysis of clinical EEG and HR data from n = 96 American Society of Anesthesiologists (ASA) Physical Status II/III AlphaMax study patients.20 

Data Collection

Thirty-two–channel EEG and single-channel electrocardiogram were collected in n = 16 healthy subjects (8 female; age, 28.6 ± 7 yr; table 1) during slowly increasing intravenous infusion of propofol up to an estimated effect-site concentration of 4 µg/ml using the Marsh pharmacokinetic model.8  The experiment was separated into four main periods: 10 min awake, 48 min induction, 10 min peak anesthesia, and 48 min emergence. Informed written consent was obtained from all participants; details of this experiment have been published previously.8 

Table 1.

Summary of Datasets Used in This Study

Summary of Datasets Used in This Study
Summary of Datasets Used in This Study

Data Preprocessing

EEG data preprocessing was carried out with BrainVision Analyzer version 2.1 (BrainProducts GmbH, Germany), custom written MATLAB code (Matlab 2019b, Math Works Inc., USA), and the EEGLAB (v2019.1) analysis toolbox. The EEG and electrocardiographic data were re-referenced to the common average of signals from all EEG channels. This was done as theoretical reasons suggest scalp average to be a robust null reference that decreases volume conduction effects.21  Independent component analysis and bad channel rejection were performed to remove EEG data with blinks and ocular movements. EEG data were band-pass–filtered with a phase-preserving third order 0.5- to 45-Hz Butterworth filter. EEG data were downsampled to 100 Hz and electrocardiogram to 500 Hz.

Time Series Electrocardiogram Analysis

HR, electrocardiogram waveform templates, and R-wave amplitudes were extracted using the BioSPPy toolbox (https://github.com/PIA-Group/BioSPPy/), which uses Hamilton segmentation.22  This was used on each subject to identify individual R-wave peaks, the HR, and electrocardiogram waveform templates, which were subsequently Spearman-correlated to the propofol effect-site concentration.

HR Variability Electrocardiogram Analysis

In order to explore correlates of autonomic activity, standard HR variability metrics were extracted for 5-min segments in each subject using the pyHRV toolbox23  and Spearman-correlated with propofol effect-site concentration at the group level. These HR variability metrics included root-mean-square successive difference between R-wave peaks as well as other frequency domain metrics. Specifically, the frequency domain metrics used were the ratio of low-frequency (largely sympathetic, 0.04- to 0.15-Hz) and high-frequency (parasympathetic, 0.15- to 0.4-Hz) HR variability and the peak frequency in the high-frequency band. The low- to high-frequency ratio is thought to index the balance between sympathetic and parasympathetic activity, with a low low-frequency/high-frequency reflecting parasympathetic dominance, although this simple interpretation has been challenged.24 

Slow-wave Analysis

Slow-wave activity was defined as the Fourier power in the 0.5- to 1.5-Hz band on the frontal Fz (international 10-20 system) channel, as is common convention.8  The correlation between slow-wave activity and HR/effect-site concentration and its significance was found using Spearman correlation in 5-min segments. Individual slow waves were then identified using standard methodology based on amplitude and duration thresholding implemented in the yasa toolbox.15,16,25  In brief, each slow wave had to have amplitude in the 99th percentile of the 0.5- to 4-Hz amplitude and negative duration between 0.25 and 1.25 s. This wider filter is commonly chosen to capture more details of the nonsinusoidal slow wave shape,26  although proposals for multiple slow-wave types exist.27  Slow-wave onset was defined as the initial downward zero crossing, and slow-wave frequency was extracted as the inverse of the slow-wave period.

Corticocardiac Coupling Analysis

Once heartbeats and slow waves were identified, we aimed to test whether heartbeats occurred at preferential times in the slow-wave cycle. For each slow wave detected, the time delay relative to the slow-wave start (initial downward zero crossing) was noted for eight heartbeats closest to it. Eight beats were chosen as this window length fully covers a slow wave. This resulted in eight R-wave to slow-wave intervals, using methodology similar to previous cardiorespiratory analyses.28,29  We wanted to know if electrocardiogram R-wave to EEG slow-wave timings were distributed randomly or in phase with the slow-wave onset. For robustness, this was tested against a surrogate null distribution in several ways. First, we utilized the same method that has previously been used to study cardiorespiratory coupling.30  This method compares the R-wave to slow-wave-1 interval (i.e., the time interval between slow-wave start and the preceding R-wave peak) to a uniformly random null distribution. Starting from the beginning of each subject’s set of R-wave to slow-wave-1 values, we used a moving window of 40 slow waves, and placed the corresponding R-wave to slow-wave-1 intervals in a 10-bin histogram with outer limits of 0 and the mean heart period for that window. From the histogram, the proportional Shannon entropy is calculated as follows:

where Pb is the histogram probability of bin b and N is the number of histogram bins.

During perfect coupling, all R-wave to slow-wave-1 intervals fall into one bin, and SHP = 0. In the absence of coupling, R-wave to slow-wave-1 intervals are distributed randomly, producing maximum entropy with SHP = 1. For each subject, the mean SHP across the whole experiment was computed. To determine a significance threshold, SHP was computed for n = 10,000 surrogate series of 200 random numbers each, drawn from a uniform distribution between 0 and 1 (mean HR of 60 beats/min). The 0.1st percentile was used to indicate significance at the P = 0.001 level (SHP = 0.970). Additional methods testing autocorrelation of the R-wave to slow-wave histogram and more complex surrogates were used to further verify robustness of this result (see Supplemental Digital Content, https://links.lww.com/ALN/D329).

Additionally, for each slow wave identified, ± 2 s of EEG and electrocardiogram activity was saved around the slow-wave start. This was then averaged across slow waves and subjects to reveal any coherent electrocardiogram patterns during a slow wave.

Clinical Dataset Analysis

In order to carry out a preliminary exploration of whether our HR results could be replicated in clinical data, we performed a post hoc analysis of HR and drug concentrations in n = 96 patients collected as part of the AlphaMax study (median age, 74 yr [range, 61 to 86 yr]; 66 men; ASA Physical Status II/III; variety of procedures; table 1). This dataset contained EEG, drug concentration, HR, and demographic data. Unfortunately, individual electrocardiogram waveforms were not available in this dataset, so we could not determine the HR variability or corticocardiac coupling analyses.

The AlphaMax study patients received a standardized desflurane- and fentanyl-based maintenance general anesthesia that was titrated to maximize the EEG alpha power in the intervention group. For each patient, HR and drug concentrations (propofol, fentanyl, desflurane) were sampled—or estimated using population based pharmacokinetic models—every 5 s. The HR was smoothed with a 2-min moving median window to suppress artifacts, and any HR greater than 250 beats/min or less than 10 beats/min was not used. A large mixed-effects general linear model was constructed with HR, drug concentrations, and demographic variables as regressors. Specifically, the fixed effects of propofol, fentanyl, and desflurane (plus their linear interaction terms), as well as age, body mass index, ASA Physical Status (II or III), and sex were studied. A random effect of each individual’s mean HR was included. In summary, the model equation was

HR [beats/min] = β0 + β1 × age + β2 × body mass index + β3 × ASA + β4 × (sex = female) + β5 × prop + β6 × fent + β7 × des + β8 × prop × fent + β9 × prop × des + β10 × fent × des + (1 | patient number),

where βi are the model coefficients, and prop = propofol [µg/ml], des = desflurane [%], and fent = fentanyl [ng/ml] effect-site concentrations. To compare possible corticocardiac interactions with the propofol dataset, individual slow waves were extracted from intraoperative EEG (from first incision to the end of closing up) and mean slow-wave frequency per subject extracted and compared with the mean HR.

Statistical Analyses

As these analyses were all post hoc analyses of previously collected and published data, no power calculation was done. Spearman correlations and their P values were used to test associations between electrocardiogram/EEG parameters (HR, R-wave amplitude, root-mean-square successive difference, low-frequency to high-frequency ratio, peak high frequency, slow-wave power) and propofol concentration. Repeated-measures ANOVA was performed on electrocardiogram/EEG-derived parameter traces in 5-min segments to further test for significant changes. For display purposes, mean ± standard error across participants is shown, except where the data were not normally distributed (as tested with D’Agostino and Pearson’s test). In these nonnormally distributed cases, median ± bootstrapped 95% CI (10,000 iterations) are shown. Mean ± SD (or median [25th, 75th percentile]) are given in the text. Significance was set at the P = 0.05 level unless otherwise specified. All custom code used in this study is available at https://gitlab.com/marcoFabus/fabus2022_brain_heart.

Time Series Electrocardiogram analysis

First, we tracked the HR and time series electrocardiogram properties across an ultraslow propofol induction and emergence in n = 16 healthy volunteers (fig. 1). At higher propofol doses, we observed a shortening of the QT segment and decrease in R-wave amplitude (fig. 1A). In every subject, the HR increased and very robustly tracked the propofol dose with Spearman correlation of ρ = 0.923, P < 0.001 (fig. 1B). Mean HR across volunteers increased from 58.2 ± 10 beats/min at baseline to 73.4 ± 8.8 beats/min at peak anesthesia, equivalent to an increase of 4.2 ± 1.5 beats/min/(μg · ml−1). The maximum effect size comparing HR at baseline and peak propofol was Cohen’s d = 1.546. A linear regression showed the HR and propofol relationship to be HR [beats/min] = 56.1 (54.9 to 57.2) + 4.23 (3.75 to 4.80) × propofol [µg/ml], where brackets show 95% CI. Similarly, the R-wave amplitude was also strongly inversely correlated with the propofol effect-site concentration (Spearman ρ = −0.902; P < 0.001; fig. 1C). R-wave amplitude decreased from 966 [707, 1,133] μV at baseline to 742 [627, 1,068] μV at peak anesthesia, equivalent to a decrease of –83 [–245, –28] μV.

Fig. 1.

Propofol increased heart rate and changed electrocardiogram (ECG) shape in healthy volunteers. (A) Mean ECG waveform for all heart beats across n = 16 subjects across propofol doses. High propofol concentration (pink) was characterized by a decrease in R-wave amplitude (top inset) and an earlier T-wave (bottom inset). (B) Propofol increased the heart rate. Group-level heart rate results (purple; mean ± SEM) and propofol effect-site concentration (black; Spearman ρ = 0.923; P < 0.001). (C) At the group level, R-wave amplitude (purple; mean ± SEM) inversely tracked propofol concentration (black; Spearman ρ = −0.902; P < 0.001).

Fig. 1.

Propofol increased heart rate and changed electrocardiogram (ECG) shape in healthy volunteers. (A) Mean ECG waveform for all heart beats across n = 16 subjects across propofol doses. High propofol concentration (pink) was characterized by a decrease in R-wave amplitude (top inset) and an earlier T-wave (bottom inset). (B) Propofol increased the heart rate. Group-level heart rate results (purple; mean ± SEM) and propofol effect-site concentration (black; Spearman ρ = 0.923; P < 0.001). (C) At the group level, R-wave amplitude (purple; mean ± SEM) inversely tracked propofol concentration (black; Spearman ρ = −0.902; P < 0.001).

Close modal

HR Variability Analysis

Next, we studied autonomic activity through HR variability (fig. 2). The root-mean-square successive difference between heartbeats, which indexes parasympathetic tone, decreased in proportion to propofol concentration, and rebounded on emergence (fig. 2A; Spearman ρ = −0.785; P < 0.001; Cohen’s d = 1.296 for baseline vs. peak concentration). This was confirmed by a repeated-measures ANOVA across subjects with significance P < 0.001.

Fig. 2.

Propofol decreased parasympathetic activity, as indexed by heart rate variability in healthy volunteers. (A) Root-mean-square successive difference between heartbeats, a measure of parasympathetic activity (blue; median ± 95% CI), inversely tracks propofol concentration and rebounds on emergence (black). (B) Low-frequency to high-frequency ratio (green; median ± 95% CI), a metric thought to index balance of sympathetic and parasympathetic activity, shows an increase with propofol (black).

Fig. 2.

Propofol decreased parasympathetic activity, as indexed by heart rate variability in healthy volunteers. (A) Root-mean-square successive difference between heartbeats, a measure of parasympathetic activity (blue; median ± 95% CI), inversely tracks propofol concentration and rebounds on emergence (black). (B) Low-frequency to high-frequency ratio (green; median ± 95% CI), a metric thought to index balance of sympathetic and parasympathetic activity, shows an increase with propofol (black).

Close modal

With regard to the frequency domain metrics, the low- to high-frequency ratio showed higher between-subject variability, but the group average confirmed the shift toward a relative predominance of sympathetic activity with increasing propofol concentration (fig. 2B; Spearman ρ = −0.763; P < 0.001). The associated repeated-measures ANOVA result also showed a significant change (P = 0.003; Cohen’s d = 0.539) between baseline and highest propofol concentration. The peak frequency in the high-frequency parasympathetic HR variability range also tracked with propofol concentrations (supplementary fig. 7, https://links.lww.com/ALN/D329; Spearman ρ = 0.885; P < 0.001; repeated-measures ANOVA P < 0.001).

Slow-wave Analysis

Cortical activity during propofol anesthesia is known to be associated with non-REM deep sleep–like slow-wave activity (fig. 3). We first confirmed the previous finding of saturation of frontal slow-wave activity with propofol dose at the Fz electrode8  (fig. 3A). However, more strikingly, this slow-wave activity increase correlated very strongly with the increasing HR (fig. 3B; Spearman ρ = 0.910; P < 0.001).

Fig. 3.

Propofol increased frontal cortical slow-wave activity, and this increase tracked heart rate increases in healthy volunteers. (A) Group-level Fz slow-wave activity results (purple; mean ± standard error) against propofol concentration (black). Slow-wave activity increases and plateaus with drug dose. (B) At the group level, increases in heart rate saw corresponding increases in slow-wave power. Each purple dot represents a 5-min segment of the experiment with standard errors across subjects shown.

Fig. 3.

Propofol increased frontal cortical slow-wave activity, and this increase tracked heart rate increases in healthy volunteers. (A) Group-level Fz slow-wave activity results (purple; mean ± standard error) against propofol concentration (black). Slow-wave activity increases and plateaus with drug dose. (B) At the group level, increases in heart rate saw corresponding increases in slow-wave power. Each purple dot represents a 5-min segment of the experiment with standard errors across subjects shown.

Close modal

Corticocardiac Coupling

The identified correlation between slow-wave activity and HR, as well as previous literature describing their coupling in sleep, led us to focus on quantifying the presence of any time-related coupling between individual slow waves and heartbeats. The methodology is illustrated on the single-subject level in figure 4. After identifying individual slow waves, we studied the distribution of the eight closest heartbeats around each slow-wave onset. If there was no coupling, it would be expected that R-wave to slow-wave intervals should follow a uniform probability distribution, and that a time-averaged electrocardiogram would converge on a horizontal line around zero. As observed in previous work on cardiorespiratory coupling,28  the distribution of time intervals between electrocardiogram R-waves and EEG slow waves (R-wave to slow-wave intervals; fig. 4B) was nonuniform and concentrated around specific phases in the slow-wave cycle. This appeared as a residual low-frequency oscillation in the electrocardiogram, after averaging around the slow-wave onset (fig. 4C), and as peaks in the distribution of heartbeat timings (fig. 4D).

Fig. 4.

Single-subject low-frequency corticocardiac coupling. (A) Example 15 s of electroencephalogram (EEG) and electrocardiogram (ECG) data at high propofol concentrations. Slow-wave starts are shown with blue dots. For each slow-wave start, duration of intervals to nearest heartbeats is determined (R-wave/slow-wave). Individual R-waves are marked with red vertical lines. (B) Example illustrative raster plot of R-wave/slow-wave intervals during 10 min of peak anesthesia. Heartbeats are not distributed randomly and show modest horizontal clustering in horizontal lines, demonstrating nonrandom coupling to slow-wave onset (proportional Shannon entropy = 0.890). (C) Single-subject average slow-wave (blue; mean ± SEM across all slow waves) and ECG amplitudes (red = broadband; purple = 0.5 to 1.5 Hz only) time-locked to slow-wave onset. The ECG pattern is not uniform random noise but shows a clear ongoing low-frequency ECG oscillation. (D) Histogram of R-wave timings relative to slow-wave onset. Individual heartbeats preferentially occur in phase with the slow wave, explaining the oscillatory appearance of C and stripes seen in B.

Fig. 4.

Single-subject low-frequency corticocardiac coupling. (A) Example 15 s of electroencephalogram (EEG) and electrocardiogram (ECG) data at high propofol concentrations. Slow-wave starts are shown with blue dots. For each slow-wave start, duration of intervals to nearest heartbeats is determined (R-wave/slow-wave). Individual R-waves are marked with red vertical lines. (B) Example illustrative raster plot of R-wave/slow-wave intervals during 10 min of peak anesthesia. Heartbeats are not distributed randomly and show modest horizontal clustering in horizontal lines, demonstrating nonrandom coupling to slow-wave onset (proportional Shannon entropy = 0.890). (C) Single-subject average slow-wave (blue; mean ± SEM across all slow waves) and ECG amplitudes (red = broadband; purple = 0.5 to 1.5 Hz only) time-locked to slow-wave onset. The ECG pattern is not uniform random noise but shows a clear ongoing low-frequency ECG oscillation. (D) Histogram of R-wave timings relative to slow-wave onset. Individual heartbeats preferentially occur in phase with the slow wave, explaining the oscillatory appearance of C and stripes seen in B.

Close modal

Importantly, this effect was present and significant at the group level (fig. 5). The group-average lag between the electrocardiogram peak and EEG slow-wave onset was 447 [392, 510] ms (fig. 5B). The slow-wave/R-wave coupling, as measured by entropy in relation to a uniform null distribution, was SHP = 0.866 ± 0.05 (P < 0.001 compared to a uniform null hypothesis). Additional tests to verify this is not a random effect were carried out and can be found in the Supplemental Digital Content (https://links.lww.com/ALN/D329). Furthermore, at the group level, the subjects’ mean HRs and slow-wave frequencies were significantly linearly correlated (Pearson r = 0.519; P = 0.0395).

Fig. 5.

Group-level corticocardiac coupling in healthy volunteers. (A) Overall group-average slow wave (blue) with each subject’s mean in gray. (B) Group-average electrocardiogram (ECG; red = broadband; purple = 0.5 to 4 Hz only) time-locked to slow-wave (black) onset. The ongoing low-frequency ECG oscillation is also present at the group level. (C) Slow-wave onset is significantly linked to the preceding heartbeat. The mean proportional Shannon entropy of intervals between the electroencephalogram (EEG) slow-wave start and previous ECG R-wave across volunteers is shown in red, and the value for each of n = 16 subjects is shown as a gray vertical line. The entropy of the surrogate distribution assuming random timings intervals between the EEG slow-wave start and previous ECG R-wave is in blue. (D) Individual heart rate and slow-wave frequency are related (Pearson r = 0.519; P = 0.0395).

Fig. 5.

Group-level corticocardiac coupling in healthy volunteers. (A) Overall group-average slow wave (blue) with each subject’s mean in gray. (B) Group-average electrocardiogram (ECG; red = broadband; purple = 0.5 to 4 Hz only) time-locked to slow-wave (black) onset. The ongoing low-frequency ECG oscillation is also present at the group level. (C) Slow-wave onset is significantly linked to the preceding heartbeat. The mean proportional Shannon entropy of intervals between the electroencephalogram (EEG) slow-wave start and previous ECG R-wave across volunteers is shown in red, and the value for each of n = 16 subjects is shown as a gray vertical line. The entropy of the surrogate distribution assuming random timings intervals between the EEG slow-wave start and previous ECG R-wave is in blue. (D) Individual heart rate and slow-wave frequency are related (Pearson r = 0.519; P = 0.0395).

Close modal

These analysis results were qualitatively unchanged when EEG data were re-referenced to linked mastoids (supplementary fig. 3, https://links.lww.com/ALN/D329) and when the electrocardiogram was time-locked to slow-wave trough instead of downward zero crossing (supplementary fig. 4, https://links.lww.com/ALN/D329).

Clinical Dataset Analysis

In order to explore whether these HR results hold in a clinical setting, we analyzed the association between effect-site drug concentration and HR using a large general linear model with n = 96 older patients collected as part of the AlphaMax trial (table 2; supplementary fig. 5, https://links.lww.com/ALN/D329). After adjusting for age, body mass index, sex, and ASA Physical Status, all agents (propofol, fentanyl, and desflurane) had a significant effect on the HR (P < 0.001). Propofol led to a mild increase in HR, on average, with a coefficient of +1.3 beats/min/(μg · ml−1) (95% CI, 1.1 to 1.5). Fentanyl, however, led to a decrease in the HR, on average −2.6 (95% CI, − 2.7 to −2.5) beats/min/(ng · ml−1), as did desflurane, with average of −1.84 (95% CI, −1.90 to −1.78) beats/min/(1% end-tidal concentration). The interaction terms were also significant, although with smaller coefficients. With mean individual HR included as a regressor, no demographic parameters were significant, suggesting that the drug effects on the HR may be independent of these demographic variables. The effect size comparing HR with propofol less than 0.5 µg/ml and greater than 3 µg/ml was Cohen’s d = 0.796.

Table 2.

General Linear Model Coefficients for Heart Rate in the AlphaMax Clinical Study

General Linear Model Coefficients for Heart Rate in the AlphaMax Clinical Study
General Linear Model Coefficients for Heart Rate in the AlphaMax Clinical Study

Interestingly, at the group level, mean intraoperative slow-wave frequency was not related to the mean HR in this dataset (appendix 1; P = 0.65). Furthermore, desflurane-fentanyl slow-wave frequency was significantly higher than propofol slow-wave frequency (propofol f = 1.01 ± 0.11 Hz; desflurane-fentanyl f = 1.26 ± 0.15Hz (mean ± SD); Welch t test P < 0.0001), as seen before.31  Unfortunately, individual electrocardiogram waveforms were not included in these data, so coupling between slow waves and individual heartbeats was not assessed.

Propofol and the HR

In this paper, we have shown that administration of propofol leads to increased mean HR. The ultraslow propofol administration in healthy volunteers led to an increase in the mean HR of roughly +4 beats/min/(μg/ml propofol concentration). In an exploratory analysis of an older patient population, the effect of propofol on heart was about threefold smaller, on average +1.3 beats/min/(µg/ml). These clear and significant mean HR increases confirmed our hypothesis but are surprising in view of the mixed existing literature.

Experimental studies have seen a HR increase across a variety of research paradigms.7,9–12,32  We contend the lack of HR increase (or HR decrease) with propofol reported in some clinical studies may be due to other drugs given, the patient population, the surgical context, or dose-/rate-dependent effects. Clinically, it is common to administer opioids and other premedication, which can decrease the HR and affect cardiovascular dynamics.32,33  The smaller HR increase observed in the older clinical population could be due to a previously proposed U-shaped relationship between propofol and HR.11  Older patients have higher anesthetic sensitivity and may be more susceptible to a HR decrease at relatively high propofol concentrations. This is supported by a previous healthy volunteer study where plasma concentrations of about 7.4 µg/ml increased the HR by approximately 30 beats/min, but excessively high concentrations up to mean plasma levels of 18.3 μg/ml reversed the effect and decreased the HR compared to lower concentrations.12  Finally, clinical procedures may provide autonomic stimulation, which could affect intrinsic HR increases with propofol.

The ultraslow induction rate used in our healthy volunteer study may also affect cardiac changes. This is supported by previous work finding rate-dependent cardiac effects of propofol with greater HR decrease in fast inductions, perhaps due to a rate-limiting central nervous system distribution process.19,34  Our HR increase is unlikely to be due to anxiety, as while HR increased from baseline to loss of responsiveness, it continued to increase at the same rate when drug concentration increased beyond the point of loss of consciousness (appendix 2). Finally, while it is plausible some of the HR increase could be due to endothelial irritation, no “pain on propofol injection” phenomenon was reported by volunteers.

Propofol and HR Variability

The biologic basis for the increase in mean HR may be due to propofol inhibiting cardioinhibitory vagal neurons in the brainstem.11  Studies of propofol’s effect on autonomic cardiac influences have also produced mixed results.

The literature agrees that propofol reduces HR variability,6,12,35–37  a result also confirmed in our experiment. The distinct sympathetic and parasympathetic contributions to this are less clear. An early study proposed that propofol mostly depresses sympathetic activity and suggested this as a mechanism for propofol bradycardia and hypotension. However, opioids were also used in that study.36  Several studies since, including this one, have concluded that propofol predominantly decreases high-frequency HR variability. This is conventionally thought to reflect a decrease in parasympathetic vagal influences.6,12,35,37  However, vagal and sympathetic activity tend to be mutually reciprocal, and thus we cannot conclusively show if the vagal decrease is a primary propofol effect or is secondary to sympathetic activation. Notably, however, the low-frequency/high-frequency ratio, traditionally taken as a measure of sympathetic/parasympathetic balance, had a much less consistent relationship with propofol concentration, indicating the sympathetic response was less consistent than the vagal response. We therefore conclude that propofol causes a shift in the autonomic balance to cause tachycardia, probably mainly via the parasympathetic branch, as this relationship was more consistent and more strongly correlated with propofol concentration. This is further supported by our R-wave/slow-wave coupling results. These would not occur through slow sympathetic influences with a lag of up to 10s of seconds.

This increased HR is unlikely to entirely be a reflex tachycardic response to vasodilation, as propofol has been shown to depress the baroreflex.38,39  Specifically, reflex tachycardic responses to hypertension are reduced by propofol in conditions of both normocapnia and hypercapnia.40  In our study, as previously reported,8  hypotension was not observed, and subjects had baseline end-tidal CO2 39.3 ± 3.3 mmHg and peak end-tidal CO2 47.1 ± 5.7 mmHg. A limitation of the current study is the lack of hemodynamic data, as tachycardia associated with decreased stroke volume due to propofol could be part of the HR increase mechanism. Future studies should focus on distinguishing between centrally and hemodynamically mediated mechanisms.41  A final potential confound is differing anesthetic susceptibility, and future studies could benefit from direct sampling of plasma propofol levels.

Propofol and Electrocardiogram Morphology

We observed dose-dependent changes in electrocardiogram morphology that may reflect both direct cardiac effects and centrally mediated changes of propofol. The R-wave amplitude decrease seen in our study might be related to previously observed propofol effects on ventricular depolarisation.42  Propofol may decrease myocardial contractility, possibly due to a direct propofol effect on myocyte ability to expel intracellular calcium,43  although this may only happen above clinical doses.44  A change in the mean electrical axis or direct vagal effect could also explain the R-wave amplitude decrease; findings have been mixed so far.42  Future studies should establish these changes in ventilated patients with more direct cardiac measurements such as cardiac output to assess the clinical significance of this finding.

Corticocardiac Coupling in Propofol

Low-frequency corticocardiac coupling has been observed in sleep.16,45  As propofol slow waves show some sleep-like properties,15  we hypothesized this effect would also be present in anesthesia. We found individual cortical slow waves and cardiac R-waves were coupled as hypothesized. A heartbeat was most likely to precede the slow-wave onset by about 450 ms, a time interval similar to that seen in sleep.16,45  We did not see any evidence of dose-dependent coupling effects (supplementary fig. 10, https://links.lww.com/ALN/D329).

Importantly, this coupling is nontrivial, as it relates to the phase relationship between individual EEG slow waves and the electrocardiogram, not just two ongoing oscillations that happen to have similar frequencies around 1 Hz. The fact that the mean electrocardiogram line in figures 4D and 5C is not zero indicates a time-linked relationship—analogous to evoked potentials in EEG work. We explored this in detail through a simulation study, showing that our proportional entropy metric is sensitive to genuine coupling (Supplemental Digital Content Section 3, https://links.lww.com/ALN/D329).

Mensen et al. proposed several hypotheses for why this coupling may occur.16  The first was a possible metabolic constraint. Overall, neurochemical tone favors hyperpolarized down states with heartbeats acting as a stimulus to evoke a down state when neuronal resources are depleted. Lower regional blood flow between heartbeats could have this effect on a few critical neurons, leading to a network change. However, this seems unlikely as the necessary time resolution of changes in metabolic energy demand seems shorter than that of a damped feeder capillary blood flow, coupled with the energy substrate diffusion time and the presence of intrinsic neuronal energy stores. The other possibility is a third generator controlling both the HR and slow-wave genesis. Knowing this effect is present both in sleep and propofol anesthesia suggests a possible nature of this generator. Sleep and anesthesia differ in noradrenaline levels, but both show low acetylcholine levels.46  Combined with the brainstem projecting both in a cephalad direction to higher brain areas and caudally to the heart, we propose it as a possible place for a common generator. For instance, the nucleus of the solitary tract or cholinergic pontine nuclei may project both to fast-spiking GABA-mediated interneurons in the thalamus and to medullar regions controlling the HR.47,48  Given that the thalamus is involved in slow-wave generation in vivo, a brainstem connection could explain corticocardiac coupling, perhaps by weak-coupling synchronization.49,50  This is supported by subjects with a faster HR also having faster slow-wave frequency. Interestingly, this frequency relationship was not observed during desflurane-fentanyl slow waves (appendix 1), suggesting volatiles may differ in their corticocardiac coupling effects. Further work is needed to explain the relationship between slow waves and cardiac activity, especially as pertains to wider coupling of autonomic and central activity.17,18  Our proposed common brainstem generator could be ruled out if patients with pacemakers also show this coupling.

In summary, in this observational study, slow propofol administration in healthy subjects robustly led to an increase in mean HR that was strongly proportional to drug concentration, and not influenced by changes in behavioral responsiveness. A preliminary analysis in a larger clinical dataset replicated this result, but with a decreased effect size. The HR increase could be explained with decreased parasympathetic inputs, as indexed by decreased high- frequency HR variability. Similar to sleep, frontal cortical slow waves preferentially occurred coupled to the heart rhythm, perhaps due to a common brainstem generator. The observational nature of the study limits causal inferences that can be made, and more work is needed to elucidate the mechanism and role of these cardiac changes and the clinical significance of their coupling to the cortex.

Thus, in the clinical management of patient hemodynamics, propofol should not be assumed to decrease the HR. In fact, particularly for slow infusions and younger patients, propofol is likely to increase the HR. Ultimately, HR will be a complex result of opioid, hypnotic, and surgical factors.

Acknowledgments

The authors thank Irene Tracey, M.A., D.Phil., F.R.C.A., F.Med.Sci., and colleagues at the University of Oxford, Nuffield Department of Clinical Neurosciences, Oxford, United Kingdom, who conducted the original propofol study used in this work.

Research Support

Supported by the Wellcome Trust (London, United Kingdom; grant No. 203139/Z/16/Z; to M.S. Fabus); the Department of Anesthesiology, University of Auckland (Auckland, New Zealand; to Dr. Sleigh); and the Medical Research Council Development Pathway Funding Scheme (London, United Kingdom; award reference No. MR/R006423/1; to Dr. Warnaby).

Competing Interests

Dr. Warnaby is listed as an inventor on patents for “Perception Loss Detection” (World Intellectual Property Organization No. 2013/179048 A1) under anesthesia. Dr. Sleigh is a handling editor for Anesthesiology. Dr. Fabus declares no competing interests.

Supplementary figures and control analyses, https://links.lww.com/ALN/D329

1.
Tramèr
MR
,
Moore
RA
,
McQuay
HJ
:
Propofol and bradycardia: Causation, frequency and severity.
Br J Anaesth
1997
;
78
:
642
51
2.
Wouters
PF
,
Van de Velde
M
,
Van Hemelrijck
J
:
Cardiovascular effects of propofol.
Semin Cardiothorac Vasc Anesth
1997
;
1
:
121
31
3.
Hug
CC
,
McLeskey
CH
,
Nahrwold
ML
,
Roizen
MF
,
Stanley
TH
,
Thisted
RA
,
Walawander
CA
,
White
PF
,
Apfelbaum
JL
,
Grasela
TH
:
Hemodynamic effects of propofol: Data from over 25,000 patients.
Anesth Analg
1993
;
77
:
S21
9
4.
McKeage
K
,
Perry
CM
:
Propofol.
CNS Drugs
2003
;
17
:
235
72
5.
Fulton
B
,
Sorkin
EM
:
Propofol.
Drugs
1995
;
50
:
636
57
6.
Kanaya
N
,
Hirata
N
,
Kurosawa
S
,
Nakayama
M
,
Namiki
A
:
Differential effects of propofol and sevoflurane on heart rate variability.
Anesthesiology
2003
;
98
:
34
40
7.
de Wit
F
,
van Vliet
AL
,
de Wilde
RB
,
Jansen
JR
,
Vuyk
J
,
Aarts
LP
,
de Jonge
E
,
Veelo
DP
,
Geerts
BF
:
The effect of propofol on haemodynamics: Cardiac output, venous return, mean systemic filling pressure, and vascular resistances.
Br J Anaesth
2016
;
116
:
784
9
8.
Ni Mhuircheartaigh
R
,
Warnaby
C
,
Rogers
R
,
Jbabdi
S
,
Tracey
I
:
Slow-wave activity saturation and thalamocortical isolation during propofol anesthesia in humans.
Sci Transl Med
2013
;
5
:
208ra148
9.
Claeys
MA
,
Gepts
E
,
Camu
F
:
Haemodynamic changes during anaesthesia induced and maintained with propofol.
Br J Anaesth
1988
;
60
:
3
9
10.
Blake
DW
,
Jover
B
,
McGrath
BP
:
Haemodynamic and heart rate reflex responses to propofol in the rabbit. Comparison with althesin.
Br J Anaesth
1988
;
61
:
194
9
11.
Wang
X
,
Huang
Z-G
,
Gold
A
,
Bouairi
E
,
Evans
C
,
Andresen
MC
,
Mendelowitz
D
:
Propofol modulates γ-aminobutyric acid–mediated inhibitory neurotransmission to cardiac vagal neurons in the nucleus ambiguus.
Anesthesiology
2004
;
100
:
1198
205
12.
Mäenpää
M
,
Penttilä
J
,
Laitio
T
,
Kaisti
K
,
Kuusela
T
,
Hinkka
S
,
Scheinin
H
:
The effects of surgical levels of sevoflurane and propofol anaesthesia on heart rate variability.
Eur J Anaesthesiol
2007
;
24
:
626
33
13.
Steriade
M
,
Nuñez
A
,
Amzica
F
:
A novel slow (< 1 Hz) oscillation of neocortical neurons in vivo: Depolarizing and hyperpolarizing components.
J Neurosci
1993
;
13
:
3252
65
14.
Rosanova
M
,
Fecchio
M
,
Casarotto
S
,
Sarasso
S
,
Casali
AG
,
Pigorini
A
,
Comanducci
A
,
Seregni
F
,
Devalle
G
,
Citerio
G
,
Bodart
O
,
Boly
M
,
Gosseries
O
,
Laureys
S
,
Massimini
M
:
Sleep-like cortical OFF-periods disrupt causality and complexity in the brain of unresponsive wakefulness syndrome patients.
Nat Commun
2018
;
9
:
4427
15.
Murphy
M
,
Bruno
M-A
,
Riedner
BA
,
Boveroux
P
,
Noirhomme
Q
,
Landsness
EC
,
Brichant
J-F
,
Phillips
C
,
Massimini
M
,
Laureys
S
,
Tononi
G
,
Boly
M
:
Propofol anesthesia and sleep: A high-density EEG study.
Sleep
2011
;
34
:
283
91A
16.
Mensen
A
,
Zhang
Z
,
Qi
M
,
Khatami
R
:
The occurrence of individual slow waves in sleep is predicted by heart rate.
Sci Rep
2016
;
6
:
29671
17.
Zambotti
M de
,
Trinder
J
,
Silvani
A
,
Colrain
IM
,
Baker
FC
:
Dynamic coupling between the central and autonomic nervous systems during sleep: A review.
Neurosci Biobehav Rev
2018
;
90
:
84
103
18.
Azzalini
D
,
Rebollo
I
,
Tallon-Baudry
C
:
Visceral signals shape brain dynamics and cognition.
Trends Cogn Sci
2019
;
23
:
488
509
19.
Obert
DP
,
Sepúlveda
P
,
Kratzer
S
,
Schneider
G
,
Kreuzer
M
:
The influence of induction speed on the frontal (processed) EEG.
Sci Rep
2020
;
10
:
19444
20.
Gaskell
A
,
Pullon
R
,
Hight
D
,
Termaat
J
,
Mans
G
,
Voss
L
,
Kreuzer
M
,
Schmid
S
,
Kratzer
S
,
Rodriguez
A
,
Schneider
G
,
Garcia
P
,
Sleigh
J
:
Modulation of frontal EEG alpha oscillations during maintenance and emergence phases of general anaesthesia to improve early neurocognitive recovery in older patients: Protocol for a randomised controlled trial.
Trials
2019
;
20
:
146
21.
Nunez
PL
:
REST: A good idea but not the gold standard.
Clin Neurophysiol
2010
;
121
:
2177
80
22.
Hamilton
P
:
Open source ECG analysis.
Comput Cardiol
2002
;
29
:
101
4
23.
Gomes
P
,
Margaritoff
P
,
Silva
H
:
pyHRV: Development and evaluation of an open-source python toolbox for heart rate variability (HRV)
,
Proceedings on the International Conference on Electrical, Electronic and Computing Engineering (IcETRAN)
.
Belgrade
,
ETRAN Society
,
2019
, pp
822
8
24.
Shaffer
F
,
Ginsberg
JP
:
An overview of heart rate variability metrics and norms.
Front Public Health
2017
;
5
:
258
25.
Vallat
R
,
Walker
MP
:
An open-source, high-performance tool for automated sleep staging.
Elife
2021
;
10
:
e70092
26.
Amzica
F
,
Steriade
M
:
Electrophysiological correlates of sleep delta waves.
Electroencephalogr Clin Neurophysiol
1998
;
107
:
69
83
27.
Siclari
F
,
Bernardi
G
,
Riedner
BA
,
LaRocque
JJ
,
Benca
RM
,
Tononi
G
:
Two distinct synchronization processes in the transition to sleep: A high-density electroencephalographic study.
Sleep
2014
;
37
:
1621
37
28.
Galletly
DC
,
Larsen
PD
:
Cardioventilatory coupling during anaesthesia.
Br J Anaesth
1997
;
79
:
35
40
29.
Galletly
DC
,
Larsen
PD
:
Relationship between cardioventilatory coupling and respiratory sinus arrhythmia.
Br J Anaesth
1998
;
80
:
164
8
30.
Tzeng
YC
,
Larsen
PD
,
Galletly
DC
:
Cardioventilatory coupling in resting human subjects.
Exp Physiol
2003
;
88
:
775
82
31.
Leslie
K
,
Sleigh
J
,
Paech
MJ
,
Voss
L
,
Lim
CW
,
Sleigh
C
:
Dreaming and electroencephalographic changes during anesthesia maintained with propofol or desflurane.
Anesthesiology
2009
;
111
:
547
55
32.
Stephan
H
,
Sonntag
H
,
Schenk
HD
,
Kettler
D
,
Khambatta
HJ
:
Effects of propofol on cardiovascular dynamics, myocardial blood flow and myocardial metabolism in patients with coronary artery disease.
Br J Anaesth
1986
;
58
:
969
75
33.
Grounds
RM
,
Twigley
A
,
Carli
F
,
Whitwam
JG
,
Morgan
M
:
The haemodynamic effects of intravenous induction.
Anaesthesia
1985
;
40
:
735
40
34.
Stokes
DN
,
Hutton
P
:
Rate-dependent induction phenomena with propofol: Implications for the relative potency of intravenous anesthetics.
Anesth Analg
1991
;
72
:
578
83
35.
Galletly
DC
,
Buckley
DHF
,
Robinson
BJ
,
Corfiatis
T
:
Heart rate variability during propofol anaesthesia.
Br J Anaesth
1994
;
72
:
219
20
36.
Deutschman
CS
,
Harris
AP
,
Fleisher
LA
:
Changes in heart rate variability under propofol anesthesia: A possible explanation for propofol-induced bradycardia.
Anesth Analg
1994
;
79
:
373
7
37.
Scheffer
GJ
,
Ten Voorde
BJ
,
Karemaker
JM
,
Ros
HH
,
Lange
JJ de
:
Effects of thiopentone, etomidate and propofol on beat-to-beat cardiovascular signals in man.
Anaesthesia
1993
;
48
:
849
55
38.
Sato
M
,
Tanaka
M
,
Umehara
S
,
Nishikawa
T
:
Baroreflex control of heart rate during and after propofol infusion in humans.
Br J Anaesth
2005
;
94
:
577
81
39.
Cullen
PM
,
Turtle
M
,
Prys-Roberts
C
,
Way
WL
,
Dye
J
:
Effect of propofol anesthesia on baroreflex activity in humans.
Anesth Analg
1987
;
66
:
1115
20
40.
Ebert
TJ
,
Muzi
M
:
Propofol and autonomic reflex function in humans.
Anesth Analg
1994
;
78
:
369
75
41.
Paventi
S
,
Santevecchi
A
,
Ranieri
R
:
Effects of sevoflurane versus propofol on QT interval.
Minerva Anestesiol
2001
;
67
:
637
40
42.
Tascanov
MB
,
Tanriverdi
Z
,
Gungoren
F
,
Besli
F
,
Bicer Yesilay
A
,
Altiparmak
HI
,
Bayram
G
,
Demir
K
:
The effect of propofol on frontal QRS-T angle in patients undergoing elective colonoscopy procedure.
J Clin Pharm Ther
2020
;
45
:
185
90
43.
Puttick
RM
,
Diedericks
J
,
Sear
JW
,
Glen
JB
,
Foëx
P
,
Ryder
WA
:
Effect of graded infusion rates of propofol on regional and global left ventricular function in the dog.
Br J Anaesth
1992
;
69
:
375
81
44.
Sprung
J
,
Ogletree-Hughes
ML
,
McConnell
BK
,
Zakhary
DR
,
Smolsky
SM
,
Moravec
CS
:
The effects of propofol on the contractility of failing and nonfailing human heart muscles.
Anesth Analg
2001
;
93
:
550
9
45.
Lechinger
J
,
Heib
DPJ
,
Gruber
W
,
Schabus
M
,
Klimesch
W
:
Heartbeat-related EEG amplitude and phase modulations from wakefulness to deep sleep: Interactions with sleep spindles and slow oscillations.
Psychophysiology
2015
;
52
:
1441
50
46.
Pal
D
,
Mashour
GA
:
Consciousness, anesthesia, and acetylcholine.
Anesthesiology
2021
;
134
:
515
7
47.
Benarroch
EE
:
Interface between the autonomic and arousal systems.
Auton Neurosci
2015
;
192
:
38
48.
Naji
M
,
Krishnan
GP
,
McDevitt
EA
,
Bazhenov
M
,
Mednick
SC
:
Timing between cortical slow oscillations and heart rate bursts during sleep predicts temporal processing speed, but not offline consolidation.
J Cogn Neurosci
2019
;
31
:
1484
90
49.
Sanchez-Vives
MV
:
Origin and dynamics of cortical slow oscillations.
Curr Opin Physiol
2020
;
15
:
217
23
50.
Glass
L
:
Synchronization and rhythmic processes in physiology.
Nature
2001
;
410
:
277
84

Appendix 1: Slow-wave Frequency in Clinical Data

Fig. A1.

Mean slow-wave (SW) frequency and heart rate (beats/min) are not related during clinical desflurane- fentanyl intraoperative general anesthesia (n = 96 patients; Pearson r = –0.05; P = 0.65).

Fig. A1.

Mean slow-wave (SW) frequency and heart rate (beats/min) are not related during clinical desflurane- fentanyl intraoperative general anesthesia (n = 96 patients; Pearson r = –0.05; P = 0.65).

Close modal

Appendix 2: Heart Rate around Loss of Responsiveness in Volunteer Data

Fig. A2.

Heart rate (beats/min) in ultra-slow propofol infusion increases from baseline to loss of behavioral response (LOBR) to auditory and laser pain stimuli and peak anesthesia. **Bonferroni-corrected paired t test P < 0.01; ***P < 0.001.

Fig. A2.

Heart rate (beats/min) in ultra-slow propofol infusion increases from baseline to loss of behavioral response (LOBR) to auditory and laser pain stimuli and peak anesthesia. **Bonferroni-corrected paired t test P < 0.01; ***P < 0.001.

Close modal
This is an open access article distributed under the Creative Commons Attribution License 4.0 (CCBY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.