Biosensors, such as microelectrode arrays that record in vitro neuronal activity, provide powerful platforms for studying neuroactive substances. This study presents a machine learning workflow to analyze drug-induced changes in neuronal biosensor data using complex network measures from graph theory. Microelectrode array recordings of neuronal networks exposed to bicuculline, a GABA \(_A\) receptor antagonist known to induce hypersynchrony, demonstrated the workflow’s ability to detect and characterize pharmacological effects. The workflow integrates network-based features with synchrony, optimizing preprocessing parameters, including spike train bin sizes, segmentation window sizes, and correlation methods. It achieved high classification accuracy (AUC up to 90%) and used Shapley Additive Explanations to interpret feature importance rankings. Significant reductions in network complexity and segregation, hallmarks of epileptiform activity induced by bicuculline, were revealed. While bicuculline’s effects are well established, this framework is designed to be broadly applicable for detecting both strong and subtle network alterations induced by neuroactive compounds. The results demonstrate the potential of this methodology for advancing biosensor applications in neuropharmacology and drug discovery.
