Comparative assessment of Ca2+ oscillations in 2- and 3-dimensional hiPSC derived and isolated cortical neuronal networks.

Human induced Pluripotent Stem Cell (hiPSC) derived neural cells offer great potential for modelling neurological diseases and toxicities and have found application in drug discovery and toxicology. As part of the European Innovative Medicines Initiative (IMI2) NeuroDeRisk (Neurotoxicity De-Risking in Preclinical Drug Discovery), we here explore the Ca2+ oscillation responses of 2D and 3D hiPSC derived neuronal networks of mixed Glutamatergic/GABAergic activity with a compound set encompassing both clinically as well as experimentally determined seizurogenic compounds. Both types of networks are scored against Ca2+ responses of a primary mouse cortical neuronal 2D network model serving as an established comparator assay. Parameters of frequency and amplitude of spontaneous global network Ca2+ oscillations and the drug-dependent directional changes to these were assessed, and predictivity of seizurogenicity scored using contingency table analysis. In addition, responses between models were compared between both 2D models as well as between 2D and 3D models. Concordance of parameter responses was best between the hiPSC neurospheroid and the mouse primary cortical neuron model (77% for frequency and 65% for amplitude). Decreases in spontaneous Ca2+ oscillation frequency and amplitude were found to be the most basic shared determinants of risk of seizurogenicity between the mouse and the neurospheroid model based on testing of clinical compounds with documented seizurogenic activity. Increases in spontaneous Ca2+ oscillation frequency were primarily observed with the 2D hIPSC model, though the specificity of this effect to seizurogenic clinical compounds was low (33%), while decreases to spike amplitude in this model were more predictive of seizurogenicity. Overall predictivities of the models were similar, with sensitivity of the assays typically exceeding specificity due to high false positive rates. Higher concordance of the hiPSC 3D model over the 2D model when compared to mouse cortical 2D responses may be the result of both a longer maturation time of the neurospheroid (84-87 days for 3D vs. 22-24 days for 2D maturation) as well as the 3-dimensional nature of network connections established. The simplicity and reproducibility of spontaneous Ca2+ oscillation readouts support further investigation of hiPSC derived neuronal sources and their 2- and 3-dimensional networks for neuropharmacological safety screening.

GABAA receptor-mediated seizure liabilities: a mixed-methods screening approach.

GABAA receptors, members of the pentameric ligand-gated ion channel superfamily, are widely expressed in the central nervous system and mediate a broad range of pharmaco-toxicological effects including bidirectional changes to seizure threshold. Thus, detection of GABAA receptor-mediated seizure liabilities is a big, partly unmet need in early preclinical drug development. This is in part due to the plethora of allosteric binding sites that are present on different subtypes of GABAA receptors and the critical lack of screening methods that detect interactions with any of these sites. To improve in silico screening methods, we assembled an inventory of allosteric binding sites based on structural data. Pharmacophore models representing several of the binding sites were constructed. These models from the NeuroDeRisk IL Profiler were used for in silico screening of a compiled collection of drugs with known GABAA receptor interactions to generate testable hypotheses. Amoxapine was one of the hits identified and subjected to an array of in vitro assays to examine molecular and cellular effects on neuronal excitability and in vivo locomotor pattern changes in zebrafish larvae. An additional level of analysis for our compound collection is provided by pharmacovigilance alerts using FAERS data. Inspired by the Adverse Outcome Pathway framework, we postulate several candidate pathways leading from specific binding sites to acute seizure induction. The whole workflow can be utilized for any compound collection and should inform about GABAA receptor-mediated seizure risks more comprehensively compared to standard displacement screens, as it rests chiefly on functional data.