Generic modeling of chemotactic based self-wiring of neural networks

Biased Random Walk Model of Neuronal Dynamics on Substrates with Periodic Geometrical Patterns

Neuronal networks are complex systems of interconnected neurons responsible for transmitting and processing information throughout the nervous system. The building blocks of neuronal networks consist of individual neurons, specialized cells that receive, process, and transmit electrical and chemical signals throughout the body. The formation of neuronal networks in the developing nervous system is a process of fundamental importance for understanding brain activity, including perception, memory, and cognition. To form networks, neuronal cells extend long processes called axons, which navigate toward other target neurons guided by both intrinsic and extrinsic factors, including genetic programming, chemical signaling, intercellular interactions, and mechanical and geometrical cues. Despite important recent advances, the basic mechanisms underlying collective neuron behavior and the formation of functional neuronal networks are not entirely understood. In this paper, we present a combined experimental and theoretical analysis of neuronal growth on surfaces with micropatterned periodic geometrical features. We demonstrate that the extension of axons on these surfaces is described by a biased random walk model, in which the surface geometry imparts a constant drift term to the axon, and the stochastic cues produce a random walk around the average growth direction. We show that the model predicts key parameters that describe axonal dynamics: diffusion (cell motility) coefficient, average growth velocity, and axonal mean squared length, and we compare these parameters with the results of experimental measurements. Our findings indicate that neuronal growth is governed by a contact-guidance mechanism, in which the axons respond to external geometrical cues by aligning their motion along the surface micropatterns. These results have a significant impact on developing novel neural network models, as well as biomimetic substrates, to stimulate nerve regeneration and repair after injury.

In silico framework to inform the design of repair constructs for peripheral nerve injury repair

Peripheral nerve injuries affect millions of people per year and cause loss of sensation and muscle control alongside chronic pain. The most severe injuries are treated through a nerve autograft; however, donor site morbidity and poor outcomes mean alternatives are required. One option is to engineer nerve replacement tissues to provide a supportive microenvironment to encourage nerve regeneration as an alternative to nerve grafts. Currently, progress is hampered due to a lack of consensus on how to arrange materials and cells in space to maximize rate of regeneration. This is compounded by a reliance on experimental testing, which precludes extensive investigations of multiple parameters due to time and cost limitations. Here, a computational framework is proposed to simulate the growth of repairing neurites, captured using a random walk approach and parameterized against literature data. The framework is applied to a specific scenario where the engineered tissue comprises a collagen hydrogel with embedded biomaterial fibres. The size and number of fibres are optimized to maximize neurite regrowth, and the robustness of model predictions is tested through sensitivity analyses. The approach provides an in silico tool to inform the design of engineered replacement tissues, with the opportunity for further development to multi-cue environments.

Mathematical models of neuronal growth

The establishment of a functioning neuronal network is a crucial step in neural development. During this process, neurons extend neurites-axons and dendrites-to meet other neurons and interconnect. Therefore, these neurites need to migrate, grow, branch and find the correct path to their target by processing sensory cues from their environment. These processes rely on many coupled biophysical effects including elasticity, viscosity, growth, active forces, chemical signaling, adhesion and cellular transport. Mathematical models offer a direct way to test hypotheses and understand the underlying mechanisms responsible for neuron development. Here, we critically review the main models of neurite growth and morphogenesis from a mathematical viewpoint. We present different models for growth, guidance and morphogenesis, with a particular emphasis on mechanics and mechanisms, and on simple mathematical models that can be partially treated analytically.

Interactions in active colloids

The past two decades have seen a remarkable progress in the development of synthetic colloidal agents which are capable of creating directed motion in an unbiased environment at the microscale. These self-propelling particles are often praised for their enormous potential to self-organize into dynamic nonequilibrium structures such as living clusters, synchronized super-rotor structures or self-propelling molecules featuring a complexity which is rarely found outside of the living world. However, the precise mechanisms underlying the formation and dynamics of many of these structures are still barely understood, which is likely to hinge on the gaps in our understanding of how active colloids interact. In particular, besides showing comparatively short-ranged interactions which are well known from passive colloids (Van der Waals, electrostatic etc), active colloids show novel hydrodynamic interactions as well as phoretic and substrate-mediated 'osmotic' cross-interactions which hinge on the action of the phoretic field gradients which are induced by the colloids on other colloids in the system. The present article discusses the complexity and the intriguing properties of these interactions which in general are long-ranged, non-instantaneous, non-pairwise and non-reciprocal and which may serve as key ingredients for the design of future nonequilibrium colloidal materials. Besides providing a brief overview on the state of the art of our understanding of these interactions a key aim of this review is to emphasize open key questions and corresponding open challenges.

A stochastic framework to model axon interactions within growing neuronal populations

The confined and crowded environment of developing brains imposes spatial constraints on neuronal cells that have evolved individual and collective strategies to optimize their growth. These include organizing neurons into populations extending their axons to common target territories. How individual axons interact with each other within such populations to optimize innervation is currently unclear and difficult to analyze experimentally in vivo. Here, we developed a stochastic model of 3D axon growth that takes into account spatial environmental constraints, physical interactions between neighboring axons, and branch formation. This general, predictive and robust model, when fed with parameters estimated on real neurons from the Drosophila brain, enabled the study of the mechanistic principles underlying the growth of axonal populations. First, it provided a novel explanation for the diversity of growth and branching patterns observed in vivo within populations of genetically identical neurons. Second, it uncovered that axon branching could be a strategy optimizing the overall growth of axons competing with others in contexts of high axonal density. The flexibility of this framework will make it possible to investigate the rules underlying axon growth and regeneration in the context of various neuronal populations.

Understanding how neuronal cells establish complex circuits with specific functions within a developing brain is a major current challenge. Over the last past years, enormous progress has been done to precisely resolve brain anatomy and to dissect the mechanisms controlling the establishment of precise neuronal networks. However, due to the extreme complexity of the brain, it is still experimentally difficult to investigate in vivo how neurons interact with each other and with their physical environments to innervate target territories during development. Here, we have developed a framework that integrates a dynamic 3D mathematical model of single axonal growth with parameters estimated from neurons grown in vivo and simulations of entire populations of growing axons. The emergent properties of our model enable the study of the mechanistic principles underlying the growth of axonal population in developing brains. Specifically, our results highlight the impact of mechanical interactions on both individual and collective axon growth, and uncover how branching regulate this process.

Neural electrical activity and neural network growth

The development of central and peripheral neural system depends in part on the emergence of the correct functional connectivity in its input and output pathways. Now it is generally accepted that molecular factors guide neurons to establish a primary scaffold that undergoes activity-dependent refinement for building a fully functional circuit. However, a number of experimental results obtained recently shows that the neuronal electrical activity plays an important role in the establishing of initial interneuronal connections. Nevertheless, these processes are rather difficult to study experimentally, due to the absence of theoretical description and quantitative parameters for estimation of the neuronal activity influence on growth in neural networks. In this work we propose a general framework for a theoretical description of the activity-dependent neural network growth. The theoretical description incorporates a closed-loop growth model in which the neural activity can affect neurite outgrowth, which in turn can affect neural activity. We carried out the detailed quantitative analysis of spatiotemporal activity patterns and studied the relationship between individual cells and the network as a whole to explore the relationship between developing connectivity and activity patterns. The model, developed in this work will allow us to develop new experimental techniques for studying and quantifying the influence of the neuronal activity on growth processes in neural networks and may lead to a novel techniques for constructing large-scale neural networks by self-organization.

A spatially resolved network spike in model neuronal cultures reveals nucleation centers, circular traveling waves and drifting spiral waves

We show that in model neuronal cultures, where the probability of interneuronal connection formation decreases exponentially with increasing distance between the neurons, there exists a small number of spatial nucleation centers of a network spike, from where the synchronous spiking activity starts propagating in the network typically in the form of circular traveling waves. The number of nucleation centers and their spatial locations are unique and unchanged for a given realization of neuronal network but are different for different networks. In contrast, if the probability of interneuronal connection formation is independent of the distance between neurons, then the nucleation centers do not arise and the synchronization of spiking activity during a network spike occurs spatially uniform throughout the network. Therefore one can conclude that spatial proximity of connections between neurons is important for the formation of nucleation centers. It is also shown that fluctuations of the spatial density of neurons at their random homogeneous distribution typical for the experiments in vitro do not determine the locations of the nucleation centers. The simulation results are qualitatively consistent with the experimental observations.

Biomaterials and computation: a strategic alliance to investigate emergent responses of neural cells

Synergistic use of biomaterials and computation allows to identify and unravel neural cell responses.

Multi-phasic bi-directional chemotactic responses of the growth cone

The nerve growth cone is bi-directionally attracted and repelled by the same cue molecules depending on the situations, while other non-neural chemotactic cells usually show uni-directional attraction or repulsion toward their specific cue molecules. However, how the growth cone differs from other non-neural cells remains unclear. Toward this question, we developed a theory for describing chemotactic response based on a mathematical model of intracellular signaling of activator and inhibitor. Our theory was first able to clarify the conditions of attraction and repulsion, which are determined by balance between activator and inhibitor, and the conditions of uni- and bi-directional responses, which are determined by dose-response profiles of activator and inhibitor to the guidance cue. With biologically realistic sigmoidal dose-responses, our model predicted tri-phasic turning response depending on intracellular Ca²⁺ level, which was then experimentally confirmed by growth cone turning assays and Ca²⁺ imaging. Furthermore, we took a reverse-engineering analysis to identify balanced regulation between CaMKII (activator) and PP1 (inhibitor) and then the model performance was validated by reproducing turning assays with inhibitions of CaMKII and PP1. Thus, our study implies that the balance between activator and inhibitor underlies the multi-phasic bi-directional turning response of the growth cone.

A Model for Improving the Learning Curves of Artificial Neural Networks

In this article, the performance of a hybrid artificial neural network (i.e. scale-free and small-world) was analyzed and its learning curve compared to three other topologies: random, scale-free and small-world, as well as to the chemotaxis neural network of the nematode Caenorhabditis Elegans. One hundred equivalent networks (same number of vertices and average degree) for each topology were generated and each was trained for one thousand epochs. After comparing the mean learning curves of each network topology with the C. elegans neural network, we found that the networks that exhibited preferential attachment exhibited the best learning curves.

Topographical strategies to control neural outgrowth

In this work a synergistic approach is used to investigate how directional anisotropic surfaces (i.e., nanogratings) control the alignment of PC12 neurites. Finite Element models were used to assess the distribution of stresses in non-spread growth cones and filopodia. The stress field was assumed to be the main triggering cause fostering the increase and stabilization of filopodia, so the local stress maxima were directly related to the neuritic orientation. Moreover, a computational framework was implemented within an open source Java environment (CX3D), and in silico simulations were carried out to reproduce and predict biological experiments. No significant differences were found between biological experiments and in silico simulations (alignment angle, p = 0.4685; tortuosity, p = 0.9075) with a standard level of confidence (95%).

Deterministic control of mean alignment and elongation of neuron-like cells by grating geometry: a computational approach

Neuron-like cells are driven by their surrounding environment through local topography. A causal mechanotransductive web of topography-force relationships influences and controls complex cellular phenomena such as growth and alignment. This work aimed to provide a computational framework able to model the behaviour of neuron-like (PC12) cells on gratings, accounting for the twofold ability of topographical cues to simultaneously align and enhance the growth of cells. In particular, starting from the mechanical behaviour of the growth cone and filopodia, the effect of grating geometry (e.g., the periodicity and the size of grooves and ridges) on the neuritic mean alignment angle and on the outgrowth rate of cells was explored through theoretical tools and combinatorial simulations, which were able to predict (R(2) > 0.9) experimental data in a time range of 72-120 hours.

A hybrid computational model to predict chemotactic guidance of growth cones

The overall strategy used by growing axons to find their correct paths during the nervous system development is not yet completely understood. Indeed, some emergent and counterintuitive phenomena were recently described during axon pathfinding in presence of chemical gradients. Here, a novel computational model is presented together with its ability to reproduce both regular and counterintuitive axonal behaviours. In this model, the key role of intracellular calcium was phenomenologically modelled through a non standard Gierer-Meinhardt system, as a crucial factor influencing the growth cone behaviour both in regular and complex conditions. This model was able to explicitly reproduce neuritic paths accounting for the complex interplay between extracellular and intracellular environments, through the sensing capability of the growth cone. The reliability of this approach was proven by using quantitative metrics, numerically supporting the similarity between in silico and biological results in regular conditions (control and attraction). Finally, the model was able to qualitatively predict emergent and counterintuitive phenomena resulting from complex boundary conditions.

Discrete stochastic model for the generation of axonal trees

In this work we propose a 2D discrete stochastic model for the simulation of axonal biogenesis. The model is defined by a third order Markov Chain. The model considers two main processes: the growth process that models the elongation and shape of the neurites and the bifurcation process that models the generation of branches. The growth process depends, among other variables, on the external attraction field generated by a chemoattractant molecule secreted by the target area. We propose an estimation scheme of the involved parameters from real fluorescent confocal microscopy images of single neurons within intact adult Drosophila fly brains. Both normal neurons and neurons in which certain genes were inactivated have been considered (two mutations). In total, 53 images (18 normal, 21 type 1 mutant and 14 type 2 mutant) were used. The model parameters allow us to describe pathological characteristics of the mutated populations.

Active Brownian agents with concentration-dependent chemotactic sensitivity

We study a biologically motivated model of overdamped, autochemotactic Brownian agents with concentration-dependent chemotactic sensitivity. The agents in our model move stochastically and produce a chemical ligand at their current position. The ligand concentration obeys a reaction-diffusion equation and acts as a chemoattractant for the agents, which bias their motion towards higher concentrations of the dynamically altered chemical field. We explore the impact of concentration-dependent response to chemoattractant gradients on large-scale pattern formation, by deriving a coarse-grained macroscopic description of the individual-based model, and compare the conditions for emergence of inhomogeneous solutions for different variants of the chemotactic sensitivity. We focus primarily on the so-called receptor-law sensitivity, which models a nonlinear decrease of chemotactic sensitivity with increasing ligand concentration. Our results reveal qualitative differences between the receptor law, the constant chemotactic response, and the so-called log law, with respect to stability of the homogeneous solution, as well as the emergence of different patterns (labyrinthine structures, clusters, and bubbles) via spinodal decomposition or nucleation. We discuss two limiting cases, where the model can be reduced to the dynamics of single species: (I) the agent density governed by a density-dependent effective diffusion coefficient and (II) the ligand field with an effective bistable, time-dependent reaction rate. In the end, we turn to single clusters of agents, studying domain growth and determining mean characteristics of the stationary inhomogeneous state. Analytical results are confirmed and extended by large-scale GPU simulations of the individual based model.

Neuronal growth as diffusion in an effective potential

Current understanding of neuronal growth is mostly qualitative, as the staggering number of physical and chemical guidance cues involved prohibit a fully quantitative description of axonal dynamics. We report on a general approach that describes axonal growth in vitro, on poly-D-lysine-coated glass substrates, as diffusion in an effective external potential, representing the collective contribution of all causal influences on the growth cone. We use this approach to obtain effective growth rules that reveal an emergent regulatory mechanism for axonal pathfinding on these substrates.

Transient epileptiform signaling during neuronal network development: regulation by external stimulation and bimodal GABAergic activity

A predominance of excitatory activity, with protracted appearance of inhibitory activity, accompanies cortical neuronal development. It is unclear whether or not inhibitory neuronal activity is solicited exclusively by excitatory neurons or whether the transient excitatory activity displayed by developing GABAergic neurons contributes to an excitatory threshold that fosters their conversion to inhibitory activity. We addressed this possibility by culturing murine embryonic neurons on multi-electrode arrays. A wave of individual 0.2-0.4 mV signals ("spikes") appeared between approx. 20-30 days in culture, then declined. A transient wave of high amplitude (>0.5 mV) epileptiform activity coincided with the developmental decline in spikes. Bursts (clusters of ≥3 low-amplitude spikes within 0.7s prior to returning to baseline) persisted following this decline. Addition of the GABAergic antagonist bicuculline initially had no effect on signaling, consistent with delayed development of GABAergic synapses. This was followed by a period in which bicuculline inhibited overall signaling, confirming that GABAergic neurons initially display excitatory activity in ex vivo networks. Following the transient developmental wave of epileptiform signaling, bicuculline induced a resurgence of epileptiform signaling, indicating that GABAergic neurons at this point displayed inhibitory activity. The appearance of transition after the developmental and decline of epileptiform activity, rather than immediately after the developmental decline in lower-amplitude spikes, suggests that the initial excitatory activity of GABAergic neurons contributes to their transition into inhibitory neurons, and that inhibitory GABAergic activity is essential for network development. Prior studies indicate that a minority (25%) of neurons in these cultures were GABAergic, suggesting that inhibitory neurons regulate multiple excitatory neurons. A similar robust increase in signaling following cessation of inhibitory activity in an artificial neural network containing 20% inhibitory neurons supported this conclusion. Even a minor perturbation in GABAergic function may therefore foster initiation and/or amplification of seizure activity, as well as perturbations in long-term potentiation.

Growth dynamics explain the development of spatiotemporal burst activity of young cultured neuronal networks in detail

A typical property of isolated cultured neuronal networks of dissociated rat cortical cells is synchronized spiking, called bursting, starting about one week after plating, when the dissociated cells have sufficiently sent out their neurites and formed enough synaptic connections. This paper is the third in a series of three on simulation models of cultured networks. Our two previous studies [26], [27] have shown that random recurrent network activity models generate intra- and inter-bursting patterns similar to experimental data. The networks were noise or pacemaker-driven and had Izhikevich-neuronal elements with only short-term plastic (STP) synapses (so, no long-term potentiation, LTP, or depression, LTD, was included). However, elevated pre-phases (burst leaders) and after-phases of burst main shapes, that usually arise during the development of the network, were not yet simulated in sufficient detail. This lack of detail may be due to the fact that the random models completely missed network topology .and a growth model. Therefore, the present paper adds, for the first time, a growth model to the activity model, to give the network a time dependent topology and to explain burst shapes in more detail. Again, without LTP or LTD mechanisms. The integrated growth-activity model yielded realistic bursting patterns. The automatic adjustment of various mutually interdependent network parameters is one of the major advantages of our current approach. Spatio-temporal bursting activity was validated against experiment. Depending on network size, wave reverberation mechanisms were seen along the network boundaries, which may explain the generation of phases of elevated firing before and after the main phase of the burst shape.In summary, the results show that adding topology and growth explain burst shapes in great detail and suggest that young networks still lack/do not need LTP or LTD mechanisms.

Mathematical modeling of axonal formation. Part I: Geometry

A stochastic model is proposed for the position of the tip of an axon. Parameters in the model are determined from laboratory data. The first step is the reduction of inherent error in the laboratory data, followed by estimating parameters and fitting a mathematical model to this data. Several axonogenesis aspects have been investigated, particularly how positive axon elongation and growth cone kinematics are coupled processes but require very different theoretical descriptions. Preliminary results have been obtained through a series of experiments aimed at isolating the response of axons to controlled gradient exposures to guidance cues and the effects of ethanol and similar substances. We show results based on the following tasks; (A) development of a novel filtering strategy to obtain data sets truly representative of the axon trail formation; (B) creation of a coarse graining method which establishes (C) an optimal parameter estimation technique, and (D) derivation of a mathematical model which is stochastic in nature, parameterized by arc length. The framework and the resulting model allow for the comparison of experimental and theoretical mean square displacement (MSD) of the developing axon. Current results are focused on uncovering the geometric characteristics of the axons and MSD through analytical solutions and numerical simulations parameterized by arc length, thus ignoring the temporal growth processes. Future developments will capture the dynamic growth cone and how it behaves as a function of time. Qualitative and quantitative predictions of the model at specific length scales capture the experimental behavior well.

Modeling of Neuronal Growth In Vitro: Comparison of Simulation Tools NETMORPH and CX3D

We simulate the growth of neuronal networks using the two recently published tools, NETMORPH and CX3D. The goals of the work are (1) to examine and compare the simulation tools, (2) to construct a model of growth of neocortical cultures, and (3) to characterize the changes in network connectivity during growth, using standard graph theoretic methods. Parameters for the neocortical culture are chosen after consulting both the experimental and the computational work presented in the literature. The first (three) weeks in culture are known to be a time of development of extensive dendritic and axonal arbors and establishment of synaptic connections between the neurons. We simulate the growth of networks from day 1 to day 21. It is shown that for the properly selected parameters, the simulators can reproduce the experimentally obtained connectivity. The selected graph theoretic methods can capture the structural changes during growth.

Mathematical modelling and numerical simulation of the morphological development of neurons

Background

The morphological development of neurons is a very complex process involving both genetic and environmental components. Mathematical modelling and numerical simulation are valuable tools in helping us unravel particular aspects of how individual neurons grow their characteristic morphologies and eventually form appropriate networks with each other.

Methods

A variety of mathematical models that consider (1) neurite initiation (2) neurite elongation (3) axon pathfinding, and (4) neurite branching and dendritic shape formation are reviewed. The different mathematical techniques employed are also described.

Results

Some comparison of modelling results with experimental data is made. A critique of different modelling techniques is given, leading to a proposal for a unified modelling environment for models of neuronal development.

Conclusion

A unified mathematical and numerical simulation framework should lead to an expansion of work on models of neuronal development, as has occurred with compartmental models of neuronal electrical activity.

SELF-WIRING IN NEURAL NETS OF POINT-LIKE CORTICAL NEURONS FAILS TO REPRODUCE CYTOARCHITECTURAL DIFFERENCES

We propose a model for description of activity-dependent evolution and self-wiring between binary neurons. Specifically, this model can be used for investigation of growth of neuronal connectivity in the developing neocortex. By using computational simulations with appropriate training pattern sequences, we show that long-term memory can be encoded in neuronal connectivity and that the external stimulations form part of the functioning neocortical circuit. It is proposed that such binary neuron representations of point-like cortical neurons fail to reproduce cytoarchitectural differences of the neocortical organization, which has implications for inadequacies of compartmental models.

Effects of random external background stimulation on network synaptic stability after tetanization: a modeling study

We constructed a simulated spiking neural network model to investigate the effects of random background stimulation on the dynamics of network activity patterns and tetanus induced network plasticity. The simulated model was a "leaky integrate-and-fire" (LIF) neural model with spike-timing-dependent plasticity (STDP) and frequency-dependent synaptic depression. Spontaneous and evoked activity patterns were compared with those of living neuronal networks cultured on multi-electrode arrays. To help visualize activity patterns and plasticity in our simulated model, we introduced new population measures called Center of Activity (CA) and Center of Weights (CW) to describe the spatio-temporal dynamics of network-wide firing activity and network-wide synaptic strength, respectively. Without random background stimulation, the network synaptic weights were unstable and often drifted after tetanization. In contrast, with random background stimulation, the network synaptic weights remained close to their values immediately after tetanization. The simulation suggests that the effects of tetanization on network synaptic weights were difficult to control because of ongoing synchronized spontaneous bursts of action potentials, or "barrages." Random background stimulation helped maintain network synaptic stability after tetanization by reducing the number and thus the influence of spontaneous barrages. We used our simulated network to model the interaction between ongoing neural activity, external stimulation and plasticity, and to guide our choice of sensory-motor mappings for adaptive behavior in hybrid neural-robotic systems or "hybrots."

Deterministic and Stochastic Elements of Axonal Guidance

An enormous literature has been developed on investigations of the growth and guidance of axons during development and after injury. In this review, we provide a guide to this literature as a resource for biomedical investigators. We first review briefly the molecular biology that is known to regulate migration of the growth cone and branching of axonal arbors. We then outline some important fundamental considerations that are important to the modeling of the phenomenology of these guidance effects and of what is known of their underlying internal mechanisms. We conclude by providing some thoughts on the outlook for future biomedical modeling in the field.

Bacterial self-organization: co-enhancement of complexification and adaptability in a dynamic environment

During colonial development, bacteria generate a wealth of patterns, some of which are reminiscent of those occurring in abiotic systems. They can exhibit rich behaviour, reflecting informative communication capabilities that include exchange of genetic materials and the fact that the colony's building blocks are biotic. Each has internal degrees of freedom, informatic capabilities and freedom to respond by altering itself and others via emission of signals in a self-regulated manner. To unravel the special secrets of bacterial self-organization, we conducted an integrative (experimental and theoretical) study of abiotic and biotic systems. Guided by the notion of general biotic motives and principles, I propose that the informative communication between individuals makes possible the enhancement of the individuals' regulated freedom, while increasing their cooperation. This process is accomplished via cooperative complexification of the colony through self-organization of hierarchical spatio-temporal patterning. The colonial higher complexity provides the degree of plasticity and flexibility required for better colonial adaptability and endurability in a dynamic environment. The biotic system can modify the environment and obtain environmental information for further self-improvement. I reflect on the potential applications of the new understanding on 'engineered self-organization of systems too complex to design' and other issues.

Formation of electrically active clusterized neural networks

Ordinarily, in vitro neurons self-organize into homogeneous networks of single neurons linked by dendrites and axons. We show that under special conditions they can also self-organize into neuronal clusters, which are linked by bundles of axons. Multielectrode array measurement reveals that the clusterized networks are also electrically active and exhibit synchronized bursting events similar to those observed in the homogeneous networks. From time-lapse recording, we deduced the features required for the neuronal clusterized versus homogeneous self-organization and developed a simple model for testing their validity.

Self-assembling of networks in an agent-based model

We propose a model to show the self-assembling of networklike structures between a set of nodes without using preexisting positional information or long-range attraction of the nodes. The model is based on Brownian agents that are capable of producing different local (chemical) information and respond to it in a nonlinear manner. They solve two tasks in parallel: (i) the detection of the appropriate nodes, and (ii) the establishment of stable links between them. We present results of computer simulations that demonstrate the emergence of robust network structures and investigate the connectivity of the network by means of both analytical estimations and computer simulations.

Morphological characterization of in vitro neuronal networks

We use in vitro neuronal networks as a model system for studying self-organization processes in the nervous system. We follow the neuronal growth process, from isolated neurons to fully connected two-dimensional networks. The mature networks are mapped into connected graphs and their morphological characteristics are measured. The distributions of segment lengths, node connectivity, and path length between nodes, and the clustering coefficient of the networks are used to characterize network morphology and to demonstrate that our networks fall into the category of small-world networks.