# Modeling the cellular mechanisms and olfactory input underlying the triphasic response of moth pheromone-sensitive projection neurons

Yuqiao Gu

► To cite this version:

Yuqiao Gu. Modeling the cellular mechanisms and olfactory input underlying the triphasic response of moth pheromone-sensitive projection neurons. PLoS ONE, 2015, 10 (5), pp.1-22. 10.1371/journal.pone.0126305. [https://doi.org/10.1371/journal.pone.0126305](https://doi.org/10.1371/journal.pone.0126305). This is an open access article distributed under a Creative Commons Attribution 4.0 International License.

HAL Id: hal-02629566

https://hal.inrae.fr/hal-02629566

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Modeling the Cellular Mechanisms and Olfactory Input Underlying the Triphasic Response of Moth Pheromone-Sensitive Projection Neurons

Yuqiao Gu¹,²*

1 UMR 1392 iEES-Paris, National Institute for Agricultural Research, Versailles, France, 2 Cognitive Neuroscience, International School for Advanced Studies, Trieste, Italy

* alice63_gu@hotmail.com

## Abstract

In the antennal lobe of the noctuid moth *Agrotis ipsilon*, most pheromone-sensitive projection neurons (PNs) exhibit a triphasic firing pattern of excitation (E₁)-inhibition (I)-excitation (E₂) in response to a pulse of the sex pheromone. To understand the mechanisms underlying this stereotypical discharge, we developed a biophysical model of a PN receiving inputs from olfactory receptor neurons (ORNs) via nicotinic cholinergic synapses. The ORN is modeled as an inhomogeneous Poisson process whose firing rate is a function of time and is fitted to extracellular data recorded in response to pheromone stimulations at various concentrations and durations. The PN model is based on the Hodgkin-Huxley formalism with realistic ionic currents whose parameters were derived from previous studies. Simulations revealed that the inhibitory phase I can be produced by a SK current (Ca²⁺-gated small conductance K⁺ current) and that the excitatory phase E₂ can result from the long-lasting response of the ORNs. Parameter analysis further revealed that the ending time of E₁ depends on some parameters of SK, Ca²⁺, nACh and Na⁺ currents; I duration mainly depends on the time constant of intracellular Ca²⁺ dynamics, conductance of Ca²⁺ currents and some parameters of nACh currents; The mean firing frequency of E₁ and E₂ depends differentially on the interaction of various currents. Thus it is likely that the interplay between PN intrinsic currents and feedforward synaptic currents are sufficient to generate the triphasic firing patterns observed in the noctuid moth *A. ipsilon*.

## Introduction

Odor coding by the olfactory system has been studied by various experimental and modeling approaches. Natural odor stimuli can be characterized not only by their molecular features but also by properties such as concentration, spatial and temporal change of chemical components. Behavioral experiments on vertebrates [1], terrestrial [2–5] and aquatic invertebrates [6–7] showed that the physical characteristics of odor stimuli condition the behavioral

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Gu Y (2015) Modeling the Cellular Mechanisms and Olfactory Input Underlying the Triphasic Response of Moth Pheromone-Sensitive Projection Neurons. PLoS ONE 10(5): e0126305. doi:10.1371/journal.pone.0126305

**Academic Editor:** Melissa J. Coleman, Claremont Colleges, UNITED STATES

**Received:** August 10, 2014

**Accepted:** March 31, 2015

**Published:** May 11, 2015

**Copyright:** © 2015 Yuqiao Gu. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was funded by the Agence Nationale de la Recherche within the French-British ANR BBSRC SysBio 006 01 "Pherosys" initiative, the European FP7-ICT 2007 STREP Bio-ICT convergence "Neurochem" and the state program "Investissements d'avenir" managed by the Agence Nationale de la Recherche (grant ANR-10-BINF-05 "Pherotaxis"). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
**Competing Interests:** The authors have declared that no competing interests exist.

response to an odorant. In moths, for example, intermittent and continuous stimulation with the same odor (the sex pheromone) evokes two distinct flight behaviors of upwind zig-zagging flight towards the odor source and cast/counter turn across the wind line, respectively [8]. The stimulation features are encoded and analyzed by individual neurons and neural networks of the olfactory system involving the antennae, the antennal lobes (ALs), the mushroom bodies (MBs) and the lateral horn in insects. Odorants are first detected and encoded by different types of olfactory receptor neurons (ORNs) situated in the antenna. Features of odorant stimuli are further analyzed in the AL, the first-order processing center. ORNs of the same type project to the same glomerulus [9] where they establish synaptic connections with multiglomerular local neurons (LNs), intrinsic to the AL, and uniglomerular projection neurons (PNs) [9]. The stimulation features/parameters influence the spatial and temporal activity patterns throughout the glomerular array and the response characteristics of individual PNs [10]. However, the neural basis of the differential response to the different physical odor stimulation features of the same odorant is poorly known. Neurons at different processing stages show different response patterns in response to the same stimulus [11–13]. Recordings in mitral/tufted cells in vertebrates and PNs in insects revealed that the odor-evoked responses of these second-order neurons are generally complex, consisting of both depolarizing and hyperpolarizing phases [12–19]. Remarkably, the temporal patterns of spike activity observed in some vertebrate mitral/tufted cells and insect PNs are very similar [12–20], suggesting common principles of cellular and/or synaptic mechanisms. In the macroglomerular complex (MGC) of the moth *A. ipsilon*, i.e. the specialist system processing pheromone information in the insect AL, a large majority of pheromone-sensitive PNs exhibit a triphasic firing pattern when the antenna is stimulated with pulses of the sex pheromone [12, 21]. Patch-clamp experiments revealed that several types of Na⁺, Ca²⁺, and K⁺ ionic currents are expressed in PNs [22–25] suggesting that they may play roles in shaping the activity patterns of PNs. Especially, it was recently found that SK channels are expressed in AL PNs both in *Drosophila* [24] and in *Agrotis* [25].

A number of biophysical models of neuron and network have been developed to investigate the cellular, synaptic, network structure and dynamical mechanisms underlying PN firing patterns and odor coding in the MGC or the AL of insects. First, a simplified Hodgkin–Huxley (HH) type neuron model and a neural network model with disinhibition mechanism were developed to simulate the low-frequency (< 10 Hz) background activity and the high-frequency (> 100 Hz) bursting capacity of pheromone-sensitive PNs in moth MGC [26–27]. In the absence of stimulation, the modeled PN is inhibited by a LN (denoted LN2), the model exhibits about 3Hz spontaneous oscillations. During odor stimulation, LN2 is inhibited by another LN (denoted LN1), the PN is released from inhibition and exhibits a burst response at frequency higher than 100 Hz. *I*<sub>Ca</sub> and *I*<sub>K(Ca)</sub> are responsible for burst and quiescent period generation whereas *I*<sub>A</sub> reduces the firing frequency. However, in this model the activation and inactivation of the ionic channels of *I*<sub>Na</sub>, *I*<sub>K</sub>, *I*<sub>A</sub>, *I*<sub>Ca</sub> and *I*<sub>K(Ca)</sub> are simplified and not biophysically realistic. Second, to simulate the temporal activity patterns induced by odor stimuli in the locust AL, neural networks with randomly connected neurons based on HH type models of PNs and LNs were developed [28–29]. In these models, *I*<sub>Ca</sub> and *I*<sub>K(Ca)</sub> for slow patterning generation were in non-spiking LNs, but not in PNs. The temporal patterns of PNs were generated through strong GABA<sub>B</sub>-mediated slow inhibition. Third, in another small neural network model of AL consisting of identical PNs and LNs of HH type, the *I*<sub>Na</sub>, *I*<sub>K</sub>, *I*<sub>Ca</sub> and *I*<sub>K(Ca)</sub> were located in both PNs and LNs [30]. It was found that *I*<sub>Ca</sub> and *I*<sub>K(Ca)</sub> in PNs are sufficient to account for the slow patterning. The authors showed that the major effect of network inhibition is to redistribute the action potentials of the PNs from bursting to one action potential per cycle of oscillation. Fourth, based on morphological and physiological data from glomerular circuitry of insect AL
and using models of PNs and LNs developed in [28–29], a cross-scale neurodynamical model of the AL was developed [31]. This model demonstrates the effects of connectivity and complex dynamics in amplifying weak odor signals, in discriminating signals, and in detecting odor similarity, difference and specialty. Simulation results also showed that the spatiotemporal patterns of the odor information emerging in the glomeruli of the AL rely on the glomerular morphology, the connectivity and the complex dynamics of the AL circuits. Fifth, a model of the MGC in the moth *Manduca sexta* with HH type neuronal models and two types of inhibitory LNs, LNs-IIa, and LNs-IIb was proposed [32]. It was shown that synaptic inhibition, intrinsic currents $I_A$ and $I_{SK}$ in PNs can account for the first and second inhibitory phases and contribute to a rapid encoding of pheromone information. Sixth, recently, using a model of AL with PNs and LNs developed in [28–29], the relationship between a structural property of a network—its colorings, $Ca^{2+}$ dynamics and the spatiotemporal activity and synchronization properties of PNs were explored [33]. Seventh, an inhibitory neural network model of MGC, which was quantitatively reduced from a HH conductance-based model to a mean field one, was recently developed [34]. It was analytically shown that the network's ability to operate on signal amplitudes across several orders of magnitude is optimal when a disinhibitory model is close to losing stability and the network dynamics are close to bifurcation.

However, most of the previous modeling work oversimplified the ORN inputs as an input current to PN. Moreover the ionic currents and parameters of the PN and LN models were not taken from insect neurons. In this work, the ORN was modeled from the experimental data recorded in the noctuid moth *A. ipsilon* from our lab [12]. In light of the availability of patch-clamp data of some ionic currents in PNs or other types of insect neurons [22, 36–38], we developed a biophysical model of PN. In a recent work, using a similar PN model with SK currents we reproduced the $E_1$ phase and I phase [21]. In the present paper, we describe the PN model and its parameters in detail and we consider a realistic ORN input modeled from experimental data. Based on the convergence rate in the moth pheromone system [35], we connected 100 ORNs to 1 PN by fast nicotinic cholinergic synapses to form a simple model of the MGC. Our model was built based on three types of experimental data: intracellular, extracellular and patch-clamp data recorded from ORNs, PNs and other neurons in insects obtained from our lab and other labs. Because ORNs do not show triphasic patterns we simply modeled each ORN firing by an inhomogeneous Poisson process. The firing frequency of the ORN model is a function of time and is fitted to the extracellular recorded data in response to pheromone stimulations varying in concentration and duration [12]. In order to better understand the cellular and synaptic mechanisms underlying the triphasic response patterns of PNs, we made a biophysical PN model taking into account the nicotinic cholinergic currents resulting from ORN synapses onto PNs and various intrinsic ionic currents found in PNs. The parameters in the voltage-dependent steady state and time-dependent functions were fitted to patch-clamp data [22, 36]. When no data were available on PN currents we utilized data from other neuron types in insects [37–38] or even vertebrates [39]. We hypothesize that the multiphasic firing patterns of PNs may be generated by the ionic currents in PNs and ORN inputs, the cholinergic synaptic currents from ORNs to PN may affect the PN response characteristics. Using this model we reproduced the recorded triphasic response patterns of PNs. Then, we investigated the ionic current mechanisms underlying these patterns. We further performed thorough analysis on how the response characteristics change with stimulation parameters and how the ORN inputs, the intrinsic and synaptic currents affect the response characteristics. In addition, we also reconstructed a model of LN and explored possible influences of LNs on the PN response characteristics through $GABA_A$- and $GABA_B$-mediated inhibition. Finally, we draw some conclusions based on our modeling study.
# Results

## Reproducing the triphasic pattern and frequency of PN responses

To understand the cellular mechanisms underlying the triphasic firing patterns (E₁/I/E₂) of MGC neurons in A. ipsilon in response to pheromone stimuli, we developed a simple biophysical MGC network model. This model (see Methods) consists of 100 Poisson ORNs connected to one PN through cholinergic synapses. The outline of the model is shown in S1 Fig and its parameter values are given in Tables 1–4. Using this model we reproduced the triphasic PN response pattern to high concentration pheromone stimuli. Results are shown in Fig 1A–1D. In the simulations, the stimulation onset is 5000ms and the ORN response latency is 140ms. The pheromone stimulus duration and dose are 500ms and 10ng respectively. Since the parameter values of various intrinsic currents were taken from different types of neurons, in order to produce the firing pattern shown in Fig 1A–1D we have modified the values of some parameters from the experimental data. The modified values are also shown in Tables 2–4 (denoted by modified value). Comparing Fig 1A and 1C shows that the spontaneous frequency of PN is higher than that of ORN; the E₁ phase in PN corresponds to the initial response of ORNs where their firing frequency is the highest; the E₂ phase in PN corresponds to the late response of ORNs where their firing rate is lower. Since the PN receives convergent inputs from 100 ORNs, the frequency of the spontaneous activity and those of E₁ and E₂ phases are higher in PN than in ORN. These results agree with the experimental findings in [12–13]. Fig 1B and 1D indicate that the I phase corresponds to the falling phase of intracellular Ca²⁺ of PN (shown by the green rectangle in Fig 1B and 1D). In order to see how the intrinsic current Iₐ affects the PN firing pattern in Fig 1E and 1F, we turned off Iₐ. The frequency of the PN spontaneous activity and that of the second excitatory phase E₂ in Fig 1E and 1F are reduced compared with Fig 1B and 1C. By contrast, in Fig 1G and 1H we reduced the half-activation voltage V₀.₅act of Iₐ from -32.7 to -40.0 mV. The frequency of the second excitatory phase E₂ was significantly increased, whereas the PN spontaneous activity was further reduced. This means that the A current affects the PN firing frequency of various phases in a parameter-dependent way. We

<table><caption>Table 1. Parameter values of ORN model fitted to extracellular recorded data.</caption><thead><tr><th>Dose (ng)</th><th>Period (ms)</th><th>f<sub>sp</sub> (Hz)</th><th>f<sub>pe</sub> (Hz)</th><th>f<sub>pl</sub> (Hz)</th><th>T<sub>lat</sub> (ms)</th><th>T<sub>d2pe</sub> (ms)</th><th>T<sub>pl</sub> (ms)</th><th>τ<sub>rise</sub> (ms)</th><th>τ<sub>fall1</sub> (ms)</th><th>τ<sub>fall2</sub> (s)</th><th>τ<sub>fall3</sub> (s)</th><th>q</th></tr></thead><tbody><tr><td>0.1</td><td>200</td><td>1.5</td><td>16</td><td>—</td><td>250</td><td>150</td><td>—</td><td>180</td><td>130</td><td>20</td><td>—</td><td>0.9</td></tr><tr><td>1.0</td><td>200</td><td>1.5</td><td>35</td><td>—</td><td>250</td><td>115</td><td>—</td><td>128.6</td><td>170</td><td>10</td><td>—</td><td>0.9</td></tr><tr><td>10.0</td><td>200</td><td>1.5</td><td>154</td><td>—</td><td>150</td><td>115</td><td>—</td><td>155</td><td>115</td><td>5</td><td>—</td><td>0.9</td></tr><tr><td>10.0</td><td>500</td><td>1.5</td><td>125</td><td>30</td><td>140</td><td>160</td><td>330</td><td>150</td><td>40</td><td>0.2</td><td>10.5</td><td>0.72</td></tr><tr><td>10.0</td><td>1000</td><td>1.5</td><td>130</td><td>30</td><td>170</td><td>110</td><td>870</td><td>140</td><td>70</td><td>0.3</td><td>11.791</td><td>0.72</td></tr></tbody></table>

doi:10.1371/journal.pone.0126305.t001

<table><caption>Table 2. Parameter values of Ca dynamics and passive parameters of the PN model.</caption><thead><tr><th></th><th>Name</th><th>Value</th><th>Modified value</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="3">Passive parameters</td><td>C<sub>m</sub>(pF)</td><td>22.9</td><td>—</td><td>[22]</td></tr><tr><td>E<sub>L</sub> (mV)</td><td>-61.4</td><td>—</td><td>[22]</td></tr><tr><td>g<sub>L</sub> (nS)</td><td>11.16<sup>a</sup></td><td>—</td><td>[22]</td></tr><tr><td rowspan="3">Parameters of Ca dynamics</td><td>f<sub>Ca</sub></td><td>1.6</td><td>1.7</td><td>[39]</td></tr><tr><td>τ<sub>Ca</sub> (ms)</td><td>656</td><td>2000</td><td>[39]</td></tr><tr><td>Ca<sub>∞</sub>(nM)</td><td>113.0</td><td>—</td><td>[39]</td></tr></tbody></table>

doi:10.1371/journal.pone.0126305.t002

<footnote><sup>a</sup>Calculated by g<sub>L</sub> = 1/R<sub>M</sub> = 1/(89.6 MΩ) = 11.16 nS.</footnote>
<table><caption>Table 3. Parameter values of the ionic currents of the PN model given or calculated from data.</caption><tbody><tr><td>Steady-state functions of I<sub>Na</sub></td><td>g&#773;<sub>Na</sub> (nS)</td><td>V<sub>0.5act</sub> (mV)</td><td>s<sub>m</sub></td><td>V<sub>0.5inact</sub> (mV)</td><td>s<sub>h</sub></td><td>E<sub>Na</sub>(mV)</td><td>Ref.</td></tr><tr><td></td><td>206 modified2500</td><td>-25.8</td><td>9.32</td><td>-41.1 modified43.0</td><td>9.75</td><td>+48.2(RP)+47.9(EP)</td><td>[38]</td></tr><tr><td rowspan="3">Time constant functions of I<sub>Na</sub></td><td>a<sub>tm,up</sub></td><td>V<sub>tm,0.5up</sub> (mV)</td><td>s<sub>tm,up</sub></td><td>a<sub>tm,dn</sub></td><td>V<sub>tm,0.5dn</sub></td><td>s<sub>tm,dn</sub></td><td>Ref.</td></tr><tr><td>0.5</td><td>-30</td><td>3.7</td><td>0.5</td><td>-15</td><td>13.7</td><td>Fitted</td></tr><tr><td>a<sub>th,up</sub></td><td>V<sub>th,0.5up</sub></td><td>s<sub>th,up</sub></td><td>a<sub>th,dn</sub></td><td>V<sub>th,0.5dn</sub></td><td>s<sub>th,dn</sub></td><td>to</td></tr><tr><td></td><td>2.1</td><td>-55</td><td>5</td><td>0.7</td><td>-10</td><td>11</td><td>[38]</td></tr><tr><td>Steady-state functions of I<sub>Ca</sub></td><td>g&#773;<sub>Ca</sub> (nS)</td><td>V<sub>0.5act</sub> (mV)</td><td>s<sub>m</sub></td><td>V<sub>0.5inact</sub> (mV)</td><td>s<sub>h</sub></td><td>E<sub>Ca</sub>(mV)</td><td>Ref.</td></tr><tr><td></td><td>16.1 modified45</td><td>-10.6</td><td>8.5</td><td>-29.6</td><td>8.4</td><td>160</td><td>[22]</td></tr><tr><td rowspan="2">Time constant functions of I<sub>Ca</sub></td><td>a<sub>tm,up</sub></td><td>s<sub>tm,up</sub></td><td>a<sub>tm,dn</sub></td><td>V<sub>tm,dn</sub></td><td>s<sub>tm,dn</sub></td><td>—</td><td>Ref.</td></tr><tr><td>0.046</td><td>20.73</td><td>0.19</td><td>19.8</td><td>10</td><td>—</td><td>[39]</td></tr><tr><td rowspan="2">Steady-state functions of I<sub>sk</sub></td><td>a<sub>msk</sub></td><td>b<sub>msk</sub></td><td>s<sub>msk</sub></td><td>—</td><td>—</td><td>—</td><td>Ref.</td></tr><tr><td>1.120</td><td>2.508</td><td>1000</td><td>—</td><td>—</td><td>—</td><td>[39]</td></tr><tr><td rowspan="2">Steady-state functions of I<sub>Kd</sub></td><td>g&#773;<sub>Kd</sub> (nS)</td><td>V<sub>0.5act</sub> (mV)</td><td>s<sub>m</sub></td><td>E<sub>K</sub> (mV)</td><td>—</td><td>—</td><td>Ref.</td></tr><tr><td>8.17<sup>a</sup> modified700</td><td>-18.5</td><td>22.5 modified 20.0</td><td>-91.6</td><td>—</td><td>—</td><td>[36]</td></tr><tr><td rowspan="3">Time constant functions of I<sub>Kd</sub></td><td>a<sub>tm,up</sub></td><td>V<sub>tm,0.5up</sub></td><td>s<sub>tm,up</sub></td><td>—</td><td>—</td><td>—</td><td>Ref.</td></tr><tr><td>0.125</td><td>-40</td><td>11.0</td><td>—</td><td>—</td><td>—</td><td>Fitted</td></tr><tr><td>a<sub>tm,dn</sub></td><td>V<sub>tm,0.5dn</sub></td><td>s<sub>tm,dn</sub></td><td>—</td><td>—</td><td>—</td><td>to</td></tr><tr><td></td><td>0.15</td><td>25</td><td>45.7</td><td>—</td><td>—</td><td>—</td><td>[37]</td></tr><tr><td>Steady-state functions of I<sub>A</sub></td><td>g&#773;<sub>A</sub> (nS)</td><td>V<sub>0.5act</sub> (mV)</td><td>s<sub>m</sub></td><td>V<sub>0.5inact</sub> (mV)</td><td>s<sub>h</sub></td><td>E<sub>K</sub>(mV)</td><td>Ref.</td></tr><tr><td></td><td>17.35<sup>b</sup> modified 500</td><td>-32.69</td><td>17.5</td><td>-53.3</td><td>7.23</td><td>-91.6</td><td>[36]</td></tr><tr><td rowspan="3">Time constant functions of I<sub>A</sub></td><td>a<sub>tm,up</sub></td><td>V<sub>tm,0.5up</sub></td><td>s<sub>tm,up</sub></td><td>a<sub>tm,dn</sub></td><td>V<sub>tm,0.5dn</sub></td><td>s<sub>tm,dn</sub></td><td>Ref.</td></tr><tr><td>0.5</td><td>-30</td><td>13.7</td><td>0.42</td><td>-15</td><td>46</td><td>Fitted</td></tr><tr><td>a<sub>th,up</sub></td><td>V<sub>th,0.5up</sub></td><td>s<sub>th,up</sub></td><td>a<sub>th,dn</sub></td><td>V<sub>th,0.5dn</sub></td><td>s<sub>th,dn</sub></td><td>to</td></tr><tr><td></td><td>0.04</td><td>-55</td><td>25</td><td>0.045</td><td>40</td><td>55</td><td>[37]</td></tr></tbody></table>

doi:10.1371/journal.pone.0126305.t003

<table><caption>Table 4. Parameter values of nACh synaptic current.</caption><thead><tr><th>Name</th><th>g&#773;<sub>nACh</sub>(µS)</th><th>E<sub>nACh</sub> (mV)</th><th>α (ms<sup>-1</sup>)</th><th>β (ms<sup>-1</sup>)</th><th>A</th><th>t<sub>max</sub> (ms)</th></tr></thead><tbody><tr><td>Value</td><td>0.3</td><td>0.0</td><td>10</td><td>0.2</td><td>0.5</td><td>0.3</td></tr><tr><td>Reference</td><td>[28]</td><td>[28]</td><td>[28]</td><td>[28]</td><td>[28]</td><td>[28]</td></tr><tr><td>Modified value</td><td>0.017</td><td>—</td><td>—</td><td>2.0</td><td>0.8</td><td>—</td></tr></tbody></table>

doi:10.1371/journal.pone.0126305.t004

further checked the influence of $I_{Ca}$ on the firing pattern. We found that decreasing $g_{Ca}$, on the one hand, enhances the firing frequency during spontaneous activity, as well as the $E_1$ and $E_2$ phases; on the other hand, it decreases the duration of the I phase. One of the results is shown in Fig 1I and 1J. In the subsection “Effects of intrinsic PN parameters on the response characteristics”, the influence of PN intrinsic currents is detailed.

In order to further reveal the mechanisms underlying the generation of the $E_1/I/E_2$ pattern, we analyzed the PN depolarizing and repolarizing currents in our simulation (Fig 1A–1D). Employing the low-pass (cut-off frequency: 5Hz) and high-pass (cut-off frequency: 5Hz) Butterworth filters (see Methods), we extracted the slow and fast components of the depolarizing currents $I_{nACh}$, $I_{Na}$, $I_{Ca}$ and the repolarizing currents $I_{SK}$, $I_A$ and $I_{Kd}$ (Fig 2). During $E_1$ and I, the kinetics of the slow component of the repolarizing $Ca^{2+}$-dependent $K^+$ current $I_{SK}$ are similar to those of the depolarizing synaptic current $I_{nACh}$ (Fig 2A). Similarly, during $E_1$, the kinetics of the slow components of $I_A$ and $I_{Kd}$ are similar to those of the $I_{Na}$ and $I_{Ca}$, respectively (Fig 2C and

<footnote>a,bCalculated from Kloppenburg et al.,1999.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Triphasic response pattern reproduced by the MGC model in response to a pheromone stimulus (500 ms, 10ng). In the results shown in A to D, for most parameters in Eqs (4–14), we used the original values from literature given in Tables 2–4 except for some values that were modified (modified values in Tables 2–4). A. Spikes of one ORN (blue lines) and mean firing frequency curve of ORNs (red line). B. Dynamics of the PN membrane potential *V*. C. Spikes of the PN (blue lines) and PN firing frequency (red line). D. Kinetics of intracellular $Ca^{2+}$ concentration of the PN. E. Spikes of the PN (blue lines) and PN firing frequency (red line) ($\bar{g}_A$ = 0 µS, $\bar{g}_{Ca}$ = 0.045 µS). F. Dynamics of the PN membrane potential *V* ($\bar{g}_A$ = 0 µS, $\bar{g}_{Ca}$ = 0.045 µS). G. Spikes of the PN (blue lines) and PN firing frequency (red line) ($\bar{g}_A$ = 0.5 µS, $V_{0.5actA}$ = -40.0 mV, $\bar{g}_{Ca}$ = 0.045 µS). H. Dynamics of the PN membrane potential *V* ($\bar{g}_A$ = 0.5 µS, $V_{0.5actA}$ = -40.0 mV, $\bar{g}_{Ca}$ = 0.045 µS). I. Spikes of the PN (blue lines) and PN firing frequency (red line) ($\bar{g}_A$ = 0.5 µS, $V_{0.5actA}$ = -37.0 mV, $\bar{g}_{Ca}$ = 0.035 µS). J. Dynamics of the PN membrane potential *V* ($\bar{g}_A$ = 0.5 µS, $V_{0.5actA}$ = -37.0 mV, $\bar{g}_{Ca}$ = 0.035 µS).</figcaption></figure>

doi:10.1371/journal.pone.0126305.g001
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Synaptic and intrinsic currents in a PN from the simulation results shown in Fig 1. Left panel: the slow components of the repolarizing currents (black lines) and depolarizing currents (magenta lines); *I*<sub>SK</sub> and-*I*<sub>nACh</sub> (A), *I*<sub>A</sub> and—*I*<sub>Na</sub> (C) and *I*<sub>Kd</sub> and—*I*<sub>Ca</sub> (E). Right panel: Fast components of the repolarizing (black lines) and depolarizing currents (magenta lines); *I*<sub>SK</sub> and *I*<sub>nACh</sub> (B), *I*<sub>A</sub> and *I*<sub>Na</sub> (D) and *I*<sub>Kd</sub> and *I*<sub>Ca</sub> (F). Note that in the left panel we draw the minus values of *I*<sub>nACh</sub>, *I*<sub>Na</sub>, *I*<sub>Ca</sub> for comparing their amplitudes, while in the right panel we draw the values of *I*<sub>nACh</sub>, *I*<sub>Na</sub>, *I*<sub>Ca</sub> directly for comparing their depolarizing and repolarizing effects). The slow components of *E*<sub>1</sub> (from 5140 to 5770ms) and I (from 5770 to 6700ms) are enlarged in the insets in A, C and E; and the fast components of *E*<sub>1</sub> and the period transiting from *E*<sub>1</sub> to I (from 5770 to 6200ms) are enlarged in the insets in B, D and F.</figcaption></figure>

doi:10.1371/journal.pone.0126305.g002

2E). In the beginning of E₁ (from 5140 to 5520ms) the amplitudes of the slow and fast compo-
nents of the depolarizing synaptic current I<sub>nACh</sub> from ORNs are higher than those of the repolar-
izing current I<sub>SK</sub> (insets in Fig 2A and 2B). With depolarization I<sub>SK</sub> increases because of the
accumulation of intracellular Ca²⁺. At 5520ms the amplitude of the slow component of I<sub>SK</sub> ex-
ceeds that of I<sub>nACh</sub> (inset in Fig 2A). From 5520ms onwards I<sub>SK</sub> competes with I<sub>nACh</sub> and slows
down the PN firing. At 5770ms the PN stops spiking and transits to the I phase. During I, I<sub>SK</sub>
flows out against I<sub>nACh</sub> to create a hyperpolarization phase until the intracellular Ca²⁺ concentra-
tion falls. During this hyperpolarization phase, the voltage-gated currents (I<sub>A</sub>, I<sub>Kd</sub>, I<sub>Na</sub> and I<sub>Ca</sub>)
cannot be activated due to the low membrane potential (Fig 2C–2F). Interestingly, both the slow and fast components of $I_A$ have the same kinetics as $I_{Na}$ and the amplitude of $I_A$ is slightly smaller than that of $I_{Na}$, especially during E₁ and E₂ phases (Fig 2C and 2D and the insets). Whenever $I_{Na}$ depolarizes the membrane to make a spike, $I_A$ flows out and repolarizes the membrane, so that the amplitude of each spike is reduced and the firing frequency is increased.

## Effects of the pheromone stimulus parameters on the response characteristics

To see how the stimulation parameters affect the PN response characteristics (such as duration and frequency) whose calculation is described in Methods section, we varied the concentration and duration of the pheromone pulses. For each parameter value, the computer simulation was repeated ten times and the mean and standard deviation of the response characteristics were calculated over the ten trials. The results are shown in Fig 3. The parameter values are the same as those in Fig 1A–1D. The duration of E₁ nearly linearly increases with the stimulus duration and the duration of I for 200ms stimulus duration is slightly lower than that for 500 and 1000ms stimulus duration (Fig 3A). The mean response frequency of E₁ decreases with stimulus duration (this is due to frequency adaptation during E₁) (Fig 3B). E₁ and I durations are independent of stimulation doses (Fig 3C). These relationships between E₁ duration and stimulation duration and dose agree with the experimental findings described in [12]. The mean response frequency of E₁ increases with stimulus concentration while that of E₂ does not (Fig 3D).

## Effects of intrinsic PN parameters on the response characteristics

We examined how the intracellular $Ca^{2+}$ dynamics and intrinsic currents in the PN affect its dynamic response characteristics. To this end we varied the value of each parameter from low to high in a range while keeping the values of other parameters the same as in Fig 1A–1D. At a given value of each parameter the computer simulation was repeated ten times and the mean and standard deviation of each response characteristics were calculated over the ten trials.

**Effects of $Ca^{2+}$ dynamics, $I_{Ca}$ and $I_{SK}$**. The simulation results indicate that the inhibitory phase following E₁ is generated by the slow component of the repolarizing current $I_{SK}$ that exceeds the slow component of the depolarizing synaptic current $I_{nACh}$. $I_{SK}$ depends on $Ca^{2+}$ concentration which in turn depends on $I_{Ca}$. In addition, $I_{Ca}$ also affects the triphasic response pattern (Fig 1I–1J). Therefore we further analyzed the effects of the parameters of the $Ca^{2+}$ dynamics, $I_{Ca}$ and $I_{SK}$ on the PN response characteristics. Fig 4 shows that the duration of the E₁ phase exponentially decreases with $\bar{g}_{SK}$ (A) while the duration of the I phase linearly increases with $\bar{g}_{Ca}$ (C) and $\tau_{Ca}$ (E). The mean frequency of the E₁ phase linearly decreases with $\bar{g}_{Ca}$ (D) and that of the E₂ phase exponentially decreases with $\bar{g}_{Ca}$ (D) and $\tau_{Ca}$ (F). Some parameters for the steady-state activation $m_{\infty}$, the time constant $\tau_m$ of $I_{Ca}$ and for the steady-state function of the gating variable $m_{sk\infty}$ of $I_{SK}$ have also clear influences on the PN response duration, especially the duration of the I phase. These effects are illustrated in S2 Fig. The I duration decreases with $a_{rm,up}$ and $S_{msk}$, while it increases with $S_{rm,up}$ and $a_{msk}$.

**Effects of INa, IA and Ikd.** In this section, by varying the parameter values of $I_{Na}$, $I_{kd}$ and $I_A$, we quantitatively investigated how the PN response characteristics depend on these parameters. The major influences of $I_{Na}$ on the response characteristics are illustrated in Fig 5. The maximal conductance and parameters for steady-state activation and inactivation function of $I_{Na}$ strongly affect the response frequency of E₁ phase: E₁ frequency increases with the half-activation parameter $V_{0.5act}$ (Fig 5D) and decreases with $\bar{g}_{Na}$ (Fig 5B) and the slope factor $S_h$ (Fig 5H); both E₁ and E₂ frequencies increase with the slope factor $S_m$ (Fig 5F). Parameters $V_{0.5act}$, $S_m$ and $S_h$ of $I_{Na}$ have also some influences on E₁ and I duration (Fig 5C, 5E and 5G). Other parameters of $I_{Na}$
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Effects of stimulation parameters on PN response characteristics. Top panel: effects of stimulus duration on E₁ and I duration (A) and mean firing frequency of E₁ and E₂ (B). Bottom panel: effects of stimulus concentration on E₁ and I durations (C) and mean firing frequency of E₁ and E₂ phases (D).</figcaption></figure>

doi:10.1371/journal.pone.0126305.g003

affect one or two PN response characteristics as shown in S3 Fig. *I*<sub>Kd</sub> has also some influences on the PN response characteristics as shown in S4 Fig. As illustrated in Fig 1 the transient potassium current *I*<sub>A</sub> affects the firing frequency of E<sub>1</sub>, E<sub>2</sub>. Here we further investigated the influences of *I*<sub>A</sub>. We found that the mean maximal conductance, the half activation, and the slope factor *S*<sub>m</sub> of *I*<sub>A</sub> have strong effects on PN response frequency. Increasing the maximal mean conductance *g*̄<sub>A</sub> and decreasing the voltage of half activation *V*<sub>0.5act</sub> and the slope factor *S*<sub>m</sub> increase the mean firing frequency of both E<sub>1</sub> and E<sub>2</sub> phases (Fig 6B, 6D and 6F). The increased frequency of spontaneous activity of the PN is also due to the higher conductance and lower voltage of half activation and smaller slope factor *s*<sub>m</sub> of this current (data not shown). In addition, increasing *g*̄<sub>A</sub> has a small effect of decreasing the duration of the I phase at low *g*̄<sub>A</sub> values (Fig 6A); increasing *V*<sub>0.5act</sub> clearly decreases the duration of the I phase and increases the duration of the E<sub>1</sub> phase (Fig 6C). Moreover, decreasing *S*<sub>m</sub> clearly decreases I duration whereas increases E<sub>1</sub> duration. Some parameters in time constant functions of *I*<sub>A</sub> have slight influences on PN response characteristics (data not shown): I duration increases with *a*<sub>rmup</sub>, *a*<sub>rlvdn</sub> and *V*<sub>rm0.5dn</sub> when it is less than -10 mV while decreases with *S*<sub>rmup</sub>, *a*<sub>rlvup</sub> and *V*<sub>rlv0.5dn</sub> when it is less than -40 mV; E<sub>1</sub> duration decreases with *a*<sub>rmup</sub> and E<sub>1</sub> frequency decreases with *V*<sub>rm0.5up</sub>, *S*<sub>rmdn</sub>, *a*<sub>rlvup</sub>, *S*<sub>rlvup</sub> and *a*<sub>rlvdn</sub>; E<sub>2</sub> frequency increases with *S*<sub>rmup</sub> and *a*<sub>rlvdn</sub> when it is less than 0.7.
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Effects of *I*<sub>Ca</sub>, *I*<sub>SK</sub> and dynamics of intracellular Ca<sup>2+</sup> on PN response characteristics. Top panel: effects of the mean maximal conductance g&#773;<sub>SK</sub> of *I*<sub>SK</sub> on E<sub>1</sub> and I durations (A) and mean firing frequency of E<sub>1</sub> and E<sub>2</sub> (B). Intermediate panel: effects of the mean maximal conductance g&#773;<sub>Ca</sub> of *I*<sub>Ca</sub> on E<sub>1</sub> and I duration (C) and mean firing frequency of E<sub>1</sub> and E<sub>2</sub> (D). Bottom panel: effects of time constant τ<sub>Ca</sub> of Ca<sup>2+</sup> dynamics on E<sub>1</sub> and I duration (E) and mean firing frequency of E<sub>1</sub> and E<sub>2</sub> phases (F).</figcaption></figure>

doi:10.1371/journal.pone.0126305.g004

Effect of nACh synaptic parameters on PN firing patterns

In the presence or absence of pheromone stimulation, PN dendrites receive feedforward cho-
linergic synaptic inputs from ORNs through nicotinic receptors. Hence, the nACh synaptic
currents are the stimulation inputs of PNs. We studied how the parameters of I<sub>nACh</sub> affect the
PN response characteristics.
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Major effects of *I*<sub>Na</sub> on PN response characteristics. Effects of *g*̄<sub>Na</sub>, *V*<sub>0.5act</sub>, *S*<sub>m</sub> and *S*<sub>h</sub> on *E*<sub>1</sub> and *I* duration (left panel) and mean firing frequency of *E*<sub>1</sub> and *E*<sub>2</sub> (right panel).</figcaption></figure>

doi:10.1371/journal.pone.0126305.g005
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Major effects of *I*<sub>A</sub> on PN response characteristics. Top panel: effects of *g*̄<sub>A</sub> on E<sub>1</sub> and I duration (A) and mean firing frequency of E<sub>1</sub> and E<sub>2</sub> (B). Intermediate panel: effects of *V*<sub>0.5act</sub> on E<sub>1</sub> and I duration (C) and mean firing frequency of E<sub>1</sub> and E<sub>2</sub> phases (D). Bottom panel: effects of *S*<sub>m</sub> on E<sub>1</sub> and I duration (E) and mean firing frequency of E<sub>1</sub> and E<sub>2</sub> phases (F).</figcaption></figure>

doi:10.1371/journal.pone.0126305.g006

First, the effects of presynaptic ACh transmitter delivered as square pulses of duration *t*<sub>max</sub> and concentration A were investigated (S5 Fig). E<sub>1</sub> duration increases while I duration decreases with Ach pulse duration *t*<sub>max</sub> (S5A Fig). The mean frequency of the E<sub>2</sub> phase *f*̄<sub>E2</sub> significantly increases with *t*<sub>max</sub>, whereas the influence of this parameter on the average frequency of the E<sub>1</sub> phase *f*̄<sub>E1</sub> is not monotonic (S5B Fig). This result is interesting because this parameter has a different effect on the average frequency of the E<sub>1</sub> and E<sub>2</sub> phases. The *f*̄<sub>E1</sub> increases at values of *t*<sub>max</sub> smaller than 0.3ms then decreases with *t*<sub>max</sub>. This is due to the fact that the firing frequency adaptation induced by SK currents takes effect when the E<sub>1</sub> duration increases with
$t_{max}$. The ACh concentration parameter A has a very different effect on the duration of E₁ and I phase and the average frequency of E₁ phase $\bar{f}_{E1}$. The duration of E₁ phase increases and that of I phase decreases with A when A is less than 1, then the duration of E₁ and I phases reaches saturation (S5C Fig). The clearest effect of A is that it can significantly increase the average frequency of E₁ phase $\bar{f}_{E1}$ (S5D Fig).

Second, we investigated the effects of the opening and closing rates of nACh postsynaptic channels on PN response characteristics (S6A–S6D Fig). Increasing the opening rate $\alpha$ slightly increases E₁ duration and decreases I duration (S6A Fig) and significantly increases the mean frequency of E₂ phase (S6B Fig purple line). Increasing the closing rate $\beta$ has exactly opposite effects (S6C and S6D Fig purple line). Interestingly, increasing $\alpha$ and $\beta$ has the same clear effect of increasing the average frequency of the E₁ phase (S6B and S6D Fig red lines). The increasing effect of $\beta$ on $\bar{f}_{E1}$ can be explained by its decreasing effect on E₁ duration. When the value of $\beta$ is small, the E₁ duration is quite long (S6C Fig). Thus the frequency adaptation mediated by SK currents is strong. As a result of frequency adaptation, the mean frequency of E₁ phase is low.

Finally, the effects of the mean maximal conductance $g_{nACh}$ on PN response characteristics were studied. Increasing $g_{nACh}$ clearly increases the duration of the E₁ phase and strongly decreases that of the I phase (S6E Fig) and significantly increases the mean frequency of both phases (S6F Fig).

Exploring possible effects of LNs on PN response patterns

In order to investigate possible network mechanisms, particularly the influences from LNs we have developed a model of type I LNs (LNIs), that generate Na⁺-driven action potentials [22–23]. The mathematical description of the LNI model and its parameter values are given in S1 Text. We have done some exploratory simulations by connecting 80 ORNs and 20 LNI with one PN. The LNIs receive inputs from ORNs and PN through fast nicotinic cholinergic synapses, while the PN receives fast nicotinic cholinergic synaptic inputs from ORNs and fast GABAergic inhibitory inputs from LNIs mediated by GABA<sub>A</sub> receptors or slow GABAergic inhibitory inputs from LNIs mediated by metabotropic GABA<sub>B</sub> receptors. Preliminary results revealed that the synaptic interactions between PN and LNIs affect the synchronization among the PN and LNIs, PN response pattern and PN response characteristics. Synchronization and I duration changes with fast GABAergic inhibition. S7 Fig qualitatively shows how the closing rate $\beta$ of the GABA<sub>A</sub> synaptic currents from LNI to PN affects the response characteristics and synchronization. Decreasing $\beta$ decreases the synchronization of PN and LNIs and also decreases the duration of I phase. As for the influence of slow GABAergic inhibition, we found that in the parameter range of the slow GABAergic inhibition given in S1 Text (Eq S3) the PN maintained the triphasic response pattern (S8A Fig). Decreasing the rate parameter $r_3$ in Eq S3 prolonged the I duration (S8E Fig). In these cases the rising kinetics of G protein concentration is close to that of the intracellular Ca²⁺ concentration in PN (S8C and S8F Fig). The LNIs showed synchronized triphasic response pattern (S8B and S8D Fig). When increasing $r_3$ and decreasing $r_4$ in Eq S3 the rising kinetics of GABA<sub>B</sub> receptor-coupled G protein concentration became faster than that of the intracellular Ca²⁺ concentration in PN. In this case the E₁ phase of PN may become an I phase (data not shown) or terminate early due to the GABAergic inhibition (S8G Fig). The shortened or disappeared E₁ is at odds with the experimental findings in [12] and S9 Fig showing that the E₁ duration lasts approximately as long as the stimulus and increases with the duration of pheromone stimuli. After an I phase, another excitatory phase appeared (S8G Fig) which corresponds to the rising period of intracellular Ca²⁺ concentration (S8H Fig). The second excitatory phase in turn is followed by another short I phase corresponding to the early falling phase of Ca²⁺ concentration. After the short I phase, the late excitatory phase of PN appeared.
Discussion

In this work, using a simple MGC model with ORNs based on the recorded frequency curves under different pheromone stimuli, PN based on patch-clamp data from PN and other neurons and a biophysical model of nACh synapses, we reproduced the triphasic response pattern of PNs at high pheromone stimulation concentrations. We investigated mechanisms generating this triphasic response pattern. Our results show that it can be shaped by intrinsic mechanisms in ORNs and PNs: in our model Ca²⁺-dependent SK current in PN is responsible for the I phase following the E₁ phase and the E₂ phase is due to the long-lasting excitatory response of ORNs. We further investigated how the external stimulation parameters and the parameters of the internal ionic currents in PN and the nACh synaptic currents from ORNs to PN affect the duration of E₁ and I, the firing frequency of E₁ and E₂ and other response characteristics of PN. In our model E₁ duration significantly increases with stimulation duration. The ending time of E₁ clearly depends on parameters: **g̅<sub>SK</sub>**: **g̅<sub>Ca</sub>**: V₀.5act, S_m and S_h of I<sub>Na</sub>; t_max, A, β and **g̅<sub>nACh</sub>** of I<sub>nACh</sub>. The external stimulation parameters have no significant influence on I duration. This implies that I phase is an intrinsic property of the network. Our results revealed that I duration linearly increases with the time constant of intracellular Ca²⁺ (τ<sub>Ca</sub>) and **g̅<sub>Ca</sub>**, decreases with t_max and **g̅<sub>nACh</sub>** of I<sub>nACh</sub>. I duration is also influenced by some other parameters such as S_m, S<sub>tm,up</sub> of I<sub>Na</sub>; a<sub>tm,up</sub>, S<sub>tm,up</sub> of I<sub>Ca</sub>; a<sub>msk</sub> and S<sub>msk</sub> of I<sub>SK</sub>; **g̅<sub>A</sub>** and V₀.5act of I<sub>A</sub>; and the concentration of the pulse of ACh transmitter delivered. The mean firing frequency of E₁ phase **f̃<sub>E1</sub>** increases with stimulation concentration and decreases with stimulation duration. **f̃<sub>E1</sub>** and **f̃<sub>E2</sub>** are also strongly affected by some intrinsic parameters. Both of **f̃<sub>E1</sub>** and **f̃<sub>E2</sub>** increase with S_m of I<sub>Na</sub>; **f̃<sub>E1</sub>** increases with V₀.5act of I<sub>Na</sub>; **g̅<sub>A</sub>**; **g̅<sub>nACh</sub>**, A, α and β of I<sub>nACh</sub>; while decreases with **g̅<sub>Na</sub>**, S_h and V₀.5inact of I<sub>Na</sub>; V₀.5act and S_m of I<sub>Kd</sub>; **g̅<sub>Ca</sub>**; and V₀.5act of I<sub>A</sub>. Besides, **g̅<sub>SK</sub>** and t_max have non monotonic effects on **f̃<sub>E1</sub>**·**f̃<sub>E2</sub>** increases with **g̅<sub>A</sub>**, t_max, A, α and **g̅<sub>nACh</sub>**, while it decreases with **g̅<sub>Ca</sub>**, τ<sub>Ca</sub> and β.

The aim of this study is to show that it is possible to explain the triphasic PN responses with
intrinsic mechanisms only (and realistic parameter values) without denying the possible impli-
cation of extrinsic mechanisms and influences. We have also explored possible network effects
on the PN response characteristics, particularly the influences from LNs. Preliminary results
show that I duration and synchronization between PN and LNs change with GABAergic inhi-
bition. Slow GABAergic inhibition does not affect the triphasic PN response pattern in the pa-
rameter range given in S1 and S2 Texts. When the rising kinetics of GABA<sub>B</sub> receptor-coupled
G protein concentration became faster than that of the intracellular Ca<sup>2+</sup> concentration in PN,
GABA<sub>B</sub> mediated inhibition may change the triphasic pattern and result in a much reduced E<sub>1</sub>
duration which is at odds with the experimental data. This might indicate that in the MGC of
A. ipsilon moths if GABA<sub>B</sub> receptor mediated slow inhibitory synapses exist GABA<sub>B</sub> receptor
may couple to Ca<sup>2+</sup> activated SK channel via G protein as reviewed in [41]. Previously we have
found that applying bicuculline (BIC), an antagonist of GABA<sub>A</sub> and a SK channel blocker, to
A. ipsilon moths abolished the inhibitory phase in all tested neurons [21]. However, applying
picrotoxin (PTX), another GABA<sub>A</sub> antagonist, to A. ipsilon led to different effects. This experi-
mental result together with our modeling results indicate that the Ca<sup>2+</sup>-activated SK channel is
likely responsible for the generation of I phase while interactions between LNs and PNs might
affect the PN response characteristics. These are very preliminary qualitative results showing
the influences of LNs on PN response pattern. The network structure of MGC may also affect
the PN response characteristics. Further experimental and modeling investigations are needed
to study how the intrinsic property of LNs, various synaptic mechanisms and network struc-
ture of MGC affect PN response. Besides, the functional roles of PN response patterns in the
dynamical representation, classification and discrimination of pheromone stimuli and in guid-
ing the moths tracking in turbulent and intermittent pheromone plumes to be elucidated.
In conclusion, our modeling study revealed that $I_{SK}$ and the long-lasting excitatory response of ORNs can be intrinsic mechanisms for the generation of triphasic response patterns of pheromone-sensitive PN. The parameters of $Ca^{2+}$, nACh synaptic, $Na^+$ and A currents have strong influence on the response patterns and the response characteristics. Preliminary results show that network interactions between PN and LNIs can also affect the PN response. Although SK channels can be responsible for the generations of I phase, parameters of $I_{Ca}$, $I_{Na}$ and $I_{A}$, as well as the synaptic currents can also affect the I phase. Therefore, experimentally blockers that affect any of these parameters might block the I phase.

## Methods

Based on various experimental findings (see S2 text, Experimental findings in ORNs and PNs), we developed models of ORN and PN and of the MGC neural network. The model parameters were fitted to the experimental data.

### The Poisson model of ORNs

To construct the Poisson model of ORNs, we fitted the extracellular recorded data [12] about the rise and fall of the mean instantaneous response frequency as a function of time following different concentrations and durations of the pheromone stimulus. At any concentration the rising phase of the frequency curve can be fitted by a single exponential function. However, the dynamics of the falling phase depend on the stimulation parameters. For short stimulation periods of 100 and 200 ms at any concentration, the falling phase can be fitted by the sum of two exponential functions, one fast with a small time constant $\tau_{fall1}$ and one slow with a larger time constant $\tau_{fall2}$ as shown in Fig 7A. The fitting function used in this case is given by Eq 1. The fitted curves are shown in Fig 7B and the blue (stimulation period: 100 ms, stimulation dose: 10 ng) and purple curves (stimulation period: 200 ms, stimulation dose: 10 ng) in Fig 7D. The fitted falling time constants decrease with pheromone concentration (Table 1). For stimulation concentration at 10 ng with long stimulation periods of 500 ms and 1000 ms, the falling phase of frequency undergoes two stages: a rapid falling stage to a plateau and a slow falling stage. The rapid falling stage can be fitted by one exponential function with a small time constant $\tau_{fall1}$ and the slow falling stage can be fitted by two exponential functions with one intermediate time constant $\tau_{fall2}$ and one larger time constant $\tau_{fall3}$ as shown in Fig 7C. The fitting function used in this case is given by Eq 2. The fitted curves are shown by the green and red curves in Fig 7D. The fitted parameter values of Eqs 1 and 2 are given in Table 1.

$$ \bar{f}(t) = \begin{cases} f_{sp}, & \text{if } t \le t_{sti} + T_{lat} \\ f_{sp} + (f_{pe} - f_{sp}) \cdot \left( 1 - e^{-\frac{t - t_{sti} - T_{lat}}{\tau_{rise}}} \right), & \text{else if } (t_{sti} + T_{lat}) \le t \le t_{sti} + T_{lat} + T_{d2pe} \\ f_{sp} + (f_{pe} - f_{sp}) \cdot \left( 1 - e^{-\frac{T_{d2pe}}{\tau_{rise}}} \right) \cdot \left( qe^{-\frac{t - (t_{sti} + T_{lat} + T_{d2pe})}{\tau_{fall1}}} + (1 - q)e^{-\frac{t - (t_{sti} + T_{lat} + T_{d2pe})}{\tau_{fall2}}} \right), & \text{otherwise} \end{cases} \quad (1) $$
<figure><img src="image_10.png" /><figcaption>**Fig 7.** Mean frequency response curves of the ORN population in response to different concentrations and durations of the pheromone stimulation. A. Response data curve (blue) and fitted curve (red) to stimulus: 10 ng and 200 ms. B. Fitted response curves to the same stimulation period 200 ms and different stimulation doses from 0.1 to 10 ng. C. Response data curve (blue) and fitted curve (red) to stimulus: 10 ng and 1000 ms. D. Fitted response curves to the same stimulation dose 10 ng and different stimulation periods from 100 to 1000 ms.</figcaption></figure>

doi:10.1371/journal.pone.0126305.g007

$$ \bar{f}(t) = \begin{cases} f_{sp}, & \text{if } t \le t_{sti} + T_{lat} \\ f_{sp} + (f_{pe} - f_{sp}) \cdot \left(1 - e^{-\frac{t - t_{sti} - T_{lat}}{\tau_{rise}}}\right), & \text{else if } (t_{sti} + T_{lat}) \le t \le t_{pe} \\ f_{pl} + \left(f_{sp} + (f_{pe} - f_{sp}) \cdot \left(1 - e^{-\frac{T_{d2pe}}{\tau_{rise}}}\right) - f_{pl}\right) e^{-\frac{t - (t_{sti} + T_{lat} + T_{d2pe})}{\tau_{fall1}}}, & (2) \\ f_{sp} + (f_{pl} - f_{sp}) \cdot \left( qe^{-\frac{t - (t_{sti} + T_{lat} + T_{d2pe} + T_{pl})}{\tau_{fall2}}} + (1 - q)e^{-\frac{t - (t_{sti} + T_{lat} + T_{d2pe} + T_{pl})}{\tau_{fall3}}} \right), & \text{otherwise} \end{cases} $$
where, $f_{sp}$, $f_{pe}$ and $f_{pl}$ are the mean spontaneous frequency of the ORN, peak frequency and plateau frequency in response to stimulation; $t_{sti}$ the time of stimulation onset; $T_{lat}$, $T_{d2pe}$ and $T_{pl}$ the response latency of the ORN population, the duration to peak frequency from $t_{sti}$ + $T_{lat}$ and duration of the plateau; $\tau_{rise}$, $\tau_{fall1}$, $\tau_{fall2}$ and $\tau_{fall3}$ the rising and falling time constants respectively; q coefficient of the fast falling component.

We modelled the ORN spike train by a Poisson process characterized by a single parameter, the mean firing rate $\bar{f}(t)$ given by Eqs (1) and (2). For sufficiently short interval $\delta t$, and a mean frequency $\bar{f}(t)$ varying slowly with respect to $\delta t$, the probability of a spike occurring during $\delta t$ is equal to the value of the instantaneous firing frequency during this time interval times the length of the interval

$$P\{1 \text{ spike during } \delta t\} \approx \bar{f}(t) \cdot \delta t. \quad (3)$$

At iteration time $t$, a random number $R[t]$, uniformly distributed between 0 and 1 is generated. If $R[t] \le \bar{f}(t) \cdot \delta t$ where $\delta t$ is the time step used in simulations, the membrane potential of ORNs is set to $V_{ORNS} = 50$ mV. Otherwise, $V_{ORNS} = -62$ mV (resting potential).

## The biophysical model of PN

The model is mathematically described by Hodgkin—Huxley type equations (Eqs (4–14)). The membrane activity of PN satisfies the following differential equation:

$$C_m \frac{dV}{dt} = -I_{Na} - I_{Ca} - I_{Kd} - g_L(V - E_L) - I_A - I_{SK} - I_{nAch}, \quad (4)$$

where $V$ is the membrane potential, $C_m$ the membrane capacitance, $g_L$ and $E_L$ the conductance and reversal potential of the leak current, respectively. The values of these passive parameters are given in Table 2.

The intrinsic currents of PN. The intrinsic inward ($I_{Na}$ and $I_{Ca}$) and outward ($I_A$, $I_{Kd}$ and $I_{SK}$) ionic currents in PN are described by

$$I_j = \bar{g}_j m^M h^N (V - E_j), \quad (5)$$

where $\bar{g}_j$ and $E_j$ are the maximal mean conductance and reversal potential for the ionic current $j$. The values of these two parameters of each current are given in Table 3. $M = 3$, $N = 1$ for $I_{Na}$ and $I_A$; $M = 1$, $N = 1$ for $I_{Ca}$; $M = 3$, $N = 0$ for $I_{Kd}$; $M = 2$, $N = 0$ for $I_{SK}$. The gating variables $m$ and $h$ in Eq (5) satisfy Eqs (6) and (7) except that $h = h_\infty$ for $I_{Ca}$.

$$\dot{m} = (m_\infty - m)/\tau_m, \quad (6)$$

$$\dot{h} = (h_\infty - h)/\tau_h, \quad (7)$$

where the steady-state activation $m_\infty$ and inactivation $h_\infty$ of the voltage-activated currents are described by Boltzmann equations as Eqs (8) and (9)

$$m_\infty = 1/\{1 + \exp[(V_{0.5act} - V)/S_m]\}, \quad (8)$$

$$h_\infty = 1/\{1 + \exp[(V - V_{0.5inact})/S_h]\} \quad (9)$$

The voltage dependency of the time constants of $m$ and $h$ of the voltage-activated currents is described by functions as Eqs (10) and (11) except that the $\tau_m$ of Ca$^{2+}$ current takes the form of
Eq (10')

$$ \tau_m(V) = \frac{1}{a_{\tau m, up} e^{(V_{\tau m, 0.5 up} - V)/S_{\tau m, up}} + a_{\tau m, dn} e^{(V - V_{\tau m, 0.5 dn})/S_{\tau m, dn}}}, \quad (10) $$

$$ \tau_{mCa}(V) = (a_{(\tau m, up)} \exp(-\frac{V}{S_{\tau m, up}}) + a_{(\tau m, dn)} \left( \frac{-V + V_{\tau m, dn}}{\exp\left(\frac{-V+V_{\tau m, dn}}{S_{\tau m, dn}} - 1\right)} \right) +)^{-1}, \quad (10) $$

$$ \tau_h(V) = \frac{1}{a_{th,up} e^{(V\tau h,0.5up-V)/Sth,up} + a_{th,dn} e^{(V-V\tau h,0.5dn)/Sth,dn}}, \quad (11) $$

For Na⁺ currents $I_{Na}$, value of parameter $V_{0.5act}$ in Eq (8) was taken from [38] (DUM cells of the cockroach *Periplaneta americana*), values of parameter $S_m$, $V_{0.5inact}$ and $S_h$ were modified from [38] and values of parameters in Eqs (10) and (11) were fitted to the data given in [38]. For Ca²⁺ currents $I_{Ca}$, values of parameters in Eqs (8) and (9) were taken from [22] (PN in *P. americana*), time constant for activation takes the form of Eq (10') described in [39], $h = h_∞$. For the sustained and transient voltage-gated K⁺ currents $I_{kd}$ and $I_A$, values of parameters in Eqs (8) and (9) were taken from [36] (MGC PN in male sphinx moth *Manduca sexta*), and values of parameters in Eqs (10) and (11) were fitted to the data given in [37] (in *M. sexta*). The values of various parameters in the voltage dependent steady-state and time constant function of $I_{Na}$, $I_{Ca}$, $I_{Kd}$ and $I_A$ are given in Table 3.

The mathematical description of $m_∞$ and current of the Ca²⁺-dependent K⁺ currents $I_{SK}$ and that of the Ca²⁺ dynamics were borrowed directly from [39] as Eqs (12–14)

$$ \frac{dCa}{dt} = -f_{Ca}I_{Ca} - (Ca - Ca_{\infty})/\tau_{Ca}, \quad (12) $$

$$ m_{SK\infty} = 1/(1 + e^{-a_{msk} - b_{msk} \log \frac{C_a - C_{a\infty}}{S_{msk}}}), \quad (13) $$

$$ I_{SK} = \bar{g}_{SK} \cdot m_{SK\infty}^2 (V - E_K), \quad (14) $$

where the values of parameters of Ca²⁺ dynamics are given in Table 2 and those of $I_{SK}$ are given in Table 3.

The cholinergic synaptic current from ORNs to PN. The fast nicotinic cholinergic synaptic currents calculated according to

$$ I_{nACh} = \sum_{i=1}^{N} \bar{g}_{nACh} \cdot [O]_{i}(t) \cdot (V - E_{nACh}), \quad (15) $$

where N is the number of ORNs, $\bar{g}_{nACh}$ the mean peak conductance, and $E_{nACh} = 0$ mV the reversal potential of the current respectively. The fraction of open channels $[O]_i$ is modeled by first-order activation scheme (see review in [40])

$$ \frac{d[O]_i}{dt} = \alpha(1 - [O]_i)[T]_i - \beta[O]_i \quad (16) $$

The release of cholinergic transmitter $[T]_i$ from ith ORN was modeled by a square pulse

$$ [T]_i = A\theta(t_0 + t_{max} - t)\theta(t - t_0) \quad (17) $$

Parameter values of the nACh synaptic current are given in Table 4.
## The MGC network model

We constructed a MGC network model by connecting 100 ORNs to one PN as shown in S1 Fig. No LNs were included in the MGC model in order to test the hypothesis that the triphasic firing patterns of PN can be generated by the ionic currents in PN and ORN inputs. Computer simulations of the model were performed in Microsoft visual studio 2008. The simulation results were analyzed with Matlab 7.5. The total computer simulation time is 25s and the pheromone stimulation started at 5s.

## Low-pass and high-pass Butterworth filters

In order to better understand how different ionic currents contribute to the generation of the PN firing pattern we separated the slow and fast components of the depolarizing and repolarizing currents in PN. We designed 10th-order lowpass and highpass Butterworth filters with cut-off frequency 5 Hz using the Matlab function "butter". By applying the designed lowpass and highpass filters the slow and fast components of each ionic currents were extracted.

## Analysis of the PN response characteristics

The PN response pattern is quantitatively characterized by duration of E₁ and I phases and frequency of E₁ and E₂ phases. These features were defined in S9 Fig and were calculated as follows and expressed as means ± standard error of the mean.

**Duration.** The E₁ durations were measured from the first spike of the PN response to the spike just preceding the inhibitory phase. The I durations were measured from the last spike of E₁ to the first spike of E₂.

**Frequency.** We first calculated the interspike intervals (ISIs) between successive spikes. Then the ISIs were averaged in 10 spikes around each spike using Matlab function smooth. The frequencies are the inverse of the ISIs.

## Supporting Information

**S1 Fig. Simplified model of moth MGC.** The model is composed of 100 Poisson ORNs and one biophysical PN. ORNs receive pheromone stimuli and the PN receive Ach synaptic inputs from ORNs through nicotinic receptors at the dendrites. (TIF)

**S2 Fig. Effects of parameters for $m_\infty$, $\tau_m$ function of $I_{Ca}$ and for $m_{sk\infty}$ function of $I_{SK}$ on PN response duration.** (TIF)

**S3 Fig. Effects of parameters for $\tau_m$, $\tau_h$ and $h_\infty$ function of $I_{Na}$ on PN response characteristics.** I duration clearly increases with $S_{rm,up}$ (A) and $V_{\tau h>0.5dn}$ (G), while it decreases with $a_{\tau h dn}$ (C). $E_1$ frequency linearly decreases with $V_{0.5inact}$ (F); $E_2$ frequency increases with $a_{rh,dn}$ (D), $S_{rh,dn}$ (B) and $V_{\tau h>0.5up}$ (E) while it decreases with $V_{0.5inact}$ (F) and $V_{\tau h>0.5dn}$ (H). (TIF)

**S4 Fig. Effects of parameters of $I_{Kd}$ on PN response characteristics.** $I_{Kd}$ clearly affects $E_1$ frequency which increases with $\bar{g}_{Kd}$ when $\bar{g}_{Kd}$ is below 0.7 µS then decreases (B) while it decreases with $V_{0.5act}$ (E) and $S_m$ of $I_{Kd}$ (H). (TIF)

**S5 Fig. Effects of Ach pulse duration $t_{max}$ and concentration A on PN response characteristics.** Top panel: effects of $t_{max}$ on $E_1$ and I duration (A) and mean firing frequency of $E_1$ and $E_2$
(B). Bottom panel: effects of A on E₁ and I duration (C) and mean firing frequency of E₁ and E₂ phases (D).

(TIF)

**S6 Fig.** Effects of nACh postsynaptic channel opening rate α, closing rate β and $\bar{g}_{nACh}$ on PN response characteristics. Top panel: effects of α on E₁ and I duration (A) and mean firing frequency of E₁ and E₂ (B). Middle panel: effects of β on E₁ and I duration (C) and mean firing frequency of E₁ and E₂ phases (D). Bottom panel: effects of $\bar{g}_{nACh}$ on E₁ and I duration (E) and mean firing frequency of E₁ and E₂ (F).

(TIF)

**S7 Fig.** Effects of inhibition mediated by fast GABA<sub>A</sub> receptors on the activity patterns of PN and three LNIs in a network with 80 ORNs, 1 PN and 20 LNIs. The postsynaptic channel closing rate β of the GABA synapses from LNI to PN is 0.1 (left panel), 2.0 (middle panel) and 3.0 (right panel) respectively.

(TIF)

**S8 Fig.** Effects of slow inhibition mediated by metabotropic GABA<sub>B</sub> receptors on the response activity of PN and three LNIs in a network with 80 ORNs, 1 PN and 20 LNIs. Panel A-D show that the GABA<sub>B</sub> mediated synaptic inhibition did not alter the triphasic response pattern of PN in the normal parameter range: A. PN potential, spikes and frequency; B. potentials of LNI8; C. normalized concentration of intracellular Ca in PN and of GABA<sub>B</sub> receptor-coupled G protein; D. potentials of LNI17. Panel E-F show that the I duration is prolonged when r<sub>3</sub> is decreased: E. PN potential, spikes and frequency; F. normalized concentration of intracellular Ca in PN and of GABA<sub>B</sub> receptor-coupled G protein. Panel G-H show that the GABA<sub>B</sub> mediated synaptic inhibition changed the triphasic response pattern when r<sub>3</sub> is increased and r<sub>4</sub> is decreased: G. PN potential, spikes and frequency; H. normalized concentration of intracellular Ca in PN and of GABA<sub>B</sub> receptor-coupled G protein.

(TIF)

**S9 Fig.** Extracellularly recorded response patterns of the moth pheromone sensitive PN in MGC. Left panel: top trace, response pattern to low dose pheromone stimulus; bottom trace, response pattern to high dose pheromone stimulus. Right panel: from top to bottom trace the duration of pheromone stimuli were increased at a given stimulation concentration.

(TIF)

**S1 Text.** Model of type I LNs with sodium spikes.

(PDF)

**S2 Text.** Experimental findings in ORNs, PNs and type I LNs.

(PDF)

## Author Contributions

Conceived and designed the experiments: YG. Performed the experiments: YG. Analyzed the data: YG. Contributed reagents/materials/analysis tools: YG. Wrote the paper: YG.

## References

1. Youngentob SL, Mozell MM, Sheehe PR, Hornung DE. A quantitative analysis of sniffing strategies in rats performing odor detection tasks. Physiol Behav. 1987; 41: 59–69. PMID: 3685154

2. Willis MA, Baker TC. Effects of intermittent and continuous pheromone stimulation on the flight behaviour of the oriental fruit moth, *Grapholita molesta*. Physiol Entomol. 1984; 9: 341–358.
3. Vickers NJ, Baker TC. Male *Heliothis virescens* sustain upwindflight in response to experimentally pulsed filaments of their sex pheromone. J Insect Behav. 1992; 5: 669-687.

4. Vickers NJ, Baker TC. Reiterative responses to single strands of odor promote sustained upwind flight and odor source location by moths. Proc Natl Acad Sci USA. 1994; 91: 5756-5760. PMID: 11607476

5. Mafra-Neto A, Cardé RT. Fine-scale structure of pheromone plumes modulates upwind orientation of flying moths. Nature. 1994; 369: 142-144.

6. Moore PA. A model of the role of adaptation and disadaptation in olfactory receptor neurons. Chem Senses. 1994; 19: 71-86. PMID: 8055260

7. Atema J. Chemical signals in the marine environment: dispersal, detection, and temporal signal analysis. Proc Natl Acad Sci USA. 1995; 92: 62-66. PMID: 7816848

8. Cardé RT, Minks AK. Insect pheromone research: new directions. New York: Chapman and Hall; 1997.

9. Hansson BS, Anton S. Function and morphology of the antennal lobe: new developments. Annu Rev Entomol. 2000; 45: 203-231. PMID: 10761576

10. Christensen TA, Waldrop BR, Hildebrand JG. Multitasking in the olfactory system: context-dependent responses to odors reveal dual GABA-regulated coding mechanisms in single olfactory projection neurons. J Neurosci. 1998; 18(15): 5999-6008. PMID: 9671685

11. Christensen TA, Heinbockel T, Hildebrand JG. Olfactory information processing in the brain: encoding chemical and temporal features of odors. J Neurobiol. 1996; 30(1): 82-91. PMID: 8727985

12. Jarriault D, Gadenne C, Lucas P, Rospars J-P, Anton S. Transformation of the sex pheromone signal in the noctuid moth *Agrotis ipsilon*: from peripheral input to antennal lobe output. Chem Senses 2010; 35: 705-715. doi: 10.1093/chemse/bjq069 PMID: 20601375

13. Rospars J-P, Grémiaux A, Jarriault D, Chaffiol A, Monsempes C, Deisig N, et al. Heterogeneity and convergence of olfactory first-order neurons account for the high speed and sensitivity of second-order neurons. PLoS Comput Biol. 2014; 10(12): e1003975. doi: 10.1371/journal.pcbi.1003975 PMID: 25474026

14. Christensen TA, Waldrop BR, Harrow ID, Hildebrand JG. Local interneurons and information processing in the olfactory glomeruli of the moth *Manduca sexta*. J Comp Physiol A. 1993; 173(4): 385-399. PMID: 8254565

15. Wellis DP, Scott JW, Harrison TA. Discrimination among odorants by single neurons of the rat olfactory bulb. J Neurophysiol. 1989; 61(6): 1161-1177. PMID: 2746317

16. Hamilton KA, Kauer JS. Patterns of intracellular potentials in salamander mitral/tufted cells in response to odor stimulation. J Neurophysiol. 1989; 62(3): 609-625. PMID: 2549211

17. Mori K, Yoshihara Y. Molecular recognition and olfactory processing in the mammalian olfactory system. Prog Neurobiol. 1995; 45(6): 585-619. PMID: 7624486

18. Laurent G. Dynamical representation of odors by oscillating and evolving neural assemblies. Trends Neurosci. 1996; 19(11): 489-496. PMID: 8931275

19. Brown SL, Joseph J, Stopfer M. Encoding a temporally structured stimulus with a temporally structured neural representation. Net Neurosci. 2005; 8(11): 1568-1576. PMID: 16222230

20. Kanzaki R, Arbas EA, Strausfeld NJ, Hildebrand JG. Physiology and morphology of projection neurons in the antennal lobe of the male moth *Manduca sexta*. J Comp Physiol A. 1989; 165: 427-453. PMID: 2769606

21. Martinez D, Chaffiol A, Voges N, Gu Y, Anton S, Rospars J-P, et al. Multiphasic On/Off pheromone signaling in moths as neural correlates of a search strategy. PLoS ONE. 2013; 8(4): e61220. doi: 10.1371/journal.pone.0061220 PMID: 23613816

22. Husch A, Paehler M, Fusca D, Paeger L, Kloppenburg P. Calcium current diversity in physiologically different local interneuron types of the antennal lobe. J Neurosci. 2009; 29(3): 716-726. doi: 10.1523/JNEUROSCI.3677-08.2009 PMID: 19158298

23. Mercer AR, Hildebrand JG. Developmental changes in the density of ionic currents in antennal-lobe neurons of the sphinx moth, *Manduca sexta*. J Neurophysiol. 2002; 87: 2664-2675. PMID: 12037169

24. Abou Tayoun AN, Li X, Chu B, Hardie RC, Juusola M, Dolph PJ. The Drosophila SK channel (dSK) contributes to photoreceptor performance by mediating sensitivity control at the first visual network. J Neurosci. 2011; 31: 13897-13910. doi: 10.1523/JNEUROSCI.3134-11.2011 PMID: 21957252

25. Defaix C, Anton S, Rospars J-P, Martinez D, Lucas P. Firing properties and ionic currents of antennal lobe neurons of a moth. Poster in the 11th Neural Coding International Workshop, Versailles, France, 2014.

26. Av-Ron E. The role of a transient potassium current in a bursting neuron model. J Math Biol 1994; 33: 71-87. PMID: 7836871
27. Av-Ron E, Rospars J-P. Modeling insect olfactory neuron signaling by a network utilizing disinhibition. BioSystems. 1995; 36: 101–108. PMID: 8573691

28. Bazhenov M, Stopfer M, Rabinovich M, Huerta R, Abarbanel HDI, Sejnowski TJ, et al. Model of transient oscillatory synchronization in the locust antennal lobe. Neuron. 2001a; 30: 553–567. PMID: 11395014

29. Bazhenov M, Stopfer M, Rabinovich M, Huerta R, Abarbanel HDI, Laurent G. Model of cellular and network mechanisms for odour-evoked temporal patterning in the locust antennal lobe. Neuron. 2001b; 30: 569–581.

30. Sivan E, Kopell N. Oscillations and slow patterning in the antennal lobe. J Comput Neurosci. 2006; 20: 85–96. doi: 10.1007/s10827-006-4087-z PMID: 16511657

31. Gu Y, Liljenström H. Modelling efficiency in insect olfactory information processing. BioSystems. 2007; 89: 236–243. PMID: 17307286

32. Belmabrouk H, Nowotny T, Rospars J-P, Martinez D. Interaction of cellular and network mechanisms for efficient pheromone coding in moths. Proc Natl Acad Sci USA. 2011; 108(49): 19790–19795. doi: 10.1073/pnas.1112367108 PMID: 22109556

33. Assisi C, Stopfer M, Bazhenov M. Using the structure of inhibitory networks to unravel mechanisms of spatiotemporal patterning. Neuron. 2011; 69: 373–386. doi: 10.1016/j.neuron.2010.12.019 PMID: 21262473

34. Buckley CL, Nowotny T. Multiscale model of an inhibitory network shows optimal properties near bifurcation. Phys Rev Lett. 2011; 106(23): 238109(4). PMID: 21770552

35. Anton S, Homberg U. Antennal lobe structure. In: Hansson BS, editor. Insect olfaction. Springer, Berlin, Heidelberg; 1999. pp. 97–124.

36. Kloppenburg P, Ferns D, Mercer AR. Serotonin enhances central olfactory neuron responses to female sex pheromone in the male sphinx moth *Manduca sexta*. J Neurosci. 1999; 19(19): 8172–8181. PMID: 10493719

37. Mercer AR, Hayashi JH, Hildebrand JG. Modulatory effects of 5-hydroxytryptamine on voltage-activated currents in cultured antennal lobe neurons of the sphinx moth *Manduca sexta*. J Exp Biol. 1995; 198: 613–627. PMID: 7714451

38. Lapied B, Malecot CO, Pelhate M. Patch-clamp study of the properties of the sodium current in cockroach single isolated adult aminergic neurons. J Exp Biol. 1990; 151: 387–404.

39. Roper P, Callaway J, Shevchenko T, Teruyama R, Armstrong W. AHP's, HAP's and DAP's: How potassium currents regulate the excitability of rat supraoptic neurones. J Comput Neurosci. 2003; 15: 367–389. PMID: 14618071

40. Destexhe A, Mainen ZF, Sejnowski TJ. Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism. J Comp Neurosci. 1994; 1: 195–230.

41. Bettler B, Kaupmann K, Mosbacher J, Gassmann M. Molecular Structure and Physiological Functions of GABA<sub>B</sub> Receptors. Physiol Rev. 2004; 84:835–867, doi: 10.1152/physrev.00036.2003 PMID: 15269338