NeuronPRM
Related Projects:  Computational Biology    Sampling-based Planning    Protein Folding    Ligand Binding    RNA Folding  

The brain has extraordinary computational power to represent and interpret complex natural environments. The topology and geometry of the structures in the brain have a decissive role in these natural computations. NeuronPRM is a framework to construct a 3D model of a cortical network, using probabilistic roadmap methods. Although not the usual motion planning problem, our objective of building a network that encodes the pathways of the cortical network is analogous to the PRM objective of constructing roadmaps that contain a representative sample of feasible paths. We represent the network as a large-scale directed graph and use L-systems and statistics data to 'grow' neurons that are morphologically indistinguishable from real neurons. Our ultimate goal is to map and understand the connectivity and geometry of the cortical network.



NeuronPRM:

A Framework for Constructing Cortical Networks

We begin our modeling by partitioning the cerebral cortex into a set of finite elements (FEs). Each FE has a set of geometrical models of recontructed neurons, which are used to generate additional anatomically realistic synthetic neurons. The cortical network is completed by making connections or synapses between the neurons. The construction process borrows from PRM (Probabilistic Roadmap) techniques by using simple, local techniques to understand a much larger unknown space viewed as a robotics problem. Here the workspace is enclosed inside the cortical model, and, within this space, points (neurons) are generated and connected. Each point is a different configuration of a tree-like robot. Unlike a traditional PRM, the tree-like robots have the ability to change shape (i.e., connectivity and number of links), so each roadmap node, like a reconfigurable robot, will represent not only a transformation of the robot's joints but also its topology (joint connectivity). In short, we are trying to solve a robotics problem with an infinite numbers of degrees of freedom.

NeuronPRM framework neurons
NeuronPRM Framework and a resulting synthetic neuron

Neuron Model

Each neuron is represented as a tree of neural components. The first components are the soma and a set of arbors. Each arbor has as children a set of segments who in turn have children that are either junctions (when they branch) or termination marks (when they are terminal segments). Each segment has a list of micro-segments, which may or may not have a spine where synapses take place.

unbounded neuron
Synthetic neuron

Neuron Placement and Abstract Synapses

Neuron-PRM (N-PRM) samples graph nodes and connects nodes using local information. Each vertex in the graph defines a configuration of a neuron and each directed edge represents an abstract connection or synapse between neurons. Here, we define an abstract synpase as a set of real synapses from segments of one neuron to another. Since actual synapses can be quite complex, these abstract synapses help us not only maintain data more easily but also provide a hierarchical representation for searching the cortical network. Following figure shows the result of N-PRM in abstract level.

cortex surface partitioned into set of finite elements cortex populated with neurons cortex populated with neurons (close-up)
Mouse Cortex. (Left) empty cortex with FE decomposition, (Center) Filled with abstract neurons and connections, (Right) A closer view.

Bounding Volumes and Synapses

Statistically each neuron has thousands of synapses. An adult human being has more than one hundred billion neurons and even the cortex of the mouse has more than seventeen billion neurons. Based on current technology, it is not feasible to compute and store such a huge amount of synaptic connections between all possible pairs of neurons. The tradition PRM connection strategy of computing all pairwise distances and attempting to connect the k-closet nodes is not feasible for roadmaps with millions or billions of nodes, each of which consists of thousands of segments. To deal with such large numbers of complex nodes, we define simple distance metrics to reject neuron pairs and find potentioal synapses very quickly. Most neuron pairs are quickly rejected by a filtering test, which checks for intersection of their bounding volumes. Thus, detailed distance computations between the neuronal segments will only be performed for those neuron pairs that pass the bounding volume test. Three bounding volume tests are implemented. They are bounding sphere, boudning box, bounding convex hull.

finite element from mouse cortex filled with abstract neural network neurons represented as bounding box abstract connections between synthetic neurons
Connections in Finite Element or FE: (Left) FE extracted from mouse cortex filled with abstract neural network, (Center) same network with neurons are represented as bounding boxes, (Right) abstract connection between synthetic neurons.

After constructing the roadmaps with abstract synaptic connections, we applied the simple synapse discovery algorithm to identify real synapses. The running time for discovery is linear in the number of edges of the input roadmap, so the convex hull is the fastest.

We determined that the approximate convex hull bounding volume was the fastest overall (considering both the cost of the abstract and the actual synapse discovery phases); its drawback is the increased storage requirements for the convex hull as compared to the bounding sphere or box.

granule neuron generated throug L-System bounding box and convex hull for granule neuron
Bounding Volumes. (Left) a Granule neuron generated through an L-System, (Right) Bounding box and Convex Hull for this neuron.

Related Publications

Neuron PRM: A Framework for Constructing Cortical Networks, Jyh-Ming Lien, Marco Morales, Nancy M. Amato, Neurocomputing, Vol: 54, pp. 191-197, Jun 2003. DOI: 10.1016/S0925-2312(02)00728-2
Keywords: Machine Learning, Sampling-Based Motion Planning
Links : [Published]

BibTex

@article{LIEN2003191,
title = "Neuron PRM: a framework for constructing cortical networks",
journal = "Neurocomputing",
volume = "52-54",
pages = "191 - 197",
year = "2003",
note = "Computational Neuroscience: Trends in Research 2003",
issn = "0925-2312",
doi = "https://doi.org/10.1016/S0925-2312(02)00728-2",
url = "http://www.sciencedirect.com/science/article/pii/S0925231202007282",
author = "Jyh-Ming Lien and Marco Morales and Nancy M. Amato",
keywords = "Cortical networks, PRM, BTS, L-system, Rectangle tree",
abstract = "The brain's extraordinary computational power to represent and interpret complex natural environments is essentially determined by the topology and geometry of the brain's architectures. We present a framework to construct cortical networks which borrows from probabilistic roadmap methods developed for robotic motion planning. We abstract the network as a large-scale directed graph, and use L-systems and statistical data to ‘grow’ neurons that are morphologically indistinguishable from real neurons. We detect connections (synapses) between neurons using geometric proximity tests."
}


Abstract

The brain's extraordinary computational power to represent and interpret complex natural environments is essentially determined by the topology and geometry of the brain's architectures. We present a framework to construct cortical networks which borrows from probabilistic roadmap methods developed for robotic motion planning. We abstract the network as a large-scale directed graph, and use L-systems and statistical data to ‘grow’ neurons that are morphologically indistinguishable from real neurons. We detect connections (synapses) between neurons using geometric proximity tests.