Data-Driven Predictive Modeling of Neuronal Dynamics using Long Short-Term Memory

Modeling brain dynamics to better understand and control complex behaviors underlying various cognitive brain functions are of interests to engineers, mathematicians, and physicists from the last several decades. With a motivation of developing computationally efficient models of brain dynamics to use in designing control-theoretic neurostimulation strategies, we have developed a novel data-driven approach in a long short-term memory (LSTM) neural network architecture to predict the temporal dynamics of complex systems over an extended long time-horizon in future. In contrast to recent LSTM-based dynamical modeling approaches that make use of multi-layer perceptrons or linear combination layers as output layers, our architecture uses a single fully connected output layer and reversed-order sequence-to-sequence mapping to improve short time-horizon prediction accuracy and to make multi-timestep predictions of dynamical behaviors. We demonstrate the efficacy of our approach in reconstructing the regular spiking to bursting dynamics exhibited by an experimentally-validated 9-dimensional Hodgkin-Huxley model of hippocampal CA1 pyramidal neurons. Through simulations, we show that our LSTM neural network can predict the multi-time scale temporal dynamics underlying various spiking patterns with reasonable accuracy. Moreover, our results show that the predictions improve with increasing predictive time-horizon in the multi-timestep deep LSTM neural network.

[1]  R. Hari,et al.  Dynamics of brain activation during picture naming , 1994, Nature.

[2]  Chethan Pandarinath,et al.  Inferring single-trial neural population dynamics using sequential auto-encoders , 2017 .

[3]  Sergey Levine,et al.  Neural Network Dynamics for Model-Based Deep Reinforcement Learning with Model-Free Fine-Tuning , 2017, 2018 IEEE International Conference on Robotics and Automation (ICRA).

[4]  Addolorata Marasco,et al.  On the mechanisms underlying the depolarization block in the spiking dynamics of CA1 pyramidal neurons , 2012, Journal of Computational Neuroscience.

[5]  Johan A. K. Suykens,et al.  Nonlinear system identification using neural networks , 1996 .

[6]  Marie-Therese Horstmann,et al.  Role of axonal NaV1.6 sodium channels in action potential initiation of CA1 pyramidal neurons. , 2008, Journal of neurophysiology.

[7]  Bartlett W. Mel,et al.  Arithmetic of Subthreshold Synaptic Summation in a Model CA1 Pyramidal Cell , 2003, Neuron.

[8]  Jürgen Schmidhuber,et al.  Long Short-Term Memory , 1997, Neural Computation.

[9]  Constantinos Siettos,et al.  Multiscale modeling of brain dynamics: from single neurons and networks to mathematical tools , 2016, Wiley interdisciplinary reviews. Systems biology and medicine.

[10]  Yoshua Bengio,et al.  Understanding the difficulty of training deep feedforward neural networks , 2010, AISTATS.

[11]  PAUL J. WERBOS,et al.  Generalization of backpropagation with application to a recurrent gas market model , 1988, Neural Networks.

[12]  Wulfram Gerstner,et al.  Firing patterns in the adaptive exponential integrate-and-fire model , 2008, Biological Cybernetics.

[13]  Wen Yu,et al.  Non-linear system modeling using LSTM neural networks , 2018 .

[14]  Michael C. Mozer,et al.  A Focused Backpropagation Algorithm for Temporal Pattern Recognition , 1989, Complex Syst..

[15]  Stefen Hui,et al.  Application of feedforward neural networks to dynamical system identification and control , 1993, IEEE Trans. Control. Syst. Technol..

[16]  Yu Wang,et al.  A new concept using LSTM Neural Networks for dynamic system identification , 2017, 2017 American Control Conference (ACC).

[17]  Terrence J. Sejnowski,et al.  Gradient Descent for Spiking Neural Networks , 2017, NeurIPS.

[18]  Razvan Pascanu,et al.  On the difficulty of training recurrent neural networks , 2012, ICML.

[19]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[20]  Petros Koumoutsakos,et al.  Data-driven forecasting of high-dimensional chaotic systems with long short-term memory networks , 2018, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[21]  Krasimira Tsaneva-Atanasova,et al.  A unified model of CA1/3 pyramidal cells: an investigation into excitability. , 2011, Progress in biophysics and molecular biology.

[22]  Tim Gollisch,et al.  Modeling Single-Neuron Dynamics and Computations: A Balance of Detail and Abstraction , 2006, Science.

[23]  David Golomb,et al.  Contribution of persistent Na+ current and M-type K+ current to somatic bursting in CA1 pyramidal cells: combined experimental and modeling study. , 2006, Journal of neurophysiology.

[24]  N. Spruston,et al.  Synapse Distribution Suggests a Two-Stage Model of Dendritic Integration in CA1 Pyramidal Neurons , 2009, Neuron.

[25]  Harald Haas,et al.  Harnessing Nonlinearity: Predicting Chaotic Systems and Saving Energy in Wireless Communication , 2004, Science.

[26]  Addolorata Marasco,et al.  Fast and accurate low-dimensional reduction of biophysically detailed neuron models , 2012, Scientific Reports.

[27]  Ross A. Knepper,et al.  DeepMPC: Learning Deep Latent Features for Model Predictive Control , 2015, Robotics: Science and Systems.

[28]  Naresh K. Sinha,et al.  Modeling Robot Dynamics Using Dynamic Neural Networks , 1997 .

[29]  M. Breakspear Dynamic models of large-scale brain activity , 2017, Nature Neuroscience.

[30]  Surya Ganguli,et al.  SuperSpike: Supervised Learning in Multilayer Spiking Neural Networks , 2017, Neural Computation.

[31]  Ganapati Panda,et al.  Identification of nonlinear dynamic systems using functional link artificial neural networks , 1999, IEEE Trans. Syst. Man Cybern. Part B.

[32]  Wulfram Gerstner,et al.  How Good Are Neuron Models? , 2009, Science.

[33]  Quoc V. Le,et al.  Sequence to Sequence Learning with Neural Networks , 2014, NIPS.

[34]  Karl J. Friston,et al.  A Hierarchy of Time-Scales and the Brain , 2008, PLoS Comput. Biol..

[35]  Maurizio Corbetta,et al.  The human brain is intrinsically organized into dynamic, anticorrelated functional networks. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[36]  Ivan Raikov,et al.  Interneuronal mechanisms of hippocampal theta oscillations in a full-scale model of the rodent CA1 circuit , 2016, eLife.

[37]  Frances K. Skinner,et al.  Simple, biologically-constrained CA1 pyramidal cell models using an intact, whole hippocampus context , 2014, F1000Research.

[38]  P J Webros BACKPROPAGATION THROUGH TIME: WHAT IT DOES AND HOW TO DO IT , 1990 .

[39]  Nelson Spruston,et al.  Dendritic sodium spikes are required for long-term potentiation at distal synapses on hippocampal pyramidal neurons , 2015, eLife.

[40]  C A Bailer-Jones,et al.  A recurrent neural network for modelling dynamical systems. , 1998, Network.

[41]  Stephen A. Billings,et al.  International Journal of Control , 2004 .

[42]  Karthik Duraisamy,et al.  Long-time predictive modeling of nonlinear dynamical systems using neural networks , 2018, Complex..

[43]  Wulfram Gerstner,et al.  Neuronal Dynamics: From Single Neurons To Networks And Models Of Cognition , 2014 .

[44]  Steven Lake Waslander,et al.  Multistep Prediction of Dynamic Systems With Recurrent Neural Networks , 2018, IEEE Transactions on Neural Networks and Learning Systems.

[45]  Erin C. McKiernan,et al.  CA1 pyramidal cells have diverse biophysical properties, affected by development, experience, and aging , 2017, PeerJ.