Predictive Modeling of Bioelectric Activity on Mammalian Multilayered Neuronal Structures in the Presence of Supraphysiological Electric Fields

Investigators
Theodore W. Berger, Gianluca Lazzi
Contact info (email)
lazzi@usc.edu; berger@usc.edu; jbouteil@usc.edu
1. Define context(s)
identify/explore new therapies
reveal new biological insights
other
Current Conformance Level / Target Conformance Level
10
Primary goal of the model/tool/database

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.

Biological domain of the model
Nervous System (hippocampus, retina)
Structure(s) of interest in the model
From biomolecular processes to neurons to network to tissue
Spatial scales included in the model
from nanometers to millimeters
Time scales included in the model
microseconds to seconds
2. Data for building and validating the model
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: ________________________
3. Validate within context(s)
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
4. Limitations
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
5. Version control
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
6. Documentation
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
7. Dissemination
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
8. Independent reviews
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?
9. Test competing implementations
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!
10. Conform to standards
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
11. (optional) Additional information to support items 1-10

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).