“Can we predict synchrony and asynchrony in networks coupled by multiple dendritic gap junctions?”

1. “Can we predict synchrony and asynchrony in networks coupled by multiple dendritic gap junctions?” Frances K. Skinner Toronto Western Research Institute University Health Network and University of Toronto New York University April 13, 2008

2. From Scholarpedia:Mathematical Biology article ofFrank Hoppensteadt • “highly interdisciplinary nature” • “barriers to collaborations between mathematicians and biologists” • “a shift from mathematical analysis to computer simulation due mostly to improvements in computer power and accessibility.. With the shift being made possible to include more information in models and still derive useful insights from them.” Especially in neuroscience with all the details being uncovered, increasingly sophisticated techniques etc. these comments are very timely. With increasing specialization and interdisciplinarity and potential moving apart of mathematical and biological sciences (or separation of viewpoints) we organized a Theoretical Neuroscience Minisymposium at 2006 Society for Neuroscience meeting, one of the aims being to help counteract this.

3. Interneuron Heterogeneity Domain-specific innervation of hippocampal interneurons apical basal purple –laminae where axonal arbor typically extends turquoise indicates that other interneurons rather than principal cells are targets Different types of interneurons containing calcium-binding proteins and neuropeptides McBain and Fisahn NRN 2001

4. General challenge – how to best consider various neurobiological details. • Specific challenge – understanding the contribution of electrical coupling in different contexts. • Outline of Talk: background, discussion of some of our previous work, and then get to question posed for this talk.

5. “Can we predict synchrony and asynchrony in networks coupled by multiple dendritic gap junctions?” Frances K. Skinner Toronto Western Research Institute University Health Network and University of Toronto New York University April 13, 2008

6. Acknowledgements Tariq Zahid Fernanda Saraga, Leo Ng NSERC of Canada Computing support – RIS of UHN

7. The hippocampus (part of medial temporal lobe) is an intensely studied region of the brain because: • It is associated with memory and learning (i.e., LTP, LTD), epileptic seizures, and neurogenesis. • It exhibits a wide range of population rhythmic activity patterns (<1 to >200 Hz) that are associated with various behavioural states. • It is amenable to experiment, retaining its synaptic circuitry and thus population activities in the slice.

8. EEG activities of mouse hippocampus Theta-Gamma SPW-ripples Electrode location Sharp wave-ripples

9. Spontaneous Rhythmic Field Potentials (SRFPs) (Liang Zhang’s lab; Wu et al. J Physiol 2002, J Neurophysiol 2005) Gillis et al., J. Neurosci. Meth. (2005)

10. Blockade of field rhythms and pyramidal IPSPs by GABA-A receptor antagonist rhythmic activities also dependent on electrical coupling (gap junctions)

11. Background Electrical coupling (i.e., gap junctions) is present in much of the mammalian brain (e.g., inferior olive, striatum, neocortex, hippocampus, retina, thalamus). In particular, gap junctions occur between inhibitory cells, often of the same type, and can be located at sites quite distant (> 200 μm) from the soma.

12. Interneurons represent 10-20% of the neuronal population but may provide the precise temporal structure necessary for ensembles of neurons to perform specific functions. - Buzsáki and Chrobak, 1995

13. Interneuron Heterogeneity Domain-specific innervation of hippocampal interneurons apical basal purple –laminae where axonal arbor typically extends turquoise indicates that other interneurons rather than principal cells are targets Different types of interneurons containing calcium-binding proteins and neuropeptides McBain and Fisahn NRN 2001

14. Background Gap junctions located far from cell bodies, at non-proximal sites (basket cells in hippocampus) Gap junctions can be modulated Inhibitory cells have active dendrites, spikes can be generated in dendrites From Fukuda & Kosaka, J Neurosci 2000

15. Dendrodendritic Gap Junctions Fukuda & Kosaka, J Neurosci 2000

16. 50 mV 100 ms Model (Hippocampal Basket Cell) Passive dendrites 372-compartment model developed in NEURON Morphology from Gulyas et al.(1999) Saraga et al., J Neurophysiol 2006

17. WB used for kinetic model basis, Martina and Jonas (1997), Martina et al (1998) used as conductance value basis and spike characteristics and electrophysiological responses from Morin et al. (1996) and van Hooft et al. (2000).

18. Vout=0.14 s d Vin =0.25 d “Reduced” 3-compartment model based on matching electrotonic length from soma (Vout) Vout Vin 0.12 1.2 Electrotonic distance from soma (L) Electrotonic distance to soma (L) 0.4 0.04 0 100 300 500 0 100 300 500 Anatomical distance to soma (mm) Anatomical distance from soma (mm)

19. s d d s d d 50 mV 10 ms (Distal) phase response curves (PRCs) Voltage along dendrite Using the reduced model geometry 0.5% PRCs calculated using XPPAUT (Ermentrout, 2002) 1% Phase Shift 1.5% 10% Phase

20. 28 26 24 23 20 18 16 14 75 50 25 0 LOW MEDIUM HIGH % Phase Lag Intrinsic Frequency (Hz) 0.1 1 10 100 % Active Predicted Network Dynamics • Weakly coupled oscillator • theory used to define three • different dynamic regions • LOW, MEDIUM, HIGH • that refer to PRCs with • particular characteristics • (e.g., negative PRCs • for MEDIUM) Phase lags determined from interaction functions calculated using XPPAUT (Ermentrout, 2002)

21. Simulations confirm theoretical predictions 75 50 25 0 LOW MED HIGH % PhaseLag Intrinsic Frequency(Hz) ggap 0.1 1 10 100 % Active % Active Results % Phase Lag • “Weak coupling” is about 10 pS (comparing predicted and simulated) • Compare full and reduced model phase lag values to • “define” synchronous and asynchronous • Synchronous is 10% or less phase lag, asynchronous otherwise

22. Cell 1: 15% basal attenuation, 2% apical attenuation Cell 2: 8% basal attenuation, 8% apical attenuation Cell 3: 6% basal attenuation, 14% apical attenuation apical apical CELL 3 basal basal

23. Cell 1, apical coupling (multistability)

24. Beyond weak coupling

25. Beyond weak coupling

26. Beyond weak coupling

27. Discussion and Conclusions • PRC skewness quantifications can be used to predict whether synchronous or asynchronous modes occur in electrically coupled basket cells. • Averaged PRCs can be used to predict modes for coupling at multiple sites. • Predictions cannot be made under all circumstances and multistability can occur. • Different apical and basal attenuation (due to different channel densities) allow more ‘robust’ asynchrony to occur with coupling on the more attenuated dendritic side. • Network couplings that produce asynchrony (as compared to synchrony) with weak coupling encompass more dynamic richness (i.e., range of possible phase lags) with gap junction conductance changes. • Thus, gap junction coupling may be able to tune networks in and out of synchronous activities if asynchrony with weak coupling is predicted.

28. the end • Thank you!