The end goal of this multiscale modeling research is to bridge the gap existing between three-dimensional, full-wave, macro-modeling of electrical and magnetic biointeractions (global modeling) and cellular-level modeling strategies. This research effort aims to predict spatio-temporal distributions of active neurons based on current densities created by multi-electrode electrical stimulation, and dependent upon a set of "core models" of molecular (receptor-channel kinetics), synaptic, neuron, and multi-neuron activity. These models and their inputs and outputs must be integrated into a global model of the extracellular media/matrix including relevant multi-electrode arrays. Successful modeling at these levels will allow hypotheses about space-time patterns of electrical stimulation to produce predictions about the number and distribution of activated inputs (based on known spatial distributions of afferent axons). The linked molecular, synaptic, neuron, multi-neuron, and global model will provide the basis for emerging predictions of the spatio-temporal distribution of active neurons and thus, the spatio-temporal distributions of spike train activity that encode all information in the nervous system.
| Data for building the model | Published? | Private? | How is credibility checked? | Current Conformance Level / Target Conformance Level |
|---|---|---|---|---|
| in vitro (primary cells cell, lines, etc.) | yes | yes | cross check with other results | good |
| ex vivo (excised tissues) | ||||
| in vivo pre-clinical (lower-level organism or small animal) | yes | yes | cross check with other results | good |
| in vivo pre-clinical (large animal) | ||||
| Human subjects/clinical | ||||
| Other: ________________________ |
| Data for validating the model | Published? | Private? | How is credibility checked? | Current Conformance Level / Target Conformance Level |
|---|---|---|---|---|
| in vitro (primary cells cell, lines, etc.) | yes | yes | cross check with other results | good |
| ex vivo (excised tissues) | ||||
| in vivo pre-clinical (lower-level organism or small animal) | yes | yes | cross check with other results | good |
| in vivo pre-clinical (large animal) | ||||
| Human subjects/clinical | ||||
| Other: ________________________ |
| Who does it? | When does it happen? | How is it done? | Current Conformance Level / Target Conformance Level | |
|---|---|---|---|---|
| Verification | student/postdoc/faculty | throughout | through systematic and/or manual testing | good |
| Validation | student/postdoc/faculty | throughout | through systematic and/or manual testing | good |
| Uncertainty quantification | ||||
| Sensitivity analysis | student/postdoc/faculty | sporadically | SA studies | good |
| Other:__________ | ||||
| Additional Comments |
| Disclaimer statement (explain key limitations) | Who needs to know about this disclaimer? | How is this disclaimer shared with that audience? | Current Conformance Level / Target Conformance Level |
|---|---|---|---|
| mechanistic models derived from fixed protocols and experiments which inherently limit prediction levels to in-range protocols/values; extrapolation may lead to higher prediction errors | models users | clearly laid out in publications | good |
| data-driven models also have strong predictive power within bounds of data used to construct them | models users | clearly laid out in publications | good |
| Current Conformance Level / Target Conformance Level |
|---|
| Good - Our group uses subversion versioning system. All changes to model structure and parameters values are recorded and tracked. The suite of tests performed to verify/validate the models are also versioned. |
| Naming Conventions? | Repository? | Code Review? | |
|---|---|---|---|
| individual modeler | no | Subversion / Github | yes |
| within the lab | no | Subversion / Github | yes |
| collaborators | no | Subversion / Github | not yet for outside members |
| Current Conformance Level / Target Conformance Level | |
|---|---|
| Code commented? | Code is commented to improve code readability |
| Scope and intended use described? | this is described globally, and within indivudual files and functions |
| User’s guide? | Currently in manuscript and on GitHub for finalized and shared models |
| Developer’s guide? | Guidance is mostly done through comments in the code |
| Current Conformance Level / Target Conformance Level |
|---|
| good |
| Target Audience(s): | “Inner circle” | Scientific community | Public |
|---|---|---|---|
| Simulations | yes | yes | yes |
| Models | yes | yes - although not all models are currently public | yes - although not all models are currently public |
| Software | yes | yes | yes |
| Results | yes | yes | yes |
| Implications of results | yes | yes | yes |
| Current Conformance Level / Target Conformance Level |
|---|
| As part of the BioMedical Simulation Resource (BMSR), our research directions and efforts have undergone third-party evaluation. More granular models evaluation will have to take place according to grant requirements and the proposed timeline. |
| Reviewer(s) name & affiliation: | |
|---|---|
| When was review performed? | |
| How was review performed and outcomes of the review? |
| Current Conformance Level / Target Conformance Level |
|---|
| We compared our IO models to multiple IO and mechanistic implementations, which helped us quantitatively determine our model's capabilities and limitations (details in publications). The starting point for our mechanistic neuron models are previously published models which are cited in our publications. In adapting these models to our platform, tests are performed to verify that our models are capable of recreating the figures associated with the original models. At the network level, there are yet to be comparable models to which we can compare. This is the same for modeling the effects of extracellular stimulation. Notably, our efforts to discretize the extracellular space led to two meshing methods (hexahedral and tetrahedral) which enabled us to mutually validate their results. |
| Yes or No (briefly summarize) | |
|---|---|
| Were competing implementations tested? | when possible (biomolecular level, neuron level), but harder to do when competing implementations do not exist |
| Did this lead to model refinement or improvement? | Yes! |
| Current Conformance Level / Target Conformance Level |
|---|
| relatively good (depending on the model of interest). Our models are developed using best coding practices and operating procedures. Notably, the neuron and network models do not yet conform to a standard mode description language such as NeuroML. We plan on converting the code to this standard by 2021. On the other hand, the mechanistic models of biomolecular processes are all coded in the SBML standard (Systems Biology Markup Language). |
| Yes or No (briefly summarize) | |
|---|---|
| Are there operating procedures, guidelines, or standards for this type of multiscale modeling? | NeuroML, SBML |
| How do your modeling efforts conform? | relatively good |
For rule 1: The goal of the project is to predict spatio-temporal distributions of active neurons based on current densities created by multi-electrode electrical stimulation, and dependent upon a set of "core models" of molecular (receptor-channel kinetics), synaptic, neuron, and multi-neuron activity. The model is intended to be used for (1) theoretical research that links multiple physical scales to network activity or a cognitive process (e.g., spatial navigation or memory), (2) electrode and stimulus design, and (3) drug discovery.
For rule 2: All model parameters were constrained using experimental data where available, e.g. unitary postsynaptic potentials for the appropriate pre-post neuron pairs, neuron electrophysiology specific for a neuron type, axon distributions, etc. The models validated against the appropriate experimental data, e.g. evoked potentials due to extracellular stimulus. These data and their references are described in all relevant publications. In certain cases, model development using prior published models, not from our research groups, as starting points and are referenced in publications. Their structure and parameters values are available in their respective publications.
For rule 3: All models are evaluated within their respective contexts and corresponding levels of granularity (See List of Planned Activities Table, section Systemic Calibration and Evaluation). For example, parameters for our abstract (data-driven) models derived from mechanistic models are estimated while minimizing error between output from the mechanistic (detailed and computationally complex) model and output from the abstract model using a 'training' dataset; validation is performed using a separate 'validation' dataset.
For rule 4: Mechanistic models drive their predictive powers from the test protocols with which they were initially calibrated and tested; significant variations in the simulation protocols may consequently affect model predictability levels. The specific protocols using which the models were constrained are described in publications; further limitations are described in the discussion sections. In addition, for models already shared with the community, these publications are explicitly mentioned on the webpages on which the models are shared.
For rule 7: Detailed models out of prototyping phase and used for development of input-output models were published and their code shared on open-source database (ModelDB). IO model structure was also described in multiple publications. Other IO models will be similarly shared as they become available. Notably, the admittance/neuron bridge model will be made accessible within the next year. The ROOTS algorithm is already freely available on GitHub and through the PyPi library.
For rule 9: We compared our IO models to multiple IO and mechanistic implementations, which helped us quantitatively determine our model's capabilities and limitations (details in publications). The starting point for our mechanistic neuron models are previously published models which are cited in our publications. In adapting these models to our platform, tests are performed to verify that our models are capable of recreating the figures associated with the original models. At the network level, there are yet to be comparable models to which we can compare. This is the same for modeling the effects of extracellular stimulation. Notably, our efforts to discretize the extracellular space led to two meshing methods (hexahedral and tetrahedral) which enabled us to mutually validate their results.
For rule 10: Our models are developed using best coding practices and operating procedures. Notably, the neuron and network models do not yet conform to a standard mode description language such as NeuroML. We plan on converting the code to this standard by 2021. On the other hand, the mechanistic models of biomolecular processes are all coded in the SBML standard (Systems Biology Markup Language).