Research Article  Open Access
Benjamin C. Yan, Johnson F. Yan, "A TreeLike Model for Brain Growth and Structure", Journal of Biophysics, vol. 2013, Article ID 241612, 5 pages, 2013. https://doi.org/10.1155/2013/241612
A TreeLike Model for Brain Growth and Structure
Abstract
The FloryStockmayer theory for the polycondensation of branched polymers, modified for finite systems beyond the gel point, is applied to the connection (synapses) of neurons, which can be considered highly branched “monomeric” units. Initially, the process is a linear growth and treelike branching between dendrites and axons of nonselfneurons. After the gel point and at the maximum “tree” size, the treelike model prescribes, on average, one pair of twin synapses per neuron. About 13% of neurons, “unconnected” to the maximum tree, migrate to the surface to form cortical layers. The number of synapses in each neuron may reach 10000, indicating a tremendous amount of flexible, redundant, and neuroplastic loopforming linkages which can be preserved or pruned by experience and learning.
1. Introduction
The molecular weight distribution (MWD) of polycondensation of branchedchain monomers of the type has been derived classically by Flory [1] and generalized by Stockmayer [2]. Here, is the number of functional groups, or “functionality” of group A. Mathematically, this widely quoted distribution function has been treated as powerseries distribution and compound distribution that provides a simple concept; that is, singleparameter expressions of number and weightaverage “degrees of polymerization” (DP) are sufficient to generate the entire MWD for branched polymers [3]. Moreover, using a cascade formulation involving functionals and probability generating functions (PGF), this distribution can be extended to finite systems [4].
Here, the previously derived properties of this finite distribution are applied to synapse formation in the brain. A neuron has multiple dendritic processes and an axon, which can also be branched. Neurons are generally three or more orders of magnitude greater in size than molecular units. However, the functionality of a neuron may be 10^{3} times larger than that of typical branched molecules (). This large functionality also means there is great accessibility to the connection sites, and the long flexible axons offer a favorable condition for connection between neurons. The “treelike,” or “ringfree,” assumption in the FloryStockmayer theory can be satisfied by the initial linkage of headtotail linear chains and followed by a “treelike branching.” A neuron itself can be considered a small tree. Similarly, the peripheral nervous system (PNS) also resembles a tree made of the nerve bundles, which can be as large as 1.5 meters. The Finite FloryStockmayer theory (FFST) deals with numbers of highly branched repeat units and their associationdissociation mechanism. Therefore it is applicable to the statistical treatment of brain growth, neuronal connectivity, and information transmission. Linear growth and subsequent treelike branching allow two simple equations to be applied up to a maximum tree size constrained by the FFST.
In practical, real and natural systems, lignin in wood is a treelike molecule [5], and a major class of corals is also treelike [6]. Thus the applicability of the FS theory in molecular, cellular or animal size scales is equally sound if the repeating branched units can be properly identified.
Early brain growth is very robust and rapid. In the human brain the crosslinking between neurons and the growth of the linked neurons (polyneurons) proceed with a rapid rate. Because of their highly branched dendrites and axons, there are three types of crosslinking: axodendritic (AD), axosomatic (AS), and axoaxonic (AA) synapses, in the order of decreasing abundance. In addition, there are at least three forms of crosslinks. The first are endlinks which are links between the tips of an axon and a dendritic spine. The second form is a splitend link, which is shaped like the Greek letter on top of (the shape of a pair of dendritic spines), designated linking. The latter conformation represents a pair of twin synapses. The third form is the most abundant “” shaped, mostly AD linkages.
Structural and functional features of neurons set the limits and modes of various crosslinks. For example, there is no dendrite to dendrite (DD) synapse because dendrites have been shown to demonstrate selfavoidance even in fruit fly brains [7]. With respect to a growing polyneuron tree, the “ringfree” assumption is also a selfavoidance rule. A conceptual illustration furnishes a visual model of a pair of twin synapses by removing the axonic “spaghetti maze” for a clearer view [8].
The rate of growth and rate of crosslinking in the brain occur exponentially in the months before and after birth in the human brain [9] and in the cortex of laboratory animals [10]. Initial linear growth and treelike branching appear to be the most efficient pathway to grow and expand rapidly in a confined or open space [11].
The following is an ideal human brain model with neurons assumed to be homogeneous, with identical size, shape, and number of crosslinking sites.
There are two basic assumptions in the FFST. First, all linkages (synapses) are formed with equal probability, or the same “extent of reaction.” Second, no ring is formed between branches of a growing chain. The second assumption is also termed the “treelike model” [4]. The brain growth process is discussed with limiting forms of equations derived previously, and upper bounds of the variables and parameters are strictly constrained.
2. Limiting Equations and Numerical Examples
Symbols and abbreviations derived previously [4]. = system size = number of neurons (in a brain or a “module”). = functionality = number of connection sites per monomer (neuron) = 10000 or in thousands [11]. = number of generations in the “transformer” probability generating functions (PGF) in the cascade formulation [4]. DP = degree of polymerization = number of connected neurons in a growing chain = . = extent of reaction or fraction of connected synapse sites. = which is approximately for large = crosslinking index. = numberaverage DP. = weightaverage DP. = 1 is value of at gel point. = 2, maximum value of allowed in treelike model. = gel fraction at .
At large values, most of the finite system equations derived previously [4] become independent of . The FFST is represented by the two averages: The last equation is derived and proven in Supplementary Material session available online at http://dx.doi.org/10.1155/2013/241612. In its application, (2a) can be written in a more compact form: In particular,
The maximum value for is , that is, a giant gel particle with neurons all connected in this “polyneuron” tree. With this value for at = 2 in (2b), the result is simply
Equation (4a) expresses the system (brain) size in terms of the number of cell division, as powers of 2. The maximum value of is also , when using a precise value [4] of .
The unconnected monomer neurons have a weight fraction of which can be approximated with the weight fraction in an infinite system [4] as Equation (5) provides a physically measurable definition for . Since this monomer fraction does not join the tree growth scheme, the maximum “treelike” value for is 2 and the attainable DP for the gel particle is , which may be applied to the lefthand side of (4a) which serves as a correction for computing . But this is only a minor correction; for example, for , , (4a) gives and (6) yields (Cf. Table 1). The gel fraction is denoted as .

The case of and in Table 1 is the highfunctionality system of the human brain [11]. The pregel distribution or the distribution in the sol fraction in Table 1 is narrow, with a dispersion ratio or dispersity (ratio of weightaverage DP to numberaverage DP) very close to 1, even for a system of 10^{11} neurons. Equation (5) indicates that the sol fraction contains mainly unconnected monomer neurons. These free, unconnected mononeurons may form other classes of aggregate structures, such as cortex layers (horizontal sheets) and minicolumns (vertical sheets) [11], or they may join the gel in a nontree fashion, as reaches hundreds or even thousands.
3. Postgel Relations with Extensive Ring Formation
As , ring formation in the gel particle cannot be avoided [1]. Indeed when , only 1 in 10^{11} neurons remains in the sol fraction. In a human neuron with available connection sites, may reach a large value of 8000. In this case, there is a nearzero chance of finding a free neuron that remains unconnected. The large values of are a measure of excessive numbers of rings or loops. Information transfer is efficient when it takes a treelike path in such a sea of ring structures.
This theory thus provides a statistical, nongeometrical model for brain growth, neuron packing, and neuron firing. Geometrical description such as calculation with cell and branch volumes tends to ignore the fact that in chemical solutions or cellular suspensions, volumes of cellular components are not additive [11].
In the postgel range of , treelike growth reaches a maximum size determined by (6). At , monomer neuron fraction can be considered as the sol fraction . Thus the approximation or the gel fraction is a good approximation for large . Indeed the plot of versus , in the entire range of , is remarkably similar to a “brain growth” curve [12]. This kind of “growth curves” can also be depicted as developmental curves of spine densities in the cortex of rat, mouse, and guineapig [13].
The fractions of free and unconnected neurons, , being , and at and 2, respectively, are in exact agreement with those obtained by using a simple binomial distribution for the probability of connecting neighboring cortical neurons at the binomial averages [10] of 1 and 2. This agreement is not surprising since the PGFs of the originating “root” (zero generation) and subsequent generations of fold compounding all take the binomial form (Supplementary Material).
The asterisk on denotes that it is a relation that holds well beyond gelation in the presence of excessive rings. In fact, denoting the sol and gel properties by and after the gel point, the size distribution becomes heterogeneous and splits into two very narrow peaks at and around , with [4]
The FFST defines an mer as having linkages without any loop structure. The maximum tree size has a DP of neurons with , or precisely , linkages. At this maximum size, , and (2a) reduces to (4a) as a constraint set by the system size. In the entire range ; however, is also the number of linkages per neuron.
If pairs of twin synapses are uniformly distributed along a linear meric chain, it would then serve as a backbone for the entire system of connected neurons. This scheme leaves all the loopforming linked connections in the branched, nonlinear chains. However, uniform distribution appears unlikely because at , there is still a fraction of or 13.5% (by weight) of neurons remained unconnected at the maximum tree size. This is approximately the same fraction as that of human cortex [14].
In a given neuron, there may be 10000 synapses that are formed by 10000 incoming postsynapses and 10000 outgoing presynapses [11]. The total number of linkages is , and the total redundancy is defined as For , there is a redundancy of linkages. Of these linkages at the maximum tree size, there are treelike linkages paired by exactly loopforming linkages. This seems contradictory, because in a strict sense of the ringfree assumption in FFST, a linkage, or a pair of splitend linkages, is a smallest loop. Such a pair of twin synapses stabilizes a crosslink or the entire network. For , almost all 8000 redundant linkages provide not only greater redundant security, but also greater neuroplasticity to the gellike network structure of the brain. At the early explosive growth stage of synaptogenesis [9], endlinking or specially the linking may be a preferred mode for linear and treelike branching chain growth (Table 1). In this sense, each treelike linkage is considered as having a pair of twin synapses.
Redundant, loopforming linkages do not contribute to the DP of a growing chain. Thus the maximumsized gel remains at a DP of , instead of the discussed above. The maximum numberaverage DP is also , rendering the dispersity at 1. Finally the unconnected or cortical fraction of neurons should also join this growing gel, approaching a final DP of , where the treelike equations (1a), (1b), and (2a) no longer hold, but – are still valid as . This is a state in which neurons in solsol, solgel, and gelgel are all linked. This is the complete gelation scheme proposed by Flory [1, 4], which differs from that of Stockmayer and others [4].
4. Neuron Wiring in the Brain
The simplest strategy for neuron wiring is a snakelike strategy as depicted in Figure 1(a) of [11], in which a long wire (axon fiber) is connected to all 209 neurons in that figure. In a “10000 nearest neighbors” model, the long wire makes a single connection to each of the 10000 light bulbs tangentially, meaning only 1 point of contact at each bulb.
Since a snake does not bite its own tail, the snaking scheme is ringfree.
The wiring model proposed here for the connection of noncortical neurons may be called concentric “hardball” and “softball” model, obtained by invoking the above treelike definitions, numbers, and equations, plus 3 types and 3 forms of synapses, and at least 2 selfavoidance rules.
4.1. The Hardball Is Where Most of the Neuron Mass Concentrated
It is centered at the soma of a given neuron, with a radius encompassing the entire soma, plus axon and dendrite branches near the soma. All 3 cell components (S, A, and D) can provide postsynapses to an incoming presynaptic axon. All 3 types (AS, AA, and AD) of synapses may be present inside or on the surface of this hardball. The hardball is a (relatively) solid, rigid, and stationary target for the incoming straightshooting presynaptic axon tip [10].
Incoming axons can be envisioned as straight arrows, each with a long string attached to its end [10]. In the early stage of synaptogenesis, an axon arrow is shot toward the hardball target, most likely with 0 or 1 hit [10].
The random hits to the hardball can be described with a binomial distribution: where is a binomial distribution function, with random variable and parameter . But , the available target sites, is a large constant. The PGF and the mean are Note that the random variable starts with 0, and is the dummy variable [4] of the PGF. At low values, for example, , the distribution is broad, with twin peaks at and . This means the most probable outcome is 0 or only 1 hit [10]. The cases for and 2, with and , have been discussed above as the fractions of “unconnected neurons” at these averages.
The simple binomial PGF in (13) and its distribution are different from FFST size distribution in other aspects, besides the similarity in fractions and . The simple binomial distribution is for the random hits because it has nothing to do with size of aggregates or polyneurons since there is no gelation. This “zeroorone” scheme, proposed for random hits on linking of mouse cortical neurons [10], actually fits the principle of cascade formulation on treelike connection of noncortical neurons [4], which is also employed in the present treatment of ringfree assumption up to the limit of maximum tree size.
A more recent estimate is a narrow peak at for linking two neighboring cortical neurons [14]. This means the distribution is peaked around , but the peak is no longer significant. That is, with an expected 5 linkages between two cortical neurons, there is practically no chance for a “no hit.” These probabilistic arguments suggest that the connection of cortical neurons is much more orderly or structured than that in the initial bulk phase of treelike connection.
Contrary to the simple binomial probability distribution, the branchedchain size distributions in Table 1, for or more, are narrow at low , broadening with increasing values and reaching a maximum at gel point [4] when .
It can be further simplified by assuming that synapses formed onto the hardball are entirely of the twin synapses of the form, for the following reasons. Firstly, the rough, spiny, and more abundant dendrite branches provide conditions required for twin synapse formation. Secondly, even with a smooth and less active surface on the soma, the AS connections may stabilize themselves with such twin linkages. All linkages in the hardball are designed to provide a “skeleton” for a brain lasting for a lifetime of the host animal [9].
4.2. The Softball Layer
In addition to the long, thin axons, the thick layer outside the hardball, of radius from the central soma, is mainly composed of flexible and distal portions of dendrite branches. Incoming axons shot through this flexible, ribbonlike cloud have nonstationary targets for contact. Connections to this softball region are made with nontree, redundant links. A given neuron has up to 10000 contact sites for 10000 axonic presynapses from other neurons. If the synaptic density (in numbers of synapses per unit volume) is uniform in both hardball and softball regions, then where . And 2 is the number of synapses, of the form, connected to the smaller hardball. The result is simply , which is easily achievable to satisfy this large number [11] of .
5. Conclusions
For comparison with other statistical models, the FFST provides treelike parameters such as (an abstract concept of fold compounding in PGF), and (as weightaverage DP and as number of generations in cell divisions) at various crosslinking stages. Crosslinking index is for pregel linear growth; at the gel point, and at the maximum tree size. When extensive ring formation occurs, as the observed range of number of synapses per neuron.
Human cortical neurons are more homogeneous, less dense in synapses, and more ordered in their connection than the wellstudied much smaller mouse cortex [10, 11]. In the human cortex, these advantages translate into a greater neuroplasticity and much greater freedom in the information transmission process. Low synapse density in noncortical regions means an early onset of treelike branching, so that the massive, lateoccurring links can have sufficient space for high neuroplasticity. In this sense, the treelike model proposed here is behaving as a “supporting” system—much like that animal skeletons providing support, with muscles and organs filling in, for their host bodies.
The medium for neurons is well dispersed with glia cells and other mixture components. The brain growth in terms of brain volume or synapse density follows a general curve shape of gel development in . Another developmental versus time curve, in months before and after birth, surges with a rapid rate to a maximum at 8 postnatal months. Thereafter, it declines and levels off at and continues to a very mature age of 70 years [9]. These surge and decline do not affect the shape of a growth curve because the values of are in thousands.
The statistical model for a maximum tree gel leads to a physical model wherein all neuron units are connected entirely by twin synapses in the hardball region, accounting for 3 types and 3 forms of observed synapses. In the outer softball region highly redundant, flexible, accessible, and neuroplasticity are the characters of the linked synapses. In a brain, twin synapses are present in a much smaller amount than the links, but the twins are obviously more visible and easily identified [8, 11, 13].
Perhaps the simplest brain is that of C. elegans [15], in which a large sensory neuron has only one pair of “horseshoe” shaped dendrite branches. Thus, the prediction of “one linkage, twin synapses” from FFST holds even for this primitive animal. In higher animals, the observation of twin synapses [8, 15] indicates that they are relics of treelike structure.
Acknowledgments
The authors thank the following friends and colleagues for their help in preparing the paper and for their fundamental and vigorous discussions on this treelike model: Dimitris Argyropoulos, Albert Chang, George Kletecka, C. C. Wang, and Alexander K. Yan.
Supplementary Materials
The molecular weight distribution (MWD) of polycondensation of branchedchain monomers of the type RA_{f}. Here, f is the number of functional groups, or “functionality” of group A. Mathematically, this widely quoted distribution function has been treated as powerseries distribution and compound distribution that provides a simple concept; that is, singleparameter expressions of number and weightaverage “degrees of polymerization” (DP) are sufficient to generate the entire MWD for branched polymers. Moreover, using a cascade formulation involving functionals and probability generating functions (PGF), this distribution can be extended to finite systems.
References
 P. J. Flory, Principles of Polymer Chemistry, chapter 9, Cornell University Press, Ithaca, NY, USA, 1953.
 W. H. Stockmayer, “Theory of molecular size distribution and gel formation in branchedchain polymers,” The Journal of Chemical Physics, vol. 11, no. 2, pp. 45–55, 1943. View at: Google Scholar
 J. F. Yan, “A new derivation of molecular size distribution in nonlinear polymers,” Macromolecules, vol. 11, no. 4, pp. 648–649, 1978. View at: Google Scholar
 J. F. Yan, “Gelation in finite polycondensation systems,” The Journal of Chemical Physics, vol. 78, no. 11, pp. 6893–6896, 1983. View at: Google Scholar
 J. F. Yan, “Kinetics of delignification: a molecular approach,” Science, vol. 215, no. 4538, pp. 1390–1392, 1982. View at: Google Scholar
 http://news.stanford.edu/news/2008/february20/coralsr022008.html.
 P. Soba, S. Zhu, K. Emoto et al., “Drosophila sensory neurons require Dscam for dendritic selfavoidance and proper dendritic field organization,” Neuron, vol. 54, no. 3, pp. 403–416, 2007. View at: Publisher Site  Google Scholar
 “Photograph by Graham Johnson,” Science, vol. 309, p. 1990, 2005, http://www.sciencemag.org/site/feature/misc/webfeat/vis2005/show/images/slide1_large.jpg. View at: Google Scholar
 P. R. Huttenlocher, “Synaptogenesis in human cerebral cortex,” in Human Behavior and the Developing Brain, G. Dawson and K. W. Fisher, Eds., chapter 4, pp. 137–153, Guilford Press, New York, NY, USA, 1994. View at: Google Scholar
 V. Braitenberg and A. Schuz, Anatomay of the Cortex, Springer, Berlin, Germany, 1991.
 J. M. J. Murre and D. P. F. Sturdy, “The connectivity of the brain: multilevel quantitative analysis,” Biological Cybernetics, vol. 73, no. 6, pp. 529–545, 1995. View at: Publisher Site  Google Scholar
 K. L. Sakai, “Language acquisition and brain development,” Science, vol. 310, no. 5749, pp. 815–819, 2005. View at: Publisher Site  Google Scholar
 F. Karube, Y. Kubota, and Y. Kawaguchi, “Axon branching and synaptic bouton phenotypes in GABAergic nonpyramidal cell subtypes,” Journal of Neuroscience, vol. 24, no. 12, pp. 2853–2865, 2004. View at: Publisher Site  Google Scholar
 D. E. Feldman and M. Brecht, “Map plasticity in somatosensory cortex,” Science, vol. 310, no. 5749, pp. 810–815, 2005. View at: Publisher Site  Google Scholar
 http://www.sfu.ca/biology/faculty/hutter/hutterlab/research/Ce_nervous_system.html.
Copyright
Copyright © 2013 Benjamin C. Yan and Johnson F. Yan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.