Pharmacokinetics and pharmacodynamics of sublingual microdosed lysergic acid diethylamide in healthy adult volunteers
In a double‑blind Phase 1 trial in 80 healthy male volunteers, a one‑compartment population pharmacokinetic model and LC‑MS/MS assay characterised 10 µg sublingual LSD (Cmax 0.20 µg/L, Tmax 1.51 h, t1/2 3.08 h, clearance 7.78 L/h/70 kg, Vc 32.9 L/70 kg). Microdosing produced minimal cardiovascular and subjective effects, no change in peripheral BDNF, and suggested qualitative influences of CYP genotypes (notably CYP2D6) on concentrations, indicating the need for larger genotype studies and more sensitive pharmacodynamic measures.
Authors
- Suresh Muthukumaraswamy
Published
Abstract
Introduction: Microdosing is the practice of taking psychedelic drugs at doses that produce no or minimal perceptible subjective or behavioural effects. This study investigated the pharmacokinetics and pharmacodynamics of microdosed lysergic acid diethylamide (LSD). Methods: This was a Phase 1 double-blind placebo-controlled parallel-groups trial with 80 healthy male volunteers (four withdrawals due to anxiety). Plasma samples were taken at 0.5, 1, 2, 4 and 6 h after 10 µg sublingual LSD and analysed with liquid chromatography-tandem mass spectrometry (LC-MS/MS). LSD pharmacokinetics were modelled. Population analyses were performed using nonlinear mixed effects models. Heart rate and a visual analogue scale (‘feel effect’) were used to describe LSD pharmacodynamics. The effect of the relevant cytochrome P450 (CYP) genotype on LSD pharmacokinetics was qualitatively assessed. Plasma and serum levels of brain-derived neurotrophic factor (BDNF) were evaluated. Results: A one-compartment model best described LSD pharmacokinetics. Mean (95% confidence interval): elimination clearance = 7.78 L/h/70 kg (6.75–8.77), central volume of distribution = 32.9 L/70 kg (30.1, 36.0). Maximal concentration (0.20 µg/L), time to maximal concentration (1.51 h) and elimination half-life (3.08 h). The maximal increase in heart rate and visual analogue scale was small (<15%) compared to baseline estimates limiting the modelling. Two of the participants withdrawn from the study due to anxiety had intermediate-weak CYP2D6 activity. CYP2D6, CYP1A6, CYP2B6 and CYP2C9 qualitatively appeared to influence concentration. No evidence of alterations of peripheral BDNF with microdosing was found. Conclusion: This study provides a population pharmacokinetic model and LC-MS/MS assay that can inform clinical and bioequivalence studies. Relevant CYP genotypes should be studied in larger samples as combined potential biomarkers of response. Microdose-sensitive and reliable pharmacodynamic measures are needed.
Research Summary of 'Pharmacokinetics and pharmacodynamics of sublingual microdosed lysergic acid diethylamide in healthy adult volunteers'
Introduction
Microdosing with LSD refers to repeated administration of sub‑perceptual doses (commonly about 10 µg) that users report may improve mood, cognition and other symptoms. Earlier experimental work in healthy volunteers has shown dose-dependent effects of low doses (5–26 µg) on physiological and behavioural measures, but pharmacokinetic (PK) characterisation at microdose levels, especially for sublingual administration, is limited. There are no published population PK models for sublingual microdoses of LSD, and uncertainty remains about concentration–effect relationships, appropriate pharmacodynamic (PD) endpoints sensitive to microdoses, and the influence of metabolic genotype (notably CYP2D6) on exposure and subjective effects. Morse and colleagues set out to describe the pharmacokinetics and pharmacodynamics of a 10 µg sublingual LSD microdose in healthy adult male volunteers, develop and validate a high‑sensitivity LC‑MS/MS assay for low plasma concentrations (including the major metabolite O‑H‑LSD and iso‑LSD), and to explore whether CYP genotype and peripheral brain‑derived neurotrophic factor (BDNF) concentrations change acutely or after a 6‑week microdosing regimen. The study aimed to provide a population PKPD model that could inform dose selection, bioequivalence studies and future clinical development of LSD microdosing.
Methods
This was a Phase I, double‑blind, placebo‑controlled, parallel‑groups trial conducted at the University of Auckland Clinical Research Centre. Eighty healthy male volunteers aged 25–60 years were randomised 1:1 to receive either 10 µg sublingual LSD or inactive treatment (placebo). Key exclusions included elevated resting blood pressure, significant medical or psychiatric disorders, recent psychedelic use (past year) or any prior microdosing. Participants provided informed consent and the protocol was registered and ethically approved. On the first treatment visit (morning sessions), an intravenous cannula was placed and participants self‑administered 1 mL syringes containing either 10 µg LSD in distilled water (0.2% ethanol) or distilled water placebo held sublingually for 30 seconds before swallowing. A titration regimen (reducing by 5 µg then increasing by 1 µg steps as tolerated) was available for participants who reported overstimulation; six entered titration (five active, one placebo) and four active participants were subsequently withdrawn due to anxiety/overstimulation. Blood samples for PK analysis were collected at 30, 60, 120, 240 and 360 minutes post‑dose; plasma and serum for BDNF and DNA for genotyping were also collected at baseline and at a final visit after 6 weeks of home microdosing (only the first‑dose PK/PD data are the focus of this report). LSD and metabolites were quantified using a validated LC‑MS/MS method developed for low plasma concentrations; assay reproducibility was demonstrated with re‑analysis of samples showing ≤15% difference. PKPD modelling used nonlinear mixed‑effects methods (NONMEM v7.5.1) with one‑ and two‑compartment structures considered, first‑order absorption/elimination, allometric scaling of body weight (exponent 3/4 for clearances, 1 for volumes), and an effect‑compartment to allow for delay between plasma concentration and effect. PD endpoints were heart rate and a visual analogue scale (VAS) item capturing a subjective 'feel effect'. Residual error models included additive and proportional components, and model selection used objective function value changes and visual predictive checks; bootstrap (100 replications) assessed parameter uncertainty. Pharmacodynamic statistical analyses used linear mixed‑effects models (lmerTest in R) with Group (LSD vs placebo) and Time as fixed effects and participant as a random effect, with interest in the group × time interaction. CYP genotyping was performed on extracted DNA using Agena MassARRAY with the Veridose panels to assess CYP2D6, CYP2C19, CYP1A2, CYP2B6, CYP2C9 and CYP3A4 and to call activity scores per contemporary CPIC guidance; copy number variants (CNVs) were also assessed. Plasma and serum BDNF were measured by ELISA (R&D Systems Quantikine), samples assayed in duplicate, with outliers and values below limit of detection retested; intra‑ and inter‑assay variability were reported ≤6.1%.
Results
Seventy‑nine participants contributed blood samples for PK analyses (one participant could not be cannulated). Using the LC‑MS/MS assay, LSD was quantified in 100% of collected samples from active participants and the intended sample collection rate for the active group was 96.5% (two missing samples at 120 and 360 minutes). A one‑compartment model provided the most parsimonious description of the sublingual 10 µg LSD PK data; a two‑compartment model reduced the objective function value (ΔOFV = 11.4, p = 0.003) but was judged over‑parameterised with poorly estimated intercompartmental parameters and no improvement in visual predictive checks. Final population parameter estimates (with reported variability) yielded a central volume of distribution around 32.9 L/70 kg and an elimination clearance approximately 7.78 L/h/70 kg. Derived metrics were: mean Cmax ≈ 0.20 µg/L (204 pg/mL; reported ±0.13), Tmax ≈ 1.51 h (±0.66), and elimination half‑life ≈ 3.08–3.09 h (±0.37). Allometric scaling by body weight was applied in the model and visual predictive checks and bootstrap resampling (100 replicates) were used to evaluate model performance and parameter uncertainty. Pharmacodynamic effects at this microdose were minimal. Heart rate and VAS 'feel effect' showed small maximal increases (<15% relative to baseline), which limited the ability to robustly characterise concentration–effect relationships; effect‑compartment modelling was used but estimated EMAX parameters remained small. Mixed‑effects analyses examined group × time interactions for PD measures but the low magnitude of changes and a relatively strong placebo response reduced discriminability between groups. The authors note a higher placebo response compared with some prior microdose studies and discuss potential reasons (parallel design, demographic differences, prior psychedelic experience). CYP genotyping in the active group (n = 40) showed 21/40 participants with CYP2D6 extensive metabolism, 13/40 intermediate metabolisers and no clearly classifiable poor or ultrarapid metabolisers based on core diplotypes; five participants had CNVs or hybrid alleles complicating phenotype assignment. Two participants were classified as intermediate‑weak (activity score 0.5) and both were among the four active participants withdrawn for overstimulation/anxiety during the home microdosing phase. Exploratory plots comparing concentration–time curves by genotype suggested qualitative effects for several enzymes: no clear CYP2C19 influence, possible lower concentrations in some CYP1A2 (rs762551) A/A carriers consistent with induction, and trends in expected directions for CYP2B6 and CYP2C9 variants though numbers were small. The dataset lacked sufficient numbers of poor or ultra‑rapid CYP2D6 phenotypes to permit formal covariate modelling of activity scores. BDNF analyses: for plasma, 214/240 (89%) observations were analysable after excluding missing samples, below‑limit values and outliers. Linear mixed models showed no group × time interaction at 6 h (t = 0.14, p = 0.89) or at 43 days (t = 0.40, p = 0.68). There was an effect of time at 6 h across groups (t = 2.15, p = 0.03) but no group effect. For serum BDNF, 224/240 (93%) samples were analysable; models showed no interaction at 6 h (t = 0.234, p = 0.81) but a statistically significant interaction at 43 days (t = 2.45, p = 0.016). The authors attribute that interaction to a low baseline mean in the placebo group rather than an increase in the LSD group; reported means were: placebo pre‑mean 26,496 (SD 6,676) pg/mL, placebo post‑mean 28,528 (6,661) pg/mL, LSD pre‑mean 28,412 (6,284) pg/mL, LSD post‑mean 28,482 (5,767) pg/mL. Overall, the study found no convincing evidence that a single 10 µg microdose or a 6‑week microdosing regimen altered peripheral BDNF in healthy male volunteers.
Discussion
Morse and colleagues interpret their findings as providing a validated population PKPD framework and a high‑sensitivity LC‑MS/MS assay suitable for microdosing studies. The one‑compartment PK description, estimated clearance and volume of distribution, and derived Cmax, Tmax and half‑life are consistent with prior microdose work and earlier reports of LSD PK. The presented assay detected LSD and O‑H‑LSD in plasma at microdose concentrations and was reproducible across re‑analysed samples. Pharmacodynamic effects at the 10 µg sublingual dose were small, with heart rate and self‑reported 'feel effect' VAS changes under 15%, limiting the capacity of the PKPD model to characterise robust concentration–effect relationships. The investigators highlight a relatively strong placebo response and note that the parallel‑groups design may have improved blinding compared with crossover designs, potentially contributing to smaller between‑group PD differences. Demographic differences and prior psychedelic experience compared with other studies are also discussed as factors that may explain differing placebo and active responses. Genotyping results suggest that CYP2D6 and other CYP enzymes may influence LSD concentrations qualitatively, but the sample contained no clearly classifiable poor or ultra‑rapid CYP2D6 metabolisers and only two intermediate‑weak participants. The two intermediate‑weak participants were among those withdrawn for overstimulation in the home‑dosing phase, a finding the authors consider potentially relevant if replicated in larger samples: CYP2D6 activity scores might be useful biomarkers to guide dosing or to stratify risk of adverse effects. The authors caution that multiple CYPs and CNVs contribute to a complex genotype profile and that larger studies are needed to quantify the average contribution of specific SNPs and CNVs to plasma exposure. Regarding BDNF, the study found no evidence that a single 10 µg microdose or a 6‑week regimen increased peripheral BDNF in healthy men. The authors contrast this null result with prior small studies reporting increases but note those earlier studies had substantial missing data and small effective sample sizes. High baseline variability in peripheral BDNF and differences between plasma and serum matrices are discussed; most peripheral BDNF is platelet‑stored and serum processing releases platelet BDNF, so plasma (with fewer platelets) might better reflect acute effects, but the evidence remains inconsistent. The investigators therefore recommend larger, within‑subject designs and careful sample handling to probe peripheral BDNF as a biomarker. Key limitations acknowledged include the exclusive inclusion of males (restricting generalisability and not accounting for menstrual cycle effects), limited sampling during the elimination phase and lack of diversity in CYP2D6 phenotypes (few or no poor/ultrarapid metabolisers and several CNV/hybrid genotypes complicating classification). The authors call for further studies to identify PD measures sensitive to microdoses, to include females with attention to menstrual phase and life stage, and to extend PKPD modelling across formulations, doses and populations so that dosing regimens can be rationally developed for potential therapeutic use.
Conclusion
This study delivers a validated LC‑MS/MS assay for low plasma concentrations of LSD, a population one‑compartment PK model for a 10 µg sublingual microdose, and preliminary PKPD and pharmacogenetic observations. No convincing acute or 6‑week effects on peripheral BDNF were detected in healthy male volunteers. The authors conclude that CYP2D6 activity score merits further investigation in larger samples as a potential biomarker of response or adverse effects, and that future work should prioritise sensitive PD metrics, inclusion of females with menstrual considerations, and larger pharmacogenomic studies to inform dose titration and personalised microdosing regimens.
View full paper sections
METHODS
A phase 1 trial with double-blind placebo-controlled parallel groups was used to generate pharmacokinetic, pharmacodynamic and BDNF observations for analysis. The trial protocol was prospectively published, and primary results were reported. The trial was registered at ANZCTR (trial ID = 381476). All participants provided written information consent, and ethical approval was provided by the New Zealand Health and Disability Ethics Committee (reference number: 19/STH/91). Lab visits in the study took place at the Clinical Research Centre at the University of Auckland.
RESULTS
Linear mixed-effects modelling of the pharmacodynamic measurements was undertaken using the lmerTest package in Rwith Group (LSD, placebo) and Timebeing treated as fixed effects with dummy coding, and participants as a random effect, with the primary estimates of interest being the group × time interaction effect.
CONCLUSION
The potential for microdoses of LSD to be developed as medicines for mental health disorders has renewed interest in the PK and PD of LSD given in different formulations. In addition, understanding the factors that may mediate response and tolerance to LSD may pave the way to precision medicine approaches to using microdosing of LSD in treating mental illness. This study validates the findings ofby providing the accompanying blood sample concentration of LSD. We successfully collected 96.5% of the samples we intended, to quantify LSD concentrations. We were able to quantify LSD in 100% of the analysed samples. We present the highsensitivity LC-MS/MS assay we used for the detection of low concentrations of LSD in plasma along with validation data. The assay is suitable for use in microdosing studies. Furthermore, the assay demonstrated the detection of OH-LSD. Although it was unlikely that OH-LSD metabolite would be detected in plasma samples following a microdose of LSD, this method could be used as a starting point for future plasma/urinary analyses that require the quantification of this metabolite. We demonstrate that a one-compartment model was found to best describe the pharmacokinetics of microdosed sublingual LSD, consistent with previous studies that used single or noncompartment modelling for microdosing with LSD. One previous study on intravenous LSD used a two-compartment model. We had insufficient samples of the elimination phase to consider a two-compartment mode as in. Unfortunately, even if we had these samples the overall lack of detectable effect at low concentrations of LSD (in our VAS and vital sign measures of PD) would still limit further inference on the nature of the PK compartments (in physiological terms). Our analysis extends those in the literature by accounting for size-related changes in LSD pharmacokinetics using theorybased allometric scaling of total body weight and quantifying the influence of inactive treatment on effect. Pharmacokinetic parameters such as clearance and volume of distribution are consistent with previous modelling studies. The C max of 204 pg/mL was consistent with, who also reported a 10 µg dose (their C max was 151 pg/mL). T max at 1.5 h sat midway between the studies led by Liechti and his group at the University Hospital of Basel(Ms = 1.1-1.7 h). By contrast, one other group has reported fast T max = 50.4 min. T 1/2 of 3 h was consistent with several previous papers. We used heart rate and VAS scales (specifically whether the participants were able to 'feel effect') to describe LSD pharmacodynamics. The maximal increase in heart rate and VAS were small (<15%) compared to baseline estimates; the lack of change is likely attributed to the low dose of LSD administered and the consequent concentrations. Furthermore, one noticeable difference between the VAS data we report compared (Figure) to the data presented inis a much stronger response in the placebo group. This could be due to several, possibly interacting factors. Firstly, our study used a parallel-groups trial design which is better able to blind participants compared to a crossover trial where participants can use their psychological responses in previous sessions to help them guess the identity of subsequent intervention sessions. Secondly, demographic factors may be important as Holze et al. (2021a) used a young (M = 23 years, SD = 3) mixed-sex cohort where all participants had previous psychedelic use experience, whereas the MDLSD trial had a wider age range (M = 36, SD = 7) of male participants of whom 30% in each group were psychedelic naïve. Incorporating covariates (such as genotypes or phenotypic traits) that can be used to predict response in the individual patient is one of the most sought-after outcomes of modelling in addition to determining PD measures that are sensitive to microdoses. CYP2D6 has been studied extensively for the potential usefulness of classifying metabolism status by activity score. A poor, intermediate or ultra-rapid status has been reported as potentially useful in guiding clinical decision-making. Individuals with CYP2D6 poor metaboliser status have increased LSD plasma concentrations and decreased LSD clearance. The only participant with a CYP2D6 activity score of 0 (indicating a poor metaboliser) from their core panel-derived diplotype also had a copy-number variation meaning they could not be accurately classified from their diplotype. We did not incorporate activity scores into the modelling because of these low numbers. However, there were two participants in our cohort who were intermediate-weak (activity score = 0.5) metabolisers and they were both among the four participants who were withdrawn from the study due to over-stimulation/anxiety (Figure). This finding is potentially relevant to the design of LSD microdosing therapeutic interventions if the tendency for overstimulation to occur in people with activity scores 0-0.5 was replicated in larger samples. In such cases, CYP profiles could be added to models. It is also possible in the future that potential patients entering into an LSD microdosing regimen could be genotyped for CYP2D6 genotype status prior to starting a treatment regimen and be started at a lower dose to reduce adverse effects if they are found to be poor metabolisers. Visual inspection of Figure(c) and (d) also reveals the importance of LSD plasma concentration in terms of generating over-stimulation effects during home microdosing protocols. This suggests that dose titration/escalation protocols should be used to reduce adverse effect incidence and allow patients to find an individually optimised dose-an approach we are taking in ongoing Phase 2 trials. Overall, it is apparent that more than CYP2D6 activity scores are needed to predict LSD concentration given the wide variability in individual concentration, including around Figure(b) intermediate-weak metabolisers. From our results, CYP2C19 clearly had no effect on the concentration following the first microdose of LSD. However, CYP2C19, CYP1A2 and CYP3A4 are induced and so in a less controlled sample (i.e. in a clinical or recreational setting), these genotypes are likely to become more important, particularly 1A2 and 3A4. CYP1A2 rs762551 variant carriers did have qualitatively higher concentrations of LSD (Figure) which may reflect some induction. CYP2B6 rs3745274 and CYP2C9 variant carriers showed either a less clear spread, or low numbers, respectively, so were less compelling, but were trending in predicted directions. CYP3A4 had two extreme values in the variant carriers (the highest AUC and the fifth lowest). As summarised in Table, individual participants had a complex profile of CYP genotypes. The data support larger studies with similar panels, where the average contribution of different SNPs to drug concentration in plasma could be quantified, and it could be determined if it was useful to inform titration regimens for LSD microdosing.report that CYP2E1 is also important to LSD metabolism and should be included in future studies. The current study found no evidence to support the hypothesis that LSD microdosing causes an increase in peripheral BDNF concentrations 6 h after a 10 µg dose nor after completing a 6-week regimen of microdosing. One previous studyfound increases in plasma BDNF concentrations in the hours following microdosing. However, despite that trial containing 24 participants, many samples were not available for assay. Statistical analyses were conducted using a complete case analysis of n = 10 for the 5 µg dose, n = 9 for the 10 µg dose and n = 8 for the 20 µg dose with a 37% data availability rate. The lack of available data may have increased the bias within the study and consequently, results may not reflect the true population time course of BDNF. Conversely, this study had a larger sample size n = 40 with high data availability (~90%). From the data, it is clear that there is high variability even at baseline (Figure)reinforcing the need for high numbers and within-subject analyses when assessing peripheral BDNF. Furthermore, more recent studies with macrodosing have largely failed to see an increase in plasma or serum BDNF following 100-200 µg doses. This study measured BDNF concentrations in both plasma and serum as it is not clear which sample type might be most sensitive to psychedelics. Most peripheral BDNF is stored in platelets (which cannot cross the blood-brain barrier) and as such plasma BDNF might be a better proxy of acute drug effect changes. The coagulation processes that occur during serum sample processing causes BDNF to be released from platelets into serum allowing this pool of BDNF to be accessed. Plasma does still contain some platelets and as yet no studies with psychedelics have measured platelet-poor plasma in which these extra platelets have been removed by additional centrifugation. The current study did observe plasma BDNF concentrations increased only at the 6-h timepoint in both groups, an effect which might be due to some non-specific study effect. Particularly because in males, plasma BDNF is expected to decrease from morning to afternoon, and serum BDNF is expected to remain stable, which it did. Despite this study not identifying the effects of LSD microdosing on BDNF concentrations, clinical populations such as those with depression (rather than healthy volunteers) may exhibit a different time course of BDNF since that population has been shown to have lower peripheral BDNF concentrations.
Study Details
- Study Typeindividual
- Populationhumans
- Characteristicsre analysisplacebo controlledrandomizedparallel groupdouble blind
- Journal
- Compounds
- Topics
- Author