Next Article in Journal
Effect of T-DNA Integration on Growth of Transgenic Populus × euramericana cv. Neva Underlying Field Stands
Previous Article in Journal
Simultaneous Irradiation with UV-A, -B, and -C Lights Promotes Effective Decontamination of Planktonic and Sessile Bacteria: A Pilot Study
Previous Article in Special Issue
A Novel Genetic Variant in MBD5 Associated with Severe Epilepsy and Intellectual Disability: Potential Implications on Neural Primary Cilia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice

by
Dorota Nieoczym
1,*,
Nancy Saana Banono
2,
Katarzyna Stępnik
3,
Agnieszka A. Kaczor
4,
Przemysław Szybkowski
5,6,
Camila Vicencio Esguerra
2,
Wirginia Kukula-Koch
7 and
Kinga Gawel
5,*
1
Department of Animal Physiology and Pharmacology, Institute of Biological Sciences, Faculty of Biology and Biotechnology, Maria Curie-Skłodowska University, Akademicka 19, 20-033 Lublin, Poland
2
Chemical Neuroscience Group, Centre for Molecular Medicine Norway, University of Oslo, Gaustadalleen 21, Forskningsparken, 0349 Oslo, Norway
3
Department of Physical Chemistry, Institute of Chemical Sciences, Faculty of Chemistry, Maria Curie-Skłodowska University, Pl. M. Curie-Skłodowskiej 3/243, 20-031 Lublin, Poland
4
Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modeling Laboratory, Faculty of Pharmacy, Medical University of Lublin, 4A Chodzki St., 20-093 Lublin, Poland
5
Department of Experimental and Clinical Pharmacology, Medical University of Lublin, Jaczewskiego St. 8b, 20-090 Lublin, Poland
6
Clinical Provincial Hospital No. 2 St. Jadwiga Krolowej in Rzeszow, Lwowska St. 60, 35-301 Rzeszow, Poland
7
Department of Pharmacognosy with Medicinal Plants Garden, Medical University of Lublin, Chodźki St. 1, 20-093 Lublin, Poland
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(16), 12949; https://doi.org/10.3390/ijms241612949
Submission received: 23 July 2023 / Revised: 16 August 2023 / Accepted: 16 August 2023 / Published: 18 August 2023
(This article belongs to the Special Issue Epilepsy: Molecular Mechanisms and Therapies)

Abstract

:
The aim of this study is to evaluate the anticonvulsant potential of schisandrin B, a main ingredient of Schisandra chinensis extracts. Schisandrin B showed anticonvulsant activity in the zebrafish larva pentylenetetrazole acute seizure assay but did not alter seizure thresholds in the intravenous pentylenetetrazole test in mice. Schisandrin B crosses the blood–brain barrier, which we confirmed in our in silico and in vivo analyses; however, the low level of its unbound fraction in the mouse brain tissue may explain the observed lack of anticonvulsant activity. Molecular docking revealed that the anticonvulsant activity of the compound in larval zebrafish might have been due to its binding to a benzodiazepine site within the GABAA receptor and/or the inhibition of the glutamate NMDA receptor. Although schisandrin B showed a beneficial anticonvulsant effect, toxicological studies revealed that it caused serious developmental impairment in zebrafish larvae, underscoring its teratogenic properties. Further detailed studies are needed to precisely identify the properties, pharmacological effects, and safety of schisandrin B.

1. Introduction

Epilepsy is a common and serious neurological disease characterized by spontaneous and unpredictable seizures. The reason for these symptoms is the abnormal excessive and synchronous activity of neurons in the brain caused by an imbalance between inhibitory and excitatory neurotransmission, ion channel dysfunctions (mainly K+ and Na+ channels), and/or genetic mutations. Approximately 75 million people worldwide are affected by epilepsy [1] with about 50–60% of them also having other neuropsychiatric ailments, mainly anxiety, depression, and cognitive disorders [2,3,4,5]. The pathophysiological mechanism between epilepsy and other neuropsychiatric problems seems to be bidirectional because some patients might suffer from psychiatric disorders prior to epilepsy development. Moreover, some clinical observations showed that epileptic disorders in particular brain regions might provoke anxiety attacks [2] and, on the other hand, some stressful life events might trigger epilepsy onset; anxiety and stressful incidents are identified as some of the most frequent factors that initiate seizures [3]. The coexistence of both disorders, as well as their close link, might be explained by the involvement of the same neurotransmitter systems, mainly γ-aminobutyric acid (GABA), as well as serotoninergic, dopaminergic, and adrenergic ones [2].
Both epilepsy itself, as well as comorbid neuropsychiatric disorders, must be appropriately treated, and emphasis should not only be placed on the achievement of seizure remission, but also on the treatment of other comorbidities to improve the quality of patients’ lives as much as possible. Despite at least the partial understanding of the mechanisms of epilepsy and the availability of numerous antiseizure medications (ASMs) on the pharmaceutical market, its therapy still encounters many difficulties. Firstly, pharmacotherapy with available anticonvulsant drugs, whether as mono- or polytherapy, does not present the expected therapeutic effects in about 30% of patients and are therefore considered pharmacoresistant [6,7]. Secondly, anticonvulsant drugs provide only symptomatic treatment; they do not prevent and/or suppress epileptogenesis—a process that converts a healthy brain into an epileptic one after the occurrence of the damage factor, i.e., infection, stroke, or trauma [7]. Thirdly, the occurrence of burdensome side effects related to the anticonvulsant medication results in the discontinuation of pharmacotherapy in up to 25% of patients and the inability to achieve fully effective therapeutic doses [8]. Considering the above-described problems, the coexistence of psychological and psychiatric disorders causes additional complications in the therapy of patients with epilepsy. There is still an urgent need to discover new medications that are safe and effective in the treatment of epilepsy as well as comorbid disorders.
One strategy to find new anticonvulsant medications and improve the therapy of epilepsy is to screen natural compounds from plants used in folk medicine [9,10]. In the present study, we investigate the potential anticonvulsant properties of schisandrin B using zebrafish larvae and mouse pentyleneterazole (PTZ) seizure tests. Schisandrin B (Figure 1) is one of the main active compounds of Schisandra chinensis—a traditional Chinese herb used for thousands of years to treat some neurobehavioral disorders, including insomnia and dysphoria [11,12]. The experimental studies confirmed the previously known traditional curative properties of Schisandra chinensis extracts [12,13] and demonstrated other biological activities, i.e., anti-inflammatory, antioxidant, and antitumor properties [14], and possibilities for the medicinal use of this plant and its constituents in the treatment of some neurodegenerative diseases. Extracts of Schisandra chinensis fruits improved cognition and learning abilities, prevented some pathological changes in the brain of mice with an intracerebroventricular injection of amyloid-β1–42 (Aβ1–42) [15,16]. Li et al. [17] revealed that a lower generation of Aβ in schisandrin B-treated mice might result from a decreased transcription and translation of β-secretase 1—an enzyme responsible for the generation of Aβ peptides. Moreover, schisandrin B prevented scopolamine-induced dementia in mice by decreasing acetylcholinoesterase activity and increasing the acetylcholine level in the brain [18] and might also alleviate cognitive dysfunctions occurring during menopause [19].
An experimental study confirmed the antidepressant and anxiolytic properties of Schisandra chinensis [20]. The antidepressant-like activity of schisandrin B alone was associated with the enhancement of serotoninergic neurotransmission (i.e., increased serum serotonin concentration, mRNA expression of serotonin receptor, and brain-derived neurotrophic factor (BDNF)) as well as the attenuation of neuroinflammatory processes (i.e., decreased indoleamine-2,3-dioxygenase expression, microglia activation, and inflammatory cytokines, TNF-α, IL-1β, and IL-6, in the brain) [21]. Schisandrin B alleviated anxiety-like behavior in mice subjected to the acute stressor of forced swimming and this action partly occurred due to its antioxidant activity [22].
The main aim of the present study was to evaluate the influence of schisandrin B on convulsant activity in the PTZ-induced seizure tests in larval zebrafish and mice. An analysis of the local field potential (LFP) recordings from the optic tectum (midbrain), as well as c-fos and bdnf mRNA expressions, were used to verify the anticonvulsant effect previously observed in behavioral seizure tests in zebrafish larvae. This was followed by the evaluation of schisandrin B safety in zebrafish embryos and larvae analyzing their development and behavior in the light–dark transition and thigmotaxis tests. The employed models are recommended for the initial screening of potentially active compounds in the preclinical studies [23,24]. In silico studies were performed before studies involving mouse models to assess the properties of schisandrin B, particularly its ability to cross the blood–brain barrier (BBB). The potential of schisandrin B to pass the BBB was also confirmed based on the determination of changes in its brain concentration at eight different time points after intraperitoneal (ip) administration in a tailored HPLC-ESI-QTOF-MS/MS method. In addition, changes in the plasma level of schisandrin B were also evaluated. The intravenous (iv) PTZ test was employed to assess its anticonvulsant effect in mice. In the final stage of the study, molecular modeling was conducted to find potential targets for schisandrin B that may be involved in its pharmacological effects.

2. Results

2.1. Effect of Schisandrin B in the PTZ Seizure Tests in Zebrafish Larvae

First, we determined maximal tolerated concentration (MTC) of schisandrin B. Four-day-old fish were exposed for 24 h to the following concentrations of the compound, i.e., 0.75, 1.5, 3, 6, and 12 µM. After this time, each fish was scored to determine abnormalities, e.g., heart oedema, yolk sac abnormalities, lack of swim bladder, body curvature, etc., or death [25]. Additionally, we assessed the touch response to evaluate muscle performance and function [26]. All larvae incubated with schisandrin B at 12 µM died. Evident morphological changes were seen in larvae incubated in a solution of 6 µM. Neither obvious morphological changes nor touch response disturbance were noted in larvae incubated in solutions of schisandrin B at concentrations up to 3 µM. Considering this result, the dose of 3 µM was selected as MTC for the PTZ-induced seizure assays.
In our study, the influence of schisandrin B at concentrations of 0.75, 1.5, and 3 µM after 24 h of incubation on the zebrafish larvae convulsant behavior was evaluated. Two-way repeated measure analysis of variance (ANOVA) revealed statistically significant differences in the locomotor activity of larvae in the studied groups at different time points (treatment: F(7, 362) = 166.5, p < 0.001; time: F(14, 1907) = 53.97, p < 0.001; treatment × time interaction: F(98, 5068) = 13.08, p < 0.001; n = 36–48 per group, Figure 2A). The total distance traveled was analyzed with one-way ANOVA and statistically significant differences between the experimental groups were also revealed (F(7, 362) = 166.5, p < 0.001; n = 36–48 per group; Figure 2B). The lowest tested concentration of schisandrin B, i.e., 0.75 µM, did not exert anticonvulsant activity in PTZ-treated larvae but shifted the peak of PTZ-induced hyperlocomotion 4 min later. This meant that the full convulsant activity of PTZ developed 4 min later in schisandrin B-pretreated larvae than in the control (vehicle + PTZ-treated) group (p > 0.05; Figure 2A). The locomotor activity values of fish pretreated with 1.5 or 3 µM schisandrin B were 25.3% and 51.6% lower (both p < 0.001; Figure 2B), respectively, compared to the control (vehicle + PTZ-treated) group. None of the tested schisandrin B concentrations affected the basic locomotion of larvae (p > 0.05, Figure 2).

2.2. The Influence of Schisandrin B on PTZ-Induced LFP Discharges in Zebrafish Larvae

To verify behavioral outcomes from the zebrafish larval PTZ assay, LFPs from the larval optic tectum (midbrain) were recorded in 7-day-old larvae. One-way ANOVA showed statistically significant differences among the experimental groups in relation to both the number of epileptiform-like events (F(3, 38) = 18.48, p < 0.001; n = 5–15/group; Figure 3A) and their mean duration (F(3, 38) = 21.20, p < 0.0001, n = 5–15/group; Figure 3B). Schisandrin B significantly reduced the number of epileptiform-like events (p < 0.01; Figure 3A) as well as their mean duration (p < 0.001; Figure 3B) in comparison to the control (vehicle + PTZ-treated) group. Treatment with schisandrin B alone did not influence the brain activity of larval zebrafish (p > 0.05; Figure 3A,B), similar to the behavioral observations. Representative LFP recordings from experimental groups are presented in Figure 3C.

2.3. The Effect of Schisandrin B on c-fos and bdnf mRNA Expressions in PTZ-Treated Zebrafish Larvae

To assess the influence of schisandrin B on c-fos and bdnf mRNA expressions, zebrafish larvae were pretreated with schisandrin B (3 µM, 24 h) and incubated with PTZ (20 mM, 90 min); then, the larvae were pooled at 10/sample. One-way ANOVA of RT-qPCR data indicated that schisandrin B significantly affected both c-fos (F(3, 28) = 40.88, p < 0.001) and bdnf mRNA expressions (F(3, 28) = 17.30, p < 0.001) in comparison to the control (vehicle + PTZ-treated) group. PTZ application upregulated the expression of both genes (p < 0.001), while schisandrin B significantly attenuated these changes (p < 0.001 for c-fos and p < 0.01 for bdnf). Moreover, schisandrin B alone did not induce any changes in the mRNA expression of those markers (p > 0.05) compared to the respective control group, i.e., vehicle + vehicle-bathed fish (Figure 4).

2.4. Effect of Schisandrin B on Developing Zebrafish

In the first step of the toxicity study, we counted the number of embryos hatched at 72 and 96 h post-fertilization (hpf) in order to test whether schisandrin B caused any disturbances in zebrafish development. We noted that schisandrin B at concentrations of 1.5 and 3 µM delayed hatching in 72 h old embryos (both p < 0.05 vs. the control group); however, there was no difference between these groups and control littermates at 96 hpf (p > 0.05, Figure 5A). Moreover, 60% and 40% of larval zebrafish exposed to 6 µM of schisandrin B were delayed in hatching at 72 and 96 hpf, respectively (both p < 0.001; Figure 5A).
Schisandrin B in a dose-dependent manner induced severe morphological abnormalities in developing zebrafish after 95 h of incubation. Only fish exposed to 0.375 µM of schisandrin B did not differ from their control counterparts (p > 0.05; Figure 5B–F). For doses of 0.75 and 1.5 µM, 40% of fish did not develop their swim bladders (p < 0.01 vs. the control group). Moreover, a concentration of 3 µM caused yolk sac abnormalities in 40% of fish incubated (p < 0.001). All or almost all zebrafish exposed to the highest doses of schisandrin B, i.e., 3 and 6 µM, experienced severe abnormalities—they did not inflate their swim bladders (p < 0.001) and yolk sac abnormalities (p < 0.001), heart oedema (p < 0.001), and hemorrhaging around the liver (p < 0.001) were observed. Moreover, fish bodies were curved (p < 0.001) and hyperpigmented at 96 hpf. For representative images, see Figure 6.

2.5. Effect of Schisandrin B in the Light–Dark Test in Zebrafish Larvae

The effects of schisandrin B on larval behavior in the light–dark test were represented in over 1 min time bins to reveal the temporal dynamics within each phase (Figure 7A) and the aggregate behavior per phase (Figure 7B). Two-way repeated measures of ANOVA revealed statistically significant differences in the distance traveled by larvae among the different treatment groups and illumination conditions (treatment: F(3, 68) = 1.722, p > 0.05; illumination condition: F(2.206, 150) = 29.58, p < 0.001; treatment × illumination condition: F(9, 204) = 1.547, p > 0.05; n = 18 larvae per group). During the dark phase, there was a tendency for all groups to travel a greater distance compared to the light phase, reaching statistical significance in most instances. The schisandrin B-treated groups, 0.75 (p < 0.01) and 1.5 µM (p < 0.05), consistently showed alternating reductions and increasing activity in the light and dark conditions, respectively. However, the treatment groups did not differ significantly during each illumination condition.
We also observed that the outer-zone preference of larvae was significantly influenced by both illumination conditions and schisandrin B treatment (two-way repeated measures ANOVA, treatment: F(3, 68) = 3.417, p < 0.05; illumination condition: F(2.067, 140.6) = 15.05, p < 0.001; treatment × illumination condition: F(9, 204) = 1.766, p > 0.05; n = 18 per group; Figure 7C). Larvae across all groups spent more time in the outer zone during the light phase than in the dark phase, and thigmotaxis in the light phase was generally higher than in the dark phase. Although there were no significant differences in the zone preference between the groups in most cases, we observed a tendency towards statistical significance in the schisandrin B-treated groups compared to the vehicle-treated group during the second dark phase (vehicle-treated group vs. schisandrin B 1.5 µM-treated group, p < 0.01; vs. schisandrin B 0.75 µM-treated group, p = 0.0506 and vs. schisandrin B 3 µM-treated group, p = 0.0758).

2.6. Analysis of BBB Descriptors

To observe if schisandrin B could effectively cross the BBB, some pharmacokinetic parameters were calculated using ACD/Percepta software (version 2012, Advanced Chemistry Development, Inc., Toronto, ON, Canada). The following parameters were calculated: logBB—the distribution of a substance in the blood–brain area (the BBB penetration descriptor), logPS—the rate of passive diffusion/permeability (the permeability–surface area product), logPS*fu,brain—the brain/plasma equilibration rate, Fu—the fraction unbound in plasma, and Fb—the fraction unbound in the brain. The results are presented in Table 1.
In further studies, other physicochemical parameters, i.e., steric, electronic, and lipophilic ones, were calculated. The following logarithms of partition coefficients, as descriptors of lipophilicity, were calculated: n-octanol/water (logPow) and cyclohexane/water (logPcw). The results are presented in Table 2.
Moreover, according to the linear free-energy relationship methodology [27,28], the excess molar refraction (E) was determined. This parameter combined with the lipophilic descriptor logPow and the hydrogen-bonding parameter ΔlogP (the difference between the logarithms of n-octanol/water and cyclohexane/water partition coefficient values) provided a promising combination to estimate the logBB values [29].
To calculate the logBB value for schisandrin B, previously established in [29], the quantitative structure–activity relationship (QSAR) model was also used (Equation (1)). The model was set up to examine the quantitative relationship between the ability of a compound to cross the BBB and the molecule structure.
logBB = −0.114 − 0.098 ΔlogP + 0.278 logPow + 0.218E
In the equation presented above, n = 40, R2 CV = 78.25%, R2 pred = 74.02%, and S = 0.436.
Immobilized artificial membrane (IAM) liquid chromatography was also used as the biomimetic system to determine the BBB permeability of schisandrin B. For this purpose, the value of the logarithm of the retention factor extrapolated to pure water (logkw) was determined. The obtained logkw value was used for the calculation of the logBB-IAM value based on the above-mentioned QSAR equation, instead of the logPow parameter.
The obtained logBB values are presented in Figure 8.

2.7. Changes in Schisandrin B Concentration in Mouse Serum and Brain Tissue

The applied chromatographic protocol was tailored to elevate the sensitivity of the mass spectrometer toward the detection of schisandrin B in biological samples. For this purpose, the settings of gas and sheath gas temperatures, gas flow rate, capillary voltage, fragmentation voltage, and collision energies were optimized using the solution of schisandrin B standard (1 mg/mL). As a result, the capillary voltage was set at 3000 V (that was selected from the range of 2800–4000 V), gas and sheath gas temperatures were 275 and 300 °C (from the range of 250–375 °C), gas flow was 12 L/min (from the range of 8–12 l/min), and the fragmentation was 80 V (from the range of 70–130 V). In the analysis, the authors were able to detect the peak of schisandrin B, which was visible in the mass chromatogram as a sodium adduct [M + Na]+ at the same retention time as the reference compound (Figure 9A,B) and with the same m/z parameter. The accurate mass measurement of the molecule and the analysis of its MS/MS spectra provided sufficient information for the identification of this natural product in comparison to the reference compound. Schisandrin B in the tested samples was detected with m/z of 423.1793 Da. The measurement showed a 2.17 ppm deviation from the theoretical mass. The molecule was characterized by the double-bond equivalent value (DBE) of 13.
The quantitative analysis of the compound of interest in the biological samples was performed based on the calibration curve equation that was obtained for the standard of schisandrin B (purity > 95%, Glentham Life Sciences Ltd., Corsham, UK) that was injected in the same way as the tested samples in the range of concentrations from 0.0005 to 0.1 mg/mL.
Fifteen min after ip administration of schisandrin B at a dose of 50 mg/kg, it was detected in the serum at a concentration of 0.025 ± 0.002 µg/mL and in the following time point, i.e., 30 min after administration, its level increased to the highest noted value of 0.073 ± 0.007 µg/mL. At later time points, the concentration of schisandrin B in plasma was gradually reduced; although, even 8 and 12 h after administration, the level of schisandrin B in the blood plasma was quite high and reached 0.017 ± 0.008 µg/mL and 0.012 ± 0.004 µg/mL, respectively. We did not observe any peak of the compound in the plasma 24 h after the injection. Levels of schisandrin B in mouse plasma are presented in Figure 9C.
A detectable concentration of schisandrin B in the brain tissue was found as early as 15 min after its administration and reached 0.095 ± 0.003 µg/g. The level of this compound increased at subsequent time points (i.e., 30 and 60 min after administration) and reached a maximum value of 0.327 ± 0.027 µg/g 2 h after its administration. Schisandrin B level gradually decreased at the next time points (i.e., 4, 8, and 12 h) and it was not detected in the samples collected 24 h after its administration. Changes in schisandrin B concentration in the brain tissue are presented in Figure 9D.

2.8. Effect of Schisandrin B on Seizure Thresholds in the iv PTZ Seizure Test in Mice

Based on the data concerning changes in schisandrin B levels in the brain tissue, we decided to evaluate its anticonvulsant effect (previously noted in the zebrafish larvae PTZ assay) in the timed PTZ test in mice at 4 different time points, i.e., 30, 60, 120, and 240 min, after its ip administration at a dose of 50 mg/kg. There were no statistically significant changes in the seizure thresholds for the myoclonic twitches (one-way ANOVA: F(4, 50) = 1.026; p > 0.05), generalized clonic seizures (one-way ANOVA: F(4, 50) = 2.16, p > 0.05), as well as tonic convulsions (one-way ANOVA: F(4, 50) = 2.005, p > 0.05) in mice (Figure 10A–C).

2.9. Molecular Docking

In order to rationalize, at the molecular level, the observed in vivo anticonvulsant activity of schisandrin B in zebrafish larvae, molecular modeling was performed. As PASS 2022 and Pharma Expert 2022 (geneXplain, Wolfenbüttel, Germany) software could not able indicate the possible molecular targets for the observed activity of schisandrin B (which confirmed the originality of its structure among potential ASMs), the respective receptors were selected based on our knowledge and the available literature. The activation of dopamine D2 or the serotonin 5-HT2A receptor, the blockade of the NMDA ion channel, and binding in the benzodiazepine (BDZ) site of the GABAA receptor were studied as the possible molecular mechanisms of the anticonvulsant activity of schisandrin B.
The results of the molecular docking of schisandrin B to the orthosteric pockets of D2 and 5-HT2A (both are G-protein-coupled receptors—GPCR) receptors are presented in Figure 11. The most striking difference between the positions of schisandrin B in the dopamine D2 receptor (Figure 11A) and serotonin 5-HT2A receptor (Figure 11C) is the shallow position of the ligand in the D2 receptor. Although the molecular docking was performed on the orthosteric site, in the case of the D2 receptor, schisandrin B was not able to fully penetrate the orthosteric pocket and its position in the receptor could be classified as allosteric (Figure 11B,E). As for the 5-HT2A receptor, schisandrin B occupied the classical orthosteric pocket (Figure 11D,F). For both receptors, the docking of schisandrin B resulted in favorable values for the Glide score (−5.341 and −4.775, respectively) and MMGBSA dG bind energy (−50.27 and −33.54 kcal/mol, respectively) (Table 3).
Our analysis revealed that schisandrin B interacted with the NMDA ion channel (Figure 12B,C) by blocking it and acting as an uncompetitive antagonist. One of its oxygen atoms formed a hydrogen bond with Asn614. Other ligand–receptor interactions were mainly of a hydrophobic nature. The Glide docking score and MMGBSA dG bind energy were also favorable for this ligand–receptor complex (−7.125 and −42.63 kcal/mol, respectively, Table 3).
Finally, the anticonvulsant activity of schizandrin B may be rationalized by its interactions with the GABAA receptor (Figure 13A). We showed here that schisandrin B could interact with the BDZ site of the GABAA receptor, similarly to alprazolam. In this case, there hydrophobic-only interactions were observed, in particular Phe77 was involved π–π stacking contact with one of the aromatic rings of the ligand (Figure 13B,C). Glide docking score and MMGBSA dG bind energy were also favorable for this ligand–receptor complex (−6.955 and −49.25 kcal/mol, respectively, Table 3) indicating the possibility of attractive interactions.

3. Discussion

One of the main goals of the present study was to evaluate the anticonvulsant properties of schisandrin B. The main purpose of anticonvulsant drugs is to control convulsions based on the attenuation of excessive excitatory processes in the central nervous system [30]. Natural products can be treated similar to an adjuvant or even, in certain cases, a primary treatment for various kinds of diseases. Some of the plant metabolites could be interesting sources of compounds with anticonvulsant properties, e.g., Melissa parviflora [31], Berberis sibrica [32], Zingiber officinale [25], or Carissa edulis [33].
Schisandrin B is one of the compounds that are the most abundant in Schisandra chinensis. It is a secondary metabolite that belongs to the class of organic compounds named hydrolysable tannins. It has proven, among others, anti-inflammatory [34], antioxidant [22,35], cardioprotective [36], liver-protective [37], as well as neuroprotective [38,39] properties. Wei et al. [40] demonstrated that Schisandra chinensis extract exhibits sedative–hypnotic effects in rats and affects levels of some neurotransmitters (i.e., glutamate, GABA, noradrenaline, dopamine, and serotonin) in the brain. Moreover, the sedative and hypnotic activities of schisandrin B in rodents were previously revealed by Li et al. [41]. This effect results from the modulation of GABAergic and glutamatergic neurotransmission systems, more precisely, from an increased GABA/glutamate ratio both in the brain structures and peripheral blood, as well as the upregulated mRNA expression of α1 and γ2 subunits of the GABAA receptor in the brain tissues [41]. Described mechanisms of action suggest that schisandrin B might also present anticonvulsant activity. We verified this hypothesis in our present study firstly using the PTZ-induced seizure assay in larval zebrafish. Zebrafish have been used for many years in biomedical studies to model various neurological and psychiatric disorders and, at the present time, they are recommended as a respective animal model for the preclinical screening of drugs with central nervous system effects [42,43].
Our results clearly show the anticonvulsant effect of schisandrin B that reduces behavioral (i.e., distance traveled) and electrographic (i.e., number and duration of epileptiform-like events) signs of seizures in PTZ-treated zebrafish larvae. Moreover, the studied compound influenced the mRNA expression of two molecular markers, i.e., c-fos and bdnf. C-Fos is a recognized indirect marker of the functional activity of neurons and numerous studies have proved its involvement in the seizure generation process both in rodent [44,45] and zebrafish [32,46] PTZ-induced seizure/epilepsy models. Similar to c-fos, epileptic insults also upregulate BDNF expression both at the mRNA and protein levels in brain structures [32,44,47]. BDNF influences, among others, the release of neurotransmitters, neuronal morphology, and synaptogenesis [48]. Moreover, BDNF binding to tropomyosin receptor kinase B (TrkB) activates various intracellular pathways and induces some transcription factors, including c-fos. Numerous studies demonstrated various and complex influences of BDNF on GABAergic neurotransmission. BDNF might enhance GABA transporter 1 activity and thus decrease presynaptic GABA concentration. Moreover, the expression of GABAA receptor subunits might also be changed and GABAA functions might be disrupted through an elevated Ca2+ intracellular concentration [49]. Considering that schisandrin B increases GABA levels and the expression of GABAA receptor subunits [41], we can speculate that a decrease in the mRNA expression of bdnf in schisandrin B-treated zebrafish larvae can enhance GABA neurotransmission and thus reduce seizure activity in the PTZ seizure test. Moreover, a reduced bdnf expression might also be responsible, at least partially, for the decrease in c-fos expression in schisandrin B-treated zebrafish larvae.
In the context of anticonvulsant properties, it is crucial to study the ability of schisandrin B to cross the BBB and its distribution to the brain—the site of its pharmacological action. The BBB is a specialized barrier created by endothelial cells of the capillaries and glia in the brain that is essential for maintaining the homeostasis of the brain and ensuring its protection from potentially harmful substances; however, it is also a major challenge for drug delivery to the central nervous system [49]. Lipophilicity is one of the most crucial features influencing the ability of substances to cross biological barriers [50]. It strongly determines the ability of a chemical compound to be distributed in the biosystem. To observe if schisandrin B can effectively cross the BBB, some pharmacokinetic parameters were calculated using computational methods (Table 1).
The scientific reports indicated that the lower the values of the BBB permeability–surface area product (PS) or the fraction unbound in the brain (Fb), the longer the time it takes to reach brain equilibrium [51]. To achieve the desired pharmacological effect, it must be emphasized that drugs in the blood are present in the unbound form but may also be bound to plasma proteins and erythrocytes. The free drug hypothesis postulates that all the distribution processes depend on the unbound drug concentration [52]. In our experiment, it was shown that the values of both the fraction unbound in the brain and the fraction unbound in the plasma were very low (0.01 and 0.11, respectively). Moreover, the brain/plasma equilibration rate was relatively high (logPSFu,brain −3.2). A high rate of penetration results from the high permeability expressed as PS as well as low binding to the brain tissue. However, most contemporary biomimetic and computational methods estimate the ability of a molecule to cross the BBB based on the logBB parameter being the brain/blood partitioning one in a steady state. It is defined as the logarithmic ratio between the concentration of a substance in the brain and its concentration in the blood [53]. As a result of the analysis of all the above-mentioned pharmacokinetic parameters, it can be concluded that schisandrin B is able to cross the BBB.
However, it should be strongly emphasized that the permeability of the BBB is determined by the physicochemical properties of a compound, mainly by its lipophilicity as well as it is influence by the active transport at the BBB [52]. Lipophilic descriptors characterize the transport and resorptive properties of a compound, including the interaction of a molecule with competing binding sites and metabolizing enzymes. Lipophilicity is most often characterized by the substance partition coefficients in various experimental systems, since hydrophobic interactions are responsible for the interaction of a molecule with the molecular target. In our case, lipohilicity was determined by the following logarithms of partition coefficients: n-octanol/water (logPow) and cyclohexane/water (logPcw). Analyzing these parameters, it can be concluded that schisandrin B is a highly lipophilic substance, capable of crossing biological barriers, including the BBB. According to Hansch’s [50,54] approach, the permeability of a substance through biological barriers was also determined by the other physicochemical parameters of a molecule, including steric and electronic ones. Schisandrin B has a molar weight equal to 400.46 g/mol, and it is polarizable. Moreover, to penetrate the BBB and act on receptors in the central nervous system, the topological polar surface area (TPSA) should be less than 90 Å2 [55]. Our in silico studies showed that the TPSA value of schisandrin B was equal to 55.38 Å2. To summarize, the analysis of the physicochemical parameters performed, among others, based on the Hansch approach [50,54] and the Lipinski “rule of five” [56] proved that schisandrin B was capable of crossing the BBB.
The quantitative structure–activity relationship (QSAR) plays an important role in contemporary drug design. The QSAR methodology is based on the theory that biological activity is a sum of various interactions of a compound during its transport through biological barriers and in the reaction with the sites of action (receptor) [57]. In this experiment, QSAR studies were used for the investigation of the dependence between the structure of schisandrin B and its capability to cross the BBB expressed as logBB. Many QSAR models for the estimation of the ability of compounds to cross the BBB were previously established. Most of them were based on both the Hansch and LFER models [29,53,57,58,59,60,61]. To calculate the logBB values for schisandrin B, the QSAR model previously obtained was used (Equation (1)) [29]. This predictive model is based on the lipophilic properties (logPow), excess molar refraction (E) taken from the Abraham model, and the hydrogen-bonding potential, expressed as the difference between the logarithms of n-octanol/water and cyclohexane/water partition coefficient values (ΔlogP). The logBB-QSAR value was calculated based on the values obtained in silico (see Table 2). Similar values of logBB in silico and logBB-QSAR may suggest that the QSAR model is applicable and predictive in assessing the permeation of schisandrin B through the BBB. To confirm this assumption, biomimetic chromatographic studies were performed using the immobilized artificial membrane column (IAM) being a more realistic model than a theoretical one.
The use of IAMs as chromatographic packing materials started the era of biomimetic chromatographic partitioning systems [38]. The permeability of a substance through biological barriers is determined by the membrane partition coefficient, which, in turn, is difficult to measure in vivo [62]. The usefulness of IAM liquid chromatography has been extensively explored since the mid-nineties [63]. The IAM column contains a monolayer of phosphatidylcholine covalently bound to a propylamino-silica core. It is used in high-performance liquid chromatography (HPLC) systems with a physiological buffer as an eluent [64]. A monolayer of phosphatidylcholine immobilized on a silica support is an imitation of the phospholipid bilayer that can mimic drug/biomembrane interactions [65]. The retention on the IAM stationary phase is often related to the specific biological activity of a given compound, including permeation through the BBB, oral absorption, transdermal transport, etc. [63,66].
In this study, IAM chromatography was used to determine the BBB permeability of schisandrin B. For this purpose, the value of the logarithm of the retention factor extrapolated to pure water (logkw) was determined. It was recognized to be an alternative to the logPow lipophilicity descriptor [67]. Therefore, the obtained logkw value was used as a lipophilic descriptor for the calculation of logBB based on the QSAR model (Equation (1)). This value was then compared with the logBB in silico value and logBB-QSAR one. A comparison of these three logBB values calculated with different methods is shown in Figure 8. To date, the literature does not present such computational and biomimetic studies with regard to schisandrin B.
In our study, we decided to confirm the results from the zebrafish larvae PTZ seizure test (i.e., anticonvulsant activity of schisandrin B) and in silico analysis (i.e., the ability of schisandrin B to cross the BBB) in in vivo studies in mice. The increasing amount of data confirms that the BBB in zebrafish shares both structural and functional similarities with the mammalian BBB, and thus zebrafish are considered useful screening tools for the estimation of drug distribution through the BBB [68]. However, the development and maturation of the BBB in the zebrafish larvae occurred between 3 and 10 days post-fertilization (dpf) [69]. Since, in the PTZ seizure assay, 6 dpf zebrafish larvae were used, there was some probability that the distribution of schisandrin B into the brain may not have exactly correlated with that in adult zebrafish, mammals, and humans. In our study, changes in schisandrin B concentration in mouse brain tissues were determined to verify its ability to penetrate the BBB. Moreover, the determination of the time point of maximal concentration in the brain was also utilized to choose appropriate time points to evaluate the time–course effect of schisandrin B in the IV PTZ test in mice. For a more comprehensive view of the pharmacokinetic profile of schisandrin B, changes in its plasma concentration were also determined.
First, our results confirm the ability of schisandrin B to cross the BBB. Moreover, it is worth noting that the plasma concentration of schisandrin B reaches the highest peak much earlier, i.e., 30 min after its administration, in comparison with brain tissue content, which is the highest 120 min after injection. Our data also indicate that schisandrin B is distributed into the brain relatively rapidly (it is recognized in the brain tissue 15 min after administration) and is slowly eliminated (a measurable concentration of schisandrin B is persistent up to 12 h). A similar conclusion was previously presented by Zhu et al. [70]. Moreover, our data on the changes in schisandrin B concentration in the brain tissue are consistent with the results of the study conducted by Wang et al. [71] in which micronized schisandrin B was administered orally at a dose of 20 mg/kg in rats. The highest concentration of schisandrin B in the brain tissue was also observed 120 min after its administration; however, what was interesting was that its level was dependent on the sex of the animals—it was much higher in female than in male rats. Since in our study we used schisandrin B, which was not changed by any kind of formulation and only male mice were used in the PTZ test, there was some probability that using female animals and any kind of formulation that increased the bioavailability of the compound might have changed the distribution of schisandrin B and resulted from the seizure test in mice.
The anticonvulsant effect noted in the zebrafish larvae model was verified in the IV PTZ seizure test in mice. This model is a very sensitive test that allows us to evaluate the effect of the studied compounds on three separate components of the seizure behavior, i.e., myoclonic, generalized clonic, and tonic seizures [72]. Although our results confirm the ability of schisandrin B to penetrate the BBB, we did not note its anticonvulsant effect in the IV PTZ test in mice.
The ability of compounds to cross the BBB is a key condition in the development of new drugs influencing the central nervous system functions; however, their effectiveness is also determined by other factors, i.e., the degree of binding by brain tissue and affinity to the target. The high total brain concentration of the compound does not reflect the real amount that is available to reach the respective target in brain tissues, and thus it does not always result in the high efficacy of the drug [73,74]. Considering the results from our in silico analysis and the seizure tests on mice, we can suppose that the lack of an anticonvulsant effect might have resulted from a very low level of the unbound fraction of schisandrin B in the brain. Although this compound revealed a significant anticonvulsant effect in the zebrafish seizure assay, it did not alleviate convulsant activity in the mice because the level of the unbound fraction may not have been sufficient for an effective influence on the cellular targets, for example, GABAA and/or glutamate receptors (two important targets for anticonvulsant compounds). Since we did not determine the unbound fraction of schisandrin B in the brain, further and more precise studies are necessary. Moreover, although Afrikanova et al. [75] validated the PTZ seizure assay in zebrafish larvae as a useful method that might be employed to screen compounds with potential anticonvulsant properties, they also noted some differences in the efficacy of the individual ASMs against PTZ-induced seizures in rodent and zebrafish larvae. These discrepancies might not only result from space, but also from the differences in the maturity of the brain as well as from the uptake and metabolism of the studied compounds [75].
To identify some possible targets for schisandrin B that were responsible or its anticonvulsant effect in the zebrafish larvae, we conducted an analysis of the molecular docking, taking into consideration the selected receptors potentially relevant for this effect. Dopamine receptor modulation affects susceptibility to seizures [76]. In particular, dopamine D2 receptor agonists are anticonvulsants and decrease pilocarpine and kindled seizures [77,78], while antagonists of this receptor are proconvulsant [79]. Similarly, the activation of the serotonin 5-HT2A receptor may also have an antiseizure effect. The 5-HT2A receptor negatively controls the expression of experimental absence seizures, and selective agonists of this GPCR can lead to the creation of innovative anti-absence drugs [80]. Molecular modeling was performed to investigate if schisandrin B could bind to these GPCRs in an active conformation to exert a possible antiseizure effect.
However, it must be emphasized that schisandrin B does not fulfill the classical pharmacophore model requirements for orthosteric ligands of aminergic GPCRs as it does not possess a protonatable nitrogen atom that can interact with the conserved Asp3.32. [81,82]. Thus, such an interaction was not observed and neither were any other polar interactions observed in the case of these two ligand–receptor complexes. There are, however, many ligands of aminergic GPCRs that either do not possess a protonatable nitrogen atom [83] or do not have a nitrogen atom at all, such as salvinorin A—a potent agonist of the κ opioid receptor [84] (also possessing aspartate in the 3.32 position). The lack of polar interactions in the ligand–receptor complexes does not need to result in weak binding only, as hydrophobic interactions are a promising driving force for biomedical applications [85]. Such interactions were observed in the case of complexes of schisandrin B with D2 and 5-HT2A receptors. The allosteric binding pocket in the D2 receptor is formed i.a. by Phe389 (6.51), His393 (6.55), Tyr408 (7.34), Ile403 (el3), Pro404 (7.30), and Pro405 (7.31). The orthosteric binding pocket in 5-HT2A is constituted, i.e., by Asp155 (3.32), Phe339 (6.51), and Phe340 (6.52). Phe340 is involved in π–π stacking contact with one of the aromatic rings of the ligand. The favorable molecular docking results of schisandrin B to D2 and the 5-HT2A receptor in active conformation indicate the agonism or positive allosteric modulation of these receptors by schisandrin B as the possible molecular mechanism of its antiseizure activity.
Another plausible explanation for the antiseizure activity of schisandrin B is its ability to block the NMDA ion channel (Figure 12A). It is well known that competitive and uncompetitive NMDA antagonists exert an anticonvulsant effect [86,87] as the activation of NMDA receptors induces oxidative stress in neurons and triggers excitotoxicity. Moreover, our analysis revealed that the manner in which schisandrin B interacted with the ion channel in the NMDA receptor was similar to that of phencyclidine—a non-competitive NMDA receptor antagonist with anticonvulsant properties [88].
Finally, the anticonvulsant activity of schisandrin B may also be rationalized by its interactions with the GABAA receptor (Figure 13A). GABA, the main inhibitory neurotransmitter in the cerebral cortex, maintains the inhibitory tone that counterbalances neuronal excitation [89]. When this balance is perturbed, it may result in seizures. The results of our analysis reveal that schisandrin B might act as a positive allosteric modulator of the GABAA receptor.
Since we noted a significant anticonvulsant effect of schisandrin B in the zebrafish larval PTZ seizure test, we used this animal model for testing of schisandrin B toxicity, i.e., its influence on zebrafish development and early life-stage swimming behavior. The zebrafish embryo acute toxicity test (Test No. 236) was employed to assess the safety of schisandrin B after 95 h-long exposure, starting at 1 hpf embryos [90], and the light–dark transition and thigmotaxis test were employed to evaluate the changes in their behavior. A zebrafish embryo acute toxicity study, endorsed by the Organization for Economic Cooperation and Development (OECD), is a guideline for the testing of chemicals [90]. The light–dark and thigmotaxis tests were not only used in toxicological studies [91], but also as valuable tools for the assessment of anxiety-like behavior [92]. An increase in locomotor activity induced by the transition from light to dark, as well as staying or moving close to the borders of the space available (“wall hugging” or thigmotaxis), are a well-validated index of anxiety in different animals, including both adult and larval zebrafish [91,92,93]. Our toxicological studies revealed that schisandrin B provoked severe developmental changes in the zebrafish larvae; however, it did not affect their behavior in the light–dark and thigmotaxis tests on mice. These results suggest that the studied compound might have teratogenic activity and might cause developmental impairments during the prenatal stage of development. However, concentrations of schisandrin B that exerted an anticonvulsant effect in the PTZ seizure zebrafish larvae assay did not cause any damage in the later stages of zebrafish development. The results obtained suggest that schisandrin B cannot be used in the prenatal period; however, it is quite safe for organisms in later periods of development.

4. Materials and Methods

4.1. Chemicals

The analytical standard of schisandrin B was purchased from Glentham Life Sciences Ltd. (Corsham, UK). Acetonitrile and dimethyl sulfoxide (DMSO) were purchased from Merck (Darmstadt, Germany; p.a.), citric acid and disodium hydrogen phosphate (Na2HPO4) from Sigma Aldrich (St. Louis, MO, USA; p.a.), and PTZ from Sanofi-Aventis (Paris, France). Distilled water was obtained from the Direct-Q3 UV apparatus (Millipore, Warsaw, Poland).
For the zebrafish experiments, schisandrin B was dissolved in DMSO to achieve a final concentration of 10 mM in stock and then diluted in zebrafish E3 medium. The final concentration of DMSO in schisandrin B solutions did not exceed 0.06%. The equivalent dose of DMSO was diluted in a E3 medium and this group served as the control group. PTZ was dissolved in E3 medium to achieve a final concentration of 60 mM in stock.
For the mouse studies, schisandrin B was suspended in 1% Tween 80 in normal saline (0.9% NaCl) and injected at a volume of 10 mL/kg. The control group received the vehicle, i.e., 1% Tween 80 in saline, at the same volume.

4.2. Animals and Ethical Approval

For the zebrafish embryo toxicity assay, embryos of the AB strain were obtained from the zebrafish facility at the Experimental Medicine Centre, Medical University of Lublin, Poland. For the remaining experiments, zebrafish embryos of the AB strain were raised at the zebrafish facility at the Centre for Molecular Medicine Norway (NCMM), Oslo, Norway. In both laboratories, the zebrafish embryos were kept under standard housing conditions (28.5 °C; 14 h/10 h light/dark cycle). They were maintained in E3 embryo water, i.e., medium (1.5 mM HEPES, pH 7.6, 17.4 mM NaCl, 0.21 mM KCl, 0.12 mM MgSO4, and 0.18 mM Ca(NO3)2).
For the experiments conducted at NCMM, larvae up to 7 dpf were used. The experiments were approved by the Norwegian Food Safety Authority experimental animal administration’s supervisory and application system (“Forsøksdyrforvatningentilsyns- og søknadssystem”, FOTS-ID 23935). For the experiments conducted in Poland, larvae up to 4 dpf were used. Ethical permission was not required for the experiments using embryos and larval zebrafish up to 120 hpf. A total of 15 µM tricaine solution was used for zebrafish euthanasia. Experiments conducted using zebrafish larvae were repeated three times and the obtained data were pooled together.
Naïve male Swiss mice weighing 23–28 g obtained from a licensed breeder (Laboratory Animals Breeding, Warszawa, Poland) were also used in the present study. The animals were housed in polycarbonate cages under standardized conditions with ambient temperature set at 21–23 °C, relative humidity at 45–55%, a 12/12 light/dark cycle (light on at 6:00 a.m.), and free access to chow pellets and tap water. Prior to the experiments, the animals were exposed to a one-week-long acclimatization period. All experiments were performed at the same time of day (between 8:00 a.m. and 3:00 p.m.) to minimize circadian influences. Control and drug experiments were always performed on the same day to avoid day-to-day variations in convulsive susceptibility. The experimental procedures and protocols were approved by the Local Ethics Committee in Lublin, Poland (license no. 40/2020).
The experiments in all laboratories were performed in compliance with the National Institute of Health Guidelines for the Care and Use of Laboratory Animals, the European Community Council Directive of November 2010 for Care, and Use of Laboratory Animals (Directive 2010/63/EU) guidelines. All efforts were made to minimize the number of animals used, as well as their stress and suffering levels.

4.3. Zebrafish Toxicity Test

To assess the safety of schisandrin B on developing organisms, a zebrafish embryo acute toxicity study was conducted in accordance with the OECD guideline for the testing of chemicals (Test No. 236), as described in detail in our previous papers [90,94]. Viable and completely transparent embryos were selected one hour after fertilization and exposed to various concentrations of schisandrin B (0.375, 0.75, 1.5, 3, or 6 µM) for 1–96 hpf, with fresh solutions provided every 24 h. The control group was exposed to 0.06% DMSO. At first, we tested the toxicity of schisandrin B at log concentrations, i.e., 1000, 100, 1, and 0.1 µM. If we found 5–4 dead fishes per group after a 23 h-long incubation period (from 1 to 24 hpf), we chose lower concentrations of schisandrin B until we found a dose that did not cause mortality until 96 hpf. Only the last experiments were replicated a few times using a higher number of zebrafish larvae. At 72 and 96 hpf, the hatching rate was monitored to assess developmental delay. At 96 hpf, observations were performed to detect abnormalities, including pericardial edema, jaw development, yolk sac necrosis, swim bladder development, body axis shape/curvature, hemorrhage, and changes in pigmentation [95]. For imaging purposes, 96 hpf-old larvae were mounted in methylcellulose and photographed (Olympus SZ61, Olympus, Tokyo, Japan).

4.4. Evaluation of MTC of Schisandrin B for PTZ-Induced Seizure Assay

Four-day-old larval zebrafish were pretreated with schisandrin B (12, 6, 3, 1.5, or 0.75 µM) solutions for 24 h. After this time, the larvae were individually checked under the microscope (Olympus SZ61, Tokyo, Japan) to detect signs of toxicity or abnormalities. The morphological traits were scored as previously described [25,26]. In addition, escape response upon a touch of the tail (to determine muscle function) and hypoactivity (to determine sedation) were checked.

4.5. PTZ-Induced Seizure Assay in Larval Zebrafish

The locomotor activity experiments were conducted as described in detail in our previous papers [25,32,96]. Zebrafish larvae at 6 dpf were pretreated with schisandrin B (concentrations of 0.75, 1.5, or 3 µM) for 24 h. The tested doses were determined taking into account the MTC and results obtained in the subsequent trials. The aim was to determine the lowest doses of schisandrin B that had an anticonvulsant effect. Prior to a vehicle (control group) or PTZ application, the larvae were habituated to the behavior tracker for 15 min. Then, the vehicle or PTZ (final concentration of 20 mM) was added to each well of the 48-well plate. Five min later, larval activity was monitored (ZebraBox, Viewpoint, Lyon, France) for 30 min with 2 min time intervals. The distance covered by each larva in millimeters (mm) was measured. The data from three independent experiments were pooled together.

4.6. LFP Recordings in Larval Zebrafish

The LFP recordings were taken as described in our previous papers [25,32,97]. Six-day-old larvae were pretreated with schisandrin B (3 µM, the highest tested dose in the zebrafish larvae PTZ assay) or vehicle for 24 h with 5 min-long follow-up exposure to 20 mM PTZ. Then, the larvae were mounted on a glass slide in a layer of low-melting-point agarose (2%). Subsequently, the glass electrode filled with artificial cerebrospinal fluid (124 mM NaCl, 2 mM KCl, 2 mM MgSO4, 2 mM CaCl2, 1.25 mM KH2PO4, 26 mM NaHCO3, 10 mM glucose) was inserted into the optic tectum (MultiClamp 700B amplifier, Digidata 1550 digitizer, Axon Instruments, Scottsdale, AZ, USA). The LFP recordings for each larva were collected within a 20 min time frame. For the analysis, Clampfit 10.2 software (Molecular Devices Corporation, San Jose, CA, USA) and custom-written R script for Windows were used.

4.7. RT-qPCR Assessment

RT-qPCR analysis was conducted according to our previously published methods [32]. Six-day-old larvae were pretreated with schisandrin B (3 µM, the highest tested dose in the zebrafish larvae PTZ assay) or vehicle for 24 h with 90 min-long follow-up exposure to 20 mM PTZ. After the treatment, the larvae were collected in a pool of n = 10/sample and RNA was extracted using TRIZOL (Invitrogen, Carlsbad, CA, USA). Subsequently, cDNA was synthesized using SuperScript™ IV First-Strand Synthesis System (Invitrogen, Carlsbad, CA, USA) and amplified using PowerUp™ SYBR™ Green Master Mix (Applied Biosystems, Waltham, MA, USA), following the manufacturer’s instructions. Relative enrichment was computed according to the 2−ΔΔt method, while expression levels were normalized against the housekeeping gene glyceraldehyde 3-phosphate dehydrogenase (gapdh) [32]. Primer sequences:
(1) gapdh_f_5′GTGGAGTCTACTGGTGTCTTC3′;
(2) gapdh_r_5′GTGCAGGAGGCATTGCTTACA3′;
(3) c-fos_f_5′CCGATACACTGCAAGCTGA 3′;
(4) c-fos_r_ 5′TGCGGCGAGGATGAACTCTA3′;
(5) bdnf_f_5′AGCTGAAGAGACAACTTGC3′;
(6) bdnf_r_5′CCATAGTAACGAACAGGAT3′.

4.8. Light–Dark Transition and Thigmotaxis Tests

The light–dark transition and thigmotaxis tests were performed, as previously described [97]. Briefly, 7 dpf larvae pretreated in schisandrin B (at the same doses as in the zebrafish larvae PTZ test, i.e., 0.75, 1.5, and 3 µM) for 24 h were placed individually in a 24-well plate. The larvae were habituated in the dark for 15 min. Then, they were tracked 10 min each in alternating light, dark, light and dark illumination conditions. For the XY graph, locomotor activity for 1 min time bins was extracted. Locomotor activity in each illumination condition was reported as the total distance, which was the sum of small and large movements. The thigmotactic behavior of the larvae in each illumination condition was calculated as a percentage of the total distance traveled in the outer zone of the arena. The data were pooled from three independent experiments performed on the same day.

4.9. Evaluation of Schisandrin B Concentration in the Plasma and Brain of Mice

The concentration of schisandrin B was measured in plasma samples and brain homogenates using an HPLC-ESI-QTOF-MS/MS platform (Agilent Technologies, Santa Clara, CA, USA) in the tailored quantitative experiment. First, the serum samples (three samples per group) were mixed with acetonitrile in the ratio of 1:1 v/v, vortexed for 10 min, centrifuged at 12,000 rpm for 30 min, and filtered through a nylon syringe filter with a pore-size diameter of 0.1 µm. The collected brains (three organs per group) were homogenized with 300 µL of acetonitrile:water (50:50 v/v) using a tissue homogenizer (PRO 200, PRO ScientificInc., Oxford, CT, USA). They were centrifuged at 12,000 rpm for 30 min and the supernatant was filtered through nylon syringe filters with a pore diameter of 0.1 µm to the autosampler vials. Directly after the preparation of the samples, they were introduced to the HPLC-MS instrument. For the analysis, the HPLC chromatograph (1200 Series) equipped with a degasser (G1322A), binary pump (G1312C), autosampler (G1329B), and mass detector (G6530B) with electrospray ionization was used (Agilent Technologies, Santa Clara, CA, USA). The chromatographic separation of the matrix was achieved on a Zorbax Eclipse Plus RP-18 (150 mm × 2.1 mm; dp = 3.5 um) column, in a gradient of acetonitrile with 0.1% of formic acid (solvent A) in a 0.1% aqueous solution of the acid (solvent B) in the following program: 0 min: 98% B; 5 min: 80% B; 40 min: 75% B; 45–49 min: 30% B, 50 min: 2%. All solvents were of LC-MS grade. The analysis time was set at 60 min, the flow rate at 0.2 µL/min, the temperature at 20 °C, and the injection volume at 10 µL. The analysis was recorded in both positive and negative ionization modes, in the following settings of the mass spectrometer: capillary voltage—3000 V, gas and sheath gas temperatures—275 and 300 °C, collision energies—10 and 20 V, skimmer voltage—65 V, and nebulizer pressure—35 psig. The spectra were registered in the m/z range of 50–1200 u, and the MS/MS spectra were recorded for every two major signals per scan. Agilent Mass Hunter Workstation (version B.10.00) was used to record and handle the data. The standard of schisandrin B was purchased from Glentham Life Sciences Ltd. (Corsham, UK) and was used for the preparation of a calibration curve equation and the following quantitative analysis of schisandrin B in the biological samples. The samples were injected 6 times (n = 6).

4.10. The Timed PTZ Infusion Test in Mice

The test was conducted according to the method described previously by Nieoczym et al. [96]. At the appropriate time after schisandrin B (50 mg/kg) or vehicle (1% Tween 80) administration, the mice were placed in the restrainer. To perform the infusion into the lateral tail vein, a needle (27G, ¾ in., Sterican®, B. Braun Melsungen AG, Melsungen, Germany) attached by a polyethylene tube (PE20RW, Plastic One Inc. Roanoke, VA, USA) with a plastic syringe containing a 1% solution of PTZ in saline and placed in the syringe pump (model Physio 22, Hugo Sachs Elektronik-Harvard Apparatus GmbH, March-Hugstetten, Germany) was used. The proper placement of the needle in the vein was confirmed by the presence of blood in the tube and, then, adhesive tape was used to secure the needle to the tail. The animal was taken out of the restrainer and the PTZ solution was infused into the vein of a freely moving mouse at a constant rate of 0.2 mL/min. The time intervals from the start of the infusion of the PTZ solution to the appearance of three separate endpoints, i.e., first myoclonic twitch, generalized clonus with loss of righting reflex, and forelimb tonus, were registered. The seizure threshold (in mg of PTZ per kg of body weight) for each endpoint was calculated according to the following formula: PTZ (mg/kg) = (infusion duration (s) × infusion rate (mL/s) × PTZ concentration (mg/mL))/weight (kg).
Experimental groups consisted of 11–12 animals.
The IV PTZ test dose of schisandrin B was the same as in the experiments evaluating the changes in its concentration in the mouse brain tissue and plasma. This was determined based on our preliminary studies.

4.11. In Silico Studies

4.11.1. Chromatographic Biomimetic Studies

The Shimadzu Vp liquid chromatographic system (Shimadzu, Kyoto, Japan) equipped with LC 10AT pump, SPD 10A UV-Vis detector, SCL 10A system controller, CTO-10 AS chromatographic oven, and Rheodyne injector valve with a 20 µL loop was applied in the HPLC measurements.
The solution of schisandrin B was prepared in DMSO at a concentration of 1 mg/mL. The chromatographic conditions were first optimized. The IAM.PC.DD2 stationary phase (100 × 4.6 mm i.d., 10 µm; Regis Chemicals Company, Morton Grove, IL, USA) was employed. Acetonitrile–phosphate-buffered solutions (0.4; 0.5; 0.6; 0.7 v/v; pH 7.4) were used as mobile phases.
The buffer was prepared from the solutions of Na2HPO4 (0.02 mol/dm3) and citric acid (0.01 mol/dm3), and the pH 7.4 value was fixed before the preparation of the mobile phases. The mobile phase flow rate was set at 1 mL/min. The analyte was detected with a UV light at 225 nm. All measurements were obtained at 25 °C. The dead time values were measured from the non-retained compound (citric acid) peaks. The retention time was the average of at least three independent measurements. The values of the peak asymmetry factor were in the acceptable range.
ACD/Percepta software (version 2012, Advanced Chemistry Development, Inc., Toronto, ON, Canada) was used in the in silico studies. The statistical analysis of the obtained results was performed using the Minitab 18 Statistical Software (Minitab Inc., State College, PA, USA).

4.11.2. Molecular Modeling

The 3D structure of schisandrin B was modeled using a LigPrep module [98] of Schrödinger suite of software, v. 2019-4, as previously reported [99,100]. In order to sample the ligand protonation state Epik module [101] of Schrödinger suite of software, v. 2019-4 was applied.
PASS 2022 and Pharma Expert software were used to obtain possible molecular targets of schisandrin B [102]. Receptor structures were downloaded from Protein Data Bank. In the case of the dopamine D2 receptor, the cryo-EM structure of the human receptor in active conformation in the complex with an agonist bromoergocriptine (PDB ID: 6VMS) [103] was used. Similarly, for the serotonin 5-HT2A receptor, the cryo-EM structure of the receptor in active conformation in a complex with an agonist 25-CN-NBOH (PDB ID: 6WHA) [104] was selected for molecular docking. In the case of the NMDA receptor, the cryo-EM structure of the human GluN1/GluN2A NMDA receptor in the glycine/CPP bound state (PDB ID: 7EOQ) [105] was used to further studies. Finally, for the GABAA receptor, the cryo-EM structure of the human full-length heteromeric α1β3γ2L GABAA in a complex with alprazolam (Xanax) and GABA (PDB ID: 6HUO) [106] was applied. The structures of the receptors were preprocessed using the Protein Preparation Wizard of Maestro Release 2019.4 [107], as previously reported [99,100]. Protonation states of the respective receptor residues were assigned. Hydrogen atoms were added, and minimization was performed, restricted to hydrogen atoms only. Yasara Structure v. 20.12.24 [108] tool for loop modeling was used to construct receptor loops, if necessary. Other missing parts of the receptors were built using Prime [109] of Schrödinger suite of software, v. 2019-4.
The Standard precision (SP) method of molecular docking with Glide [110] from Schrödinger release 2019-4 was applied. For GPCRs and GABAA receptors, the grid files were obtained based on co-crystallized ligands (alprazolam in the case of GABAA receptor). For the NMDA receptor, grid file was based on the phencyclidine position in the ion channel taken from the rat structure of NMDA receptor in complex with phencyclidine (PDB ID: 7SAB) [111]. The selected receptor hydroxyl groups in the active sites were made flexible for GPCRs, as previously reported [112,113]. In the case of the NMDA receptor, the hydroxyl groups of the following residues were made flexible, Thr646A, Thr648B, Thr646C, and Thr648D, while for the GABAA receptor: Tyr58C, Tyr160D, Ser205D, and Ser206D. A total of 100 poses were generated in the case of each ligand–receptor complex. Final poses were selected from the highest-scoring poses based on the visual inspection. Maestro Release 2019.4 [107] and PyMol 2.0.4 [114] software were used for the visualization of the molecular modeling results. Prime MM/GBSA [115] of Schrödinger suite of software v. 2019-4 was used to estimate the relative binding affinity of the protein–ligand complex, without taking into account any simulation process. The flexible residue distance was defined for 5 Å.

4.12. Statistical Analysis

One-way or two-way ANOVAs with repeated measures, with Tukey’s or Bonferroni’s post hoc tests, respectively, were used to analyze locomotor activity, thigmotaxis data, LFP recordings, and RT-qPCR from the zebrafish larvae, as well as the results obtained from the IV PTZ test in mice. The results are presented as mean ± SEM, ± SD or individual measurements with mean ± SEM. In the toxicological assay, the chi-squared and Fisher’s exact tests were used for analyses. The criterion for statistical significance was set at p < 0.05. For statistical purposes and figure generation, GraphPad Prism 8.4.3 or 9.3.1 versions (San Diego, CA, USA) were used.

5. Conclusions

Our study revealed the anticonvulsant properties of schisandrin B. However, the anticonvulsant effect was observed only in the larval zebrafish PTZ seizure assay. Although schisandrin B could cross the BBB, which was verified both in in silico and in vivo studies in mice, it did not influence convulsant activity in the timed PTZ infusion test in mice. The lack of anticonvulsant action in the mouse seizure test might have resulted from a very low level of unbound fraction of the compound in brain tissue, which significantly prevented its interaction with the respective molecular targets, i.e., with GABAA and/or NMDA receptors. The toxicological studies reveal that schisandrin B has a detrimental impact on zebrafish embryonic development, which suggests its teratogenic potential. Further detailed studies to clarify the properties, safety, and mechanisms of action of schisandrin B are required.

Author Contributions

Conceptualization, D.N. and K.G.; methodology, D.N. and K.G.; software, D.N. and K.G.; validation, K.G.; formal analysis, D.N. and K.G.; investigation, D.N., N.S.B., K.S., A.A.K., P.S., W.K.-K. and K.G.; resources, D.N and K.G.; data curation, K.G.; writing—original draft preparation, D.N., K.S., A.A.K., W.K.-K.; writing—review and editing, D.N., N.S.B., C.V.E., W.K.-K. and K.G.; visualization, D.N., A.A.K., K.G.; supervision, D.N. and K.G.; project administration, D.N. and K.G.; funding acquisition, D.N. and K.G. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by statutory funds recieved from the Maria Skłodowska-Curie University (MPK 502-11-30-07; ZFIN 00000010) and the Medical University of Lublin (DS 448, K.G.).

Institutional Review Board Statement

The experimental protocols and procedures are in accordance with the Guide for the Care and Use of Laboratory Animals and were approved by the Local Ethics Committee in Lublin, Poland (approval numbers: 40/2020). Experiments conducted at NCMM, University of Oslo, Norway, were approved by the Norwegian Food Safety Authority experimental animal administration’s supervisory and application system (“Forsøksdyrforvatningentilsynsog søknadssystem”, FOTS-ID 23935). All laboratories followed EU Directive 2010/63/EU for animal experiments.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

P.S. is a volunteer at the Department of Experimental and Clinical Pharmacology, Medical University of Lublin, Poland. The authors declare no conflict of interest.

References

  1. Weltha, L.; Reemmer, J.; Boison, D. The role of adenosine in epilepsy. Brain Res. Bull. 2019, 151, 46–54. [Google Scholar] [CrossRef] [PubMed]
  2. Beyenburg, S.; Mitchell, A.J.; Schmidt, D.; Elger, C.E.; Reuber, M. Anxiety in patients with epilepsy: Systematic review and suggestions for clinical management. Epilepsy Behav. 2005, 7, 161–171. [Google Scholar] [CrossRef] [PubMed]
  3. Hingray, C.; McGonigal, A.; Kotwas, I.; Micoulaud-Franchi, J.A. The Relationship Between Epilepsy and Anxiety Disorders. Curr. Psychiatry Rep. 2019, 21, 40. [Google Scholar] [CrossRef]
  4. Rao, G.; Mashkouri, S.; Aum, D.; Marcet, P.; Borlongan, C.V. Contemplating stem cell therapy for epilepsy-induced neuropsychiatric symptoms. Neuropsychiatr. Dis. Treat. 2017, 13, 585–596. [Google Scholar] [CrossRef]
  5. Kanner, A.M.; Helmstaedter, C.; Sadat-Hossieny, Z.; Meador, K. Cognitive disorders in epilepsy I: Clinical experience, real-world evidence and recommendations. Seizure 2020, 83, 216–222. [Google Scholar] [CrossRef] [PubMed]
  6. Fattorusso, A.; Matricardi, S.; Mencaroni, E.; Dell’Isola, G.B.; Di Cara, G.; Striano, P.; Verrotti, A. The Pharmacoresistant Epilepsy: An Overview on Existant and New Emerging Therapies. Front. Neurol. 2021, 12, 674483. [Google Scholar] [CrossRef]
  7. Wahab, A. Difficulties in Treatment and Management of Epilepsy and Challenges in New Drug Development. Pharmaceuticals 2010, 3, 2090–2110. [Google Scholar] [CrossRef]
  8. Perucca, P.; Gilliam, F.G. Adverse effects of antiepileptic drugs. Lancet Neurol. 2012, 11, 792–802. [Google Scholar] [CrossRef]
  9. Lin, C.H.; Hsieh, C.L. Chinese Herbal Medicine for Treating Epilepsy. Front. Neurosci. 2021, 15, 682821. [Google Scholar] [CrossRef]
  10. Aghdash, S.N. Herbal Medicine in the Treatment of Epilepsy. Curr. Drug Targets 2021, 22, 356–367. [Google Scholar] [CrossRef]
  11. Zhang, M.; Xu, L.; Yang, H. Schisandra chinensis Fructus and Its Active Ingredients as Promising Resources for the Treatment of Neurological Diseases. Int. J. Mol. Sci. 2018, 19, 1970. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, M.; Li, N.; Jing, S.; Wang, C.; Sun, J.; Li, H.; Liu, J.; Chen, J. Schisandrin B exerts hypnotic effects in PCPA-treated rats by increasing hypothalamic 5-HT and γ-aminobutyric acid levels. Exp. Ther. Med. 2020, 20, 142. [Google Scholar] [CrossRef]
  13. Zhang, C.; Mao, X.; Zhao, X.; Liu, Z.; Liu, B.; Li, H.; Bi, K.; Jia, Y. Gomisin N isolated from Schisandra chinensis augments pentobarbital-induced sleep behaviors through the modification of the serotonergic and GABAergic system. Fitoterapia 2014, 96, 123–130. [Google Scholar] [CrossRef] [PubMed]
  14. Fan, X.; Elkin, K.; Shi, Y.; Zhang, Z.; Cheng, Y.; Gu, J.; Liang, J.; Wang, C.; Ji, X. Schisandrin B improves cerebral ischemia and reduces reperfusion injury in rats through TLR4/NF-κB signaling pathway inhibition. Neurol. Res. 2020, 42, 693–702. [Google Scholar] [CrossRef] [PubMed]
  15. Jeong, E.J.; Lee, H.K.; Lee, K.Y.; Jeon, B.J.; Kim, D.H.; Park, J.H.; Song, J.H.; Huh, J.; Lee, J.H.; Sung, S.H. The effects of lignan-riched extract of Shisandra chinensis on amyloid-β-induced cognitive impairment and neurotoxicity in the cortex and hippocampus of mouse. J. Ethnopharmacol. 2013, 146, 347–354. [Google Scholar] [CrossRef] [PubMed]
  16. Xu, M.; Yan, T.; Fan, K.; Wang, M.; Qi, Y.; Xiao, F.; Bi, K.; Jia, Y. Polysaccharide of Schisandra Chinensis Fructus ameliorates cognitive decline in a mouse model of Alzheimer’s disease. J. Ethnopharmacol. 2019, 237, 354–365. [Google Scholar] [CrossRef]
  17. Li, Z.; Liu, X.C.; Li, R.; Chang, J. Reduction of Aβ Generation by Schisandrin B through Restraining Beta-Secretase 1 Transcription and Translation. Med. Sci. Monit. 2018, 24, 1219–1224. [Google Scholar] [CrossRef]
  18. Giridharan, V.V.; Thandavarayan, R.A.; Sato, S.; Ko, K.M.; Konishi, T. Prevention of scopolamine-induced memory deficits by schisandrin B, an antioxidant lignan from Schisandra chinensis in mice. Free Radic. Res. 2011, 45, 950–958. [Google Scholar] [CrossRef]
  19. Jiang, Z.J.; Wang, C.Y.; Xie, X.; Yang, J.F.; Huang, J.N.; Cao, Z.P.; Xiao, P.; Li, C.H. Schizandrin ameliorates ovariectomy-induced memory impairment, potentiates neurotransmission and exhibits antioxidant properties. Br. J. Pharmacol. 2015, 172, 2479–2492. [Google Scholar] [CrossRef]
  20. Song, Y.; Shan, B.; Zeng, S.; Zhang, J.; Jin, C.; Liao, Z.; Wang, T.; Zeng, Q.; He, H.; Wei, F.; et al. Raw and wine processed Schisandra chinensis attenuate anxiety like behavior via modulating gut microbiota and lipid metabolism pathway. J. Ethnopharmacol. 2021, 266, 113426. [Google Scholar] [CrossRef]
  21. Yeh, C.-H.; Lin, M.-W.; Lin, H.-C. Anti-neuroinflammaiton and antidepressant effects of Schisandrin B in mice. FASEB J. 2013, 27, 1099.5. [Google Scholar] [CrossRef]
  22. Wu, Y.; Li, Z.C.; Yao, L.Q.; Li, M.; Tang, M. Schisandrin B alleviates acute oxidative stress via modulation of the Nrf2/Keap1-mediated antioxidant pathway. Appl. Physiol. Nutr. Metab. 2019, 44, 1–6. [Google Scholar] [CrossRef] [PubMed]
  23. Norton, W. Measuring Larval Zebrafish Behavior: Locomotion, Thigmotaxis, and Startle. In Zebrafish Protocols for Neurobehavioral Research; Springer: Berlin/Heidelberg, Germany, 2012; pp. 3–20. [Google Scholar] [CrossRef]
  24. Basnet, R.M.; Zizioli, D.; Taweedet, S.; Finazzi, D.; Memo, M. Zebrafish Larvae as a Behavioral Model in Neuropharmacology. Biomedicines 2019, 7, 23. [Google Scholar] [CrossRef] [PubMed]
  25. Gawel, K.; Kukula-Koch, W.; Banono, N.S.; Nieoczym, D.; Targowska-Duda, K.M.; Czernicka, L.; Parada-Turska, J.; Esguerra, C.V. 6-Gingerol, a Major Constituent of Zingiber officinale Rhizoma, Exerts Anticonvulsant Activity in the Pentylenetetrazole-Induced Seizure Model in Larval Zebrafish. Int. J. Mol. Sci. 2021, 22, 7745. [Google Scholar] [CrossRef]
  26. Gawel, K.; Turski, W.A.; van der Ent, W.; Mathai, B.J.; Kirstein-Smardzewska, K.J.; Simonsen, A.; Esguerra, C.V. Phenotypic Characterization of Larval Zebrafish (Danio rerio) with Partial Knockdown of the cacna1a Gene. Mol. Neurobiol. 2020, 57, 1904–1916. [Google Scholar] [CrossRef]
  27. Abraham, M.H.; Takács-Novák, K.; Mitchell, R.C. On the partition of ampholytes: Application to blood-brain distribution. J. Pharm. Sci. 1997, 86, 310–315. [Google Scholar] [CrossRef]
  28. Abraham, M.H. The factors that influence permeation across the blood-brain barrier. Eur. J. Med. Chem. 2004, 39, 235–240. [Google Scholar] [CrossRef]
  29. Stępnik, K.; Kukula-Koch, W. In Silico Studies on Triterpenoid Saponins Permeation through the Blood-Brain Barrier Combined with Postmortem Research on the Brain Tissues of Mice Affected by Astragaloside IV Administration. Int. J. Mol. Sci. 2020, 21, 2534. [Google Scholar] [CrossRef]
  30. Hu, T.; Zhang, J.; Wang, J.; Sha, L.; Xia, Y.; Ortyl, T.C.; Tian, X.; Chen, L. Advances in Epilepsy: Mechanisms, Clinical Trials, and Drug Therapies. J. Med. Chem. 2023, 66, 4434–4467. [Google Scholar] [CrossRef]
  31. Bhat, J.U.; Nizami, Q.; Asiaf, A.; Parray, S.A.; Ahmad, S.T.; Aslam, M.; Khanam, R.; Mujeeb, M.; Umar, S.; Siddiqi, A. Anticonvulsant activity of methanolic and aqueous extracts of Melissa parviflora in experimentally induced Swiss albino mice. EXCLI J. 2012, 11, 1–6. [Google Scholar]
  32. Gawel, K.; Kukula-Koch, W.; Nieoczym, D.; Stepnik, K.; Ent, W.V.; Banono, N.S.; Tarabasz, D.; Turski, W.A.; Esguerra, C.V. The Influence of Palmatine Isolated from Berberis sibiricaRadix on Pentylenetetrazole-Induced Seizures in Zebrafish. Cells 2020, 9, 1233. [Google Scholar] [CrossRef] [PubMed]
  33. Ya’u, J.; Yaro, A.H.; Malami, S.; Musa, M.A.; Abubakar, A.; Yahaya, S.M.; Chindo, B.A.; Anuka, J.A.; Hussaini, I.M. Anticonvulsant activity of aqueous fraction of Carissa edulis root bark. Pharm. Biol. 2015, 53, 1329–1338. [Google Scholar] [CrossRef] [PubMed]
  34. Nasser, M.I.; Zhu, S.; Chen, C.; Zhao, M.; Huang, H.; Zhu, P. A Comprehensive Review on Schisandrin B and Its Biological Properties. Oxid. Med. Cell. Longev. 2020, 2020, 2172740. [Google Scholar] [CrossRef] [PubMed]
  35. Zhu, N.; Cai, C.; Zhou, A.; Zhao, X.; Xiang, Y.; Zeng, C. Schisandrin B Prevents Hind Limb from Ischemia-Reperfusion-Induced Oxidative Stress and Inflammation via MAPK/NF-κB Pathways in Rats. Biomed Res. Int. 2017, 2017, 4237973. [Google Scholar] [CrossRef]
  36. Chen, N.; Ko, M. Schisandrin B-induced glutathione antioxidant response and cardioprotection are mediated by reactive oxidant species production in rat hearts. Biol. Pharm. Bull. 2010, 33, 825–829. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, Q.; Bao, L.; Lv, L.; Xie, F.; Zhou, X.; Zhang, H.; Zhang, G. Schisandrin B regulates macrophage polarization and alleviates liver fibrosis via activation of PPARγ. Ann. Transl. Med. 2021, 9, 1500. [Google Scholar] [CrossRef]
  38. Ba, Q.; Cui, C.; Wen, L.; Feng, S.; Zhou, J.; Yang, K. Schisandrin B shows neuroprotective effect in 6-OHDA-induced Parkinson’s disease via inhibiting the negative modulation of miR-34a on Nrf2 pathway. Biomed. Pharmacother. 2015, 75, 165–172. [Google Scholar] [CrossRef]
  39. Hu, X.-L.; Guo, C.; Hou, J.-Q.; Feng, J.-H.; Zhang, X.-Q.; Xiong, F.; Ye, W.-C.; Wang, H. Stereoisomers of Schisandrin B Are Potent ATP Competitive GSK-3β Inhibitors with Neuroprotective Effects against Alzheimer’s Disease: Stereochemistry and Biological Activity. ACS Chem. Neurosci. 2019, 10, 996–1007. [Google Scholar] [CrossRef]
  40. Wei, B.; Li, Q.; Fan, R.; Su, D.; Chen, X.; Jia, Y.; Bi, K. Determination of monoamine and amino acid neurotransmitters and their metabolites in rat brain samples by UFLC-MS/MS for the study of the sedative-hypnotic effects observed during treatment with S. chinensis. J. Pharm. Biomed. Anal. 2014, 88, 416–422. [Google Scholar] [CrossRef]
  41. Li, N.; Liu, J.; Wang, M.; Yu, Z.; Zhu, K.; Gao, J.; Wang, C.; Sun, J.; Chen, J.; Li, H. Sedative and hypnotic effects of Schisandrin B through increasing GABA/Glu ratio and upregulating the expression of GABAA in mice and rats. Biomed. Pharmacother. 2018, 103, 509–516. [Google Scholar] [CrossRef]
  42. Gawel, K.; Langlois, M.; Martins, T.; van der Ent, W.; Tiraboschi, E.; Jacmin, M.; Crawford, A.D.; Esguerra, C.V. Seizing the moment: Zebrafish epilepsy models. Neurosci. Biobehav. Rev. 2020, 116, 1–20. [Google Scholar] [CrossRef] [PubMed]
  43. Gawel, K.; Banono, N.S.; Michalak, A.; Esguerra, C.V. A critical review of zebrafish schizophrenia models: Time for validation? Neurosci. Biobehav. Rev. 2019, 107, 6–22. [Google Scholar] [CrossRef] [PubMed]
  44. Malhi, S.M.; Jawed, H.; Hanif, F.; Ashraf, N.; Zubair, F.; Siddiqui, B.S.; Begum, S.; Kabir, N.; Simjee, S.U. Modulation of c-Fos and BDNF protein expression in pentylenetetrazole-kindled mice following the treatment with novel antiepileptic compound HHL-6. Biomed Res. Int. 2014, 2014, 876712. [Google Scholar] [CrossRef]
  45. Yang, H.; Shan, W.; Zhu, F.; Yu, T.; Fan, J.; Guo, A.; Li, F.; Yang, X.; Wang, Q. C-Fos mapping and EEG characteristics of multiple mice brain regions in pentylenetetrazol-induced seizure mice model. Neurol. Res. 2019, 41, 749–761. [Google Scholar] [CrossRef] [PubMed]
  46. Baraban, S.C.; Taylor, M.R.; Castro, P.A.; Baier, H. Pentylenetetrazole induced changes in zebrafish behavior, neural activity and c-fos expression. Neuroscience 2005, 131, 759–768. [Google Scholar] [CrossRef]
  47. Skupien-Jaroszek, A.; Walczak, A.; Czaban, I.; Pels, K.K.; Szczepankiewicz, A.A.; Krawczyk, K.; Ruszczycki, B.; Wilczynski, G.M.; Dzwonek, J.; Magalska, A. The interplay of seizures-induced axonal sprouting and transcription-dependent Bdnf repositioning in the model of temporal lobe epilepsy. PLoS ONE 2021, 16, e0239111. [Google Scholar] [CrossRef]
  48. Binder, D.K.; Croll, S.D.; Gall, C.M.; Scharfman, H.E. BDNF and epilepsy: Too much of a good thing? Trends Neurosci. 2001, 24, 47–53. [Google Scholar] [CrossRef]
  49. Kim, J.; Lee, S.; Kang, S.; Kim, S.H.; Kim, J.C.; Yang, M.; Moon, C. Brain-derived neurotropic factor and GABAergic transmission in neurodegeneration and neuroregeneration. Neural Regen. Res. 2017, 12, 1733–1741. [Google Scholar] [CrossRef]
  50. Hansch, C. Quantitative approach to biochemical structure-activity relationships. Acc. Chem. Res. 1969, 2, 232–239. [Google Scholar] [CrossRef]
  51. Liu, X.; Smith, B.J.; Chen, C.; Callegari, E.; Becker, S.L.; Chen, X.; Cianfrogna, J.; Doran, A.C.; Doran, S.D.; Gibbs, J.P.; et al. Use of a physiologically based pharmacokinetic model to study the time to reach brain equilibrium: An experimental analysis of the role of blood-brain barrier permeability, plasma protein binding, and brain tissue binding. J. Pharmacol. Exp. Ther. 2005, 313, 1254–1262. [Google Scholar] [CrossRef]
  52. Hammarlund-Udenaes, M.; Fridén, M.; Syvänen, S.; Gupta, A. On the rate and extent of drug delivery to the brain. Pharm. Res. 2008, 25, 1737–1750. [Google Scholar] [CrossRef] [PubMed]
  53. Platts, J.A.; Abraham, M.H.; Zhao, Y.H.; Hersey, A.; Ijaz, L.; Butina, D. Correlation and prediction of a large blood-brain distribution data set—An LFER study. Eur. J. Med. Chem. 2001, 36, 719–730. [Google Scholar] [CrossRef] [PubMed]
  54. Hansch, C.; Fujita, T. p-σ-π Analysis. A Method for the Correlation of Biological Activity and Chemical Structure. J. Am. Chem. Soc. 1964, 86, 1616–1626. [Google Scholar] [CrossRef]
  55. Hitchcock, S.A.; Pennington, L.D. Structure-brain exposure relationships. J. Med. Chem. 2006, 49, 7559–7583. [Google Scholar] [CrossRef]
  56. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 1997, 23, 3–25. [Google Scholar] [CrossRef]
  57. Gozalbes, R.; Pineda-Lucena, A. QSAR-based solubility model for drug-like compounds. Bioorganic Med. Chem. 2010, 18, 7078–7084. [Google Scholar] [CrossRef]
  58. Young, R.C.; Mitchell, R.C.; Brown, T.H.; Ganellin, C.R.; Griffiths, R.; Jones, M.; Rana, K.K.; Saunders, D.; Smith, I.R.; Sore, N.E.; et al. Development of a new physicochemical model for brain penetration and its application to the design of centrally acting H2 receptor histamine antagonists. J. Med. Chem. 1988, 31, 656–671. [Google Scholar] [CrossRef]
  59. Rose, K.; Hall, L.H.; Kier, L.B. Modeling Blood-Brain Barrier Partitioning Using the Electrotopological State. J. Chem. Inf. Comput. Sci. 2002, 42, 651–666. [Google Scholar] [CrossRef]
  60. Lombardo, F.; Blake, J.F.; Curatolo, W.J. Computation of brain-blood partitioning of organic solutes via free energy calculations. J. Med. Chem. 1996, 39, 4750–4755. [Google Scholar] [CrossRef]
  61. Kaliszan, R.; Markuszewski, M. Brain/blood distribution described by a combination of partition coefficient and molecular mass. Int. J. Pharm. 1996, 145, 9–16. [Google Scholar] [CrossRef]
  62. Ong, S.; Liu, H.; Qiu, X.; Bhat, G.; Pidgeon, C. Membrane Partition Coefficients Chromatographically Measured Using Immobilized Artificial Membrane Surfaces. Anal. Chem. 1995, 67, 755–762. [Google Scholar] [CrossRef]
  63. Russo, G.; Grumetto, L.; Szucs, R.; Barbato, F.; Lynen, F. Screening therapeutics according to their uptake across the blood-brain barrier: A high throughput method based on immobilized artificial membrane liquid chromatography-diode-array-detection coupled to electrospray-time-of-flight mass spectrometry. Eur. J. Pharm. Biopharm. 2018, 127, 72–84. [Google Scholar] [CrossRef] [PubMed]
  64. Darrouzain, F.; Dallet, P.; Dubost, J.-P.; Ismaili, L.; Pehourcq, F.; Bannwarth, B.; Matoga, M.; Guillaume, Y.C. Molecular lipophilicity determination of a huperzine series by HPLC: Comparison of C18 and IAM stationary phases. J. Pharm. Biomed. Anal. 2006, 41, 228–232. [Google Scholar] [CrossRef]
  65. Ong, S.; Liu, H.; Pidgeon, C. Immobilized-artificial-membrane chromatography: Measurements of membrane partition coefficient and predicting drug membrane permeability. J. Chromatogr. A 1996, 728, 113–128. [Google Scholar] [CrossRef] [PubMed]
  66. Kotecha, J.; Shah, S.; Rathod, I.; Subbaiah, G. Relationship between immobilized artificial membrane chromatographic retention and human oral absorption of structurally diverse drugs. Int. J. Pharm. 2007, 333, 127–135. [Google Scholar] [CrossRef] [PubMed]
  67. Jandera, P.; Fischer, J.; Effenberger, H. Characterisation of retention in micellar high-performance liquid chromatography and in micellar electrokinetic chromatography using lipophilicity and polarity indices. J. Chromatogr. A 1998, 807, 57–70. [Google Scholar] [CrossRef]
  68. Kim, S.S.; Im, S.H.; Yang, J.Y.; Lee, Y.R.; Kim, G.R.; Chae, J.S.; Shin, D.S.; Song, J.S.; Ahn, S.; Lee, B.H.; et al. Zebrafish as a Screening Model for Testing the Permeability of Blood-Brain Barrier to Small Molecules. Zebrafish 2017, 14, 322–330. [Google Scholar] [CrossRef]
  69. Fleming, A.; Diekmann, H.; Goldsmith, P. Functional characterisation of the maturation of the blood-brain barrier in larval zebrafish. PLoS ONE 2013, 8, e77548. [Google Scholar] [CrossRef]
  70. Zhu, H.; Zhang, X.; Guan, J.; Cui, B.; Zhao, L.; Zhao, X. Pharmacokinetics and tissue distribution study of schisandrin B in rats by ultra-fast liquid chromatography with tandem mass spectrometry. J. Pharm. Biomed. Anal. 2013, 78–79, 136–140. [Google Scholar] [CrossRef]
  71. Wang, Z.; You, L.; Cheng, Y.; Hu, K.; Wang, Z.; Cheng, Y.; Yang, J.; Yang, Y.; Wang, G. Investigation of pharmacokinetics, tissue distribution and excretion of schisandrin B in rats by HPLC–MS/MS. Biomed. Chromatogr. 2018, 32, e4069. [Google Scholar] [CrossRef]
  72. Mandhane, S.N.; Aavula, K.; Rajamannar, T. Timed pentylenetetrazol infusion test: A comparative analysis with s.c.PTZ and MES models of anticonvulsant screening in mice. Seizure 2007, 16, 636–644. [Google Scholar] [CrossRef] [PubMed]
  73. Spreafico, M.; Jacobson, M.P. In silico prediction of brain exposure: Drug free fraction, unbound brain to plasma concentration ratio and equilibrium half-life. Curr. Top. Med. Chem. 2013, 13, 813–820. [Google Scholar] [CrossRef] [PubMed]
  74. Mateus, A.; Matsson, P.; Artursson, P. A high-throughput cell-based method to predict the unbound drug fraction in the brain. J. Med. Chem. 2014, 57, 3005–3010. [Google Scholar] [CrossRef]
  75. Afrikanova, T.; Serruys, A.S.; Buenafe, O.E.; Clinckers, R.; Smolders, I.; de Witte, P.A.; Crawford, A.D.; Esguerra, C.V. Validation of the zebrafish pentylenetetrazol seizure model: Locomotor versus electrographic responses to antiepileptic drugs. PLoS ONE 2013, 8, e54166. [Google Scholar] [CrossRef] [PubMed]
  76. Brodovskaya, A.; Kapur, J. Anticonvulsant dopamine type 2 receptor agonist activates inhibitory parvalbumin interneurons. Epilepsia 2021, 62, e147–e152. [Google Scholar] [CrossRef] [PubMed]
  77. Turski, L.; Cavalheiro, E.A.; Bortolotto, Z.A.; Ikonomidou-Turski, C.; Kleinrok, Z.; Turski, W.A. Dopamine-sensitive anticonvulsant site in the rat striatum. J. Neurosci. 1988, 8, 4027–4037. [Google Scholar] [CrossRef]
  78. Wahnschaffe, U.; Löscher, W. Anticonvulsant effects of ipsilateral but not contralateral microinjections of the dopamine D2 agonist LY 171555 into the nucleus accumbens of amygdala-kindled rats. Brain Res. 1991, 553, 181–187. [Google Scholar] [CrossRef]
  79. Turski, W.A.; Cavalheiro, E.A.; Ikonomidou, C.; Bortolotto, Z.A.; Klockgether, T.; Turski, L. Dopamine control of seizure propagation: Intranigral dopamine D1 agonist SKF-38393 enhances susceptibility to seizures. Synapse 1990, 5, 113–119. [Google Scholar] [CrossRef]
  80. Venzi, M.; David, F.; Bellet, J.; Cavaccini, A.; Bombardi, C.; Crunelli, V.; Di Giovanni, G. Role for serotonin2A (5-HT2A) and 2C (5-HT2C) receptors in experimental absence seizures. Neuropharmacology 2016, 108, 292–304. [Google Scholar] [CrossRef]
  81. Kaczor, A.A.; Rutkowska, E.; Bartuzi, D.; Targowska-Duda, K.M.; Matosiuk, D.; Selent, J. Computational methods for studying G protein-coupled receptors (GPCRs). Methods Cell Biol. 2016, 132, 359–399. [Google Scholar] [CrossRef]
  82. Bartuzi, D.; Kaczor, A.A.; Targowska-Duda, K.M.; Matosiuk, D. Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors. Molecules 2017, 22, 340. [Google Scholar] [CrossRef] [PubMed]
  83. Staroń, J.; Bugno, R.; Pietruś, W.; Satała, G.; Mordalski, S.; Warszycki, D.; Hogendorf, A.; Hogendorf, A.S.; Kalinowska-Tłuścik, J.; Lenda, T.; et al. Rationally designed N-phenylsulfonylindoles as a tool for the analysis of the non-basic 5-HT(6)R ligands binding mode. Eur. J. Med. Chem. 2021, 209, 112916. [Google Scholar] [CrossRef] [PubMed]
  84. Cunningham, C.W.; Rothman, R.B.; Prisinzano, T.E. Neuropharmacology of the naturally occurring kappa-opioid hallucinogen salvinorin A. Pharmacol. Rev. 2011, 63, 316–347. [Google Scholar] [CrossRef] [PubMed]
  85. Xiao, F.; Chen, Z.; Wei, Z.; Tian, L. Hydrophobic Interaction: A Promising Driving Force for the Biomedical Applications of Nucleic Acids. Adv. Sci. 2020, 7, 2001048. [Google Scholar] [CrossRef]
  86. Löscher, W.; Rundfeldt, C.; Hönack, D. Low doses of NMDA receptor antagonists synergistically increase the anticonvulsant effect of the AMPA receptor antagonist NBQX in the kindling model of epilepsy. Eur. J. Neurosci. 1993, 5, 1545–1550. [Google Scholar] [CrossRef]
  87. Nasiri-Boroujeni, S.; Rahimi-Madiseh, M.; Lorigooini, Z.; Piroti, K.; Rafieian-Koupaei, M.; Amini-Khoei, H. NMDA Receptor Mediates the Anticonvulsant Effect of Hydroalcoholic Extract of Artemisia persica in PTZ-Induced Seizure in Mice. Evid.-Based Complement. Alternat Med. 2021, 2021, 6422451. [Google Scholar] [CrossRef]
  88. Hayes, B.A.; Balster, R.L. Anticonvulsant properties of phencyclidine-like drugs in mice. Eur. J. Pharmacol. 1985, 117, 121–125. [Google Scholar] [CrossRef]
  89. Treiman, D.M. GABAergic mechanisms in epilepsy. Epilepsia 2001, 42 (Suppl. 3), 8–12. [Google Scholar] [CrossRef]
  90. Nakonieczna, S.; Grabarska, A.; Gawel, K.; Wróblewska-Łuczka, P.; Czerwonka, A.; Stepulak, A.; Kukula-Koch, W. Isoquinoline Alkaloids from Coptis chinensis Franch: Focus on Coptisine as a Potential Therapeutic Candidate against Gastric Cancer Cells. Int. J. Mol. Sci. 2022, 23, 330. [Google Scholar] [CrossRef]
  91. Haigis, A.C.; Ottermanns, R.; Schiwy, A.; Hollert, H.; Legradi, J. Getting more out of the zebrafish light dark transition test. Chemosphere 2022, 295, 133863. [Google Scholar] [CrossRef]
  92. Rock, S.; Rodenburg, F.; Schaaf, M.J.M.; Tudorache, C. Detailed Analysis of Zebrafish Larval Behaviour in the Light Dark Challenge Assay Shows That Diel Hatching Time Determines Individual Variation. Front. Physiol. 2022, 13, 827282. [Google Scholar] [CrossRef]
  93. Schnörr, S.J.; Steenbergen, P.J.; Richardson, M.K.; Champagne, D.L. Measuring thigmotaxis in larval zebrafish. Behav. Brain Res. 2012, 228, 367–374. [Google Scholar] [CrossRef]
  94. Stępnik, K.; Kukula-Koch, W.; Plazinski, W.; Gawel, K.; Gaweł-Bęben, K.; Khurelbat, D.; Boguszewska-Czubara, A. Significance of Astragaloside IV from the Roots of Astragalus mongholicus as an Acetylcholinesterase Inhibitor-From the Computational and Biomimetic Analyses to the In Vitro and In Vivo Studies of Safety. Int. J. Mol. Sci. 2023, 24, 9152. [Google Scholar] [CrossRef] [PubMed]
  95. Di Paola, D.; Iaria, C.; Capparucci, F.; Arangia, A.; Crupi, R.; Cuzzocrea, S.; Spanò, N.; Gugliandolo, E.; Peritore, A.F. Impact of Mycotoxin Contaminations on Aquatic Organisms: Toxic Effect of Aflatoxin B1 and Fumonisin B1 Mixture. Toxins 2022, 14, 518. [Google Scholar] [CrossRef] [PubMed]
  96. Nieoczym, D.; Socała, K.; Gawel, K.; Esguerra, C.V.; Wyska, E.; Wlaź, P. Anticonvulsant Activity of Pterostilbene in Zebrafish and Mouse Acute Seizure Tests. Neurochem. Res. 2019, 44, 1043–1055. [Google Scholar] [CrossRef] [PubMed]
  97. Banono, N.S.; Gawel, K.; De Witte, L.; Esguerra, C.V. Zebrafish Larvae Carrying a Splice Variant Mutation in cacna1d: A New Model for Schizophrenia-Like Behaviours? Mol. Neurobiol. 2021, 58, 877–894. [Google Scholar] [CrossRef] [PubMed]
  98. Schrödinger Release 2019-4: LigPrep; Schrödinger, LLC: New York, NY, USA, 2019.
  99. Sparaco, R.; Kędzierska, E.; Kaczor, A.A.; Bielenica, A.; Magli, E.; Severino, B.; Corvino, A.; Gibuła-Tarłowska, E.; Kotlińska, J.H.; Andreozzi, G.; et al. Synthesis, Docking Studies and Pharmacological Evaluation of Serotoninergic Ligands Containing a 5-Norbornene-2-Carboxamide Nucleus. Molecules 2022, 27, 6492. [Google Scholar] [CrossRef]
  100. Kondej, M.; Wróbel, T.M.; Targowska-Duda, K.M.; Leandro Martínez, A.; Koszła, O.; Stępnicki, P.; Zięba, A.; Paz, A.; Wronikowska-Denysiuk, O.; Loza, M.I.; et al. Multitarget Derivatives of D2AAK1 as Potential Antipsychotics: The Effect of Substitution in the Indole Moiety. ChemMedChem 2022, 17, e202200238. [Google Scholar] [CrossRef]
  101. Schrödinger Release 2019-4: Epik; Schrödinger, LLC: New York, NY, USA, 2019.
  102. Poroikov, V.V.; Filimonov, D.A.; Borodina, Y.V.; Lagunin, A.A.; Kos, A. Robustness of Biological Activity Spectra Predicting by Computer Program PASS for Noncongeneric Sets of Chemical Compounds. J. Chem. Inf. Comput. Sci. 2000, 40, 1349–1355. [Google Scholar] [CrossRef]
  103. Yin, J.; Chen, K.-Y.M.; Clark, M.J.; Hijazi, M.; Kumari, P.; Bai, X.-C.; Sunahara, R.K.; Barth, P.; Rosenbaum, D.M. Structure of a D2 dopamine receptor–G-protein complex in a lipid membrane. Nature 2020, 584, 125–129. [Google Scholar] [CrossRef]
  104. Kim, K.; Che, T.; Panova, O.; DiBerto, J.F.; Lyu, J.; Krumm, B.E.; Wacker, D.; Robertson, M.J.; Seven, A.B.; Nichols, D.E.; et al. Structure of a Hallucinogen-Activated Gq-Coupled 5-HT2A Serotonin Receptor. Cell 2020, 182, 1574–1588.e1519. [Google Scholar] [CrossRef] [PubMed]
  105. Wang, H.; Lv, S.; Stroebel, D.; Zhang, J.; Pan, Y.; Huang, X.; Zhang, X.; Paoletti, P.; Zhu, S. Gating mechanism and a modulatory niche of human GluN1-GluN2A NMDA receptors. Neuron 2021, 109, 2443–2456.e5. [Google Scholar] [CrossRef] [PubMed]
  106. Masiulis, S.; Desai, R.; Uchański, T.; Serna Martin, I.; Laverty, D.; Karia, D.; Malinauskas, T.; Zivanov, J.; Pardon, E.; Kotecha, A.; et al. GABAA receptor signalling mechanisms revealed by structural pharmacology. Nature 2019, 565, 454–459. [Google Scholar] [CrossRef] [PubMed]
  107. Schrödinger Release 2019-4: BioLuminate; Schrödinger, LLC: New York, NY, USA, 2019.
  108. Ozvoldik, K.; Stockner, T.; Rammner, B.; Krieger, E. Assembly of Biomolecular Gigastructures and Visualization with the Vulkan Graphics API. J. Chem. Inf. Model. 2021, 61, 5293–5303. [Google Scholar] [CrossRef]
  109. Schrödinger Release 2019-4: Prime; Schrödinger, LLC: New York, NY, USA, 2019.
  110. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide:  Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef]
  111. Chou, T.H.; Epstein, M.; Michalski, K.; Fine, E.; Biggin, P.C.; Furukawa, H. Structural insights into binding of therapeutic channel blockers in NMDA receptors. Nat. Struct. Mol. Biol. 2022, 29, 507–518. [Google Scholar] [CrossRef]
  112. Magli, E.; Kędzierska, E.; Kaczor, A.A.; Severino, B.; Corvino, A.; Perissutti, E.; Frecentese, F.; Saccone, I.; Massarelli, P.; Gibuła-Tarłowska, E.; et al. Synthesis, docking studies, and pharmacological evaluation of 5HT2C ligands containing the N′-cyanoisonicotinamidine or N′-cyanopicolinamidine nucleus. Arch. Pharm. 2019, 352, e1800373. [Google Scholar] [CrossRef]
  113. Magli, E.; Kędzierska, E.; Kaczor, A.A.; Bielenica, A.; Severino, B.; Gibuła-Tarłowska, E.; Kotlińska, J.H.; Corvino, A.; Sparaco, R.; Esposito, G.; et al. Synthesis, docking studies, and pharmacological evaluation of 2-hydroxypropyl-4-arylpiperazine derivatives as serotoninergic ligands. Arch. Pharm. 2021, 354, e2000414. [Google Scholar] [CrossRef]
  114. Schrödinger, LLC. The PyMOL Molecular Graphics System, version 2.0; Schrödinger, LLC: New York, NY, USA, 2019. [Google Scholar]
  115. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert. Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef]
Figure 1. The chemical structure of schisandrin B.
Figure 1. The chemical structure of schisandrin B.
Ijms 24 12949 g001
Figure 2. The effect of schisandrin B (0.75–3 µM) on PTZ-induced seizure-like activity in zebrafish larvae. The data are expressed as distance traveled (in mm) by larvae in 2 min-long time bins (A); total distance covered by larvae during 30 min of the assay (B). Data are analyzed using two-way or one-way ANOVAs with Bonferroni’s or Tukey’s post hoc tests, respectively. Data are presented as mean ± standard error of the mean (SEM). Group schisandrin B 3 µM + vehicle consisted of 36 zebrafish larvae, group schisandrin B 0.75 µM + PTZ consisted of 46 larvae, and the remaining experimental groups had 48 zebrafish larvae. *** p < 0.001 vs. vehicle + vehicle-treated group; ^^^ p < 0.001 vs. vehicle + PTZ-treated group. PTZ, pentylenetetrazole.
Figure 2. The effect of schisandrin B (0.75–3 µM) on PTZ-induced seizure-like activity in zebrafish larvae. The data are expressed as distance traveled (in mm) by larvae in 2 min-long time bins (A); total distance covered by larvae during 30 min of the assay (B). Data are analyzed using two-way or one-way ANOVAs with Bonferroni’s or Tukey’s post hoc tests, respectively. Data are presented as mean ± standard error of the mean (SEM). Group schisandrin B 3 µM + vehicle consisted of 36 zebrafish larvae, group schisandrin B 0.75 µM + PTZ consisted of 46 larvae, and the remaining experimental groups had 48 zebrafish larvae. *** p < 0.001 vs. vehicle + vehicle-treated group; ^^^ p < 0.001 vs. vehicle + PTZ-treated group. PTZ, pentylenetetrazole.
Ijms 24 12949 g002
Figure 3. The effect of schisandrin B on epileptiform-like activity in the optic tectum of larval zebrafish treated with PTZ. Six-day-old fish were incubated with schisandrin B (3 µM) for 24 h and subsequently exposed to PTZ (20 mM). After a 5 min delay, LFP recordings from the larval optic tectum were obtained. Data are shown as: number of events (nr/20 min) (A); mean duration of events (s/20 min) (B); and representative LFP recordings (C). Data were analyzed using one-way ANOVA followed by Tukey’s post hoc test. Data are depicted as individual measurements (as symbols) and mean values. The experimental groups consisted of the following number of zebrafish larvae: vehicle + vehicle − 9, vehicle + PTZ − 15, schisandrin B + vehicle − 5, and schisandrin B + PTZ − 13. *** p < 0.001 vs. vehicle + vehicle; ^^^ p < 0.001, ^^ p < 0.01 vs. vehicle + PTZ. PTZ, pentylenetetrazole.
Figure 3. The effect of schisandrin B on epileptiform-like activity in the optic tectum of larval zebrafish treated with PTZ. Six-day-old fish were incubated with schisandrin B (3 µM) for 24 h and subsequently exposed to PTZ (20 mM). After a 5 min delay, LFP recordings from the larval optic tectum were obtained. Data are shown as: number of events (nr/20 min) (A); mean duration of events (s/20 min) (B); and representative LFP recordings (C). Data were analyzed using one-way ANOVA followed by Tukey’s post hoc test. Data are depicted as individual measurements (as symbols) and mean values. The experimental groups consisted of the following number of zebrafish larvae: vehicle + vehicle − 9, vehicle + PTZ − 15, schisandrin B + vehicle − 5, and schisandrin B + PTZ − 13. *** p < 0.001 vs. vehicle + vehicle; ^^^ p < 0.001, ^^ p < 0.01 vs. vehicle + PTZ. PTZ, pentylenetetrazole.
Ijms 24 12949 g003
Figure 4. The effect of schisandrin B on c-fos (A) and bdnf (B) mRNA expressions in PTZ-bathed larval zebrafish. Six-day-old fish were incubated with schisandrin B (3 µM) for 24 h and subsequently exposed to PTZ (20 mM). After a 90 min delay, fish were collected and pooled (10/sample). Data were analyzed using one-way ANOVA with the Tukey’s post hoc test. Data are depicted as individual measurements (as symbols) and mean + SEM (n = 8/group). *** p < 0.001 vs. vehicle + vehicle; ^^^ p < 0.001, ^^ p < 0.01 vs. vehicle + PTZ. PTZ − pentylenetetrazole.
Figure 4. The effect of schisandrin B on c-fos (A) and bdnf (B) mRNA expressions in PTZ-bathed larval zebrafish. Six-day-old fish were incubated with schisandrin B (3 µM) for 24 h and subsequently exposed to PTZ (20 mM). After a 90 min delay, fish were collected and pooled (10/sample). Data were analyzed using one-way ANOVA with the Tukey’s post hoc test. Data are depicted as individual measurements (as symbols) and mean + SEM (n = 8/group). *** p < 0.001 vs. vehicle + vehicle; ^^^ p < 0.001, ^^ p < 0.01 vs. vehicle + PTZ. PTZ − pentylenetetrazole.
Ijms 24 12949 g004
Figure 5. The effect of schisandrin B on hatching at 72 and 96 hpf (A); swim bladder inflation (B); yolk sac abnormalities (C); heart malformation (D); hemorrhage (E) and body shape of 96 hpf old larvae (F). Zebrafish larvae were exposed to schisandrin B at a concentration range of 0.375–6 µM for 95 h. All scores were performed in triplicates and the results were pooled together. The experimental groups consisted of the following number of zebrafish larvae: vehicle + vehicle − 36, schisandrin B 0.375 µM − 32, schisandrin B 0.75 µM − 33, schisandrin B 1.5 µM − 36, schisandrin B 3 µM − 36, schisandrin B 6 µM − 34. * p < 0.05; ** p < 0.01; *** p < 0.001 vs. control group, Fisher’s exact test.
Figure 5. The effect of schisandrin B on hatching at 72 and 96 hpf (A); swim bladder inflation (B); yolk sac abnormalities (C); heart malformation (D); hemorrhage (E) and body shape of 96 hpf old larvae (F). Zebrafish larvae were exposed to schisandrin B at a concentration range of 0.375–6 µM for 95 h. All scores were performed in triplicates and the results were pooled together. The experimental groups consisted of the following number of zebrafish larvae: vehicle + vehicle − 36, schisandrin B 0.375 µM − 32, schisandrin B 0.75 µM − 33, schisandrin B 1.5 µM − 36, schisandrin B 3 µM − 36, schisandrin B 6 µM − 34. * p < 0.05; ** p < 0.01; *** p < 0.001 vs. control group, Fisher’s exact test.
Ijms 24 12949 g005
Figure 6. Representative images of the zebrafish larvae from the toxicity test. One hpf zebrafish embryo was exposed to different concentrations of schisandrin B (from 0.375 to 6 µM) for 95 h and, afterwards, pictures of the representative zebrafish larvae were taken. Larvae are shown to the same scale (bar = 1 mm).
Figure 6. Representative images of the zebrafish larvae from the toxicity test. One hpf zebrafish embryo was exposed to different concentrations of schisandrin B (from 0.375 to 6 µM) for 95 h and, afterwards, pictures of the representative zebrafish larvae were taken. Larvae are shown to the same scale (bar = 1 mm).
Ijms 24 12949 g006
Figure 7. The effect of schisandrin B (3 µM) on zebrafish larval behavior in the light–dark and thigmotaxis tests. The data are shown as: distance traveled (in mm) by larvae when exposed to alternating light (100%) and dark (0%) conditions in 1 min-long time bins (A), aggregate data over the entire stimulus duration (10 min) (B), the percentage of time that the larvae are located in the outer zone of the swimming area when exposed to alternating light (100%) and dark (0%) conditions (C). Data were analyzed using two-way repeated measures with Tukey’s post hoc test. Data are presented as individual measurements (as symbols) and mean ± SEM. All experimental groups consisted of 18 zebrafish larvae. Letters indicate significant differences within treatment groups: a. comparison from light to dark phase, b. comparison from dark to light phase; ** p < 0.01.
Figure 7. The effect of schisandrin B (3 µM) on zebrafish larval behavior in the light–dark and thigmotaxis tests. The data are shown as: distance traveled (in mm) by larvae when exposed to alternating light (100%) and dark (0%) conditions in 1 min-long time bins (A), aggregate data over the entire stimulus duration (10 min) (B), the percentage of time that the larvae are located in the outer zone of the swimming area when exposed to alternating light (100%) and dark (0%) conditions (C). Data were analyzed using two-way repeated measures with Tukey’s post hoc test. Data are presented as individual measurements (as symbols) and mean ± SEM. All experimental groups consisted of 18 zebrafish larvae. Letters indicate significant differences within treatment groups: a. comparison from light to dark phase, b. comparison from dark to light phase; ** p < 0.01.
Ijms 24 12949 g007
Figure 8. LogBB values calculated for schisandrin B. logBB in silico—the value calculated using ACD/Percepta software, logBB-QSAR—the value calculated based on Equation (1), and logBB-IAM—the value determined from the biomimetic IAM system.
Figure 8. LogBB values calculated for schisandrin B. logBB in silico—the value calculated using ACD/Percepta software, logBB-QSAR—the value calculated based on Equation (1), and logBB-IAM—the value determined from the biomimetic IAM system.
Ijms 24 12949 g008
Figure 9. Mass chromatogram recorded in the positive ionization mode for the reference solution of schisandrin B (concentration 1 mg/mL, injection volume: 10 µL) (A); extracted ion chromatogram of schisandrin B in brain tissue recorded in the positive ionization mode (B), and changes in schisandrin B concentration in the blood plasma (C) and brain tissue (D) at 8 time points, i.e., 15 and 30 min as well as 1, 2, 4, 8, 12, and 24 h, after its ip administration at a dose of 50 mg/kg in mice. Plasma and brain concentrations of schisandrin B are presented as mean ± standard deviation (SD) (n = 6).
Figure 9. Mass chromatogram recorded in the positive ionization mode for the reference solution of schisandrin B (concentration 1 mg/mL, injection volume: 10 µL) (A); extracted ion chromatogram of schisandrin B in brain tissue recorded in the positive ionization mode (B), and changes in schisandrin B concentration in the blood plasma (C) and brain tissue (D) at 8 time points, i.e., 15 and 30 min as well as 1, 2, 4, 8, 12, and 24 h, after its ip administration at a dose of 50 mg/kg in mice. Plasma and brain concentrations of schisandrin B are presented as mean ± standard deviation (SD) (n = 6).
Ijms 24 12949 g009
Figure 10. The time–course effect of schisandrin B (50 mg/kg, ip) on the threshold for myoclonic twitches (A), generalized clonic seizures with the loss of righting reflex (B), and tonic seizures (C) in the iv PTZ test in mice. Control group and groups tested 30 and 240 min after schisandrin B treatment consisted of 11 mice while groups tested 60 and 120 min after schisandrin B treatment had 12 animals. Data were analyzed using one-way ANOVA with the Tukey’s post hoc test. Data are depicted as mean ± SEM. White bars represent the control groups while colored represent schisandrin B-treated groups.
Figure 10. The time–course effect of schisandrin B (50 mg/kg, ip) on the threshold for myoclonic twitches (A), generalized clonic seizures with the loss of righting reflex (B), and tonic seizures (C) in the iv PTZ test in mice. Control group and groups tested 30 and 240 min after schisandrin B treatment consisted of 11 mice while groups tested 60 and 120 min after schisandrin B treatment had 12 animals. Data were analyzed using one-way ANOVA with the Tukey’s post hoc test. Data are depicted as mean ± SEM. White bars represent the control groups while colored represent schisandrin B-treated groups.
Ijms 24 12949 g010
Figure 11. Molecular interactions of schisandrin B with D2 receptor (A,B,E) and 5-HT2A receptor (C,D,F) based on molecular docking. (A,C) Overview of ligand pose in the receptor. Receptor shown in gray in cartoon representation, ligand shown as blue balls; (B,D) 3D details of the binding site. The most important protein residues are shown as sticks with gray carbon atoms. Ligand shown as sticks with blue carbon atoms; (E,F) 2D overview of the binding site.
Figure 11. Molecular interactions of schisandrin B with D2 receptor (A,B,E) and 5-HT2A receptor (C,D,F) based on molecular docking. (A,C) Overview of ligand pose in the receptor. Receptor shown in gray in cartoon representation, ligand shown as blue balls; (B,D) 3D details of the binding site. The most important protein residues are shown as sticks with gray carbon atoms. Ligand shown as sticks with blue carbon atoms; (E,F) 2D overview of the binding site.
Ijms 24 12949 g011
Figure 12. Molecular interactions of schisandrin B with NMDA receptor based on molecular docking. (A) Overview of ligand pose in the receptor. Receptor shown in gray in cartoon representation, ligand shown as blue balls; (B) 3D details of the binding site. The most important protein residues are shown as sticks with gray carbon atoms. Ligand shown as sticks with blue carbon atoms. Polar interactions shown as red dashed lines; (C) 2D overview of the binding site.
Figure 12. Molecular interactions of schisandrin B with NMDA receptor based on molecular docking. (A) Overview of ligand pose in the receptor. Receptor shown in gray in cartoon representation, ligand shown as blue balls; (B) 3D details of the binding site. The most important protein residues are shown as sticks with gray carbon atoms. Ligand shown as sticks with blue carbon atoms. Polar interactions shown as red dashed lines; (C) 2D overview of the binding site.
Ijms 24 12949 g012
Figure 13. Molecular interactions of schisandrin B with GABAA receptor based on molecular docking. (A) Overview of ligand pose in the receptor. Receptor shown in gray in cartoon representation, ligand shown as blue balls; (B) 3D details of the binding site; (C) 2D overview of the binding site.
Figure 13. Molecular interactions of schisandrin B with GABAA receptor based on molecular docking. (A) Overview of ligand pose in the receptor. Receptor shown in gray in cartoon representation, ligand shown as blue balls; (B) 3D details of the binding site; (C) 2D overview of the binding site.
Ijms 24 12949 g013
Table 1. The BBB pharmacokinetic descriptors of schisandrin B calculated in silico (ACD/Percepta).
Table 1. The BBB pharmacokinetic descriptors of schisandrin B calculated in silico (ACD/Percepta).
Log BBLog PSlogPS*fu,brainFuFb
1.15−1.1−3.20.110.01
Table 2. Physicochemical parameters calculated in silico.
Table 2. Physicochemical parameters calculated in silico.
LogPowLogPcwMW [g/mol]TPSA [Å2]PolarizabilityE
4.7913.753400.4655.3843.582.23
Table 3. Glide scores and MMGBSA results for the selected complexes of schisandrin B with the studied receptors.
Table 3. Glide scores and MMGBSA results for the selected complexes of schisandrin B with the studied receptors.
No.ReceptorGlide ScoreMMGBSA dG Bind [kcal/mol]
1D2−5.341−50.27
25-HT2A−4.775−33.54
3NMDA−7.125−42.63
4GABAA−6.955−49.25
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nieoczym, D.; Banono, N.S.; Stępnik, K.; Kaczor, A.A.; Szybkowski, P.; Esguerra, C.V.; Kukula-Koch, W.; Gawel, K. In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice. Int. J. Mol. Sci. 2023, 24, 12949. https://doi.org/10.3390/ijms241612949

AMA Style

Nieoczym D, Banono NS, Stępnik K, Kaczor AA, Szybkowski P, Esguerra CV, Kukula-Koch W, Gawel K. In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice. International Journal of Molecular Sciences. 2023; 24(16):12949. https://doi.org/10.3390/ijms241612949

Chicago/Turabian Style

Nieoczym, Dorota, Nancy Saana Banono, Katarzyna Stępnik, Agnieszka A. Kaczor, Przemysław Szybkowski, Camila Vicencio Esguerra, Wirginia Kukula-Koch, and Kinga Gawel. 2023. "In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice" International Journal of Molecular Sciences 24, no. 16: 12949. https://doi.org/10.3390/ijms241612949

APA Style

Nieoczym, D., Banono, N. S., Stępnik, K., Kaczor, A. A., Szybkowski, P., Esguerra, C. V., Kukula-Koch, W., & Gawel, K. (2023). In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice. International Journal of Molecular Sciences, 24(16), 12949. https://doi.org/10.3390/ijms241612949

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop