Review and application of group theory to molecular systems biology
- Edward A Rietman^{1, 2, 6},
- Robert L Karp^{3} and
- Jack A Tuszynski^{4, 5}Email author
https://doi.org/10.1186/1742-4682-8-21
© Rietman et al; licensee BioMed Central Ltd. 2011
Received: 5 February 2011
Accepted: 22 June 2011
Published: 22 June 2011
Abstract
In this paper we provide a review of selected mathematical ideas that can help us better understand the boundary between living and non-living systems. We focus on group theory and abstract algebra applied to molecular systems biology. Throughout this paper we briefly describe possible open problems. In connection with the genetic code we propose that it may be possible to use perturbation theory to explore the adjacent possibilities in the 64-dimensional space-time manifold of the evolving genome.
With regards to algebraic graph theory, there are several minor open problems we discuss. In relation to network dynamics and groupoid formalism we suggest that the network graph might not be the main focus for understanding the phenotype but rather the phase space of the network dynamics. We show a simple case of a C_{6} network and its phase space network. We envision that the molecular network of a cell is actually a complex network of hypercycles and feedback circuits that could be better represented in a higher-dimensional space. We conjecture that targeting nodes in the molecular network that have key roles in the phase space, as revealed by analysis of the automorphism decomposition, might be a better way to drug discovery and treatment of cancer.
1. Introduction
In 1944 Erwin Schrödinger published a series of lectures in What is Life?[1]. This small book was a major inspiration for a generation of physicists to enter microbiology and biochemistry, with the goal of attempting to define life by means of physics and chemistry. Though an enormous amount of work has been done, our understanding of "Life Itself" [2] is still incomplete. For example, the standard way in which biology textbooks list the necessary characteristics of life--in order to delineate it from nonliving matter--includes metabolism, self-maintenance, duplication involving genetic material and evolution by natural selection. This largely descriptive approach does not address the real complexity of organisms, the dynamical character of ecological systems, or the question of how the phenotype emerges from the genotype (e.g., for disease processes [3]).
The universe can be viewed as a large Riemannian resonator in which evolution takes place through energy dispersal processes and entropy reduction. Life can be thought of as some of the machinery the universe uses to diminish energy gradients [4]. This evolution consists of a step-by-step symmetry breaking process, in which the energy density difference relative to the surrounding is diminished. When the universe was formed via the Big Bang 13.7 billion years ago, a series of spontaneous symmetry-breaking events took place, which evolved the uniform quantum vacuum into the heterogenous structure we observe today. In fact the quantum fluctuations of the early universe got blown up to cosmological scales, through a process known as cosmic inflation, and the remnants of these quantum fluctuations can be observed directly in the variation of the cosmic microwave background radiation in different directions. At each stage along the evolution of the universe--from quantum gravity, to fundamental particles, atoms, the first stars, galaxies, planets--there was a further breaking of symmetry. These cosmological, stellar, and atomic particle abstractions can be powerfully expressed in terms of group theory [5].
It also turns out that the very foundation of all of modern physics is based on group theory. There are four fundamental interactions (or forces) in Nature: strong (responsible for the stability of nuclei despite the repulsion of the positively charged protons), weak (manifested in beta-decay), electromagnetic and gravitational. The first three are described by quantum theories: an SU(3) gauge group for the quarks, and an SU(2) × U(1) theory for the unified electro-weak interactions [6–8]. From these theories one can derive, for example, Maxwell's theory of electromagnetism, which is the basis of contemporary electrical engineering and photonics, including laser action. Group theory provides a framework for constructing analogies or models from abstractions, and for the manipulation of those abstractions to design new systems, make new predictions and propose new hypotheses.
The motivation of this paper is to examine an alternative set of mathematical abstractions applied to biology, and in particular systems biology. Symmetry and symmetry breaking play a prominent role in developmental biology, from bilaterians to radially symmetric organisms. Brooks [9], Woese [10] and Cohen [11] have all called for deeper analyses of life by applying new mathematical abstractions to biology. The aim of this paper is not so much to address the hard question raised by Schrödinger, but rather to enlarge the set of mathematical techniques potentially applicable to integrating the massive amounts of data available in the post-genomic era, and indirectly contribute to addressing the hard question. Here we will focus on questions of molecular systems biology using mathematical techniques in the domain of abstract algebra which heretofore have been largely overlooked by researchers. The paper will encompass a review of the literature and also offer some new work. We begin with an introduction to group theory, then review applications to the genetic code, and the cell cycle. The last section explores ideas expanding group theory into contemporary molecular systems biology.
2. Introduction to Group Theory
Group theory is a branch of abstract algebra developed to study and manipulate abstract concepts involving symmetry [12]. Before defining group theory in more specific terms, it will help to start with an example of one such abstract concept, a rotation group.
Cayley table for the rotation example (see Figure 1).
e | r _{1} | r _{2} | r _{3} | |
---|---|---|---|---|
e | e | r _{1} | r _{2} | r _{3} |
r _{ 1 } | r _{1} | e | r _{3} | r _{2} |
r _{ 2 } | r _{2} | r _{3} | e | r _{1} |
r _{ 3 } | r _{3} | r _{2} | r _{1} | e |
The symmetry about the diagonal in the Cayley table tells us that the group is abelian: when the rotations are performed in pairs, they are commutative, so that r_{m} r_{n} = r_{n} r_{m}.
- 1.
Associativity: for all a,b,c ∈ G, (a *b) * c = a * (b* c).
- 2.
Identity: There is an identity element e ∈ G, such that a *e = e* a = a for all a ∈ G.
- 3.
Inverse: For any a ∈ G there is an element b ∈ G such that a*b = b* a = e.
- 1.
Rotations by multiples of 2π/n.
- 2.
The group of integers module n.
These will be discussed later.
Infinite groups are harder to study, but those that have additional structure--like the structure of a topological space or of a manifold--where this additional structure is compatible with the group structure, have also been classified. Of particular interest are the Lie groups, which are simultaneously groups and topological spaces, and the group multiplication and inverse operation are both continuous functions. Lie groups are completely classified, many of them arising as matrix groups. The matrix representation allows us to use conventional matrix algebra to manipulate the group objects, but does not play any special role. In fact any group, finite or infinite, is isomorphic to a subgroup of matrix groups. This is the realm of group representation theory.
The special orthogonal group SO(n) consists of those orthogonal matrices whose determinant is +1, and they form a subgroup of the orthogonal group: SO(n) ⊂ O(n). Geometrically, the special orthogonal group SO(n) is the group of rotations in n dimensional Euclidian space, while the orthogonal group O(n) in addition contains the reflections as well.
form a group (where H means complex conjugation of each matrix element together with transposition). Special unitary matrices, SU(n), satisfy the additional det(U) = +1 constraint, and also form groups.
Finally, we mention the "symplectic" or Sp(2n) groups, but given the fact that these are harder to define, we will not give a formal definition here. As will be shown later, these matrix groups are used in describing the "condensation" of the genetic code.
Another important definition which we will encounter later involves groupoids. A groupoid is more general than a group, and consists of a pair (G,μ), where G is a set of elements, for example, the set of integers Z, and μ is a binary operation--again usually referred to as "multiplication," but not to be confused with arithmetic multiplication--however, the binary operation μ is not defined for every pair in G. We will see that groupoids are useful in describing networks, and thus transcriptome and interactome networks.
3. The Genetic Code
In this section we review some work describing the genetic code in groupoid and group theory terms. One could easily imagine genetic codes based only on RNA or protein, or combinations thereof [13]. When the genetic code "condensed" from the "universe of possibilities" there were many potential symmetry-breaking events.
The triple cross product has 4^{3} = 64 possible triplets. As is known, the full three-way product table contains redundancies in the code. This was all worked out in the '60s, without group theory, using empirical knowledge of the molecular structure of the bases [14].
Klein group table for genetic code exchange operators.
e | α | β | γ | |
---|---|---|---|---|
e | e | α | β | γ |
α | α | e | γ | β |
β | β | γ | e | α |
γ | γ | β | α | e |
The corners of this hypercube form two octets of dinucleotides, the two sets M_{1} and M_{2} . The vertices of each octet lie at the planes of a continuously connected region. One such region M_{1} is shown in the shading of Figure 2. The octets are neither subgroups nor cosets of a subgroup. They are both unchanged under the operations (e, e) and (β,e). These two octets can also be interchanged by acting on one of them with (α,α) and/or (γ,α).
Here, we have written the corresponding dinucleotide next to the operator in the format (e,e):AA, etc.; the bar over some dinucleotides indicates membership in a different octet of completely degenerate codons, while the other dinucleotides are ambiguous codons.
to arrive at a 64-dimensional hypercube as the general genetic code. But of course multiple vertices of this hypercube code for the same amino acid. This is said to be a surjective map, because more than one nucleotide triplet codes for the same amino acid. In 1982 Findley et al. [17] describe further symmetry breakdown of the group D, and show various isomorphic subgroups including the Klein group and describe alternative coding schemes in this hyperspace.
and from a phylogenetic tree to recreate trajectories in this space. It should be possible to thus see regions of the information-spacetime that have not been explored by evolution. One may speculate on the code-trajectory by bringing in Stuart Kauffman's theory on the adjacent possible [19–21] by a perturbation theory. Further, the curves on this manifold should map, in a complex way, to the symmetry breaking described below, or bifurcation, and thus give a second route to the differential geometry of Findley et al. [18].
The initial sp(6) symmetry breaks into 6 subspaces sp(4) and su(2). Sp(4) then splits into su(2) ⊗ su(2) while the second su(2) factors into u(1). Details are given in Hornos and Hornos [22] and Forger et al. [23] on how this maps to the existing genetic code.
4. Cell Cycle and Multi-Nucleated Cells
Cell cycle is an example of a natural application of group theory because of the cyclic symmetry governing the process. The steps in the cell cycle include G1 → S → G2 →M, and back to G1. In some cases G0 is essentially so brief as to be nonexistent so we will ignore that state.
Group table for the cell cycle.
G1 | S | G2 | M | |
---|---|---|---|---|
G1 | G1 | S | G2 | M |
S | S | G2 | M | G1 |
G2 | G2 | M | G1 | S |
M | M | G1 | S | G2 |
Group table for the converged stated of binucleated cells (see Figure 4).
G1 | S | G2 | M | |
---|---|---|---|---|
G1 | G1 | S | G1/G2 | M |
S | S | S | S | M |
G2 | G1/G2 | S | G2 | M |
M | M | M | M | M |
On the other hand, it is clear from the multiplication table the we cannot have a group structure on the set {G1, G2, S, M}. Namely, in a group G any row or column of the multiplication table will contain the elements of G precisely once, hence will be a permutation of elements. This property fails for the rows of S and M. Furthermore, the product of G1 and G2 is undefined. Nevertheless, the set {G1, G2, S, M} carries the structure of a groupoid--which is discussed below.
Similar considerations apply if we fuse cells of different type, or differentiation state. These types of experiments were carried out for different stem cells, as reviewed in Hanna [28]. Another fusion-type experiment involves nuclear transfer from one type of somatic cell to another, and determining the identity of the outcome. A variant of this is to transfer RNA populations between cells, and observe the change in the cell's phenotype [29].
5. Algebraic Graph Theory: Graph Morphisms
With modern datasets it is possible to begin exploring molecular systems dynamics on a network level by using morphism concepts and algebraic graph theory. For example, using these datasets we may be able to impute missing connections in PPI networks, or build vector-matrix-based models representing the dynamics of changing PPI networks. In other cases we may be able to prove algebraic graph theory concepts using the PPI-data. Our focus here will be to continue exploring the cell cycle by including transcription data and protein-protein interaction data from high-throughput screenings. We will first review a few algebraic graph theorems. Godsil and Royle [34] will be our primary reference for algebraic graph theory.
Mathematically a network is a graph G = G(V, E) of a set of n vertices {V} (also called nodes), and a set of e edges {E}, or links. Graphs can be represented using the adjacency matrix A. The adjacency matrix of a finite graph on n vertices is the n × n matrix where the non-diagonal i-j th entry A_{ ij } is the number of edges from vertex i to vertex j, while the diagonal entry A_{ ii } , depending on the convention, is either once or twice the number of edges (loops) from vertex i to itself.
Where tr(A) represents the trace of the matrix, n is the number of edges, and t represents the number of triangles in the graph. An excellent review of spectral graph theory is given by Chung [35].
Another important matrix is the incidence matrix, which has some very useful properties. The incidence matrix B(G) of a graph G, is a matrix having one row for each vertex and a column for each edge, with nonzero elements for those node-edge pairs for which the node is an end-node of the edge. This matrix is therefore not square. An interesting property is that if we let G be a graph with n vertices, c_{0} its bipartite connected components, and B the incidence matrix of G, then its rank is given by rk(B) n - c_{0}.
The concept of automorphism of a graph is an important one, and as we will see it has applicability to subgraphs within more complex graphs. Automorphisms of a graph are permutations of the vertices that preserve the adjacency of the graph, i.e., if (u, v) is an edge, and P is the graph automorphism, then (P_{ u } , P_{ v } ) is also an edge. As a result, an automorphism maps a vertex of valence m to a vertex of valence m. Whole graph automorphisms applied to asymmetric graphs, similar to the yeast PPI network shown in Figure 5, detect core symmetric regions.
β_{ G } in this case is 1.02693.
As we will see later, this may be applicable to molecular interactomes. A full molecular interactome (not just PPI) is a directed graph, and describes an underlying dynamical system in terms of ordinary differential equations: dx_{ i } /dt = f_{ i } (A,x_{ j } ) where x_{ i } is the state of molecular species i, and A is the full interactome adjacency matrix, an asymmetric matrix. Golubitsky and Stewart [41] point out that the symmetry groups determine the dynamics of the network. When the symmetry changes in one or more factors of the automorphism group, because of a protein mutation or misfolding, for example, this will affect the overall symmetry and thus the dynamics. A catalog of the automorphism groups for interactomes is thus a list of the dynamic behaviors allowed. It might be possible to map these automorphism group elements to disease states. Incidentally, a neural network technique to perform automorphism partitioning is described in Jain and Wysotzki [42].
A proof for this theorem is given for example in Godsil and Royle [34].
We can use this theorem to examine the effects of removing a vertex. If we let e = uv be an edge of G, then the graph G\e is obtained by deleting the edge e from G. The existing PPI network is an extreme case in which a set of unknown edges E and unknown vertices V have been removed from the actual interactome to give us the observed graph P = G \(E,V).
It would be interesting to see how far these deletion theorems can be extended as one approaches graphs with current density. One should be able to test these new theorems empirically with real world data from a manufacturing plant, say an integrated circuit fab. One could start with the full manufacturome and begin deleting edges or vertices and evaluating the theorems observing the effects on the automorphism groups. We know the full interactome should be a directed graph. With the manufacturome, which is of course a directed graph, it should be possible to evaluate and extend other algebraic graph theorems to directed and undirected graphs.
This suggests it may be possible to asymptotically approach an optimal embedding for nteractomes by a gradient descent algorithm to minimize the energy of the embedding.
A number of questions then arise such as: What is the biological, or evolutionary, significance of the embedding space? How does it relate to the automorphism group and the actual molecular network dynamics? Are patterns noticeable for disease trajectories in this higher-dimensional space, or even simple cell cycle trajectories in this space? Are there routes from differentiated cells to pluripotent states? Are there noticeable automorphism group differences between normal cells and polyploidy cells? Is there an isomorphism between the automorphism group and the motifs of Alon [43], and an isomorphism between the order of the automorphism group | Aut(G) | and the average degree distribution < k > or other network statistics? These are all open research questions and some methods described below may be applicable to efforts aimed at answering these questions.
6. Network Dynamics and the Groupoid Formalism
In the above section we described group theory formalism applied to graphs. Here we step up in symmetry, and describe another algebraic object, groupoids; this will allow us to bring more dynamics into the study [41, 44, 45]. Obviously this has importance for understanding the dynamics of molecular interactome networks.
Recall that a directed graph encodes the dynamics given by dx_{ i } /dt = f_{ i } (A,x_{ j } ) where x_{ i } is the state of molecular species i, and A_{ ij } is the full interactome adjacency matrix. More precisely the automorphism group of the network implicitly encodes the dynamics. Further, we know that interactome-like network graphs are composed of multiple copies of a few basic components, e.g. . Groupoids are algebraic objects that resemble groups but the conventional group operation is undefined. In other words, we recognize symmetry but automorphisms are nontrivial. This formalism will allow us to apply group-theory methods to network graphs. Most of this will be focused on small subnets within the larger interactome, where we observe permutation-type automorphisms.
The notion of groupoid is most transparent if we approach it from a categorical angle [46]. The standard definition of a category C involves a collection of objects, A, B,..., and a set of morphisms (which could be the empty set) for each pair of objects; Hom(A, B) for objects A and B. The composition of morphisms is defined and is associative, and there is an identity element in each Hom(A, A), therefore Hom(A, A) is never empty.
But a category C can be viewed as an algebraic structure in itself, endowed with a binary operation, making it similar to a group or semigroup. We call this associated algebraic structure G(C). The "elements" (since the collection of objects do not necessarily form a set) of G(C) are the morphism of C, and the "product" is the composition, which is an associative partial binary operation with identity elements. If C has only one object, then any two morphisms can be composed, and we have only one identity element. The axioms of a category guarantee that G(C) is a semigroup. Furthermore, if we insist on the invertibility of each morphism in C, then G(C) is a group.
Now it is natural to extend the notion of a group by requiring that the objects of C form a set, i.e., C is a small category, and also ask that each morphism of C is invertible. This is the categorical definition of a groupoid. It is easy to translate this definition into the algebraic language, and get a notion similar to the definition of a group [47]. But perhaps it is the categorical definition that illuminates the power of groupoids. Namely, while groups are ideally suited to describe the symmetries of an object, groupoids can similarly capture the symmetries of collections of objects. This is perfectly illustrated in modern algebraic geometry, when one tries to form classifying space, known as moduli spaces, but the algebraic varieties one wants to classify (say elliptic curves) have different symmetries. This problem is solved using the language of stacks and groupoids [48]. The necessity for the same powerful generalization arises in string theory, where symmetries of the physical theory cannot be mathematically realized in terms of topological spaces and groups, only in terms of stacks and groupoids [49].
In the groupoid approach we will examine not the symmetry of the small subnetworks and motifs, but rather the dynamics of these small networks, when they are directed graphs, and in particular when these small nets are wired together to make larger networks (circuits). The symmetries we will observe are not the network symmetry but the symmetries in the phase space or the space of the dynamics.
The interactome, and indeed the full chemical reaction network comprising a cell is a complicated network with numerous feedback loops and feedforward circuits. Its dynamics is no doubt complicated and the details of the full network are only now being elucidated; but we can begin to speculate on some of the possible dynamics by exploiting work from a slightly more mature field--neuronal networks.
We know that biological neuronal nets comprise two- and three-dimensional arrays of frequency-controlled oscillators, voltage controlled oscillators, and logic gates. Engineers have constructed random and non-random networks of these components and discovered not only that the network is capable of memory storage in the form of dynamic patterns and limit cycles (for example memorizing a Bach minuet) but initially random pulse patterns coursing through the network will, after a time delay for component integration, entrain other components and produce continual limit cycles. In large arrays of these networks the limit cycles interact with each other to produce emergent dynamics. In the following we draw on work of Hasslacher and Tilden [50], Rietman et al. [51] and Rietman and Hillis [52]. We argue that by analogy similar dynamics would occur in molecular interaction network of the cell.
Using these dynamics a four-node ring, for example, can exhibit the following three states: (0000), (0101), (0001). Here we employ a permutation-like notation, where, for example (0001) → (0010) → (0100) → (1000) are equivalent to (0001). (Known as the 1 equivalent class, where the underscore is to remind us that this is not a number but a groupoid.) Interestingly, the three states shown here are isomorphic to the nuclei-fusion group: (0000) → dead cell; (0001) → normal healthy cell; (0101) → G1/G2 (equivalence class 5).
One can debate whether or not this is a good model of the cell cycle, but feedforward nets of similar central pattern generators (CPGs) are able to rapidly adapt to changing external stimuli to maintain some entrainment or global stability [53], and from a molecular perspective this is exactly what is required of biological cells. The molecular network in living cells consists of a highly complex interconnected feedback and feedforward chemical reaction system. Walhout and colleagues [54, 55] and others [56, 57] have been discerning some of these details. They have found that feedback and ring circuits, often with inhibitory connections, are common in transcription regulatory networks (protein-DNA interaction networks). One could envision that the basic cell cycle is the primary limit cycle in the dynamics of the cell and the transcription regulator dynamics are used to control and simultaneously be controlled by the cell cycle.
where d_{ i } are the divisors of n with d_{1} = 1, d_{v(n)}= n ; v(n) is the number of divisors of n; ϕ(n) is the totient function, and F(.) is the Fibonacci sequence (where F_{n} = (F_{n-1}) + (F_{n-2}) [58]). The totient function, also called the Euler totient function, is the number of positive integers less than n which are relatively prime to n[51].
Consequently, even small rings of only a dozen nodes can maintain a large number of stable states. Coupling these motifs into networks can produce overall global stability. As Golubitsky and Stewart [41, 44] point out--and as is apparent in the large network of Figure 5--the overall network has very low global symmetry.
where the numbers are a decimal representation of the bit string; they are underlined to remind us that these are group symbols, and are not to be manipulated as numbers. This string of elements interspersed with a rotation operation represents the elements for the group and the main operation. We represent this group by where the superscript reminds us that the group is for six-node rings and the subscript is the lowest decimal equivalent of the bit string in this group.
The above set of mappings shows cyclic permutations from rotation operations on the individual states s represented as decimal equivalent. The and groups are said to be of order 6. The group is of third order and the group is of second order. The similarities between group theory and conventional dynamics are now obvious. The two 6 order groups are 6-cycles. The one third-order group is a three-cycle and the second-order group is a two-cycle.
As the number of rotations needed to return to the starting state decreases for a given group, the periodicity increases--e.g. a two-cycle is faster than a 6-cycle. Similarly, as the number of rotations needed to return to the starting state decreases, the order of the group decreases and the symmetry increases. As we point out later, a symmetry phase transition occurs during signal input and ring coupling.
where ρ is the operator that adds two elements in the group a, b and then applies the modulus operation. The identity element of C_{6} is 0, and the inverse of each element a is b = 6 - a.
The C_{6} group table.
0 | 1 | 2 | 3 | 4 | 5 | |
---|---|---|---|---|---|---|
0 | 0 | 1 | 2 | 3 | 4 | 5 |
1 | 1 | 2 | 3 | 4 | 5 | 0 |
2 | 2 | 3 | 4 | 5 | 0 | 1 |
3 | 3 | 4 | 5 | 0 | 1 | 2 |
4 | 4 | 5 | 0 | 1 | 2 | 3 |
5 | 5 | 0 | 1 | 2 | 3 | 4 |
where k | p means k divides p; φ(k) is the totient function (as discussed above) and Z is the set of integers.
This maps the CPG group to the first nonnegative integers in the cyclic group C_{6} .
There are four subgroups in , and there are four subgroups in C_{6} :{(0),(0,3),(0,2,4),(0,1,2,3,4,5)}.
In order to use these ideas with concepts such as signal input (known as sensor fusion in the control community) and network (ring) coupling to make larger networks, we need to define operators, Φ _{r} , that transform one group into another group. Let the subscript on the operator represent the number of rotations when the signal is injected. Then we can write all of the allowed operations on the groups and their results.
To consider sensor input and/or coupling to two or more of these dynamic rings we consider the example . This equation says that when the CPG circuit has one, 1 cycling through the ring and if a pulse of duration equal to the time constant of the nodes is injected at rotation 2 (subscript to operator), this will be the equivalent of initializing the ring circuit with (000101) or decimal 5. Hence, the circuit is transformed to the group. Explicitly this would be written as (000100) + (000001) → (000101).
These equations are interpreted as follows. From Figure 8 we see a ring in state (000001). If we inject a pulse of short duration (i.e., less then the response time of the logic gates with the associated components) into that ring at position 0 while in that configuration, it will have no effect, . If injected into position 1 while the ring is in this (000001) state it will force the system to transition to state (000101), following the operation . If we inject a short pulse into the network at position 2 it will also transition to (000101) . On the other hand, a short pulse injected at position 3 will cause the ring circuit to exhibit the stable state (001001), according to the operation . In this case, the subscript on the operator indicates the node distance from node 0 in state 1 (000001), while the superscript and subscript on the symbol remind us that the ring is a 6-node ring and it is in state (001001). These transition rules apply for either injected pulses, such as from external sensors, or for internal pulses, such as from rings coupled to make larger networks. The number of states the ring can sustain is still dictated by the ring size as given by the above 2-ary necklace function.
7. Cellular Dynamics Models via Graph Morphisms
Our interpretation of the protein-protein interaction (PPI) network, shown in Figure 5, needs to be considered carefully. The first problem it represents is a biophysical interaction of two proteins as observed in yeast 2-hybrid experiments [59, 60]. These biophysical interactions do not necessarily occur in the actual organism. Second, the PPI networks for most organisms represent only about 10% of the actual possible protein-protein interactions. Third, it is a static network, or time invariant, which is an almost meaningless concept for life forms. We also know that to include the catalytic set for self-replication, the full interactome should include small molecules, large biopolymers, DNA, RNA, oligobiopolymers, etc.
Given these caveats, we will now proceed to parse the PPI in time. We can do this by conducting a relational join between transcription data and PPI data. We start with the expression data as a function of time. Several expression data sets exist; here we mention only the more recent ones by Pramila et al. [61] and Granovskaia et al. [62]. Both of these teams conducted experiments collecting transcription microarray data at five-minute intervals for the yeast S. cerevisiae. The Pramila data (accession number GSE4987) was from cDNA-spotted arrays, and therefore consists of data in the range (-2,2), where zero represents not expressed, below zero represents down regulated, and above zero, up regulated, respectively. The Granovskaia data consist of Affymetrix RNA data (PN 520055) and the numerical data are in the range (-3,2), where data above zero are considered expressed and those below zero not expressed, while the discrimination between up regulated and down regulated is not provided.
The Granovskaia data set is described in their technical paper [62]. They distribute, via links, both the full set of expression data for 6378 gene IDs and a parsed set consisting of 588 genes associated with the cell cycle which clearly show oscillations. Here we ask: What are the large-scale protein-protein interactome changes as a function of time during the cell cycle?
Before we address this question, we note that both the Pramila et al. and the Granovskaia papers show heat maps for the major several hundred genes expressed during the cell cycle. These heat maps show periodic structure and represent periodicity in the transcriptome. Lastly, a paper by de Lichtenbert et al. [63] examines the yeast cell cycle with particular emphasis on parsing the proteome into molecular machines during the cell cycle. Our method differs, as our emphasis is on graph morphisms.
This figure shows for the first time some of the interactome details as a function of time. Each graph represents the changes in the interactome, as represented by the transcriptome, during the indicated time period. Each time period is a 5-minute segment. The red nodes represent those proteins whose expression has disappeared in the time period, and the blue nodes represent those proteins whose expression has appeared in that time period. In a later publication we will be analyzing these graphs in more detail, along with the graph statistical metrics and the automorphism group.
where ||A|| represents the norm of the sum of the two matrices and γ is a Lagrange multiplier called the regularization coefficient. The first term on the RHS is the least mean square of the difference between the target, T, and the learning machine response, R. This regularization technique effectively forces the values for the adjustable parameters in the nonlinear fit, the weight matrices, to very small numerical values, often near zero. Their magnitude is proportional to the regularization coefficient.
An intuitive argument for this regularization may be found in the analogy of fitting 40 data points to a 6000-order polynomial in 2-space. With 6000 adjustable parameters and using a conventional polynomial-fitting algorithm, the plot of the function with the fitted data points would show wild oscillations in the function, with every data point perfectly intercepted by the function. If we fit the same 40 data points to a third-order polynomial, we would find that many of the points were not intercepted by the curve, and there would be an error associated with the fitting. But comparing interpolation on this third-order polynomial and the 6000-dimensional polynomial, we find that our interpolated errors are much lower and the interpolation is more reliable. Now if we again fit our 40 data points to a 6000-dimensional polynomial, but we also force the magnitude of any of the coefficients to be very small, the net effect will resemble a low-order polynomial. There will be an error associated with the fitting, and a much lower error associated with the interpolation. The regularization algorithm does much the same thing; it forces the magnitude of the weight terms to be small, even very small [66, 67].
Table of the gene IDs with error greater than 2 standard deviations (2 × 0
gene ID | seq. no. | mse |
---|---|---|
YKL164C | 2 | 0.5977 |
YNL327W | 5 | 1.0649 |
YNR067C | 8 | 0.6106 |
YOR264W | 10 | 0.4415 |
YKR077W | 36 | 0.3798 |
YDR146C | 397 | 0.2492 |
YGR108W | 399 | 0.5018 |
YMR001CA | 412 | 0.4618 |
YJL051W | 420 | 0.5781 |
YHL028W | 433 | 0.7127 |
YOR049C | 524 | 0.3274 |
YPL158C | 571 | 0.4270 |
YDL179W | 575 | 0.4269 |
YOL101C | 587 | 0.2593 |
With knowledge of the sensitivity analysis we can plot a Pareto chart showing the importance of each of the individual inputs with respect to the output. One could imagine also conducting multi-way digital knockout experiments with this system and comparing it with known experimental results.
Conclusions
In this review we have touched on a few mathematical ideas that may expand our understanding of the boundary between living and non-living systems. We recognize that there are other important works, including category theory [2, 69], genetic networks [70], complexity theory and self-organization [20, 69–71], autopoiesis [72], Turing machines and information theory [73], and many others. It would take a full-length book to review the many subjects that already come into play in discussing the boundaries between living and nonliving. Here we focused on mostly group theory and abstract algebra applied to molecular systems biology. Throughout this paper we have briefly described possible open problems. Here we collect them with respect to the subsections of the paper.
In the section on the genetic code we proposed that it may be possible to use perturbation theory to explore the adjacent possibilities in the 65-dimensional space-time manifold of the evolving genome. One could start by using phylogenetic mappings as historical data on this manifold and compute distances in this space. The statistics of these distances may then be fed back via the perturbation theory to study the trajectory. Of course, we recognize that the existing state-of-the-art bioinformatics makes this proposal mostly unfeasible at this time. But crude outlines of the technique could be developed.
With regards to algebraic graph theory, there are several minor open problems we discussed. Here we reiterate only a few. First, it may be possible to map the automorphism group to disease state through an isomorphism with the phase space. Second, it may be possible to use some of the graph deletion concepts G \(E,V) to evaluate existing protein-protein interactomes. One would start with a known complete graph, say a manufacturing plant, remove edges (vertices) and compute conventional network graph statistics. Third, we pose the following questions. What is the minimal embedding space for something like a protein-protein network? Are there patterns in that space? What biological significance is there to these observed patterns?
In the section on network dynamics and groupoid formalism we suggested that the network graph might not be the main focus for understanding the phenotype but rather the phase space of the network dynamics. We showed a simple case of a C_{6} network and its phase space network. We envision that the molecular network of a cell is actually a complex network of hypercycles and feedback circuits that could be better represented in a higher-dimensional space. Targeting one protein may not have much effect on the overall phenotype or the overall phase space dynamics. For example, a planar array of 6-cycles would give rise to frustration points, not unlike a spin-glass [74]. The overall dynamics would then give rise to not only emergent phase space dynamics but also emergent patterns in the phase space that would not be computable from the molecular reaction network graph [75–79]. Targeting one or two proteins in this network, based on molecular interaction maps, may prove to be futile in many cases. We conjecture that targeting nodes in the molecular network that play key roles in the phase space, as revealed by analysis of the automorphism decomposition might be a better way to carry out drug discovery and treatment of cancer [80].
Declarations
Acknowledgements
We thank Arthur B. Pardee for discussions involving the multi-nucleated cell group table. We thank David Hill and Michael Cusick for helpful discussions. We thank Patrick Rourke for improving the readability at an early stage in the manuscript. EAR thanks Marc Vidal for the free time to pursue this research. We thank Philip Winter for his assistance.
Authors’ Affiliations
References
- Schrödinger E: What Is Life?: The Physical Aspect of the Living Cell ; with, Mind and Matter ; & Autobiographical Sketches. 1992, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Rosen R: Life Itself: A Comprehensive Inquiry into the Nature, Origin, and Fabrication of Life. 1991, New York: Columbia University PressGoogle Scholar
- Goh K-I, Cusick ME, Valle D, Childs B, Vidal M, Barabási A-L: The human disease network. Proc Natl Acad Sci USA. 2007, 104: 8685-8690. 10.1073/pnas.0701361104.PubMed CentralView ArticlePubMedGoogle Scholar
- Annila A: Space, time and machines. 2009, arXiv:0910.2629v1 [physics.gen-ph]Google Scholar
- Penrose R: The Road to Reality: A Complete Guide to the Laws of the Universe. 2004, London: Jonathan CapeGoogle Scholar
- Weinberg S: The Quantum Theory of Fields. 1995, Cambridge: Cambridge University Press, 1:View ArticleGoogle Scholar
- Kaku M: Quantum Field Theory: A Modern Introduction. 1993, New York: Oxford University PressGoogle Scholar
- Rosen J: Symmetry Rules: How Science and Nature Are Founded on Symmetry. 2008, New York: Springer, 1View ArticleGoogle Scholar
- Brooks R: The relationship between matter and life. Nature. 2001, 409: 409-411. 10.1038/35053196.View ArticlePubMedGoogle Scholar
- Woese CR: A New Biology for a New Century. Microbiol Mol Biol Rev. 2004, 68: 173-186. 10.1128/MMBR.68.2.173-186.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Cohen JE: Mathematics Is Biology's Next Microscope, Only Better; Biology Is Mathematics' Next Physics, Only Better. PLoS Biol. 2004, 2: e439-10.1371/journal.pbio.0020439.PubMed CentralView ArticlePubMedGoogle Scholar
- Goodman FM: Algebra: Abstract and Concrete. 2006, Iowa City, IA: SemiSimple Press, 2.5Google Scholar
- Woolfson A: Life Without Genes. 2000, London: HarperCollinsGoogle Scholar
- Crick FHC: Codon--anticodon pairing: The wobble hypothesis. J Mol Biol. 1966, 19: 548-555. 10.1016/S0022-2836(66)80022-0.View ArticlePubMedGoogle Scholar
- Danckwerts HJ, Neubert D: Symmetries of genetic code-doublets. J Mol Evol. 1975, 5: 327-332. 10.1007/BF01732219.View ArticlePubMedGoogle Scholar
- Bertman MO, Jungck JR: Group graph of the genetic code. J Hered. 1979, 70: 379-384.PubMedGoogle Scholar
- Findley AM, Findley GL, McGlynn SP: Genetic coding: approaches to theory construction. J Theor Biol. 1982, 97: 299-318. 10.1016/0022-5193(82)90106-0.View ArticlePubMedGoogle Scholar
- Findley AM, McGlynn SP, Findley GL: Applications of differential geometry to molecular genetics. J Biol Phys. 1985, 13: 87-94. 10.1007/BF01878385.View ArticleGoogle Scholar
- Kauffman SA: At Home in the Universe: The Search for Laws of Self-Organization and Complexity. 1995, New York: Oxford University PressGoogle Scholar
- Kauffman SA: Investigations. 2000, Oxford: Oxford University PressGoogle Scholar
- Kauffman S: Molecular autonomous agents. Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences. 2003, 361: 1089-1099. 10.1098/rsta.2003.1186.View ArticleGoogle Scholar
- Hornos JEM, Hornos YMM: Algebraic model for the evolution of the genetic code. Phys Rev Lett. 1993, 71: 4401-10.1103/PhysRevLett.71.4401.View ArticlePubMedGoogle Scholar
- Forger M, Hornos YMM, Hornos JEM: Global aspects in the algebraic approach to the genetic code. Phys Rev E. 1997, 56: 7078-10.1103/PhysRevE.56.7078.View ArticleGoogle Scholar
- Stone M, Goldbart PM: Mathematics for Physics: A Guided Tour for Graduate Students. 2009, Cambridge, UK: Cambridge University PressView ArticleGoogle Scholar
- Gallian JA: Contemporary Abstract Algebra. 2006, Boston, MA: Houghton Mifflin, 6Google Scholar
- Rao PN, Johnson RT: Mammalian Cell Fusion : Studies on the Regulation of DNA Synthesis and Mitosis. Nature. 1970, 225: 159-164. 10.1038/225159a0.View ArticlePubMedGoogle Scholar
- Johnson RT, Rao PN: Mammalian Cell Fusion : Induction of Premature Chromosome Condensation in Interphase Nuclei. Nature. 1970, 226: 717-722. 10.1038/226717a0.View ArticlePubMedGoogle Scholar
- Hanna JH, Saha K, Jaenisch R: Pluripotency and Cellular Reprogramming: Facts, Hypotheses, Unresolved Issues. Cell. 2010, 143: 508-525. 10.1016/j.cell.2010.10.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim J, Eberwine J: RNA: state memory and mediator of cellular phenotype. Trends Cell Biol. 2010, 20: 311-318. 10.1016/j.tcb.2010.03.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Junker BH, Schreiber F: Analysis of Biological Networks. 2008, Hoboken, N.J: Wiley-InterscienceView ArticleGoogle Scholar
- Mason O, Verwoerd M: Graph theory and networks in Biology. IET Syst Biol. 2007, 1: 89-119. 10.1049/iet-syb:20060038.View ArticlePubMedGoogle Scholar
- Newman MEJ, Barabási A-L, Watts DJ: The Structure and Dynamics of Networks. 2006, Princeton: Princeton University PressGoogle Scholar
- Albert R, Barabási A-L: Statistical mechanics of complex networks. Rev Mod Phys. 2002, 74: 47-97. 10.1103/RevModPhys.74.47.View ArticleGoogle Scholar
- Godsil CD, Royle G: Algebraic Graph Theory. 2001, New York: SpringerView ArticleGoogle Scholar
- Chung FRK: Spectral Graph Theory. 1997, Providence, R.I: Published for the Conference Board of the mathematical sciences by the American Mathematical SocietyGoogle Scholar
- MacArthur BD, Anderson JW: Symmetry and Self-Organization in Complex Systems. 2006, arXiv:cond-mat/0609274v1 [cond-mat.dis-nn]Google Scholar
- Xiao Y, MacArthur BD, Wang H, Xiong M, Wang W: Network quotients: Structural skeletons of complex systems. Phys Rev E. 2008, 78: 046102-View ArticleGoogle Scholar
- MacArthur BD, Sánchez-García RJ, Anderson JW: Symmetry in complex networks. Discrete Appl Math. 2008, 156: 3525-3531. 10.1016/j.dam.2008.04.008.View ArticleGoogle Scholar
- McKay BD: Practical Graph Isomorphism. Congressus Numerantium. 1981, 30: 45-87.Google Scholar
- Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, Califano A: Reverse engineering of regulatory networks in human B cells. Nat Genet. 2005, 37: 382-390. 10.1038/ng1532.View ArticlePubMedGoogle Scholar
- Golubitsky M, Stewart I: Nonlinear dynamics of networks: the groupoid formalism. Bulletin of the American Mathematical Society. 2006, 43: 305-364. 10.1090/S0273-0979-06-01108-6.View ArticleGoogle Scholar
- Jain BJ, Wysotzki F: Automorphism Partitioning with Neural Networks. Neural Processing Letters. 2003, 17: 205-215. 10.1023/A:1023657727387.View ArticleGoogle Scholar
- Alon U: An Introduction to Systems Biology: Design Principles of Biological Circuits. 2007, Boca Raton, FL: Chapman & Hall/CRCGoogle Scholar
- Golubitsky M, Stewart I: The Symmetry Perspective: From Equilibrium to Chaos in Phase Space and Physical Space. 2002, Boston, MA: BirkhäuserView ArticleGoogle Scholar
- Brown R: From Groups to Groupoids: a Brief Survey. Bulletin of the London Mathematical Society. 1987, 19: 113-134. 10.1112/blms/19.2.113.View ArticleGoogle Scholar
- Higgins PJ: Notes on Categories and Groupoids. 1971, London: Van Nostrand Reinhold CoGoogle Scholar
- Dicks W, Ventura E: The group fixed by a family of injective endomorphisms of a free group. 1996, Providence, R.I.: American Mathematical SocietyView ArticleGoogle Scholar
- Laumon G, Moret-Bailly L: Champs algebriques. 1999, Berlin: SpringerGoogle Scholar
- Karp RL: Quantum Symmetries and Exceptional Collections. Communications in Mathematical Physics. 2010, 301: 1-21.View ArticleGoogle Scholar
- Hasslacher B, Tilden MW: Living machines. Robotics and Autonomous Systems. 1995, 15: 143-169. 10.1016/0921-8890(95)00019-C.View ArticleGoogle Scholar
- Rietman EA, Tilden MW, Askenazi M: Analog computation with rings of quasiperiodic oscillators: the microdynamics of cognition in living machines. Robotics and Autonomous Systems. 2003, 45: 249-263. 10.1016/j.robot.2003.08.002.View ArticleGoogle Scholar
- Rietman EA, Hillis RW: Neural Computation with Rings of Quasiperiodic Oscillators. 2006, arXiv:cs/0611136v1 [cs.RO]Google Scholar
- Pikovsky A, Rosenblum M, Kurths J: Synchronization: A Universal Concept in Nonlinear Sciences. 2001, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Vermeirssen V, Barrasa MI, Hidalgo CA, Babon JAB, Sequerra R, Doucette-Stamm L, Barabási A-L, Walhout AJM: Transcription factor modularity in a gene-centered C. elegans core neuronal protein-DNA interaction network. Genome Res. 2007, 17: 1061-1071. 10.1101/gr.6148107.PubMed CentralView ArticlePubMedGoogle Scholar
- Arda HE, Taubert S, MacNeil LT, Conine CC, Tsuda B, Van Gilst M, Sequerra R, Doucette-Stamm L, Yamamoto KR, Walhout AJM: Functional modularity of nuclear hormone receptors in a Caenorhabditis elegans metabolic gene regulatory network. Mol Syst Biol. 2010, 6: 367-PubMed CentralView ArticlePubMedGoogle Scholar
- Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, Zeitlinger J, Jennings EG, Murray HL, Gordon DB, Ren B, Wyrick JJ, Tagne J-B, Volkert TL, Fraenkel E, Gifford DK, Young RA: Transcriptional Regulatory Networks in Saccharomyces cerevisiae. Science. 2002, 298: 799-804. 10.1126/science.1075090.View ArticlePubMedGoogle Scholar
- Xu Z, Wei W, Gagneur J, Perocchi F, Clauder-Munster S, Camblong J, Guffanti E, Stutz F, Huber W, Steinmetz LM: Bidirectional promoters generate pervasive transcription in yeast. Nature. 2009, 457: 1033-1037. 10.1038/nature07728.PubMed CentralView ArticlePubMedGoogle Scholar
- Wolfram MathWorld: The Web's Most Extensive Mathematics Resource.http://mathworld.wolfram.com/http://mathworld.wolfram.com/
- Bartel PL, Fields S: The Yeast Two-Hybrid System. 1997, New York: Oxford University PressGoogle Scholar
- Dreze M, Monachello D, Lurin C, Cusick ME, Hill DE, Vidal M, Braun P: High-Quality Binary Interactome Mapping. Methods in Enzymology. 2010, Academic Press, 470: 281-315.Google Scholar
- Pramila T, Wu W, Miles S, Noble WS, Breeden LL: The Forkhead transcription factor Hcm1 regulates chromosome segregation genes and fills the S-phase gap in the transcriptional circuitry of the cell cycle. Genes Dev. 2006, 20: 2266-2278. 10.1101/gad.1450606.PubMed CentralView ArticlePubMedGoogle Scholar
- Granovskaia MV, Jensen LJ, Ritchie ME, Toedling J, Ning Y, Bork P, Huber W, Steinmetz LM: High-resolution transcription atlas of the mitotic cell cycle in budding yeast. Genome Biol. 2010, 11: R24-R24. 10.1186/gb-2010-11-3-r24.PubMed CentralView ArticlePubMedGoogle Scholar
- de Lichtenberg U, Jensen LJ, Brunak S, Bork P: Dynamic Complex Formation During the Yeast Cell Cycle. Science. 2005, 307: 724-727. 10.1126/science.1105103.View ArticlePubMedGoogle Scholar
- Gershenfeld NA: The Nature of Mathematical Modeling. 1999, Cambridge: Cambridge University PressGoogle Scholar
- Rumelhart DE, McClelland JL, University of California, San Diego: Parallel Distributed Processing: Explorations in the Microstructure of Cognition. 1986, Cambridge, Mass: MIT Press, 1:Google Scholar
- Vapnik VN: Statistical Learning Theory. 1998, New York: WileyGoogle Scholar
- Vovk V, Gammerman A, Shafer G: Algorithmic Learning in a Random World. 2005, New York: SpringerGoogle Scholar
- Shawe-Taylor J, Cristianini N: Kernel Methods for Pattern Analysis. 2004, Cambridge, UK: Cambridge University PressView ArticleGoogle Scholar
- Kampis G: Self-Modifying Systems in Biology and Cognitive Science: A New Framework for Dynamics, Information, and Complexity. 1991, Oxford: Pergamon Press, 1Google Scholar
- Kauffman SA: The Origins of Order: Self-Organization and Selection in Evolution. 1993, New York: Oxford University PressGoogle Scholar
- Feistel R, Ebeling W: Evolution of Complex Systems: Self-Organization, Entropy, and Development. 1989, Dordrecht, Holland: Kluwer Academic PublishersGoogle Scholar
- Maturana HR, Varela FJ: Autopoiesis and Cognition: The Realization of the Living. 1980, Boston, MA: Kluwer Academic PublishersView ArticleGoogle Scholar
- Chaitin GJ: Metabiology. 2010,http://www.umcs.maine.edu/~chaitin/http://www.umcs.maine.edu/~chaitin/Google Scholar
- Fischer KH, Hertz J: Spin Glasses. 1991, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Anderson PW: More Is Different. Science. 1972, 177: 393-396. 10.1126/science.177.4047.393.View ArticlePubMedGoogle Scholar
- Wolfram S: Undecidability and intractability in theoretical physics. Phys Rev Lett. 1985, 54: 735-738. 10.1103/PhysRevLett.54.735.View ArticlePubMedGoogle Scholar
- Moore C: Unpredictability and undecidability in dynamical systems. Phys Rev Lett. 1990, 64: 2354-2357. 10.1103/PhysRevLett.64.2354.View ArticlePubMedGoogle Scholar
- Laughlin RB, Pines D: The Theory of Everything. Proc Natl Acad Sci USA. 2000, 97: 28-31. 10.1073/pnas.97.1.28.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu M, Weedbrook C, Perales Á, Nielsen MA: More really is different. Physica D. 2009, 238: 835-839. 10.1016/j.physd.2008.12.016.View ArticleGoogle Scholar
- Huang S, Eichler G, Bar-Yam Y, Ingber DE: Cell Fates as High-Dimensional Attractor States of a Complex Gene Regulatory Network. Phys Rev Lett. 2005, 94: 128701-View ArticlePubMedGoogle Scholar
- Yu H, Braun P, Yıldırım MA, Lemmens I, Venkatesan K, Sahalie J, Hirozane-Kishikawa T, Gebreab F, Li N, Simonis N, Hao T, Rual J-F, Dricot A, Vazquez A, Murray RR, Simon C, Tardivo L, Tam S, Svrzikapa N, Fan C, de Smet A-S, Motyl A, Hudson ME, Park J, Xin X, Cusick ME, Moore T, Boone C, Snyder M, Roth FP, Barabási A-L, Tavernier J, Hill DE, Vidal M: High-Quality Binary Protein Interaction Map of the Yeast Interactome Network. Science. 2008, 322: 104-110. 10.1126/science.1158684.PubMed CentralView ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.