Exposure to pesticides across species has been associated with cognitive and motor impairments. As the problem impacts ecosystem stability, food production and public health, it is urgent to develop multifactorial solutions, from regulatory legislation to pharmacological alternatives that ameliorate the impairments. Fipronil, a commonly used insecticide, acts as a GABAA receptor (GABAAR) antagonist and induces motor impairments in vertebrates and invertebrates. Here, we hypothesized that kaempferol, a secondary metabolite derived from plants, acting as an allosteric modulator of GABAARs, would protect against the negative effects induced by the administration of fipronil in adults of the fruit fly Drosophila melanogaster. We further evaluated our hypothesis via co-administration of flumazenil, a competitive antagonist on the GABAAR, and through in silico analyses. We administered kaempferol prophylactically at three concentrations (10, 30 and 50 µmol l−1) and evaluated its protective effects against motor impairments induced by fipronil. We then used a single dose of kaempferol (50 µmol l−1) to evaluate its protective effect while administering flumazenil. We found that oral administration of fipronil impaired motor control and walking ability. In contrast, kaempferol was innocuous and protected flies from developing the motor-impaired phenotype, whereas the co-administration of flumazenil counteracted these protective effects. These results are supported by the binding of the ligands with the receptor. Together, our results suggest that kaempferol exerts a protective effect against fipronil via positive allosteric modulation of GABAARs, probably within brain areas such as the central complex and the mushroom bodies. These findings further support current attempts to use metabolites derived from plants as protectors against impairments produced by pesticides.

For almost a century, the subtle, yet significant effects of exposure to sublethal levels of pesticides have been of concern because of their impacts ranging from food security to pollinator declines and public health. A renowned example is the association between neuroactive insecticides and the decline of bee populations (Goulson et al., 2015; Holder et al., 2018; Sánchez-Bayo and Wyckhuys, 2019). Neonicotinoids, insecticides acting as agonists of acetylcholine receptors (AChRs), are broadly used for pest control because of their presumed low affinity for human AChRs. However, the exposure to neonicotinoids impairs learning and memory, motor control and sensory sensitivity of beneficial species such as honey bees, even at very low dosages (Charreton et al., 2015; Decourtye et al., 2004; Li et al., 2019; Williamson et al., 2014). Similarly, in humans, exposure to pesticides is associated with a higher risk of developing neurocognitive diseases, such as Parkinson's disease (Elbaz et al., 2009; Godinho et al., 2016; Narayan et al., 2017; Ratner et al., 2014). These concerns have led to regulations and bans that must be supported by further understanding of the mechanisms of action of the pesticides and by strategies to protect animals, including humans, before and after exposure.

Among pesticides, fipronil indiscriminately affects invertebrates and vertebrates through the blockage of inhibitory pathways. Fipronil acts as a non-competitive antagonist of the GABAA receptor (GABAAR), blocking inward currents of chloride ions (Ikeda et al., 2004; Zhao et al., 2003). Hence, not surprisingly, exposure to fipronil leads to neuronal hyperexcitation, loss of motor control and a decrease in memory retention (Bharatiya et al., 2020; Godinho et al., 2016; Park et al., 2016; Stehr et al., 2006; Suzuki et al., 2021; Zaluski et al., 2015). Moreover, exposure to fipronil leads to reduced dopamine and tyrosine hydroxylase levels, thus deteriorating motor activity (Bharatiya et al., 2020; Park et al., 2016), characterizing the so-called ‘parkinsonian phenotype’ in rats (Park et al., 2016). In fact, in humans, this impact on the GABAergic and dopaminergic pathways partially explains the association between the onset of parkinsonian traits such as tremors, dyskinesia and in some cases seizures following exposure to fipronil (Mohamed et al., 2004).

Over recent decades, flavonoids, secondary metabolites derived from plants (Jan and Abbas, 2018), have been recognized as inducers of physiological protection (Jung and Kim, 2018; Maher, 2017), antioxidants (Williams et al., 2004), anti-inflammatory agents (Hämäläinen et al., 2007), modulators of GABAergic activity (Benkherouf et al., 2019; Eghorn et al., 2014; Hall et al., 2014) and inducers of detoxification against xenobiotics (Bernklau et al., 2019; Liao et al., 2017). For example, in honey bees, prophylactic administration of quercetin and p-coumaric acid increases longevity in individuals exposed to low concentrations of imidacloprid (Wong et al., 2018), propiconazole and chlorantraniliprole (Liao et al., 2020). Also, in bumble bees, the prophylactic administration of rutin induces protection against the cognitive impairment produced by oral administration of imidacloprid and fipronil (Riveros and Gronenberg, 2022). However, the potential effect of flavonoids in the regulation of GABAergic pathways supporting motor control has not been addressed. Can prophylactic administration of a flavonoid protect against impairments to motor control produced by exposure to fipronil?

Here, our goal was threefold. First, we aimed to test the impairment produced by the administration of sublethal doses of fipronil on motor control in the fruit fly Drosophila melanogaster. We used the fruit fly because it has been a central model in the understanding of behavioral ecotoxicology and in the development of new drugs (Nagoshi, 2018; Tasman et al., 2021a,b). In D. melanogaster, exposure to pesticides negatively impacts reproduction (Tasman et al., 2021b), immunity (Daisley et al., 2017), memory (Tasman et al., 2021a) and locomotion (Chaudhuri et al., 2020), and induces neuronal death (Martelli et al., 2020). Moreover, D. melanogaster has been used to understand the effects of exposure to pesticides such as paraquat and rotenone in connection with the development of Parkinson's disease (Bastías-Candia et al., 2019; Coulom and Birman, 2004; Maher, 2019; Nagoshi, 2018; Pallanck and Whitworth, 2005; Park et al., 2012).

Second, we aimed to test the protective effect against fipronil of the prophylactic administration of the flavonol kaempferol, primarily recognized by its antiviral (Mitrocotsa et al., 2000), anti-inflammatory (Hämäläinen et al., 2007) and antioxidant (Grünz et al., 2012) effects. However, several accounts suggest an impact on the nervous system because of its role as anxiolytic (Grundmann et al., 2009; Grünz et al., 2012) and inducer of neuroprotective effects (Filomeni et al., 2012). Remarkably, kaempferol might act as a positive allosteric modulator of the GABAAR (Liu et al., 2018), as suggested by the abolishment of its protective effects following the co-administration of flumazenil (Grundmann et al., 2009), a competitive antagonist of the benzodiazepine site (Whitman and Amrein, 1995). Thus, we hypothesized that the allosteric modulation of the GABAergic pathway by kaempferol enhances the release of GABA, thus counteracting the impairment produced by fipronil.

As a third and final goal, we evaluated whether any protective effect produced by kaempferol was dose dependent and whether it would be absent in the presence of flumazenil. As a part of our third goal, we analyzed in silico the effect of administration of kaempferol, flumazenil and fipronil by focusing on the interaction between kaempferol, flumazenil and fipronil with the GABAAR in flies and humans.

Subject maintenance

We commercially acquired wild-type fruit flies (Drosophila melanogaster Meigen 1830) from a local provider (Universidad Nacional de Colombia, Bogotá, Colombia) and reproduced them and maintained them under laboratory conditions (25°C, 12 h:12 h light:dark cycle). Flies were fed using a diet modified from the recipe proposed by the Bloomington Drosophila Stock Center (Indiana University at Bloomington). The diet is composed of agar (1.6 g), sucrose (5 g), corn starch (6 g), yeast (2 g) and water (to a final volume of 100 ml). Diet was supplemented with propionic acid (0.3% v/v) and replaced once a week if needed to avoid mold growth. For all experiments, we arbitrarily selected females (most likely not virgins as they were collected 3 days after eclosion; see below), aiming to decrease a potential confounding effect of sex.

Selection of concentrations of fipronil, kaempferol and flumazenil

The information available on doses of fipronil is scarce and whereas concentrations are reported, individual dosages are not available. Thus, we utilized our empirically determined dose of 1 ng per individual based on our separate accounts in honey bees and bumble bees indicating sublethal cognitive impairment (Riveros and Gronenberg, 2022; L. M. Garcia, V. Caicedo-Garzón and A.J.R., unpublished). In contrast, the concentrations of kaempferol were selected based on the protective and antioxidant effect demonstrated in a model of Parkinson's disease in D. melanogaster (Rahul et al., 2020). All solutions of kaempferol (CAS: 520-18-3; Sigma, St Louis, MO, USA) were obtained through dilutions of a sonicated (90% of maximum intensity) stock (170 µmol l−1) and varied depending on the experiment between 10 and 50 µmol l−1. Finally, as information on in vivo concentrations of flumazenil in D. melanogaster and other insects is scarce, we decided to adjust the concentration from the dose that reverses the anesthetic effect of midazolam, a benzodiazepine, and has no behavioral effects on the invertebrate Daphnia pulex (Dong et al., 2013).

Experiment 1: dose-dependent effect of kaempferol against fipronil

Three days after hatching, we collected females and individually kept them for 24 h of starvation in the wells of culture plates (2 ml volume per well). Following this period of starvation, we randomly assigned each fly to one of four feeding treatments (Table 1): sucrose water (Control: 146 mmol l−1 sucrose solution), 10 µmol l−1 kaempferol (K10), 30 µmol l−1 kaempferol (K30) or 50 µmol l−1 kaempferol (K50). After 22 h, we randomly assigned each fly to one of two treatments: 146 mmol l−1 sucrose water or 2.06 µmol l−1 fipronil (Fip; Astuto 200 scv; ∼0.9 ng of pesticide per fly). Thus, each fly belonged to one of eight treatments depending upon their feeding schedule: Control, Fip, K10, K30, K50, K10/Fip, K30/Fip, K50/Fip. In all cases, solutions were prepared by adding 1 µl of the main component (e.g. 1 µl of 30 µmol l−1 kaempferol) to a sucrose solution (final concentration of sucrose: 146 mmol l−1). Lastly, we conducted evaluation of motor activity 4 h after the administration of fipronil (see below).

Table 1.

Experimental design for the evaluation of cognitive protection against a commercial form of fipronil

Experimental design for the evaluation of cognitive protection against a commercial form of fipronil
Experimental design for the evaluation of cognitive protection against a commercial form of fipronil

Experiment 2: effect of co-administration of flumazenil on kaempferol protection against fipronil

Following our results from the first experiment, we tested whether the effect of kaempferol was due to its action as a positive allosteric modulator of GABAARs. We selected the K50 treatment because its administration led to protection against the deleterious effects of fipronil in most behavioral assays.

Thus, 3 days after hatching, we collected females and individually kept them for 24 h of starvation in 24 well culture plates. Following starvation, we randomly assigned each fly to one of four feeding treatments: sucrose water (Control: 146 mmol l−1 sucrose solution), 50 µmol l−1 kaempferol (K50), 20 µmol l−1 flumazenil (Flum; Diazenil) or 50 µmol l−1 kaempferol+20 µmol l−1 flumazenil (K50/Flum). After 22 h, we randomly assigned each fly to one of two treatments: 146 mmol l−1 sucrose water or 2.06 µmol l−1 fipronil (Fip; Astuto 200 scv; ∼0.9 ng of pesticide per fly). Thus, each fly belonged to one of eight treatments depending upon their feeding schedules across the three phases: Control, K50, Flum, K50/Flum, Fip, Flum/Fip, K50/Fip, K50/Flum/Fip. In all cases, solutions were prepared by adding 1 µl of the main component (e.g. 1 µl of 50 µmol l−1 kaempferol) to a sucrose solution (final concentration of sucrose: 146 mmol l−1). Lastly, we conducted the evaluation of motor activity 4 h after the administration of fipronil.

Motor evaluation

Open-field tests

We evaluated the tendency of the flies to walk and remain on the border of a container by recording their individual trajectories within the cell for 1 min (30 frame s−1, iPhone 6, Apple Inc.). The distance and trajectory of each fly inside the cell were analyzed using the plugin Mtrack and the Scholl's analysis from Fiji (sampling rate: 5 Hz; Schindelin et al., 2012). We defined the ‘edge’ as a region 0.3 cm from the wall (one body length) and counted the total number of intersections outside this range (Fig. 1A). We counted the number of flies that completed at least one lap inside the edge region (nonelap). Then, the border preference index (BPI) was defined as [(nonelap/ntotal)×100]. The time spent walking was recorded manually from the videos using a chronometer. All treatments were represented during each run of the test.

Fig. 1.

Experimental setup for motor evaluation. (A) The behavior of individual flies inside the well was recorded for 60 s. We calculated the distance traveled, border preference index (BPI) and number of intersections. (B) Each fly was transferred to an individual lane and, after 10 s, the total distance traveled and percentage of flies below 5 cm were calculated.

Fig. 1.

Experimental setup for motor evaluation. (A) The behavior of individual flies inside the well was recorded for 60 s. We calculated the distance traveled, border preference index (BPI) and number of intersections. (B) Each fly was transferred to an individual lane and, after 10 s, the total distance traveled and percentage of flies below 5 cm were calculated.

Negative geotaxis test

We transferred each fly to an individual lane of a 3D printed arena and allowed them to acclimate for 5 min before conducting the test (Fig. 1B). The arena consisted of 15 lanes (each lane 7.5×0.4×0.2 cm) covered with three microscope slides. Once flies were in the arena, we taped it to displace the flies to the bottom. We recorded the flies using a camera (30 frames s−1, iPhone 6) and measured the walking distance after 10 s using Fiji (Schindelin et al., 2012). We further calculated the percentage of flies below a threshold of 5 cm. All treatments were represented during each run of the test and each fly was recorded once.

Experiment 3: bioinformatic analysis of the interaction between ligands and the GABAAR

Homology modeling

We aligned the D. melanogaster GABAAR (DGABAAR) amino acid sequence obtained from NCBI (NP_523991.2) against the PDB database to find the most suitable template (PDB 5TIN). We generated 3D homology modeling of DGABAAR using the program MODELLER 9.3 (Webb and Sali, 2016). We then scored 10 models using MODELLER's DOPE scoring function and selected the lowest for further analysis. Finally, we selected the best model after validation of the alignment using the stereochemical quality (Ramachandran plot) and by comparing the Z-score of the modeled GABAAR with the Z-scores of the experimental structures of proteins of the same size.

Molecular dockings

We conducted molecular dockings of the interaction between kaempferol, flumazenil, fipronil and fipronil sulfone, the principal metabolite of fipronil, which is also toxic (Suzuki et al., 2021), with the allosteric site of GABAAR from D. melanogaster and humans using the open-source software Chimera UCSF (v.1.15 for Mac; Pettersen et al., 2004). We performed the analyses using Autodock vina (v.1.1.2 for Mac; Trott and Schroer, 2010). The binding site of flumazenil and kaempferol was located between the residues Ser205 and Tyr58 in D. melanogaster. The binding site for fipronil and fipronil sulfone both in D. melanogaster and humans in GABAA was located around the centroid of the residues −2′ and 2 from the five chains.

Data analyses

We analyzed the normality and homoscedasticity parameters of continuous variables (distance, time, velocity) using a Kolmogorov–Smirnov and Levene test, respectively (Massey, 1951; Glass, 1966). Data exhibiting a normal distribution and homoscedasticity were analyzed using an ANOVA followed by multiple t-test planned comparisons (Student, 1908). Of interest were the following comparisons (depending on the experiment, see above): Control versus Fip to evaluate the effect of the insecticide, Control versus K (10, 30, 50) to evaluate the innocuousness of kaempferol, K (10, 30, 50) versus K (10, 30, 50)/Fip, to test whether there was full protection against Fip, and Fip versus K (10, 30, 50)/Fip, Fip versus K50/Flum/Fip to test whether there was any significant protection. If the data did not exhibit a normal distribution, we compared the treatments using a Kruskal–Wallis test (Hoffman, 2019) followed by the Wilcoxon signed-rank test (for planned comparisons, see above; Woolson, 2008). Nominal variables (e.g. yes/no climbed above threshold or inside the edge region) were analyzed using a chi-square test (Tallarida and Murray, 1987). We identified and excluded outliers after a single run of the Mahalanobis distance test (Rousseeuw and van Zomeren, 1990). Analyses of behavioral results were performed in JMP v.14.2 (SAS Institute). In all cases, error associated with multiple comparisons was corrected with the false discovery rate (FDR) method (Benjamini and Hochberg, 1995) and the corrected P-values (q-values) are presented.

Experiment 1: dose-dependent effect of kaempferol against fipronil

Open field-test

We collected, maintained and evaluated 206 flies. Some flies were excluded following an outlier analysis (see above). Thus, for our final analysis, we included 194 flies distributed across eight treatments: Control (N=26), K10 (N=23), K30 (N=28), K50 (N=19), Fip (N=23), K10/Fip (N=26), K30/Fip (N=27) and K50/Fip (N=22). Fig. 2 shows the results of the open field test for experiment 1, as trajectories, number of intersections, BPI, distance traveled, time walking and average speed.

Fig. 2.

Effects of kaempferol and fipronil on walking control. (A) Overlapping trajectories within the circular cell for flies in the four feeding treatments [control and 10, 30 or 50 μmol l−1 kaempferol (K10, K30, K50)] with or without fipronil treatment (Fip). (B) Number of intersections of trajectories with an imaginary circular region one body length away from the border. (C) BPI. (D) Distance traveled. (E) Time walking. (F) Average speed. Data are means±s.e.m. *P<0.05, **P<0.01 compared with flies in the Fip group; P<0.05, ‡‡P<0.01, ‡‡‡P<0.001 compared with flies in the Control group. P-values are reported after false discovery rate (FDR) correction. Only significant planned comparisons are shown (see Materials and Methods). Sample sizes (N) are given within the bars.

Fig. 2.

Effects of kaempferol and fipronil on walking control. (A) Overlapping trajectories within the circular cell for flies in the four feeding treatments [control and 10, 30 or 50 μmol l−1 kaempferol (K10, K30, K50)] with or without fipronil treatment (Fip). (B) Number of intersections of trajectories with an imaginary circular region one body length away from the border. (C) BPI. (D) Distance traveled. (E) Time walking. (F) Average speed. Data are means±s.e.m. *P<0.05, **P<0.01 compared with flies in the Fip group; P<0.05, ‡‡P<0.01, ‡‡‡P<0.001 compared with flies in the Control group. P-values are reported after false discovery rate (FDR) correction. Only significant planned comparisons are shown (see Materials and Methods). Sample sizes (N) are given within the bars.

We found that the trajectory patterns differed among treatments in terms of BPI and the number of intersections (Kruskal–Wallis test: χ32=34.21, P≤0.0001). Flies fed with fipronil (Fip group) changed their trajectory patterns relative to Control flies (Fig. 2A–C). For instance, Control flies followed circular trajectories, whereas Fip flies typically exhibited trajectories often distanced from the wall. Thus, the trajectories of Fip flies had more intersections than those of Control flies (mean±s.e.m. number of intersections: Control 1.6±0.47, Fip 6.5±1.05; Wilcoxon test: Z=−4.292, P=0.039; Fig. 2B) and shorter BPI (χ21,47=19.67, P=0.012; Fig. 2C). Additionally, we found that the administration of fipronil and kaempferol affected the distance traveled (ANOVA: F7,186=5.46, P<0.0001), time (Kruskal–Wallis test: χ32=57.18, P<0.0001) and speed (Kruskal–Wallis test: χ32=21.57, P=0.003). For instance, flies in the Fip group spent a significantly shorter time walking and, when moving, traveled significantly shorter distances than flies fed with sucrose solution (Fip versus Control: mean±s.e.m. time: Control 36.92±2.7 s, Fip 15.13±3.21 s; Wilcoxon test: Z=−4.77, P=0.0005; Fig. 2E; mean±s.e.m. distance: Control 20.45±1.4 cm, Fip 12.10±1.11 cm; t182=−4.92, P=0.0005; Fig. 2D). Interestingly, Fip flies were significantly faster than Control flies (Control 0.58±0.03 cm s−1, Fip 1.02±1.02 cm s−1; Wilcoxon test: Z=2.57, P=0.025; Fig. 2F).

In contrast, administration of kaempferol was innocuous to the flies at all concentrations. Flies across all concentrations of kaempferol exhibited trajectory patterns that did not differ from those of Control flies (mean±s.e.m. number of intersections: K10 2.78±0.65; Wilcoxon test: Z=−0.50, P=0.61; K30 2.73±0.7; Wilcoxon test: Z=−1.15, P=0.31; K50 2.5±0.6; Wilcoxon test: Z=−1.17, P=0.31; Fig. 2A,B; BPI: K10: χ21,51=0.38, P=0.56; K30: χ21,56=0.487, P=0.60; K50: χ21,47=0.02, P=0.89; Fig. 2A,C). Further, flies in the K10, K30 and K50 groups exhibited a performance in terms of distance, time and speed of trajectories that did not differ from that of Control flies (mean±s.e.m. distance: Control 20.45±1.4 cm, K10 19.2±1.04 cm; t182: 0.71, P=0.096; K30 18.37±1.45 cm; t182: 1.65, P=0.14; K50 18.33±0.67 cm; t182: 1.19, P=0.24; Fig. 2D; mean±s.e.m. time: Control 36.92±2.7 s, K10 36.17±2.21 s; Wilcoxon test: Z=0.06, P=0.952; K30 37.32±2.57 s; Wilcoxon test: Z=0.06, P=0.952; K50 37.75±2.34 s; Wilcoxon test: Z=−0.66, P=0.64; Fig. 2E; mean±s.e.m. speed: Control 0.58±0.03 cm s−1, K10 0.65±0.68 cm s−1; Wilcoxon test: Z=0.27, P=0.23; K30 0.50±0.03 cm s−1; Wilcoxon test: Z=2.13, P=0.062; K50 0.51±0.04 cm s−1; Wilcoxon test: Z=2.33, P=0.05; Fig. 2F).

Most remarkably, the prophylactic administration of kaempferol led to an overall dose-dependent performance improvement in flies exposed to fipronil. Flies in the K50/Fip treatment followed trajectories that did not differ from those of Control flies in the total number of intersections or BPI (mean±s.e.m. number of intersections: K50 2.5±0.6, K50/Fip 4.32±1.35; Wilcoxon test: Z=0.7221, P=0.522; BPI: χ21,41=1.420, P=0.17; Fig. 2A–C), but significantly differed from those of Fip flies (mean±s.e.m. number of intersections K50/Fip 4.32±1.35; Wilcoxon test: Z=0.72, P=0.045; BPI: χ21,41=1.42, P=0.007; Fig. 2A–C), demonstrating full protection. However, this protection was not observed when we prophylactically administered the lower concentrations of kaempferol in the K10/Fip or the K30/Fip flies (Fip versus 10/30K/Fip: mean±s.e.m. number of intersections: K10/Fip 5.42±1.03; Wilcoxon test: Z=−0.867, P=0.31; K30/Fip 7.90±1.29; Wilcoxon test: Z=0.722, P=0.31; Fig. 2A,B; BPI: K10/Fip: χ21,49=11.31, P=0.08; K30/Fip: χ21,55=8.20, P=0.13; Fig. 2A,C). Importantly, flies in the K10/Fip, K30/Fip and K50/Fip groups had walking distances that did not significantly differ from those of Fip flies (mean±s.e.m. distance: K10/Fip 14.61±1.28 cm; t182=−148, P=0.12; K30/Fip 15.06±1.27 cm; t182=−0.85, P=0.22; K50/Fip 13.86±1.29 cm; t182=−1.0, P=0.2; Fig. 2D), but significantly differed from those of Control flies (K10/Fip: t182=3.56, P=0.02; K30/Fip: t182=4.10, P=0.02; K50/Fip: t182=3.84, P=0.02; Fig. 2D). Further, flies that were fed with kaempferol spent a significantly longer time walking than flies in the Fip group (mean±s.e.m. time: Fip 7.37±0.62 s, K10/Fip 24.58±3.09 s; Wilcoxon test: Z=−2.0, P=0.004; K30/Fip 23.04±1.68 s; Wilcoxon test: Z=−3.0, P=0.003; K50/Fip 26.95±2.66 s; Wilcoxon test: Z=−3.32, P=0.002; Fig. 2E), but a significantly shorter time than their respective controls (K10/Fip: Wilcoxon test: Z=2.77, P=0.012; K30/Fip: Wilcoxon test: Z=3.96, P=0.002; K50/Fip: Wilcoxon test: Z=2.44, P=0.027; Fig. 2E). We further observed protection when evaluating the speed of K30/Fip and K50/Fip flies (30/50K/Fip versus Fip: mean±s.e.m. speed K30/Fip 0.75±0.11 cm s−1; Wilcoxon test: Z=2.10, P=0.049; K50/Fip 0.68±0.15 cm s−1; Wilcoxon test: Z=2.67, P=0.025; Fig. 2F). Additionally, protected flies walked at speeds that did not differ significantly from those of Controls (K10/Fip: Wilcoxon test: Z=−0.70, P=0.232; K30/Fip: Wilcoxon test: Z=−0.57, P=0.062; K50/Fip: Wilcoxon test: Z=0.55, P=0.232; Fig. 2F). These results suggest that K30 and K50 are the most protective concentrations, while K10 seems to lead to an intermediate level of protection.

Negative geotaxis test

We evaluated geotaxis in a total of 192 flies distributed across eight treatments: Control (N=25), Fip (N=25), K10 (N=23), K30 (N=25), K50 (N=19), K10/Fip (N=22), K30/Fip (N=28), K50/Fip (N=25). We found that feeding flies with fipronil or kaempferol had an effect on the distance climbed (Kruskal–Wallis test: χ32=48.70, P<0.0001). Flies in the Fip group climbed shorter distances compared with flies in the Control group (mean±s.e.m. distance: Control 5.43±0.40 cm, Fip 1.26±0.25 cm; Wilcoxon test: Z=−5.292, P=0.003; Fig. 3A). Also, 78.9% of flies in the Fip group remained below the 5 cm threshold while only 14.3% of the flies in the Control group remained below 5 cm (χ21,50=4.348, P=0.003; Fig. 3B). Prophylactic treatment with K10 and K50 increased the distance traveled but flies in the K30/Fip treatment did not differ from those in the Fip group (mean±s.e.m. distance: K10/Fip 3.23±0.54 cm; Wilcoxon test: Z=2.59, P=0.01; K30/Fip 2.88±0.52 cm; Wilcoxon test: Z=1.67, P=0.08; K50/Fip 3.22±0.44 cm; Wilcoxon test: Z=3.44, P=0.003; Fig. 3A). Moreover, all groups receiving kaempferol showed a decrease in the percentage of flies below the 5 cm threshold compared with the Fip group (K10/Fip: χ21,46=9.83, P=0.003; K30/Fip: χ21,53=9.68, P=0.003; K50/Fip: χ21,49=7.12, P=0.01; Fig. 3B). Notably, only the K50/Fip group behaved similar to its control in distance (mean±s.e.m. distance: K10/Fip 3.23±0.54 cm; Wilcoxon test: Z=−2.53, P=0.012; K30/Fip 2.88±0.52 cm; Wilcoxon test: Z=−2.59, P=0.012; K50/Fip 3.22±0.43 cm; Wilcoxon test: Z=−1.24, P=0.15; Fig. 3A) and the percentage of flies below 5 cm (K10/Fip: χ21,46=6.88, P=0.02; K30/Fip: χ21,53=8.40, P=0.02; K50/Fip: χ21,49=10.8, P=0.09; Fig. 3B), suggesting partial protection by kaempferol at the lowest concentrations. Importantly, flies exposed to kaempferol did not significantly differ in distance traveled or percentage of flies below the 5 cm threshold from those in the Control group (mean±s.e.m. distance: K10 5.14±0.40 cm; Wilcoxon test: Z=−0.730, P=0.47; K30 4.93±0.40 cm; Z=−1.24, P=0.24; K50 4.22±0.62 cm; Z=−1.44, P=0.19; percentage of flies below 5 cm: K10: χ21,49=0.16, P=0.69; K30: χ21,50=0.37, P=0.60; K50: χ21,44=2.76, P=0.12; Fig. 3).

Fig. 3.

Effects of kaempferol on the negative geotactic response after exposure to fipronil. (A) Distance traveled. (B) Percentage of flies below the 5 cm threshold. Data are means±s.e.m. *P<0.05, **P<0.01 compared with flies in the Fip group; P<0.05, ‡‡P<0.01 compared with flies in the Control group. P-values are reported after FDR correction. Sample sizes (N) are given within the bars.

Fig. 3.

Effects of kaempferol on the negative geotactic response after exposure to fipronil. (A) Distance traveled. (B) Percentage of flies below the 5 cm threshold. Data are means±s.e.m. *P<0.05, **P<0.01 compared with flies in the Fip group; P<0.05, ‡‡P<0.01 compared with flies in the Control group. P-values are reported after FDR correction. Sample sizes (N) are given within the bars.

Experiment 2: effect of co-administration of flumazenil on kaempferol protection against fipronil

Open field

We collected, maintained and evaluated 266 flies. Some flies were excluded following an outlier analysis (see above). Thus, for our final analysis, we included 225 flies distributed across eight treatments: Control (N=31), K50 (N=33), Flum (N=31), K50/Flum (N=28), Fip (N=24), K50/Fip (N=30), K50/Flum/Fip (N=20), Flum/Fip (N=28). Fig. 4 shows the results of the open field test for experiment 2, as trajectories, number of intersections, BPI, distance traveled, time walking and average speed.

Fig. 4.

Effects of co-administration of flumazenil, kaempferol and fipronil on walking control. (A) Overlapping trajectories within the circular cell for flies in the four treatment groups [control, 50 μmol l−1 kaempferol (K50), 20 μmol l−1 flumazenil (Flum) and 50 μmol l−1 kaempferol+20 μmol l−1 flumazenil (K50/Flum)], with or without fipronil treatment. (B) Number of intersections of trajectories with an imaginary circular region one body length away from the border. (C) BPI. (D) Distance traveled. (E) Time walking. (F) Average speed. Data are means±s.e.m. *P<0.05, **P<0.01, ***P<0.001 compared with flies in the Fip group; §P<0.05, §§P<0.001 relative to flies in the K50/Fip group; and P<0.05, ‡‡P<0.01, ‡‡‡P<0.001 compared with flies in the Control group. P-values are reported after FDR correction. Only significant planned comparisons are shown (see Materials and Methods). Sample sizes (N) are given within the bars.

Fig. 4.

Effects of co-administration of flumazenil, kaempferol and fipronil on walking control. (A) Overlapping trajectories within the circular cell for flies in the four treatment groups [control, 50 μmol l−1 kaempferol (K50), 20 μmol l−1 flumazenil (Flum) and 50 μmol l−1 kaempferol+20 μmol l−1 flumazenil (K50/Flum)], with or without fipronil treatment. (B) Number of intersections of trajectories with an imaginary circular region one body length away from the border. (C) BPI. (D) Distance traveled. (E) Time walking. (F) Average speed. Data are means±s.e.m. *P<0.05, **P<0.01, ***P<0.001 compared with flies in the Fip group; §P<0.05, §§P<0.001 relative to flies in the K50/Fip group; and P<0.05, ‡‡P<0.01, ‡‡‡P<0.001 compared with flies in the Control group. P-values are reported after FDR correction. Only significant planned comparisons are shown (see Materials and Methods). Sample sizes (N) are given within the bars.

We found that the administered compounds affected the distance traveled (ANOVA: F7,217=19.11, P<0.0001), time spent walking (Kruskal–Wallis test: χ32=120.4, P<0.0001), speed (Kruskal–Wallis test: χ32=55.50, P<0.0001) and trajectory patterns (intersections: Kruskal–Wallis test: χ32=24, P=0.001). We found that the co-administration of flumazenil to the flies counteracted the protective effect of administration of kaempferol (K50/Fip versus K50/Flum/Fip) in terms of the number of intersections (mean±s.e.m. number of intersections: K50Flum/Fip 4.15±0.89, K50/Fip 1.80±0.34; Wilcoxon test: Z=1.95, P=0.04; Fig. 4A,B), walking time (mean±s.e.m. time: K50Flum/Fip 8.65±0.98 s, K50/Fip 15.67±1.58 s; Wilcoxon test: Z=3.12, P=0.002; Fig. 4E) and speed (mean±s.e.m. speed K50Flum/Fip 0.90±0.05 cm s−1, K50/Fip 0.88±0.004 cm s−1; Wilcoxon test: Z=2.68, P=0.02; Fig. 4F). Interestingly, the BPI and the total distance walked by flies in the K50/Flum/Fip group were not significantly different from those of flies in the K50/Fip group (BPI: χ21,50=3, P=0.08; Fig. 4D; mean±s.e.m. distance: K50/Flum/Fip 8.77±0.96 cm, K50/Fip 12.11±1.12 cm; Wilcoxon test: Z=−1.73, P=0.08; Fig. 4E). Thus, flies receiving co-treatment with flumazenil seemed to emulate the performance of flies receiving fipronil alone (mean±s.e.m. number of intersections: Fip 4.04±0.77; Wilcoxon test: Z=0.097, P=0.55; Fig. 4A,B; BPI: χ21,44=1.13, P=0.33; Fig. 4A,C; mean±s.e.m. distance: Fip 8.38±0.74 cm; Wilcoxon test: Z=0.80, P=0.43; Fig. 4D; mean±s.e.m. time: Fip 7.37±0.62 s; Wilcoxon test: Z=1.40, P=0.19;Fig. 4E; mean±s.e.m. speed: Fip 1.15±0.07 cm s−1; Wilcoxon test: Z=−0.91, P=0.22; Fig. 4F).

Co-treatment with flumazenil led to a shorter time walking and a lower speed but not to differences in trajectory patterns or distance compared with flies in the Control group (mean±s.e.m. number of intersections: Control 1.065±0.29, K50/Flum 1.71±0.40; Wilcoxon test: Z=−1.18, P=0.38; Fig. 4A,B; BPI: χ21,60=0.35, P=0.55; Fig. 4C; mean±s.e.m. distance: Control 20.03±1.31 cm, K50/Flum 17.58±1.33 cm; Wilcoxon test: Z=1.34, P=0.40; Fig. 4D; mean±s.e.m. time: Control 30.06±1.90 s, K50/Flum 22.96±2.14 s; Z=2.06, P=0.02; Fig. 4E; mean±s.e.m. speed: Control 0.69±0.03 cm s−1, K50/Flum 0.83±0.04 cm s−1; Wilcoxon test: Z=2.22, P=0.03; Fig. 4F). Importantly, flumazenil alone did not reverse the effects of fipronil; in fact, flies in the Flum/Fip group did not perform differently from flies in the Fip group (mean±s.e.m. number of intersections: Flum/Fip 3.36±0.76; Z=−1.35, P=0.15; Fig. 4A,B; BPI: χ21,52=1.82, P=0.12; Fig. 4A,C; mean±s.e.m. distance: Flum/Fip9.38±0.66 cm; Z=0.80, P=0.29; mean±s.e.m. time: Flum/Fip 8.96±0.71 s; Z=1.40, P=0.1; Fig. 4D; mean±s.e.m. speed: Flum/Fip 1.14±0.08 cm s−1; Z=−0.91, P=0.40; Fig. 4F). Additionally, flies in the Flum group behaved similar to flies in the Control group, except for time walking and speed (mean±s.e.m. number of intersections: Control 1.06±0.29, Flum 1.13±0.31; Z=−1.68, P=0.35; Fig. 4A,B; BPI: χ21,63=2.91, P=0.12; Fig. 4A,C; mean±s.e.m. distance: Control 20.03±1.31 cm, Flum 18.10±1.24 cm; Z=0.23, P=0.40; Fig. 4D; mean±s.e.m. time: Control 30.06±1.90 s, Flum 23.45±1.98 s; Z=2.52, P=0.02; Fig. 4E; mean±s.e.m. speed: Control 0.69±0.03 cm s−1, Flum 0.84±0.04 cm s−1; Z=−2.65, P=0.02; Fig. 4F).

Negative geotaxis test

We evaluated geotaxis in a total of 247 flies distributed across eight treatments: Control (N=28), K50 (N=37), Flum (N=32), K50/Flum (N=28), Fip (N=31), K50/Fip (N=30), K50/Flum/Fip (N=30), Flum/Fip (N=31). We found that these compounds affected the distance climbed (Kruskal–Wallis test: χ32=38.48, P<0.0001). For instance, the co-administration of flumazenil and kaempferol significantly reduced the distance climbed and increased the percentage of flies under the 5 cm threshold compared with that in the K50/Fip group (mean±s.e.m. distance: K50/Fip 3.42±0.48 cm, K50/Flum/Fip 2.08±0.42 cm; Z=−2.08, P=0.04; Fig. 5A; percentage of flies below 5 cm: χ21,58=8.73, P=0.03; Fig. 5B). Moreover, flies in the K50/Flum/Fip group did not differ from those in the Fip group in terms of distance or percentage of flies below the 5 cm threshold (mean±s.e.m. distance: K50/Flum/Fip 2.077±0.22 cm, Fip 2.14±0.42 cm; Z=0.62, P=0.36; percentage of flies below 5 cm: χ21,62=0.1, P=0.40; Fig. 5A,B). Importantly, administration of flumazenil alone did not have any effect against fipronil (Flum/Fip versus Fip: mean±s.e.m. distance: Flum/Fip 2.47±0.45 cm; Z=0.62, P=0.42;Fig. 5A; percentage of flies below the 5 cm threshold: χ21,62=0.10, P=0.47; Fig. 5B). Finally, the performance of flies in the Flum group did not significantly differ from that of flies in the Control group (mean±s.e.m. number of intersections: Flum 1.13±0.31, Control 1.06±0.29; Z=−0.05, P=0.959; Fig. 5A; BPI: χ21,60=0.58, P=0.50; Fig. 5B).

Fig. 5.

Effects of co-administration of flumazenil and kaempferol on the negative geotactic response after exposure to fipronil. (A) Distance traveled. (B) Percentage of flies below the 5 cm threshold. Data are means±s.e.m. *P<0.05 relative to flies in the Fip group; §P<0.05 relative to flies in the K50/Fip group; P<0.05 relative to the corresponding control. P-values are reported after FDR correction. Sample sizes (N) are given within the bars.

Fig. 5.

Effects of co-administration of flumazenil and kaempferol on the negative geotactic response after exposure to fipronil. (A) Distance traveled. (B) Percentage of flies below the 5 cm threshold. Data are means±s.e.m. *P<0.05 relative to flies in the Fip group; §P<0.05 relative to flies in the K50/Fip group; P<0.05 relative to the corresponding control. P-values are reported after FDR correction. Sample sizes (N) are given within the bars.

Experiment 3: bioinformatic analysis of the interaction between ligands and the GABAAR

Homology modeling

We found that 88.3% of residues from the modeled DGABAAR were in the most favorable area, 10.8% were in the allowed regions, 0.5% were in the generously allowed regions and just 0.4% were in the disallowed region (Ramachandran plot; Fig. 6A). The Z-scores of the modeled DGABAAR were similar to those obtained by NMA or X-ray data of protein of similar size (Fig. 6B). We considered that the model was of good quality as the Z-scores relied on the range of similar proteins even though the percentage of residues in the most favorable area in the Ramachandran plot was below 90%. However, we took special care in subsequent steps for docking preparation.

Fig. 6.

Structure validation of the modeled GABAARs of Drosophila melanogaster. (A) Ramachandran plot showing 88.3% of residues from the modeled D. melanogaster GABAA receptors (DGABAAR) were in the most favorable area, 10.8% were in the allowed regions, 0.5% were in the generously allowed regions and just 0.4% were in the disallowed region. (B) Estimation of the Z-score of chain A of DGABAAR compared with experimental data. The calculated value is shown as a black dot in the plot.

Fig. 6.

Structure validation of the modeled GABAARs of Drosophila melanogaster. (A) Ramachandran plot showing 88.3% of residues from the modeled D. melanogaster GABAA receptors (DGABAAR) were in the most favorable area, 10.8% were in the allowed regions, 0.5% were in the generously allowed regions and just 0.4% were in the disallowed region. (B) Estimation of the Z-score of chain A of DGABAAR compared with experimental data. The calculated value is shown as a black dot in the plot.

Molecular dockings

We docked kaempferol and flumazenil to the allosteric site, and fipronil and fipronil sulfone to the antagonist site on DGABAAR (Fig. 7A,B). We found that flumazenil and kaempferol bond to the DGABAAR in the canonical allosteric site in the extracellular domain. Flumazenil, the competitive antagonist of benzodiazepines at the GABAAR (Whitman and Amrein, 1995), binds to the DGABAAR through the interaction between the Y254, L249, Y90, Y109, F206 and S205 residues with a free energy (ΔG) of −8.5 kcal mol−1 (Fig. 7C). Kaempferol forms hydrogen bonds with the S205 and S176 residues with the OH-group of the B ring and with the oxygen from the C ring (ΔG=−8.3 kcal mol−1; Fig. 7D). We corroborated that fipronil and its principal metabolite fipronil sulfone binds to the transmembrane domain (TMD) of the five chains of the receptor (Zheng et al., 2014). Fipronil and fipronil sulfone bind with the trifluoromethyl group directly to the intracellular domain and interact with the 2′ Ala, 6′ Thr and 9′ Leu residues from the five chains (Fipronil: ΔG=−8.1 kcal mol−1, fipronil sulfone: ΔG=−7.9 kcal mol−1; Fig. 7E,F).

Fig. 7.

Allosteric and antagonist binding site on the DGABAAR. (A) Structural formulae of flumazenil, kaempferol, fipronil and fipronil sulfone. (B) Top and side view of the allosteric and antagonist binding site in the DGABAAR. TMD, transmembrane domain. (C,D) Detailed architecture of the allosteric binding pocket in complex with (C) flumazenil and (D) kaempferol. (E,F) Detailed architecture of the antagonist binding pocket in complex with (E) fipronil and (F) fipronil sulfone. Black dotted lines indicate hydrogen bonds.

Fig. 7.

Allosteric and antagonist binding site on the DGABAAR. (A) Structural formulae of flumazenil, kaempferol, fipronil and fipronil sulfone. (B) Top and side view of the allosteric and antagonist binding site in the DGABAAR. TMD, transmembrane domain. (C,D) Detailed architecture of the allosteric binding pocket in complex with (C) flumazenil and (D) kaempferol. (E,F) Detailed architecture of the antagonist binding pocket in complex with (E) fipronil and (F) fipronil sulfone. Black dotted lines indicate hydrogen bonds.

Interestingly, we found that kaempferol binds to the allosteric site in human GABAAR and corroborated that fipronil binds to the TMD of human GABAAR (Ci et al., 2007). Also, we found that fipronil sulfone binds to the GABAAR. Kaempferol forms hydrogen bonds with the S205 residue with oxygen from the A and C rings, and interacts with Y160, Y58, F100, R144 and F77 from the binding pocket (ΔG=−8.7 kcal mol−1; Fig. 8A). Fipronil and fipronil sulfone bind to the TMD via the trifluoromethyl group directly to the intracellular domain through interaction with L259, T256, A252 and P253 residues (fipronil: ΔG=−8.4 kcal mol−1, fipronil sulfone: ΔG=−8.0 kcal mol−1; Fig. 8B,C).

Fig. 8.

Detailed binding of kaempferol, fipronil and fipronil sulfone to human GABAARs. (A) Allosteric binding pocket in complex with kaempferol. (B,C) Antagonist site in the TMD in complex with (B) fipronil and (C) fipronil sulfone.

Fig. 8.

Detailed binding of kaempferol, fipronil and fipronil sulfone to human GABAARs. (A) Allosteric binding pocket in complex with kaempferol. (B,C) Antagonist site in the TMD in complex with (B) fipronil and (C) fipronil sulfone.

Recent decades have seen a rise in concerns over the negative effects of exposure to sublethal doses of pesticides in invertebrates and vertebrates. Current approaches include restrictions and bans, yet the continuous development of new pesticides calls for the development of strategies that counteract the deleterious effects in animals. Fipronil, a commonly used insecticide, acts as a GABAA antagonist, promoting neuronal hyperexcitation, dopamine depletion and motor loss (Stehr et al., 2006; Zhao et al., 2003). Here, we evaluated the protective effects of the flavonoid kaempferol against motor impairment induced by the administration of fipronil. We used kaempferol because its neuroprotective properties have been studied in neurodegenerative diseases and because of its action as a potential regulator of the GABAergic system. Our results confirm the detrimental effects of exposure to fipronil on motor activity. The impairments resemble parkinsonian features, characterized by fast walking, shorter walking activity and shorter distances climbed. In contrast, we found that kaempferol was innocuous and conferred protection to the flies against the impairing effects of fipronil. Importantly, the use of co-administration of flumazenil and the bioinformatic analyses enabled us to infer interactions between ligands and GABAAR that help explain the observed protection. Thus, together, our results highlight the motor impairment produced by fipronil, the innocuous effect of kaempferol and the prophylactic protection induced by its administration.

First, we observed that flies exposed to fipronil were faster but spent a shorter time walking, a trait resembling the involuntary hastening and propulsion-like gait in those with Parkinson's disease (PD; Hohler and de Leon, 2011; Koller, 1984; Pallis, 1971), albeit this effect decreased over time (Hohler and de Leon, 2011). Interestingly, whereas the faster pace in people with PD contrasts with a freezing phenotype (Nonnekes et al., 2019), the higher speed in the flies exposed to fipronil reflected bursts of fast walking episodes followed by longer periods of immobility. These higher speeds may be caused by the release of the inhibitory action imposed by the GABAergic pathways within the central complex (CC; a major area within the brain associated with motor control), controlling limb coordination and therefore speed (Martin et al., 1999). The irregular oscillation between periods of immobility and walking was perhaps initiated by impairment of the regulatory control of walk termination supported by the mushroom bodies (MBs; Martin et al., 1998; Poeck et al., 2008). This argument is further supported by the recent discovery of a tract connecting the central complex and the MBs (Li et al., 2020).

A second key feature of the Fip flies was erratic and shorter trajectories distanced from the border, probably explained by an impairment of visual orientation. The exposure to fipronil is associated with decreases in dopamine levels produced indirectly by the impairment of GABAergic pathways (Bharatiya et al., 2020; Błaszczyk, 2016; Kottler et al., 2019). This is consistent with the cognitive impairments in PD patients exhibiting visuo-spatial orientation and attentional deficits that predate motor symptoms and, in some cases, worsen it (Azulay et al., 1999; Davidsdottir et al., 2005). In fruit flies as in other animals, coordinated execution of walking and climbing depends on the integration of proprioceptive and exteroceptive (mostly visual) information (Bausenwein et al., 1994; Strauß and Heisenberg, 1990; Weir and Dickinson, 2015). We hypothesize that fipronil induced changes in the GABAergic system and an indirect reduction in dopamine levels, causing impairment of optomotor responses. Although GABAARs are widely distributed in the fruit fly brain (Enell et al., 2007), the impairment of optomotor responses may result from damage to at least two structures: the CC and the MBs. CC mutants lose their ability to orient in Buridan's paradigm (Strauss and Heisenberg, 1993). Also, orientation memory relies on GABAergic neurons in the ellipsoid body, a sub-region within the CC (Neuser et al., 2008). Further, in fruit flies, border preference is related to spatial orientation, visual processing and dopamine signaling (Kottler et al., 2019; Soibam et al., 2012; White et al., 2010). However, the neuronal circuit is not yet fully understood. Although border preference has been primarily associated with the MBs (Besson and Martin, 2005), a role of the CC cannot be disregarded because of its coordination of movements based on visual inputs (Martin et al., 1999).

In flies, the hallmark for protection in models of motor impairment focuses on a better performance in the climbing assay (Nagoshi, 2018) rather than cognitive features such as attention deficit and visual impairments. Here, we found that kaempferol had protective effects on the detrimental impacts of fipronil in walking features associated with coordination and visual orientation. The protective effect of kaempferol may be explained by its direct action on GABAARs, affecting primarily but not exclusively the visual processing and motor circuits within the CC and the MBs. Kaempferol improves cognitive function via GABAAR modulation (Grundmann et al., 2009); however, the importance of the modulation in improving motor functions is unknown. Interestingly, we found that the co-administration of flumazenil, a competitive antagonist of the allosteric site in GABAARs (Whitman and Amrein, 1995), reversed the protective effects of kaempferol against fipronil. The counteraction was evident through three features of the phenotype: lower walking activity, trajectory patterns separated from the border and shorter distances climbed. This is consistent with our hypothesis that kaempferol acts as a positive allosteric modulator (PAM) of the GABAARs. The effect of PAMs is dependent on ligand concentration, mostly exhibiting a biphasic response (Sigel and Baur, 1988). This supports our results of a dose-dependent protective effect of kaempferol against fipronil.

That kaempferol has a protective effect against fipronil via the GABAergic system may be explained by its action through at least three mechanisms: enhancement of GABAA transmission, blocking the binding of fipronil and regulation of antioxidant activity. First, kaempferol allosterically enhances GABAAR transmission, which prevents hyperexcitation in the system and halts the oxidative stress induced by fipronil and later neuronal loss. This agrees with previous accounts showing that a balance between excitation and inhibition is essential to a normal transmission in the brain. Second, it is possible that the binding of kaempferol in the allosteric site inhibits the binding of fipronil and thus prevents the blockade of chloride influx. For example, PAMs inhibit the binding of picrotoxin, a GABAAR antagonist that binds in the same site as fipronil (Leeb-Lundberg et al., 1981). Also, benzodiazepines (recognized PAMs) have been used to treat fipronil intoxication in humans (Mohamed et al., 2004). Third, in PD models, protective effects of kaempferol are primarily attributable to its antioxidant activity (Li and Pu, 2011; Rahul et al., 2020). We observed that the protective effects were reverted by flumazenil; therefore, we suggest that the antioxidant activity is not the primary mechanism of protection against fipronil, although it cannot be ruled out as an indirect course of protection. Hence, kaempferol may maintain the integrity of neuronal circuits within the CC and the MBs via GABAAR enhancement, allowing an improvement in motor performance and, probably, visual orientation in flies.

Thus, we conclude that acute exposure to fipronil impairs motor control in the fruit fly D. melanogaster while kaempferol protects against motor impairment induced by fipronil in a dose-dependent fashion. The protective effects of kaempferol are presumably the result of its positive allosteric activity on GABAARs. Hence, our results highlight the potential of kaempferol as a prophylactic nutraceutical to protect against motor impairment produced by exposure to fipronil, a pesticide broadly used around the globe.

We thank Julieta Troncoso (Universidad Nacional de Colombia) and Wulfila Gronenberg (University of Arizona) Mauricio Nava (Universidad del Rosario) and Jhon Jairo Sutachán (Pontificia Universidad Javeriana) for comments on early versions of the manuscript. We also thank members of the Riveros' Lab (Grupo de Investigaciones CANNON) for comments on the experimental design. We thank Juan Manuel Quiroga, Lida Triviño and Amparo Duarte for continued technical support.

Author contributions

Conceptualization: D.M.R., K.F.L., A.J.R.; Methodology: D.M.R., A.J.R.; Formal analysis: D.M.R., K.F.L., A.J.R.; Investigation: D.M.R., K.F.L.; Resources: A.J.R.; Data curation: D.M.R.; Writing - original draft: D.M.R., A.J.R.; Writing - review & editing: D.M.R., K.F.L., A.J.R.; Supervision: A.J.R.; Project administration: A.J.R.; Funding acquisition: A.J.R.

Funding

This project was funded by a Big Grant of the Dirección de Investigación e Innovación (DIeI; Universidad del Rosario) to A.J.R. Open Access funding provided by University of Arizona. Deposited in PMC for immediate release.

Azulay
,
J. P.
,
Mesure
,
S.
,
Amblard
,
B.
,
Blin
,
O.
,
Sangla
,
I.
and
Pouget
,
J.
(
1999
).
Visual control of locomotion in Parkinson's disease
.
Brain
122
,
111
-
120
.
Bastías-Candia
,
S.
,
Zolezzi
,
J. M.
and
Inestrosa
,
N. C.
(
2019
).
Revisiting the paraquat-induced sporadic Parkinson's disease-like model
.
Mol. Neurobiol.
56
,
1044
-
1055
.
Bausenwein
,
B.
,
Müller
,
N. R.
and
Heisenberg
,
M.
(
1994
).
Behavior–dependent activity labeling in the central complex of Drosophila during controlled visual stimulation
.
J. Comp. Neurol.
340
,
255
-
268
.
Benjamini
,
Y.
and
Hochberg
,
Y.
(
1995
).
Controlling the false discovery rate: a practical and powerful approach to multiple testing
.
J. R. Stat. Soc. Ser. B
57
,
289
-
300
.
Benkherouf
,
A. Y.
,
Soini
,
S. L.
,
Stompor
,
M.
and
Uusi-Oukari
,
M.
(
2019
).
Positive allosteric modulation of native and recombinant GABAA receptors by hops prenylflavonoids
.
Eur. J. Pharmacol
852
,
34
-
41
.
Bernklau
,
E.
,
Bjostad
,
L.
,
Hogeboom
,
A.
,
Carlisle
,
A.
and
Arathi
,
H. S.
(
2019
).
Dietary phytochemicals, honey bee longevity and pathogen tolerance
.
Insects
10
,
14
.
Besson
,
M.
and
Martin
,
J. R.
(
2005
).
Centrophobism/thigmotaxis, a new role for the mushroom bodies in Drosophila
.
J. Neurobiol
62
,
386
-
396
.
Bharatiya
,
R.
,
Bratzu
,
J.
,
Lobina
,
C.
,
Corda
,
G.
,
Cocco
,
C.
,
De Deurwaerdere
,
P.
,
Argiolas
,
A.
,
Melis
,
M. R.
and
Sanna
,
F.
(
2020
).
The pesticide fipronil injected into the substantia nigra of male rats decreases striatal dopamine content: A neurochemical, immunohistochemical and behavioral study
.
Behav. Brain Res.
384
,
112562
.
Błaszczyk
,
J. W.
(
2016
).
Parkinson's disease and neurodegeneration: GABA-collapse hypothesis
.
Front. Neurosci.
10
,
269
.
Charreton
,
M.
,
Decourtye
,
A.
,
Henry
,
M.
,
Rodet
,
G.
,
Sandoz
,
J. C.
,
Charnet
,
P.
and
Collet
,
C.
(
2015
).
A locomotor deficit induced by sublethal doses of pyrethroid and neonicotinoid insecticides in the honeybee Apis mellifera
.
PLoS One
10
,
e0144879
.
Chaudhuri
,
A.
,
Johnson
,
R.
,
Rakshit
,
K.
,
Bednářová
,
A.
,
Lackey
,
K.
,
Chakraborty
,
S. S.
,
Krishnan
,
N.
and
Chaudhuri
,
A.
(
2020
).
Exposure to Spectracide® causes behavioral deficits in Drosophila melanogaster: Insights from locomotor analysis and molecular modeling
.
Chemosphere
248
,
126037
.
Ci
,
S.
,
Ren
,
T.
and
Su
,
Z.
(
2007
).
Modeling the interaction of fipronil-related non-competitive antagonist with the GABA Beta3-receptor
.
J. Mol. Model.
13
,
457
-
464
.
Coulom
,
H.
and
Birman
,
S.
(
2004
).
Chronic exposure to rotenone models sporadic Parkinson's disease in Drosophila melanogaster
.
J. Neurosci.
24
,
10993
-
10998
.
Daisley
,
B. A.
,
Trinder
,
M.
,
McDowell
,
T. W.
,
Welle
,
H.
,
Dube
,
J. S.
,
Ali
,
S. N.
,
Leong
,
H. S.
,
Sumarah
,
M. W.
and
Reid
,
G.
(
2017
).
Neonicotinoid-induced pathogen susceptibility is mitigated by Lactobacillus plantarum immune stimulation in a Drosophila melanogaster model
.
Sci. Rep.
7
,
2703
.
Davidsdottir
,
S.
,
Cronin-Golomb
,
A.
and
Lee
,
A.
(
2005
).
Visual and spatial symptoms in Parkinson's disease
.
Vision Res.
45
,
1285
-
1296
.
Decourtye
,
A.
,
Armengaud
,
C.
,
Renou
,
M.
,
Devillers
,
J.
,
Cluzeau
,
S.
,
Gauthier
,
M.
and
Pham-Delègue
,
M. H.
(
2004
).
Imidacloprid impairs memory and brain metabolism in the honeybee (Apis mellifera L.)
.
Pestic. Biochem. Physiol
78
,
83
-
92
.
Dong
,
C.
,
Hu
,
A.
,
Ni
,
Y.
,
Zuo
,
Y.
and
Li
,
G. H.
(
2013
).
Effects of midazolam, pentobarbital and ketamine on the mRNA expression of ion channels in a model organism Daphnia pulex
.
BMC Anesthesiol
13
,
32
.
Eghorn
,
L. F.
,
Hoestgaard-Jensen
,
K.
,
Kongstad
,
K. T.
,
Bay
,
T.
,
Higgins
,
D.
,
Frølund
,
B.
and
Wellendorph
,
P.
(
2014
).
Positive allosteric modulation of the GHB high-affinity binding site by the GABAA receptor modulator monastrol and the flavonoid catechin
.
Eur. J. Pharmacol.
740
,
570
-
577
.
Elbaz
,
A.
,
Clavel
,
J.
,
Rathouz
,
P. J.
,
Moisan
,
F.
,
Galanaud
,
J. P.
,
Delemotte
,
B.
,
Alpérovitch
,
A.
and
Tzourio
,
C.
(
2009
).
Professional exposure to pesticides and Parkinson disease
.
Ann. Neurol
66
,
494
-
504
.
Enell
,
L.
,
Hamasaka
,
Y.
,
Kolodziejczyk
,
A.
and
Nässel
,
D. R.
(
2007
).
γ-Aminobutyric acid (GABA) signaling components in Drosophila: Immunocytochemical localization of GABAB receptors in relation to the GABAA receptor subunit RDL and a vesicular GABA transporter
.
J. Comp. Neurol
505
,
18
-
31
.
Filomeni
,
G.
,
Graziani
,
I.
,
de Zio
,
D.
,
Dini
,
L.
,
Centonze
,
D.
,
Rotilio
,
G.
and
Ciriolo
,
M. R.
(
2012
).
Neuroprotection of kaempferol by autophagy in models of rotenone-mediated acute toxicity: Possible implications for Parkinson's disease
.
Neurobiol. Aging
33
,
767
-
785
.
Glass
,
G. V.
(
1966
).
Testing Homogeneity of Variances
.
Am. Educ. Res. J.
3
,
187
-
190
.
Godinho
,
A. F.
,
Chagas
,
A. C. S.
,
Carvalho
,
C. C.
,
Horta
,
D. F.
,
De Fraia
,
D.
,
Anselmo
,
F.
,
Chaguri
,
J. L.
and
Faria
,
C. A.
(
2016
).
Memory impairment due to fipronil pesticide exposure occurs at the GABAA receptor level, in rats
.
Physiol. Behav
165
,
28
-
34
.
Goulson
,
D.
,
Nicholls
,
E.
,
Botías
,
C.
and
Rotheray
,
E. L.
(
2015
).
Bee declines driven by combined stress from parasites, pesticides, and lack of flowers
.
Science (80–.)
347
,
1255957
.
Grundmann
,
O.
,
Nakajima
,
J. I.
,
Kamata
,
K.
,
Seo
,
S.
and
Butterweck
,
V.
(
2009
).
Kaempferol from the leaves of Apocynum venetum possesses anxiolytic activities in the elevated plus maze test in mice
.
Phytomedicine
16
,
295
-
302
.
Grünz
,
G.
,
Haas
,
K.
,
Soukup
,
S.
,
Klingenspor
,
M.
,
Kulling
,
S. E.
,
Daniel
,
H.
and
Spanier
,
B.
(
2012
).
Structural features and bioavailability of four flavonoids and their implications for lifespan-extending and antioxidant actions in C. elegans
.
Mech. Ageing Dev.
133
,
1
-
10
.
Hall
,
B. J.
,
Karim
,
N.
,
Chebib
,
M.
,
Johnston
,
G. A. R.
and
Hanrahan
,
J. R.
(
2014
).
Modulation of ionotropic GABA receptors by 6-methoxyflavanone and 6-methoxyflavone
.
Neurochem. Res
39
,
1068
-
1078
.
Hämäläinen
,
M.
,
Nieminen
,
R.
,
Vuorela
,
P.
,
Heinonen
,
M.
and
Moilanen
,
E.
(
2007
).
Anti-inflammatory effects of flavonoids: Genistein, kaempferol, quercetin, and daidzein inhibit STAT-1 and NF-κB activations, whereas flavone, isorhamnetin, naringenin, and pelargonidin inhibit only NF-κB activation along with their inhibitory effect on iNOS expression and NO production in activated macrophages
.
Mediators Inflamm.
2007
,
1
-
10
.
Hoffman
,
J. I. E
. (
2019
).
Chapter 25 - Analysis of Variance. I. One-Way
. In
Basic Biostatistics for Medical and Biomedical Practitioners
, 2nd edn (ed.
J. I. E.
Hoffman
), pp.
391
-
417
.
Academic Press
.
Hohler
,
A. D.
and
de Leon
,
M. P.
(
2011
).
Festination
. In
Encyclopedia of Clinical Neuropsychology
(ed.
J. S.
Kreutzer
,
J.
DeLuca
and
B.
Caplan
), pp.
1033
-
1034
.
Springer
.
Holder
,
P. J.
,
Jones
,
A.
,
Tyler
,
C. R.
and
Cresswell
,
J. E.
(
2018
).
Fipronil pesticide as a suspect in historical mass mortalities of honey bees
.
Proc. Natl. Acad. Sci. U. S. A
115
,
13033
-
13038
.
Ikeda
,
T.
,
Nagata
,
K.
,
Kono
,
Y.
,
Yeh
,
J. Z.
and
Narahashi
,
T.
(
2004
).
Fipronil modulation of GABAA receptor single-channel currents
.
Pest Manag. Sci
60
,
487
-
492
.
Jan
,
S.
and
Abbas
,
N
. (
2018
).
Chemistry of himalayan phytochemicals
. In
Himalayan Phytochemicals, JMP®, Version v.14.2
, pp.
1989
-
2021
.
Cary, NC
:
SAS Institute Inc
.
Jung
,
U. J.
and
Kim
,
S. R.
(
2018
).
Beneficial effects of flavonoids against Parkinson's disease
.
J. Med. Food
21
,
421
-
432
.
Koller
,
W. C.
(
1984
).
The diagnosis of Parkinson's disease
.
Arch. Intern. Med.
144
,
2146
-
2147
.
Kottler
,
B.
,
Faville
,
R.
,
Bridi
,
J. C.
and
Hirth
,
F.
(
2019
).
Inverse control of turning behavior by Dopamine D1 receptor signaling in columnar and ring neurons of the central complex in Drosophila
.
Curr. Biol
29
,
567
-
577.e6
.
Leeb-Lundberg
,
F.
,
Napias
,
C.
and
Olsen
,
R. W.
(
1981
).
Dihydropicrotoxinin binding sites in mammalian brain: Interaction with convulsant and depressant benzodiazepines
.
Brain Res.
216
,
399
-
408
.
Li
,
S.
and
Pu
,
X. P.
(
2011
).
Neuroprotective effect of kaempferol against a 1–methyl-4-phenyl-1,2,3,6- tetrahydropyridine-induced mouse model of Parkinson's disease
.
Biol. Pharm. Bull.
34
,
1291
-
1296
.
Li
,
Z.
,
Yu
,
T.
,
Chen
,
Y.
,
Heerman
,
M.
,
He
,
J.
,
Huang
,
J.
,
Nie
,
H.
and
Su
,
S.
(
2019
).
Brain transcriptome of honey bees (Apis mellifera) exhibiting impaired olfactory learning induced by a sublethal dose of imidacloprid
.
Pestic. Biochem. Physiol.
156
,
36
-
43
.
Li
,
F.
,
Lindsey
,
J.
,
Marin
,
E. C.
,
Otto
,
N.
,
Dreher
,
M.
,
Dempsey
,
G.
,
Stark
,
I.
,
Bates
,
A. S.
,
Pleijzier
,
M. W.
,
Schlegel
,
P.
et al. 
(
2020
).
The connectome of the adult Drosophila mushroom body provides insights into function
.
Elife
9
,
e62576
.
Liao
,
L. H.
,
Wu
,
W. Y.
,
Berenbaum
,
M. R.
(
2017
).
Impacts of dietary phytochemicals in the presence and absence of pesticides on longevity of honey bees (Apis mellifera)
.
Insects
8
,
22
.
Liao
,
L. H.
,
Pearlstein
,
D. J.
,
Wu
,
W. Y.
,
Kelley
,
A. G.
,
Montag
,
W. M.
,
Hsieh
,
E. M.
and
Berenbaum
,
M. R.
(
2020
).
Increase in longevity and amelioration of pesticide toxicity by natural levels of dietary phytochemicals in the honey bee, Apis mellifera
.
PLoS One
15
,
e0243364
.
Liu
,
Z.
,
Lindemeyer
,
A. K.
,
Liang
,
J.
,
Wallner
,
M.
,
Shao
,
X. M.
,
Shao
,
Y.
,
Tao
,
Y.
and
Olsen
,
R. W.
(
2018
).
Flavonoids isolated from Tibetan medicines, binding to GABAA receptor and the anticonvulsant activity
.
Phytomedicine
50
,
1
-
7
.
Maher
,
P.
(
2017
).
Protective effects of fisetin and other berry flavonoids in Parkinson's disease
.
Food Funct
8
,
3033
-
3042
.
Maher
,
P.
(
2019
).
The potential of flavonoids for the treatment of neurodegenerative diseases
.
Int. J. Mol. Sci.
20
,
3056
.
Martelli
,
F.
,
Zhongyuan
,
Z.
,
Wang
,
J.
,
Wong
,
C. O.
,
Karagas
,
N. E.
,
Roessner
,
U.
,
Rupasinghe
,
T.
,
Venkatachalam
,
K.
,
Perry
,
T.
,
Bellen
,
H. J.
et al. 
(
2020
).
Low doses of the neonicotinoid insecticide imidacloprid induce ROS triggering neurological and metabolic impairments in Drosophila
.
Proc. Natl. Acad. Sci. U. S. A
117
,
25840
-
25850
.
Martin
,
J. R.
,
Ernst
,
R.
and
Heisenberg
,
M.
(
1998
).
Mushroom bodies suppress locomotor activity in Drosophila melanogaster
.
Learn. Mem
5
,
179
-
191
.
Martin
,
J. R.
,
Raabe
,
T.
and
Heisenberg
,
M.
(
1999
).
Central complex substructures are required for the maintenance of locomotor activity in Drosophila melanogaster
.
J. Comp. Physiol. A Sensory, Neural, Behav. Physiol
185
,
277
-
288
.
Massey
,
F. J.
(
1951
).
The Kolmogorov-Smirnov Test for Goodness of Fit
.
J. Am. Stat. Assoc.
46
,
68
-
78
.
Mitrocotsa
,
D.
,
Mitaku
,
S.
,
Axarlis
,
S.
,
Harvala
,
C.
and
Malamas
,
M.
(
2000
).
Evaluation of the antiviral activity of kaempferol and its glycosides against human cytomegalovirus
.
Planta Med.
66
,
377
-
379
.
Mohamed
,
F.
,
Senarathna
,
L.
,
Percy
,
A.
,
Abeyewardene
,
M.
,
Eaglesham
,
G.
,
Cheng
,
R.
,
Azher
,
S.
,
Hittarage
,
A.
,
Dissanayake
,
W.
,
Sheriff
,
M. H. R.
et al. 
(
2004
).
Acute human self-poisoning with the N-phenylpyrazole insecticide fipronil - A GABAA-gated chloride channel blocker
.
J. Toxicol. Clin. Toxicol.
42
,
955
-
963
.
Nagoshi
,
E.
(
2018
).
Drosophila models of sporadic parkinson's disease
.
Int. J. Mol. Sci.
19
,
3343
.
Narayan
,
S.
,
Liew
,
Z.
,
Bronstein
,
J. M.
and
Ritz
,
B.
(
2017
).
Occupational pesticide use and Parkinson's disease in the Parkinson Environment Gene (PEG) study
.
Environ. Int
107
,
266
-
273
.
Neuser
,
K.
,
Triphan
,
T.
,
Mronz
,
M.
,
Poeck
,
B.
and
Strauss
,
R.
(
2008
).
Analysis of a spatial orientation memory in Drosophila
.
Nature
453
,
1244
-
1247
.
Nonnekes
,
J.
,
Giladi
,
N.
,
Guha
,
A.
,
Fietzek
,
U. M.
,
Bloem
,
B. R.
and
Růžička
,
E.
(
2019
).
Gait festination in parkinsonism: introduction of two phenotypes
.
J. Neurol.
266
,
426
-
430
.
Pallanck
,
L. J.
and
Whitworth
,
A. J.
(
2005
).
Drosophila Models of Parkinson Disease
.
Movement Disorders
(Second Edition) (ed. M. S. LeDoux), pp. 345-354.
Academic Press
.
Pallis
,
C. A.
(
1971
).
Parkinsonism: natural history and clinical features
.
Br. Med. J
3
,
683
-
690
.
Park
,
J. H.
,
Jung
,
J. W.
,
Ahn
,
Y. J.
and
Kwon
,
H. W.
(
2012
).
Neuroprotective properties of phytochemicals against paraquat-induced oxidative stress and neurotoxicity in Drosophila melanogaster
.
Pestic. Biochem. Physiol.
104
,
118
-
125
.
Park
,
J. H.
,
Park
,
Y. S.
and
Koh
,
H. C.
(
2016
).
Progressive loss of nigrostriatal dopaminergic neurons induced by inflammatory responses to fipronil
.
Toxicol. Lett.
258
,
36
-
45
.
Pettersen
,
E. F.
,
Goddard
,
T. D.
,
Huang
,
C. C.
,
Couch
,
G. S.
,
Greenblatt
,
D. M.
,
Meng
,
E. C.
and
Ferrin
,
T. E.
(
2004
).
UCSF Chimera - A visualization system for exploratory research and analysis
.
J. Comput. Chem.
25
,
1605
-
1612
.
Poeck
,
B.
,
Triphan
,
T.
,
Neuser
,
K.
and
Strauss
,
R.
(
2008
).
Locomotor control by the central complex in Drosophila - An analysis of the tay bridge mutant
.
Dev. Neurobiol.
68
,
1046
-
1058
.
Rahul
,
N.
,
Jyoti
,
F.
,
and Siddique
,
S.
and
H
,
Y.
(
2020
).
Effect of kaempferol on the transgenic Drosophila model of Parkinson's disease
.
Sci. Rep.
10
,
13793
.
Ratner
,
M. H.
,
Farb
,
D. H.
,
Ozer
,
J.
,
Feldman
,
R. G.
and
Durso
,
R.
(
2014
).
Younger age at onset of sporadic Parkinson's disease among subjects occupationally exposed to metals and pesticides
.
Interdiscip. Toxicol.
7
,
123
-
133
.
Riveros
,
A. J.
and
Gronenberg
,
W.
(
2022
).
The flavonoid rutin protects the bumble bee Bombus impatiens against cognitive impairment by imidacloprid and fipronil
.
J. Exp. Biol
225
,
jeb244526
.
Rousseeuw
,
P. J.
and
van Zomeren
,
B. C.
(
1990
).
Unmasking multivariate outliers and leverage points
.
Null
85
,
633
-
639
.
Sánchez-Bayo
,
F.
and
Wyckhuys
,
K. A. G.
(
2019
).
Worldwide decline of the entomofauna: A review of its drivers
.
Biol. Conserv.
232
,
8
-
27
.
Schindelin
,
J.
,
Arganda-Carreras
,
I.
,
Frise
,
E.
,
Kaynig
,
V.
,
Longair
,
M.
,
Pietzsch
,
T.
,
Preibisch
,
S.
,
Rueden
,
C.
,
Saalfeld
,
S.
,
Schmid
,
B.
et al. 
(
2012
).
Fiji: An open-source platform for biological-image analysis
.
Nat. Methods
9
,
676
-
682
.
Sigel
,
E.
and
Baur
,
R.
(
1988
).
Allosteric modulation by benzodiazepine receptor ligands of the GABA(A) receptor channel expressed in Xenopus oocytes
.
J. Neurosci.
8
,
289
-
295
.
Soibam
,
B.
,
Mann
,
M.
,
Liu
,
L.
,
Tran
,
J.
,
Lobaina
,
M.
,
Kang
,
Y. Y.
,
Gunaratne
,
G. H.
,
Pletcher
,
S.
and
Roman
,
G.
(
2012
).
Open-field arena boundary is a primary object of exploration for Drosophila
.
Brain Behav.
2
,
97
-
108
.
Stehr
,
C. M.
,
Linbo
,
T. L.
,
Incardona
,
J. P.
and
Scholz
,
N. L.
(
2006
).
The developmental neurotoxicity of fipronil: Notochord degeneration and locomotor defects in zebrafish embryos and larvae
.
Toxicol. Sci.
92
,
270
-
278
.
Strauss
,
R.
and
Heisenberg
,
M.
(
1993
).
A higher control center of locomotor behavior in the Drosophila brain
.
J. Neurosci
13
,
1852
-
1861
.
Strauß
,
R.
and
Heisenberg
,
M.
(
1990
).
Coordination of legs during straight walking and turning in Drosophila melanogaster
.
J. Comp. Physiol. A
167
,
403
-
412
.
Student
(
1908
).
The probable error of a mean
.
Biometrika
6
,
1
-
25
.
Suzuki
,
T.
,
Hirai
,
A.
,
Khidkhan
,
K.
,
Nimako
,
C.
,
Ichise
,
T.
,
Takeda
,
K.
,
Mizukawa
,
H.
,
Nakayama
,
S. M. M.
,
Nomiyama
,
K.
,
Hoshi
,
N.
et al. 
(
2021
).
The effects of fipronil on emotional and cognitive behaviors in mammals
.
Pestic. Biochem. Physiol.
175
,
104847
.
Tallarida
,
R. J.
and
Murray
,
R. B
. (
1987
).
Chi-square test
. In
Manual of Pharmacologic Calculations: With Computer Programs
(ed.
R. J.
Tallarida
and
R. B.
Murray
), pp.
140
-
142
.
Springer New York
.
Tasman
,
K.
,
Hidalgo
,
S.
,
Zhu
,
B.
,
Rands
,
S. A.
and
Hodge
,
J. J. L.
(
2021a
).
Neonicotinoids disrupt memory, circadian behaviour and sleep
.
Sci. Rep
11
,
2061
.
Tasman
,
K.
,
Rands
,
S. A.
and
Hodge
,
J. J. L.
(
2021b
).
The power of Drosophila melanogaster for modeling neonicotinoid effects on pollinators and identifying novel mechanisms
.
Front. Physiol
12
,
659440
.
Trott
,
A. J. O. O.
and
Schroer
,
A.
(
2010
).
AutoDock vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading
.
J. Comput. Chem
31
,
455
-
461
.
Webb
,
B.
and
Sali
,
A.
(
2016
).
Comparative protein structure modeling using MODELLER
.
Curr. Protoc. Bioinforma
54
,
5.6.1
-
5.6.37
.
Weir
,
P. T.
and
Dickinson
,
M. H.
(
2015
).
Functional divisions for visual processing in the central brain of flying Drosophila
.
Proc. Natl. Acad. Sci. U. S. A
112
,
E5523
-
E5532
.
White
,
K. E.
,
Humphrey
,
D. M.
and
Hirth
,
F.
(
2010
).
The dopaminergic system in the aging brain of Drosophila
.
Front. Neurosci.
4
,
205
.
Whitman
,
J. G.
and
Amrein
,
R.
(
1995
).
Pharmacology of flumazenil
.
Acta Anaesthesiol. Scand
108
,
3
-
14
.
Williams
,
R. J.
,
Spencer
,
J. P. E.
and
Rice-Evans
,
C.
(
2004
).
Flavonoids: Antioxidants or signalling molecules?
Free Radic. Biol. Med.
36
,
838
-
849
.
Williamson
,
S. M.
,
Willis
,
S. J.
and
Wright
,
G. A.
(
2014
).
Exposure to neonicotinoids influences the motor function of adult worker honeybees
.
Ecotoxicology
23
,
1409
-
1418
.
Wong
,
M. J.
,
Liao
,
L. H.
and
Berenbaum
,
M. R.
(
2018
).
Biphasic concentration-dependent interaction between imidacloprid and dietary phytochemicals in honey bees (Apis mellifera)
.
PLoS One
13
,
e0206625
.
Woolson
,
R. F
. (
2008
).
Wilcoxon Signed-Rank Test.. In: Encyclopedia of Biostatistics. Vol 6 (ed. P. Armitage and T. Colton), pp. 4739-4740, England: John Wiley & Sons
.
Zaluski
,
R.
,
Kadri
,
S. M.
,
Alonso
,
D. P.
,
Martins Ribolla
,
P. E.
and
de Oliveira Orsi
,
R.
(
2015
).
Fipronil promotes motor and behavioral changes in honey bees (Apis mellifera) and affects the development of colonies exposed to sublethal doses
.
Environ. Toxicol. Chem.
34
,
1062
-
1069
.
Zhao
,
X.
,
Salgado
,
V. L.
,
Yeh
,
J. Z.
and
Narahashi
,
T.
(
2003
).
Differential actions of fipronil and dieldrin insecticides on GABA-gated chloride channels in cockroach neurons
.
J. Pharmacol. Exp. Ther.
306
,
914
-
924
.
Zheng
,
N.
,
Cheng
,
J.
,
Zhang
,
W.
,
Li
,
W.
,
Shao
,
X.
,
Xu
,
Z.
,
Xu
,
X.
and
Li
,
Z.
(
2014
).
Binding Difference of Fipronil with GABAARs in Fruitfly and Zebrafish: Insights from Homology Modeling, Docking, and Molecular Dynamics Simulation Studies
.
J. Agric. Food Chem.
62
,
10646
-
10653
.

Competing interests

Part of the results presented here are included in Patent Cooperation Treaty (PCT) application W02021046388A1.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.