Pacemaker neuron and network oscillations depend on a neuromodulator-regulated linear current

Linear leak currents have been implicated in the regulation of neuronal excitability, generation of neuronal and network oscillations, and network state transitions. Yet, few studies have directly tested the dependence of network oscillations on leak currents or explored the role of leak currents on network activity. In the oscillatory pyloric network of decapod crustaceans neuromodulatory inputs are necessary for pacemaker activity. A large subset of neuromodulators is known to activate a single voltage-gated inward current I(MI), which has been shown to regulate the rhythmic activity of the network and its pacemaker neurons. Using the dynamic clamp technique, we show that the crucial component of I(MI) for the generation of oscillatory activity is only a close-to-linear portion of the current-voltage relationship. The nature of this conductance is such that the presence or the absence of neuromodulators effectively regulates the amount of leak current and the input resistance in the pacemaker neurons. When deprived of neuromodulatory inputs, pyloric oscillations are disrupted; yet, a linear reduction of the total conductance in a single neuron within the pacemaker group recovers not only the pacemaker activity in that neuron, but also leads to a recovery of oscillations in the entire pyloric network. The recovered activity produces proper frequency and phasing that is similar to that induced by neuromodulators. These results show that the passive properties of pacemaker neurons can significantly affect their capacity to generate and regulate the oscillatory activity of an entire network, and that this feature is exploited by neuromodulatory inputs.

section_content: hormonally or by descending projection neurons (Nusbaum and Beenhakker, 2002). When neuromodulators are removed, e. g., by decentralization, the bursting oscillations are usually disrupted and the pyloric neurons either spike tonically or become quiescent (Marder and Bucher, 2006). Using the dynamic clamp technique, we show that in a decentralized preparation, a reduction in the level of linear (leak) conductance in a single pyloric pacemaker neuron recovers oscillatory activity in that neuron, and interestingly, brings back the entire network oscillation. 

Many of peptidergic and other modulators converge to unmask a voltage-dependent inward current (I MI ) characterized by an inverted bell-shaped current-voltage relationship (Golowasch and Marder, 1992b; Swensen and Marder, 2000). An inverted bellshaped I-V curve has traditionally been thought to confer a region of negative-slope conductance to the overall I-V curve of a neuron, enhancing their excitability and their tendency to produce regenerative membrane potential changes, including slow oscillations. It could be argued that an inverted bell-shaped I-V curve could be approximated by two linear segments, one with negative and one with positive slope. The positive-conductance component can depolarize the cell but reduces excitability whereas the negativeconductance component both increases excitability and reduces the stability of the membrane potential near the resting potential. Our dynamic clamp results indicate that regulation of a linear conductance with a reversal potential near the pacemaker neuron's

section_content: Recent studies have shown that linear leak currents play an important role in the generation of network oscillations (Cymbalyuk et al., 2002; Blethyn et al., 2006; Koizumi and Smith, 2008; Pang et al., 2009) and regulating excitability (Rekling et al., 2000; Brickley et al., 2007). Neuromodulatory substances are known to modify neuronal excitability by activating metabotropic receptors that can regulate ionic currents and pumps (Bayliss et al., 1992; Cantrell and Catterall, 2001; LeBeau et al., 2005; Tobin and Calabrese, 2005). In particular, the negative regulation of leak currents by several neurotransmitters has been proposed to play an important role in the regulation of neuronal excitability (Bayliss et al., 1992; Lee and Boden, 1997; Aller et al., 2005; Kim, 2005; Pratt and Aizenman, 2007; Weber et al., 2008), oscillations (Xu et al., 2007) and sleepwake transitions (Goldstein et al., 2001). Nevertheless, few studies have directly tested the dependence of pacemaker and network oscillatory activity on leak currents and the role of neuromodulators in this process. 

In this study we examine the role of linear conductances on the generation of oscillations by the pyloric network of the crab Cancer borealis. The pyloric network, located in the stomatogastric ganglion (STG), generates a rhythmic pattern of bursting activity driven by a pacemaker group of one anterior burster (AB) and two pyloric dilator (PD) neurons. The pyloric oscillations are dependent on and regulated by a variety of neuromodulators that are released Pacemaker neuron and network oscillations depend on a neuromodulator-regulated linear current Shunbing Zhao 1, Jorge Golowasch 1,2 and Farzan Nadim 1,2 * 1 Federated Department of Biological Sciences, Rutgers University and New Jersey Institute of Technology, Newark, NJ, USA 2 Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, NJ, USA Linear leak currents have been implicated in the regulation of neuronal excitability, generation of neuronal and network oscillations, and network state transitions. Yet, few studies have directly tested the dependence of network oscillations on leak currents or explored the role of leak currents on network activity. In the oscillatory pyloric network of decapod crustaceans neuromodulatory inputs are necessary for pacemaker activity. A large subset of neuromodulators is known to activate a single voltage-gated inward current I MI, which has been shown to regulate the rhythmic activity of the network and its pacemaker neurons. Using the dynamic clamp technique, we show that the crucial component of I MI for the generation of oscillatory activity is only a close-to-linear portion of the current-voltage relationship. The nature of this conductance is such that the presence or the absence of neuromodulators effectively regulates the amount of leak current and the input resistance in the pacemaker neurons. When deprived of neuromodulatory inputs, pyloric oscillations are disrupted; yet, a linear reduction of the total conductance in a single neuron within the pacemaker group recovers not only the pacemaker activity in that neuron, but also leads to a recovery of oscillations in the entire pyloric network. The recovered activity produces proper frequency and phasing that is similar to that induced by neuromodulators. These results show that the passive properties of pacemaker neurons can significantly affect their capacity to generate and regulate the oscillatory activity of an entire network, and that this feature is exploited by neuromodulatory inputs. 

where E MI = -10 mV and τ m = 4 ms. The values for g MI _ varied depending on the experiment as described in the Results. The standard leak current is described by

) where E L = -68 mV is the leak current reversal potential and g L is the leak conductance. When the value of leak conductance was set to be negative in the dynamic clamp software it resulted in a reduction of the total neuronal conductance. We describe this current as

The voltage-dependent current I MI was divided into regions approximating a linear I-V curve with positive linear conductance (MI-PL) or negative linear conductance (MI-NL). The linear conductance regions were restricted in their voltage range so that they represented the negative and positive-conductance regions of I MI :

where τ = 4 ms. Although I MI is non-inactivating, to build I MI-NL, an inactivation variable h N with very fast kinetics is used in dynamic clamp as a method of restricting the voltage range of activation below the upper limit of -46 mV (see Figure 6A, red trace). A similar use of the activation variable m P was made to produce I MI-PL (Figure 6A, gray trace). 

It should be noted that dynamic clamp conductance manipulations introduce artificial ionic currents into the point of penetration, in this case the soma. The ionic currents in biological neurons are distributed over the entire structure of the cell and it is possible therefore that the size of the currents injected by dynamic clamp does not properly represent the equivalent biological ionic current. In particular, the currents injected in our dynamic clamp experiments may be much larger than what is actually needed for the pacemaker neurons to produce oscillations. 

We measured neuronal input conductance before and after injecting negative conductances with dynamic clamp to avoid reducing the total conductance of the neuron to less than zero. When simultaneous dynamic clamp current and constant current injections were required, a Brownlee Precision Amplifier (Brownlee Precision) was used to add the currents before input to the Axoclamp amplifier. The acquired data were saved as binary files and were analyzed with the software Readscope (http://stg. rutgers. edu/software) and Clampfit 9. 2 (Molecular Devices). 

For experiments performed to compare the effect of neuromodulatory input, the STG was isolated from anterior neuromodulatory inputs (i. e. decentralized) by replacing the saline solution within a resting potential, controls the ability of the pyloric pacemaker neuron to generate oscillations. Furthermore, we show that the only essential element of the non-linear I MI that contributes to the enhancement of membrane excitability and promoting membrane potential oscillations is the negative-conductance region spanning a voltage range that is aligned with the resting potential of the pacemaker neuron.
section_name: IntroductIon

section_content: Experiments were carried on adult male crabs (Cancer borealis) purchased from local distributors (Newark, NJ, USA). The animals were maintained in artificial seawater tanks at 12-15ºC until they were used, when they were first anesthetized by cooling on ice for 15-30 min prior to each dissection. The stomatogastric nervous system (including the stomatogastric ganglion, STG, the esophageal ganglion, OG, and the commissural ganglia, CoGs, Figure 1A ) was removed using standard methods (Selverston et al., 1976; Harris-Warrick et al., 1992) and each preparation was pinned down in a Sylgard-coated Petri dish. The STG was desheathed to expose the cell bodies for microelectrode impalement, and to allow effective superfusion. Preparations were superfused using normal saline containing (in mM): 11 KCl, 440 NaCl, 13 CaCl 2, 26 MgCl 2, 11. 2 Trizma base, 5. 1 Maleic acid, pH 7. 4-7. 5, and kept at ∼12ºC. 

Microelectrodes were pulled using a Flaming-Brown micropipette puller (P97, Sutter Instruments) and filled with 0. 6 M K 2 SO 4 + 0. 02 M KCl (resistance of 15-25 MΩ). Extracellular recordings were made with two stainless steel pins, one placed in the bath and the other within a petroleum jelly well built around the nerve of interest. Neuronal identification was accomplished by matching the intracellular action potential recordings to their corresponding extracellular nerve recordings (Selverston et al., 1976; Harris-Warrick et al., 1992). Intracellular recordings were made from the soma of the neurons using Axoclamp 2B amplifiers (Molecular Devices) and extracellular recordings were amplified using a differential AC amplifier model 1700 (A-M Systems). 

We used the dynamic clamp technique (Sharp et al., 1992) to introduce ionic conductances into, or subtract them from, the biological neurons. A NI PCI-6070-E board (National Instruments, Austin, TX, USA) was used for current injection in dynamic clamp experiments. Data acquisition was performed using the Digidata 1332A data acquisition board and pClamp 9. 2 software (Molecular Devices). The dynamic clamp software was developed in our laboratories (available for download at http://stg. rutgers. edu/software) in the LabWindows/CVI software environment (National Instruments) on a Windows platform. In dynamic clamp, ionic currents are calculated using Hodgkin-Huxley-type equations as described below and continuously updated by recording the membrane potential (V) of the neurons in real time. I MI was described by:

petroleum jelly well built around the single stomatogastric nerve (stn, Figure 1A ) with 0. 75-M sucrose and 1-μM TTX (Biotium) to block action potential transmission. In experiments to measure the input resistance of neurons 1-μM TTX was used to superfuse the ganglion. SAS (SAS Institute), SigmaStat 2. 3 (Systat Software) and Origin 7 (OriginLab) software packages were used for statistical and graphical analysis and all final figures were made in CorelDraw 12 (Corel Corp. ). Reported statistical significance indicated a significance level p ≤ 0. 05. In all statistical tests the data have been checked for normality by the SAS software package. All error bars shown and error values reported are standard deviations of the mean (SD).
section_name: MaterIals and Methods

section_content: The pyloric network is a central pattern generator located in the stomatogastric ganglion, one of the four ganglia in the stomatogastric nervous system, an extension of the crab CNS (Figure 1A ). The pyloric network is an extremely well-studied network and consists, primarily, of six types of neurons that are connected with inhibitory synapses as well as gap junctions (Figure 1B ) (Nusbaum and Beenhakker, 2002; Marder and Bucher, 2006). Oscillations in the pyloric network are relatively stereotypical and consist of three phases with a cycle frequency of ∼1 Hz (Figure 1C ). All pyloric neurons take part in the network oscillation by producing bursts of action potentials on top of slow wave oscillations of their membrane potential. These oscillations are generated by a pacemaker group consisting of the AB and two PD neurons which are strongly coupled through gap junctions and produce synchronous bursting activity. 

Pyloric oscillations are conditional upon the presence of neuromodulators that are released hormonally and in paracrine fashion from descending neurons that project to the STG via the stomatogastric nerve (stn; Figure 1A ). Blocking action potential conduction in the stn (decentralization) in vitro results in the disruption of the pyloric oscillations (Figure 1C ). Neuromodulatory inputs have been shown to target ionic currents in all pyloric neurons and modify the strengths of chemical synapses and gap junctions in the network (Johnson et al., 1994; Dickinson, 2006; Marder and Bucher, 2006). Although distinct modulators diverge in their target neurons, many modulators of the pyloric network, particularly neuropeptides and muscarinic agonists, have been shown to converge to activate a single voltage-gated inward current (I MI ) (Golowasch Oscillations regulated by leak current and Marder, 1992b; Swensen and Marder, 2000). It should be noted that not all modulatory actions converge on this current and a variety of other ionic currents are also regulated by neuromodulators, particularly neuroamines (Harris-Warrick et al., 1995; Johnson et al., 2003). 

The intact pyloric network is presumed to be subject to modulation by a variety of substances, each targeting a subset of neurons and synapses, the combined actions of which result in the observed triphasic pyloric oscillations. Here we examine the hypothesis that the primary action of neuromodulators that results in the triphasic pyloric oscillation is a simple modification of the leakiness of the pacemaker neurons. 

During pyloric oscillations the pacemaker neurons AB and PD produce synchronous oscillations and follower neurons such as the LP neuron burst out of phase with the pacemakers (Figures 1C and 2A ). When the preparation was decentralized, bursting oscillations in all the neurons ceased and the slow wave oscillation of their membrane potentials were suppressed and replaced with a almost completely quiescent membrane potential of -60 to -50 mV (Figure 2B ). Using the dynamic clamp technique, we reduced the magnitude of the leak current (i. e. increased the input resistance R IN ) in individual pyloric neurons in a way that the resting potential was not significantly changed. This was done by adding a negative linear conductance (Eq. 2) with a reversal potential near the resting potential of the neuron. A small value of g NL in any of the neurons did not affect the output of a decentralized preparation. When a sufficiently large g NL value was added to the AB neuron (N = 3), bursting activity in the AB neuron resumed and other pyloric neurons also resumed their oscillations with a triphasic pyloric activity pattern (Figure 2C ). All oscillatory activity ceased immediately as soon as the dynamic clamp current was removed (not shown). Similarly, if a sufficiently large g NL was added to either one of the two PD neurons, it resulted in the recovery of the triphasic pyloric activity (Figure 2D ; N = 8). In contrast, addition of an I NL to the LP neuron never resulted in the recovery of oscillations in either the LP neuron or the network no matter what conductance value g NL was used (Figure 2E ; N = 3). 

These results indicated that reducing the amount of leak current in a single pacemaker neuron, but not a follower neuron, was sufficient to recover the oscillations of the entire pyloric network. There appeared to be no difference in the ability to recover the oscillations by the dynamic clamp application of I NL in either the AB or PD neurons. Thus, from this point on we focused our experiments on the PD neurons because the single small AB neuron is often difficult to find and harder to record from than the PD neurons. Henceforth any reference to I NL is in the PD neuron unless otherwise specified. 

We examined the effect of the conductance of the I NL on the recovered pyloric oscillation. Although in different preparations different amounts of g NL led to the recovery of pyloric oscillations, increasingly negative g NL in a recovered rhythm always led to a faster oscillation frequency (Figure 3A ). On average, a g NL value of -70 to -80 nS was required to bring back the oscillation frequency to the average control value observed before decentralization (horizontal line). We also examined the effect of I NL on the phase of activity of the individual pyloric neurons. The phase of individual neurons in the pyloric oscillation is defined as the fraction of each cycle duration (T) during which that neuron is active (Figure 3B ). As seen in Figure 3C, the I NL not only recovered network oscillations but produced activity phases (open boxes) that were statistically similar to those in the control pyloric oscillation recorded before decentralization (filled boxes; two-way RM ANOVA; All pair-wise p values > 0. 1; N = 6). This was observed in all preparation in which I NL led to a recovery of rhythmic pyloric activity. 

The fact that pyloric oscillations can be recovered in a decentralized preparation by reducing the leak conductance in the pacemaker neurons, leads to the hypothesis that modulatory inputs result in a conductance decrease in these neurons. We examined this hypothesis by measuring the leak conductance in the PD neurons in Oscillations regulated by leak current control and in a decentralized preparation in which action potential transmission along the stn was blocked with only 0. 75-M sucrose to allow reversal of the block. These experiments were done in twoelectrode voltage clamp in which the PD membrane potential was ramped from a holding value of -60 mV to a final value of -80 mV in 500 ms and the resulting I-V curve was fit with a linear function to measure the slope (Figure 4A ). The current that is blocked by decentralization can then be measured as a difference current (Figure 4B ). The average reversal potential for the difference current was -63. 6 ± 9. 5 mV. 

To prevent any possible contamination of the leak conductance measurement by activation of voltage-gated currents the neuron was kept in two-electrode voltage clamp while the preparation was decentralized (see Materials and Methods) and while it recovered from decentralization after sucrose was washed from around the stn. These experiments showed that the leak conductance in the PD neuron was significantly increased when the preparation was decentralized, compared to control and recovery (Figure 4C ; N = 6; p = 0. 04), supporting the hypothesis that the conductance in the decentralized system is reduced when neuromodulators are restored. 

As mentioned above, several neuromodulators of the pyloric network are known to activate a single voltage-gated inward current I MI in pyloric neurons (Swensen and Marder, 2000). I MI was first characterized through actions of proctolin (Golowasch and Marder, 1992b), an endogenously released neuromodulator that is known Oscillations regulated by leak current to activate the pyloric oscillations in a decentralized preparation (Nusbaum and Marder, 1989). If neuromodulators could affect the rhythm by affecting the pacemaker's leak current or by targeting I MI it is important to know if these effects are similar or related. We therefore compared the effect of I NL with adding a dynamic clamp version of I MI in the PD neuron and with bath application of 10 -6 M proctolin in decentralized preparations. We found that the dynamic clamp injections of I NL or I MI and bath application of proctolin all had very similar effects in recovering the pyloric oscillation. A comparison of these effects on the PD neuron is shown in Figure 5. This result led us to the idea that I NL and I MI act on the PD neuron through similar mechanisms. How could the actions of a non-linear inward current and a linear current with negative conductance be similar? We explored this question by examining the I-V curve of I MI. As shown in Figure 6A, this typical inverted bell-shaped I-V curve (black trace) can be divided into two regions, one with negative-slope conductance (red trace) and the other with positive slope conductance (gray trace). The section of the I-V curve with negative-slope conductance can be argued to be well approximated with I NL (dashed blue curve in Figure 6A ). We denoted the negative (positive) slope conductance section of I MI as I MI-NL (I MI-PL ) and compared the effect of dynamic clamp injection of I NL, I MI, I MI-NL and I MI-PL. These currents were modeled in dynamic clamp with the appropriate parameters so that I MI was well approximated by I MI-NL + I MI-PL (Eq. 3) and I NL matched the negative slope of I MI-NL in the appropriate voltage range as shown in Figure 6A (see Materials and Methods). When I MI-NL, I NL or or I MI were added to the PD neuron using dynamic clamp the result was always the recovery of oscillatory activity (Figures 6D, F, G ). In contrast, adding I MI-PL had little effect on the PD neuron and never resulted in oscillatory activity (Figure 6E ). These results indicated that the actions of I MI on the PD neuron were mainly through the limited region of negativeslope conductance which effectively resulted in a reduction of the leak current in the PD neuron.
section_name: results

section_content: Regulation of leak currents has been proposed to play an important role in the control of neuronal excitability, oscillatory activity and sleep-wake transitions (Bayliss et al., 1992; Lee and Boden, 1997; Goldstein et al., 2001; Aller et al., 2005; Kim, 2005; Pratt and Aizenman, 2007; Xu et al., 2007; Weber et al., 2008). A number of neuromodulators acting on metabotropic receptors have been shown to negatively regulate leak currents by decreasing their conductance; these include glutamate (Blethyn et al., 2006; Brickley et al., 2007), thyrotropin-releasing hormone (TRH) (Bayliss et al., 1992) and serotonin (Weber et al., 2008). We show that the action of the neuromodulatory inputs in producing pyloric oscillations can be mimicked by a simple reduction in the leak current of the pyloric pacemaker neurons. When a leak current with negative conductance is injected, using the dynamic clamp technique, in any one of the neurons that comprise the pyloric pacemaker group, the rhythmic pattern of the entire network is recovered. This recovered rhythm is essentially identical to the ongoing pyloric rhythm prior to decentralization and to the rhythm produced by bath application of a modulatory neuropeptide in a decentralized preparation. In contrast, injecting the negative-conductance leak current into non-pacemaker neurons in the pyloric network fails to activate the pyloric rhythm or any rhythmic activity in that neuron. 

How does the recovery of the pyloric rhythm by reducing the leak conductance in the pacemaker neurons relate to the known effects of neuromodulators in this network? Much of the neuromodulatory actions on individual pyloric neurons converge, through unknown second messenger mechanisms, to activate a voltage-gated inward current I MI (Swensen and Marder, 2000). We show that the I-V curve of this non-linear inward current can basically be decomposed into two linear components, one with positive-conductance I MI-PL and one with negative conductance I MI-NL. When we introduced one or the other of these two components into the PD neuron using dynamic clamp, we found that I MI-NL led to the recovery of oscillations whereas I MI-PL did not. The recovered oscillation was identical to oscillations produced with the full I-V range of I MI or that produced by bath application of the neuropeptide proctolin. 

Is the linear decomposition of a continuously voltage-dependent current valid? The red and gray traces shown in Figure 6 suggest that the sum of these two linear components, restricted in their voltage range, very closely match the total current I MI over the full voltage range. However, we observed that the linear negative-conductance component (I MI-NL or I NL ) effectively activates the pyloric rhythm only when its reversal potential was negative relative to the membrane potentials trajectory of the pacemaker neurons' slow wave. Reversal potential values more positive than that led to currents that were outward for most of this trajectory and thus did not enhance their excitability (not shown). Conversely, with its positive reversal potential, the positive-conductance linear component (I MI-PL ) only led to a slight depolarization of the baseline voltage but decreased the input resistance of the cells and did not enhance their excitability. 

Electrical loading of a neuron can affect its membrane potential amplitude changes and temporal properties such as spiking and bursting activity. Such electrical loading can occur via two mechanisms: electrical coupling to other cells or changes in the passive membrane conductance. In the pyloric network, electrical coupling of the intrinsically oscillatory AB neuron to non-oscillatory neurons, such as the PD neuron, affects the oscillatory activity of the whole ensemble (Hooper and Marder, 1987; Kepler et al., 1990; Soto-Trevino et al., 2005). If the load produced by the nonoscillatory neurons is large, it would kill the oscillations. A negative conductance in any of the coupled neurons effectively reduces the electric load on the oscillatory cell, leading to an enhanced excitability. Any reduction in leak current has this effect as long as the membrane potential of the cell is more positive than the current's reversal potential (e. g. Sirois et al., 2002; Blethyn et al., 2006; Koizumi and Smith, 2008). Thus, the reduction of the leak current in the PD neuron may lead to the recovery of oscillatory activity because of the strong electrical coupling of the PD neurons to the AB neuron. However, our experiments do not exclude the possibility that the PD neuron itself becomes an intrinsic oscillator under these conditions. 

Previous studies of the pyloric network have shown that blocking outward currents can result in a recovery of oscillations in the AB and PD neurons, although a recovery of the full network activity was not reported (Harris-Warrick and Johnson, 1987). Additionally, modulatory actions by dopamine, that do not target I MI, have also been reported to result in a significant increase of the input resistance in the AB neuron and a simultaneous decrease in the PD neuron (Johnson et al., 1993) and dopamine has been shown to induce oscillatory activity in the isolated AB neuron in lobsters (Johnson et al., 1992). It is therefore possible that recovery of oscillations may occur through mechanisms that include, among other effects, a decrease in leak conductance, but that are independent of I MI. Changes in leak conductance have been shown to be the likely mechanism underlying the excitatory action of several substances acting on metabotropic receptors (Bayliss et al., 1992; Sirois et al., 2002; Blethyn et al., 2006; Brickley et al., 2007; Koizumi and Smith, 2008; Weber et al., 2008). In many if not all these cases, it appears that the consequence of a reduced leak current is to move the voltage of the cell to the activation range of voltage-dependent pacemaker currents. Of course, excitability cannot result from a simple reduction in the endogenous leak current; it is necessary to produce a negative-slope conductance in the I-V curve that creates a regenerative mechanism for membrane depolarization. In the case of the pacemaker pyloric network neurons, we suggest that the negative-conductance component of I MI is itself the pacemaker current since no other current has been reported to be active in the voltage range of this conductance component (Graubard and Hartline, 1991; Golowasch and Marder, 1992a; Soto-Trevino et al., 2005). Although this hypothesis needs to be tested experimentally, it effectively implies that the fundamental pacemaker properties in this network result from the actions of a linear portion of a non-linear current that is activated by neuromodulatory inputs. In other words, it appears that pacemaker activity in this network depends on a linear current, whose essential properties are implemented by the modulator-activated inward current I MI.
section_name: dIscussIon

section_content: We thank Dr. Vahid Tohidi for his initial experiments on this work which were published in abstract form. This research was supported by NIH grants MH64711 (Jorge Golowasch) and MH60605 (Farzan Nadim).
section_name: acknowledgMents