Toward blueprints for network architecture, biophysical dynamics, and signal transduction

We review mathematical aspects of biophysical dynamics, signal transduction and network architecture that have been used to uncover functionally significant relations between the dynamics of single neurons and the networks they compose. We focus on examples that combine insights from these three areas to expand our understanding of systems neuroscience. These range from single neuron coding to models of decision making and electrosensory discrimination by networks and populations, as well as coincidence detection in pairs of dendrites and the dynamics of large networks of excitable dendritic spines. We conclude by describing some of the challenges that lie ahead as the applied mathematics community seeks to provide the tools that will ultimately underpin systems neuroscience.

Release currents of IP₃ receptor channel clusters and concentration profiles

We simulate currents and concentration profiles generated by Ca2+ release from the endoplasmic reticulum (ER) to the cytosol through IP3 receptor channel clusters. Clusters are described as conducting pores in the lumenal membrane with a diameter from 6 nm to 36 nm. The endoplasmic reticulum is modeled as a disc with a radius of 1–12 mm and an inner height of 28 nm. We adapt the dependence of the currents on the trans Ca2+ concentration (intralumenal) measured in lipid bilayer experiments to the cellular geometry. Simulated currents are compared with signal mass measurements in Xenopus oocytes. We find that release currents depend linearly on the concentration of free Ca2+ in the lumen. The release current is approximately proportional to the square root of the number of open channels in a cluster. Cytosolic concentrations at the location of the cluster range from 25 μM to 170 μM. Concentration increase due to puffs in a distance of a few micrometers from the puff site is found to be in the nanomolar range. Release currents decay biexponentially with timescales of < 1 s and a few seconds. Concentration profiles decay with timescales of 0.125–0.250 s upon termination of release.