The following sections first outline the theoretical aspects of the model and describe the mathematical approach taken. Then, support for the validity of the mathematical approach is provided by identifying the thermodynamic, network and immunological phenomena explained by the model.
3.1. Compartments of the System
Throughout the article, we examine the adaptive immune system in an abstract space called configuration space. Configuration here refers to the arrangement of all the constituent Abs in a system of interacting molecules. A general Ab structure is modified by the arrangement of atoms that contribute to the binding site of the Abs, creating a huge, conformational landscape. Since the formation of non-covalent bonds between interaction partners requires shape congruency, arrangement in this space determines both Ab and Ag structure and binding specificity and strength. The manifold enclosing components of the immune system in this configuration space thus define both Abs and Ags. In this space, cells and molecules of the immune system, corresponding to clones bearing or being a particular Ab, respectively, are positioned according to their potential to interact with a target and in the direction of the target. The interaction potential is the chemical potential of the Ab in body fluids. Targets are Ag molecules that drive the evolution of the system. Considering a three-dimensional Euclidean space as the configuration space of the system, target Ag shapes—epitopes—form a continuity on an imaginary, spherical canvas enclosing the system (
Figure 1).
From a physical point of view, the system is embedded in a reservoir, which serves as a heat and particle bath; this is the host organism that maintains constant temperature and pressure of the system. The system is in a stationary state in the sense that it is constantly dispersing energy while maintaining its state functions. Overall, we can regard the system as an ensemble of overlapping and hierarchically arranged thermodynamic ensembles.
We consider three compartments in our configuration space model: 1. the core of the system, surrounded by the core surface; 2. interaction space around the core, extending from the core surface to the system surface; and 3. the system environment (surroundings) beyond the system surface. We use two scales for examining and visualizing the system: the energy scale (
Figure 1A), which is the logarithmic form of the interaction probability, or the thermodynamic activity scale (
Figure 1B). The system core contains the conceptual energy source of the system that is located in the origin of the configuration space. Events within the core, therefore, can be thermodynamically unfavorable, consuming energy to maintain a continuous supply of particles in the interaction space. The core surface is spherical and is characterized by system components with average molecular free energy μ
0.
The system is composed of Ab and Ag molecules. B cells, which carry sequence information, at various stages of differentiation, determine Ab structure and define clonotypes. Soluble Abs released by B cells mediate energy transfer by binding and releasing Ags. Energy corresponds to Ag molecules, which are transferred from lower- to higher-affinity Abs, with a corresponding decrease in the system free energy. These constituents are present both in the core and in the interaction space. Interaction space is defined and confined by the totality of chemical potential vectors, each of which represents the interaction energy of a particular Ab clonotype. Radial geometrical distance in this configuration space is measured in units of chemical energy per molecule or chemical potential, while vector directions identify Ag structures and molecular shapes in the surroundings. From a network point of view, in the humoral adaptive immune system, nodes are basically Ab molecules, while links represent Ab structural similarities and, thereby, represent potential pathways of Ag transport.
3.2. Principle of Self-Organization
Self-organization is the ability of a system to arrange its constituents so as to adjust to its environment and also maintain a responsive state. Self-organizing systems tend to reach an optimal state associated with minimal interaction and dissipation [
21,
22,
23]. The immune system adjusts the quality (intensive physical property) and number of molecules (extensive physical property) in the host by the coupling of a sensor and an effector mechanism [
20]. The sensor mechanism is a B cell that survives only in the presence of signals triggered by the target molecule via an Ab displayed on the cell surface. The effector mechanism is the generation of molecules, soluble Abs, that bind both to the target molecule and to cells that remove the resulting complexes. Coupling means the balanced adjustment of the chemical potential μ of the Abs on the sensor surface and of the effector molecules against that molecule species. Biologically, it means that the sensor and effector B cells are genetically closely clonally related, expressing identical or similar Abs either on their surface or as secreted components. Relations in terms of network connectivity are determined not only by similarity but also by the similarity of interaction partners [
24]. In the case of the humoral immune system, this means that Abs that are clonally not or only distantly related may contribute to the binding and elimination of a given Ag molecule. In other words, in a subnetwork of the system dedicated for a given Ag, distinct Ab clonotypes may be linked together in the network of Ag transfer pathways. The concept of clone sizes and corresponding Ab chemical potential and Ag thermodynamic activity, with corresponding chemical potentials, can be visualized by activity maps and energy diagrams (
Figure 2). The thermodynamic activity of a particular Ab clone is related to the frequency of the cells belonging to that clone (clone size) and the serum concentration of Abs. The chemical potential of an Ab is the logarithm of its activity. Similarity is represented by the topology of the activity map, with closely related clones being neighbors. The thermodynamic activity of an Ag is related to its ability to stimulate the immune system and the corresponding chemical potential is its logarithm. Again, structural similarity is represented by location in the activity map.
Effector cells (plasma cells) secrete Abs, which bind Ags with an efficiency determined by chemical potential and Ab concentration. During the active phase of the immune response, the chemical potential of the effector molecule is raised, and more bound Ag is removed from the system. At the same time, the coupled increase in the sensor sensitivity allows the cell to survive in the presence of less Ag, and a steady state with higher chemical potential is established. The active phase of the immune response is followed by contraction, wherein memory cells, the long-lived sensors and effectors of the system, are selected and maintained. We assume that, in the resting phase of the immune response, a steady state is maintained regarding Ab and Ag concentrations. Since several different Abs with distinct chemical potentials may contribute to Ag free-energy adjustment, the system selects clones with the chemical potential and concentration required. Our argument here is that this selection process is governed by thermodynamic rules.
The cells of the adaptive immune system, a tissue that penetrates the whole organism, possesses a unique property that no other cells have; lymphocytes are capable of directing the evolution of proteins within the organism [
15]. These proteins, called Ag receptors, are generated by genetic recombination and mutation and selection events, resulting in a diversity that well exceeds the diversity of all other proteins in the organism. Depending on the size of the host, lymphocyte numbers are in the range of 10
8–10
12 cells [
25,
26]. Generation of diversity is a random molecular process, while the function of the adaptive immune system is to adjust, by cycles of division and checkpoints of selection, lymphocyte specificity and affinity. These two features correspond to the direction and length of the vector in the configuration space. A steady state, which is immunological rest in our system, is reached by the adjustment of an intensive and an extensive property. The intensive property is the chemical potential of Abs, determined by a cell’s genes and the Ab sequence and structure. The extensive property is diversity and numerosity of Abs, the coverage of the conformational space of the system surface (
Figure 1). Thus, in a broader sense, balancing is between molar free energy and system entropy.
Self-organization also maintains responsiveness in the system. In our case, it is maintained by the sensors, B cells displaying membrane-bound Abs, that are able to initiate an immune response and, thereby, reset the chemical potential of both sensors and effectors and lead to the reorganization of local hierarchy in the system. Activity maps and energy diagrams in this paper, therefore, represent a particular point in the lifetime of a highly dynamic system.
3.3. Architecture of the System
The immune system is theoretically capable of interacting with any molecular structure. This means it takes up volume in all directions of Ag space, creating a spherical core (
Figure 3). It has to allocate its limited resources for the generation of particles in the interaction space in the most energy-efficient way. Therefore, it is split up into a number of thermodynamic ensembles, distributing its resources between these. A particular ensemble is directed against a particular target epitope, adjusting its properties to attain a steady state. Because resources of the system are finite, ensembles compete for them and may overlap, cluster and form hierarchies. At any one moment in time, the distribution of the state of the ensembles is the distribution of chemical potentials and their way of sharing system energy. The configuration space of the self-organizing system is a systemic ensemble in the sense that it is a collection of coexisting grand canonical ensembles, each ensemble responding to a different antigenic component of the environment. In the vector space of the chemical potentials of the system, the location of each ensemble identifies a particular direction, representing the potential energy of the interaction with the particular part of the environment.
The probability density function of absolute thermodynamic activity λ of elements of the system, as determined by chemical potentials μ
Ab~ln(λ
Ab), can be modeled by a mixture of exponential and lognormal distributions, as proposed in the RADARS model [
27]. The model is based on the assumptions that:
Chemical potentials of Ags in the interaction space, corresponding to the molecular free energy of the interactions of the system with the surroundings, are distributed exponentially, similar to a Boltzmann distribution;
Equilibrium binding constants of Abs are distributed lognormally [
28], with a corresponding normal distribution of molar free energies of binding;
The immune system arranges and adjusts Ab chemical potentials to maintain a steady state.
The mixture of exponential and lognormal distributions, which we use here with the following parametrization, was described by Reed [
29] and Mitzenmacher [
30] as:
where:
α is the rate of the exponential distribution of chemical potential of Ag;
λ is the absolute thermodynamic activity of Abs;
μ is the chemical potential of Ag, the distance from average chemical potential μ0.
Solving the integral (see
Appendix A), the result is a power law function for λ ≤ 1, that is, for the system core,
Whereas, for λ ≥ 1, that is, in the interaction space,
These functions are the probability densities of the absolute thermodynamic activity of Abs in the system and also define a network of interaction pathways. Thermodynamic activity corresponds to the degree of nodes in the sense that it represents the probability of interactions. The exponent of network node degree distribution, −(), is thus determined by the rate of the exponential distribution of the chemical potential of Ag molecules, α, generating the interaction network.
3.4. Architecture of Interaction Space: A Hierarchical, Scale-Free Network
Interaction space spans the range of interaction energies between μ
0 and the upper limit of non-covalent binding energy. This is the space where effector Abs establish a network of interactions for Ag transport. For the description of interaction space, we can use another mathematical approach, the combination of exponentials, as suggested by Reed [
31] and Newman [
32]. Reed proposed that, when observing exponential growth after an exponentially distributed time, the process exhibits power law behavior. Here, we consider the exponential relationship between Ag chemical potential and the absolute thermodynamic activity of Abs combined with the exponential distribution of the chemical potential of Ags. In terms of molecular interactions, the former is the relationship between the free energy of the molecule and the number of different states of binding. Since a greater surface area, or more non-covalent bonds available for binding (greater buried surface area in the bound molecule), corresponds to greater free energy, absolute thermodynamic activity is a summation of the number of ways of interaction. For a given molecule, therefore, it is the number of links to other molecular structures with shared binding properties, an expression of similarity in reactivity. The combination of exponentials with the generation of networks can be illustrated by the overlay of activity maps (
Figure 4) and the mapping of the pathways of Ag flow.
Again, we assume that the chemical potential μ of Ags is distributed exponentially when the immune system is in a resting, steady state.
If the absolute thermodynamic activity λ of Abs is related to Ag chemical potential as
then the combination of these exponential functions, the distribution of Ag chemical potential and the thermodynamic activity of Abs is given by
In other words, the combination of the exponential distribution of the chemical potential of Ag and the exponential relationship between Ag chemical potential and Ab thermodynamic activity determines properties of the energy transduction network and is defined by the power law distribution of Ab thermodynamic activity. Of note, this expression is an alternative to the system equilibrium binding constant K
sys used in a previous publication [
27]. The value of λ
Ab reflects the number of pathways a single bond can evolve into a full binding surface and is, therefore, the number of links to nodes lower in hierarchy. Absolute thermodynamic activity thus corresponds to node degree, and its distribution determines network node degree distribution in the system.
Beyond relating Ab and Ag thermodynamic properties, the value of β influences how nodes are linked during renormalization. Increasing α/β is accompanied by an increasing hierarchy of sub-nodes (
Figure 4). From Equation (6), we can see that, in order to keep the value of degree exponent between 2 and 3, the scale-free network regimen, α/β should be in the range 1 < α/β < 2. Comparing Equation (6) and Equation (3), we can see that
Therefore,
and
meaning that the assumption of normally distributed energies with a variance of 1 restricts the possible values of β, and the value of α alone, via determining the distribution of Ag chemical potentials, determines the distribution of Ab thermodynamic activity and the degree distribution of the Ag transport network.
3.5. Distribution of Energy by Links in the System
The humoral adaptive immune system is a transport system for Ag molecules. The flow of Ag molecules in the interaction space of the system, resulting from a concatenated series of interactions, is an energy transduction process. As an Ag molecule is transferred from a lower-affinity Ab to a higher-affinity Ab, the free energy of the system decreases. We may examine in the system the pathways of energy transduction by assigning direction and weight to links. The direction of Ag transport is from lower to higher chemical potential, so links are always directed to nodes with a higher degree. We can assign to each link the chemical potential difference between the two nodes connected by the link. The molecular basis of this assignment is the capability of Abs to take over Ag from weaker binders (lower chemical potential). As
Figure 4 suggests, renormalization results in the joining of links with lower weight to higher-energy nodes, with hierarchy being inversely related to link weight. The number of links assigned to a chemical potential energy level is the sum of in-degrees of the nodes of that level. The average number of interchangeable links to nodes of a given chemical potential can be obtained from the relative number of direct links β and the rate variable of the distribution of Ag chemical potentials α as
This value is the degeneracy
g of the energy level. Links of identical chemical potential leading to the same node are interchangeable; they represent degeneracy in binding states. The total number of accessible states or molecular partition function Z is given by the integral of the product of Equation (9) and the exponential distribution of Ab chemical potential obtained from the distribution of thermodynamic activity in Equation (6) as
In other words, with increasing node chemical potential, the increase in activity of the nodes is limited by the number and weight of incoming links, as determined by the hierarchy of sharing interaction energies. The logarithm of the weighted probabilities of interactions in the energy shells, of the number of microstates the system ranges over, is a thermodynamic potential, entropy (S). For a system consisting of N elements, the entropy can be obtained from the partition function as
where U is internal energy, T is thermodynamic temperature, k
B is the Boltzmann constant and N is the number of particles in the system. An ideal state of the system is described by
, that is, when α/β = 1.618…, the golden ratio (19). The number of accessible states,
, as a function of Ab chemical potential in this configuration, is given by
and describes energy transport in an ideal, stationary state, supported by a golden network with an optimal balance between increasing node energy and the sharing of transport.
3.6. Thermodynamic Validation of the Model
The mathematical constructs obtained from the model so far can be corroborated by understanding the physical meaning of the variables α and β. Zheng and Wang proposed [
28] that the distribution of the equilibrium binding constant is related to the conformational space available at the given temperature and thermal fluctuation, which is, in turn, related to the heat capacity and flexibility of the molecules [
33,
34,
35]. A more flexible interaction allows for the binding of Ag to distinct Ab clonotypes, exploring greater conformation surface, while a rigid Ag epitope shows stronger preference for a single clonotype. On the other hand, distinct, flexible Ab clonotypes may bind to the same Ag, their probability-weighted chemical potentials determining the average binding energy of Ag. A higher value of α/β thus corresponds to greater flexibility and binding promiscuity, and the tail of the distribution becomes flatter; the probability of high-affinity interactions in a steady state decreases (
Figure 5).
In analogy to the heat capacity ratio, which relates heat capacity at constant pressure to heat capacity at constant volume,
The ratio of α/β defines a binding flexibility ratio in configuration space where
This relates binding flexibility to chemical potential μ, kept constant to flexibility with number of system elements N kept constant, in the following sense: The partial differential expression in the numerator describes the extent of conformational surface S explored by Ag molecules when these molecules are provided by the surroundings so as to keep chemical potential μ constant. The expression in the denominator describes the potential to explore the conformational surface due to thermal fluctuation T; when no extra molecules are supplied, N is invariant. The immune system is, indeed, capable of fine thermodynamic tuning of the interactions of its Abs by selecting the structure with appropriate flexibility by isotype switching and affinity maturation and by targeting selected Ag epitopes. The regulation of these Ab properties for each and every Ag is the essence of controlling Ag concentrations.
If we regard blood plasma as an open, single-phase, multicomponent system, we can use a thermodynamic potential for the description of energies related to Ab–Ag interactions. Following the nomenclature used by Emmerich [
36] for a thermodynamic potential expressed for extensive entropy,
Thermodynamic potential Λ summarizes the essence of the adaptive thermodynamic system: create a large conformational landscape (conformational entropy S) and generate N elements with chemical potential μ expressed as the ability to cover conformation space,
Equation (16), in our model, states that the addition of molecules to the system changes the conformational entropy according to the molecular free energy μ per unit temperature T. We can replace the summation in Equation (15) with the probability-weighted sum of chemical potentials, which is the average chemical potential
, and combine with Equation (11) in the form:
Thus, the thermodynamic potential of Abs depends on their number N, the number of interactions they can engage in, summarized in Z, and their average chemical potential. The configuration space model introduced here arranges constituents of the system based on their “power” to change other constituents’ chemical potential. This power is reflected in the topology and network linkage of constituents. In a stationary state, the distribution, and, thereby, the average chemical potential of all Ag molecules, is kept constant, adjusting binding energy that tickles but does not activate sensors, memory or naïve B cells. This is achieved by setting a chemical potential that immunologically suits the host: low chemical potential for harmless self molecules, high chemical potential for dangerous non-self molecules. Concentrations of targets are then adjusted accordingly, meaning that Ags with high chemical potential are eliminated with higher efficiency. The process leads to a wide range of chemical potentials of Ags and Abs in the system. A theoretical, system-wide steady state is reached when Ab free energies are distributed exponentially. In the ideal configuration of Equation (12), the partition function Z is equal to one; thus, Equation (17) simplifies to the product of element number N and the average chemical potential.
Self-organization of the system refers to how many Ab molecules are generated and how Ab chemical potentials are distributed. In a thermodynamic steady state, the system minimizes free energy and maximizes entropy; this state can be described by statistical mechanics functions. The adaptive, self-organizing biological ensemble is a co-existing collection of large numbers of copies of states, with a regulated total extensive free energy G and a composition that may fluctuate with the distribution of chemical potentials adjusted via network formation (
Table 1). In immunological terms, in order to reach a steady state, the chemical potential of Abs in the system has to be tuned with regards to all potential binding partners, not just the Ag that triggered a response. The pathway to reaching the steady state is the selection of Abs (cell clones) that fit into a previously organized network. Steady state is reached when the system, again, possesses an extensive free energy G that continuously drives the binding of Abs to Ags.
The biological, self-organizing system diversifies its tool of energy transfer, Ab molecules in our case, adjusting the distribution of chemical potentials so as to maintain the system against forces of change in the environment. By deploying mechanisms to sense the free energy of constituents in the environment, the system can adjust, readjust and evolve with the environment.
3.7. Immunological Validation of the Model
The humoral immune system must regulate, over a very wide range, the concentration of a vast number of molecules that are found in and constitute an organism. How it is technically possible to eliminate certain molecules while leaving others unharmed has led to the long-standing question of tuning recognition specificity and affinity [
37], breadth and depth [
38]. The solution provided by an adaptive immune system apparently requires a scale of numerosity, diversity and affinity comparable to that found in the organism—or, more accurately, the supraorganism [
39]—itself. A mechanistic approach to the feedback and tuning procedure, based on saturation of Abs and Ags, was proposed in a quantitative model of B-cell development and Ag removal [
8,
40]. The model treated distinct molecules independently, neglecting the effect of cross-reactivity but, nevertheless, providing a general framework for understanding the system. On the level of individual cells, the level of engagement of Ag receptors determines the cell’s fate: programmed death, survival, proliferation or differentiation. An initial repertoire of naïve cells, which are constantly generated, displays receptors produced by the random rearrangement of gene segments and subsequent selection steps. The immunologically controlled encounters of these cells with the Ags of the inner and outer environment lead to the expansion and differentiation of clones according to the immunological ranking of the target Ag. These processes are happening continuously, manifesting as the adaptation of the system to the antigenic environment or, in other words, as self-organization. As a result, a repertoire of memory B cells and long-lived plasma cells is produced, which reflects past adaptation and serves as the basis of further evolution in response to the environment [
20,
41].
The present model identifies variables that characterize network formation and define the system on a thermodynamic basis. Both Abs and Ags are cross-reactive, being able to bind multiple partners. The average binding flexibility or cross-reactivity on the system level manifests as the slope of the distribution curves when a logarithmic scale is used (
Figure 6). Changing the slope of either influences the other via self-organization, as described above. Complexes of Abs and Ags are constantly removed from the system via cells with Ab Fc receptors. Changes in the composition and availability of Ags trigger the activation of sensor cells, leading to readjustment of the properties of the repertoire (
Figure 6). It is tempting to speculate that the local hierarchy of binding has important immunological consequences; minimal hierarchy would result in a “powerful” Ab response, in the sense of high-affinity binding and appropriately switched Ig isotypes, such as an allergic IgE response. With the application of novel tools for the deeper analysis of serological reactions, such questions can be experimentally addressed [
18].
Immunoassays probe the system of Abs via measuring Ab binding to selected Ags. These assays are also called serological assays because they characterize serum Ab reactivity against medically relevant Ag targets. An immunoassay that uses Ag titration, a gradient of Ag concentrations, can be interpreted as the measurement of the changes of chemical potential in the interaction space and can be modeled by the Richards differential equation [
20]. In the sigmoid Richards curve, the point of inflection corresponds to the sum of probability-weighted chemical potentials and is, therefore, a measure of Ab affinity. Meanwhile, the asymmetry parameter of the Richards function is related to α/β in the system, which corresponds to the limiting activity coefficient of Ags. Thus, the variable that describes the hierarchy and network organization of Abs in the system appears in a biochemical measurement as a thermodynamic variable.
A theoretically important message of the model is the continuity of self. The interaction space belongs to the system; it is self. However, it also incorporates elements of the surroundings in a regulated way. Ag-binding energies in the interaction space range from system average to very far from average, and elements are distributed according to an exponential distribution. Therefore, the frequency or mole fraction contribution to the system also covers a very wide range (
Figure 7). All binding events belong to the system, becoming incorporated into the architecture as imprints in the network. This is consistent with the liquid hypothesis of self, which states that immunological identity is continuous and dynamic [
42]. Even though we have no exact physical models describing the functioning of the immune system, we have more and more experimental data on the structural and molecular mechanisms of Ab function and networks of Abs based on sequencing. It is, thus, possible to relate our model to these observations. The total immunoglobulin concentration ([Ab]
t) in adult human plasma is in the high micromolar range (~10
−4 M). We can express Ab quantities as a fraction of the system by using this reference as [Ab]/[Ab]
t (
Figure 7).
Assuming that the system is in a steady state when all its components are present at their equilibrium dissociation concentrations K
D, their distribution as a function of chemical potential is determined by the standard free energy of their interactions. Gibbs free energy has two components:
Namely, an entropic (S) and an enthalpic (H) component. We can relate the free energy of binding to these components using the relationship between the distributions of chemical potentials and interactions, combining Equations (10) and (18) as
The difference between using only the enthalpic component or both components is illustrated by the solid and dashed lines, respectively, of
Figure 7. The concentration of an Ab clone with a particular binding energy in the system in steady state is higher than that of the same Ab alone in equilibrium with its single binding partner. Simply stated, because interaction energies are “spread out” over the network, equilibrium concentrations in a system are higher than in a solution of two components, such as a monoclonal Ab and a single target Ag.
3.8. Network Validation of the Model
A theoretical network of Abs can be generated on the basis of sequence and structure similarity [
27], which corresponds to experimentally determined, Ag-binding correlation networks [
43]. On the molecular level, these networks allow the transfer of Ags from an Ab node to another Ab node in the direction of increasing affinity. The configuration space model of this article organizes such links into a hierarchy of chemical potential and renders structural similarity into topological relationships. Increasing chemical potential in the configuration space means shifting in the hierarchy and increasing node degree, which is related to the probability of receiving an Ag molecule from a node lower in hierarchy. All nodes high in the hierarchy are “supported” by a large number of nodes with lower chemical potential. This arrangement causes the spreading of hubs over the interaction space and manifests as the repulsion of hubs and results in a disassortative network.
The power law degree distribution and scale-independent property of the network were previously derived from a lognormal distribution of equilibrium binding constants of random interactions in combination with an exponential distribution of the standard deviation of binding energy [
27]. Here, a power law distribution is also obtained from a Boltzmannian distribution of molar free energy of Ags and the relationship between the absolute thermodynamic activity of Abs and the chemical potential of Ag. Cross-reactivity is the result of partial interactions with lower-than-maximal energy. Higher-energy interactions allow more cross-reactivity, since a greater molecular surface area is available for binding. This cross-reactivity appears in the renormalization strategy (
Figure 4), where nodes with lower degree are joined to nodes with a higher degree. This renormalization is similar to the box-counting method of determining fractal dimension [
44] in the sense that grids with increasing side lengths are applied to reveal network architecture. Disassortativity and fractality have been recognized as features of biological networks [
45,
46].
From a network perspective, the ratio of our variables α/β corresponds to the ratio of fractal dimension d
B and degree exponent of boxes dk, as described by Song et al. [
45,
47]. The network degree distribution exponent γ can be calculated from these as
This is in agreement with our findings that α/β is the rate parameter of the exponential distribution of chemical potential, and α/β+1 is the exponent of the power law distribution of thermodynamic activity and network degree distribution (Equation (6)) (also see
Appendix B). The evolution of the network in time is reflected in configuration space, in as much as the nodes with greater chemical potential that develop in the time course of immunological reactions are located more distantly from the system surface. Therefore, these indices can be interpreted as factors of renormalization in energy levels.
According to Caetano-Anollés et al. [
48], scale-free networks can follow a trajectory in network morphospace from homogenous, non-modular towards heterogenous, modular network organization. This trajectory can correspond to a mesh-like structure of random binding events at the core surface of the system, evolving into a highly hierarchical and modular network of specific and high-affinity interactions. This corresponds to a star-like network distribution of energy when α/β approaches 1 and multiple transfers of energy when α/β increases (
Figure 5). Whereas the probability density of binding energy is associated with a Boltzmannian distribution, scale-free networks are known for the power law distribution of network node degree. Our configuration space model of a system suggests that these two phenomena are two sides of the same coin; an identical, complementary cumulative distribution results in the exponential probability density function of chemical potential and power law probability function of node degrees of elements of the system (see
Appendix B).
3.9. Interpretation of the Model as a Complex System
In their maximum-entropy model for Ab diversity, Mora et al. associated Ab sequences with effective energies as if the sequences represented a particular state in a system in equilibrium [
16]. In this article, it is suggested that Ab structures, which are defined by their amino acid sequences, are, indeed, distinct states in the configuration space of the system. The molecular free energy of these structures is determined by the composition of the system, which is itself regulated by the system.
Following the active, expansive stage of the immune response, a contraction stage establishes a balance between the thermodynamic activity of Abs and of Ags by selecting only a subset of the cells that evolved in the active stage. This process leads to a steady state, whereby secreted Abs are continuously removed along with bound Ags, keeping both free Ab and free Ag concentrations at the immunologically adjusted values. A two-way feedback mechanism of sensor and effector cells adjusts chemical potential so that cells are poised between under- and overactivation (see [
8,
20,
40] for details). The effector cell secretes Abs capable of reducing target Ag concentration; the extent of reduction is determined by the chemical potential of the Ab. The sensor cell receives an amount of energy in the form of survival signals by the target Ag—this is called tickling in immunological jargon. This signal is adjusted by the effector cells, secreting soluble Abs. From an immunological point of view, this means the parallel generation of memory B cells (sensors) and plasma cells (effectors) in germinal centers with closely related if not identical binding properties [
49,
50,
51]. Perturbations at any point in the interaction space trigger the rearrangement of this hierarchy and network. In this sense, the system is poised at criticality, a phenomenon suggested to be present in all biological systems [
52]. Critical events represent the reorganization of hierarchy in a system. These events are the coalescence of sub-ensembles or subnetworks. The size of these events often shows power law distribution over time (strength of earthquakes, size of forest fires, avalanches of sandpiles) [
53]. The immune system presumably reorganizes itself constantly, adapting to the antigenic environment, via such events. Occasionally, triggered by infections, massive reorganization is necessary, which may correspond to a huge critical event.
A variety of phenomena has been shown to follow power law, including natural events and systems and artificial, man-made systems. In these systems, the power law applies only for values greater than minimum value x
min. The distributions described here suggest that the minimum value of x corresponds to activity with the average energy or reference chemical potential μ
0 of the system. By our definition, this reference value is zero, and the corresponding activity is one.
It also means that, in our directed network, all isolated nodes possess an in-degree of 1 (see
Figure 5). Nodes with a higher degree connect to form the networks, which represent the organization of the interactions of the system with its environment.
Power law relationships can be broadly assigned into two main categories: distributions of frequencies and of magnitudes [
54]. Complementary cumulative distributions (cCDF) of frequencies, rank–frequency plots, Zipf’s law and the Pareto distributions are examples of the first type [
30,
55]. Protein interaction and metabolic networks are examples of the second type [
56,
57].
Magnitudes represent an intensive physical property of the adaptive system. The distribution of an intensive physical property describes a hierarchy and often represents a hierarchical network that organizes the system. It follows that these magnitudes are related to the network degree of the entities. The configuration space model presented here suggests that these magnitudes are absolute thermodynamic activities and are related to molar free energy or chemical potential, an intensive physical property.
Frequency is an extensive property and corresponds to the number of links (interactions) of a given energy in our model. If the energy of interactions is distributed exponentially (Equations (4) and (12)), then the logarithm of probability (Equation (9)) decreases linearly (
Figure 6,
Figure 7 and
Figure 8). The corresponding linear growth of entropy (for cumulative distribution) is characteristic of criticality in thermodynamic systems [
16,
58]. In unchanging, simple surroundings, the system can strongly adapt (
Figure 8, vertical arrow) and distribute its energies over a wider range of chemical potentials. On the other hand, changing, diverse surroundings prohibit effective adaptation (
Figure 8, horizontal arrow), allowing large numbers of weak interactions. The system can develop and grow, maintaining its organization by parallel growth of its intensive and extensive properties, chemical potentials and the number of interactions as a function of those potentials, respectively.
Power law distributions, in general, characterize critical points and phase transition. In this respect, self-organization aims to maintain the state of phase transition. In an abstract sense, transition is between the two phases of organization of matter in the system and in the surroundings. Self-organization thus counteracts the effects of the surrounding environment by generating an interaction space and, therein, sustaining phase transition, maintaining steady-state reactions towards all environmental components that would otherwise disintegrate the system.