#### Abstract

Boolean networks are used for modeling and analysis of complex systems of interacting entities. Classical Boolean networks are binary and they are relevant for modeling systems with complex switch-like causal interactions. More descriptive power can be provided by the introduction of gradation in this model. If this is accomplished by using conventional fuzzy logics, the generalized model cannot secure the Boolean frame. Consequently, the validity of the model’s dynamics is not secured. The aim of this paper is to present the Boolean consistent generalization of Boolean networks, interpolative Boolean networks. The generalization is based on interpolative Boolean algebra, the -valued realization of Boolean algebra. The proposed model is adaptive with respect to the nature of input variables and it offers greater descriptive power as compared with traditional models. For illustrative purposes, IBN is compared to the models based on existing real-valued approaches. Due to the complexity of the most systems to be analyzed and the characteristics of interpolative Boolean algebra, the software support is developed to provide graphical and numerical tools for complex system modeling and analysis.

#### 1. Introduction

Approaches based on logic are used for modeling systems with complex causal interactions between the components [1–3]. Obtained models usually have the form of directed graphs, with a logical state variable and a logical function assigned to each node. The state of a node evolves in continuous or discrete time, in accordance with the value realization of the logical function. In comparison to traditional models of dynamical systems, logical models are less detailed and requires less information about the process of interest. Determination of an appropriate level of detail is considered as the major challenge in modeling. If the model is very detailed, it might be extremely difficult to gain insight; thus, the essential nature of the system may be obscured. On the other hand, if too little detail is included, the model may not be a reliable representation of reality. The appropriate level of detail is determined with respect to the purpose of the model and the amount of information about the process of interest. The logical modeling is aimed at capturing the nonlinear dynamics of possibly large-scale systems. The modeling of such systems is typically hampered by the lack of information. According to the principle of incompatibility [4], as the complexity of a system increases our ability to obtain both relevant and precise model of the system diminishes, and, eventually, these two model characteristics become incompatible. In the logical approaches, the precision suffers for the relevance. This enables modeling of systems that are too complex for traditional approaches.

Boolean networks (BNs) are the common logical model for modeling complex systems based on binary logic. They may be seen as discrete dynamical systems where each node logical function works dynamically to compute the node state at the next time step depending on the input values coming from other nodes [5, 6]. Therefore, the research in this field mainly addresses the structural and dynamical aspect of BN. The BN structure can be specified by an expert [7], identified from the data [8] or randomly constructed [5, 9]. On the other hand, the validity of the dynamics is secured by the laws of Boolean algebra (BA). According to the principle of truth functionality, the values of logical functions in such models are calculated by the values of their arguments in parallel or serial manner [10]. Classical models are widely used in biology and medicine [11–14], geosciences [15–18], and social sciences [19–22].

BNs are relevant for modeling systems which exhibit threshold behavior. If more details are available, BNs can be generalized by the introduction of gradation. In this case, logical variables and functions are treated as -valued. In [9, 23–28], BNs are generalized using conventional, truth functional fuzzy logics. Since such fuzzy logics are not a BA [29, 30], the generalized models cannot secure the Boolean frame. Consequently, the validity of the dynamics of the model cannot be secured.

In this paper, we proposed the Boolean consistent -valued generalization of BN. This approach, called the interpolative Boolean networks (IBN), is based on interpolative Boolean algebra (IBA), the -valued realization of BA [31]. Unlike conventional fuzzy logics which are based on the truth functionality principle, IBA is based on the structural functionality principle to secure Boolean consistent generalization. According to this principle, logical functions have vector nature, and thus, the Boolean consistent calculations of values can be accomplished by the immanent structure vectors. This enables IBA, and consequently the IBN, to be used for modeling and analysis of complex dynamical systems in a Boolean consistent manner. Further, IBN is the general -valued approach, while Boolecube, Hillcube, and similar approaches may be seen as the special case of IBN. The validity of the dynamics of IBN is illustrated on two examples from the literature using a software support developed for this purpose.

The paper is organized as follows. In Section 2, we describe the classical binary BNs and the existing real-valued realizations of BNs. Section 3 provides the basic concepts of interpolative Boolean algebra. In Section 4, the proposed generalization of BN is introduced. Furthermore, the software for complex system modeling and analysis, implementing IBN, is presented. In Section 5, the proposed approach is illustrated on two examples. The obtained results are compared with both conventional fuzzy approach and the Boolecube/Hillcube approach. The main conclusions and directions for future work are given in Section 6.

#### 2. Boolean Networks

Boolean networks are the least detailed mathematical models of causally interacting entities. The model has the form of a directed graph, whose nodes represent the entities, and the edges represent the interactions between the entities:where is a set of nodes and is a set of direct edges [5].

##### 2.1. Binary Boolean Networks

Dynamics of BN is introduced by assigning a Boolean state variable and a Boolean function to each node. It is assumed that all nodes change state at each time step, with a unit time-delay. The dynamics of a BN reflects the topology of the directed graph. The state of node takes the value of its Boolean function whose arguments are the states of the node’s predecessors in the previous time step: where and are the state and the Boolean function of node , are the states of node predecessors in the previous time step, and is the number of predecessors of the node .

Dynamics of BN is defined by equations which update the state of all nodes in a system:where , , is the state of a Boolean network at time , and is the number of nodes.

Transition functions in a classic BN, as any binary logical function, follow the principle of truth functionality [32, 33]: the value of a complex formula is uniquely determined by the truth values of its subformulas [34]. In other words, the value of a logical function is calculated using the values of functions’ arguments. Since the algebra of binary BNs is Boolean, the validity of the BN dynamics is secured. The global dynamics is characterized by the existence of one or more point or cycle attractors. This is a state or a set of states where the dynamic of the network stabilizes. For each attractor there is a basin of attraction which contains states that converge to the attractor. The dynamics of a BN depends on which basin of attraction the initial states are located. The number and size of attractors, and the corresponding basins of attraction, are the state space properties which reveal the global dynamics of BNs.

Classical Boolean networks are adequate for modeling systems with complex switch-like causal interactions. In such systems, cause and effect relationships between the components are governed by thresholds [35]. Below (or above) a certain level, a component has a little or no influence on the behavior of others, while above (or below) this level the influence saturates rapidly to a constant level. Such a switch-like causal relationships can be represented using sigmoid functions. Binary representation in the classical BNs is an idealization of such relationships.

##### 2.2. Real-Valued Realization of Boolean Networks

When more detailed description of a system is required and/or possible, the Boolean networks can be generalized by the introduction of gradation. In this case, logical state variables and logical functions are realized as real-valued. This provides the drastic increases of the descriptive power of the model [36–38]. The behavior of a system can be quantified with much more precision, and thus, deeper insight into the system’s dynamics is possible. In [9, 23–28], Gödel, product, and Lukasiewicz conventional fuzzy logics are proposed for generalization of BNs. However, the validity of the dynamics of the proposed generalized models cannot be secured, since the Boolean laws of excluded middle and contradiction are not followed in fuzzy logic [29]. This implies that the behavior of complex systems may not be accurately predicted by such models.

The Boolecube/Hillcube model is able to predict and explain quantitative experimental data [3]. This approach is the generalization of BNs based on multilinear interpolation of Boolean functions. Boolean functions are replaced by their real-valued realizations called the Boolecubes or the Hillcubes. The Boolecube is defined by the following expression:where , , and is a Boolean function that indicates the inclusion of the variables/functions in the Boolecube.

The Hillcube approach is the generalization of the Boolecube approach since it uses the sigmoid function for mapping variables. This function is called the Hill function and it is defined by the following expression: where is the Hill coefficient which defines the steepness of the function and is the threshold parameter. Therefore, the Hillcube is defined by the following expression: where , , is a Boolean function, and is the Hill function.

#### 3. Interpolative Boolean Algebra

Interpolative Boolean algebra (IBA) is the Boolean consistent -valued generalization of finite (atomic) Boolean algebra (BA) [30, 39]. Formally, IBA is based on the principle of structural functionality [39–41]. It implies the necessity to distinctly detach the structure of the logical function from its value realization. Therefore, IBA is two-leveled algebra consisting of symbolic and valued level. The structural transformation on a symbolic level enables us to evaluate any logical function in the Boolean consistent way. Once the transformations have been conducted, the value level is introduced. The values do not have an influence on the preservation of Boolean laws. Therefore, IBA is used as a natural frame for Boolean consistent fuzzy logic [42], logical aggregation [40], similarity measure [43], and computational cognition [44].

##### 3.1. The Principle of Structural Functionality

Unlike the truth functional fuzzy logic and many-valued algebras, IBA is based on the principle of structural functionality [39–41]. This principle implies that the immanent structure of any logical function can be directly derived from the structures of the function’s arguments. In other words, any element of finite Boolean algebra uniquely corresponds to its structure (content, relation of inclusion). Relying on the principle of structural functionality is essentially different from relying on the principle of truth functionality—the structure of the observed function is at the forefront rather than the values of its arguments [45].

The main motive of introducing the principle of structural functionality is that the principle of truth functionality is binary in its essence, that is, the truth value of an expression may be true or false. Therefore, the truth functionality principle is valid in the sense that it preserves all Boolean axioms and theorems for a two-valued realization of Boolean algebra, while it is not valid in many-valued case and fuzzy logic [29]. On the other hand, according to the principle of structure functionality, the structure of the logical expression and the intensity of variables should both be taken into account in -valued case. The structure of an element determines which atomic elements of IBA are included in it and/or which are not included in it. An atomic element of IBA is the simplest elements of the Boolean algebra domain of attributes BA in the sense that they do not include in themselves anything except for a trivial Boolean constant 0. Atomic element expresses the intensities of variable within the observed structure. Detailed mathematical formalization of IBA structural functionality is given further in the text.

In [40], it is shown that the truth functional principle is a value realization of the structure functional principle which is valid (in the sense that it preserves all Boolean axioms and theorems) only for a two-valued realization of Boolean algebra. On the other hand, the structure functional principle is valid in general case [40]. Therefore, the principle of truth functionality is a special case of the principle of structural functionality from the aspect of the Boolean laws.

##### 3.2. IBA Symbolic Level

From the mathematical point of view, any logical function on IBA symbolic level is represented in a vector form, structure vector. The value of any logical function can be defined as a scalar product of two vectors: where is the structure vector of and is atomic vector. The structure vector is binary and it represents the shape/structure of the logical function, while the atomic vector expresses the intensities of variable. In classical binary case, the intensity of variables is , so it is clearly a special case of IBA approach.

The structure vector is a binary vector defined by the following expression [40]: where are the structure function of logical function , is an atom of , and is the power set of a set of free variables .

The structure vector contains information about which atoms are relevant/included in a logical function. Such information is given by the structural function , which is defined by the following expression [40]: An atom in IBA is defined by the following expression [40]:where . For instance, atomic elements of BA generated by two attributes are , , , and .

In accordance with the structural functionality principle, any logical function can be uniquely represented by the following disjunctive canonical form [39]: In order to obtain the real value realization of any logical function, IBA uses generalized Boolean polynomials (GBPs). GBP is a sum of the relevant atomic Boolean polynomials:where and is a generalized product.

##### 3.3. IBA Structural Transformation and Valued Level

IBA is technically based on the mapping procedure of logical function to GBP. Any Boolean function can be mapped into corresponding GBP using the following rules [40]:The generalized product is any function that maps and satisfies all four axioms of -norms (commutativity, associativity, monotonicity, and boundary condition) and the condition which ensures that the values of atomic Boolean polynomials are nonnegative [39, 40]. Therefore, the generalized product for may be realized as any -norm that produce the value from the following interval: Three cases may be distinguished as particularly important realizations of the generalized product: minimum function, ordinary product, and Lukasiewicz operator.

#### 4. Interpolative Boolean Networks

In this section, we present the Boolean consistent -valued generalization of Boolean networks, interpolative Boolean network. In the proposed model, states of nodes take value from the unit interval, while the Boolean functions are mapped into generalized Boolean polynomials according to (13). Such polynomials are the -valued realizations of Boolean functions with respect to the structure functionality principle. The proposed generalization provides a drastic increase of the descriptive power of BNs. It allows the model to be used for gaining quantitative insight into the behavior of complex systems. For the purpose of this paper, we assumed that interpolative Boolean networks evolve in discrete time, although continuous time models can be generalized in the same way. Formally, in the IBN, the state of node takes the value of its generalized Boolean polynomial whose arguments are the states of node predecessors in the previous time step: where and are the state and the Boolean function of node , are the states of node predecessors in the previous time step, is the structure function of , is the generalized product, and is the number of predecessors of the node .

The state of node can be defined as the scalar product of two vectors: where is the structure vector of and is a vector of atomic Boolean polynomials of at time .

Dynamics of interpolative Boolean network is defined by the following equations:where , , is the state of IBN at time and is the number of nodes.

The interpolative Boolean network evolves according to the principle of structure functionality. This secures the Boolean frame and the logical validity of the model’s dynamics. The global dynamics is characterized by the existence of point, cycle, and/or chaotic attractors, with the corresponding basins of attraction. The behavior of IBN can be fundamentally different depending on which basin of attraction the initial states are located. If the initial state of each node is two-valued, IBN reduces to the classical model; that is, the structure functionality principle reduces to the truth functionality principle. IBN can be used for modeling complex systems which may or may not exhibit the threshold behavior.

The behavior of IBN depends on the way how generalized product operators are treated. The adequate generalized product operator for the two variables is determined with respect to the nature of the variables. Two logical variables and are of the same nature if they represent the same quality. If this is the case, the adequate generalized product operator is minimum function; that is, . If two variables are of the same nature, but negatively correlated, the Lukasiewicz operator is used; thus, . On the other hand, if two variables are not correlated, that is, these variables are of the different nature, the proper operator is ordinary product, and thus, .

In IBN, different operators for generalized product can be used for any pair of variables. In this way, the nature of the variables can be fully taken into consideration. Comparing to the Boolecube/Hillcube approach, in which it is assumed that all variables are of the different nature, IBN offers more descriptive power. The Boolecube/Hillcube implicitly take into account the structure of elements; its Boolean indicator functions have the role of IBN structure vector, but the only ordinary product is applied as the -norm. In other words, if each generalized product operator in IBN is treated as ordinary product, the proposed approach reduces to the Boolecube/Hillcube approach (see (4), (6), and (15)). Further, IBN utilize comprehensible and direct IBA transformation procedure from logical function to GBP given in (13), which provides an easier interpretation of the model.

##### 4.1. Software Support

For the task of modeling and analysis of complex dynamical systems, the software tool implementing interpolative Boolean networks is developed. The architecture of the software is shown in Figure 1. The software is based on jFuzzyIBATranslator, the very first software realization of IBA, which can also process conventional fuzzy expressions [45]. jFuzzyIBATranslator is an application written in the Java programming language using NetBeans 6.9 development environment. The task of jFuzzyIBATranslator is to transform logical functions into corresponding generalized Boolean polynomials according to the transformation rules given in (13). Logical functions are validated (syntax validation, nested parentheses, etc.) and then forwarded to jFuzzyIBATranslator. The minimal GBPs are obtained after application of the distributivity law and the rule of idempotency. The developed software tool is covered by jUnit tests, a popular testing framework for Java [46], to verify software correctness.

The procedure for complex systems modeling and analysis using the software tool is the following. The input to the software is the number of nodes of IBN, where each node represents a system’s component, and an initial state and a logical function for each node. The states take initial values from the unit interval. Logical functions are validated and transformed into GBPs. The nodes evolve in discrete time, governed by the GBPs. Since the initial states and the GBPs for each node are obtained, IBN can be simulated. The software tool predicts the future behavior of a complex system. After a number of iterations, an attractor is reached, while the dynamics of each node are presented in the graphical form. The software can be used for analysis of classical Boolean networks. In this case, it is sufficient to enter the binary initial states for each node.

#### 5. Illustrative Examples

To illustrate the proposed approach, two classical BN models are taken from the literature as examples and generalized using IBA. The obtained IBN models are then simulated with the software tool presented in Section 4.1. Further, the obtained IBN results are compared to the results of conventional fuzzy BN and Boolecube approach.

##### 5.1. Comparing the IBN Approach to the Conventional Fuzzy BN Approach

To compare the proposed IBN approach with the conventional fuzzy BN approach, the example of WNT5A gene regulatory network (Figure 2) is used. The WNT5A gene regulatory network is related to melanoma and thoroughly studied as BN in [47, 48]. First, the basic model is generalized to obtain both fuzzy BN and IBN models. Further, these real-valued networks are simulated using two -norms: minimum and ordinary product. Finally, the obtained results are compared and discussed.

The basic WNT5A gene regulatory network dynamics is described with the following logical expressions [48]: where , .

The corresponding IBN model, presented in (A.1) of Appendix A, is obtained by mapping Boolean functions from (18) into the generalized Boolean polynomials. To translate the obtained GBPs, two operators of generalized product are used and their results are compared. In the first case scenario, a minimum function is used as an operator of generalized product; that is, , where . In the second case scenario, generalized product is replaced with an ordinary product, , where . For both of selected operators, the proposed IBN model is reduced to the classical fuzzy network model and results are compared.

*Case 1. *In the first case scenario, we investigate the differences between IBN and classical fuzzy BN approach in the case of the Gödel (minimum) -norm. The following equations (19) and (20) present, respectively, IBN and conventional fuzzy BN generalizations of the classical BN given in (18), obtained with minimum operator.

IBN model with minimum operator () for WNT5A gene regulation is defined with the following set of equations: where and .

The conventional fuzzy BN model with minimum operator for WNT5A gene regulation is defined with the following set of equations: Both IBN and fuzzy BN models are simulated for the same initial conditions. The obtained dynamics are presented in Figures 3 and 4.

The final results of the experiment are presented in Table 1.

Although the dynamics of the state variables are similar in two approaches, state variables , , and reach significantly different point attractors in IBN model as compared with fuzzy BN model (values are bold in Table 1). These differences are due to inconsistency (in a Boolean sense) of the conventional fuzzy logic approach.

*Case 2. *In the second case scenario, we investigate the differences between IBN and classical fuzzy BN approach in the case of the ordinary product -norm. The following equations (21) and (22) present, respectively, IBN and conventional fuzzy BN generalizations of the classical BN given in (18), obtained with product operator.

The IBN model with product operator for WNT5A gene regulation is defined by the following set of equations: where and .

On the other hand, the conventional fuzzy BN model with the same product operator for WNT5A gene regulation is defined as Notice that in this case the only difference between the two models is in the second equation. The Boolean function for state variable from (18) is too complex to be generalized in the Boolean consistent way with the conventional fuzzy approach.

Models are simulated for the same initial conditions, and obtained dynamics are presented in Figures 5 and 6.

The final results of the experiment are presented in Table 2.

In this case scenario, the dynamics of the state variables of the two models are much more similar than in the previous case. There is only one state that significantly differs in the value of the point attractor reached, and it is . It is a direct consequence of the difference between IBN and fuzzy BN model for state space equations. However, this difference is important since it indicates inconsistency in the classical fuzzy BNs. It is also important to mention that the presented WNT5A gene regulation models based on product operator need more time to reach its point attractors comparing to the previous models based on minimum operator. This is due to the characteristics of the operator itself.

##### 5.2. Comparing the IBN Approach to the Boolecube Approach

To compare the proposed IBN approach with the Boolecube approach, the example of modeling gene interaction is used. The starting point is the classical BN model used in [7] to describe interactions between the* Drosophila* segment polarity genes. The following equations specify its dynamic:where * *, is the network state at time , and .

The IBN model, presented in (B.1) of Appendix B, is obtained by mapping the Boolean functions into generalized Boolean polynomials. Operators of generalized product in this model are treated in the following way. Since the variables and , and , and are of the same nature, the adequate generalized product operator for these pairs of variables is minimum function. Thus,The remaining variables are of different nature and, therefore, the adequate operator of generalized product is the ordinary product. The final IBN model is presented in (B.2) of Appendix B.

On the other hand, in the Boolecube/Hillcube approach, each generalized product in (B.1) is replaced with an ordinary product. The final Boolecube/Hillcube model is presented in (B.3).

Both IBN and Boolecube/Hillcube models are simulated using the following initial conditions , . The obtained models’ dynamics are presented in Figures 7 and 8.

The final results of the experiment are presented in Table 3.

It is easy to see differences in IBN and Boolecube models dynamics and their point attractors reached. This difference is a direct consequence of the nature of the model variables that is taken into account when selecting an appropriate operator in the IBN approach, unlike the Boolecube/Hillcube approach. For this reason, the IBN offers superior descriptive power and it is adaptive with respect to nature of inputs. In fact, the IBN approach can be considered as the generalization of the Boolecube/Hillcube approach since it can be reduced to it when only the product -norm is used as an operator of the generalized product.

#### 6. Conclusion

Logical approaches are widely used for modeling and analysis of complex dynamical systems. Classical models are within the Boolean frame, which secures the validity of their dynamics. Such models are adequate for modeling system which exhibit threshold behavior. If more detail is required, the Boolean networks should be generalized to include real-valued case.

In this paper, interpolative Boolean networks are proposed. They are based on interpolative Boolean algebra and represent a Boolean consistent real-valued generalization of classical Boolean networks. The validity of their dynamics is secured thanks to the principle of structure functionality. Furthermore, the proposed approach is a generalization of the current state-of-the-art approaches (e.g., Boolecube/Hillcube) and offers a guideline for the appropriate realization of generalized product. Also, IBN utilizes comprehensible and direct IBA transformation procedure from logical function to GBP.

Interpolative Boolean networks can be used to predict behavior of complex systems. For the purpose of their practical application, a simple software support tool is developed. The proposed approach is applied on two examples from the literature and compared against other real-valued logical approaches. Interpolative Boolean networks offer superior descriptive power compared to other approaches. This approach facilitates modeling and analysis of systems too complex for quantitative mathematical approaches.

#### Appendix

#### A. WNT5A Gene Regulatory Network Realized Using IBN

The WNT5A gene regulatory network realized as IBN is defined by the following equations:where and .

#### B. *Drosophila* Segment Polarity Genes Regulatory Network Realized Using IBN

The* Drosophila* segment polarity genes network realized as IBN is defined by the following equations:where * *, is the network state at time , and .

Taking into account the nature of the variables, this IBN is realized as follows: