Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and 120∼330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks.
more »
« less
Deep learning of material transport in complex neurite networks
Abstract Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and $$120 \sim 330$$ 120 ∼ 330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks.
more »
« less
- Award ID(s):
- 1804929
- NSF-PAR ID:
- 10293562
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Neurons exhibit remarkably complex geometry in their neurite networks. So far, how materials are transported in the complex geometry for survival and function of neurons remains an unanswered question. Answering this question is fundamental to understanding the physiology and disease of neurons. Here, we have developed an isogeometric analysis (IGA) based platform for material transport simulation in neurite networks. We modeled the transport process by reaction-diffusion-transport equations and represented geometry of the networks using truncated hierarchical tricubic B-splines (THB-spline3D). We solved the Navier-Stokes equations to obtain the velocity field of material transport in the networks. We then solved the transport equations using the streamline upwind/Petrov-Galerkin (SU/PG) method. Using our IGA solver, we simulated material transport in three basic models of the network geometry: a single neurite, a neurite bifurcation, and a neurite tree with three bifurcations. In addition, the robustness of our solver is illustrated by simulating material transport in three representative and complex neurite networks. From the simulation we discovered several spatial patterns of the transport process. Together, our simulation provides key insights into how material transport in neurite networks is mediated by their complex geometry.more » « less
-
Modeling material transport regulation and traffic jam in neurons using PDE-constrained optimizationAbstract The intracellular transport process plays an important role in delivering essential materials throughout branched geometries of neurons for their survival and function. Many neurodegenerative diseases have been associated with the disruption of transport. Therefore, it is essential to study how neurons control the transport process to localize materials to necessary locations. Here, we develop a novel optimization model to simulate the traffic regulation mechanism of material transport in complex geometries of neurons. The transport is controlled to avoid traffic jam of materials by minimizing a pre-defined objective function. The optimization subjects to a set of partial differential equation (PDE) constraints that describe the material transport process based on a macroscopic molecular-motor-assisted transport model of intracellular particles. The proposed PDE-constrained optimization model is solved in complex tree structures by using isogeometric analysis (IGA). Different simulation parameters are used to introduce traffic jams and study how neurons handle the transport issue. Specifically, we successfully model and explain the traffic jam caused by reduced number of microtubules (MTs) and MT swirls. In summary, our model effectively simulates the material transport process in healthy neurons and also explains the formation of a traffic jam in abnormal neurons. Our results demonstrate that both geometry and MT structure play important roles in achieving an optimal transport process in neuron.more » « less
-
more » « less
Abstract The intracellular transport process plays an important role in delivering essential materials throughout branched geometries of neurons for their survival and function. Many neurodegenerative diseases have been associated with the disruption of transport. Therefore, it is essential to study how neurons control the transport process to localize materials to necessary locations. Here, we develop a novel optimization model to simulate the traffic regulation mechanism of material transport in three-dimensional complex geometries of neurons. The transport is controlled to avoid traffic jams of materials by minimizing a predefined objective function. The optimization subjects to a set of partial differential equation (PDE) constraints that describe the material transport process based on a macroscopic molecular-motor-assisted transport model of intracellular particles. The proposed PDE-constrained optimization model is solved in complex tree structures by using the isogeometric analysis. Different simulation parameters are used to introduce traffic jams and study how neurons handle the transport issue. Specifically, we successfully model and explain the traffic jam caused by the reduced number of microtubules (MTs) and MT swirls. In summary, our model effectively simulates the material transport process in healthy neurons and also explains the formation of a traffic jam in abnormal neurons. Our results demonstrate that both geometry and MT structure play important roles in achieving an optimal transport process in neurons.
-
Biological regulatory networks depend upon chemical interactions to process information. Engineering such molecular computing systems is a major challenge for synthetic biology and related fields. The chemical reaction network (CRN) model idealizes chemical interactions, abstracting away specifics of the molecular implementation, and allowing rigorous reasoning about the computational power of chemical kinetics. Here we focus on function computation with CRNs, where we think of the initial concentrations of some species as the input and the eventual steady-state concentration of another species as the output. Specifically, we are concerned with CRNs that are rate-independent (the computation must be correct independent of the reaction rate law) and composable (𝑓∘𝑔 can be computed by concatenating the CRNs computing f and g). Rate independence and composability are important engineering desiderata, permitting implementations that violate mass-action kinetics, or even “well-mixedness”, and allowing the systematic construction of complex computation via modular design. We show that to construct composable rate-independent CRNs, it is necessary and sufficient to ensure that the output species of a module is not a reactant in any reaction within the module. We then exactly characterize the functions computable by such CRNs as superadditive, positive-continuous, and piecewise rational linear. Our results show that composability severely limits rate-independent computation unless more sophisticated input/output encodings are used.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41598-021-90724-3
