1
|
Staii C. Biased Random Walk Model of Neuronal Dynamics on Substrates with Periodic Geometrical Patterns. Biomimetics (Basel) 2023; 8:267. [PMID: 37366862 DOI: 10.3390/biomimetics8020267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Revised: 06/07/2023] [Accepted: 06/16/2023] [Indexed: 06/28/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Cristian Staii
- Department of Physics and Astronomy, Tufts University, Medford, MA 02155, USA
| |
Collapse
|
2
|
Abstract
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.
Collapse
Affiliation(s)
- Hadrien Oliveri
- Mathematical Institute, University of Oxford, Oxford, OX2 6GG, UK
| | - Alain Goriely
- Mathematical Institute, University of Oxford, Oxford, OX2 6GG, UK.
| |
Collapse
|
3
|
Hartley C, Farmer S, Berthouze L. Temporal ordering of input modulates connectivity formation in a developmental neuronal network model of the cortex. PLoS One 2020; 15:e0226772. [PMID: 31923200 PMCID: PMC6953763 DOI: 10.1371/journal.pone.0226772] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2019] [Accepted: 12/06/2019] [Indexed: 12/13/2022] Open
Abstract
Preterm infant brain activity is discontinuous; bursts of activity recorded using EEG (electroencephalography), thought to be driven by subcortical regions, display scale free properties and exhibit a complex temporal ordering known as long-range temporal correlations (LRTCs). During brain development, activity-dependent mechanisms are essential for synaptic connectivity formation, and abolishing burst activity in animal models leads to weak disorganised synaptic connectivity. Moreover, synaptic pruning shares similar mechanisms to spike-timing dependent plasticity (STDP), suggesting that the timing of activity may play a critical role in connectivity formation. We investigated, in a computational model of leaky integrate-and-fire neurones, whether the temporal ordering of burst activity within an external driving input could modulate connectivity formation in the network. Connectivity evolved across the course of simulations using an approach analogous to STDP, from networks with initial random connectivity. Small-world connectivity and hub neurones emerged in the network structure—characteristic properties of mature brain networks. Notably, driving the network with an external input which exhibited LRTCs in the temporal ordering of burst activity facilitated the emergence of these network properties, increasing the speed with which they emerged compared with when the network was driven by the same input with the bursts randomly ordered in time. Moreover, the emergence of small-world properties was dependent on the strength of the LRTCs. These results suggest that the temporal ordering of burst activity could play an important role in synaptic connectivity formation and the emergence of small-world topology in the developing brain.
Collapse
Affiliation(s)
- Caroline Hartley
- Centre for Mathematics and Physics in the Life Sciences and Experimental Biology, University College London, London, United Kingdom
- Department of Paediatrics, University of Oxford, Oxford, United Kingdom
- * E-mail:
| | - Simon Farmer
- Institute of Neurology, University College London, London, United Kingdom
| | - Luc Berthouze
- Centre for Computational Neuroscience and Robotics, University of Sussex, Falmer, Brighton, United Kingdom
| |
Collapse
|
4
|
Marinov T, López Sánchez HA, Yuchi L, Adewole DO, Cullen DK, Kraft RH. A computational model of bidirectional axonal growth in micro-tissue engineered neuronal networks (micro-TENNs). In Silico Biol 2020; 14:85-99. [PMID: 32390612 PMCID: PMC7505002 DOI: 10.3233/isb-180172] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Micro-Tissue Engineered Neural Networks (Micro-TENNs) are living three-dimensional constructs designed to replicate the neuroanatomy of white matter pathways in the brain and are being developed as implantable micro-tissue for axon tract reconstruction, or as anatomically-relevant in vitro experimental platforms. Micro-TENNs are composed of discrete neuronal aggregates connected by bundles of long-projecting axonal tracts within miniature tubular hydrogels. In order to help design and optimize micro-TENN performance, we have created a new computational model including geometric and functional properties. The model is built upon the three-dimensional diffusion equation and incorporates large-scale uni- and bi-directional growth that simulates realistic neuron morphologies. The model captures unique features of 3D axonal tract development that are not apparent in planar outgrowth and may be insightful for how white matter pathways form during brain development. The processes of axonal outgrowth, branching, turning and aggregation/bundling from each neuron are described through functions built on concentration equations and growth time distributed across the growth segments. Once developed we conducted multiple parametric studies to explore the applicability of the method and conducted preliminary validation via comparisons to experimentally grown micro-TENNs for a range of growth conditions. Using this framework, the model can be applied to study micro-TENN growth processes and functional characteristics using spiking network or compartmental network modeling. This model may be applied to improve our understanding of axonal tract development and functionality, as well as to optimize the fabrication of implantable tissue engineered brain pathways for nervous system reconstruction and/or modulation.
Collapse
Affiliation(s)
- Toma Marinov
- Penn State Computational Biomechanics Group, Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, PA, USA
| | - Haven A. López Sánchez
- The Laboratory of Physicochemistry and Engineering of Proteins, Department of Biochemistry, Facultad de Medicina, National Autonomous University of Mexico, Mexico
| | - Liang Yuchi
- Penn State Computational Biomechanics Group, Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, PA, USA
| | - Dayo O. Adewole
- Center for Brain Injury & Repair, Department of Neurosurgery, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Department of Bioengineering, School of Engineering & Applied Science, University of Pennsylvania, Philadelphia, PA, USA
- Center for Neurotrauma, Neurodegeneration & Restoration, Corporal Michael J. Crescenz Veterans Affairs Medical Center, Philadelphia, PA, USA
| | - D. Kacy Cullen
- Center for Brain Injury & Repair, Department of Neurosurgery, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Department of Bioengineering, School of Engineering & Applied Science, University of Pennsylvania, Philadelphia, PA, USA
- Center for Neurotrauma, Neurodegeneration & Restoration, Corporal Michael J. Crescenz Veterans Affairs Medical Center, Philadelphia, PA, USA
| | - Reuben H. Kraft
- Penn State Computational Biomechanics Group, Mechanical and Nuclear Engineering, Pennsylvania State University, University Park, PA, USA
- Department of Biomedical Engineering, Pennsylvania State University, University Park, PA, USA
| |
Collapse
|
5
|
Nguyen H, Dayan P, Pujic Z, Cooper-White J, Goodhill GJ. A mathematical model explains saturating axon guidance responses to molecular gradients. eLife 2016; 5:e12248. [PMID: 26830461 PMCID: PMC4755759 DOI: 10.7554/elife.12248] [Citation(s) in RCA: 89] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2015] [Accepted: 12/18/2015] [Indexed: 11/13/2022] Open
Abstract
Correct wiring is crucial for the proper functioning of the nervous system. Molecular gradients provide critical signals to guide growth cones, which are the motile tips of developing axons, to their targets. However, in vitro, growth cones trace highly stochastic trajectories, and exactly how molecular gradients bias their movement is unclear. Here, we introduce a mathematical model based on persistence, bias, and noise to describe this behaviour, constrained directly by measurements of the detailed statistics of growth cone movements in both attractive and repulsive gradients in a microfluidic device. This model provides a mathematical explanation for why average axon turning angles in gradients in vitro saturate very rapidly with time at relatively small values. This work introduces the most accurate predictive model of growth cone trajectories to date, and deepens our understanding of axon guidance events both in vitro and in vivo.
Collapse
Affiliation(s)
- Huyen Nguyen
- Queensland Brain Institute, The University of Queensland, St. Lucia, Australia.,School of Mathematics and Physics, The University of Queensland, St. Lucia, Australia
| | - Peter Dayan
- Gatsby Computational Neuroscience Unit, University College London, London, United Kingdom
| | - Zac Pujic
- Queensland Brain Institute, The University of Queensland, St. Lucia, Australia
| | - Justin Cooper-White
- Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, St. Lucia, Australia
| | - Geoffrey J Goodhill
- Queensland Brain Institute, The University of Queensland, St. Lucia, Australia.,School of Mathematics and Physics, The University of Queensland, St. Lucia, Australia
| |
Collapse
|
6
|
Roccasalvo IM, Sergi PN, Tonazzini I, Cecchini M, Micera S. Topographical strategies to control neural outgrowth. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2016; 2015:7147-50. [PMID: 26737940 DOI: 10.1109/embc.2015.7320040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
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%).
Collapse
|
7
|
Roccasalvo IM, Micera S, Sergi PN. A hybrid computational model to predict chemotactic guidance of growth cones. Sci Rep 2015; 5:11340. [PMID: 26086936 PMCID: PMC4471899 DOI: 10.1038/srep11340] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2014] [Accepted: 05/15/2015] [Indexed: 11/08/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Silvestro Micera
- The BioRobotics Institute, Scuola Superiore Sant’Anna, Pisa, Italy
- Bertarelli Foundation Chair in Translational NeuroEngineering Laboratory, Institute of Bioengineering, School of Engineering, Ecole Polytechnique Federale de Lausanne, Lausanne, Switzerland
- Center for Neuroprosthetics, Ecole Polytechnique Federale de Lausanne, Lausanne, Switzerland
| | | |
Collapse
|
8
|
Sergi PN, Morana Roccasalvo I, Tonazzini I, Cecchini M, Micera S. Cell guidance on nanogratings: a computational model of the interplay between PC12 growth cones and nanostructures. PLoS One 2013; 8:e70304. [PMID: 23936404 PMCID: PMC3735603 DOI: 10.1371/journal.pone.0070304] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2013] [Accepted: 06/17/2013] [Indexed: 12/05/2022] Open
Abstract
Background Recently, the effects of nanogratings have been investigated on PC12 with respect to cell polarity, neuronal differentiation, migration, maturation of focal adhesions and alignment of neurites. Methodology/Principal Findings A synergistic procedure was used to study the mechanism of alignment of PC12 neurites with respect to the main direction of nanogratings. Finite Element simulations were used to qualitatively assess the distribution of stresses at the interface between non-spread growth cones and filopodia, and to study their dependence on filopodial length and orientation. After modelling all adhesions under non-spread growth cone and filopodial protrusions, the values of local stress maxima resulted from the length of filopodia. Since the stress was assumed to be the main triggering cause leading to the increase and stabilization of filopodia, the position of the local maxima was directly related to the orientation of neurites. An analytic closed form equation was then written to quantitatively assess the average ridge width needed to achieve a given neuritic alignment (R2 = 0.96), and the alignment course, when the ridge depth varied (R2 = 0.97). A computational framework was implemented within an improved free 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, p = 0.3571; tortuosity, p = 0.2236) with a standard level of confidence (95%). Conclusions/Significance A mechanism involved in filopodial sensing of nanogratings is proposed and modelled through a synergistic use of FE models, theoretical equations and in silico simulations. This approach shows the importance of the neuritic terminal geometry, and the key role of the distribution of the adhesion constraints for the cell/substrate coupling process. Finally, the effects of the geometry of nanogratings were explicitly considered in cell/surface interactions thanks to the analytic framework presented in this work.
Collapse
Affiliation(s)
- Pier Nicola Sergi
- Neural Engineering Area, The BioRobotics Institute, Scuola Superiore Sant'Anna, Pisa, Italy.
| | | | | | | | | |
Collapse
|
9
|
Tracking Growing Axons by Particle Filtering in 3D + t Fluorescent Two-Photon Microscopy Images. ACTA ACUST UNITED AC 2013. [DOI: 10.1007/978-3-642-37431-9_21] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2023]
|
10
|
van Ooyen A. Using theoretical models to analyse neural development. Nat Rev Neurosci 2011; 12:311-26. [DOI: 10.1038/nrn3031] [Citation(s) in RCA: 95] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
11
|
Pearson YE, Castronovo E, Lindsley TA, Drew DA. Mathematical modeling of axonal formation. Part I: Geometry. Bull Math Biol 2011; 73:2837-64. [PMID: 21390561 DOI: 10.1007/s11538-011-9648-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2010] [Accepted: 02/18/2011] [Indexed: 12/28/2022]
Abstract
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.
Collapse
|
12
|
Staii C, Viesselmann C, Ballweg J, Williams JC, Dent EW, Coppersmith SN, Eriksson MA. Distance dependence of neuronal growth on nanopatterned gold surfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2011; 27:233-9. [PMID: 21121598 DOI: 10.1021/la102331x] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Understanding network development in the brain is of tremendous fundamental importance, but it is immensely challenging because of the complexity of both its architecture and function. The mechanisms of axonal navigation to target regions and the specific interactions with guidance factors such as membrane-bound proteins, chemical gradients, mechanical guidance cues, etc., are largely unknown. A current limitation for the study of neural network formation is the ability to control precisely the connectivity of small groups of neurons. A first step in designing such networks is to understand the "rules" central nervous system (CNS) neurons use to form functional connections with one another. Here we begin to delineate novel rules for growth and connectivity of small numbers of neurons patterned on Au substrates in simplified geometries. These studies yield new insights into the mechanisms determining the organizational features present in intact systems. We use a previously reported atomic force microscopy (AFM) nanolithography method to control precisely the location and growth of neurons on these surfaces. By examining a series of systems with different geometrical parameters, we quantitatively and systematically analyze how neuronal growth depends on these parameters.
Collapse
Affiliation(s)
- Cristian Staii
- Department of Physics and Astronomy, Tufts University, Medford, Massachusetts 02155, United States.
| | | | | | | | | | | | | |
Collapse
|
13
|
A hybrid approach for the control of axonal outgrowth: preliminary simulation results. Med Biol Eng Comput 2010; 49:163-70. [DOI: 10.1007/s11517-010-0687-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2009] [Accepted: 09/15/2010] [Indexed: 12/17/2022]
|
14
|
Kobayashi T, Terajima K, Nozumi M, Igarashi M, Akazawa K. A stochastic model of neuronal growth cone guidance regulated by multiple sensors. J Theor Biol 2010; 266:712-22. [PMID: 20688081 DOI: 10.1016/j.jtbi.2010.07.036] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2009] [Revised: 07/25/2010] [Accepted: 07/28/2010] [Indexed: 01/24/2023]
Abstract
Neuronal growth cones migrate directionally under the control of axon guidance molecules, thereby forming synapses in the developing brain. The signal transduction system by which a growth cone detects surrounding guidance molecules, analyzes the detected signals, and then determines the overall behavior remains undetermined. In this study, we describe a novel stochastic model of this behavior that utilizes multiple sensors on filopodia to respond to guidance molecules. Overall growth cone behavior is determined by using only the concentration gradients of guidance molecules in the immediate vicinity of each sensor. The detected signal at each sensor, which is treated as a vector quantity, is sent to the growth cone center and then integrated to determine axonal growth in the next step by means of a simple vector operation. We compared the results of computer simulations of axonal growth with observations of actual axonal growth from co-culture experiments using olfactory bulb and septum. The probabilistic distributions of axonal growth generated by the computer simulation were consistent with those obtained from the culture experiments, indicating that our model accurately simulates growth cone behavior. We believe that this model will be useful for elucidating the as yet unknown mechanisms responsible for axonal growth in vivo.
Collapse
Affiliation(s)
- Taichiro Kobayashi
- Division of Information Science and Biostatistics, Niigata University Graduate School of Medical and Dental Sciences, 1 Asahimachi, Niigata 951-8520, Japan
| | | | | | | | | |
Collapse
|
15
|
Haines C, Goodhill GJ. Analyzing neurite outgrowth from explants by fitting ellipses. J Neurosci Methods 2010; 187:52-8. [DOI: 10.1016/j.jneumeth.2009.12.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2009] [Revised: 12/16/2009] [Accepted: 12/16/2009] [Indexed: 11/29/2022]
|
16
|
Zubler F, Douglas R. A framework for modeling the growth and development of neurons and networks. Front Comput Neurosci 2009; 3:25. [PMID: 19949465 PMCID: PMC2784082 DOI: 10.3389/neuro.10.025.2009] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2009] [Accepted: 10/19/2009] [Indexed: 01/21/2023] Open
Abstract
The development of neural tissue is a complex organizing process, in which it is difficult to grasp how the various localized interactions between dividing cells leads relentlessly to global network organization. Simulation is a useful tool for exploring such complex processes because it permits rigorous analysis of observed global behavior in terms of the mechanistic axioms declared in the simulated model. We describe a novel simulation tool, CX3D, for modeling the development of large realistic neural networks such as the neocortex, in a physical 3D space. In CX3D, as in biology, neurons arise by the replication and migration of precursors, which mature into cells able to extend axons and dendrites. Individual neurons are discretized into spherical (for the soma) and cylindrical (for neurites) elements that have appropriate mechanical properties. The growth functions of each neuron are encapsulated in set of pre-defined modules that are automatically distributed across its segments during growth. The extracellular space is also discretized, and allows for the diffusion of extracellular signaling molecules, as well as the physical interactions of the many developing neurons. We demonstrate the utility of CX3D by simulating three interesting developmental processes: neocortical lamination based on mechanical properties of tissues; a growth model of a neocortical pyramidal cell based on layer-specific guidance cues; and the formation of a neural network in vitro by employing neurite fasciculation. We also provide some examples in which previous models from the literature are re-implemented in CX3D. Our results suggest that CX3D is a powerful tool for understanding neural development.
Collapse
Affiliation(s)
- Frederic Zubler
- Institute of Neuroinformatics, University of Zurich/Swiss Federal Institute of Technology Zurich Zurich, Switzerland.
| | | |
Collapse
|
17
|
|
18
|
Causin P, Facchetti G. Autocatalytic loop, amplification and diffusion: a mathematical and computational model of cell polarization in neural chemotaxis. PLoS Comput Biol 2009; 5:e1000479. [PMID: 19714204 PMCID: PMC2722090 DOI: 10.1371/journal.pcbi.1000479] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2008] [Accepted: 07/21/2009] [Indexed: 12/11/2022] Open
Abstract
The chemotactic response of cells to graded fields of chemical cues is a complex process that requires the coordination of several intracellular activities. Fundamental steps to obtain a front vs. back differentiation in the cell are the localized distribution of internal molecules and the amplification of the external signal. The goal of this work is to develop a mathematical and computational model for the quantitative study of such phenomena in the context of axon chemotactic pathfinding in neural development. In order to perform turning decisions, axons develop front-back polarization in their distal structure, the growth cone. Starting from the recent experimental findings of the biased redistribution of receptors on the growth cone membrane, driven by the interaction with the cytoskeleton, we propose a model to investigate the significance of this process. Our main contribution is to quantitatively demonstrate that the autocatalytic loop involving receptors, cytoplasmic species and cytoskeleton is adequate to give rise to the chemotactic behavior of neural cells. We assess the fact that spatial bias in receptors is a precursory key event for chemotactic response, establishing the necessity of a tight link between upstream gradient sensing and downstream cytoskeleton dynamics. We analyze further crosslinked effects and, among others, the contribution to polarization of internal enzymatic reactions, which entail the production of molecules with a one-to-more factor. The model shows that the enzymatic efficiency of such reactions must overcome a threshold in order to give rise to a sufficient amplification, another fundamental precursory step for obtaining polarization. Eventually, we address the characteristic behavior of the attraction/repulsion of axons subjected to the same cue, providing a quantitative indicator of the parameters which more critically determine this nontrivial chemotactic response.
Collapse
Affiliation(s)
- Paola Causin
- Department of Mathematics F Enriques, Università degli Studi di Milano, Milano, Italy.
| | | |
Collapse
|
19
|
Abstract
Proper wiring up of the nervous system is critical to the development of organisms capable of complex and adaptable behaviors. Besides the many experimental advances in determining the cellular and molecular machinery that carries out this remarkable task precisely and robustly, theoretical approaches have also proven to be useful tools in analyzing this machinery. A quantitative understanding of these processes can allow us to make predictions, test hypotheses, and appraise established concepts in a new light. Three areas that have been fruitful in this regard are axon guidance, retinotectal mapping, and activity-dependent development. This chapter reviews some of the contributions made by mathematical modeling in these areas, illustrated by important examples of models in each section. For axon guidance, we discuss models of how growth cones respond to their environment, and how this environment can place constraints on growth cone behavior. Retinotectal mapping looks at computational models for how topography can be generated in populations of neurons based on molecular gradients and other mechanisms such as competition. In activity-dependent development, we discuss theoretical approaches largely based on Hebbian synaptic plasticity rules, and how they can generate maps in the visual cortex very similar to those seen in vivo. We show how theoretical approaches have substantially contributed to the advancement of developmental neuroscience, and discuss future directions for mathematical modeling in the field.
Collapse
|
20
|
Aletti G, Causin P. Mathematical characterisation of the transduction chain in growth cone pathfinding. IET Syst Biol 2008; 2:150-61. [PMID: 18537455 DOI: 10.1049/iet-syb:20070059] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Axon guidance by graded diffusible ligands plays an important role in the developing nervous system. Concentration gradients induce an asymmetric localisation of molecules in the axon tip, the growth cone, and the consequent internal polarised signalling pathway leads to rearrangement of the growth cone cytoskeleton and, ultimately, to motility. Here the authors provide a mathematical description of the growth cone transduction chain as a series of functional boxes characterised by input/output relations. The model relies on the assumption that the characteristic time of independent concentration measures by growth cone receptors, the characteristic time of growth cone internal reorganisation preceding motion and the characteristic time needed for a discernible axon turning belong to separated scales. The results give insight into the deterministic against stochastic regime of internal growth cone functions that are not readily accessible from experimental observations, pointing out a substantial equilibrium of the two contributions. The mathematical model predicts the decrease of the coefficient of variation of the signal moving down the functional chain leading to motion. Moreover, possible mechanisms that allow for buffering against noise are highlighted. These results have an interest also for the more experimentally minded reader, since they can be used to predict sample sizes for detecting significant differences in benchmark gradient assays.
Collapse
Affiliation(s)
- G Aletti
- Università degli Studi di Milano, Dipartimento di Matematica F. Enriques, Milano, Italy.
| | | |
Collapse
|
21
|
Borisyuk R, Cooke T, Roberts A. Stochasticity and functionality of neural systems: Mathematical modelling of axon growth in the spinal cord of tadpole. Biosystems 2008; 93:101-14. [DOI: 10.1016/j.biosystems.2008.03.012] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2008] [Revised: 03/10/2008] [Accepted: 03/20/2008] [Indexed: 12/24/2022]
|
22
|
Mortimer D, Fothergill T, Pujic Z, Richards LJ, Goodhill GJ. Growth cone chemotaxis. Trends Neurosci 2008; 31:90-8. [DOI: 10.1016/j.tins.2007.11.008] [Citation(s) in RCA: 141] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2007] [Revised: 11/29/2007] [Accepted: 11/29/2007] [Indexed: 12/16/2022]
|
23
|
Abstract
A filopodium is a cytoplasmic projection, exquisitely built and regulated, which extends from the leading edge of the migrating cell, exploring the cell's neighborhood. Commonly, filopodia grow and retract after their initiation, exhibiting rich dynamical behaviors. We model the growth of a filopodium based on a stochastic description which incorporates mechanical, physical, and biochemical components. Our model provides a full stochastic treatment of the actin monomer diffusion and polymerization of each individual actin filament under stress of the fluctuating membrane. We investigated the length distribution of individual filaments in a growing filopodium and studied how it depends on various physical parameters. The distribution of filament lengths turned out to be narrow, which we explained by the negative feedback created by the membrane load and monomeric G-actin gradient. We also discovered that filopodial growth is strongly diminished upon increasing retrograde flow, suggesting that regulating the retrograde flow rate would be a highly efficient way to control filopodial extension dynamics. The filopodial length increases as the membrane fluctuations decrease, which we attributed to the unequal loading of the membrane force among individual filaments, which, in turn, results in larger average polymerization rates. We also observed significant diffusional noise of G-actin monomers, which leads to smaller G-actin flux along the filopodial tube compared with the prediction using the diffusion equation. Overall, partial cancellation of these two fluctuation effects allows a simple mean field model to rationalize most of our simulation results. However, fast fluctuations significantly renormalize the mean field model parameters. The biological significance of our filopodial model and avenues for future development are also discussed.
Collapse
|
24
|
Li WC, Cooke T, Sautois B, Soffe SR, Borisyuk R, Roberts A. Axon and dendrite geography predict the specificity of synaptic connections in a functioning spinal cord network. Neural Dev 2007; 2:17. [PMID: 17845723 PMCID: PMC2071915 DOI: 10.1186/1749-8104-2-17] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2007] [Accepted: 09/10/2007] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND How specific are the synaptic connections formed as neuronal networks develop and can simple rules account for the formation of functioning circuits? These questions are assessed in the spinal circuits controlling swimming in hatchling frog tadpoles. This is possible because detailed information is now available on the identity and synaptic connections of the main types of neuron. RESULTS The probabilities of synapses between 7 types of identified spinal neuron were measured directly by making electrical recordings from 500 pairs of neurons. For the same neuron types, the dorso-ventral distributions of axons and dendrites were measured and then used to calculate the probabilities that axons would encounter particular dendrites and so potentially form synaptic connections. Surprisingly, synapses were found between all types of neuron but contact probabilities could be predicted simply by the anatomical overlap of their axons and dendrites. These results suggested that synapse formation may not require axons to recognise specific, correct dendrites. To test the plausibility of simpler hypotheses, we first made computational models that were able to generate longitudinal axon growth paths and reproduce the axon distribution patterns and synaptic contact probabilities found in the spinal cord. To test if probabilistic rules could produce functioning spinal networks, we then made realistic computational models of spinal cord neurons, giving them established cell-specific properties and connecting them into networks using the contact probabilities we had determined. A majority of these networks produced robust swimming activity. CONCLUSION Simple factors such as morphogen gradients controlling dorso-ventral soma, dendrite and axon positions may sufficiently constrain the synaptic connections made between different types of neuron as the spinal cord first develops and allow functional networks to form. Our analysis implies that detailed cellular recognition between spinal neuron types may not be necessary for the reliable formation of functional networks to generate early behaviour like swimming.
Collapse
Affiliation(s)
- Wen-Chang Li
- School of Biological Sciences, University of Bristol, Woodland Road, Bristol BS8 1UG, UK
| | - Tom Cooke
- Centre for Theoretical and Computational Neuroscience, University of Plymouth, Plymouth PL4 8AA, UK
| | - Bart Sautois
- Department of Applied Mathematics and Computer Science, Ghent University, Krijgslaan 281-S9, B-9000 Ghent, Belgium
| | - Stephen R Soffe
- School of Biological Sciences, University of Bristol, Woodland Road, Bristol BS8 1UG, UK
| | - Roman Borisyuk
- Centre for Theoretical and Computational Neuroscience, University of Plymouth, Plymouth PL4 8AA, UK
| | - Alan Roberts
- School of Biological Sciences, University of Bristol, Woodland Road, Bristol BS8 1UG, UK
| |
Collapse
|
25
|
Abstract
Filopodia are actin-based structures composed of parallel bundles of actin filaments and various actin-associated proteins, and they play important roles in cell-cell signaling, guidance toward chemoattractants, and adhesion to the extracellular matrix. Two mechanisms for the formation of filopodia have been suggested, each using different sets of actin-regulating proteins, creating some controversy in the field. New molecules, some of unknown functions, have also been implicated in filopodium formation, suggesting that other possible mechanisms of filopodium formation exist. We discuss established and novel proteins that mediate the formation and dynamics of filopodia, different mechanisms of filopodium formation, and the various functions that distinct filopodia perform.
Collapse
Affiliation(s)
- Stephanie L Gupton
- Department of Biology, Center for Cancer Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.
| | | |
Collapse
|
26
|
Thivierge JP, Balaban E. Getting into shape: Optimal ligand gradients for axonal guidance. Biosystems 2007; 90:61-77. [PMID: 17187926 DOI: 10.1016/j.biosystems.2006.07.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2005] [Revised: 07/02/2006] [Accepted: 07/02/2006] [Indexed: 12/01/2022]
Abstract
During neural development, neurons from downstream, presynaptic regions of the nervous system (such as the retina) send spatially patterned axonal projections to upstream, target regions (the tectum or superior colliculus). A servomechanism model has been proposed to explain the pattern and time-course of axonal growth between these two regions [Honda, H., 1998. Topographic mapping in the retinotectal projection by means of complementary ligand and receptor gradients: a computer simulation study. J. Theor. Biol., 192, 235-246]. Here, we show that a modification of this model incorporating a different criterion for axonal decision-making, called the local optimum rule, is guaranteed to converge to a topographic map under a wide range of conditions encountered during neural development. A theoretical investigation of these conditions leads to new hypotheses regarding map formation.
Collapse
|
27
|
Gafarov FM. SELF-WIRING IN NEURAL NETS OF POINT-LIKE CORTICAL NEURONS FAILS TO REPRODUCE CYTOARCHITECTURAL DIFFERENCES. J Integr Neurosci 2006; 5:159-69. [PMID: 16783866 DOI: 10.1142/s0219635206001136] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2006] [Accepted: 04/03/2006] [Indexed: 11/18/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Fail M Gafarov
- Department of Theoretical Physics, Tatar State University of Humanity and Pedagogic, 420021 Kazan, Mezhlauk Street, 1, Russia.
| |
Collapse
|
28
|
Carvalho RF, Beutler M, Marler KJM, Knöll B, Becker-Barroso E, Heintzmann R, Ng T, Drescher U. Silencing of EphA3 through a cis interaction with ephrinA5. Nat Neurosci 2006; 9:322-30. [PMID: 16491080 DOI: 10.1038/nn1655] [Citation(s) in RCA: 146] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2005] [Accepted: 01/27/2006] [Indexed: 11/09/2022]
Abstract
EphAs and ephrinAs are expressed in multiple areas of the developing brain in overlapping countergradients, notably in the retina and tectum. Here they are involved in targeting retinal axons to their correct topographic position in the tectum. We have used truncated versions of EphA3, single-amino acid point mutants of ephrinA5 and fluorescence resonance energy transfer technology to uncover a cis interaction between EphA3 and ephrinA5 that is independent of the established ligand-binding domain of EphA3. This cis interaction abolishes the induction of tyrosine phosphorylation of EphA3 and results in a loss of sensitivity of retinal axons to ephrinAs in trans. Our data suggest that formation of this complex transforms the uniform expression of EphAs in the nasal part of the retina into a gradient of functional EphAs and has a key role in controlling retinotectal mapping.
Collapse
Affiliation(s)
- Ricardo F Carvalho
- Medical Research Council, Centre for Developmental Neurobiology, London SE1 1UL, UK
| | | | | | | | | | | | | | | |
Collapse
|
29
|
Xu J, Rosoff WJ, Urbach JS, Goodhill GJ. Adaptation is not required to explain the long-term response of axons to molecular gradients. Development 2005; 132:4545-52. [PMID: 16176951 DOI: 10.1242/dev.02029] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
It has been suggested that growth cones navigating through the developing nervous system might display adaptation, so that their response to gradient signals is conserved over wide variations in ligand concentration. Recently however, a new chemotaxis assay that allows the effect of gradient parameters on axonal trajectories to be finely varied has revealed a decline in gradient sensitivity on either side of an optimal concentration. We show that this behavior can be quantitatively reproduced with a computational model of axonal chemotaxis that does not employ explicit adaptation. Two crucial components of this model required to reproduce the observed sensitivity are spatial and temporal averaging. These can be interpreted as corresponding, respectively, to the spatial spread of signaling effects downstream from receptor binding, and to the finite time over which these signaling effects decay. For spatial averaging, the model predicts that an effective range of roughly one-third of the extent of the growth cone is optimal for detecting small gradient signals. For temporal decay, a timescale of about 3 minutes is required for the model to reproduce the experimentally observed sensitivity.
Collapse
Affiliation(s)
- Jun Xu
- Department of Neuroscience, Georgetown University Medical Center, 3900 Reservoir Road NW, Washington, DC 20007, USA
| | | | | | | |
Collapse
|
30
|
Sakumura Y, Tsukada Y, Yamamoto N, Ishii S. A molecular model for axon guidance based on cross talk between rho GTPases. Biophys J 2005; 89:812-22. [PMID: 15923236 PMCID: PMC1366631 DOI: 10.1529/biophysj.104.055624] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2004] [Accepted: 05/02/2005] [Indexed: 12/11/2022] Open
Abstract
To systematically understand the molecular events that underlie biological phenomena, we must develop methods to integrate an enormous amount of genomic and proteomic data. The integration of molecular data should go beyond the construction of biochemical cascades among molecules to include tying the biochemical phenomena to physical events. For the behavior and guidance of growth cones, it remains largely unclear how biochemical events in the cytoplasm are linked to the morphological changes of the growth cone. We take a computational approach to simulate the biochemical signaling cascade involving members of the Rho family of GTPases and examine their potential roles in growth-cone motility and axon guidance. Based on the interactions between Cdc42, Rac, and RhoA, we show that the activation of a Cdc42-specific GEF resulted in switching responses between oscillatory and convergent activities for all three GTPases. We propose that the switching responses of these GTPases are the molecular basis for the decision mechanism that determines the direction of the growth-cone expansion, providing a spatiotemporal integration mechanism that allows the growth cone to detect small gradients of external guidance cues. These results suggest a potential role for the cross talk between Rho GTPases in governing growth-cone movement and axon guidance and underscore the link between chemodynamic reactions and cellular behaviors.
Collapse
Affiliation(s)
- Yuichi Sakumura
- Graduate School of Information Science, Nara Institute of Science and Technology, Ikoma, Japan
| | | | | | | |
Collapse
|
31
|
Goodhill GJ, Xu J. The development of retinotectal maps: a review of models based on molecular gradients. NETWORK (BRISTOL, ENGLAND) 2005; 16:5-34. [PMID: 16353341 DOI: 10.1080/09548980500254654] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Information about the world is often represented in the brain in the form of topographic maps. A paradigm example is the topographic representation of the visual world in the optic tectum/superior colliculus. This map initially forms during neural development using activity-independent molecular cues, most notably some type of chemospecific matching between molecular gradients in the retina and corresponding gradients in the tectum/superior colliculus. Exactly how this process might work has been studied both experimentally and theoretically for several decades. This review discusses the experimental data briefly, and then in more detail the theoretical models proposed. The principal conclusions are that (1) theoretical models have helped clarify several important ideas in the field, (2) earlier models were often more sophisticated than more recent models, and (3) substantial revisions to current modelling approaches are probably required to account for more than isolated subsets of the experimental data.
Collapse
Affiliation(s)
- Geoffrey J Goodhill
- Queensland Brain Institute, Department of Mathematics, and Institute for Molecular Bioscience, University of Queensland, St. Lucia, Australia.
| | | |
Collapse
|