Journal of Complex Systems

Volume 2014 (2014), Article ID 693938, 17 pages

http://dx.doi.org/10.1155/2014/693938

Research Article

Conception of Biologic System: Basis Functional Elements and Metric Properties

Garri Davydyan

Appletree Medical Group, 1902 Robertson Road, Ottawa, ON, Canada, K2H 7B3

Received 21 May 2013; Accepted 16 October 2013; Published 29 April 2014

Academic Editor: Yang Tang

Copyright © 2014 Garri Davydyan. 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.

Abstract.

A notion of biologic system or just a system implies a functional wholeness of comprising system components. Positive and negative feedback are the examples of how the idea to unite anatomical elements in the whole functional structure was successfully used in practice to explain regulatory mechanisms in biology and medicine. There are numerous examples of functional and metabolic pathways which are not regulated by feedback loops and have a structure of reciprocal relationships. Expressed in the matrix form positive feedback, negative feedback, and reciprocal links represent three basis elements of a Lie algebra of a special linear group . It is proposed that the mathematical group structure can be realized through the three regulatory elements playing a role of a functional basis of biologic systems. The structure of the basis elements endows the space of biological variables with indefinite metric. Metric structure resembles Minkowski's space-time (+, −, −) making the carrier spaces of biologic variables and the space of transformations inhomogeneous. It endows biologic systems with a rich functional structure, giving the regulatory elements special differentiating features to form steady autonomous subsystems reducible to one-dimensional components.

1. Introduction

The concept “system,” introduced in the middle of the last century, was proposed to describe self-organizing properties of biologic matter. The essence of it is the ability of biologic objects to maintain their own anatomical and functional structure [1–3].

Until now positive and negative feedback are used to describe the regulatory mechanisms of biologic system. Existing data show that positive and negative feedback do not explain some of the physiological and medical data. It was an attractive idea to find out whether some regulatory components can form a functional basis of the system. In challenging the existing hypotheses it was found that reciprocal links have the same importance as positive and negative feedback in functional regulation of biological systems. In fact, positive and negative feedback and reciprocal links can form a functional basis of biological systems.

Regulatory functions of the system comprise numerous metabolic pathways and biochemical reactions, which, in fact, realize only a number of core mechanisms. Until now a great deal of effort has been undertaken to find these functional invariants. Positive and negative feedback are the only known functional structures with universal properties [4–8].

Widely presented as common regulatory mechanisms in general biology and clinical practice to explain normal physiological and pathophysiological mechanisms, feedback loops do not correspond to counter-directional (reciprocal) pathways.

There are some examples showing broad representation of reciprocal relationships in regulation of metabolism and basic biologic functions. The system of hemostasis consists of two “antagonistic” subsystems, clot formation and clot degradation cascades; concentration of the glucose in the blood is regulated by insulin and glucagon releasing mechanisms leading to the opposite results. Parathyroid hormone increases and calcitonin decreases calcium concentration in the blood. The natriuretic peptide is opposite to the renin in vasopressor function. It diminishes the blood pressure, whereas renin increases it. Morphogenesis and apoptosis also follow the reciprocity as their leading regulatory mechanisms. Less trivial examples provide subsystems, which are not opportunistic or, figuratively speaking, not lying on the line of plus-minus relationships but rather have complementary properties and are “orthogonal.” These are androgen and estrogen promoted biochemical pathways leading to the different phenotypical and functional properties [5].

Biological meaning of the reciprocal structures is that each pair of functions, linked reciprocally, seems to be originated from a common morphological and functional root (or stem). Reciprocity implies the link between two opportunistic, differentiated from common functional root, functions.

Negative feedback is well known in general and reproductive endocrinology because of its quite demonstrative mechanism balancing activities of central and peripheral endocrine organs. For example, lack of peripheral hormones signals central organs to produce more stimulating substances. Stimulating factors activate peripheral systems, which increase production of peripheral hormones. If endocrine glands produce an excessive amount of peripheral hormones, it blocks central stimulation, which, in turn, diminishes peripheral activity.

Positive feedback shows its action in clot formation cascades. Another example of positive feedback is pulse therapy. Increasing the doses of hormones administered at short intervals improves functional capabilities of both peripheral and central systems. Positive feedback loops also help to explain quickly developing emotional reactions like acute phobias. In reproductive endocrinology, in the first stage of the delivery, release of the oxytocin stimulates uterus contractions, cervix dilatation, and effacement, which, in turn, stimulates oxytocin release. This vicious circle is maintained by positive feedback until the process switches to another regulatory mechanism corresponding to the second stage of delivery [5, 9].

Based on these observations, it is proposed that “reciprocal links” (reciprocity) as well as positive and negative feedback are the major determinants and separable functional elements of the internal self-regulatory structure of the system. Treated as separable functional elements, positive and negative feedback and reciprocal links expressed in the matrix form are identical to the basis of elements of Lie algebra of special linear group ; represents real numbers.

Lie algebra and the related group are closed sets of elements endowing positive and negative feedback and reciprocal links with some additional functional properties. It includes the ability to form integrative units through the linear combinations of matrices representing these elements [10].

Formally from the structure of the chosen basis elements of the Lie algebra it follows that the space of biologic variables has indefinite metric, meaning that the space of possible regulatory scenarios is inhomogeneous. Indefinite metric gives a broad spectrum of regulatory actions determining system’s behavior.

2. Graph Representation of Functional Structure of a System

According to the classical system’s theory, a system is a set of interrelated elements, acting together as a whole to achieve a specific outcome for the given system. As a whole a system has its own global or external function and in this sense it acts as a machine, transforming input into output. (Further in the text a Biologic System is a general term applied to different morphologically and functionally separable objects of biological origin. A simplest two-element system is defined as a pair , where is a space of biologic variables, realized as , and is a space of transformations on realized as matrices over . can be represented either by , which are (2 times 2) matrices over with determinant one, or which are linear approximations of , represented by matrices of trace (the sum of diagonal elements) zero.) The input-output relationship does not provide any information about internal regulatory elements and their structure. The simplest way to express internal structure of the system is to show it as a directed cycle graph where points are internal elements and arrows are functional links among them.

Biomolecules, cells, tissues, organs, and so forth are some examples of the elements of biologic systems. Relations among the elements can be “visualized” through the chemical affinity of molecules triggering biochemical reactions, or even through the generated physical forces applied to the tissues and organs. Detailed functional structure of a system shown by a web (simple graph) of linked to each other elements usually is not observable, while generalizations may hide and lose some substantial functional properties. Common functional or regulatory characteristics inherent to different biological objects can be presented, for instance, by closed loops indicating self-regulatory functional structure.

Most of the proposed regulatory mechanisms and theories of functional organization of biologic systems are based on negative feedback loops [4], which can be demonstrated in the following example. For instance, digestive systems in protozoa and animals are utterly different, but, at the same time, they have some common features that aimed to capture sources of energy, digest, and move partially digested substances until they can be excreted. In protozoa digestion is rudimental, proceeded not in a specialized tract but in the cytosol. During digestion captured substances form electrochemical and mechanical gradients, which create forces moving the content within the cell. This process seems to be continually adjusted by feedback loops through the sensory stimuli coming from every movement of being digested particles. By feedback loops the system forms an “optimal” trajectory, which mimics temporarily created digestive tract.

Contrary to the primitive unicellular organisms, the passage of chime in animals is implemented by specialized organ-gastrointestinal tract. Despite the complex functional and anatomical organization, motility mechanism of the intestine seems to also involve negative feedback as a key regulatory component. It certainly follows from the existence of closed filling-emptying cycles as possible regimes of regulation of the motor function of intestinal segment [11].

By now there are no conceptual models describing regulatory structure of the system realized in the space of functional elements [8]. We still apply to the classical image of a system, a “black box,” because of undetermined system’s internal functional structure (Figure 1).

Figure 1: The “black box” represents undetermined internal structure of a system.

The graph representation gives an initial approach in understanding of the basic structural differences among positive feedback, negative feedback, and reciprocal links. Schematically, the simplest internal functional structure is given by a two-point closed graph: points are two subsystems, and arrows, having plus or minus signs, are the links. Each pair of elements and arrows with fixed signs gives one functional (regulatory) unit. For example, negative feedback has arrows having plus and minus signs, while positive feedback has arrows with only positive signs. This scheme shows how each subsystem responds on exiting (+), inhibiting (−), or neutral (0) stimuli. Formally, there exists over sixteen scenarios, assuming that each subsystem is also sending the signals to itself. Consider two-element graph, when each element is directly linked only to itself and signs of the loops are opposite. If “neutral,” indirect, links are assumed to be between elements, this graph will represent reciprocal (antipodal) functional structure. In the reciprocal graph elements do not affect each other directly, so the links cannot be shown by +/− arrows and are denoted as neutral (Figure 2). However, information received through the neutral pathways eventually will counterbalance functional states of the elements and the whole system. For example, functional system responsible for hemostasis consists of two parts: clot formation and clot degradation subsystems. They belong to the same hierarchical level and targeted at the only one physiological parameter—viscosity of the blood. Each subsystem contributes in an opposite way. Together they provide a wide spectrum of regulatory responses to maintain the most favorable rheology for metabolic needs.

(a)

(b)

(c)

Figure 2: Graph representation of positive feedback (a), negative feedback (b), and reciprocal links (c).

Most probably, clot formation and clot degradation subsystems were originated from the single anatomofunctional root or stem. From the stem of inseparable characteristics the process of differentiation was followed by a long way of functional and morphological splitting. As a result, some functional components became partially independent (autonomous), comprising their own cascades of biochemical reactions aimed to the opposite goals—the clot formation and the clot degradation. On the other hand, together they provide an efficient mechanism in achieving the optimal viscosity of the blood. Logically, clot formation and clot degradation biochemical pathways should be functionally linked.

Two-element graphs showing positive and negative feedback mechanisms have found practical applications mostly in endocrinology, because of the simple and convincing demonstration of two-directional relationships between central and peripheral endocrine organs. On the contrary, the graph structure of the reciprocity does not give a satisfactory picture. The logical schemes have one substantial disadvantage—they do not show regulatory structure as a dynamic process. In addition, regulatory elements should possess some integrative properties and be able to maintain functional stability of the system.

3. Classical Dynamic Modeling of Internal Functional Elements of a System

Ordinary differential equations (ODE) give a unique means to demonstrate dynamic properties of the objects described by variables changing in time. In ODE positive and negative feedback and reciprocal links are represented by the operators (Linear operators and matrices of linear operators are isomorphic objects with almost the same properties. A matrix is an operator expression relative to some basis. For simplicity, they will be used as synonyms, if the differences are not emphasized. All matrices, if not specially mentioned, are expressed relative to the standard basis.) (or matrices) relating velocities of variables with variables themselves. Qualitative solutions of ODE can be shown by curves, which are evolutions of the system’s conditions in time.

For example, three simple differential equations

may represent positive and negative feedback and reciprocal links (Figures 3–5) in two-element systems, where the system’s elements , , and (vectors) are variables of ODE. Because two variables , , , are used for each system, the spaces of variables are two-dimensional, and the operators , , and are represented by 2 × 2 matrices with real coefficients , , and . The matrix view directly follows from the graph representation of the functional elements. The operators , , and serve as the symbolic equivalents of positive and negative feedback and reciprocal relations, because the described processes are entirely determined by these operators. Each ODE has two variables , , ; hence, functional relations between variables can be shown on the plane , . Conditions of the system are points on the plane, evolving in time, along the third, the time, and axis, so that three-dimensional processes (integral curves) are shown as the projections of integral curves on the plane (phase curves). The curves characteristics are determined by the types of the matrices, so that each matrix establishes its own dynamic law of variables behavior.

Figure 3: Phase curves of differential equation . Matrix represents negative feedback.

Figure 4: Phase curves of differential equation . Matrix represents positive feedback.

Figure 5: Phase curves of differential equation . Matrix represents reciprocal links.

Trajectories corresponding to the negative feedback (operator ) (Figure 3) show fluctuations of the variables, which remain within some restricted areas. Closed curves clearly demonstrate self-regulatory process [11, 12].

On the contrary, positive feedback (operator ) changes conditions hyperbolically (Figure 4). This behavior, if not corrected, will destroy the system. Reciprocal links (operator ) are also represented by hyperbolas, where the axes of variables are asymptotes (Figure 5). This type of regulations is not self-destructing. Interestingly, each point on the hyperbolas gives the same product , meaning that the increase of the value of any variable simultaneously decreases the value of another, so that the product remains unchanged. It can be hypothesized that evolution of the system's conditions determined by operator C leaves intrinsic characteristic of the system (subsystem) unchanged. which also leaves intrinsic characteristic of the system (subsystem) unchanged, if they measured by the scalar product. For example, components of input can be entirely used in output, or there is no dissipation of the inner energy during transformations and the energy is preserved in the final product and so on.

Unlike operators and , has a diagonal matrix. Related to this matrix operator has a special property; it divides 2-dimensional space of variables on two one-dimensional invariant subspaces, which in this case are coordinate axes.

Operator makes the subspace of each variable separable and autonomous. It acts on the vectors of invariant subspaces (eigenvectors) by changing their lengths or directions to opposite, leaving them in the same subspaces. If initial condition of the system lies on coordinate axis, it will evolve along this axis.

All linear second order ODE can be grouped on four classes of equivalence. Each class is represented by the matrix and has, related to this matrix, a phase portrait. Matrices and belong to the same class, because they can be transformed to each other. Their phase portraits are the saddles turned on 45° to each other. Matrix corresponds to the phase trajectories termed the center. The center and the saddle cannot be transformed to each other by continuous transformations, so that is not similar to or . There are no means to be smoothly transformed to or . From the matrix properties it follows that negative feedback is not typically (physiologically) transformable to positive feedback or reciprocal links, suggesting some bypassing mechanisms bridging these elements. Matrices , and are not singular and have inverse ones. It means that describing processes can be reversed.

Properties of dynamic systems based on the described three matrices are discussed in detail in numerous publications devoted to biologic and medical problems [12, 13]. However, by now it is not understood whether regulatory processes involve all three (or only positive and negative feedback) functional components simultaneously or there is a consequence of alternating regulatory commands depending on the current conditions. Another question is how many functional units are involved in regulatory process and whether their superposition is possible [6].

4. Functional Structure of a Biologic System Follows Lie Group and Lie Algebra Properties

It is known that every physiologic variable fluctuates around some equilibrium states. It seems logical to assume that regulatory mechanisms preventing structural deterioration should provide reversibility of the current nonequilibrium conditions in order that the process may evolve towards the equilibrium state.

Generally, reversibility is a result of the counter regulatory mechanisms, which makes functional states fluctuate around the equilibrium points. These fluctuations are observed for the most of physiologic variables. Based on these observations, it is proposed that regulatory process reflected in fluctuations of physiological parameters possesses mathematical group properties. (A mathematical group is a set of elements of some origin satisfying the following conditions: composition of any two elements of a group belongs to the group ; each element has inverse one ; there exists a unique element of the group 1, which is a neutral element . Under the law of composition the group is closed structure.) A formal structure of the group requires a set of elements of the group and a law of composition of the elements.

Many functional models use positive and negative feedback as regulatory mechanisms of two component biologic systems [6, 7]. However, positive and negative feedback have never been discussed in terms of functional basis elements of a system.

If the basis is proposed, it infers the existence of the space of functional elements, composed of the basis. It also implies a rule how to compose elements of this space, including the basis. Biological meaning of group structure is that at each time point the regulatory process is determined by the functional unit formed of the basis regulatory elements.

In the standard basis , of the carrier space of biologic variables , positive and negative feedback and reciprocal links are realized as matrices corresponding to the basis elements

(Matrices of operators , , and are the same as , , and .) of the Lie algebra , denoted , of special linear group .

The Lie algebra and the Lie group are closely related structures. is a multiplicative group of second order () matrices over , with determinant one. It represents a group of transformations of the states of a biologic system. As a mathematical object is a smooth surface (manifold), biologic transformations of the current states of the system can be described in terms of the smooth curves rather than discrete points separating one state from another. Unit matrix is a neutral element of the group. Multiplication of any two matrices of gives element of the group. Each matrix has its inverse. Because each matrix over belongs to , determinant one restricts the four-dimensional space of transformations to three dimensions, so that is a three-dimensional manifold.

Lie algebra is a linear approximation of . It is an additive subgroup of all matrices of the trace (the sum of diagonal elements) zero. The Lie algebra is also a linear vector space over of the second order traceless matrices. This vector space is a tangent space to each point (matrix) of . Tangent vectors can be moved to the identity point of the group, so that the tangent space at the identity contains all the algebra elements. Matrices , , and lie on the tangent space at the identity to the Lie group and form a basis of three-dimensional vector space . Linear combinations of basis elements with real coefficients , (, , ), give elements of , , . The , , and matrices are not singular and have the inverse ones.

Neighborhood of to each point gives a linear approximation of the surface surrounding this point. Locally, real objects are qualitatively similar to their linear approximations. If there is a curve on its behavior can be described in the neighborhood of each point of this curve by elements of the Lie algebra lying on the tangent surfaces to these points (Figure 6).

Figure 6: The space of transformations of biologic variables and its linear approximation tangent 3-space .

As a Lie algebra, has a special bilinear operation—not symmetrical (depending on the order of elements) Lie bracket defined as , where is a matrix multiplication. For each two elements of the algebra and , result of the bracket operation is also an element of the algebra , so is closed under the Lie bracket. Bracket operation is a criterion for the matrices to belong to the algebra. It also measures symmetry (commutativity) of matrix products ( versus ). If matrices commute, then the result of bracket operation is zero; . is not a commutative group; that is, multiplication in the algebra is not a symmetrical operation, which means that the result depends on the order of multiplying matrices.

Lie brackets of basis matrices are , , . In terms of transformations of functional states noncommutativity means the following: the resultant point of transformation ( following ) is not the same as for . Elements of the group can be obtained by exponential mapping of the algebra elements, exp: , , and (Figure 4). For , , and from follows . That is, because consists of matrices with determinant =1. Due to , each matrix of the group has its inverse. Exponents of the basis elements are , , [14–17].

is a noncommutative multiplicative group. Noncommutativity of the group elements compositions may reflect nonreversibility of metabolic pathways, biochemical reactions, and so forth. The commutator gives an effective approximation of the noncommutative operations of the group. If the group elements are curves on , composition of two elements of , (left translation by or multiplication from the left) will transform the curve to the curve . Right translation by , will not move right backwards to the place coincident to the initial trajectory . . (). Lie bracket substitutes the gaps, when results of the group operations are interpreted through the elements of the algebra.

For example, composition of two one-parameter subgroups of diffeomorphisms (action of on from the left) displaces the initial integral curve related to the tangent vector (let it be vector ) at identity point. Because and do not commute, there is a difference in the results of the actions of and , depending on the order of the composition: . The “gap” between two pathways and can be approximated by the Lie bracket of the derivations of and corresponding to the tangent vectors and , respectively:

In the small neighborhood exponent mapping transforms the structure of the Lie algebra into the structure of the group, or, more precisely, transforms an additive function and bracket operation of the algebra into multiplicative operation of the group

Basis elements , , and have special for biologic applications properties related to their matrix structure. In the standard basis , only operator has diagonal matrix form . It has two one-dimensional invariant subspaces and direct sum of which is the entire two-dimensional space . is invariant under means , (). Vectors from are transformed by into the vectors of the same subspace . So operator action on the elements of invariant subspace leaves elements in the same subspace. This is also true for , because is also an invariant subspace.

Vectors of the invariant subspace are eigenvectors, and each value of the matrix diagonal element is an eigenvalue of the eigenvectors of an eigenspace (invariant subspace). Real eigenvalue gives -times elongation of the applied vectors. If sign of eigenvalue is negative, it changes vectors directions to the opposite. Basis vectors , are eigenvectors of and , respectively, with eigenvalues and −1. will leave vectors of unchanged, because of the eigenvalue and transform vectors of to the opposite (reflection) due to negative eigenvalue sign. Applied to the basis vectors (eigenvectors), will transform to the same vector and change direction of to the opposite. Because divides on two invariant subspaces, can be considered as a direct sum of two operators acting separately on the corresponding subspaces and , so that, if , .

Biological meaning of invariant spaces suggests that some parts of the system have autonomous regulations. The term subsystem implies separable structure, which found its functional correlates in invariant subspaces. The existence of regulatory mechanisms (operators) acting separately on subsystems makes the subsystems anatomically distinguishable and functionally closed. Eigenvectors of invariant spaces can be interpreted as some conditions of the system, and eigenvalues—as some commands how to change the current condition’s magnitude or direction, not their qualitative characteristics. Positive eigenvalues increase existing functional activity. Negative eigenvalues invert the process to opposite.

In the basis the operator , responsible for positive feedback, has matrix form . Its eigenvectors , have ±1 eigenvalues (obtained from characteristic polynomial equation), so that can be transformed to the diagonal form by similarity transformations using orthogonal matrix U, . belongs to the same class of equivalence as the matrix . This transformation changes the initial coordinate system to which corresponds to the eigenspaces , . The new coordinates could be obtained either by the rotation of coordinate axes clockwise on 45°, or just by the rotations of the plane without coordinates (Figures 4 and 5). Invariant subspaces and are two lines inclined by 45° to the main coordinate axes (Figure 4).

Together operators and divide on four one-dimensional subspaces: , , , and . and transform points of living vectors from , and , , in the same subspaces , , , and . Although is transformed to the same diagonal form as , and are different objects, because similarity transformations change just the matrix forms of the operator , which is different from the operator . Unlike and , is irreducible to one-dimensional invariant subspaces, because it has two complex eigenvalues . endows with the complex structure, which can be represented on the plane of real variables. Thus, the only invariant subspace for is itself, which is a two-dimensional space over . The operator just rotates vectors from leaving the vectors lengths unchanged.

If a smooth curve on is given as a prototype of some biochemical, metabolic, or physiologic transformations, the tangent vectors to the curve at each point will give linear approximations of the real processes. The value of the tangent vector to some point on the curve is the velocity at which the curve is coming through this point. This is how matrices of , used in ODE, create vector fields or fields of forces on , the spaces of biologic variables. The smooth curves on (or ) are obtained as solutions of ODE. In other words, each (operator) matrix in ODE (, , ) creates its own vector field by assigning to each point of the space of variables a vector (derivative), which is a velocity of the curve coming through this point. A Lie algebra is a derivation of a Lie group, thus at each point of the group there is a tangent space of the elements of a Lie algebra. Each element of tangent space is a matrix corresponding to the vector field. Hence, every element of the algebra has a corresponding vector field. The value of this field, related, for example, to the matrix ) at arbitrary point is defined as . Matrix or related to vector field determines evolution of along integral curves , which are solutions of the differential equation .

Obtained integral curves lie on the surface .

Correspondence between matrices of the algebra and vector fields tells us that matrix structure determines specific for each matrix physical (physiological) action.

In the standard basis , , the values of the vector fields, related to , , and at some point are

Written in the differential form , , and these vector fields will serve as infinitesimal generators of the group action on the space of biologic variables . It can be shown that ; that is, vector fields with the basis elements form a Lie algebra isomorphic to the algebra .

Because is a vector space over , linear combinations of with real coefficients will give matrices from ; . They will correspond to the integrated vector field generated by on (Figure 7).

Figure 7: The basis matrices , , and of represent corresponding regulatory elements of a biologic system. Integrated functional unit is a linear combination of the basis.

Being a linear combination of , assumes a simultaneous action of three basis elements, so that each condition of two-element system is determined by the states of the three basis regulatory subsystems, functioning together. These subsystems are positive feedback and negative feedback and reciprocal links. They link morphologic elements in a whole functional structure. The three basis elements of the Lie algebra are “infinitesimal” analogs of positive feedback and negative feedback and reciprocal links.

Points of the space of two variables , determine all possible conditions of the system. Variables change their values in time, and these changes for and are points on the curves , . The characteristics of the curves are determined by the matrix , which also determines functional relations between and . Due to nonsingularity of , these relations can be expressed by a function . determines a family of trajectories on . Any of them refers the process to some initial condition . The phase trajectories lie on the 2-dimensional surface of two variables. They are projections of the integral curves , , which are solutions of the system of the two ordinary differential equations , , where , , ,  . Transforming to one of the Jordan forms we obtain eigenvalues and and eigenvectors as new variables , . Integral curves and are one-parameter group of diffeomorphisms . Each function transforms corresponding variable from , satisfying group properties, hence determining two-directional process. A curve is a direct product of two one-parameter groups of diffeomorphisms , on and .

The Lie algebra can act not only in the space of biologic variables , but also on its own elements and elements of the Lie group . Thus we can apply some regulatory structures to the primary regulatory elements themselves, considering it as the next hierarchical level of regulations.

Using matrices one can create vector fields on elements of . For the space of matrix elements and fixed real matrix , , let be a vector field on generated by whose values on each point (matrix) are

We relate the defined vector field to the ODE , solution of which is integral curves on :

The curve is a result of an action of the one-parameter subgroup of diffeomorphisms ()  generated by on element of the group . The vector is a tangent vector to the surface at the point , and it is a velocity at which is coming through that point . The points of the curve created by the vector field are obtained by multiplying the matrix exponent on .

By applying matrix to the vector field from the left an integral curve is transformed to the new point :

This is called left-invariance of the field. Because is an element of the Lie algebra , lies on the tangent surface to the group at the identity point, and due to the left invariance, creates vector fields related to on the entire group . There are also commutation rules for vector fields created from the elements of the algebra

so that vector fields on have the structure of a Lie algebra .

If , , and are the basis vectors of , related vector fields , , and are linearly independent, and at each point they present the three basis elements , , and of three-dimensional vector space tangent to at [18, 19].

5. Actions of the Algebra on the Space of Biologic Variables and on Its Own Elements

It is proposed that the operators (matrices) , , and represent regulatory basis of the two elements systems. The matrices over transform points of a two-dimensional carrier space of biologic variables. For the genetically fixed regulatory mechanisms, which were assumed to be positive and negative feedback and reciprocal links, the corresponding morphological structures should also be fixed, such as, for example, central and peripheral endocrine glands, realizing negative feedback regulatory “algorithm,” and antipodal structural elements, which are ready to release counter acting substances—insulin versus glucagon, vWF (von Willebrand Factor) versus anticoagulant proteins C and so on. Although each of the three basis regulatory elements is acting on two element structures and links only two variables, it is not evident that substantially different mechanisms simultaneously or intermittently can be applied to the same two variables. It seems unconvincing that, for example, central and peripheral endocrine glands, normally operating through the feedback loops, are also under the direct control of reciprocal links and positive feedback mechanisms. Negative feedback implies hierarchy; that is, corresponding variables should always be on different functional levels. On the contrary, the reciprocity links variables, which are on the same line of functional capabilities, so they are on the same hierarchical level. It seems to be also true for positive feedback in the sense that each of the basis regulatory elements fixes some pairs of variables becoming a basis for corresponding two-dimensional subsystems. It is assumed, that each of the operators , , and has corresponding two-dimensional space , and the sum of these spaces is a six-dimensional space of variables . On the other hand, , and, due to correspondence between and , it follows that , and the dimension of is less than six. Nevertheless, we will consider as a direct sum of , , which is a six-dimensional space, because as the basis elements of a Lie algebra are linearly independent. Additional, physiological reason for is that each is an “irreducible” unit that operates on “its own” two-dimensional subspace .

From this it follows that the “integrated” operator on is a direct sum of , , , . In the standard basis of the matrix of has block-diagonal form and the blocks are matrices , , and (Figure 8).

Figure 8: Direct sum of corresponds to space of biologic variables.

If pairs of variables are fixed for each basis regulatory mechanism ( matrix), then the six-space of biologic variables could be stemmed by three pairs of antipodal, reciprocally related, variables linked by negative and positive feedback loops. This hypothetical network is demonstrated in Figure 9. There are examples of the mechanisms maintaining rheology of the blood, which support this functional structure. According to the current understanding of these mechanisms fibrinolysis and coagulation are regulated by positive and negative feedback loops acting simultaneously on each of these components [20]. In order to comprise closed regulatory structure, elements of fibrinolytic and coagulation cascades are supposed to be connected by reciprocal links. That is, the set of either , , , or , , variables shown in Figure 9 can be considered as a clot formation or clot degradation subsystems linked by (reciprocal) types of relations. Physiologic meaning of (= reciprocal links) is that, if, in the normal conditions, coagulation begins to be a prevailing component, it will stimulate a fibrinolytic activity and simultaneously slow down the clot formation processes. Another example, demonstrating possibility for each subsystem to have different types of simultaneously acting regulatory loops with other subsystems is a hypothalamic-hypophyseal-thyroid axis. A thyroid gland linked by negative feedback loops with hypothalamus, at the same time, is connected by reciprocal relations with parathyroids in the regulation of calcium metabolism. Production of calcitonin by C-cells of thyroid glands is also regulated by negative feedback through the concentration of calcium in the blood. On hypophyseal level ultrashort feedback loops are found as additional to the negative feedback with thyroids, means of regulation of thyroid stimulating hormone (TSH) secretion [21, 22]. These data support the existence of different regulatory mechanisms applied to the same anatomical structure (or subsystem).

Figure 9: Schematic representation of the functional regulatory structure of a biologic system.

The graph shown in Figure 9 is a single regulatory unit. It is closed and contains minimal number of links corresponding to the three basis operators. If only one variable is considered in a stem, then functional unit including all three regulatory components will correspond to four-dimensional carrier space. It is clearly seen for the pairs of variables and links that originated from only one ( or ) variable (Figure 9).

If are given as the basis for , the linear span is a three-dimensional linear space of elements over , . In this case the carrier space remains two-dimensional, because is (traceless) matrix. represents the action on two-dimensional carrier space embedded in the six-dimensional space.

The algebra itself, the space, can be considered as a carrier space for some transformations . In fact, adjoint transformations of elements (discussed further in the text) have the same algebra structure (it is ) as transformations of biologic variables.

As soon as fixes the basis elements of the space of biologic variables , the spaces of biologic variables and the space of transformations, which are the space of the algebra elements spanned by , , and , are related.

There is still an ambiguity in determining the dimension of a carrier space of biologic variables. Since is a linear combination of basis transformations , , and , composition of these elements gives integrated matrix structure that predisposes simultaneous actions of . Therefore, the all three basis elements , , and as an integrated unit, should be linked through the same variable or relate the same two variables. On the other hand, each , as emphasized before, acts on its own, corresponding to , space of variables , thus making the whole space not two, or even four-dimensional, but, because of symmetry of reciprocal links, a six-dimensional (Figure 9).

To clarify this uncertainty two separate time intervals in the course of phylogenetic development should be considered, early stage, consisting of rudimental functions, and advanced stage, when mature structures have already been formed. In the first (early) stage of functional isotropy conditions of the system were, for instance, points of -dimensional space . In this stage there was no means to regulate each of parameters separately, because is the only invariant space, and it was irreducible to the invariant subspaces of lesser dimensions. For simplicity, , so we can describe internal regulatory structure through the phase trajectories showing changes in the system’s conditions. In chosen for basis the matrices of transformations of will have common form , where is a group element. Because reflects regulatory structure of mature (not growing) biologic objects, , meaning is volume preserving transformations, preventing systems from being expanded or collapsed. Unit determinant restricts 4-space to the 3-space of transformations.

Depending on the values of the matrix entries, can be transformed to the simpler forms. There are three possibilities determined by the roots , of a characteristic polynomial equation:

. are complex numbers, , ; then is an elliptic element.

Consider that . are +1 or −1, , or , where is a parabolic element.

Consider that . are different real numbers, , where is a hyperbolic element.

Similarity transformations resulted in any of , , change the initial basis in which was expressed. And, in these new bases and related matrix forms , , , may have one-dimensional (, ) or two-dimensional () invariant subspaces. Importance of this for morphological and functional development is that the initial space becomes reducible to the subspaces of lesser dimensions. These subspaces can be maintained and developed further as autonomous systems. Thus functional properties of operators of implies the existence of three types of regulatory structures. In biological terms, regulatory structure of the system, imitating the group structure of gives three autonomous branches of transformations from the primary stem of not yet separable functions.

During phylogenetic development the changes in the bases of the carrier spaces eventually are becoming fixed genetically in new anatomical structures having autonomous features corresponding to invariant subspaces. Primary morphological and functional stem has been split on two distinguishable branches with respect to the invariant subspaces. Thus new anatomical structures have become supported by steady regulatory mechanisms. They were assumed to be negative feedback, positive feedback, and reciprocal links.

The three basis operators comprising elements of the Lie algebra are linear analogues of these phylogenetic mechanisms.

Matrix , , , , of transformations on is a linear combination of the elements of basis. The “integrated” matrix can be transformed to the diagonal (Jordan) form and be expressed through the basis of eigenvectors. To obtain this matrix, we have to find the roots of the characteristic equation . The roots , of the characteristic polynomial equation are eigenvalues of the vectors of invariant subspaces: . Consider . Depending on the values of , , , there are three possibilities: can be real, zero, or complex numbers. For a more demonstrative picture, can be viewed as an operator of .

Scenario 1. are real, if ; . Therefore, operator may have two invariant one-dimensional subspaces containing eigenvectors with real eigenvalues. Formally, that means that can be transformed to the diagonal form by similarity transformations , where for this case is an idempotent orthogonal matrix . Eigenvectors for the operator are , . In new coordinates , , and eigenvectors are orthogonal and can be taken as coordinate axes for the whole space (Figures 4 and 5).

Scenario 2. are complex numbers, ; . The entire two-dimensional space is an invariant space (Figure 3).

Scenario 3. , , , decays on two cases with the singular matrices , or . In the first case all the points of the space are stationary points. The second one, related to the nilpotent operator, is a marginal situation for the improper node when two equal positive roots aimed to zero. Its eigenspaces are stationary points of horizontal coordinate axis. Nilpotent operator has only one invariant one-dimensional subspace. All these scenarios are examples of similarity transformations , .

Each value of a physiologic variable reflects some condition of a system, which is the result of interactions of comprising the system subsystems. Superposition of, for example, subsystems organized in 2-element structures is equivalent to the interactions of or variables, depending on whether subsystems have common elements. It is postulated that for two-element subsystems only three types of regulations, negative feedback (), positive feedback (), and reciprocal links (), are genetically and anatomically fixed. Other forms are “linear combinations” of these three.

A system described by variables can be reduced to the group of simpler two-component structures; it is easier to trace a behavior of two-element system as well. Thus, we consider two-dimensional space as a space of pairs of biologic variables , of a two-element system. For example, these components are {  = clot formation system, = clot degradation system}, {  = insulin-dependent system lowering the glucose of the blood, = glucagon-dependent system facilitating glucose production and increase of its concentration in the blood}, {  = system, eliminating from the blood (decreasing its concentration), system retaining in the blood}, and so forth.

Regulatory process keeping pairs of variables within optimal levels is determined by three operators , which are the basic functional components of a SYSTEM as a regulatory machine.

The following summarizes the main features, the “logic,” of three-component regulatory system. Consider a linear two-dimensional space of biologic variables and . Regulatory actions change conditions of the system according to some phase trajectories , which are projections of on . Consider that at . determines the vector field on . is an arbitrary three-dimensional vector over , . In standard basis on operator has a form , , , , which is a linear combination of the basis vectors , , of the Lie algebra , tangent to the space of transformations . Basis vectors on are considered as biologic variables , so that the regulatory process relating and is viewed according to the matrix tangent to . (obtained, e.g., from experimental data and then linearized) is in general view, which does not show regulatory structure of . There are three possible types of similar to simpler matrix forms. They are diagonal, symplectic, and nilpotent matrices. For example, , because of the values of coefficients, is transformed to the diagonal form . According to that, initial basis changes, say, to , where the basis elements are linear combinations of . This means that in the new variables , operator has the matrix form so that the same regulatory process, but in terms of new variables, has become viewed as if consisting of two reciprocal (antipodal) subsystems. Relative to that, the new basis vectors will represent invariant subspaces; that is, only in a new variables , it is possible to distinguish autonomous subsystems, regulated by . If intrinsic characteristics of were related to or structure, the new variables will be linked by negative feedback (), or some degenerative form (), needed special consideration.

Similarity transformation transforms the basis to a new one . Geometrically this transformation is related to the rotations of the three-dimensional surface to the point in which , a two-dimensional surface of 4-space of linear transformations embedded in , coincides with or will be parallel to , , or . The curve in (or phase trajectory in ) related to makes an additional movement in in order to match with corresponding pair of variables (Figure 10).

Figure 10: Similarity transformations change initial basis of a carrier space as well as the basis of the space of transformations.

It can be summarized by the following.

Proposition 1. For each basis of a carrier linear space there exist an orthonormal basis of transformation matrices on , which can be transformed to the basis, having the view , , , . This transformation is coherent with the transformation of the basis elements of the carrier space .

6. Automorphisms of Lie Algebra Imply the Existence of Uniformly Organized Regulatory Elements of Biologic System on Different Hierarchical Levels

Inner automorphisms of the group , , , and automorphisms of the algebra , are not commutative operations, so and for the algebra elements Lie bracket is a precise means to measure the “commutation gaps,” which are elements of the algebra as well; .

The Lie bracket can be viewed in the other way, as linear transformations on elements of the algebra. By definition, adjoint action of on algebra , , is a map, which with some fixed element of the algebra , sends elements to . Thus . Each has its adjoint representative in ; that is . Consider that follows from Jacobi identity, so is a homomorphism of and linear transformations on . The basis gives three adjoint matrices of transformations on , , (Empty spots of the 3x3 matrices are zeros), which form a basis of a Lie algebra of endomorphisms on .

Adjoint representation of the algebra are related to adjoint representation of the group through the exponent mapping exp: , and differentiating returns it back to , .

Unlike similarity transformations of elements of linear spaces, inner automorphisms of , , are not precise: for example, transforms to the different point . If elements of the algebra are to be transformed, for instance, , the results of inner automorphisms will permute elements of the tangent space at the neutral element of the group, sending to , . is defined as a derivation of [18].

Restrictions of () operators on two-dimensional spaces give and , which differ from the basis elements of    and only by integer 2 and negative sign. Thus, the structure of the basis operators on and , which is considered as the next, higher, and regulatory level, remains in general the same.

7. Metrics of the Space of Biologic Variables and the Space of Transformations

An interesting idea is to define quantitative characteristics of the functional elements of the system in the context of geometrical structure of the space of regulatory functions. The metric gives such a measure of intrinsic properties of geometrical objects expressed through the lengths of vectors and norms of matrices. Another meaning of the metric is related to the characteristics of the spaces as geometrical surfaces in which matrices of transformations operate. The significance of metric for biologic objects is that it tells us what functional structure and regulatory characteristics a biologic system are preserved. The metric gives a natural way to find invariants of regulatory mechanisms.

In general, metric gives some quantitative measure to the geometrical objects. Formally, metric is a symmetrical bilinear form . It relates pairs of vectors from to real numbers by the Gramm’s matrices . In Euclidean space , a bilinear form is given by a scalar (inner) product of two vectors and . For example, the square of a rectangle is a scalar product of its sides. In general, it can be written in the form , where is a metric matrix. If is related to Euclidean space, , , then . The length or the norm of the vector is defined as a square root of the scalar product of the vector with itself . This gives the Pythagorean formula to calculate the lengths through the vector’s coordinates . This is an example of positive definite metric since all signs of the diagonal elements are positive numbers. The lengths of the vectors related to this metric are always >0.

If some of the signs of diagonal elements of a metric matrix are negative, metric is called indefinite, and the objects may have positive, negative, or zero lengths, depending on the values of components. For example, gives and minus sign at the first component of the sum gives all three possibilities. The relativistic Minkowski space-time has indefinite metric , which makes a “curved” space.

By definition, is a metric tensor of type (0, 2), which determines a tensor field on a manifold by assigning to each point of the manifold a scalar product ,, on the tangent space to this point. If , , and are the basis of , the metric of the space is indefinite with the signature (1, 2); that is, the metric matrix has one positive and two negative elements, .

To show it, consider transformations , where is a one-parameter subgroup of diffeomorphisms on , and is an element of the tangent space at the identity element of the group . If , , differentiation on will give , . Due to , endomorphisms on associated with , , will preserve quadratic form , thus metric on . Determinant is a quadratic function on matrices associated with their intrinsic properties and it is a metric invariant of adjoint transformations , because . For arbitrary , over , . This metric is pseudoeuclidean with signature (1, 2) or (+, −, −) [14, 15, 19]. The norm (or length) of the vectors of is defined as .

The Killing form is another means for determining the metric on . It is defined as . For , gives indefinite metric of the same index 1, if the signs of elements change to opposite. Its matrix form is . Metrics (+, −, −) and (−, +, +) describe similar geometrical objects, because of the same index of these forms. The index is the odds of numbers of positive and negative signs.

Every algebra element , a matrix of linear transformations on the carrier space , , is related to a vector field on . can be transformed to the three of four possible canonical Jordan matrices. Two of them have the same view as , and differ only by integers. There are two basic types of phase curves obtained from the matrices —the center (for ) and the saddle (for , ). Recall, these matrices determine functional relationships between variables shown as phase curves. The curves characteristics are related to the two types of quadratic forms: and . They are the constants of the motions along the curves hence isometries preserving some metric properties on expressed as quadratic forms—the radii of the circles in the first case, and the constant products () of the values of variables related to points on hyperbolas, in the second case. Because the describing processes are volume preserving. induces transformations of points on one-dimensional sphere; hence, is metric preserving. is related to the positive definite (Euclidean) metric with signature (0, 2). are related to the transformations on hyperbolas and determine indefinite metric with signature . They are also metric preserving. Thus, depending on the Jordan forms to which can be transformed, , considered to be a two-dimensional, becomes equipped with positive definite or indefinite metrics. There is one more class related to the nilpotent matrices, which is the result of similarity transformations of . Because nilpotent matrices are singular, the constants of motions are stationary points on the plane, which all are the points of axis . In some sense, the space of nilpotent matrices connects the first two cases.

Suppose to be a 2-dimensional linear space over of biologic variables relative to the standard basis . We described two types of metric on depending on Jordan forms obtained through the similarity transformations of . These metrics are indefinite , related to the diagonal form , and Euclidean , related to . Depending on , resembles hyperbolas or spheres. is also a two-dimensional surface, points of which are realized as conditions of two element system. These conditions are related to the points on a phase curve, which in turn, is determined by a matrix of transformations on . This matrix is an element of a three-dimensional manifold , . At the identity element of the group , , there is a tangent space , which is the Lie algebra of matrices of zero traces. Elements of the algebra form a three-dimensional vector space . Each matrix , is a linear combination of three matrices , , , considered as a basis of the algebra. At each point in the space of variables matrix creates a velocity vector, and through the exponential function is related to the integral curves, lying on .

Recall, the norm of the elements of is defined as . Let the scalar product on be corresponding to the Killing form. Henceforth, is a metric space. For the basis elements of the algebra (tangent space ) we have , where if , , , . Computations show , ( is a Kronecker’s delta), so that is orthogonal to if . Geometrical surface related to this metric is a hyperboloid. Indefinite metric of makes the space of regulatory elements inhomogeneous.

The value of the matrix determinant, , , is associated with the value of the scalar product . It can be positive, , negative , or zero .

According to terminology of the special relativity the vectors related to negative space-time interval or negative scalar product, , are time-like vectors , positive scalar product, ; space-like vectors and zero square lengths, ; light-like or null vectors . Hence, there are three types of subspaces determined by metric. The subspace to which vector belongs called its causal character (Figure 11) [14, 16].

Figure 11: Indefinite metric (+, −, −) divides the space of functional elements on three subspaces: time-like, , null, , and space-like, .

While indefinite metric gives all three possibilities, positive definite metric is always related to the space-like vectors, because the sum of the squared lengths is never negative. Given a span , we have . Restrictions of for two-dimensional spaces will give indefinite metric with signature and Euclidean metric , depending on the pairs of the basis elements for two-dimensional subspaces. Let be subspaces spanned by and .

Then and are indefinite metrics with signature . Vectors from the subspaces and are time-like vectors. Each of these subspaces contains two null vectors: for , , and for , , where , are basis vectors for and , respectively.

Consider that is a positive definite (Euclidean) and vectors from are space-like.

Consider as representatives of subspaces , , and . We have and in the new notations correspondence between the metrics and the subspaces that are is metric of time-like subspaces and . Time-like subspaces and will be denoted . Time-like subspaces contain null vectors. The light-like subspace is denoted . Consider that is a metric of a space-like subspace .

According to the three previously described scenarios, vector over , being in general position, can be transformed to one of the following forms: , , or . This transformation rotates the tangent space to have a correspondence between and one of the three subspaces , , and (Figure 10).

Vectors from the subspace , , have negative squared lengths, because of the condition for the diagonal form , equivalent to . is spanned by vectors with positive squared lengths due to . These vectors have the form . Vectors from are nilpotent matrices that have zero square lengths.

Direct computations of scalar products using normalized , , and vectors as representatives of , , and show that . (here and further in the text in the parentheses are the matrix columns)). Consider that is indefinite; —is a space of time-like vectors. Consider that is not a degenerate metric and can be transformed to the diagonal form . Its signature is , so vectors of are time-like. Consider that is degenerate; —is a space of light-like vectors. Consider that is nonsingular and can be transformed over by to the diagonal form , thus leaving unchanged. Consider that , .

Consider that , , and , . is similar to , to , and to . Hence is space-like. Thus, the metric is congruent to , and is an isometry.

Analogous to the first, initial, partition, we have , , and , where is either or . is a space of light-like vectors.

It can be shown that Lorentz vector space , being the sum of time-like (), light-like , and space-like () subspaces, can be presented as a direct sum of just space-like and time-like subspaces [14].

In biologic terms, indefinite metric with the signature and positive definite metric are the most consistent (stable) structural components of relationships among basis regulatory elements and the most distinguishable feature of the integrated regulatory structure of the system. and are the functional invariants under similarity transformations. It should be noted that metric matrices and matrices of transformations are different objects.

Indefinite metric of the space of regulatory elements is a consequence of the structural properties of the basis elements , , and of . We described indefinite metric also on the space of biologic variables. It was shown that can be furnished by Euclidean metric as well.

The carrier space inherits the metric structure of regulatory elements. Each type of the basis regulatory element is fixed by some anatomical (or molecular) pairs of variables, which groups elements of the whole anatomical and molecular structure into the families of negative feedback, positive feedback, and reciprocal links. Different basis regulatory elements comprise different families; that is, a pair of anatomical structures linked by feedback cannot be regulated by reciprocal relations, and vice versa. More specifically, for example, TSH (thyroid stimulating hormone) and thyroxine (hormone of the thyroid gland) are related through the negative feedback loops, not through the reciprocal links. Therefore, it is logical to assume that there are families of 2-space , isomorphic as linear spaces , but are different as biologic and geometrical objects. In fact, we have a direct sum of two-dimensional spaces representatives of each family: , . There are convincing examples for two-dimensional spaces of variables, related by negative feedback and reciprocal links. It is proposed that positive feedback has also representation in “its own carrier space.”

The matrices have dual features. As elements of a Lie algebra , considered as a “one-dimensional” basis of . If are the elements of four-dimensional space of transformations on whose matrix structure is related to the type of operator , then each is a restriction on and the direct sum of , will correspond to space of . This correspondence takes place because of induced structure of . In this case each , being a linear combination of the unit matrix elements (a standard matrix basis of ), belongs to . Two bases of , , where and , are isomorphic, and we have correspondences , , and . Direct sum of the operators representing negative feedback and reciprocal links will certainly make sense, because and as well as their subspaces are disjunctive . Although it was proposed that positive feedback has its own representation in the carrier space , (nilpotent operator ) and subspace (related to positive feedback) is transversal to , but not disjunctive, . On the other hand, nilpotent operators , and subspaces of light-like vectors related to are the physiologic background of distinguishing three orthogonal, disjunctive carrier subspaces.

The sum of the subspaces of transformations corresponding to the three basis matrices presents the same 3-space as a direct sum of subspaces representing negative feedback and reciprocal links: , , and . It means that and play the role of “pure” regulatory states that span : the space containing vectors are orthogonal to the space of and complements it to . Thus, the whole space can be entirely decomposed on two subspaces, having no common regulatory components in the sense that and are disjunctive as 4-vector. In other words, regulatory actions can be a superposition of negative feedback and reciprocal links. Thus the metric of subspaces of biologic variables , related to space of transformations , can be presented as having just indefinite or positive definite metric structures. Since can be continually transformed to , it seems that its functional role is to alleviate relationships between negative feedback and reciprocal links. Due to the same sign of determinants, and can be transformed to each other by continuous transformations. Unlike , and cannot be connected by a continuous pathway. In this representation matrices, representing , , , are embedded in space of transformation matrices of space of variables as diagonal blocks. Each diagonal block is a separable component of an integrated functional unit. Matrices of the tangent space being in general position relative to the basis of can be transformed to the Jordan forms, so that in a new basis to match one of the diagonal blocks, where blocks are , , and . Related to the new pair of variables (, ) the integrated subunit will have a block-diagonal view, where the blocks are or or and other diagonal elements are . Consider that , , and , whose entry, .

One of the most important functional determinants of the system is stability of its components and the whole structure. In the context of the presented materials stability is discussed in two ways, with respect to the properties of regulatory elements to be a mathematical group, and in the connection with specific features of mathematical operators and matrices themselves. The group structure of regulatory elements explains why current conditions of the system fluctuate around equilibrium states. Besides being a group, these operators have additional characteristics. For instance, an apoptosis, which is the maintenance of the equilibrium between cell proliferation and elimination, can be viewed as a distinct function fixed in reciprocal links between two opposite processes. The property of an apoptosis to prevent tissue overgrowth can formally be derived from matrix having . For example, uncontrollable growth, known as one of the peculiar characteristics of malignancy, can be related to the indefinite metrical structure of with the only but crucial difference; there should be . The matrix form could be , . represents negative feedback, which is similar to the reciprocal links in keeping variables within certain functional limits. This is also due to . needs a confined area, because the variables have the tendency to increase their values simultaneously and almost linearly when they come up to the asymptote. Because the phase curves of are hyperbolas the metric properties of the phase space are preserved, if the asymptotes serve as the coordinates.

The phase curves related to , , and are the examples of isometries preserving inner characteristics of subsystems to make regulatory process autonomous. Thus, stability of biologic systems can be described in terms of isometries of linear transformations. An isometry is a mapping preserving some geometrical properties of transforming objects given by the norm of vectors and operators or by the scalar product. In other words, isometries give additional characteristics to the regulatory elements depending on the type of geometrical surfaces of the carrier space and the space of transformations. For Euclidean spaces isometries are obtained by a group of orthogonal transformations, preserving, for instance, the radiuses of rotations. Isometries of the spaces with indefinite metric are related to the group of Lorentz transformations.

An isometry can be defined as a transformation leaving metric tensor unchanged, so that to satisfy . For two-dimensional spaces equipped with indefinite metric, , : , the matrix of isometry, or Lorentz transformations is found from the above equation. It has the view . This matrix represents the group of hyperbolic rotations over the angle . It transforms points lying on hyperbola to the points of hyperbola, thus leaving them on the same geometrical surface. Hyperbolas are pseudo spheres (dimension = 1) with real or imaginary radii satisfying or . For , correspondent surfaces are one or two chamber hyperboloids or . In other words, hyperbolic rotations preserve indefinite metric of the above quadratic forms given as the formulas of hyperbolas or hyperboloids.

Consider that gives a matrix of isometry , because is preserved by Lorentz transformations . It is known that is not a connected group and, depending on the combinations of signs of its elements, consists of four components , represented by simpler forms , , , and , respectively. It should be noted that and matrices have the same view as , considered previously as matrices of transformations on . Thus the main feature of is reciprocal transformations of subsystems conditions, preserving metric properties of the whole space.

Isometries related to , the matrix of a negative feedback, are a group of orthogonal transformations, satisfying ; , : . The matrix of isometry satisfying the above equation is . The matrix () belongs to the group of proper rotations of the plane through the angle preserving Euclidean metric that givesconcentric circles on the phase space of two variables. Consider that gives an orthogonal isometry matrix with the property . determines rotations preserving Euclidean metric.

Consider to be a space of transformations on of biologic variables, as a metric space. Its metric tensor, defined at each point of , has a signature (1, 2) inherited from the indefinite metric of the tangent spaces to these points. Tangent spaces are the algebra spanned by the basis . Because for each , and inner endomorphisms of , , give elements of the group , , preserve determinant function , thus can be considered as the volume preserving isometries.

8. Regulatory Functions of a Biologic System as Inhomogeneous Space of Metric Matrices with Preserved Isometry

Regulatory structure of biologic system can also be considered in a way, when regulatory elements are represented by metric matrices. For these purposes is realized as a three-dimensional space of metric matrices. Matrices of transformations , , and become Gramm’s matrices, denoted , , and of bilinear forms. Each has its own metric structure: , are indefinite metric matrices; is symplectic matrix. , quadratic forms on have the same features as in cases where , are transformation matrices, whereas , the symplectic metric, can be considered as a generalization of rotations in Euclidean spaces. Symplectic metric does not “separate” variables, so regulatory function, for example, in , can be applied only to pairs of variables. Scalar product of one-dimensional objects related to symplectic metric is zero; . On the contrary, indefinite metric distinguishes the class of the one-dimensional objects, so that the division of a carrier space can be reduced to the invariant “autonomous” subspaces of lesser dimensions.

Metric can be determined equivalently through the determinant function or the Killing form. The basis endows with the indefinite metric , so that is a three-dimensional Lorentz vector space (−, +, +).

Metric matrices can be transformed either to or forms, so that resulting hyperbolic or symplectic spaces become projections of initial metric structure on corresponding two-dimensional carrier spaces. Depending on the basis of , hyperbolic space can be presented by or matrices.

In general, isometries of hyperbolic and symplectic spaces can be characterized by transformations preserving the forms of corresponding metric matrices: gives matrix of isometry , because is preserved by ; is related to the symplectic matrix of isometry with the property . It has the same view as . Transformations preserving follow from , where can be represented by the matrix of the same view as .

Metric properties of the carrier spaces of biologic variables are inherited from metrics of corresponding regulatory structure. and , related to reciprocal links and positive feedback, have indefinite metric structure, whereas , related to negative feedback, has symplectic structure.

Metric properties of can be summarized in the law of conservation of energy, expressed through the ordinary differential equation , where , , and are real numbers. This equation shows that change in the state of biologic system in a time interval is determined by the integrated vector field . This equation describes autonomous system not explicitly depending on , and it is a linear approximation of “real” processes.

Each corresponds to the vector field related to the constant first integral (or constant energy level); therefore integrated field over will determine a conservative, closed system with the constant inner energy [16, 19, 23]. Depending on the values of coefficients , , , can be transformed to , , or matrices giving three types of solutions corresponding to the new variables.

9. Conclusion

Initially the purpose of this paper was to validate and describe the basis functional components of the biologic system, which makes the system a closed structure under the existing regulatory mechanisms. Later on, the Lie algebra formalism led to the metric properties of the space of regulatory elements of the biologic system. It seemed logical to include two closely related results in this paper.

Presented conception of a biologic system is based on a general approach describing regulatory structure of the whole system through the quantification of its functional parts. It is proposed that the basis functional elements resembling universal parts of regulatory structure (subsystems) have become separable during phylogenetic and ontogenetic development, when a manifestation of new functional properties and appearance of new morphological structures were accompanied by three (specific for structural pairs) types of regulatory mechanisms. In two element systems the three regulatory basis elements, positive feedback, negative feedback, and reciprocal links, form an integrative regulatory unit. Acting simultaneously, they determine current conditions of the system.

This conception is based on the structure of the basis elements of a Lie algebra and the related group , having the properties to be closed structures under compositions of their elements. Three basis elements of the algebra, , , and , represent two well-known regulatory mechanisms, negative and positive feedback, and the newly described reciprocal links. Together they form integrative vector field over on the space of biological variables. Each regulatory component “phylogenetically” is fixed on its own 2-space of biologic variables, so that together they form a morphofunctional six-dimensional integrative unit. In practice we understand it as if the whole 6-space is reduced to the two-dimensional subspaces (subsystems) regulated either by positive feedback, negative feedback, or reciprocal links.

This provides the clue to new understanding of the fact that the regulatory components of the system should form a closed functional structure. This is the key point for keeping the system stable.

Adjoint representations of the group and algebra elements ( and ) show the universal character of regulatory elements on different functional levels—regulation of metabolic processes and autoregulation of functional elements themselves. The basis elements of and related vector fields represent their characteristics in the ability of the corresponding subsystem to maintain their internal structure and autonomy. Thus each regulatory subsystem can be viewed as conservative and stable, capable of maintaining its inner energy, that is, to be independent of the external sources of energy. This simplification helps in understanding the fact that, despite the differences in functional properties of regulatory components, the group properties endow the integrated structure with the same features as the components. This makes biologic systems uniformly organized on different hierarchical levels.

From the structure of basis of the Lie algebra it formally follows that all events regarding regulatory processes of normal biologic systems, geometrically speaking, should lie on some geometrical surface, which in our case is hyperboloid.

Primary regulatory functions comprise three-dimensional space having indefinite metric structure with signature or (+, −, −). It resembles relativistic Minkowski space-time restricted to two spatial coordinates. Geometrically, regulatory processes are trajectories on the surface resembling hyperboloids. Scalar product, considered as a measure of functional activities of regulatory elements themselves or their interactions, divides the 3-space of regulatory elements on three subspaces containing “time-like,” “null,” and “space-like” vectors. Self-preservation of this metric and its components is an essential feature of biologic organization.

Indefinite metric structure provides functional and morphological flexibility to biologic systems, which could be reflected in the system's ability to be highly adapted to the constantly changing environment.

The meaning of indefinite metric for biologic objects is based on reciprocity, which is an efficient means in achieving equilibrium states. Formally, the reciprocal links directly reflect indefinite metric in the regulatory structure.

The basis elements can be transformed into each other, which can be extrapolated on the ability of real systems to combine potent regulatory mechanisms in order to obtain alternative pathways to compensate deteriorated function. Structural “deterioration” of any of the three basis components may affect the system's steadiness and the ability to maintain genetically fixed relations among regulatory components. Practically, it may cause either an inadequate stimulation or suppression of regulatory functions, leading to excessive activity or degeneration. It can be presented in the form of uncontrollable cell proliferation. Most probably, this process has a more complex underlying mechanism, related to the deterioration of all three basis regulatory components.

This concept can be generalized to -dimensional spaces. Almost all the statements, discussed in mathematical terms, can be easily interpreted using physiological language. This approach has a broad spectrum of possibilities for further theoretical and practical investigations. For example, it may boost to find out how the reciprocal links are related to already known mechanisms of regulations, what structural changes in the basis regulatory mechanisms determine functional and morphological deterioration, leading, for instance, to the uncontrollable and invasive cell proliferation. The latter is a clue to understand the functional phenomenon of cancer growth as a “deviation” from the uniformly organized regulatory structure of biologic systems.

Highlights

The mathematical model of organization and self regulation of biologic systems gives insight to the most common aspects of biologic organisation. Positive feedback, negative feedback, and reciprocal links expressed in a matrix form are identical to the basis elements of Lie algebra ; thus, they may represent a functional basis of biologic system. Three basis elements form closed functional structure maintaining stability of the system. The space of biologic variables is inhomogeneous and has indefinite metric structure.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The author would like to thank Dr. Janetta Kourzenkova for discussions and useful suggestions made while preparing this material. This work was not supported by grants.

References

L. von Bertalanffy, General System Theory: Foundations, Development, Applications, George Braziller, New York, NY, USA, 1973.

W. R. Ashby, An Introduction to Cybernetics, Taylor & Francis, New York, NY, USA, 1956.

M. Eigen and P. Schuster, The Hypercycle: A Principle of Natural Self-Organization, Springer, Berlin, Germany, 1979.

P. K. Anochin, Theory of Functional System, Science, Moscow, Russia, 1980.

L. Goldman and D. Ausiello, Cecil Textbook of Medicine, Saunders, 22nd edition, 2004.

B. Pfeuty and K. Kaneko, “The combination of positive and negative feedback loops confers exquisite flexibility to biochemical switches,” Physical Biology, vol. 6, no. 4, Article ID 046013, 2009.

O. Cinquin and J. Demongeot, “Positive and negative feedback: striking a balance between necessary antagonists,” Journal of Theoretical Biology, vol. 216, no. 2, pp. 229–241, 2002.

G. Grussing, Quantum Cybernetics: Toward a Unification of Relativity and Quantum Theory via Circulatory Causal Modelling, Springer, New York, NY, USA, 2000.

P. Felig and L. A. Frohman, Endocrinology and Metabolism, McGraw-Hill, New York, NY, USA, 4th edition, 2001.

G. Davydyan, “Functional basic elements of biologic systems,” Journal of Critical Care, vol. 21, no. 4, p. 360, 2006.

G. Davydyan and S. Regirer, “Modelling of the intestinal segment motions, annals of the Russian Academy of Science,” Fluid Mechanics, vol. 1, pp. 36–42, 1994.

D. K. Arrowsmith and C. M. Place, Ordinary Differential Equations. A Qualitative Approach with Applications, Chapman & Hall, London, UK, 1992.

V. I. Arnold, Ordinary Differential Equations, Nauka, Moscow, Russia, 1984.

B. O'Neil, Semi-Reimannian Geomery with Applications to Relativity, Elsevier, San Diego, Calif, USA, 1983.

S. Helgason, Differential Geometry, Lie Groups, and Symmetric Spaces, Academic Press, New York, NY, USA, 1978.

A. Kostrikin and Y. Manin, Linear Algebra and Geometry, Nauka, Moscow, Russia, 1986.

J. Stillwell, Naive Lie Theory, Springer, New York, NY, USA, 2008.

B. Schutz, Geometrical Methods of Mathematical Physics, Cambridge University Press, Cambridge, UK, 1980.

B. Dubrovin, S. Novikov, and A. Fomenko, Modern Geometry: Methods and Applications, Nauka, Moscow, Russia, 1986.

J. P. Greer, J. Foerster, G. M. Rodgers, et al., Wintrobe’s Clinical Hematology, vol. 1, Lippincott Williams & Wilkins, Philadelphia, Pa, USA, 12th edition, 2009.

M. F. Prummel, L. J. S. Brokken, and W. M. Wiersinga, “Ultra short-loop feedback control of thyrotropin secretion,” Thyroid, vol. 14, no. 10, pp. 825–829, 2004.

G. Gelwane, N. De Roux, D. Chevenne, J. C. Carel, and J. Léger, “Pituitary-thyroid feedback in a patient with a sporadic activating thyrotropin (TSH) receptor mutation: implication that thyroid-secreted factors other than thyroid hormones contribute to serum TSH levels,” Journal of Clinical Endocrinology and Metabolism, vol. 94, no. 8, pp. 2787–2791, 2009.

H. Weil, The Theory of Groups and Quantum Mechanics, Dover, New York, NY, USA, 1950.