Dynamic Substates of the Default Mode Network Are Associated with Mystical Experiences and Clinical Outcomes after Magnesium-Ibogaine Treatment in Veterans with Traumatic Brain Injury
This open-label observational study (n=30) in US veterans with traumatic brain injury (TBI) found that magnesium-ibogaine did not change overall default mode network connectivity, but it did alter a dynamic brain state linked to mystical experiences and clinical improvement one month later. The affected brain pattern was related to a shift towards more internally focused thought.
1 linked clinical trial·23 references indexed in Blossom
Authors
- Shinozuka, K.
- Olash, C.
- Han, T.
Published
Abstract
Preliminary evidence suggests that the combination of magnesium and ibogaine, an atypical psychedelic, may be a promising treatment for post-traumatic stress disorder (PTSD), opioid use disorder, and traumatic brain injury (TBI). The "mystical" experience elicited by ibogaine, which is characterized by feelings of awe, selflessness, and transcendence, is correlated with improvements in PTSD symptoms. Mystical experiences with other psychedelics are associated with acute decreases in the activity and connectivity of the default mode network (DMN), which mediates self-related cognition. However, the dynamic effects of ibogaine on the DMN have not yet been studied. At baseline, immediately (3-4 days) after ibogaine, and one month after ibogaine, we acquired resting-state functional magnetic resonance imaging data in an open-label, observational trial of magnesium-ibogaine treatment for 30 U.S. veterans with TBI. Magnesium-ibogaine did not significantly alter static DMN connectivity at either the immediate-post or one-month timepoint. Since static measures cannot capture time-evolving changes in connectivity, we next used Hidden Markov Models (HMM) to measure the post-acute dynamics of DMN activity. Magnesium-ibogaine was associated with significant, sustained decreases in the switching rate (i.e., increases in the duration of) a dynamic DMN substate, which was significantly correlated with the mean score on the Revised Mystical Experience Questionnaire and clinical improvements at one month-post treatment. This DMN substate exhibited a lateral-medial spatial gradient, which was significantly associated with a gradient of externally oriented (i.e., directed to the environment) to internally oriented (i.e., self-related) perception and cognition. Taken together, our results indicate that magnesium-ibogaine alters specific dynamic substates of the DMN, which correlate with its subjective and therapeutic effects.
Research Summary of 'Dynamic Substates of the Default Mode Network Are Associated with Mystical Experiences and Clinical Outcomes after Magnesium-Ibogaine Treatment in Veterans with Traumatic Brain Injury'
βBlossom's Take
Introduction
Psychedelics have shown promise for a range of psychiatric and neurological conditions, and earlier neuroimaging research has linked their effects to changes in the default mode network (DMN), a network involved in self-related thought, autobiographical memory, and mind-wandering. However, previous findings have been mixed, and most studies have relied on static connectivity measures that average brain activity across an entire scan. The authors also note that ibogaine, an atypical psychedelic, has only recently begun to show preliminary therapeutic potential, including in veterans with traumatic brain injury (TBI), and that its neural effects on the DMN had not yet been studied dynamically. Shinozuka and colleagues set out to test whether magnesium-ibogaine treatment alters DMN connectivity in veterans with TBI, and whether any neural changes relate to mystical experience and later clinical outcomes. They aimed to compare static and dynamic measures of DMN function across baseline, immediately after treatment, and one month later, with the broader goal of identifying more precise neural correlates of ibogaine’s subjective and therapeutic effects. The study builds on the same MISTIC cohort by focusing on temporal patterns within the DMN rather than treating the network as uniform across time.
Methods
The researchers used data from the prior MISTIC open-label observational study. The sample comprised 30 male U.S. veterans with mild-to-moderate TBI, each of whom received a single dose of ibogaine (12.1 ± 1.2 mg/kg) together with intravenous magnesium sulphate given before and after ibogaine (1 g each time) to reduce cardiovascular risk, particularly QT interval prolongation. Resting-state fMRI was acquired at three main timepoints: baseline, immediately post-treatment (3-4 days after ibogaine), and one month post-treatment. The extracted text reports usable fMRI data for 29 participants at baseline, 28 at the immediate post-treatment timepoint, and 25 at one month, with some missing data because of MRI safety issues or other logistical reasons. The researchers also used clinical and subjective measures collected in the parent study, including the Revised Mystical Experience Questionnaire (MEQ-30), CAPS-5 for PTSD, MADRS for depression, HAM-A for anxiety, and WHODAS-2.0 for disability. Long-term clinical follow-up was also mentioned at 3, 6, 9, and 12 months, and baseline TBI burden was characterised using the Boston Assessment of TBI-Lifetime. For static functional connectivity, the brain was parcellated with the Schaefer-200 atlas and grouped into the seven Yeo networks. The researchers calculated Fisher-transformed Pearson correlations between regions, averaged them within and between networks, and tested timepoint effects with linear mixed-effects models with subject as a random effect. They corrected timepoint contrasts using Benjamini-Hochberg false discovery rate (FDR) correction. They also correlated MEQ-30 scores with changes in within-DMN connectivity, and explored associations between clinical outcomes and static DMN connectivity. For dynamic analyses, the researchers restricted the HMM analysis to the 37 DMN parcels, fitting Gaussian Hidden Markov Models across K = 2 to K = 10 states and selecting the number of states using the elbow of the free-energy curve; K = 5 was chosen. They quantified fractional occupancy and switching rate, where switching rate refers to how often the brain transitions between substates and therefore inversely reflects how long a state persists. Mixed-effects models were again used to test timepoint effects, with additional models adjusting for head motion, age, and number of TBIs. They also correlated session-normalised changes in DMN substate 3 switching rate with MEQ-30 scores and with clinical outcomes, applying FDR correction. In a further exploratory analysis, the researchers fit HMMs to the whole-brain time series after principal component reduction to 15 dimensions. To interpret one DMN substate, the researchers constructed an internal-to-external cognition gradient using Cognitive Atlas terms, Neurosynth-derived activation maps, and a semantic embedding approach. They then tested whether the spatial pattern of DMN substate activity aligned with this gradient, using a null model that preserved spatial autocorrelation.
Results
Static connectivity analysis did not show a significant mean change in within-DMN functional connectivity after ibogaine at either immediate post-treatment or one month. The reported estimates were β = -0.016, pFDR = 0.438 at the immediate timepoint and β = -0.028, pFDR = 0.175 at one month. The only significant static within-network change was a reduction in salience network connectivity immediately after treatment (β = -0.064, pFDR = 0.016), which was not present at one month. There were no significant between-network connectivity changes. Although mean within-DMN connectivity did not shift significantly, individual differences mattered. MEQ-30 mystical experience scores were positively correlated with changes in within-DMN connectivity at one month (r = 0.478, pFDR = 0.036), with a weaker trend immediately after treatment (r = 0.336, pFDR = 0.087). However, static within-DMN connectivity was not significantly associated with any of the clinical outcomes reported in the Methods at either post-treatment timepoint. The dynamic HMM analysis identified 5 DMN substates as optimal. Ibogaine was not associated with a significant change in fractional occupancy of any substate, but it was associated with reduced switching rate, meaning longer duration, of DMN substate 2 immediately after treatment (β = -0.010, pFDR = 0.020) and DMN substate 3 at both immediate post-treatment (β = -0.006, pFDR = 0.038) and one month (β = -0.007, pFDR = 0.020). These findings remained significant after adjusting for head motion, age, and number of TBIs. Head motion itself was not correlated with switching rates. For DMN substate 3, the session-normalised switching rate at the immediate post-treatment timepoint was negatively correlated with MEQ-30 score (r = -0.457, pFDR = 0.033), indicating that stronger mystical experiences were associated with a greater reduction in switching rate, or longer state duration. The one-month correlation remained negative but was not significant (r = -0.271, pFDR = 0.200). This substate showed higher activation in lateral than medial DMN regions, and this pattern aligned significantly with a gradient from externally oriented to internally oriented cognition (r = 0.346, pspatial = 0.045). The significance of this gradient was reported to persist at a 15% selection quantile but not at 5% or 20%. Exploratory analyses found that the change in session-normalised switching rate of DMN substate 3 at one month was positively correlated with all four clinical outcomes at that timepoint: CAPS-5 (r = 0.460, pFDR = 0.048), MADRS (r = 0.567, pFDR = 0.015), HAM-A (r = 0.700, pFDR = 0.001), and WHODAS-2.0 (r = 0.548, pFDR = 0.018). The authors interpret these positive correlations as meaning that larger increases in the duration of DMN substate 3 were associated with greater clinical improvement. The same switching-rate change also correlated with later improvements at 12 months in HAM-A (r = 0.627, pFDR = 0.015) and WHODAS (r = 0.551, pFDR = 0.035). Baseline TBI burden was not significantly correlated with changes in this switching rate. Whole-brain HMMs identified states that resembled canonical networks, including the DMN, but none of these whole-brain states showed significant switching-rate changes after ibogaine, and the whole-brain DMN-like state was not significantly correlated with MEQ-30 or clinical outcomes. Overall, the dynamic DMN substate 3 was the main neural feature linked to both mystical experience and clinical improvement.
Discussion
The authors argue that this is the first study to describe post-acute effects of ibogaine on DMN dynamics, and they present three main interpretations. First, dynamic measures were more sensitive than static measures for detecting ibogaine-related changes and for relating brain changes to clinically meaningful outcomes. In their view, static connectivity still captured some individual variability, especially its positive association with mystical experience at one month, but it did not identify the treatment-related temporal reorganisation that dynamic analysis revealed. Second, they interpret the reduction in switching rate of DMN substate 3 as a specific post-treatment alteration that tracked both mystical experience and clinical improvement. They contrast this with earlier psychedelic studies, which often measured the brain during the acute drug experience and reported decreases in within-DMN connectivity. The authors note that their fMRI scans were collected 3-4 days and one month after treatment, so the post-acute pattern they observed may differ from acute psychedelic effects. They also note that the salience network showed a transient reduction in static connectivity immediately after treatment, which they suggest could relate to reduced hypervigilance in PTSD, but this did not persist at one month. Third, they argue that DMN substate 3 is spatially meaningful because it showed relatively higher activity in lateral DMN regions and lower activity in medial regions, and this pattern aligned with a gradient from externally to internally oriented cognition. They suggest, cautiously, that this may fit the phenomenology of mystical experience, in which self-focused awareness diminishes and engagement with the surrounding environment increases. However, they state that the functional meaning of this gradient was not directly experimentally verified in the study and requires further validation. The authors acknowledge several limitations. The treatment was open-label, so expectancy and context effects cannot be separated from ibogaine itself, especially because the treatment setting included other wellness activities such as sweat lodge, massage, yoga, reiki, breathwork, and meditation. The sample was all male and almost entirely White, limiting generalisability to women, non-White populations, non-veteran patients, and other veteran groups. They also did not test alternative methods for estimating DMN dynamics, such as leading eigenvector dynamics analysis, co-activation patterns, or sliding-window connectivity. Finally, they note that the gradient result was sensitive to the choice of selection quantile, and that future studies should collect fMRI both during and after the ibogaine experience to compare acute and post-acute neural effects directly. In terms of implications, the authors suggest that psychedelic neuroimaging should not treat the DMN as a single temporally uniform system. They propose that specific dynamic substates may be more informative than static averages for understanding both subjective experiences and therapeutic change after ibogaine.
Conclusion
The authors conclude that, in veterans with TBI, dynamic analysis revealed reorganisation of DMN activity after ibogaine that was not captured by static connectivity measures. They state that ibogaine specifically lengthened the duration of one DMN substate, and that this change was associated with mystical experience and improvements in PTSD, depression, anxiety, and functional disability. They also conclude that psychedelic neuroimaging is likely to be more informative when it examines specific DMN substates rather than treating the DMN as temporally uniform.
Expert Research Summary
Go Pro to unlock section-by-section summaries, Blossom's editorial take, and the complete research toolkit.
View full paper sections
SECTION 2.1. PARTICIPANTS AND TREATMENTS
The data were acquired as part of the previous MISTIC study.Details are summarized in Supplementary Methods 1-2. In brief, 30 male veterans with mild-to-moderate TBI received a single dose of ibogaine (12.1 ± 1.2 mg/kg) and were pre-and post-treated with intravenous magnesium sulfate (1 g) to mitigate cardiovascular risks, particularly QT interval prolongation.
SECTION 2.2. SUBJECTIVE AND CLINICAL MEASURES
Mystical experiences on ibogaine were measured retrospectively at the IP timepoint with the Revised Mystical Experience Questionnaire (MEQ-30).Responses to the questions were averaged into one aggregate mean score. Additionally, various clinical measures, including the Clinician-Administered PTSD Scale for DSM-5 (CAPS-5), Montgomery-Asberg Depression Rating Scale (MADRS), Structured Interview Guide for the Hamilton Anxiety Scale (HAM-A), and World Health Organization Disability Assessment Schedule 2.0 (WHODAS-2.0, a measure of functional disability), were administered at baseline, IP, and 1M and four long-term timepoints, which included 3 monthspost, 6 months-post, 9 months-post, and 12 months-post.Number of TBIs and blast exposure, assessed with the Boston Assessment of TBI -Lifetime (BATL), were measured at baseline. Details on each of these measures are provided in the previous MISTIC study, as well as a separate study reporting the long-term follow-up data.Section 2.3. fMRI data acquisition and preprocessing fMRI data were acquired at three timepoints: baseline (n = 29), IP (n = 28), and 1M (n = 25). Data were missing from some participants because of failure to meet MRI safety criteria or other logistical issues. fMRI data acquisition and preprocessing were previously reportedand are summarized in Supplementary Methods 3.
SECTION 2.4. DATA ANALYSIS
A list of primary vs. exploratory analyses is provided in Table. Section 2.4.1. Static FC fMRI data were parcellated with the Schaefer-200 atlas, which is itself divided into the seven Yeo networks: visual network (VN), somatomotor network (SMN), dorsal attention network (DAN), salience network (SN), limbic network (LN), central executive network (CEN), and default mode network (DMN).We then measured the Fisher-transformed Pearson correlation of the FC between all pairs of regions. For static within-network FC, these correlations were averaged across pairs of regions that belonged to the same network. For static between-network FC, correlations were averaged across pairs of regions that belonged to different networks. A linear mixed-effects model (LME) was used to determine statistical significance, with timepoint as a fixed effect and subject as a random effect (FC ~ Timepoint + (1 | Subject)). p-values for timepoint contrasts were corrected for multiple comparisons across timepoints and networks using the Benjamini-Hochberg false discovery rate (FDR) procedure. Pearson correlations were conducted between MEQ-30 ratings and changes in static within-DMN FC at the IP and 1M timepoints. (The data for all correlations in this study generally met the assumptions required for Pearson correlations, including linearity, lack of extreme outliers, independence of observations, and lack of severe skew.) FDR was used to correct for multiple comparisons across the two timepoints. Exploratory Pearson correlations were also conducted between clinical measurestotal scores on the CAPS-5, MADRS, HAM-A, and WHODAS-2.0 scales at IP and 1M -and static within-DMN FC. Due to the exploratory nature of these correlations, only uncorrected p-values are reported.
SECTION 2.4.2. DYNAMIC FC ANALYSIS
Hidden Markov Models (HMMs) are probabilistic models that represent a timeseries as a sequence of discrete, recurring latent states, each associated with a characteristic pattern of observed activity.The transitions between states are governed by a Markov process in which the state at a given timepoint depends only on the state at the preceding timepoint. In applications to fMRI and other human brain data, the observed regional timeseries are modeled as 'emissions' from these hidden brain states.Different types of HMMs use different distributions to describe the probability of the observed timeseries given the hidden state, with the Gaussian distribution being the most common.This framework allows the data to be decomposed into temporally dynamic brain states that can recur across time and across participants, while also estimating the probability of transitioning from one state to another. In contrast to static FC analyses, HMMs explicitly capture moment-to-moment changes in large-scale brain dynamics. The Schaefer-200 atlas was restricted to the 37 parcels assigned to the DMN. Using the GLHMM toolbox in Python,Gaussian HMMs with full covariance ranging from K = 2 to K = 10 states were fit to the concatenated DMN parcel timeseries across baseline, IP, and 1M sessions, while preserving session boundaries. Prior to fitting, parcel timeseries were z-scored within sessions. K was selected using the elbow of the free-energy curve. For each DMN substate, we quantified fractional occupancy and switching rate using GLHMM outputs. Fractional occupancy was defined as the proportion of time spent in a DMN substate within a session, and switching rate as the rate of transitions between DMN substates within a session. Timepoint effects of summary metrics were tested using LMEs with timepoint as a fixed effect and subject as a random intercept (Summary Metric ~ Timepoint + (1 | Subject)); post-baseline contrasts were corrected for multiple comparisons using FDR across DMN substate-by-timepoint tests. To ensure that head motion (measured as mean framewise displacement [FD]), age, and number of TBIs did not account for the results, the LME was rerun with these three covariates added as fixed effects (Summary Metric ~ Timepoint + FD + Age + #TBIs + (1 | Subject)). Pearson correlations were also performed between summary metrics and mean FD at each timepoint. Pearson correlations were conducted between MEQ-30 ratings and session-normalized changes in switching rate of DMN substate 3 at the IP and 1M timepoints. Note that we focused on the switching rate of DMN substate 3 because it was the only summary metric that exhibited significant changes across both timepoints (see Section 3.2). For session normalization, the switching rate of DMN substate 3 per subject per session was divided by the mean of the switching rate across DMN substates for that subject and session. FDR was used to correct for multiple comparisons across the two timepoints. Exploratory Pearson correlations were also conducted between clinical measurestotal scores on the CAPS-5, MADRS, HAM-A, and WHODAS-2.0 scales -and switching rate of DMN substate 3. To reduce the number of multiple comparisons, clinical data from only one longterm follow-up visit, i.e., 12 months-post (12M), were correlated with summary metrics. In total, FDR was used to correct for 12 comparisons (3 timepoints [IP, 1M, 12M] x 4 clinical outcomes). Whole-brain HMMs were also fit to the full Schaefer-200 parcel time series. For these analyses, parcel timeseries were mean-centered across regions at each timepoint to reduce global variance, then standardized within-session and reduced to 15 dimensions using principal component analysis. Gaussian HMMs with full covariance shared across whole-brain states, ranging from K = 2 to K = 10 whole-brain states, were trained. K was once again selected using the elbow of the free-energy curve. Fractional occupancy and switching rate were then compared across timepoints using the same mixedeffects framework described above.
SECTION 2.4.3. GRADIENT OF INTERNAL TO EXTERNAL PERCEPTION
Mystical experiences are characterized by a feeling of oneness with the environment, or a dissolution of the boundaries between the internal and external worlds. We therefore determined whether DMN substates exhibited increased activation in DMN regions involved in externally oriented cognition and decreased activation in DMN regions involved in internally oriented cognition. To do this, we embedded terms in the Cognitive Atlas, a set of words and phrases describing various cognitive phenomena (e.g., "action perception," "motor planning," "agency"), along a semantic axis representing external cognition on one end and internal cognition on the other.Frequency (TF-IDF) matrix, which assigns weights to words based on how strongly they are related to a Cognitive Atlas term.This matrix is then reduced to 100 dimensions via singular value decomposition (SVD). We averaged the vector embeddings of the anchor terms in this 100dimensional space, yielding centroids corresponding to the two poles of the internal-external semantic axis. Then, we computed the cosine similarity of the embedding of each Cognitive Atlas term to the 'internal' and 'external' centroids, thereby scoring each term along the internal-external semantic axis. Positive scores indicated that a term is more related to externally oriented cognition, whereas negative scores indicated that a term is more related to internally oriented cognition. To relate the internal-external semantic axis to DMN regions, we leveraged prior applications of Neurosynth to identify patterns of brain activity encoding each of the Cognitive Atlas terms. Neurosynth is a tool that automatically performs meta-analyses on fMRI papers about various topics, generating an aggregate brain activation map for each topic.Previous studies have applied Neurosynth to Cognitive Atlas terms, resulting in a collection of brain activation maps for each term (identifiers.org/neurovault.collection:1274).Each of these maps was parcellated to the 37 DMN regions in the Schaefer-200 atlas. Maps corresponding to terms with the top 10% score on the internalexternal semantic axis (i.e., the most 'external' terms) were averaged, and the same was done for the terms with the bottom 10% score. (Note that our result remained significant when the selection quantile was set to 15%, but not at 5% or 20%.) The difference between these two averages was then computed for each DMN region, resulting in an 'internal-external gradient' across DMN areas. Finally, we performed Pearson correlations between the DMN substate means and this internalexternal gradient. Statistical significance of the correlation was assessed using a null model that preserves spatial autocorrelation, as implemented in neuromaps.A custom volumetric parcellation containing the 37 DMN parcels was resampled to MNI152 2mm space, and Moran spectral randomization was used to generate surrogate maps with matched spatial autocorrelation structure.
SECTION 3.1. STATIC FC ANALYSIS
A schematic of the study design is shown in Figure. Ibogaine did not significantly alter static within-DMN FC at either IP (β = -0.016, SE = 0.012, pFDR = 0.438) or 1M (β = -0.028, SE = 0.012, pFDR = 0.175) (Figure). The only significant change in static within-network FC occurred in the salience network at IP (β = -0.064, SE = 0.020, pFDR = 0.016), but not at 1M. There were also no significant differences in between-network connectivity (Supplementary Table). The intensity of the mystical experience, as quantified by the mean score on the MEQ-30, was significantly and positively correlated with changes in within-DMN FC at 1M (r = 0.478, pFDR = 0.036) (Figure). The correlation trended towards significance at the IP timepoint (r = 0.336, pFDR = 0.087). Thus, although within-DMN FC did not exhibit a significant mean shift across the cohort, inter-individual differences in within-DMN FC were positively associated with MEQ. However, there were no significant correlations between any of the clinical outcomes (see Methods section) at either timepoint and within-DMN FC (Supplementary Figure), or FC within the salience network.
SECTION 3.2. DYNAMIC FC ANALYSIS
Since static FC analysis averages connectivity over the entire scanning session, we next determined whether a dynamical analysis, namely HMMs, could better differentiate DMN activity before and after ibogaine. K = 5 was identified as the optimal number of DMN substates based on the free energy of the model (Supplementary Figure). Each DMN substate was characterized by a different pattern of mean activation across DMN regions (Figure) and covariance between DMN regions (Supplementary Figure). While ibogaine was not associated with a significant change in the fractional occupancy of any of the DMN substates (Supplementary Figure), it was associated with a significant decrease in the switching rate -in other words, it increased the duration -of DMN substate 2 in this analysis at IP (β = -0.010, SE = 0.003, pFDR = 0.020) (Figure). (Note, however, that fractional occupancy and switching rate were generally correlated with each other, across states and timepoints [r = 0.59].) Ibogaine was also associated with a significant reduction in the switching rate of DMN substate 3 in this analysis at both IP (β = -0.006, SE = 0.002, pFDR = 0.038) and 1M (β = -0.007, SE = 0.003, pFDR = 0.020). The decrease in DMN substate 3 switching rate remained significant after including head motion, age, and number of TBIs as covariates in the LME (IP vs BL: β = -0.007, SE = 0.003, z = -2.80, pFDR = 0.005; 1M vs BL: β = -0.008, SE = 0.003, z = -2.86, pFDR = 0.005) (Supplementary Figure). Head motion was not correlated with switching rates (Supplementary Figure). In summary, the DMN exhibited a significant increase in the duration of a particular dynamic substate, state 3, at both timepoints. The change in the session-normalized switching rate of DMN substate 3 at IP was significantly and negatively correlated with the mean MEQ-30 score (r = -0.457, pFDR = 0.033) (Figure; correlations between all session-normalized dynamic metrics and all subjective/clinical outcomes are shown in Supplementary Table). The correlation at 1M was still negative but not significant (r = -0.271, pFDR = 0.200). The mean activation of DMN substate 3 was higher in lateral than medial regions. This lateral-medial axis was also significantly aligned with a gradient of externally oriented to internally oriented cognition (r = 0.346, pspatial = 0.045); that is, activation of DMN substate 3 tends to be higher in DMN regions that regulate cognitive processes related to perception of the external environment (e.g., visual attention) (Figure). Meanwhile, activation of DMN substate 3 is lower in DMN regions involved in mediating internal cognitive phenomena (e.g., autobiographical memory). Exploratory analyses revealed that the change in the session-normalized switching rate of DMN substate 3 at the 1M timepoint was significantly and positively correlated with all clinical outcomes at 1M, i.e., scores on the CAPS-5 (r = 0.460, pFDR = 0.048), MADRS (r = 0.567, pFDR = 0.015), HAM-A (r = 0.700, pFDR = 0.001), and WHODAS scales (r = 0.548, pFDR = 0.018), after FDR correction (Figure). (Note that the four clinical outcomes may themselves be highly inter-correlated.) In other words, the larger the increase in session-normalized duration of DMN substate 3 (conversely, the greater the decrease in session-normalized switching rate of DMN substate 3) following ibogaine, the greater the clinical improvement on all four scales. None of the other four DMN substates exhibited significant correlations with clinical outcomes at 1M. The change in this switching rate was also significantly and positively correlated with improvements in HAM-A (r = 0.627, pFDR = 0.015) and WHODAS (r = 0.551, pFDR = 0.035) scores at 12M (Supplementary Figure). There was no significant correlation between TBI symptoms at baseline with changes in the session-normalized switching rate of DMN substate 3 at either IP or 1M (number of TBIs: rIP = -0.02, r1M = 0.01; BATL: rIP = -0.16, r1M = -0.29). HMMs trained on the whole-brain timeseries revealed states that were aligned with canonical restingstate networks, including the DMN (Supplementary Figure). However, none of these whole-brain states exhibited significant changes in switching rate at either IP or 1M (Supplementary Figure). The switching rate of the whole-brain state that most strongly resembled the DMN did not correlate significantly with MEQ-30 mean score or scores on any of the four clinical scales at either IP or 1M (Supplementary Figure).
SECTION 4. DISCUSSION
To our knowledge, this is the first study to characterize the post-acute effects of ibogaine on the dynamics of the DMN. Our findings make three central contributions. First, dynamic measures of DMN activity were more sensitive to ibogaine than static measures, both in differentiating pre-from post-treatment brain activity and in correlating with clinically meaningful outcomes. While static measures can provide valuable insights into time-averaged connectivity, dynamic analyses provide complementary information by revealing the evolution of DMN activity over time. Second, we identified a specific dynamical substate of the DMN, state 3 of this analysis, whose switching rate was significantly reduced after ibogaine. This switching rate was correlated with both the intensity of the mystical experience and clinical improvement across multiple symptom domains at one month and at long-term follow-up. Third, this same DMN substate exhibited a lateral-medial spatial gradient, which may be aligned with an axis of externally-to-internally oriented perception and cognition.
SECTION 4.1. DYNAMIC MEASURES ARE COMPLEMENTARY TO STATIC MEASURES OF IBOGAINE'S POST-ACUTE EFFECTS ON THE DMN
Ibogaine was not associated with significant changes in static within-DMN FC or static betweennetwork FC at either IP or 1M. This null result converges with our group's prior work on the same cohort: Sridhar and colleagues (2026), using exploratory static FC analyses between pairs of regions, observed widespread changes in static FC across the DMN and other brain networks, yet these did not survive correction for multiple comparisons.However, the present study is the first to show that static within-DMN FC at 1M was significantly associated with mean scores on the MEQ-30. The absence of post-acute changes in static within-DMN FC, as well as its positive correlation with mystical experiences, may appear to be in tension with prior literature. Indeed, previous studies have shown that psychedelics reduce FC within the DMN and that mystical experiences are negatively correlated with within-DMN connectivity. 11,60-63 However, these studies acquired fMRI data during the psychedelic experience, whereas the fMRI data in the MISTIC trial was collected 3-4 days and one month after ibogaine dosing. (Given the 28-to 49-hour half-life of noribogaine, the primary active metabolite of ibogaine, it is possible that noribogaine was still present in the brain at the IP timepoint.) Furthermore, there is evidence that the post-acute effects of psychedelics on static FC can be the opposite of their acute effects. For example, one day after psilocybin treatment for depression, static within-DMN FC rises, yet it decreases during the treatment.(Note that mystical experiences were measured retrospectively, at the IP timepoint.) Notably, the one significant change in static FC was not in the DMN but in the salience network, which showed a reduction in within-network static FC at the IP, but not 1M, timepoint. The salience network is a core hub of threat detection, interoceptive awareness, and hyperarousal, and is hyperconnected in patients with PTSD.Its acute disruption following ibogaine may suggest that ibogaine relaxes the hypervigilance that characterizes PTSD, which is supported by the decreases in peak alpha frequency that were observed after ibogaine treatment in electroencephalography (EEG) data from the same veterans.However, salience network FC returned to baseline levels at 1M, so it cannot be a marker of the long-term effects of ibogaine. By contrast, dynamic FC analysis using HMMs revealed robust, significant effects of ibogaine treatment on the DMN. Ibogaine was associated with significant decreases in the switching rate of DMN substate 3 at both the IP and 1M timepoints, an effect that survived FDR correction across all five DMN substates and two post-treatment timepoints and was robust to inclusion of head motion, age, and number of TBIs as covariates. Crucially, after FDR correction, change in the switching rate of DMN substate 3 was also correlated with the intensity of the mystical experience and with clinical improvements at 1M on the CAPS-5, MADRS, HAM-A, and WHODAS-2.0, as well as with longterm changes in HAM-A and WHODAS-2.0 scores at 12M. Static within-DMN FC, by contrast, was not significantly correlated with any clinical outcome at any timepoint, although it did correlate with the mystical experience at 1M. Dynamic FC patterns often diverge strongly from static FC, since resting-state FC cycles through different states over time while static measures treat brain networks as if they occupied a single configuration across the entire scanning sessionKey hubs of the DMN in particular, such as the posterior cingulate cortex, exhibit multiple dynamical co-activation patterns, each one engaging a distinct set of brain regions.Although HMMs are typically applied to fMRI data of the entire brain, prior work has also fitted HMMs just to the timeseries of DMN regions and recovered discrete DMN connectivity states with distinct dynamical properties.Most studies of FC in the psychedelic literature have used static measures, but dynamic FC measures are capable of capturing more information.For instance, only dynamical analyses can reveal whether the 'repertoire' of FC states -specific patterns of connectivity between selected regions -expands under psychedelics; static measures only identify a single set of connections over the entire duration of the scan.As shown in this study, dynamical measures can also indicate changes in the switching rate or duration of certain states, which could occur even when the time-averaged signal remains the same.Post-acute fMRI activity on ibogaine only appears to exhibit significant changes in FC when measured dynamically rather than statically, highlighting the complementary value of dynamic measures. In summary, static FC appears to capture meaningful individual variability related to mystical experiences, whereas dynamic analyses capture the temporal organization of DMN activity that relates to both mystical experiences and clinical improvement.
SECTION 4.2. SUBSTATE-SPECIFIC CHANGES WITHIN THE DMN
A second contribution of this work is both methodological and conceptual: the prior psychedelic neuroimaging literature has, with few exceptions, treated the DMN as a single, monolithic entity. That is, studies typically report just one within-DMN FC value or a single DMN-component time series; hence, the DMN as a whole goes 'up' or 'down' in connectivity, which obscures more nuanced, spatially selective changes in FC. In our study, only one of the five DMN substates, state 3, showed a significant change in switching rate both immediately and one month after ibogaine. Were the DMN analyzed only as a single composite signal, the ibogaine-induced change in DMN substate 3 would have been diluted by the unchanged dynamics of the other four DMN substates, plausibly contributing to the null static-FC findings. The whole-brain HMM analysis reinforces this point. When the HMM was fit to the full 200-parcel time series rather than restricted to DMN parcels, the whole-brain state that best resembled the entire DMN did not show significant changes in switching rate after ibogaine, and none correlated with the mystical experience or clinical outcomes. Identifying DMN substates with a DMN-restricted HMM was necessary to identify dynamics that differentiated brain activity before and after ibogaine. The gradient analysis offers a principled explanation of why DMN substate 3, specifically, should be the substate that tracks mystical experience. DMN substate 3 was characterized by relatively higher activity in lateral DMN regions and lower activity in medial DMN regions. While the core, medial hubs of the DMN are associated with internally oriented, self-referential cognition, lateral regions outside these hubs encode different behaviors.In particular, more lateral regions, including lateral temporal cortex, the inferior parietal lobule, and middle temporal gyrus, are recruited by more externally directed semantic, attentional, and social tasks, while medial DMN areas deactivate in response to these tasks.Decreased switching rate of, or increased duration of, DMN substate 3 may therefore imply that ibogaine recruits lateral DMN regions that are involved in externally oriented cognitive processes (e.g., visual attention, action perception, auditory perception) while reducing activation of medial DMN regions that mediate internally oriented processes (e.g., autobiographical memory, introspection, mental imagery). This reconfiguration of the DMN is consistent with the phenomenology of mystical experiences, in which self-focused awareness recedes and people feel "one" with the surrounding environment.Although the acute effects of ibogaine involve internal, dreamlike experiences, the integration in the days and months afterwards, when we recorded fMRI, often leads to heightened awareness of and re-engagement with the external environment.However, the involvement of lateral DMN regions in externally-oriented functions and medial DMN regions in internally-oriented processes was not directly verified by experiment in this study. Instead, the gradient analysis relied on Neurosynth, which performs a meta-analysis of the literature to determine brain regions that are correlated with various behaviors. The functional relevance of the medial-lateral gradient of DMN substate 3 requires further validation. Beyond its apparent association with the mystical experience, the change in switching rate of DMN substate 3 was correlated with clinical improvement across multiple symptom domains: PTSD, depression, anxiety, and functional disability. Critically, no such associations were observed for static within-DMN FC, nor for any of the other four DMN substates or for the DMN-resembling state of the whole-brain model. In other words, the clinical relevance of ibogaine's post-acute effects on the DMN was specific to the temporal dynamics of a single DMN substate, rather than the network's time-averaged connectivity or its dynamics as a whole. DMN substate 3's correlation with both the mystical experience and clinical outcomes is consistent with prior work showing that the subjective effects of psychedelics are related to their therapeutic benefits.That is, the more intense the subjective experience, the greater the clinical improvement. Although PTSD, depression, and anxiety are typically associated with aberrations in static within-DMN FC, our results indicate that increases in the duration of a particular DMN substate are more predictive of symptom improvement after ibogaine treatment than static measures averaged across the entire DMN.Overall, the substatespecific dynamics of the DMN appear to be a robust neural correlate of the therapeutic and mystical effects of ibogaine.
SECTION 4.3. LIMITATIONS AND FUTURE DIRECTIONS
Several limitations should be acknowledged. First, ibogaine treatment was open-label; future studies should consider double-blinded or dose-controlled study designs. Currently, we cannot disentangle the effects of ibogaine from expectancy, set, and setting, especially in light of the fact that the ibogaine treatment center offered other wellness activities such as sweat lodge, massage, yoga, reiki, breathwork and meditation. Second, our sample was all-male and almost completely White. Findings may not generalize to women, to non-White populations, to non-veteran patients, or even other types of veterans (including non-SOVs). Thirdly, this study did not consider other methods of measuring DMN dynamics, such as leading eigenvector dynamics analysis, co-activation patterns, or sliding window FC. Fourthly, the gradient analysis is not robust to different selection quantiles (i.e., thresholds at which DMN regions were evaluated for their involvement in external and internal processes) and may not be robust to different anchor terms. Finally, future studies should acquire fMRI data both during and after the ibogaine experience so that the post-acute and acute neural signatures can be directly compared.
SECTION 5. CONCLUSION
In a cohort of veterans with TBI, dynamic measures captured the reorganization of DMN activity after ibogaine treatment. Ibogaine was specifically associated with an increase in the duration of a single DMN substate whose duration correlated with the intensity of the mystical experience and clinical improvement across PTSD, depression, anxiety, and functional disability. These results argue against treating the DMN as a temporally and spatially uniform network in psychedelic neuroimaging, and suggest that the post-acute neural signature of ibogaine treatment may be best resolved at the level of specific dynamical DMN substates.
DATA AND CODE AVAILABILITY
• Data: The deidentified human participant data reported in this study cannot be deposited in a public repository because they contain sensitive clinical and neuroimaging data from a vulnerable veteran population and are subject to Institutional Review Board-mandated controlled access. To request access, please contact the Stanford University Institutional Review Board and the corresponding author (Cammie Rolle, crolle@stanford.edu) with a research proposal and data security plan consistent with Stanford data-sharing requirements. • Code: Original code for the manuscript comprises the internal-external gradient analysis, which can be found here:. Otherwise, code used existing functions. (b) This correlation is still negative but not significant at 1M.
Full Paper Sections
Go Pro to read the extracted full-text sections.
Full Text PDF
Full Paper PDF
Create a free account to open full-text PDFs.
Study Details
- Study Typeindividual
- Populationhumans
- Characteristicsobservationalopen labelsurveybrain measures
- Journal
- Compound
- Topics
- APA Citation
- Citation FormatsExport citation
Linked Clinical Trial
Registry record linked to this paper.
Linked Clinical Trials
Create a free account to see the clinical-trial registry records linked to this paper.
References (23)
References cited by this study and indexed in Blossom.
Gattuso, J. J., Perkins, D., Ruffell, S. G. D. et al. · International Journal of Neuropsychopharmacology (2022)
Palhano-Fontes, F., Andrade, K. C., Tófoli, L.F. et al. · PLOS ONE (2015)
Timmermann, C., Roseman, L., Haridas, S. et al. · PNAS (2023)
Girn, M., Doss, M. K., Roseman, L. et al. · Nature Medicine (2026)
Barrett, F. S., Griffiths, R. R. · Current Topics in Behavioral Neurosciences (2017)
Carhart-Harris, R. L., Friston, K. J. · Pharmacological Reviews (2019)
Cherian, K. N., Keynan, J. N., Anker, L. et al. · Nature Medicine (2024)
Brown, R. E., Lissemore, J. I., Shinozuka, K. F. et al. · Journal of Affective Disorders (2026)
Coleman, J. A., Yang, D., Zhao, Z. et al. · Nature (2019)
MacLean, K. A., Leoutsakos, J. S., Johnson, M. W. et al. · Journal for the Scientific Study of Religion (2012)
Show all 23 referencesShow fewer
Carhart-Harris, R. L., Muthukumaraswamy, S., Roseman, L. et al. · PNAS (2016)
Carhart-Harris, R. L., Erritzoe, D., Williams, T. et al. · PNAS (2012)
Smigielski, L., Scheidegger, M., Kometer, M. et al. · NeuroImage (2019)
Glue, P., Lockhart, M., Lam, F. et al. · Journal of Clinical Pharmacology (2014)
Carhart-Harris, R. L., Roseman, L., Bolstridge, M. et al. · Scientific Reports (2017)
Lissemore, J. I., Chaiken, A., Keller, C. J. et al. · Nature Mental Health (2025)
Barrett, F. S., Krimmel, S. R., Griffiths, R. R. et al. · NeuroImage (2020)
Madsen, M. K., Stenbaek, D. S., Arvidsson, A. et al. · European Neuropsychopharmacology (2021)
Shinozuka, K., Jerotic, K., Mediano, P. A. M. et al. · Translational Psychiatry (2024)
Mortaheb, S., Fort, L. D., Mason, N. L. et al. · Biological Psychiatry (2023)
Lord, L. D., Expert, P., Atasoy, S. et al. · NeuroImage (2019)
Schenberg, E. E., de Castro Comis, M. A., Alexandre, J. F. M. et al. · Journal of Psychedelic Studies (2017)
Ross, S., Bossis, A. P., Guss, J. et al. · Journal of Psychopharmacology (2016)
Your Personal Research Library
Go Pro to save papers, add notes, rate research, and organize custom shelves.