CN112582069A - Neuron dendritic spine development mode establishing method based on reaction diffusion model - Google Patents
Neuron dendritic spine development mode establishing method based on reaction diffusion model Download PDFInfo
- Publication number
- CN112582069A CN112582069A CN202011573380.XA CN202011573380A CN112582069A CN 112582069 A CN112582069 A CN 112582069A CN 202011573380 A CN202011573380 A CN 202011573380A CN 112582069 A CN112582069 A CN 112582069A
- Authority
- CN
- China
- Prior art keywords
- dendritic
- dendritic spine
- rate
- spine
- width
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000009792 diffusion process Methods 0.000 title claims abstract description 48
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 39
- 210000002569 neuron Anatomy 0.000 title claims abstract description 26
- 230000024573 dendritic spine development Effects 0.000 title claims abstract description 18
- 210000003520 dendritic spine Anatomy 0.000 claims abstract description 178
- 238000004088 simulation Methods 0.000 claims abstract description 55
- 210000001787 dendrite Anatomy 0.000 claims abstract description 31
- 210000004556 brain Anatomy 0.000 claims abstract description 19
- 230000008569 process Effects 0.000 claims abstract description 17
- 230000002159 abnormal effect Effects 0.000 claims abstract description 7
- 230000008506 pathogenesis Effects 0.000 claims abstract description 7
- 210000001320 hippocampus Anatomy 0.000 claims abstract description 3
- 239000012190 activator Substances 0.000 claims description 31
- 239000000758 substrate Substances 0.000 claims description 29
- 239000003112 inhibitor Substances 0.000 claims description 28
- 230000015556 catabolic process Effects 0.000 claims description 15
- 238000006731 degradation reaction Methods 0.000 claims description 15
- 235000001674 Agaricus brunnescens Nutrition 0.000 claims description 13
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 210000004292 cytoskeleton Anatomy 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 208000012902 Nervous system disease Diseases 0.000 claims description 5
- 208000025966 Neurological disease Diseases 0.000 claims description 3
- 230000003197 catalytic effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000004913 activation Effects 0.000 claims 1
- 230000000977 initiatory effect Effects 0.000 claims 1
- 241001465754 Metazoa Species 0.000 abstract description 5
- 201000010099 disease Diseases 0.000 abstract description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 5
- 238000002474 experimental method Methods 0.000 abstract description 5
- 238000011160 research Methods 0.000 abstract description 5
- 210000005036 nerve Anatomy 0.000 abstract 1
- 238000007792 addition Methods 0.000 description 21
- 241000700159 Rattus Species 0.000 description 19
- 208000032612 Glial tumor Diseases 0.000 description 9
- 206010018338 Glioma Diseases 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 230000008602 contraction Effects 0.000 description 6
- 230000000971 hippocampal effect Effects 0.000 description 6
- 210000004295 hippocampal neuron Anatomy 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 210000005013 brain tissue Anatomy 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005653 Brownian motion process Effects 0.000 description 2
- 238000005537 brownian motion Methods 0.000 description 2
- 238000006555 catalytic reaction Methods 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 210000003710 cerebral cortex Anatomy 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000010172 mouse model Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 239000012192 staining solution Substances 0.000 description 2
- 238000005406 washing Methods 0.000 description 2
- VHUUQVKOLVNVRT-UHFFFAOYSA-N Ammonium hydroxide Chemical compound [NH4+].[OH-] VHUUQVKOLVNVRT-UHFFFAOYSA-N 0.000 description 1
- 208000001914 Fragile X syndrome Diseases 0.000 description 1
- 108010010803 Gelatin Proteins 0.000 description 1
- 235000011114 ammonium hydroxide Nutrition 0.000 description 1
- 238000010171 animal model Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000000170 cell membrane Anatomy 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003413 degradative effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002964 excitative effect Effects 0.000 description 1
- 229920000159 gelatin Polymers 0.000 description 1
- 239000008273 gelatin Substances 0.000 description 1
- 235000019322 gelatine Nutrition 0.000 description 1
- 235000011852 gelatine desserts Nutrition 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 210000001879 hippocampal ca1 region Anatomy 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006386 memory function Effects 0.000 description 1
- 238000001000 micrograph Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000007261 regionalization Effects 0.000 description 1
- 230000003248 secreting effect Effects 0.000 description 1
- 230000028327 secretion Effects 0.000 description 1
- 230000018448 secretion by cell Effects 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 210000004092 somatosensory cortex Anatomy 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 210000000225 synapse Anatomy 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Landscapes
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种基于反应扩散模型的神经元树突棘发育模式建立方法，该方法以大鼠大脑海马区神经元树突棘形态制定树突棘形状分类标准，并利用反应扩散模型模拟单一树突棘在不同条件下的生长过程，使用形状分类标准对仿真结果进行分类，进一步利用反应扩散模型模拟一段树突在不同条件下生长不同密度的树突棘的过程，继而得到不同模式的神经元树突棘。与现有的研究方法相比，不仅有效避免了大量的动物实验；同时本发明首次将反应扩散模型应用到神经元树突棘发育模式的建立过程中，为后续研究树突棘模式异常引发疾病的发病机制奠定了基础。The invention provides a method for establishing a developmental pattern of neuron dendritic spines based on a reaction diffusion model. The method formulates a dendritic spine shape classification standard based on the morphology of neuron dendritic spines in the hippocampus of the rat brain, and uses the reaction diffusion model to simulate a single dendritic spine shape. The growth process of dendritic spines under different conditions, the shape classification standard was used to classify the simulation results, and the reaction diffusion model was further used to simulate the process of a section of dendrites growing dendritic spines with different densities under different conditions, and then different patterns of nerves were obtained. Metadendritic spines. Compared with the existing research methods, it not only effectively avoids a large number of animal experiments; at the same time, the present invention applies the reaction diffusion model to the establishment of the neuron dendritic spine development pattern for the first time, which is for the follow-up study of abnormal dendritic spine pattern causing diseases. The pathogenesis has laid the foundation.
Description
Technical Field
The invention relates to the field of neuron dendritic spine biological patterns, in particular to a neuron dendritic spine development pattern establishing method based on a reaction diffusion model and application of the neuron dendritic spine development pattern establishing method in disease mechanism analysis.
Background
Dendritic spines are minute processes on the dendrites of neurons, widely present in the dendrites of higher animals, and play an important role in the formation of the dendritic synapses of most excitatory axes. Generally dendritic spines are divided into four patterns: elongated, columnar, mushroom, and branched. Wherein, the slender dendritic spines have high plasticity and are related to learning; mushroom-type dendritic spines exhibit poor plasticity and are involved in memory function. Abnormal patterns of dendritic spines are associated with a variety of neurological disorders. For example, there are too many elongated dendritic spines in the somatosensory cortex of fragile X syndrome model mice. The density of dendritic spines in the dendrites of neurons of glioma model mice was significantly reduced. Therefore, it is important to understand the pathogenesis of these diseases and find a therapeutic method by studying the morphology and density of the dendritic spines.
Currently, the main approaches to dendritic spine pattern studies are: the cerebral cortex of the experimental animals and the cerebral cortex of the control animals were compared by brain section. Such studies have major drawbacks: (1) a large number of anatomical animal experiments need to be performed; (2) only one factor influencing the dendritic spine pattern is proposed at a time, and the formation of the dendritic spine is a dynamic process, involves various chemical reactions and is regulated and controlled by various factors.
Therefore, in order to systematically study the mechanism of development of various neurological diseases associated with the dendritic spine patterns, it is necessary to introduce numerical simulation methods to simulate the growth process of the different patterns of dendritic spines.
There are two main types of numerical simulation methods currently used to study the dendritic spine pattern: (1) the volume of dendritic spines is simulated using a brownian motion model, but the brownian motion model describes a random phenomenon, and the pattern formation of dendritic spines is a gene-and environment-regulated process, rather than a random process; (2) the relationship between the forces generated by the cytoskeleton and the deformation of the cell membrane was studied, but this method lacks an explanation for the reason why the cytoskeleton generates different forces. In addition, in the process of carrying out numerical simulation, the method distinguishes different patterns of dendritic spines in a manual identification mode, and subjective errors are inevitably introduced. Therefore, developing a method for establishing a neuron dendritic spine development mode based on a reaction diffusion model, and making morphological classification indexes based on morphological differences of different dendritic spines are very important for researching pathogenesis of diseases caused by abnormal dendritic spine modes.
Disclosure of Invention
The invention provides a method for establishing a neuron dendritic spine development mode based on the defects of the existing mathematical simulation method simulation neuron dendritic spine mode.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for establishing a neuron dendritic spine development mode based on a reaction diffusion model comprises the following steps:
s1, identifying dendritic spines in four modes of mushroom type, column type, long and thin type and branch type in the rat brain hippocampal neuron microscopic image, and measuring the head width w of the dendritic spines in three modes of mushroom type, column type and long and thin typeheadWidth w of neckneckHeight h;
s2, according to the width w of the head of the dendritic spineheadWidth w of neckneckAnd a height h, calculating a relative average width RAW and a relative shrinkage width RCW of the dendritic spines in the three modes in step S1;
s3, establishing classification standards of dendritic spines in different modes;
s4, establishing a reaction diffusion model, setting simulation initial conditions, simulating the growth process of a single dendritic spine by setting variable parameters of the model, applying the classification standard of the step S3 to the reaction diffusion model, and classifying the shape of the dendritic spine in the simulation result;
s5, performing a simulation of a plurality of dendritic spines: and (5) simulating the process of growing a plurality of dendritic spines on a section of dendrite by using the reaction diffusion model in the step S4, and setting variable parameters of the model to obtain different patterns of dendritic spines.
Further, in step S1, the head width w of the dendritic spineheadRefers to the width of the widest part of the dendritic spine away from the dendrites; width w of neckneckRefers to the width of the dendritic spine at its narrowest point in the half near the dendrites.
Further, in step S2, the relative average width RAW of the dendritic spines is calculated as follows:
the relative contraction width RCW of the dendritic spines is calculated as follows:
wherein, wheadWidth of head of dendritic spine, wneckIs the neck width of the dendritic spine and h is the dendritic spine height.
Further, in step S3, the dendritic spine classification indexes in different patterns are as follows:
(a) branched dendritic spines, if they comprise a branched structure; (b) if the relative average width RAW of the dendritic spines is less than 0.4, the dendritic spines are slender dendritic spines; (c) if the relative average width RAW of the dendritic spines is more than 0.4 and the relative contraction width RCW is less than 0.25, the dendritic spines are columnar dendritic spines; (d) a mushroom-type dendritic spine if the relative average width RAW is greater than 0.4 and the relative shrinkage width RCW is greater than 0.25.
Further, the reaction diffusion models in the steps S4 and S5 are as follows:
wherein A, H, S, Y represents the concentration of activator, inhibitor, substrate and cytoskeleton at each point in the simulation region; c is the catalytic rate; μ is the degradation rate of the activator;ρAis the rate of production of the activator; dAIs the diffusion rate of the activator; upsilon is the degradation rate of the inhibitor; rhoHIs the rate of production of the inhibitor; dHIs the diffusion rate of the inhibitor; c. C0A rate of substrate production for the environment; gamma is the degradation rate of the substrate; ε is the rate of substrate consumption; dSIs the diffusion rate of the substrate; e is the degradation rate of the labeling substance; f is a positive feedback coefficient; deltaAIndicates the rate of exogenous activator addition, deltaHIndicating the rate of exogenous inhibitor addition.
Further, in step S4, the simulation initial condition refers to a simulation area size of 100 × 100 pixels, and an area with 10 × 5 pixels disposed in the middle above the area is used as a growth start state of the dendritic spines.
Further, in step S4, the dendritic spine region variable initial state is a ═ 2, H ═ 0.02, S ═ 1, Y ═ 1, the background region variable initial state is a ═ 0.001, H ═ 0.001, S ═ 1, Y ═ 0; at exogenous activator addition rate deltaAExogenous inhibitor addition rate deltaHSubstrate consumption rate, epsilon, as a model variable parameter, and the other parameters were fixed parameters set at c 0.002, μ 0.16, upsilon 0.04, ρA＝0.01，ρH＝0.00005，c0＝0.02，γ＝0.02，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10; the values of the variable parameters of the model are respectively epsilon 0.01-0.90 and deltaA＝0，0.005，0.01，0.015，0.02，δH＝0，0.000025，0.00005，0.000075，0.0001。
Further, in step S5, a dendrite simulation is performed by using a reaction diffusion model, the size of the dendrite simulation area is 150 × 200 pixels, an area with 5 × 10 pixels in the middle of the left side of the area is used as the growth start state of the dendrite, and the values of model parameters are: c is 0.002, μ is 0.16, upsilon is 0.04, ρA＝0.03，ρH＝0.0001，δA＝0，δH＝0，c0＝0.2，γ＝0.02，ε＝0.017，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10; secondly, the first step is to carry out the first,taking the result of the dendritic simulation as the initial condition of the plurality of dendritic spine simulations, carrying out the plurality of dendritic spine growth simulations, and adding the exogenous activator at a rate deltaAExogenous inhibitor addition rate deltaHThe substrate consumption rate epsilon as a model variable parameter, and the other parameters as fixed parameters, c 0.002, μ 0.16, upsilon 0.04, ρA＝0.02，ρH＝0.00005，c0＝0.05，γ＝0.02，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10, and the values of the model variable parameters are respectively 0.5, 1, 1.5, 2, 2.5, and δA＝0，0.01，0.02，0.03，0.04，δH＝0，0.00005，0.0001，0.00015，0.0002。
Furthermore, the method for establishing the development pattern of the neuron dendritic spines can be applied to the analysis and research of pathogenesis of nervous system diseases caused by abnormal dendritic spine patterns.
Compared with the prior art, the invention has the following advantages and effects:
1. according to the method for establishing the neuron dendritic spine development mode, the dendritic spine shape division standard is established through the measurement data of the neuron dendritic spine microscopic image in the hippocampal region of the brain of the rat of the control group, the standard is applied to the reaction diffusion model, and the dendritic spine shapes in the simulation result are automatically classified.
2. According to the method, simulation initial conditions are set according to the actual growth environment of the neuron dendritic spines, the shape of a single dendritic spine and the density of a plurality of dendritic spines are researched by utilizing a reaction diffusion model, and the neuron dendritic spines in different modes are obtained by setting fixed parameters of the model and adjusting variable parameters of the model; meanwhile, the reaction diffusion model is applied to the establishment process of the development mode of the dendritic spines of the neurons for the first time, and a foundation is laid for the subsequent research of the pathogenesis of diseases caused by abnormal dendritic spine modes.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a flowchart of a method for establishing a neuron dendritic spine development pattern based on a reaction diffusion model according to an embodiment of the present invention.
FIG. 2 is a microscopic image (a) of dendritic spines in four shapes on neurons in hippocampal regions of brain of a rat in a control group according to the present invention and a schematic diagram (b) of data to be measured;
FIG. 3 is a diagram of the initial condition setting of the single dendritic spine simulation of the present invention.
FIG. 4 is a diagram of the simulation results of the formation of the dendrite spine shape pattern under different variable model parameter settings according to the present invention.
FIG. 5 is a diagram (a) illustrating initial condition setting of dendrite simulation according to the present invention; a simulation initial condition setting map (b) in which a plurality of dendritic spines are grown on dendrites.
FIG. 6 is a graph of simulation results of formation of dendritic spine density patterns under different variable model parameter settings according to the present invention.
FIG. 7 is a brain hippocampal neuron microscopic image (a) of a control rat according to the present invention; glioma group rat brain hippocampal neuron microscopic image (b); and (c) a glioma rat brain hippocampus neuron dendritic spine mode simulation result graph obtained by utilizing the neuron dendritic spine development mode establishing method disclosed by the invention.
Detailed Description
The present invention will be described in further detail with reference to examples, which are illustrative of the present invention and are not to be construed as being limited thereto.
Example 1: as shown in the flow chart of fig. 1: a method for establishing a neuron dendritic spine development mode based on a reaction diffusion model comprises the following steps:
s1, identificationDendritic spines in four modes of mushroom type, column type, long and thin type and branch type in the microscopic image of the hippocampal neurons of the brain of the other rat are measured, and the head width w of the dendritic spines in the three modes of the mushroom type, the column type and the long and thin type is measuredheadWidth w of neckneckHeight h;
s2, according to the width w of the head of the dendritic spineheadWidth w of neckneckAnd a height h, calculating a relative average width RAW and a relative shrinkage width RCW of the dendritic spines in the three modes in step S1;
s3, establishing classification standards of dendritic spines in different modes;
s4, establishing a reaction diffusion model, setting simulation initial conditions, simulating the growth process of a single dendritic spine by setting variable parameters of the model, applying the classification standard of the step S3 to the reaction diffusion model, and classifying the shape of the dendritic spine in the simulation result;
s5, performing a simulation of a plurality of dendritic spines: and (5) simulating the process of growing a plurality of dendritic spines on a section of dendrite by using the reaction diffusion model in the step S4, and setting variable parameters of the model to obtain different patterns of dendritic spines.
Specifically, in step S1, the neurons used were SD rat hippocampal CA1 region neurons.
The experimental method is as follows: taking out brain tissue, immediately placing the brain tissue into Golgi-Cox staining solution, storing the brain tissue in dark for two weeks, wherein the staining solution is changed every 48 hours, then cutting the brain tissue into brain slices with the thickness of 150 mu m by using a vibratory microtome, placing the brain slices on a 2% gelatin-coated glass slide, staining the brain slices with ammonia water for 60 minutes, washing the brain slices for three times, soaking the brain slices for 30 minutes by using a kodak fixing solution, washing, dehydrating, cleaning, flaking, observing and shooting three-level dendrites of neurons in a hippocampal area of a rat brain under a 100 x objective lens, and finding out 3 mushroom-type dendrites, 3 columnar dendrites, 2 elongated spines and 1 branched dendrite spines in a microscopic image (as shown in figure 2(a), wherein iii, vii and ix represent mushroom types, ii, vi and viii represent columnar types, iv and v represent branched slender spines, and i represents branched spines); the head width of each type of dendritic spine except for the branched type was measured as shown in FIG. 2(b)wheadWidth w of neckneckAnd a height h. Wherein the width w of the head of the dendritic spineheadRefers to the width of the widest part of the dendritic spine away from the dendrites; width w of neckneckRefers to the width of the dendritic spine at its narrowest point in the half near the dendrites.
In step S2, defining a relative average width RAW for measuring the overall width of the dendritic spines, defining a relative contracted width RCW for measuring the difference in width between the head and neck of the dendritic spines, calculating the Relative Average Width (RAW) and the Relative Contracted Width (RCW) using the following formulas;
according to the measurement result of step S1, the relative average widths RAW of the 3 mushroom-type dendritic spines are calculated to be 0.54, 0.51, 0.65, respectively, and the relative contraction widths RCW are calculated to be 0.38, 0.48, 0.74, respectively; the relative average widths RAW of the 3 columnar dendritic spines are respectively 0.52, 0.96 and 0.93, and the relative contraction widths RCW are respectively 0.09, -0.41 and-0.42; the relative average width RAW of the 2 elongated dendritic spines is 0.30 and 0.32 respectively, and the relative contraction width RCW is 0.16 and 0.30 respectively.
In step S3, the index range corresponding to each dendritic spine shape is divided according to the calculation result in step S2: elongated dendritic spines RAW less than 0.4; the RAW of the columnar dendritic spines is more than 0.4, and the RCW is less than 0.25; mushroom type dendritic spines RAW is greater than 0.4 and RCW is greater than 0.25. From this, a classification criterion for the shape of the dendritic spines can be established: (a) branched dendritic spines, if they comprise a branched structure; (b) if the relative average width RAW of the dendritic spines is less than 0.4, the dendritic spines are slender dendritic spines; (c) if the relative average width RAW of the dendritic spines is more than 0.4 and the relative contraction width RCW is less than 0.25, the dendritic spines are columnar dendritic spines; (d) a mushroom-type dendritic spine if the relative average width RAW is greater than 0.4 and the relative shrinkage width RCW is greater than 0.25.
In steps S4 and S5, the reaction diffusion model used is as follows:
a, H, S, Y represents the activator concentration, inhibitor concentration, substrate concentration, and cytoskeleton concentration at each point in the simulated region. Equation (1) is used to calculate the rate of change of the concentration of the activator, and four factors are responsible for the change in the concentration of the activator: autocatalytic, degradative, secretory and diffusive effects of the activator. Equation (2) is used to calculate the rate of change of inhibitor concentration, and there are four factors that cause the change in inhibitor concentration: cross-catalysis, degradation, cellular secretion and diffusion of activators. Equation (3) is used to calculate the rate of change of substrate concentration, and there are four factors that cause the change in substrate concentration: environmental secretion, degradation, depletion of cellular activity and diffusion. Equation (4) is used to calculate the rate of change of the cytoskeleton concentration, and there are three factors that cause the change of the cytoskeleton concentration: cross-catalysis, degradation of the activator and positive feedback with saturation suppression. The parameters in formula (1), formula (2), formula (3) and formula (4) are explained as follows: c is the catalytic rate; μ is the degradation rate of the activator; rhoAIs the rate of production of the activator; dAIs the diffusion rate of the activator; upsilon is the degradation rate of the inhibitor; rhoHIs the rate of production of the inhibitor; dHIs the diffusion rate of the inhibitor; c. C0A rate of substrate production for the environment; gamma is the degradation rate of the substrate; epsilon is the rate at which the cell consumes the substrate; dSIs the diffusion rate of the substrate; e is the degradation rate of the labeling substance; f is a positive feedback coefficient; deltaAIndicates the rate of exogenous activator addition, deltaHIndicates the rate of addition of exogenous inhibitor, delta in the above parametersA、δHEpsilon is a variable parameter of the model, other parameters are fixed parameters of the model, and the substrate consumption rate epsilon and the exogenous activator addition rate delta of the cells are controlled in the simulation processAAnd rate of addition of exogenous inhibitorRate deltaHDifferent patterns of the neuron dendritic spines are obtained.
Specifically, in step S4, the specific method for simulating the growth process of a single dendritic spine is as follows:
as shown in fig. 3, the simulation initial conditions are set, the simulation area size is 100 × 100 pixels, and the area with 10 × 5 pixels in the middle above the area is used as the growth start state of the dendritic spines. The initial state of the dendrite spine region variable is a ═ 2, H ═ 0.02, S ═ 1, Y ═ 1, the initial state of the background region variable is a ═ 0.001, H ═ 0.001, S ═ 1, Y ═ 0; selection of deltaA(rate of exogenous activator addition), deltaH(rate of exogenous inhibitor addition), epsilon (rate of substrate consumption) as model variable parameters, and other parameters were fixed parameters set at c 0.002, μ 0.16, upsilon 0.04, ρA＝0.01，ρH＝0.00005，c0＝0.02，γ＝0.02，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10; the values of the variable parameters of the model are respectively epsilon 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, … …, 0.90 and deltaA＝0，0.005，0.01，0.015，0.02，δ H0, 0.000025, 0.00005, 0.000075, 0.0001, simulation is implemented using C + + language and runs in Visual Studio2015 environment, space step size is 0.3, and all variables are dimensionless quantities. Measuring the width w of the head of the dendritic spine in the simulation resultheadWidth w of neckneckAnd h, calculating the relative average width RAW and the relative shrinkage width RCW of each simulated dendritic spine, and applying the dendritic spine classification standard in the step S3 to classify the simulated dendritic spine, wherein partial simulation results, corresponding indexes and classification results are shown in FIG. 4.
The analysis was as follows: in step S4, by the reaction model variable parameter setting, the influence of the model variable parameter on the shape of the single dendritic spine can be found as follows: as the substrate consumption rate epsilon increases, the dendritic web shape simulation results are mushroom type, column type, long and thin type, and branch type in sequence (as shown in fig. 4 (a)); rate of addition of exogenous activator deltaAPart of the original pattern of elongated dendritic spines is converted into columnar dendritic spines, and part of the original pattern of elongated dendritic spines is converted into columnar dendritic spinesThe dendritic spines whose pattern is columnar are converted into mushroom-type dendritic spines (as shown in fig. 4 (b)); with exogenous inhibitor deltaHThe rate of addition was increased, and part of the dendritic spines, which were originally in a branched form, were converted into elongated dendritic spines (as shown in fig. 4 (c)). Therefore, during the simulation, the parameter δ is varied by the modelA(rate of exogenous activator addition), deltaHThe settings of (exogenous inhibitor adding rate) and epsilon (substrate consumption rate) can realize the conversion between dendritic spines with different shapes.
Specifically, in step S5, the specific method for simulating the growth process of the plurality of dendritic spines is as follows:
firstly, performing dendrite simulation by using a reaction diffusion model, as shown in fig. 5(a), the size of a dendrite simulation area is 150 × 200 pixels, an area with 5 × 10 pixels in the middle of the left side of the area is used as a growth starting state of dendrites, and the values of model parameters are respectively as follows: c is 0.002, μ is 0.16, upsilon is 0.04, ρA＝0.03，ρH＝0.0001，δA＝0，δH＝0，c0＝0.2，γ＝0.02，ε＝0.017，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10;
next, as shown in fig. 5(b), a plurality of dendritic spine growth simulations were performed using the result of the dendritic spine simulation as an initial condition for the plurality of dendritic spine simulations, and δ was selectedA(rate of exogenous activator addition), deltaH(rate of exogenous inhibitor addition), epsilon (rate of substrate consumption) as model variable parameters, and other parameters were fixed parameters set at c 0.002, μ 0.16, upsilon 0.04, ρA＝0.02，ρH＝0.00005，c0＝0.05，γ＝0.02，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10, and the values of the model variable parameters are respectively 0.5, 1, 1.5, 2, 2.5, and δA＝0，0.01，0.02，0.03，0.04，δH＝0，0.00005，0.0001，0.00015，0.0002。
The simulation results are shown in fig. 6, and the analysis shows the following: by setting variable parameters of the reaction model, the variable parameters of the model can be obtained for a plurality of dendritic spinesThe effect of density is: rate of addition of exogenous activator deltaAThe density of dendritic spines is increased in sequence, and the dominant type is changed from a long and thin type to a column type and then to a mushroom type (as shown in figure 6 (a)); with exogenous inhibitor deltaHThe adding rate is increased, the density of the dendritic spines is reduced in sequence, and the dominant type is changed from mushroom type to column type and then to slender type (as shown in fig. 6 (b)); as the substrate consumption rate ∈ increases, the density of dendritic spines decreases in turn, and changes from being predominantly long-length to branched (as shown in fig. 6 (c)). Therefore, during the simulation, the parameter δ is varied by the modelA(rate of exogenous activator addition), deltaHThe setting of (exogenous inhibitor adding rate) and epsilon (substrate consumption rate) can realize the conversion between dendritic spine patterns with different density and shape distribution.
Specifically, as shown in FIG. 6(a), when δH0.00005, ∈ 1; rate of addition of exogenous activator deltaAGradually increasing, part of the dendritic spines with original pattern of slender type is converted into columnar dendritic spines, and part of the dendritic spines with original pattern of columnar type is converted into mushroom type dendritic spines, so that when delta is increasedAWhen the value is 0, a dendrite spine pattern mainly having a long and thin shape is obtained, and when the value is δAWhen the number of the grooves is increased to 0.01, a main columnar dendritic spine pattern is obtained, and when the number of the grooves is increased to deltaAWhen the number of the grooves is increased to 0.04, a mushroom-type main dendritic spine pattern can be obtained.
When delta is shown in FIG. 6(b)A0.01,. epsilon.1, with exogenous inhibitor deltaHThe density of branched dendritic spines gradually decreases with increasing addition rate, therefore, when delta is increasedHWhen δ is 0.00005, a dendrite spine pattern mainly having a long and thin shape can be obtainedHWhen the thickness is increased to 0.0001, the dendritic spine pattern which is all long and thin can be obtained.
When delta is shown in FIG. 6(c)A＝0.01，δHWhen the substrate consumption rate epsilon is increased, the dendritic spine pattern is sequentially changed from a column type to a long and thin type and a branch type, when the epsilon is 1.0, the dendritic spine pattern mainly with the long and thin type can be obtained, and when the epsilon is increased to 2.0, the dendritic spine pattern mainly with the branch type can be obtained.
Further, as shown in fig. 7, the effectiveness of the method is verified by comparing a control group rat and a glioma rat dendritic spine imaging experiment with a glioma rat hippocampal neuron dendritic spine mode simulation graph obtained by using the mode establishing method of the invention.
A group of a plurality of normal SD rats are taken as a control group, a group of a plurality of SD rats suffering from glioma is taken as a glioma group, two groups of rat hippocampal brain slices are prepared, and a dendritic spine pattern on a neuron tertiary dendrite of a rat brain hippocampal area is observed under a microscope. FIGS. 7(a), (b) show microscope images of neuronal dendritic spines in hippocampal regions of normal rats and glioma rats.
The graph (c) is a glioma rat brain hippocampal neuron dendritic spine mode simulation graph obtained by utilizing the mode establishing method of the invention, and the variable parameter delta of the model is obtained at the momentA＝0.01，δH0.0001, ∈ 1, mode fixed parameter c 0.002, μ0.16, ν 0.04, ρA＝0.02，ρH＝0.00005，c0＝0.05，γ＝0.02，DA＝0.02，DH＝0.26，DS0.06, 0.0035, 0.1, 10.
Example 2: the method for establishing the development pattern of the neuron dendritic spines, which is described in the embodiment 1, is applied to the analysis and research of pathogenesis of nervous system diseases caused by abnormal dendritic spine patterns. Compared with the existing research and analysis method, a large number of animal experiments can be avoided.
In addition, it should be noted that the specific embodiments described in the present specification may differ in the shape of the components, the names of the components, and the like. All equivalent or simple changes of the structure, the characteristics and the principle of the invention which are described in the patent conception of the invention are included in the protection scope of the patent of the invention. Various modifications, additions and substitutions for the specific embodiments described may be made by those skilled in the art without departing from the scope of the invention as defined in the accompanying claims.
Claims (9)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011573380.XA CN112582069B (en) | 2020-12-25 | 2020-12-25 | A method for establishing the developmental pattern of neuronal dendritic spines based on the reaction-diffusion model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011573380.XA CN112582069B (en) | 2020-12-25 | 2020-12-25 | A method for establishing the developmental pattern of neuronal dendritic spines based on the reaction-diffusion model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112582069A true CN112582069A (en) | 2021-03-30 |
| CN112582069B CN112582069B (en) | 2023-01-17 |
Family
ID=75140121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011573380.XA Active CN112582069B (en) | 2020-12-25 | 2020-12-25 | A method for establishing the developmental pattern of neuronal dendritic spines based on the reaction-diffusion model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112582069B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106754669A (en) * | 2016-11-23 | 2017-05-31 | 河海大学常州校区 | The preparation method and preparation system of the multi-cellular structure based on Reaction-Diffusion Models |
| CN110904114A (en) * | 2019-12-16 | 2020-03-24 | 东北师范大学 | A coding sequence of Gnb2 gene and application of coding sequence of Gnb2 gene in constructing mouse model |
| CN111465852A (en) * | 2017-12-28 | 2020-07-28 | 国际商业机器公司 | Neural model generation for personalized drug therapy |
-
2020
- 2020-12-25 CN CN202011573380.XA patent/CN112582069B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106754669A (en) * | 2016-11-23 | 2017-05-31 | 河海大学常州校区 | The preparation method and preparation system of the multi-cellular structure based on Reaction-Diffusion Models |
| CN111465852A (en) * | 2017-12-28 | 2020-07-28 | 国际商业机器公司 | Neural model generation for personalized drug therapy |
| CN110904114A (en) * | 2019-12-16 | 2020-03-24 | 东北师范大学 | A coding sequence of Gnb2 gene and application of coding sequence of Gnb2 gene in constructing mouse model |
Non-Patent Citations (2)
| Title |
|---|
| MOHAMMAD NAZRUL ISHLAM PATOARY: "Parallel Stochastic Discrete Event Simulation of Calcium Dynamics in Neuron", 《ACM》 * |
| 郭艺娜: "基于气态反应模型的三维肺分枝形态发生的研究", 《中国优秀博硕士学位论文全文数据库（博士）医药卫生科技辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112582069B (en) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hartenstein et al. | Development of adult sensilla on the wing and notum of Drosophila melanogaster | |
| CN106202997B (en) | A kind of cell division detection method based on deep learning | |
| Levin | Bioelectric mechanisms in regeneration: unique aspects and future perspectives | |
| Li et al. | Fibroblasts and myofibroblasts in wound healing: force generation and measurement | |
| Watanabe et al. | Real-time analysis of ureteric bud branching morphogenesis in vitro | |
| Shefi et al. | Morphological characterization of in vitro neuronal networks | |
| Zhang et al. | Modeling neurological disorders using brain organoids | |
| Jackson et al. | Topological analysis of multicellular complexity in the plant hypocotyl | |
| Buchmann et al. | Sizing it up: the mechanical feedback hypothesis of organ growth regulation | |
| Buxhoeveden et al. | Quantitative analysis of cell columns in the cerebral cortex | |
| Panfilio et al. | High plasticity in epithelial morphogenesis during insect dorsal closure | |
| Polackwich et al. | Traction force and tension fluctuations in growing axons | |
| Scheibel | Some structural and developmental correlates of human speech | |
| Eglen et al. | Lateral cell movement driven by dendritic interactions is sufficient to form retinal mosaics | |
| Al-Absi et al. | Layers II/III of prefrontal cortex in Df (h22q11)/+ mouse model of the 22q11. 2 deletion display loss of parvalbumin interneurons and modulation of neuronal morphology and excitability | |
| Nyaanga et al. | Changes in body shape implicate cuticle stretch in C. elegans growth control | |
| Aćimović et al. | Modeling of neuronal growth in vitro: comparison of simulation tools NETMORPH and CX3D | |
| CN112582069A (en) | Neuron dendritic spine development mode establishing method based on reaction diffusion model | |
| Rahman et al. | Biomechanical modeling of cell chirality and symmetry breaking of biological systems | |
| Bertling et al. | Methods for three-dimensional analysis of dendritic spine dynamics | |
| Wang et al. | Experience-dependent structural plasticity of pyramidal neurons in the developing sensory cortices | |
| Reiner et al. | Modeling human neuronal migration deficits in 3D | |
| van Ooyen | Network formation through activity-dependent neurite outgrowth: a review of a simple model of homeostatic structural plasticity | |
| Rabinowicz et al. | Neocortical maturation during adolescence: change in neuronal soma dimension | |
| Casanova et al. | Minicolumnar width: comparison between supragranular and infragranular layers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |