In silico screening for clinical efficacy of antiseizure medications: Not all central nervous system drugs are alike
Roa’a Hamed, Amit David Eyal, Erez Berman, Sara Eyal
Abstract
A recent review by Odi et al. provided important information on the physicochemical and biopharmaceutical properties of antiseizure medications (ASMs).1 The authors included marketed ASMs, ASMs in development, and those that reached clinical trials but did not make it to the market. Based on that analysis, a set of properties was suggested as a screening tool to be used in the early developmental stages of small molecule ASMs intended for oral dosing. These properties influence the permeability of drugs across membranes, which is critical for their absorption, elimination, and distribution across the blood–brain barrier (BBB), and their interactions with drug-metabolizing enzymes. The analysis referred to parameters that may be used individually or as components of desirability scoring systems, namely the drug's molecular weight (MW), logP (the logarithm of a drug n-octanol:water partition coefficient), logD (logP at pH = 7.4), topological polar surface area (TPSA; described in more detail below), and pKa (the negative log of the acid dissociation constant). Among the scoring systems were the “Lipinski rule of five” and the central nervous system (CNS) multiparameter optimization (MPO) desirability tool. The former relates to passive permeation across membranes and predicts oral absorption based on logP, MW, and the number of hydrogen bond donors and hydrogen bond acceptors.2, 3 The CNS MPO tool, used for designing CNS drugs, is a function that maps the molecule's logP, logD, MW, TPSA, number of hydrogen bond donors, and pKa to desirability scores (Table 1).4, 5 Here, we argue that scores that have been developed for CNS drugs as a broader group are not optimal for identification of clinically-efficacious small molecule ASM candidates. Moreover, we suggest that these comparisons might also provide some insights into drug disposition within the epileptogenic brain. These notions are based on our evaluation of the characteristics of ASMs approved by the US Food and Drug Administration (FDA) and/or the European Medicines Agency for at least one epilepsy indication (other than tuberous sclerosis, as described below). We compared those drugs to ASM candidates that reached clinical trials (≥phase 2) but did not make it to the market due to efficacy considerations, and to other CNS drugs.6 We chose MW and TPSA, which are key determinants of small molecule distribution across the BBB,7, 8 as examples of parameters that require refinement. This does not exclude the need for further assessment of other parameters. In addition, the presented analysis is based solely on clinical efficacy data, and not on results from preclinical seizure models. One of the structural prerequisites of successful ASMs that was identified by Odi et al.1 is MW < 375 Da. This value is close to the upper limit of the CNS MPO's desirable MW (≤360 Da).5 The CNS MPO tool was based on an analysis of 119 marketed CNS drugs and 129 Pfizer CNS candidates.4 The drug set used to build the tool included 16 marketed ASMs but also compounds such as aripiprazole and buprenorphine whose MWs are considerably larger than 375 Da. Accordingly, the upper MW cutoff value of this tool for CNS drugs is 500 Da, in agreement with the upper boundary of Lipinski's rule of five (Table 1).2, 9 Since the publication of the analysis by Odi et al., ganaxolone (333 Da) was approved by the FDA for the treatment of seizures associated with cyclin-dependent kinase-like 5 deficiency disorder.10 In contrast, padsevonil (433 Da) showed no advantage over placebo in its primary end points, despite its efficacy in experimental models of drug-resistant epilepsy.11 This adds to the failure of the prodrug VX-765 (MW of the active metabolite = 481 Da) to perform better than placebo in treatment-resistant focal epilepsy.12 Bumetanide (364 Da), which was suggested to act through both cerebral and extracerebral targets,10 did not improve seizure control in newborns with hypoxic ischemic encephalopathy13 and showed no effectiveness in children with autism spectrum disorder, leading to early termination of the trials.10 Although the lack of VX-765 and padsevonil efficacy may not relate solely to the MW of these molecules, it is a missed opportunity to demonstrate the success of compounds larger than the currently used ASMs. It further calls for between-species comparisons of ASM permeability into epileptogenic brain tissue. Besides being a predictor of ASM marketability, the 375-Da cutoff value distinguishes clinically-used ASMs from other classes of CNS drugs, particularly psycholeptics (Figure 1, Table 2). The MW of many CNS drugs is higher than 375 yet they are below the 500-Da upper boundary (Figure 1). An outlier is bromocriptine (655 Da), which is both a substrate of uptake transporters at the BBB15 and a highly potent compound that can exert pharmacological effect at nanomolar concentrations.16 The observed difference between ASMs and other CNS drugs might be associated with distinct structural requirements for receptor binding, but is more likely to provide hints into drug disposition within epileptogenic brain tissue, given the multitude of cerebral ASM targets. In the absence of direct comparisons of BBB permeability across neurological and psychiatric disorders, the reasons for the higher MW cutoff value for CNS drugs other than ASMs have yet to be clarified. Overexpression of the BBB efflux transporter P-glycoprotein (P-gp),17 for which drugs with MW > 400 Da are good substrates,18 could have provided an elegant mechanistic explanation for the low MW cutoff in epilepsy. However, some other brain disorders are also characterized by P-gp upregulation.19, 20 Moreover, recent proteomic analyses suggest the P-gp is globally underexpressed and not overexpressed in the epileptogenic human brain.21, 22 Likewise, key tight junction proteins are downregulated in the epileptogenic brain,22, 23 supporting looser rather than tighter barrier than in the healthy brain. Alternatively, seizures might theoretically alter the porosity of endothelial cell membranes, thereby slowing down or preventing drug “jumping” through transitory holes in the membrane.7, 24 The MW limitation for ASMs appears to conflict with the parenchymal distribution of much larger molecules, for example, albumin (>2000 Da) within epileptogenic tissue.17, 23, 25, 26 The paradox may be explained by heterogeneity of the vasculature in epilepsy21; some vessels, for example, those newly formed through vascular remodeling and aberrant angiogenesis27, 28 could be “leaky,” whereas others are restrictive. That is, both albumin and large molecule ASMs may extravasate permeable capillaries. However, neither would uniformly distribute across diseased brain tissue, such that ASM concentrations near neurons distant from the permeable vessels could be subtherapeutic. ASMs additionally differ from other CNS drugs by generally having more surface area over polar atoms, as indicated by the higher TPSA values of the marketed ASM group relative to the nonmarketed ASMs and other CNS drugs (Figure 1B, Table 2). The TPSA measures the polarity of the compound and is an indicator of the molecule's ability to cross membranes by passive diffusion. TPSA larger than 140 Å2 usually restricts a compound from entering the brain,1, 29 and the optimal range on the CNS MPO score is 40–90 Å2.5 Yet simply enhancing the brain penetration of candidate molecules by decreasing their TPSA (or increasing the logP) may be associated with increased risk of in vivo toxicity,5, 30, 31 particularly for basic molecules.31 In addition, higher lipophilicity leads to unnecessary nonspecific partitioning into brain tissue.32, 33 In two recent indirect comparisons of ASM efficacy in patients with focal seizures,34, 35 the relative responder rate was highest for molecules with TPSA value > 90 (Figure 2). The exact nature of the nonspecific drug disposition sites within epileptogenic brain tissue requires further investigation. The MW and the TPSA values of the mTOR inhibitors everolimus and sirolimus, used in the treatment of tuberous sclerosis, are well beyond those of all other CNS drugs: 958.2 Da and 914.2 Da (as compared to a maximum of 654.6 Da for CNS drugs), and 205 Å2 and 195 Å2 (vs. 124 Å2), respectively. In mice, the distribution of everolimus and sirolimus across the BBB is poor (brain:plasma ratio = .0057 for sirolimus and .016 for everolimus).36 If this is also the case in humans, the antiseizure activity of those ASMs in tuberous sclerosis might be explained by extracerebral sites of action, permeable BBB, or both. Protective effects of sirolimus on the BBB were demonstrated in experimental models of temporal lobe epilepsy, suggesting the BBB as a target.37, 38 The activity of these drugs in tuberous sclerosis might also involve targets beyond the BBB; imaging studies with gadolinium-based contrast agents showed enhancement of white matter lesions, subependymal nodules, subependymal giant cell astrocytomas, and some tubers, particularly cerebellar tubers.39-41 Because gadolinium complexes are large and polar molecules (gadolinium–diethylenetriamine pentaacetic acid, MW = 1031 Da, TPSA = 315 Å2; gadopentetate dimeglumine: MW = 548 Da, TPSA = 205 Å2), it is plausible that their distribution might correlate with those of sirolimus and everolimus and might even serve as a biomarker for such treatment. Similar considerations apply to other seizure etiologies associated with considerable BBB disruption, such as neurocysticercosis.42, 43 Sirolimus and everolimus are also being evaluated for epilepsies with etiologies other than tuberous sclerosis, some of which are without overt BBB disruption (e.g., focal cortical dysplasia44 – NCT03198949, NCT0245169645). The outcomes of these trials can add to our understanding of how mTOR inhibitors exert their activity in patients with epilepsy. Currently, the majority of efficacy data for such epilepsies are based on case series and case reports (Table 3). Current screening tools lack specificity and selectivity to effectively prioritize small molecules for development as ASMs. For instance, the CNS MPO score and TPSA would likely exclude cannabidiol and fenfluramine from further development as ASMs and place topiramate at a lower priority (Figure 3). Accordingly, these parameters are not sensitive enough for ASM candidates. On the other hand, the desirable values on commonly used scores for CNS drugs such as Lipinski's rule of five are of low specificity for ASMs and encourage the testing of many compounds that are less likely to become systemically administered ASMs. Exemptions to the aforementioned scores or a refined scoring system could include (1) candidates for seizure prevention in tuberous sclerosis, neurocysticercosis, and other seizure etiologies characterized by overt contrast enhancement; (2) highly potent compounds that exert pharmacological effects even at nanomolar or submicromolar concentrations near their targets, or compounds with little systemic toxicity even at high doses; (3) drugs, prodrugs, and drug formulations that can cross the BBB by using specialized carriers and receptors; and (4) treatments that target extracerebral sites, including the BBB itself and inflammation of the neurovascular unit. Examples include interleukin-1 and interleukin-6 antagonists (anakinra and tocilizumab, respectively) for treating febrile infection-related epilepsy syndrome55 and microbiome-based therapies, taking advantage of microbiotal effects on the brain through secreted metabolites and other communication routes.56-58 In addition, drugs might better distribute into the brain during epileptogenesis if the epileptogenic insult involves BBB disruption, as has been reported for traumatic brain injury.59-61 Antiepileptogenic drugs might also treat processes that do not necessarily require their distribution across the BBB (e.g., inflammation). Accordingly, a wider assortment of drugs might be repurposed as antiepileptogenic versus antiseizure medications, but this would require thorough consideration of the epileptogenic etiology. When interpreting CNS MPO scores for potential drug repurposing in the treatment of epilepsy, researchers should consider that the model was developed for the screening of novel compounds. The score may penalize or reward drugs based on an estimation of properties, such as safety and absorption, distribution, metabolism, and excretion (ADME) profiles, that are already well established for those candidates. For example, the TPSA desirability score hump function used in the CNS MPO model assigns low desirability scores to molecules with low TPSA values such as fenfluramine, due to correlations with increased toxicity, not because of decreased CNS activity. Toxicity and ADME screening are desirable in assessment of novel drug candidates but may exclude repurposing candidates for which the aforementioned screening is redundant. In addition, they may promote others with properties favorable for successful safety and ADME profiles but of reduced relevance in prediction of CNS activity. Importantly, even refinement of efficacy-related parameters alone can minimize animal testing and increase the success of ASM candidates at preclinical stages of development or during repurposing. We acknowledge support of the Israel Research Fund (2054/18; student scholarships). S.E. is the Dame Susan Garth Chair of Cancer Research, David R. Bloom Center for Pharmacy and Dr. Adolf and Klara Brettler Center for Research in Molecular Pharmacology and Therapeutics at the Hebrew University of Jerusalem, Israel. S.E has served as a paid consultant for BioPass, Israel. The other authors have no conflicts of interest. Abbreviations: ASM, antiseizure medication; DEPDC5, pleckstrin domain-containing protein 5; FCD, focal cortical dysplasia; GATOR1, GAP activity toward RAGs 1 complex; LYK5, protein kinase LYK5; MTOR, mechanistic target of rapamycin; NPRL, nitrogen permease regulator-like; PMSE, polyhydramnios, megalencephaly, and symptomatic epilepsy syndrome; STRADα, STE20-related kinase adaptor α. Table S1. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.