Efficient Integration of Coupled Electrical-Chemical Systems in Multiscale Neuronal Simulations.
Multiscale modeling and simulations in neuroscience is gaining scientific attention due to its growing importance and unexplored capabilities. For instance, it can help to acquire better understanding of biological phenomena that have important features at multiple scales of time and space. This includes synaptic plasticity, memory formation and modulation, homeostasis. There are several ways to organize multiscale simulations depending on the scientific problem and the system to be modeled. One of the possibilities is to simulate different components of a multiscale system simultaneously and exchange data when required. The latter may become a challenging task for several reasons. First, the components of a multiscale system usually span different spatial and temporal scales, such that rigorous analysis of possible coupling solutions is required. Then, the components can be defined by different mathematical formalisms. For certain classes of problems a number of coupling mechanisms have been proposed and successfully used. However, a strict mathematical theory is missing in many cases. Recent work in the field has not so far investigated artifacts that may arise during coupled integration of different approximation methods. Moreover, in neuroscience, the coupling of widely used numerical fixed step size solvers may lead to unexpected inefficiency. In this paper we address the question of possible numerical artifacts that can arise during the integration of a coupled system. We develop an efficient strategy to couple the components comprising a multiscale test problem in neuroscience. We introduce an efficient coupling method based on the second-order backward differentiation formula (BDF2) numerical approximation. The method uses an adaptive step size integration with an error estimation proposed by Skelboe (2000). The method shows a significant advantage over conventional fixed step size solvers used in neuroscience for similar problems. We explore different coupling strategies that define the organization of computations between system components. We study the importance of an appropriate approximation of exchanged variables during the simulation. The analysis shows a substantial impact of these aspects on the solution accuracy in the application to our multiscale neuroscientific test problem. We believe that the ideas presented in the paper may essentially contribute to the development of a robust and efficient framework for multiscale brain modeling and simulations in neuroscience.
Science for Life Laboratory, Computational Science and Technology, School of Computer Science and Communication, KTH Royal Institute of TechnologyStockholm, Sweden; National Centre for Biological SciencesBangalore, India; Manipal UniversityManipal, India.
Ekaterina
E
Brocke
National Centre for Biological Sciences Bangalore, India.
Upinder S
US
Bhalla
PDC Center for High-Performance Computing, KTH Royal Institute of TechnologyStockholm, Sweden; International Neuroinformatics Coordinating Facility, Karolinska InstituteStockholm, Sweden.
Mikael
M
Djurfeldt
Science for Life Laboratory, Computational Science and Technology, School of Computer Science and Communication, KTH Royal Institute of TechnologyStockholm, Sweden; Science for Life Laboratory, Department of Numerical Analysis and Computer Science, Stockholm UniversityStockholm, Sweden; Department of Neuroscience, Karolinska InstituteStockholm, Sweden.
Jeanette
J
Hellgren Kotaleski
Department of Mathematics, School of Engineering Sciences, KTH Royal Institute of Technology Stockholm, Sweden.
Michael
M
Hanke
adaptive time step integration; backward differentiation formula; co-simulation; coupled integration; coupled system; multiscale modeling; multiscale simulation; parallel numerical integration
Prog Mol Biol Transl Sci. 2014;123:351-86
24560151
Int J Biomed Comput. 1989 Mar;24(1):55-68
2714879
Biol Cybern. 2008 Nov;99(4-5):427-41
19011929
Prog Biophys Mol Biol. 2011 Oct;107(1):169-82
21729716
PLoS One. 2013 Jul 03;8(7):e66811
23843966
Neural Netw. 2011 Nov;24(9):943-9
21616637
J Physiol. 1952 Aug;117(4):500-44
12991237
Int J Biomed Comput. 1984 Jan-Feb;15(1):69-76
6698635
Biophys J. 2007 Sep 15;93(6):1847-57
17496042