Serotonin syndrome by drug interactions with linezolid: clues from pharmacovigilance-pharmacokinetic/pharmacodynamic analysis

Purpose To characterize the post-marketing reporting of serotonin syndrome (SS) due to drug-drug interactions (DDIs) with linezolid and investigate the relationship with pharmacokinetic/pharmacodynamic (PK/PD) properties of serotonergic agents. Methods We queried the worldwide FDA Adverse Event Reporting System to extract SS records due to DDIs where linezolid was reported as suspect. For each serotonergic agent concomitantly reported, proportion of SS reports and mean number of DDIs were calculated and three different “SS reporting zones” were created. Relevant PK (peak concentration, area under plasma concentration curve, volume of distribution (VD), and lipophilicity) and PD (values of binding affinity (Ki) and IC50 for serotonin reuptake transporter (SERT) and 5-HT2A) parameters were extracted for each serotonergic agent, and relevant PK/PD indexes were calculated to assess correlation with mean number of DDIs (PV index). Results Six hundred sixty-nine reports of SS mentioning linezolid were found, being linezolid-citalopram (N = 69; 10.3%) the most frequently DDI reported. Citalopram and methadone showed respectively the highest proportion of SS reports (0.28%) and the lowest mean number of DDIs (1.41). Citalopram, escitalopram, and methadone emerged as red (i.e., alert)-zone medications: they exhibited high lipophilicity and large VD (proxies of excellent central nervous system penetration) coupled with high potency. Among PK/PD indexes, a significant correlation with PV index was found for VD/Ki SERT ratio (p = 0.05). Discussion Our integrated approach suggests that linezolid is more likely to cause SS when co-administered with citalopram, escitalopram, and methadone, as inferred from their pharmacological properties. Proper management of SS should be tailored on a case-by-case basis. Electronic supplementary material The online version of this article (10.1007/s00228-020-02990-1) contains supplementary material, which is available to authorized users.


Introduction
Linezolid is an oxazolidinone antibiotic with activity against multidrug-resistant Gram-positive organisms [1], showing lipophilic features, excellent tissue penetration including the central nervous system (CNS), and weak reversible non-selective monoamine oxidase (MAO) inhibitory effects at therapeutic serum concentrations (according to an inhibitory binding affinity constant (Ki) of 56 μM and 0.71 μM, respectively for MAO-A and MAO-B) [2,3]. MAO is involved in the metabolism of the monoamine neurotransmitters, and its inhibition may potentially lead to excess of serotonin (5-hydroxytriptamine (5-HT)) in the CNS and occurrence of serotonin syndrome (SS), a potential lifethreatening condition [4,5]. Inhibition of the serotonin reuptake transporter (SERT), causing accumulation of serotonin in the synaptic cleft and overstimulation of 5-HT 2A receptors, is supposed to be the main mechanism involved in SS [5].
Several drugs may cause serotonin excess through different pathways [6], thus posing concerns when concomitant treatments with linezolid are needed. Although reports describing occurrence of SS due to drug-drug interactions (DDIs) with linezolid were reported [3,7,8], and the Food and Drug Administration (FDA) warned against co-administration of linezolid with serotonergic drugs recommending a 2-week washout period in patients already receiving these medications [9], treatment of lifethreatening infections may not be delayed, especially if secondline alternatives are not effective or available. However, it remains unclear whether all serotonergic agents carry similar risk of SS, and the infectious disease consultant faces a "dilemma," with relevant implications on rational drug use in specific clinical scenarios, as exemplified in Fig. 1.
In order to address this aspect and provide clues to help in decision-making, we investigated the relationship between real-world data, namely, spontaneous reports submitted to the US FDA adverse event reporting system (FAERS) database, and pharmacokinetic/pharmacodynamic (PK/PD) properties of serotonergic agents.

Materials and methods
We carried out a "two-step" approach to investigate whether the reporting pattern of SS induced by linezolid is potentially related to the burden of co-medications resulting in clinically significant DDIs, and relevant PK/PD properties of serotonergic agents, an emerging approach recently proposed to assess the pharmacological basis of adverse events (AEs) observed in a worldwide pharmacovigilance (PV) database [10][11][12].
First, the FAERS public dashboard database, the US repository of AEs and medication errors comprising more than 18 million reports gathered worldwide, was queried to retrieve SS reports recorded until December 2019, using "serotonin syndrome" as preferred term, to extract only those records where linezolid was reported as suspect. In order to select potential DDIs, co-administered serotonergic agents reported as suspect were identified according to the lists proposed by Bower et al. [6] and Woytowish et al. [13]. Post-marketing studies using the FAERS are effective for continuous monitoring of old medications [14], including the analysis of potential DDIs [15], with noteworthy performance (i.e., the capacity to discriminate true from positive drug-event associations) [16], especially for side effects with low/rare background incidence and a likely drugattributable component such as SS.
For each serotonergic agent concomitantly reported with linezolid in ≥ 5 cases of SS, two indexes were created: (a) proportion of SS reports, as compared with the overall number of reports, and (b) mean number of DDIs, based on concomitant serotonergic agents recorded in SS reports, being 1 the value corresponding to the single interaction between linezolid and the drug of interest. A scatterplot showing the relationship between the two indexes for each serotonergic agent was The ID consultant "dilemma": START LINEZOLID? created. Threshold values of ≥ 0.1% and ≤ 1.5 were respectively selected for proportion of SS reports (according to 0.09% incidence of SS in a retrospective study including patients exposed to serotonergic medications [17]) and mean number of DDIs (based on the estimation that in almost half of records, no other concomitant serotonergic agent was reported, except for drug of interest), thus identifying three different "SS reporting zones": & "red (i.e., alert)-zone," including medications with high proportion of SS reports and low mean number of DDIs & "yellow-zone," including medications with high proportion of SS reports coupled with high mean number of DDIs or low proportion of SS reports coupled with low mean number of DDIs & "green-zone," including medications with low proportion of SS reports and high mean number of DDIs In the second step, physiochemical/PK (peak concentration (C max ), area under plasma concentration curve (AUC), volume of distribution (V D ), and lipophilicity (LogP)) and PD (values of Ki and IC 50 (concentration corresponding to 50% inhibition of activity in vitro) for SERT and 5-HT 2A ) data of each selected serotonergic agents were extracted from MEDLINE/PubMed and DrugBank databases [18] (searches performed on April 22, 2020). For PK parameters, the highest mean value reported for each medication when administered at therapeutic dosage was selected. The following PK/PD indexes were created: C max /K i SERT, C max /IC 50 SERT, C max /K i 5-HT 2A , AUC/K i SERT, AUC/IC 50 SERT, AUC/K i 5-HT 2A , V D /K i SERT, V D /IC 50 SERT, V D /K i 5-HT 2A , LogP/K i SERT, LogP/IC 50 SERT, and LogP/K i 5-HT 2A .
Correlation between PK/PD indexes and mean number of DDIs (PV index) was assessed by calculating Spearman's rank correlation coefficient, considering significant a p value < 0.05.

Results
Overall, 11,429 reports of SS were found, 669 (5.9%) of which mentioning linezolid as suspect agent. Mean age was 55 ± 26.3 years, with no gender preponderance; 99.1% of reports were serious (i.e., resulting in death, hospitalization, or life-threatening event), and death was reported in 41 cases (18 with citalopram co-administration).
A summary of PK/PD properties of serotonergic medications is shown in Supplementary Table 3. Among antidepressants, venlafaxine and sertraline showed respectively the highest V D and the lowest IC 50 for SERT, while methadone exhibited both the highest V D and LogP, as well as the lowest IC 50 for SERT among opioids.

Discussion
Our real-world study, by investigating the worldwide reporting of SS with linezolid and concomitant serotonergic drugs and analyzing their potential contributing role in the light of relevant PK/PD properties, provides some clues towards safer prescribing. Although previously proposed [10][11][12], to the best of our knowledge, this is the first study in which this PV-PK/PD approach combined multiple sources of data and also accounted for the contributing role of DDIs.
Citalopram, escitalopram, and methadone emerged as "redlight drugs" according to PV indexes, given the high proportion of reports coupled with their potential to precipitate SS when mainly administered with linezolid. Notably, these findings are closely correlated with PK/PD features, thus confirming the role of SERT inhibition as key molecular pathway involved in SS [5,11] and suggesting that agents with high binding affinity (low Ki or IC 50 ) for this transporter and large V D or LogP (high lipophilicity to effectively cross the blood-brain barrier) are more likely to precipitate the disease [19]. Consequently, it is expected that agents showing a low mean number of DDIs should exhibit a large V D /Ki SERT ratio (high CNS penetration and high potency), as found in our analysis. This is the case of citalopram, escitalopram, and methadone, which, by virtue of their PK/PD properties, are more likely to cause SS when used individually with linezolid, in contrast with other serotonergic agents, which were reported in SS cases only under polypharmacy regimen [20][21][22]. Our results are in line with those of recently proposed systematic bioinformatics approach [23], which also indicated SERT inhibition by methadone as a key mechanistic link to SS. Conversely, our findings showed non-significant correlation between C max , AUC, and PV indices, as well as 5-HT 2A seams to play a minor role in comparison to SERT in the occurrence of SS. Consequently, our analysis of the different proposed PK/PD indexes likely suggested that penetration in the CNS of a specific serotonergic agent, coupled with its potency on SERT inhibition, could be more important than plasma exposure in precipitating SS when co-administered with linezolid.
Additionally, different implications could be suggested for the two "yellow-zones." Particularly, the one including agents with low proportion of SS reports and low mean number of DDIs, namely, paroxetine, amitriptyline, and mirtazapine, could be considered as a zone exhibiting an intermediate reporting risk between the other "yellow-zone" (including agents with high proportion of SS reports and high mean number of DDIs) and the "red-zone", due to the ability of these serotonergic medications in precipitating SS mainly when co-administered with linezolid in absence of other serotonergic drugs.
The regulatory recommendation of a reasonable washout period when starting linezolid after serotonergic agents poses a clinical challenge when prompt treatment of serious infections is needed. Given its potent activity against multidrug-resistant Gram-positive pathogens and excellent tissue penetration [24], linezolid represents the first-line agent for management of several deep-seated infections (i.e., pneumonia, meningitis, osteomyelitis) for which potential alternatives (namely, daptomycin, ceftaroline, ceftobiprole) may exhibit PK disadvantages [25,26].
We put forward four main determinants to be considered in tailored risk-benefit assessment depending on clinical scenario: (1) color zone of the different serotonergic medications according to PV indexes and PK/PD properties, (2) severity of infection, (3) availability of effective therapeutic alternatives, and (4) risk of worsening of underlying mental disorder.
In the setting of an intubated and deeply sedated critical patient requiring treatment of a life-threatening infection, antidepressant therapy could be abruptly withdrawn, and sedation with extemporaneous or long-term administration of fentanyl or   Spearman's rank correlation coefficient was calculated, and a p value < 0.05 was considered significant N number of serotonergic medications for which specific PK/PD parameters are retrieved, PK pharmacokinetic, PD pharmacodynamic, PV pharmacovigilance, C max peak concentration, AUC area under plasma concentration curve, V D volume of distribution, Ki binding affinity, IC 50 concentration corresponding to 50% inhibition of activity in vitro, SERT serotonin reuptake transporter, DDIs drug-drug interactions, SS serotonin syndrome meperidine (yellow-zone medications) could be replaced with morphine. Conversely, discontinuing citalopram/escitalopram in a patient affected by severe mood disorder requiring prolonged treatment with linezolid for a spondylodiscitis could be challenging. In this scenario, exacerbation of mental illness and suicide risk could exceed the likelihood of SS occurrence, thus requiring close monitoring while maintaining both treatments. Similarly, linezolid should be promptly initiated under strict monitoring in treatment of severe pneumonia despite concomitant methadone replacement therapy, according to severity of infection and the risk for opiate withdrawal syndrome in case of methadone discontinuation.
We acknowledge the limitations of this study, related to both pharmacovigilance analyses, including FAERS data (e.g., quality of reports, potential existence of remaining duplicates, reporting biases, lack of exposure data, inability in establishing firm causality between drug exposure and occurrence of SS, absence of data to calculate the time to onset from drug administration in the public FAERS dashboard, and limited verification of events through clinical features [27,28]), and pharmacokinetic-pharmacodynamic assessment. Additionally, absolute incidence of SS cannot be inferred by analyzing FAERS data. Consequently, although we calculated proportion of SS reports for each serotonergic agent, our index was not intended to assess absolute incidence of SS, but it would only provide a relative proportion as a percentage of overall reports for a specific agent. Moreover, SS may be underdiagnosed due to lack of recognition, thus likely increasing the extent of under-reporting, a well-known phenomenon. Furthermore, heterogeneity of PK/PD studies in assessing the different parameters have to be considered. Particularly, the value of V D may be affected by serotonergic agents subjected to high pre-systemic extraction after oral dosage. However, a trend to significant correlation was also found for indexes correlating lipophilicity (LogP) and binding affinity to SERT (Ki or IC 50 ), thus supporting the hypothesis that agents with high lipophilicity and penetration in CNS coupled with high potency on SERT are more likely to precipitate the disease when used individually with linezolid.
In summary, citalopram, escitalopram, and methadone emerged as "red-light medications," according to their PV indices. High SERT inhibitory potency coupled with large V D and high lipophilicity, proxies for excellent CNS penetration, may explain the role of "red-light drugs" in precipitating SS when individually co-administered with linezolid. Proper management of SS, including drug discontinuation and/or strict monitoring, should be tailored on a case-by-case basis.
Acknowledgements Open access funding provided by Alma Mater Studiorum -Università di Bologna within the CRUI-CARE Agreement.
Author's contributions MG and ER made substantial contributions to conception and design of the study. MG made substantial contributions to acquisition and analysis of data. MG, ER, and FDP made substantial contributions in interpretation of data. MG was involved in drafting the manuscript. ER and FDP made substantial contributions in revising the manuscript critically for important intellectual content. All authors approved the final version of the manuscript.
Funding ER and FDP at the University of Bologna are supported by institutional research funds (Ricerca Fondamentale Orientata).

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Ethics approval Not applicable.
Consent to participate Not applicable.

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