This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Comparing biological models and theories of life with process-enablement graphs

Emmy Browna and Sean T. Vittadelloa,b,∗
aSchool of Mathematics and Statistics, The University of Melbourne,
Parkville VIC 3010, Australia
bSchool of BioSciences, The University of Melbourne,
Parkville VIC 3010, Australia
Corresponding author: [email protected]
Abstract

There are many perspectives through which biologists can study a particular living system. As a result, models of biological systems are often quite different from one another, both in form and size. Thus, in order for us to generate reliable knowledge of a particular system, we need to understand how the models that represent it are related. In previous work, we constructed a general model comparison framework to compare models representing any physical system. Here, we develop an alternative methodology that focuses on a fundamental feature of living systems, namely self-organisation. We employ a graph theoretic formalism which captures self-organising processes as cycles within particular kinds of graphs: process-enablement graphs. We then build the mathematical tools needed to compare biological models and their corresponding descriptions of self-organisation in a consistent and rigorous manner. We apply our formalism to a range of classical theories of life to show how they are similar and where they differ. We also investigate examples of putatively abiotic systems which nonetheless still realise primitive forms of self-organisation. While our current framework does not demarcate living systems from nonliving ones, it does allow us to better study the grey area surrounding life’s edge.

Keywords: Self-organisation, model comparison, perspectival realism, autopoiesis, (M,R)(M,R)-systems, autocatalytic sets, constraints, viruses, graph theory, organisational closure.

1 Introduction

Several attempts have been made throughout the past century to construct theories of life, including (M,R)(M,R)-systems Rosen (\APACyear1991), autopoiesis Maturana \BBA Varela (\APACyear1980), the chemoton Gánti (\APACyear2003), autocatalytic sets Kauffman (\APACyear1986), and constraint closure Montévil \BBA Mossio (\APACyear2015). These accounts are each slightly different from one another, and each emphasises different characteristics of living systems. However, in one way or another they all explore the concept of organisational closure, which is the idea that biological systems have to continuously produce and maintain the conditions of their own existence in order to stay alive. That is, organisms are self-determining Montévil \BBA Mossio (\APACyear2015). The precise way in which authors articulate organisational closure varies greatly, which is hardly surprising given the complexity of the biosphere, but it leaves the path towards developing future theories of life unclear. Some have suggested that we should simply incorporate all of the best features from our current models into a hybrid theory Cornish-Bowden \BBA Cárdenas (\APACyear2020). Others say that it will be nearly impossible to construct a general theory of life without first finding extraterrestrial life Cleland (\APACyear2019).

Desire for a synthesis of these different theories of life is understandable. Modern biology seems more fragmented than ever before, so the allure of unification is very tempting Vittadello \BBA Stumpf (\APACyear2022\APACexlab\BCnt2); van Hemmen (\APACyear2007). How nice it would be to have universal laws of biology! However, many biologists and philosophers alike have argued that expecting a synthesis or a grand unified theory of such varied accounts of life is unreasonable due to life’s complexity Rosen (\APACyear1991); Mitchell (\APACyear2002); Wimsatt (\APACyear2007); Gyllingberg \BOthers. (\APACyear2023). The more promising way forward is a pluralistic approach, and in particular, one aligned with the philosophy of perspectival realism. This philosophical view, as developed by \citeAmassimi2022perspectival, argues that we should embrace the many ways in which we can model a (biological) system and acknowledge that any single approach on its own is necessarily limited and imperfect. This is not to say that any given model of life is worthy of consideration, a position known as flabby pluralism Lapid (\APACyear2003). There are indeed model-independent lawlike dependencies in nature that ultimately act as the adjudicators of scientific knowledge claims Massimi (\APACyear2023). However, there is no objective vantage point, or view from nowhere, from which we can see these lawlike dependencies. We can only come to understand reality through the lens of particular scientific models, therefore from somewhere Massimi (\APACyear2018).

But it is not enough to simply create many different models of a biological system. Scientific perspectives need to intersect and interlace with one another in order for us to model reality more accurately. Thus, reliable scientific knowledge is obtained when epistemic (knowledge-producing) communities enter into dialogue with another and ask questions like: What phenomena does your model describe that ours does not? What explanation does your model give for phenomenon X? Is that a complementary or contradictory explanation to ours? Dataset Y fits our model, is it compatible with yours? Our model predicts Z, does yours predict the same? The goal here is not to find a winner or to arrive at a ‘true’ theory. This is because new results and new assumptions down the line can revive models that are currently inconsistent with observations at the moment Mitchell \BBA Dietrich (\APACyear2006). For instance, Jean-Baptiste Lamarck’s theory of inheritance of acquired characteristics was long dismissed throughout the 20th century, but it has now been partially vindicated through the modern study of transgenerational epigenetic inheritance Jablonka \BBA Raz (\APACyear2009); McGuigan \BOthers. (\APACyear2021); Ashe \BOthers. (\APACyear2021). So instead, we should keep all useful perspectives around, even if they are incompatible Chang (\APACyear2012). The coexistence of such perspectives can allow for insightful interactions as time goes on, if we commit to remaining in dialogue with one another.

Modern authors have already begun to compare and contrast the different accounts of life with each other in order to develop new ideas Montévil \BBA Mossio (\APACyear2015); Hofmeyr (\APACyear2021); Kauffman (\APACyear2019). It is this spirit of productive cross-dialogue that we wish to further cultivate here with the help of a mathematical formalism. In \citeAvittadello2021model, vittadello2022group, the authors provide a formalism for model comparison using the language of simplicial complexes, valid also for general hypergraphs. Components of each model are represented as labelled vertices and interconnections between model components are represented as higher-dimensional simplices. Once models are expressed in this abstract framework they can be compared using two methods, comparison by distance and comparison by equivalence. In this paper, we introduce a framework to compare theories of life based on key characteristics of living organisms, motivated by the approach to model comparison in \citeAvittadello2021model with two significant modifications.

First, we only consider the processes described by each model and disregard other features. The value of viewing the living world as being made up of processes, not substances, has become increasingly recognised among philosophers of biology Nicholson \BBA Dupré (\APACyear2018); Meincke (\APACyear2022). Organisms are constantly taking in new nutrients, breaking them down, and excreting the waste products back into the environment. This metabolic turnover never ceases within an organism, or as Heraclitus put it, everything flows. From ecosystems all the way down to the macromolecules that make up cells, life is a hierarchy of processes. Whenever it seems as though we are looking at a static entity, look closer and we will see that it is no more than a temporary stability which is, in fact, actively flowing. This metaphysical thesis is quite strong and we will not provide a detailed argument to defend it here (see \citeAnicholson2018everything for more details). However, at a minimum, it is surely an interesting perspective that can lead to new insights. Indeed, the rest of this paper highlights some of the advantages of taking a process perspective, especially with regard to biological systems.

Second, we shift from simplicial complexes to directed graphs. This allows us to model asymmetrical relationships between processes as opposed to symmetrical interconnections between concepts. In particular, we consider how processes enable one another within a particular model. The notion of enablement has been developed as an important way to understand how the dynamics of the living biosphere differ from the world modelled by physics Longo \BBA Montévil (\APACyear2013); Kauffman (\APACyear2019). Here, we define enablement precisely and show how cycles of processes embody a primitive form of organisational closure. Admittedly, cycles of processes can appear in all sorts of systems, including nonliving ones like candles or the water cycle. Hence, the identification of a cycle of processes is by no means a sufficient condition for identifying life Mossio \BBA Bich (\APACyear2017). However, by studying how these cycles are captured by particular models and by making comparisons using the tools of graph theory, we can explicitly compare how different theories of life articulate self-organisational structures.

More generally, we define a particular type of graph homomorphism that preserves cycles of processes between models, allowing for fine-grained perspectives to be consistently compared to coarse-grained ones. This contextualises the problem of self-organising features seeming to appear or disappear depending on the level of analysis Cusimano \BBA Sterner (\APACyear2020); Nahas \BBA Sachs (\APACyear2023). Accordingly, our framework allows lower-level models of biochemical interactions to be compared to more abstract models of the cell, and we can explicitly see when and how organisational closure is achieved in these systems. Such integrative model comparison can therefore bring vastly different perspectives into contact with one another, enabling the development of more informed theories of life without arbitrarily forcing unification Mitchell \BBA Dietrich (\APACyear2006); Gyllingberg \BOthers. (\APACyear2023).

The remainder of this paper is organised as follows. In Section 2, we develop the concepts of process and (direct) enablement in relation to the biosphere. In Section 3, we then use these concepts to define the central object of study in this paper, namely the process-enablement graph. We also put forward a mathematical conception of organisational closure in terms of cycles in process-enablement graphs, and develop mathematical tools to consistently compare such cycles between graphs. Section 4 explores a range of applications of our formalism and explains how to translate abstract models of living systems into the language of process-enablement graphs. We use each example to highlight features and possible extensions of our mathematical framework. In particular, in Section 4.4.2 we note that our current formalism is static, not dynamic, and will need to be extended appropriately to model evolution, development, and biological agency. Finally, in the Conclusion we comment on how process-enablement graphs allow us to advance biology by grounding the science in the philosophy of perspectival realism.

2 Process and enablement

The concept at the core of this paper is that of a process-enablement graph. We consider a process as an organised sequence of consecutive physicochemical changes that can involve production, destruction, alteration, or movement.

We call a set of processes interacting contemporaneously with each other a system. That is, for any process pp within a given system SS, there exists a distinct process qq in SS that will be altered by the removal of pp. The interacting condition here prevents us from placing two very distant processes together into a single system. For instance, we could not say that a particular tree in the Amazon rainforest together with a particular earthworm in the soil here in Melbourne constitute a single system, since neither one will be affected by the absence of the other. However, both organisms are part of the global network of processes within the system of Earth’s biosphere, and they may interact (albeit very indirectly) through this network. A subset of the processes and interactions occurring within a system is called a subsystem.

Although this is quite abstract, processes at all levels of the biosphere will be familiar to biologists: biochemical reactions, protein folding, transcription, translation, cell division, action potentials, predation, natural selection, inhalation of oxygen, and digestion of food. Examples of systems should also be familiar, such as: cells, organs, organisms, and ecosystems. Systems themselves are also processes as they are organised sequences of consecutive physicochemical changes. But by using the term system we are emphasising that we are thinking of a network of interacting processes. Since it is always possible to view a system in different ways from different perspectives, there is not necessarily a unique way to enumerate all of the processes within a particular system.

Given any process pp, we can associate it with the region of space XX which is occupied by pp. The space XX is generated by pp and is inherently tied to it. That is, as the process pp interacts with its environment, the space XX associated with pp will also change. In any case, this assignment allows us to say that pp is happening within the region XX and is not happening outside of XX. When required, we write pXp_{{}_{X}} for a process pp with corresponding spatial region XX.

To capture the hierarchical nature of processes in biology we introduce the following notion. Let pXp_{{}_{X}} and qYq_{{}_{Y}} be two processes. We say that pXp_{{}_{X}} happens within qYq_{{}_{Y}} if pXp_{{}_{X}} and qYq_{{}_{Y}} are occurring contemporaneously and XX is contained within YY. For example, within an adult human, a single red blood cell passing through a capillary is a process that happens within the broader process of blood circulation. The concept of happening within will actually become quite important for our formalism later on, as it will severely constrain which kinds of graphs can be compared. Indeed, we will only make a comparison between two processes if one is happening within the other. This restricts our framework to studying models of the same system, rather than exploring more abstract comparisons between process-enablement structures in different systems.

We now consider the concept of enablement. Within a system SS, we say that a process pp enables another distinct process qq when pp has provided some of the necessary conditions for qq to occur. This necessity should be understood counterfactually. In other words, if pp had not happened, then qq would not be happening. Note that enablement is defined with respect to the given system SS, so that pp may not be necessary for qq to occur within some other system TT. For example, within the cell, ribosomes are necessary to produce polypeptides from RNA sequences via translation, so the diffusion of ribosomes throughout the cell enables polypeptide synthesis. But in the lab we can make polypeptides synthetically via chemical ligation, without the need for ribosomes, so in this case polypeptide synthesis is not enabled by ribosome diffusion.

If pXp_{{}_{X}} enables qYq_{{}_{Y}} and there also exists a time interval τ\tau over which pXp_{{}_{X}} and qYq_{{}_{Y}} are interacting with each other, then we say that pXp_{{}_{X}} enables qYq_{{}_{Y}} directly and we write this as pXqYp_{{}_{X}}\to q_{{}_{Y}} or pqp\to q. To be more explicit, by ‘interacting’ we mean that XX and YY intersect over τ\tau, and within this spatiotemporal region there is some physical interaction between the two processes. Using the example from the previous paragraph, ribosome production directly enables ribosome diffusion through the cell which in turn directly enables translation. Conversely, in a particular forest, a tree releasing a seed might enable a new plant to grow once it rains. However, this enablement is not direct since there is an intermediate period where the tree is not releasing the seed and the new plant is also not growing. Direct enablements are a very useful concept in biology since the temporal overlap between processes is a key feature of life. It is one of the ways in which biological systems remain tightly organised, instead of simply waiting for the rain to come.

In summary, to investigate whether a direct enablement pXqYp_{{}_{X}}\to q_{{}_{Y}} within a system SS is genuine or not, we need to check two things. First, we need to see whether the hypothetical removal of pXp_{{}_{X}} from SS would result in qYq_{{}_{Y}} not occurring. This ensures that pXp_{{}_{X}} enables qYq_{{}_{Y}}. Second, we need to check that there is a time interval τ\tau over which pXp_{{}_{X}} and qYq_{{}_{Y}} are interacting with each other. This ensures that the enablement is direct.

3 Process-enablement graphs

In this section we introduce and develop our mathematical formalism based on our central notion of a process-enablement graph. We begin with the required background on graph theory. A directed graph G=(V(G),E(G))G=(V(G),E(G)) is a pair of finite sets where V(G)V(G) is a nonempty set of elements called vertices and E(G)E(G) is a set of ordered pairs of elements of V(G)V(G) called (directed) edges. For notational simplicity, we denote an edge (u,v)E(G)(u,v)\in E(G) by uvuv. To visualise a directed graph, we draw the vertices as points and the edges as arrows. That is, the edge uvuv is drawn as uvu\to v. This definition of a directed graph allows for loops which are edges from a vertex to itself, for example uuu\to u, but not parallel edges which are distinct edges between the same two vertices pointed in the same direction, for example uvu\rightrightarrows v. For the remainder of this paper we refer to directed graphs simply as graphs.

A subgraph of GG is a graph HH such that V(H)V(G)V(H)\subseteq V(G) and E(H)E(G)E(H)\subseteq E(G). We say that GG contains HH and write HGH\subseteq G. In particular, if V(H)V(G)V(H)\neq V(G) or E(H)E(G)E(H)\neq E(G) then we say that HH is a proper subgraph of GG. We can obtain a subgraph from a subset of vertices by considering all possible edges that can be drawn between the vertices. Explicitly, if SS is a nonempty subset of V(G)V(G) then the subgraph G[S]G[S] induced by SS has vertex set SS and edge set

E(G[S])={uvE(G)u,vS}.E(G[S])=\{\,uv\in E(G)\mid u,v\in S\,\}.

Moreover, we say that a subgraph HH of GG is an induced subgraph if there exists a nonempty subset AV(G)A\subseteq V(G) such that H=G[A]H=G[A]. A cycle C=(V(C),E(C))C=(V(C),E(C)) is a graph with vertex and edge sets

V(C)\displaystyle V(C) ={v0,,vn1},\displaystyle=\{v_{0},\ldots,v_{n-1}\},
E(C)\displaystyle E(C) ={v0v1,v1v2,,vn1v0}.\displaystyle=\{v_{0}v_{1},v_{1}v_{2},\ldots,v_{n-1}v_{0}\}.

The cycle CC has length nn and a cycle of length 11 is a loop. These concepts are illustrated in Figure 1.

AABBCCDDG1G_{1}
BBCCDDG2G_{2}
AABBCCDDG3G_{3}
Figure 1: G2G_{2} and G3G_{3} are subgraphs of G1G_{1}. G2G_{2} is an induced subgraph of G1G_{1}, but G3G_{3} is not since the edge CBCB is missing. G2G_{2} is also a cycle since V(G2)={B,D,C}V(G_{2})=\{B,D,C\} and E(G2)={BD,DC,CB}E(G_{2})=\{BD,DC,CB\}.

To compare the structures of two graphs GG and HH we can consider a graph homomorphism ϕ:GH\phi\colon G\to H, which is a structure preserving map given by a vertex function ϕ:V(G)V(H)\phi\colon V(G)\to V(H) such that if uvE(G)uv\in E(G) then ϕ(u)ϕ(v)E(H)\phi(u)\phi(v)\in E(H). We call GG the source graph and HH the target graph of ϕ\phi. For a subgraph G0GG_{0}\subseteq G the induced image of G0G_{0} under ϕ\phi is the induced subgraph H[ϕ(V(G0))]H[\phi(V(G_{0}))] of HH, which we denote by H[ϕ(G0)]H[\phi(G_{0})].

For our purposes, we want to compare graphs with slightly less structure preservation using the more general notion of a weak graph homomorphism ϕ:GH\phi\colon G\to H, which is a vertex function ϕ:V(G)V(H)\phi\colon V(G)\to V(H) such that if uvE(G)uv\in E(G) and ϕ(u)ϕ(v)\phi(u)\neq\phi(v) then ϕ(u)ϕ(v)E(H)\phi(u)\phi(v)\in E(H). We use weak graph homomorphisms rather than graph homomorphisms because the former allow us to contract edges into vertices Hell \BBA Nesetril (\APACyear2004). Within the context of our formalism, this feature of weak graph homomorphisms will allow us to zoom out from lower-level models of a system to view the system more coarsely.

Importantly, weak homomorphisms are composable, as shown by the following standard result.

Theorem 3.1.

Let G1G_{1}, G2G_{2}, and G3G_{3} be graphs, and let ϕ:G1G2\phi\colon G_{1}\to G_{2} and ψ:G2G3\psi\colon G_{2}\to G_{3} be weak graph homomorphisms. Then ψϕ:G1G3\psi\circ\phi\colon G_{1}\to G_{3} is a weak graph homomorphism.

Proof.

Let uvE(G1)uv\in E(G_{1}) and assume ψ(ϕ(u))ψ(ϕ(v))\psi(\phi(u))\neq\psi(\phi(v)). Then ϕ(u)ϕ(v)\phi(u)\neq\phi(v) so, since ϕ\phi is a weak graph homomorphism, ϕ(u)ϕ(v)E(G2)\phi(u)\phi(v)\in E(G_{2}). Further, since ψ\psi is a weak graph homomorphism, ψ(ϕ(u))ψ(ϕ(v))E(G3)\psi(\phi(u))\psi(\phi(v))\in E(G_{3}). Therefore, ψϕ\psi\circ\phi is a weak graph homomorphism. ∎

We now define our notion of a process-enablement graph for modelling physical systems from a process perspective.

Definition 3.2 (Process-enablement graph).

A process-enablement graph, or pepe-graph, is a graph GG where V(G)V(G) is a set of processes occurring contemporaneously in a particular system SS over a particular time interval τ\tau, and E(G)E(G) is the set of direct enablements between the processes.

3.1 Organisational closure

For the remainder of this paper we use the term organisational closure to refer to the existence of a cycle in a relevant pepe-graph. This definition of organisational closure is weaker than definitions put forward in other contemporary accounts. In particular, we have defined organisational closure in terms of a closure of processes which differs from \citeAmontevil2015biological’s characterisation of biological organisation as a closure of constraints. There, constraints are conceived of as “contingent causes, exerted by specific structures of dynamics, which reduce the degrees of freedom of the system on which they act.”

Processes and constraints differ from each other in a few important ways. Processes, in general, often contribute energy and matter to other downstream processes. Constraints, on the other hand, merely canalise processes at lower levels of causation, without directly contributing energy or matter Montévil \BBA Mossio (\APACyear2015). This allows constraints to do work, in the strict sense of physics, which we consider here as a constrained release of energy Kauffman (\APACyear2000). Constraints also remain unaltered as they act on systems, when viewed from an appropriate perspective. For example, our vasculature acts as a constraint on the flow of blood through our body since, without blood vessels, blood would flow in all directions instead of being neatly transported to our cells. Further, over short time periods our vasculature remains largely unaltered. But by virtue of existing in a world where everything flows, constraints are themselves processes. Using the present example, over longer timescales our vasculature will change its physical structure through processes like angiogenesis. However, processes are not necessarily constraints, as not all processes reduce the degrees of freedom of a particular system.

It is arguably important to make a distinction between the processes themselves and the constraints that impinge on them Mossio \BBA Bich (\APACyear2017). We have not made a distinction between processes and constraints in our current framework and will therefore not be able to identify where a system is doing work. Further, we also lose the powerful concept of self-constraint in our analysis which has become a fruitful framework to understanding biological organisation Juarrero (\APACyear2023); Moreno \BBA Mossio (\APACyear2015). We make this sacrifice for means of mathematical simplicity, since it is far easier to make mathematical comparisons when there is only one kind of vertex (processes) and one kind of arrow (enablements) as opposed to, for example, the multiple types of arrows used by \citeAmontevil2015biological.

Organisational closure is nevertheless biologically interesting since cycles in pepe-graphs embody a more general form of self-organisation. To see why, consider the following two perspectives. First, whenever we have a closure of constraints, we can also see such closure as a corresponding cycle in a pepe-graph. Cycles of processes are called self-organising by \citeAmossio2017makes, but are distinguished from self-determining systems which realise a closure of constraints El-Hani \BBA Nunes-Neto (\APACyear2020). We can therefore understand self-determination as a special case of self-organisation, so studying the latter can inform our understanding of the former (Juarrero, \APACyear2023, p. 101). Second, a cycle in a pepe-graph is collectively impredicative, meaning that it can only exist if all of its constituent processes are already unfolding and generating one another Jaeger \BOthers. (\APACyear2024). It is unclear how such a system could have come to exist in the first place, since none of the processes in this network can begin unless all the other processes have already begun. This is what \citeAhofmeyr2021biochemically refers to as the starting problem, analogous to \citeAturing1936computable’s halting problem. So cycles of processes possess an interesting property which is not shared by noncyclical process-enablement structures.

3.2 Model comparison

Every pepe-graph represents a particular perspective of a system. As a result, for any given system we could draw many distinct pepe-graphs that each represent a different way of looking at the same system. For instance, if two scientists observe a bird in a forest, a biochemist would likely draw a very different pepe-graph to describe it compared to an ecologist. This is a fundamental feature of the nature of scientific modelling. Every model represents a particular point of view and a particular way of thinking. Every model leaves things out and distorts reality in one way or another. It is only by comparing many models together that we can really begin to understand nature Massimi (\APACyear2022).

This doesn’t mean that every pepe-graph is a valid representation of a system. Each time you draw an arrow in a pepe-graph you are claiming the existence of a genuine lawlike dependency that could be investigated. Thus, whether the pepe-graph is accurate or not is ultimately up to empirical investigation. If such an investigation shows that the enablement is wrong then the graph should be updated accordingly.

We now want to develop the mathematical tools required to conduct such a comparative analysis of pepe-graphs. We want to enable the biochemist and ecologist, for example, to compare their points of view with each other and see where and how their models differ, but also to investigate where there are similarities. However, as we conduct these perspective shifts we want to ensure that the appropriate structural features of pepe-graphs remain intact. This desire motivates our following notions of preservation and reflection of closure.

Definition 3.3 (Preservation of closure).

Let GG and HH be graphs and let ϕ:GH\phi\colon G\to H be a weak graph homomorphism. Then ϕ\phi preserves closure if for every cycle CGC\subseteq G there exists a cycle DH[ϕ(C)]D\subseteq H[\phi(C)].

Note that standard graph homomorphisms always preserve closure, but this is not necessarily true for weak graph homomorphisms. Figure 2 illustrates the notion of preservation of closure.

AABB
ϕ\phi
A1A_{1}A2A_{2}A3A_{3}BB
ψ\psi
AABB
Figure 2: Both ϕ\phi and ψ\psi are weak graph homomorphisms, but ψ\psi preserves closure whilst ϕ\phi does not. Both ϕ\phi and ψ\psi map A1AA_{1}\mapsto A, A2AA_{2}\mapsto A, and A3AA_{3}\mapsto A. The induced image of the cycle in the centre graph, under ϕ\phi and ψ\psi, is highlighted in blue in their respective target graphs.

The following result shows that weak graph homomorphisms that preserve closure are composable.

Theorem 3.4.

Let G1,G2G_{1},G_{2}, and G3G_{3} be graphs and ϕ:G1G2\phi\colon G_{1}\to G_{2} and ψ:G2G3\psi\colon G_{2}\to G_{3} be weak graph homomorphisms that preserve closure. Then ψϕ:G1G3\psi\circ\phi\colon G_{1}\to G_{3} also preserves closure.

Proof.

Suppose C1G1C_{1}\subseteq G_{1} is a cycle. Since ϕ\phi preserves closure there exists a cycle C2G2[ϕ(C1)]C_{2}\subseteq G_{2}[\phi(C_{1})], and since ψ\psi preserves closure there exists a cycle C3G3[ψ(C2)]C_{3}\subseteq G_{3}[\psi(C_{2})]. It suffices to show C3G3[(ψϕ)(C1)]C_{3}\subseteq G_{3}[(\psi\circ\phi)(C_{1})]. For this, note that V(ψ(C2))V((ψϕ)(C1))V(\psi(C_{2}))\subseteq V((\psi\circ\phi)(C_{1})), so C3G3[ψ(C2)]G3[(ψϕ)(C1)]C_{3}\subseteq G_{3}[\psi(C_{2})]\subseteq G_{3}[(\psi\circ\phi)(C_{1})]. ∎

The next result establishes useful conditions for determining the preservation of closure based on structural properties of the induced images of cycles under the weak graph homomorphism.

Theorem 3.5 (Preservation test).

Let GG and HH be graphs and let ϕ:GH\phi\colon G\to H be a weak graph homomorphism. Then ϕ\phi preserves closure if and only if for every cycle CGC\subseteq G either

  1. (i)

    |V(H[ϕ(C)])|>1|V(H[\phi(C)])|>1, or

  2. (ii)

    H[ϕ(C)]H[\phi(C)] is a graph with a single vertex and a single loop.

Proof.

We prove the contrapositive in both directions, that is we show that ϕ\phi does not preserve closure if and only if there exists a cycle CGC\subseteq G such that H[ϕ(C)]H[\phi(C)] is a graph with a single vertex and no edges. If ϕ\phi does not preserve closure then there exists a cycle CGC\subseteq G such that the graph H[ϕ(C)]H[\phi(C)] does not contain a cycle, so H[ϕ(C)]H[\phi(C)] must consist of a single vertex and no edges. Conversely, if there exists a cycle CGC\subseteq G such that the graph H[ϕ(C)]H[\phi(C)] has a single vertex and no edges then H[ϕ(C)]H[\phi(C)] contains no cycles, hence ϕ\phi does not preserve closure. ∎

When considering pepe-graphs we want homomorphisms to preserve additional structure with respect to the processes.

Definition 3.6 (Homomorphism of pepe-graphs).

Consider a system SS, and let GG and HH be pepe-graphs which model subsystems of SS over the same time interval. Let ϕ:GH\phi\colon G\to H be a weak graph homomorphism that preserves closure. Then ϕ\phi is a homomorphism of pepe-graphs if for every process pXV(G)p_{{}_{X}}\in V(G) we have that pXp_{{}_{X}} is happening within ϕ(pX)\phi(p_{{}_{X}}). That is, XX is contained within YY where ϕ(pX)=qYV(H)\phi(p_{{}_{X}})=q_{{}_{Y}}\in V(H).

We also refer to homomorphisms of pepe-graphs as pepe-graph homomorphisms, or simply homomorphisms if it is sufficiently clear that we are referring to pepe-graphs. We will often abbreviate the notation when referring to a pepe-graph homomorphism and say ‘pp is happening within ϕ(p)\phi(p)’ when we really mean ‘pXp_{{}_{X}} is happening within qYq_{{}_{Y}}’ where qY=ϕ(pX)q_{{}_{Y}}=\phi(p_{{}_{X}}).

We now show that pepe-graph homomorphisms are composable.

Theorem 3.7.

Consider a system SS. Let G1,G2G_{1},G_{2}, and G3G_{3} be pepe-graphs representing subsystems of SS over the same time interval, and let ϕ:G1G2\phi\colon G_{1}\to G_{2} and ψ:G2G3\psi\colon G_{2}\to G_{3} be pepe-graph homomorphisms. Then ψϕ:G1G3\psi\circ\phi\colon G_{1}\to G_{3} is also a pepe-graph homomorphism.

Proof.

By Theorem 3.4, ψϕ\psi\circ\phi preserves closure. If pXp_{{}_{X}} is a process in SS with pXV(G1)p_{{}_{X}}\in V(G_{1}) then, since ϕ\phi is a pepe-graph homomorphism, there exists a process qYq_{{}_{Y}} in SS such that ϕ(pX)=qYV(G2)\phi(p_{{}_{X}})=q_{{}_{Y}}\in V(G_{2}) and pXp_{{}_{X}} happens within qYq_{{}_{Y}}. Further, since ψ\psi is a pepe-graph homomorphism, there exists a process rZr_{{}_{Z}} in SS such that ψ(ϕ(pX))=ψ(qY)=rZV(G3)\psi(\phi(p_{{}_{X}}))=\psi(q_{{}_{Y}})=r_{{}_{Z}}\in V(G_{3}) and ϕ(pX)=qY\phi(p_{{}_{X}})=q_{{}_{Y}} happens within rZr_{{}_{Z}}. It follows that pXp_{{}_{X}} happens within ψ(ϕ(pX))\psi(\phi(p_{{}_{X}})), hence ψϕ\psi\circ\phi is a pepe-graph homomorphism. ∎

Homomorphisms of pepe-graphs allow us to translate between two perspectives representing the same system, mapping all of the relevant structural properties of one pepe-graph to another. However, it is possible for the target graph to contain more structure than the source graph. If every cycle in the target graph is reflected appropriately in the source graph then we say that a homomorphism reflects closure.

Definition 3.8 (Reflection of closure).

Let GG and HH be graphs and let ϕ:GH\phi\colon G\to H be a weak graph homomorphism. Then ϕ\phi reflects closure if for every cycle DHD\subseteq H there exists a cycle CGC\subseteq G such that DH[ϕ(C)]D\subseteq H[\phi(C)].

We now show that weak graph homomorphisms that reflect closure are composable.

Theorem 3.9.

Let G1,G2G_{1},G_{2}, and G3G_{3} be graphs, and let ϕ:G1G2\phi\colon G_{1}\to G_{2} and ψ:G2G3\psi\colon G_{2}\to G_{3} be weak graph homomorphisms that reflect closure. Then ψϕ:G1G3\psi\circ\phi\colon G_{1}\to G_{3} reflects closure.

Proof.

If C3G3C_{3}\subseteq G_{3} is a cycle then, since ψ\psi reflects closure, there exists a cycle C2G2C_{2}\subseteq G_{2} such that C3G3[ψ(C2)]C_{3}\subseteq G_{3}[\psi(C_{2})]. Further, since ϕ\phi reflects closure there exists a cycle C1G1C_{1}\subseteq G_{1} such that C2G2[ϕ(C1)]C_{2}\subseteq G_{2}[\phi(C_{1})]. It suffices to show C3G3[(ψϕ)(C1)]C_{3}\subseteq G_{3}[(\psi\circ\phi)(C_{1})]. For this, note that V(ψ(C2))V((ψϕ)(C1))V(\psi(C_{2}))\subseteq V((\psi\circ\phi)(C_{1})), so C3G3[ψ(C2)]G3[(ψϕ)(C1)]C_{3}\subseteq G_{3}[\psi(C_{2})]\subseteq G_{3}[(\psi\circ\phi)(C_{1})]. Hence ψϕ\psi\circ\phi reflects closure. ∎

The following result establishes a useful condition for determining the reflection of closure.

Theorem 3.10 (Reflection test).

Let GG and HH be graphs and let ϕ:GH\phi\colon G\to H be a weak graph homomorphism. If there exists a cycle CGC\subseteq G such that H[ϕ(C)]=HH[\phi(C)]=H then ϕ\phi reflects closure.

Proof.

If CGC\subseteq G is a cycle such that H[ϕ(C)]=HH[\phi(C)]=H then for every cycle DD in HH we have DH=H[ϕ(C)]D\subseteq H=H[\phi(C)]. ∎

A homomorphism which reflects closure provides a correspondence between all of the key structural features of the two pepe-graphs, namely the processes and the cycles. As such, these maps indicate a greater similarity between two pepe-graphs. We therefore call these maps homorheisms to distinguish them from homomorphisms which do not necessarily reflect closure.

Definition 3.11 (Homorheism).

Consider a system SS, and let GG and HH be pepe-graphs which model subsystems of SS over the same time interval. Let ϕ:GH\phi\colon G\to H be a homomorphism of pepe-graphs. If ϕ\phi also reflects closure then ϕ\phi is a homorheism. We will denote homorheisms using the modified arrow \mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}.

The notion of a homorheism is illustrated in Figure 3. The following result shows that homorheisms are composable.

A1A_{1}H1H_{1}
ϕ\phi
A1A_{1}A2A_{2}A3A_{3}H2H_{2}
ψ\psi
AAH3H_{3}
Figure 3: H1H_{1}, H2H_{2}, and H3H_{3} are pepe-graphs. Since pepe-graph homomorphisms preserve closure, they exist only from H1H2H_{1}\to H_{2}, H1H3H_{1}\to H_{3}, and H2H3H_{2}\to H_{3}, but not from H2H1H_{2}\to H_{1}, H3H2H_{3}\to H_{2}, and H3H1H_{3}\to H_{1}. Of the three pepe-graph homomorphisms only ψ:H2H3\psi\colon H_{2}\to H_{3}, where ψ\psi maps each vertex in H2H_{2} to AA in H3H_{3}, reflects closure and it is therefore a homorheism.
Theorem 3.12.

Consider a system SS. Let G1,G2G_{1},G_{2}, and G3G_{3} be pepe-graphs representing subsystems of SS over the same time interval. Let ϕ:G1G2\phi\colon G_{1}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}G_{2} and ψ:G2G3\psi\colon G_{2}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}G_{3} be homorheisms. Then ψϕ:G1G3\psi\circ\phi\colon G_{1}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}G_{3} is also a homorheism.

Proof.

Follows from Theorems 3.7 and 3.9. ∎

Homorheisms allow us to identify two vantage points that convey roughly the same description of the system under investigation, though the source graph may be a more finely-grained perspective than the target graph. Note, however, that homorheisms are not symmetric so they are not an equivalence relation on pepe-graphs. That is, just because there exists a homorheism GHG\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}H does not necesarily mean there is a homorheism HGH\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}G.

But if there are no homorheisms between two pepe-graphs that represent the same system, then there is at least one fundamentally different organisational feature between them. Namely one must contain a cycle that the other does not capture in any way. Graphs like these are therefore different vantage points that offer significantly different views of the same system, with respect to their organisational structure. It is possible, however, to use intermediate pepe-graphs to translate between perspectives when there are no homorheisms directly between the graphs available. We will demonstrate this process in Section 4.3.

As a linguistic note to conclude the discussion here, the term homorheism translates to same flow. It is inspired by C. H. Waddington’s neologism homeorhesis to describe dynamical systems that return to the same dynamical trajectory, following a disturbance Waddington (\APACyear1957).

3.3 Self-enablement

Our formalism allows for a pepe-graph to have loops, but we have not yet described what this means conceptually. We wish to avoid paradoxes of self-enablement, so we do not want a loop in a pepe-graph to indicate that a process is directly enabling itself. Instead, we use loops as a convenient shorthand. Consider a system SS of processes and direct enablements (the prokaryotic cell is a good example). If this system SS contains at least one cycle then instead of writing all of the processes and enablements within SS we can simply write SSS\to S. Therefore, a loop within a pepe-graph should always indicate to the reader that the loop in question is really a finer-grained network of processes containing at least one cycle. There are, therefore, really two kinds of arrows in pepe-graphs: direct enablements (not loops) and loops.

In contrast, many models often choose to ignore lower-level organisation and treat systems as if they were not self-organising. For instance, organisms within many evolutionary models are treated as objects that can be mostly reduced down to their genes Lewontin (\APACyear1983); Walsh (\APACyear2018); Chiu (\APACyear2022). For this reason, population genetics has been referred to pejoratively as ‘beanbag genetics’ Rao \BBA Nanjundiah (\APACyear2011). We can see this disappearance of the self-organising organism in, for example, the Price equation where the self-organising properties of individuals are ignored in favour of a simple statistical equation Price (\APACyear1970); Frank (\APACyear2012). Accordingly, if we were to treat such models as reality and then constructed the resulting pepe-graphs, loops like ‘Organism \to Organism’ would never appear.

With the use of loops as a shorthand we need to keep in mind that more detailed descriptions of a system may result in the loss of a homorheism. For instance, suppose GG is the graph with a single process ‘Hepatocyte’ and a single loop (since cells are self-organising), and HH is a similar graph with a single process ‘Liver’ and a single loop. Then there is a natural homorheism ϕ:GH\phi\colon G\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}H, but clearly GG and HH model quite different systems. If we unpack the loops in GG and HH and express these graphs as more detailed networks of the systems in question then we lose the homorheism: once we unpack the process ‘Liver’ we will see how intricate and complex the network here really is, with hundreds of processes and cycles between many different cell types. Indeed, there are far more process-enablement cycles within the liver compared to those just occurring within a single hepatocyte.

4 Applications

We have now built up enough of our formalism to begin comparing a wide array of perspectives from across the biosphere. In practice, our analysis in the following subsections proceeds by first identifying the processes and enablements within a system, from a particular perspective, and then drawing corresponding pepe-graphs. Since multiple perspectives are always possible there are bound to be local differences between these graphs. For instance, we may view two processes pp and qq where pp enables qq as distinct, so pqp\to q, or we may see pp and qq as a single compound process. Homomorphisms and homorheisms will allow us to compare these differences precisely and concretely.

4.1 Autocatalytic sets

The concept of autocatalytic sets emerged during the 1980s Dyson (\APACyear1982); Kauffman (\APACyear1986). Dyson and Kauffman were both interested in the origins of life, and wanted to explore how random sets of chemical reactions could come together to form an organised network. Their ideas were then refined via a mathematical formalism known as RAF sets: reflexively autocatalytic systems generated by a food set Hordijk \BBA Steel (\APACyear2004); Hordijk \BOthers. (\APACyear2012). Here, reflexive autocatalysis means that every reaction in the system is catalysed by at least one molecule already in the system, and food-generated means that all the molecules in the system can be ultimately generated from the food set through chains of reactions within the system (Figure 4).

F1F_{1}F2F_{2}F3F_{3}F4F_{4}B1B_{1}B2B_{2}B3B_{3}Reaction 1Reaction 2Reaction 3
Figure 4: A RAF set. The food set is 𝒮={F1,F2,F3,F4}\mathcal{S}=\{F_{1},F_{2},F_{3},F_{4}\}, while B1B_{1}, B2B_{2}, and B3B_{3} are intermediate molecules produced by the network. Dotted arrows indicate catalysis and solid arrows represent molecules entering/exiting a reaction. Modified from Figure 4 in Hordijk, Steel, and Kauffman (2012).

To study autocatalytic sets using the language of pepe-graphs, we first need to consider the kinds of processes and enablement relationships going on within these chemical systems. The processes in each case are relatively clear: they are the chemical reactions of the autocatalytic set, together with the process of food diffusing freely into the system. Since at least one food molecule is required for each reaction displayed in Figure 4, food diffusion must directly enable each reaction. However, whether reactions directly enable each other depends on the nature of molecule degradation and the rate constants of the reactions. For instance, if molecule B3B_{3} degrades rapidly within the system, then Reaction 2 will directly enable Reaction 3 because Reaction 2 will need to continuously produce B3B_{3} in order to keep Reaction 3 going (Figure 4). On the other hand, if B3B_{3} is more stable and remains abundant within the system then this direct enablement could be undermined. For example, if B3B_{3} has been stockpiled by the system in large quantities, then the occurrence of Reaction 2 may be irrelevant to Reaction 3. For now, we will assume that this RAF set operates in a high degradation environment, so that whenever a reaction ceases then its respective products disappear almost instantaneously from the system111We could also allow for a low-degradation environment, but add the diffusion processes of each molecule into the pepe-graph. For simplicity we instead assume a high-degradation environment, yielding direct enablements between reactions and ignoring the diffusion processes.. We can then draw two pepe-graphs to describe the RAF sets (Figure 5).

Food DiffusionReaction1Reaction2Reaction3R1R_{1}
ϕ1\phi_{{}_{1}}
Food DiffusionReactions1+3Reaction2R2R_{2}
Figure 5: Two pepe-graph representations of the RAF set from Figure 4.

Interestingly, the pepe-graph R1R_{1} does not have a cycle, whilst R2R_{2} does (Figure 5). This demonstrates that organisational closure is partially perspective dependent. Unlike \citeAcusimano2020objectivity, we do not see this failure of objectivity as a disadvantage of the pepe-graph framework. Rather, we see it is a significant advantage because it introduces a way to identify qualitatively different analyses of the system under investigation. A comparison of these perspectives can then inform future empirical investigations and help us to better understand the system. For instance, if we see Reactions 1 and 3 as really two parts of the broader process ‘Reactions 1+31+3’, then the analysis provided in R2R_{2} will make more sense. The system contains a cycle and this part of the system therefore achieves organisational closure. Thus, depending on where this system is located (e.g. whether it be in a Petri dish or a cell) we may even expect the system to evolve structures, or act in a particular way, to preserve this cycle (Moreno \BBA Mossio, \APACyear2015, p. 71).

In a similar manner, using the homomorphism ϕ1\phi_{1}, the scientist can come to understand how a collectively impredicative system like R2R_{2} may have evolved. This is because, in R1R_{1}, Reaction 3 does not enable Reactions 1 or 2. Thus, we could have had a reaction network with Reactions 1 and 2, but not Reaction 3. We therefore have a solution to \citeAhofmeyr2021biochemically’s starting problem for this case. From one point of view, we may see a cycle that cannot seem to start without all the processes already existing (R2R_{2}). But from another point of view, related via an appropriate homomorphism (ϕ1\phi_{1}), we can see how the system may have actually come to exist, without any collective impredicativity (R1R_{1}). It is not clear to us if all collectively impredicative systems can be dissolved like this using a perspective shift, however our example here shows that it is certainly possible in particular situations.

4.2 Autopoiesis

Humberto Maturana and Francisco Varela together developed the concept of autopoiesis, which proved to be a hugely influential account of life Maturana \BBA Varela (\APACyear1980). The concept went through a number of iterations, but we will focus on a later definition given by \citeA[p. 34]varela2000fenomeno and discussed in relation to cognition in \citeA[pp. 99–103]thompson2007mind. Here, a system is defined as being autopoietic if the following criteria all hold:

  1. 1.

    The system has a semipermeable boundary constructed from molecular components. This boundary allows us to discriminate between the inside and outside of the system.

  2. 2.

    The components of the boundary are produced by a metabolic reaction network (MRN) located within the boundary.

  3. 3.

    There is an interdependence between Criteria 1 and 2. That is, the MRN is regenerated and maintained because of the existence of the boundary which keeps all of the molecules in the system in close proximity to each other.

For an autopoietic system to be maintained, much like an autocatalytic set, there must also be a ready supply of external nutrients that can flow into the system. Although this final criterion is not made explicit by Maturana and Varela, it is implicitly assumed throughout their work Varela (\APACyear2000). The prototypical example of an autopoietic system is the prokaryotic cell. It is still debated whether multicellular organisms, eusocial insect colonies, or even human societies are also autopoietic. The discussion around these trickier cases typically centres on how one defines a semipermeable boundary Thompson (\APACyear2007).

From this basic outline we can see that a minimal autopoietic system has three key processes: the MRN, boundary maintenance, and the passive diffusion of food into the system. Food diffusion directly enables the MRN since it provides the basic reactants for the system, and boundary maintenance and the MRN directly enable each other because of Criteria 2 and 3 (Figure 6).

FoodDiffusionMRNBoundaryMaintenanceA1A_{1}
FoodDiffusionMRNBoundaryMaintenanceA2A_{2}
TransportMRNBoundaryMaintenanceA3A_{3}
Figure 6: Three different kinds of autopoietic systems represented as pepe-graphs. MRN = metabolic reaction network. A1A_{1}: minimal autopoiesis. A2A_{2}: autopoiesis with an organisationally closed metabolic reaction network. A3A_{3}: autopoiesis with facilitated diffusion, active transport, and endocytosis of relevant molecules.

There is no direct requirement that the MRN itself be organisationally closed, but in practice life has incorporated many cycles into its metabolic reaction networks (for example, the Krebs cycle). For these more complex autopoietic systems, the MRN itself achieves organisational closure (Figure 6). Some autopoietic systems also play an active role in bringing nutrients and other molecules into the system. For instance, cells often use processes like facilitated diffusion, primary and secondary active transport, and endocytosis to control which molecules enter into their internal metabolic reaction networks Jeckelmann \BBA Erni (\APACyear2020). In these autopoietic systems the MRN directly enables transport by producing the relevant transport proteins (for example GLUTs, Na+/K+-ATPase, and Clathrin), and the maintenance of the boundary directly enables these transport processes since they would not occur without a cell membrane (Figure 6).

So using the language of pepe-graphs we can see at least three distinct kinds of autopoietic systems, with differing levels of organisation (Figure 6). We can also connect this perspective of life to autocatalytic sets by identifying homomorphisms between the two accounts. To do this, we first need to consider a system that could be viewed either through the lens of autopoiesis or autocatalytic sets. For example, we can imagine an extended version of the RAF set in Figure 4 where molecule B2B_{2} is a lipid which is produced in sufficient quantity to form a semipermeable boundary. This boundary then keeps all of the relevant molecules of the RAF set together, and is completely permeable to the food molecules F1F_{1}, F2F_{2}, F3F_{3}, and F4F_{4}. Depending on how we view the RAF set as discussed in the previous section, and depending on how much of the system we would like to study, this analysis yields four related pepe-graphs (Figure 7).

Food DiffusionReaction1Reaction2Reaction3R1R_{1}
ϕ1\phi_{{}_{1}}
Food DiffusionReactions1+3Reaction2R2R_{2}
ϕ2\phi_{{}_{2}}
ϕ3\phi_{{}_{3}}
FoodDiffusionMRNBoundaryMaintenanceA1A_{1}
ϕ4\phi_{{}_{4}}
FoodDiffusionMRNBoundaryMaintenanceA2A_{2}
Figure 7: Process-enablement graph homomorphisms between autocatalytic set representations and autopoietic representations. Each arrow between graphs is a homomorphism and ϕ2\phi_{2} and ϕ3\phi_{3} are homorheisms onto their induced images. The map ϕ2\phi_{2} (respectively ϕ3\phi_{3}) sends each of the reactions in R1R_{1} (respectively R2R_{2}) to the MRN in A1A_{1} (respectively A2A_{2}). MRN = metabolic reaction network.

Figure 7 illustrates quite well how loops in pepe-graphs should be understood. At first glance, ϕ4\phi_{4} can be analysed in the same way as ϕ1\phi_{1}, as discussed in the previous section. But in fact, from this zoomed out perspective, we really have no understanding as to how closure is being realised in A2A_{2} compared to A1A_{1}: it is all hidden away inside the loop. Once we contextualise what we mean by the ‘metabolic reaction network’, via either ϕ2\phi_{2} or ϕ3\phi_{3}, it then becomes clear what is actually happening at the more detailed level.

4.3 (F,A)-systems

If there is a homomorphism between two pepe-graphs, then the relationship between them is relatively easy to analyse. We can quickly check which processes and enablements match up, what gets abstracted away from the source graph, and what extra phenomena the target graph may be able to capture that are invisible to the source graph. However, homomorphisms often do not exist between pepe-graphs, so it can be useful to have some more general tools for comparison. In particular, when there are no homomorphisms between two pepe-graphs we can sometimes develop an intermediary graph which indirectly shows how the two models are related.

To illustrate this process we will compare autopoietic systems, as described in Section 4.2, to a biochemically realisable version of Robert Rosen’s metabolism-repair (M,R)(M,R)-systems Rosen (\APACyear1991). Rosen’s theory of biological organisation and (M,R)(M,R)-systems is highly abstract and has essentially no reference to biochemistry. As a result, Rosen’s framework has been largely ignored by biologists. However, his work has recently been refined by \citeAhofmeyr2019basic, hofmeyr2021biochemically to fit it onto real biochemical processes. The result is what Hofmeyr calls fabrication-assembly (F,A)(F,A)-systems. As considered within the context of the cell, an (F,A)(F,A)-system has three classes of processes: (1) covalent metabolic chemistry (fabrication); (2) supramolecular chemical processes (assembly) and; (3) maintenance of the intracellular milieu via the transporting of ions.

Covalent molecular chemistry or fabrication involves the transport of nutrients into the cell, the degradation of nutrients into the building blocks of metabolism (for example amino acids, lipids, and nucleotides) and the synthesis of macromolecules like polypeptides and nucleic acids from these building blocks. Supramolecular chemistry or assembly involves the folding of proteins, the self-assembly of ribosomes, and the maintenance of the cell membrane via noncovalent interactions. This kind of chemistry is facilitated by the regulated ionic composition, temperature, and pH of the intracellular milieu, and by chaperone proteins. Finally, the maintenance of the intracellular milieu is primarily driven by the transport of ions into and out of the cell via membrane-bound transporters.

The enablements between these three classes of processes within an (F,A)(F,A)-system can be summarised via three cycles (Figure 8).

Nutrient TransportMembraneMaintenanceCovalentChemistryProteinFoldingIonTransportI𝒫I_{\mathcal{P}}
ϕ5\phi_{{}_{5}}
ϕ6\phi_{{}_{6}}
FabricationAssemblyMaintenance ofthe intracellular milieuFAF_{A}
TransportMRNBoundaryMaintenanceA3A_{3}
Figure 8: A comparison of the (F,A)(F,A)-system model of the cell to an autopoietic perspective using pepe-graphs. FAF_{A}: (F,A)(F,A)-system representation. I𝒫I_{\mathcal{P}}: intermediate pepe-graph. A3A_{3}: autopoietic representation. There exist homorheisms ϕ5:I𝒫FA\phi_{5}\colon I_{\mathcal{P}}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}F_{A} and ϕ6:I𝒫A3\phi_{6}\colon I_{\mathcal{P}}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}A_{3}, which are described using a colouring to match vertices in I𝒫I_{\mathcal{P}} with those in FAF_{A} or A3A_{3} (see Figure A.1). There are no homomorphisms of pepe-graphs between FAF_{A} and A3A_{3}, as explained in the main text. MRN = metabolic reaction network.

Fabrication enables assembly by providing the raw materials that need to be assembled. This enablement is direct because as macromolecules are being produced within the cell, they immediately begin to self-assemble as they interact with the intracellular milieu. Conversely, assembly enables fabrication by providing the enzymes needed to catalyse the covalent chemistry of the cell. This enablement is direct because enzymes do not have a single static conformation, they are highly dynamic entities that have an ensemble of equilibrium configurations Yang \BOthers. (\APACyear2003); Nicholson (\APACyear2019). So we have the first cycle: Fabrication \rightleftarrows Assembly.

In turn, the stochastic fluctuations of proteins are modulated by the constantly shifting, local composition of the intracellular milieu. Thus, the maintenance of the intracellular milieu directly enables assembly. Conversely, the maintenance of the cell membrane directly enables ions to be transported across it via membrane transporters, since without a membrane there would be no ion channels. Thus, assembly directly enables the maintenance of the intracellular milieu. This enablement is direct because the cell membrane is constantly shifting and undergoing repair. This gives us a second cycle: Assembly \rightleftarrows Maintenance of the intracellular milieu.

Finally, assembly has a loop because protein folding and cell membrane maintenance, both processes occuring within assembly, directly enable each other. Maintenance of the cell membrane keeps ions within the cell, directly enabling protein folding. Conversely, the folding of enzymes like flippase and scramblase is needed to keep the membrane regulated and intact.

Both the (F,A)(F,A)-system and autopoiesis are models of the biological cell, and we can model each with the pepe-graphs FAF_{A} and A3A_{3}, respectively, as shown in Figure 8. But it is impossible to directly compare FAF_{A} and A3A_{3} via a pepe-graph homomorphism. This is because the two perspectives partition the processes of the cell differently. Explicitly, assembly involves both the maintenance of the cell’s boundary and protein folding; the latter is a part of the MRN as viewed from the perspective of autopoiesis. So if we wanted a pepe-graph homomorphism ϕ:FAA3\phi\colon F_{A}\to A_{3}, we would need to map assembly to both boundary maintenance and the MRN. Conversely, the MRN also involves the synthesis of nucleic acids from building blocks, which is a part of fabrication in FAF_{A}. So if there were a pepe-graph homomorphism ψ:A3FA\psi\colon A_{3}\to F_{A}, we would need to map the MRN onto both fabrication and assembly. Since a pepe-graph homomorphism cannot map a single vertex to two vertices, there are no pepe-graph homomorphisms between FAF_{A} and A3A_{3}. Yet, both FAF_{A} and A3A_{3} model the same system, so we would expect some kind of connection between them. Indeed, we can find one, if we construct an intermediate perspective as follows.

First, we need to choose a finer partition 𝒫\mathcal{P} of the processes described in each of the models. There are many choices we could make here, and different sets of finer-grained processes will reveal different similarities between the models. As an illustrative example we will choose the set 𝒫:=\mathcal{P}:= {nutrient transport, membrane maintenance, protein folding, ion transport, covalent chemistry}. Here, covalent chemistry is equivalent to (F,A)(F,A)-system fabrication, excluding nutrient transport, and is also equivalent to the autopoietic MRN, excluding protein folding. Membrane maintenance refers to the trafficking of lipids and proteins to the cell membrane, the lateral diffusion of membrane components, the flipping of lipids from one side of the membrane to the other, membrane repair processes, and the containment of cytosolic molecules within the cell as they bounce off the inner plasma membrane. Ion transport refers to the transport of ions into and out of the cell via facilitated diffusion and active transport.

Second, we match each of the processes in 𝒫\mathcal{P} to processes in the target graphs in question. We have illustrated this matching for our example in Figure 8 using two colourings. We must ensure that the choice of 𝒫\mathcal{P} is fine enough such that every process p𝒫p\in\mathcal{P} is only happening within a single process in each of the target graphs. This will ultimately allow us to construct well-defined homomorphisms.

Third, we construct the corresponding pepe-graph, I𝒫I_{\mathcal{P}}, by drawing in the direct enablements between processes in 𝒫\mathcal{P}. If the arrows are drawn correctly then, using the correspondence from the previous step, we should have weak graph homomorphisms from I𝒫I_{\mathcal{P}} to the two target graphs. Finally, we need to check whether these weak graph homomorphisms preserve or reflect closure. If this is successful then we will have indeed constructed pepe-graph homomorphisms, or possibly homorheisms, from the intermediate pepe-graph to the two target graphs. We prove this for our case in Theorems A.1 and A.2 in the Appendix. If, however, the weak graph homomorphisms do not preserve closure then we can choose a different set 𝒫\mathcal{P} and repeat the steps above. Alternatively, if we cannot construct homomorphisms from I𝒫I_{\mathcal{P}} to the target graphs then it may be an indication that there is a misplaced or missing arrow in these target graphs. Accordingly, this process of constructing an intermediate pepe-graph can facilitate model development, as well as improve existing models.

The existence of the graph I𝒫I_{\mathcal{P}} and corresponding homorheisms ϕ5:I𝒫FA\phi_{5}\colon I_{\mathcal{P}}\to F_{A} and ϕ6:I𝒫A3\phi_{6}\colon I_{\mathcal{P}}\to A_{3} in Figure 8 demonstrates that lower-level perspectives of the cell, with more molecular detail, can be organisationally compatible with more abstract models like autopoiesis. We have also shown that autopoietic and (F,A)(F,A)-system models can both arise from the same set of biological processes.

Additionally we have shown that there are no pepe-graph homomorphisms between FAF_{A} and A3A_{3}, so they must have at least one different organisational feature: and indeed they do, as the graphs FAF_{A} and A3A_{3} have different cycle structures. Thus, the processes prioritised as being central to the self-organisation of the cell are different in each case, since each perspective views the self-organisation of the cell itself in different ways. From the (F,A)(F,A)-system perspective, the key process in the model is assembly as it sits within all three cycles in FAF_{A}. The way in which the cell arranges its parts from disorganised building blocks into functional assemblies is the central process that keeps the cell alive. But from the autopoietic perspective given in A3A_{3}, the key process is the metabolic reaction network as it also sits within all three cycles.

If we analyse the constraints on the central processes in each model, we can also flesh out more of the differences between these two perspectives. The key constraint on assembly is the maintenance of the intracellular milieu, but the key constraint on the MRN is the maintenance of a boundary. These are similar processes, though they are ultimately different. The maintenance of the intracellular milieu involves keeping the ionic composition of the cytoplasm stable, whilst the maintenance of a boundary refers to keeping all of the molecules in the cell within spatial proximity, thereby enabling reactants to interact and form a reaction network. This highlights a subtle but important distinction between these two models.

Another important remark is that because ϕ5\phi_{5} and ϕ6\phi_{6} are homorheisms, not just homomorphisms, all of the cycles in FAF_{A} (respectively A3A_{3}) can be found in I𝒫I_{\mathcal{P}} and are therefore embodied in A3A_{3} (respectively FAF_{A}) via ϕ6\phi_{6} (respectively ϕ5\phi_{5}). This means that we have not ‘lost’ any of the self-organisational features when we move from FAF_{A} to A3A_{3} or vice versa, as they are just rearticulated in a different way and into different cycles. In short, we could simply say that FAF_{A} and A3A_{3} model the same self-organising processes in the cell, but they each have a different focus as to what they consider important. In the Appendix we consider the pepe-graph I𝒫I_{\mathcal{P}} in more detail, showing that each of the edges in I𝒫I_{\mathcal{P}} is a direct enablement. We also prove that ϕ5\phi_{5} and ϕ6\phi_{6} are homorheisms.

4.4 Problem cases

There are some natural problem cases that arise whenever philosophers of biology ask themselves the question ‘what is life?’. In this section we have no intention of resolving demarcation problems surrounding what counts as a living organism or not. We will, however, translate some of these interesting edge cases into the language of pepe-graphs, allowing for them to be studied from a new angle. We will also focus less on homomorphisms between pepe-graphs in this section and instead look at when cycles are present or absent from the graphs. This will facilitate our study of organisational closure in these lifelike systems.

4.4.1 Viruses

Viruses pose many puzzles for biologists and philosophers. They have genes but no internal metabolism, they can replicate but only by infecting a cell, they can evolve but have no ability to control their movement. As a virus floats around outside the confines of a cell it is difficult to characterise it as ‘alive’. All that the virus can do at this stage of its life cycle is float around until it hopefully bumps into a cell that it can infect. This perspective is shown in the graph V1V_{1} in Figure 9.

VirionDiffusionV1V_{1}
NativeMRNViral Gene ExpressionUpregulatedGlycolysisV2V_{2}
InfectionViral GeneExpressionReleaseVirionDiffusionV3V_{3}
Figure 9: Process-enablement graph representations of viruses from three perspectives. V1V_{1}: Virion diffusion alone. V2V_{2}: A ribovirocell, which is an infected host cell that has been partially remodelled to produce virions. V3V_{3}: A population-level perspective of viruses. Native MRN refers to the basal metabolic reaction network of the cell before it has been infected by a virus. Release refers to the release of virions from an infected cell via lysis, budding or exocytosis.

Here, it is useful to make a distinction between the virus as a process and the virus as a particle, also known as a virion.

As a process the virus has several stages to its life cycle, some more active than others. For instance, if a virus infects a cell then multiple processes enter the picture beyond just the diffusion of virions. At this point, the virus begins its active phase and co-opts the cell’s native cellular processes to transcribe and translate its own viral genes (Figure 9). The infected host cell is thus turned into a manufacturer of virion particles. During this transformation, viruses will also alter the host cell’s own metabolism for their own benefit. For instance, most viruses induce aerobic glycolysis as a necessary process to enable viral replication Sanchez \BBA Lagunoff (\APACyear2015). \citeAforterre2011manipulation argues that this system, which they call a virocell, constitutes the living form of the virus. Explicitly, a virocell is one that completely remodels a cell to produce virion particles, and inhibits native cellular processes like cell division. On the other hand, an infected cell that can both divide by cell division and produce virion particles is called a ribovirocell Forterre (\APACyear2016). If we conceptualise such a ribovirocell with a pepe-graph like V2V_{2}, we can see some of the organisational features present within this system (Figure 9). A ribovirocell inherits the organisational closure from the uninfected cell since it also has a metabolic reaction network (see Section 4.2). But additionally, by inducing new metabolic pathways, such as aerobic glycolysis, the virocell also forms new cycles of processes and enablements. The production of viral proteins is necessary for the induction of these new metabolic pathways and the induction of these pathways is necessary to continue the expression of viral genes. In this way, the virocell or ribovirocell embodies a unique form of organisational closure, distinct from either the isolated virion particle or the uninfected cell.

If we zoom out even further and look at a whole population of viruses and their corresponding life cycle then we can generate yet another perspective to consider. Across a population of viruses we can reasonably assume that each stage of the viral life cycle will occur in one place or another. For instance, if we consider the global population of influenza A viruses, then a patient with the flu in Australia would have infected cells actively transcribing and translating influenza A genes as well as other cells in the middle of releasing virion particles, a child in China may have just inhaled a sneeze from an infected patient beginning an infection, and there may be influenza A virion particles floating around the air of a hospital in the US. We express these contemporaneous processes, extended over a global spatial region, in the pepe-graph V3V_{3} (Figure 9). Importantly, if we restrict our attention to a single infected cell then the processes in V3V_{3} will not be contemporaneous, in fact only a subset of them will be occurring so V3V_{3} would not be a valid pepe-graph.

It is perhaps a little strange to consider a whole viral population extended out over various stages of its life cycle as a single organisationally closed system as in V3V_{3} (Figure 9). Nonetheless, this is the perspective traditionally taken by developmental systems theory, or DST Oyama (\APACyear2000). From the perspective of DST, we can study the evolution of life not through the shuffling of genes, or even the action of the environment on the organism (or vice versa), but instead as whole developmental systems and life cycles changing over time and being jointly determined by multiple causes Oyama \BOthers. (\APACyear2001). A common critique of DST is that there are no clear boundaries to a biological system from this perspective Griesemer (\APACyear2006). For instance, to study the evolution of influenza A through the lens of DST it would be just as important to study the mutation of genes within the viral genome as well as to study the introduction of new vaccines targeting the virus: both contribute to the evolution of the influenza A developmental system. But once again using the philosophy of perspectival realism, we can contextualise DST perspectives as just one way of studying life among many other possibilities. Indeed, we could recover many of the boundaries lost from the DST perspective by constructing several different pepe-graphs of the system under investigation, such that each one would have a particular boundary. We could then compare and contrast each of these perspectives using pepe-graph homomorphisms. Such an in-depth investigation is beyond the scope of this paper, but at a minimum we wish to highlight the possibility of untangling the blurriness of DST using pepe-graphs.

4.4.2 The water cycle and the candle

We include the following two examples in our discussion to highlight that pepe-graphs can really be used to model any physical system whether living or nonliving.

First, we consider the global222Much like the viral life cycle in Section 4.4, over a wide enough spatial region all of the processes in W1W_{1} will be contemporaneous, and over too narrow a spatial region only a subset will be occurring. To avoid this issue we simply consider the global water cycle. water cycle in W1W_{1} (Figure 10).

Evaporation and TranspirationCloud FormationPrecipitationRunoff and InfiltrationW1W_{1}
VaporisationCombustionMelting and Capillary ActionW2W_{2}
Figure 10: Process-enablement graph representations of the global water cycle (W1W_{1}) and a candle (W2W_{2}). We use a broad definition of precipitation to include all processes that involve liquid water contacting the ground, such as glacial melting and agricultural irrigation and denote this as ‘Precipitation’. We also use the term vaporisation to refer to the local vaporisation of liquid wax in the wick of the candle.

The direct enablement structure of W1W_{1} is as follows. Evaporation and transpiration produce water vapour which leads to cloud formation. Once the water droplets within clouds grow large enough in size, precipitation begins to occur. As rain hits the ground it immediately begins to infiltrate the ground and runoff occurs. Admittedly, other processes are also responsible for infiltration and runoff such as glacial melting and agricultural irrigation. So to ensure that we have a direct enablement we expand our definition of precipitation to include all processes that involve water contacting the ground, and refer to this broader definition as ‘Precipitation’ (Figure 10). This way, Precipitation is indeed a necessary process for infiltration and runoff to occur. Finally, the infiltration of water into the soil is necessary for it to be absorbed by plant roots, commencing transpiration. The global water cycle thus achieves a cycle of processes and enablements.

Second, we consider a simplified model of the candle, given by the graph W2W_{2} in Figure 10, which achieves organisational closure. Here, the heat produced by combustion melts the solid wax near the base of the candle’s flame and also vaporises the liquid wax in the wick. The wax vapours react with oxygen in the air around the candle, maintaining the combustion reaction. Liquid wax is also continually pulled up into the wick via capillary action.

Both of the graphs in Figure 10 achieve organisational closure, but neither are typically considered living systems. The global water cycle certainly achieves a closure of processes, as demonstrated in W1W_{1}, but it is nonetheless maintained via externally determined constraints like the sun Mossio \BBA Bich (\APACyear2017). That said, life can actually play an active role in the maintenance of hydrological cycles by forming cloud condensation nuclei (CCN) via the release and subsequent oxidation of dimethylsulfide El-Hani \BBA Nunes-Neto (\APACyear2020). From this perspective, the coupled life-cloud system realises a closure of constraints where cloud formation is dependent on the production of dimethylsulfide by marine phytoplankton, and the marine microbiota are dependent on precipitation and runoff to return sulfur, a necessary component of their metabolism, to the ocean.

Of particular interest to \citeAel2020life is the transition from an externally-constrained abiotic hydrological cycle, such as the one described in W1W_{1}, to a life-constrained ecological system. To conduct an analysis of this transition would require a theory of dynamic pepe-graphs. That is, we would need to study how pepe-graphs change over time. This would facilitate the application of pepe-graphs to fields like developmental and evolutionary biology where systems, and the processes that constitute them, change over time. Indeed, by using our framework to investigate how and when organisational closure is realised in a sequence of pepe-graphs we could study the emergence of life from abiotic systems.

On a shorter timescale, systems can also alter their process-enablement structure to adapt to their environments. This is what \citeAmoreno2015biological refer to as regulation. A good example of this is the regulation of lactose metabolism in E. coli via activation of the lac operon. When the concentration of lactose is low in the environment, E. coli metabolises glucose. However, if the concentration of glucose drops and lactose is present, E. coli will activate the lac operon and switch to lactose metabolism to continue producing energy. Thus, living systems like E. coli have the ability to fundamentally restructure the processes within themselves in order to maintain organisational closure.

Candles, on the other hand, have no ability to regulate their processes. If there is a gust of wind, or if the wax in a candle is running low, the candle has no lac operon to switch on to stay alight: the flame will simply extinguish. Thus, candles realise organisational closure, but have no second-order ability to regulate this closure when it is threatened. Again, to better study regulation, and by extension biological agency and autonomy, we would need to develop a dynamic theory of pepe-graphs in more detail, but this lies beyond the scope of this paper Moreno \BBA Mossio (\APACyear2015); Jaeger \BOthers. (\APACyear2024).

5 Conclusions

In order to develop robust theories of life we will need to compare, contrast, and build upon the models that have already been developed Cornish-Bowden \BBA Cárdenas (\APACyear2020). In this paper we have developed a tool, the process-enablement graph or pepe-graph, which can express any model of a physical system in an abstract and rigorous form. This allows us to compare and contrast vastly different models of a system using a common language. In particular, we focus on the cycles within pepe-graphs as these represent the parts of the system that are self-organising. Since all biological systems must continuously reproduce the physical conditions that enable their own existence, cycles of processes and enablements represent the lifelike features of a system Nicholson \BBA Dupré (\APACyear2018); Moreno \BBA Mossio (\APACyear2015); Montévil \BBA Mossio (\APACyear2015). We develop the mathematical theory needed to compare pepe-graphs and the cycles found within them by defining pepe-graph homomorphisms and homorheisms. These two mathematical tools allow us to compare pepe-graphs, and therefore models of a system, in a rigorous and explicit manner.

One advantage of our framework is that pepe-graphs can be used to analyse any model of a physical system. This is not to say that every model of a system is worthy of consideration, nor that every pepe-graph is accurate. Some models may describe the phenomena within a system more accurately than others, and drawing a good pepe-graph requires an intimate knowledge of the system under investigation. But importantly, there are no ‘fundamentally true’ models: every model of a system is imperfect in its own way Massimi (\APACyear2022). Thus, to improve our understanding of a system we can never look at a single model, rather we need to take many different models into consideration and see how they compare. We have conducted such an analysis in this paper by comparing autopoiesis to autocatalytic sets and (F,A)(F,A)-systems using pepe-graphs Maturana \BBA Varela (\APACyear1980); Kauffman (\APACyear1986); Hofmeyr (\APACyear2021).

As a final point, we would like to highlight that our framework of pepe-graphs and the homomorphisms between them is a kind of universal biology Mariscal \BBA Fleming (\APACyear2017). To do universal biology is to study life as it must be, rather than studying life as it exists in this or that system. In this light we argue that all life must be self-organising Jaeger \BOthers. (\APACyear2024); Rosen (\APACyear1991). Nonetheless, our approach is not exclusively focused on living matter, nor does it give a clear demarcation between biotic and abiotic systems. This is a significant advantage as it can help us to understand ‘life’s edge’ and to study the transitions between life and nonlife Maynard Smith \BBA Szathmary (\APACyear1997). That is, our framework allows us to study how self-organisation may arise in any system. Thus, pepe-graphs are a powerful tool both to advance our theories of life and to better understand self-organising systems.

Acknowledgements

EB would like to thank Léo Diaz (University of Melbourne) and Adriana Zanca (University of Melbourne) for insightful discussions.

References

  • Ashe \BOthers. (\APACyear2021) \APACinsertmetastarashe2021does{APACrefauthors}Ashe, A., Colot, V.\BCBL \BBA Oldroyd, B\BPBIP.  \APACrefYearMonthDay2021. \BBOQ\APACrefatitleHow does epigenetics influence the course of evolution? How does epigenetics influence the course of evolution?\BBCQ \APACjournalVolNumPagesPhilosophical Transactions of the Royal Society B376182620200111. {APACrefDOI} 10.1098/rstb.2020.0111 \PrintBackRefs\CurrentBib
  • Chang (\APACyear2012) \APACinsertmetastarchang2012water{APACrefauthors}Chang, H.  \APACrefYear2012. \APACrefbtitleIs water H2O?: Evidence, realism and pluralism Is water H2O?: Evidence, realism and pluralism. \APACaddressPublisherSpringer Dordrecht. {APACrefDOI} 10.1007/978-94-007-3932-1 \PrintBackRefs\CurrentBib
  • Chiu (\APACyear2022) \APACinsertmetastarchiu2022extended{APACrefauthors}Chiu, L.  \APACrefYear2022. \APACrefbtitleExtended Evolutionary Synthesis: a Review of the Latest Scientific Research Extended evolutionary synthesis: a review of the latest scientific research. \APACaddressPublisherJohn Templeton Foundation. {APACrefDOI} 10.15868/socialsector.40950 \PrintBackRefs\CurrentBib
  • Cleland (\APACyear2019) \APACinsertmetastarcleland2019quest{APACrefauthors}Cleland, C\BPBIE.  \APACrefYear2019. \APACrefbtitleThe Quest for a Universal Theory of Life: searching for Life as we don’t know it The quest for a universal theory of life: searching for life as we don’t know it. \APACaddressPublisherCambridge University Press. {APACrefDOI} 10.1017/9781139046893 \PrintBackRefs\CurrentBib
  • Cornish-Bowden \BBA Cárdenas (\APACyear2020) \APACinsertmetastarcornish2020contrasting{APACrefauthors}Cornish-Bowden, A.\BCBT \BBA Cárdenas, M\BPBIL.  \APACrefYearMonthDay2020. \BBOQ\APACrefatitleContrasting theories of life: Historical context, current theories. In search of an ideal theory Contrasting theories of life: Historical context, current theories. In search of an ideal theory.\BBCQ \APACjournalVolNumPagesBiosystems188104063. {APACrefDOI} 10.1016/j.biosystems.2019.104063 \PrintBackRefs\CurrentBib
  • Cusimano \BBA Sterner (\APACyear2020) \APACinsertmetastarcusimano2020objectivity{APACrefauthors}Cusimano, S.\BCBT \BBA Sterner, B.  \APACrefYearMonthDay2020. \BBOQ\APACrefatitleThe objectivity of organizational functions The objectivity of organizational functions.\BBCQ \APACjournalVolNumPagesActa Biotheoretica682253–269. {APACrefDOI} 10.1007/s10441-019-09365-9 \PrintBackRefs\CurrentBib
  • Dyson (\APACyear1982) \APACinsertmetastardyson1982model{APACrefauthors}Dyson, F\BPBIJ.  \APACrefYearMonthDay1982. \BBOQ\APACrefatitleA model for the origin of life A model for the origin of life.\BBCQ \APACjournalVolNumPagesJournal of Molecular Evolution18344–350. {APACrefDOI} 10.1007/BF01733901 \PrintBackRefs\CurrentBib
  • El-Hani \BBA Nunes-Neto (\APACyear2020) \APACinsertmetastarel2020life{APACrefauthors}El-Hani, C\BPBIN.\BCBT \BBA Nunes-Neto, N.  \APACrefYearMonthDay2020. \BBOQ\APACrefatitleLife on earth is not a passenger, but a driver: explaining the transition from a physicochemical to a life-constrained world from an organizational perspective Life on earth is not a passenger, but a driver: explaining the transition from a physicochemical to a life-constrained world from an organizational perspective.\BBCQ \BIn L. Baravalle \BBA L. Zaterka (\BEDS), \APACrefbtitleLife and Evolution: Latin American Essays on the History and Philosophy of Biology Life and evolution: Latin American essays on the history and philosophy of biology (\BPGS 69–84). \APACaddressPublisherSpringer Cham. {APACrefDOI} 10.1007/978-3-030-39589-6_5 \PrintBackRefs\CurrentBib
  • Forterre (\APACyear2011) \APACinsertmetastarforterre2011manipulation{APACrefauthors}Forterre, P.  \APACrefYearMonthDay2011. \BBOQ\APACrefatitleManipulation of cellular syntheses and the nature of viruses: the virocell concept Manipulation of cellular syntheses and the nature of viruses: the virocell concept.\BBCQ \APACjournalVolNumPagesComptes Rendus Chimie144392–399. {APACrefDOI} 10.1016/j.crci.2010.06.007 \PrintBackRefs\CurrentBib
  • Forterre (\APACyear2016) \APACinsertmetastarforterre2016or{APACrefauthors}Forterre, P.  \APACrefYearMonthDay2016. \BBOQ\APACrefatitleTo be or not to be alive: How recent discoveries challenge the traditional definitions of viruses and life To be or not to be alive: How recent discoveries challenge the traditional definitions of viruses and life.\BBCQ \APACjournalVolNumPagesStudies in History and Philosophy of Science Part C: Studies in History and Philosophy of Biological and Biomedical Sciences59100–108. {APACrefDOI} 10.1016/j.shpsc.2016.02.013 \PrintBackRefs\CurrentBib
  • Frank (\APACyear2012) \APACinsertmetastarfrank2012natural{APACrefauthors}Frank, S\BPBIA.  \APACrefYearMonthDay2012. \BBOQ\APACrefatitleNatural selection. IV. The price equation Natural selection. IV. The price equation.\BBCQ \APACjournalVolNumPagesJournal of Evolutionary Biology2561002–1019. {APACrefDOI} 10.1111/j.1420-9101.2012.02498.x \PrintBackRefs\CurrentBib
  • Gánti (\APACyear2003) \APACinsertmetastarganti2003principles{APACrefauthors}Gánti, T.  \APACrefYear2003. \APACrefbtitleThe Principles of Life The principles of life. \APACaddressPublisherOxford University Press. {APACrefDOI} 10.1093/acprof:oso/9780198507260.001.0001 \PrintBackRefs\CurrentBib
  • Griesemer (\APACyear2006) \APACinsertmetastarGriesemer+2006+199+237{APACrefauthors}Griesemer, J.  \APACrefYearMonthDay2006. \BBOQ\APACrefatitleGenetics from an Evolutionary Process Perspective Genetics from an evolutionary process perspective.\BBCQ \BIn E\BPBIM. Neumann-Held \BBA C. Rehmann-Sutter (\BEDS), \APACrefbtitleGenes in Development: re-reading the molecular paradigm Genes in development: re-reading the molecular paradigm (\BPGS 199–237). \APACaddressPublisherDuke University Press. {APACrefDOI} 10.1215/9780822387336-009 \PrintBackRefs\CurrentBib
  • Gyllingberg \BOthers. (\APACyear2023) \APACinsertmetastargyllingberg2023lost{APACrefauthors}Gyllingberg, L., Birhane, A.\BCBL \BBA Sumpter, D\BPBIJ.  \APACrefYearMonthDay2023. \BBOQ\APACrefatitleThe lost art of mathematical modelling The lost art of mathematical modelling.\BBCQ \APACjournalVolNumPagesMathematical Biosciences362109033. {APACrefDOI} 10.1016/j.mbs.2023.109033 \PrintBackRefs\CurrentBib
  • Hell \BBA Nesetril (\APACyear2004) \APACinsertmetastarhell2004graphs{APACrefauthors}Hell, P.\BCBT \BBA Nesetril, J.  \APACrefYear2004. \APACrefbtitleGraphs and Homomorphisms Graphs and homomorphisms. \APACaddressPublisherOxford University Press. {APACrefDOI} 10.1093/acprof:oso/9780198528173.001.0001 \PrintBackRefs\CurrentBib
  • Hofmeyr (\APACyear2019) \APACinsertmetastarhofmeyr2019basic{APACrefauthors}Hofmeyr, J\BHBIH\BPBIS.  \APACrefYearMonthDay2019. \BBOQ\APACrefatitleBasic biological anticipation Basic biological anticipation.\BBCQ \BIn R. Poli (\BED), \APACrefbtitleHandbook of anticipation: theoretical and Applied Aspects of the Use of Future in Decision Making Handbook of anticipation: theoretical and applied aspects of the use of future in decision making (\BPGS 219–233). \APACaddressPublisherSpringer Cham. {APACrefDOI} 10.1007/978-3-319-91554-8_51 \PrintBackRefs\CurrentBib
  • Hofmeyr (\APACyear2021) \APACinsertmetastarhofmeyr2021biochemically{APACrefauthors}Hofmeyr, J\BHBIH\BPBIS.  \APACrefYearMonthDay2021. \BBOQ\APACrefatitleA biochemically-realisable relational model of the self-manufacturing cell A biochemically-realisable relational model of the self-manufacturing cell.\BBCQ \APACjournalVolNumPagesBiosystems207104463. {APACrefDOI} 10.1016/j.biosystems.2021.104463 \PrintBackRefs\CurrentBib
  • Hordijk \BBA Steel (\APACyear2004) \APACinsertmetastarhordijk2004detecting{APACrefauthors}Hordijk, W.\BCBT \BBA Steel, M.  \APACrefYearMonthDay2004. \BBOQ\APACrefatitleDetecting autocatalytic, self-sustaining sets in chemical reaction systems Detecting autocatalytic, self-sustaining sets in chemical reaction systems.\BBCQ \APACjournalVolNumPagesJournal of Theoretical Biology2274451–461. {APACrefDOI} 10.1016/j.jtbi.2003.11.020 \PrintBackRefs\CurrentBib
  • Hordijk \BOthers. (\APACyear2012) \APACinsertmetastarhordijk2012structure{APACrefauthors}Hordijk, W., Steel, M.\BCBL \BBA Kauffman, S.  \APACrefYearMonthDay2012. \BBOQ\APACrefatitleThe structure of autocatalytic sets: Evolvability, enablement, and emergence The structure of autocatalytic sets: Evolvability, enablement, and emergence.\BBCQ \APACjournalVolNumPagesActa Biotheoretica60379–392. {APACrefDOI} 10.1007/s10441-012-9165-1 \PrintBackRefs\CurrentBib
  • Jablonka \BBA Raz (\APACyear2009) \APACinsertmetastarjablonka2009transgenerational{APACrefauthors}Jablonka, E.\BCBT \BBA Raz, G.  \APACrefYearMonthDay2009. \BBOQ\APACrefatitleTransgenerational epigenetic inheritance: prevalence, mechanisms, and implications for the study of heredity and evolution Transgenerational epigenetic inheritance: prevalence, mechanisms, and implications for the study of heredity and evolution.\BBCQ \APACjournalVolNumPagesThe Quarterly Review of Biology842131–176. {APACrefDOI} 10.1086/598822 \PrintBackRefs\CurrentBib
  • Jaeger \BOthers. (\APACyear2024) \APACinsertmetastarjaegernaturalizing{APACrefauthors}Jaeger, J., Riedl, A., Djedovic, A., Vervaeke, J.\BCBL \BBA Walsh, D.  \APACrefYearMonthDay2024. \BBOQ\APACrefatitleNaturalizing relevance realization: why agency and cognition are fundamentally not computational Naturalizing relevance realization: why agency and cognition are fundamentally not computational.\BBCQ \APACjournalVolNumPagesFrontiers in Psychology151362658. {APACrefDOI} 10.3389/fpsyg.2024.1362658 \PrintBackRefs\CurrentBib
  • Jeckelmann \BBA Erni (\APACyear2020) \APACinsertmetastarjeckelmann2020transporters{APACrefauthors}Jeckelmann, J\BHBIM.\BCBT \BBA Erni, B.  \APACrefYearMonthDay2020. \BBOQ\APACrefatitleTransporters of glucose and other carbohydrates in bacteria Transporters of glucose and other carbohydrates in bacteria.\BBCQ \APACjournalVolNumPagesPflügers Archiv-European Journal of Physiology4721129–1153. {APACrefDOI} 10.1007/s00424-020-02379-0 \PrintBackRefs\CurrentBib
  • Juarrero (\APACyear2023) \APACinsertmetastarjuarrero2023context{APACrefauthors}Juarrero, A.  \APACrefYear2023. \APACrefbtitleContext changes everything: How constraints create coherence Context changes everything: How constraints create coherence. \APACaddressPublisherMIT Press. {APACrefDOI} 10.7551/mitpress/14630.001.0001 \PrintBackRefs\CurrentBib
  • Kauffman (\APACyear1986) \APACinsertmetastarkauffman1986autocatalytic{APACrefauthors}Kauffman, S\BPBIA.  \APACrefYearMonthDay1986. \BBOQ\APACrefatitleAutocatalytic sets of proteins Autocatalytic sets of proteins.\BBCQ \APACjournalVolNumPagesJournal of Theoretical Biology11911–24. {APACrefDOI} 10.1016/S0022-5193(86)80047-9 \PrintBackRefs\CurrentBib
  • Kauffman (\APACyear2000) \APACinsertmetastarkauffman2000investigations{APACrefauthors}Kauffman, S\BPBIA.  \APACrefYear2000. \APACrefbtitleInvestigations Investigations. \APACaddressPublisherOxford University Press. \PrintBackRefs\CurrentBib
  • Kauffman (\APACyear2019) \APACinsertmetastarkauffman2019world{APACrefauthors}Kauffman, S\BPBIA.  \APACrefYear2019. \APACrefbtitleA world beyond physics: the emergence and evolution of life A world beyond physics: the emergence and evolution of life. \APACaddressPublisherOxford University Press. \PrintBackRefs\CurrentBib
  • Lapid (\APACyear2003) \APACinsertmetastarlapid2003through{APACrefauthors}Lapid, Y.  \APACrefYearMonthDay2003. \BBOQ\APACrefatitleThrough dialogue to engaged pluralism: the unfinished business of the third debate Through dialogue to engaged pluralism: the unfinished business of the third debate.\BBCQ \APACjournalVolNumPagesInternational Studies Review51128–131. {APACrefDOI} 10.1111/1521-9488.501019_3 \PrintBackRefs\CurrentBib
  • Lewontin (\APACyear1983) \APACinsertmetastarlewontin1983organism{APACrefauthors}Lewontin, R\BPBIC.  \APACrefYearMonthDay1983. \BBOQ\APACrefatitleThe organism as the subject and object of evolution The organism as the subject and object of evolution.\BBCQ \APACjournalVolNumPagesScientia771865. \PrintBackRefs\CurrentBib
  • Longo \BBA Montévil (\APACyear2013) \APACinsertmetastarlongo2013extended{APACrefauthors}Longo, G.\BCBT \BBA Montévil, M.  \APACrefYearMonthDay2013. \BBOQ\APACrefatitleExtended criticality, phase spaces and enablement in biology Extended criticality, phase spaces and enablement in biology.\BBCQ \APACjournalVolNumPagesChaos, Solitons & Fractals5564–79. {APACrefDOI} 10.1016/j.chaos.2013.03.008 \PrintBackRefs\CurrentBib
  • Mariscal \BBA Fleming (\APACyear2017) \APACinsertmetastarmariscal2018we{APACrefauthors}Mariscal, C.\BCBT \BBA Fleming, L.  \APACrefYearMonthDay2017. \BBOQ\APACrefatitleWhy we should care about universal biology Why we should care about universal biology.\BBCQ \APACjournalVolNumPagesBiological Theory13121–130. {APACrefDOI} 10.1007/s13752-017-0280-8 \PrintBackRefs\CurrentBib
  • Massimi (\APACyear2018) \APACinsertmetastarmassimi2018four{APACrefauthors}Massimi, M.  \APACrefYearMonthDay2018. \BBOQ\APACrefatitleFour kinds of perspectival truth Four kinds of perspectival truth.\BBCQ \APACjournalVolNumPagesPhilosophy and Phenomenological Research962342–359. {APACrefDOI} 10.1111/phpr.12300 \PrintBackRefs\CurrentBib
  • Massimi (\APACyear2022) \APACinsertmetastarmassimi2022perspectival{APACrefauthors}Massimi, M.  \APACrefYear2022. \APACrefbtitlePerspectival realism Perspectival realism. \APACaddressPublisherOxford University Press. {APACrefDOI} 10.1093/oso/9780197555620.001.0001 \PrintBackRefs\CurrentBib
  • Massimi (\APACyear2023) \APACinsertmetastarmassimi2023epistemic{APACrefauthors}Massimi, M.  \APACrefYearMonthDay2023. \BBOQ\APACrefatitleEpistemic communities and their situated practices: perspectival realism—a primer Epistemic communities and their situated practices: perspectival realism—a primer.\BBCQ \APACjournalVolNumPagesAnnals of the New York Academy of Sciences152315–10. {APACrefDOI} 10.1111/nyas.14981 \PrintBackRefs\CurrentBib
  • Maturana \BBA Varela (\APACyear1980) \APACinsertmetastarmaturana1980autopoiesis{APACrefauthors}Maturana, H\BPBIR.\BCBT \BBA Varela, F.  \APACrefYear1980. \APACrefbtitleAutopoiesis and cognition: the realization of the living Autopoiesis and cognition: the realization of the living. \APACaddressPublisherSpringer Dordrecht. {APACrefDOI} 10.1007/978-94-009-8947-4 \PrintBackRefs\CurrentBib
  • Maynard Smith \BBA Szathmary (\APACyear1997) \APACinsertmetastarsmith1997major{APACrefauthors}Maynard Smith, J.\BCBT \BBA Szathmary, E.  \APACrefYear1997. \APACrefbtitleThe major transitions in evolution The major transitions in evolution. \APACaddressPublisherOxford University Press. {APACrefDOI} 10.1093/oso/9780198502944.001.0001 \PrintBackRefs\CurrentBib
  • McGuigan \BOthers. (\APACyear2021) \APACinsertmetastarmcguigan2021epigenetics{APACrefauthors}McGuigan, K., Hoffmann, A\BPBIA.\BCBL \BBA Sgrò, C\BPBIM.  \APACrefYearMonthDay2021. \BBOQ\APACrefatitleHow is epigenetics predicted to contribute to climate change adaptation? What evidence do we need? How is epigenetics predicted to contribute to climate change adaptation? What evidence do we need?\BBCQ \APACjournalVolNumPagesPhilosophical Transactions of the Royal Society B376182620200119. {APACrefDOI} 10.1098/rstb.2020.0119 \PrintBackRefs\CurrentBib
  • Meincke (\APACyear2022) \APACinsertmetastarmeincke2022one{APACrefauthors}Meincke, A\BPBIS.  \APACrefYearMonthDay2022. \BBOQ\APACrefatitleOne or two? A Process view of pregnancy One or two? A process view of pregnancy.\BBCQ \APACjournalVolNumPagesPhilosophical Studies17951495–1521. {APACrefDOI} 10.1007/s11098-021-01716-y \PrintBackRefs\CurrentBib
  • Mitchell (\APACyear2002) \APACinsertmetastarmitchell2002integrative{APACrefauthors}Mitchell, S\BPBID.  \APACrefYearMonthDay2002. \BBOQ\APACrefatitleIntegrative pluralism Integrative pluralism.\BBCQ \APACjournalVolNumPagesBiology and Philosophy1755–70. {APACrefDOI} 10.1023/A:1012990030867 \PrintBackRefs\CurrentBib
  • Mitchell \BBA Dietrich (\APACyear2006) \APACinsertmetastarmitchell2006integration{APACrefauthors}Mitchell, S\BPBID.\BCBT \BBA Dietrich, M\BPBIR.  \APACrefYearMonthDay2006. \BBOQ\APACrefatitleIntegration without unification: an argument for pluralism in the biological sciences Integration without unification: an argument for pluralism in the biological sciences.\BBCQ \APACjournalVolNumPagesThe American Naturalist168S6S73–S79. {APACrefDOI} 10.1086/509050 \PrintBackRefs\CurrentBib
  • Montévil \BBA Mossio (\APACyear2015) \APACinsertmetastarmontevil2015biological{APACrefauthors}Montévil, M.\BCBT \BBA Mossio, M.  \APACrefYearMonthDay2015. \BBOQ\APACrefatitleBiological organisation as closure of constraints Biological organisation as closure of constraints.\BBCQ \APACjournalVolNumPagesJournal of Theoretical Biology372179–191. {APACrefDOI} 10.1016/j.jtbi.2015.02.029 \PrintBackRefs\CurrentBib
  • Moreno \BBA Mossio (\APACyear2015) \APACinsertmetastarmoreno2015biological{APACrefauthors}Moreno, A.\BCBT \BBA Mossio, M.  \APACrefYear2015. \APACrefbtitleBiological autonomy: a Philosophical and Theoretical Enquiry Biological autonomy: a philosophical and theoretical enquiry. \APACaddressPublisherSpringer Dordrecht. {APACrefDOI} 10.1007/978-94-017-9837-2 \PrintBackRefs\CurrentBib
  • Mossio \BBA Bich (\APACyear2017) \APACinsertmetastarmossio2017makes{APACrefauthors}Mossio, M.\BCBT \BBA Bich, L.  \APACrefYearMonthDay2017. \BBOQ\APACrefatitleWhat makes biological organisation teleological? What makes biological organisation teleological?\BBCQ \APACjournalVolNumPagesSynthese19441089–1114. {APACrefDOI} 10.1007/s11229-014-0594-z \PrintBackRefs\CurrentBib
  • Nahas \BBA Sachs (\APACyear2023) \APACinsertmetastarnahas2023s{APACrefauthors}Nahas, A.\BCBT \BBA Sachs, C.  \APACrefYearMonthDay2023. \BBOQ\APACrefatitleWhat’s at stake in the debate over naturalizing teleology? An overlooked metatheoretical debate What’s at stake in the debate over naturalizing teleology? An overlooked metatheoretical debate.\BBCQ \APACjournalVolNumPagesSynthese2014142. {APACrefDOI} 10.1007/s11229-023-04147-w \PrintBackRefs\CurrentBib
  • Nicholson (\APACyear2019) \APACinsertmetastarnicholson2019cell{APACrefauthors}Nicholson, D\BPBIJ.  \APACrefYearMonthDay2019. \BBOQ\APACrefatitleIs the cell really a machine? Is the cell really a machine?\BBCQ \APACjournalVolNumPagesJournal of Theoretical Biology477108–126. {APACrefDOI} 10.1016/j.jtbi.2019.06.002 \PrintBackRefs\CurrentBib
  • Nicholson \BBA Dupré (\APACyear2018) \APACinsertmetastarnicholson2018everything{APACrefauthors}Nicholson, D\BPBIJ.\BCBT \BBA Dupré, J.  \APACrefYear2018. \APACrefbtitleEverything flows: Towards a processual philosophy of biology Everything flows: Towards a processual philosophy of biology. \APACaddressPublisherOxford University Press. {APACrefDOI} 10.1093/oso/9780198779636.001.0001 \PrintBackRefs\CurrentBib
  • Oyama (\APACyear2000) \APACinsertmetastaroyama2000ontogeny{APACrefauthors}Oyama, S.  \APACrefYear2000. \APACrefbtitleThe ontogeny of information: developmental systems and evolution The ontogeny of information: developmental systems and evolution. \APACaddressPublisherDuke University Press. {APACrefDOI} 10.1215/9780822380665 \PrintBackRefs\CurrentBib
  • Oyama \BOthers. (\APACyear2001) \APACinsertmetastaroyama2003cycles{APACrefauthors}Oyama, S., Gray, R\BPBID.\BCBL \BBA Griffiths, P\BPBIE.  \APACrefYear2001. \APACrefbtitleCycles of contingency: developmental systems and evolution Cycles of contingency: developmental systems and evolution. \APACaddressPublisherMIT Press. \PrintBackRefs\CurrentBib
  • Price (\APACyear1970) \APACinsertmetastarprice1970selection{APACrefauthors}Price, G\BPBIR.  \APACrefYearMonthDay1970. \BBOQ\APACrefatitleSelection and covariance Selection and covariance.\BBCQ \APACjournalVolNumPagesNature227520–521. {APACrefDOI} 10.1038/227520a0 \PrintBackRefs\CurrentBib
  • Rao \BBA Nanjundiah (\APACyear2011) \APACinsertmetastarrao2011jbs{APACrefauthors}Rao, V.\BCBT \BBA Nanjundiah, V.  \APACrefYearMonthDay2011. \BBOQ\APACrefatitleJ.B.S Haldane, Ernst Mayr and the beanbag genetics dispute J.B.S Haldane, Ernst Mayr and the beanbag genetics dispute.\BBCQ \APACjournalVolNumPagesJournal of the History of Biology44233–281. {APACrefDOI} 10.1007/s10739-010-9229-5 \PrintBackRefs\CurrentBib
  • Rosen (\APACyear1991) \APACinsertmetastarrosen1991life{APACrefauthors}Rosen, R.  \APACrefYear1991. \APACrefbtitleLife Itself: A Comprehensive Inquiry into the Nature, Origin, and Fabrication of Life Life itself: A comprehensive inquiry into the nature, origin, and fabrication of life. \APACaddressPublisherColumbia University Press. \PrintBackRefs\CurrentBib
  • Sanchez \BBA Lagunoff (\APACyear2015) \APACinsertmetastarsanchez2015viral{APACrefauthors}Sanchez, E\BPBIL.\BCBT \BBA Lagunoff, M.  \APACrefYearMonthDay2015. \BBOQ\APACrefatitleViral activation of cellular metabolism Viral activation of cellular metabolism.\BBCQ \APACjournalVolNumPagesVirology479–480609–618. {APACrefDOI} 10.1016/j.virol.2015.02.038 \PrintBackRefs\CurrentBib
  • Thompson (\APACyear2007) \APACinsertmetastarthompson2007mind{APACrefauthors}Thompson, E.  \APACrefYear2007. \APACrefbtitleMind in life: biology, phenomenology, and the sciences of mind Mind in life: biology, phenomenology, and the sciences of mind. \APACaddressPublisherHarvard University Press. \PrintBackRefs\CurrentBib
  • Turing (\APACyear1936) \APACinsertmetastarturing1936computable{APACrefauthors}Turing, A\BPBIM.  \APACrefYearMonthDay1936. \BBOQ\APACrefatitleOn computable numbers, with an application to the Entscheidungsproblem On computable numbers, with an application to the Entscheidungsproblem.\BBCQ \APACjournalVolNumPagesProceedings of the London Mathematical Societys2–42230–265. {APACrefDOI} 10.1112/plms/s2-42.1.230 \PrintBackRefs\CurrentBib
  • van Hemmen (\APACyear2007) \APACinsertmetastarvan2007biology{APACrefauthors}van Hemmen, J\BPBIL.  \APACrefYearMonthDay2007. \BBOQ\APACrefatitleBiology and mathematics: a fruitful merger of two cultures Biology and mathematics: a fruitful merger of two cultures.\BBCQ \APACjournalVolNumPagesBiological Cybernetics9711–3. {APACrefDOI} 10.1007/s00422-007-0163-3 \PrintBackRefs\CurrentBib
  • Varela (\APACyear2000) \APACinsertmetastarvarela2000fenomeno{APACrefauthors}Varela, F.  \APACrefYear2000. \APACrefbtitleEl fenómeno de la vida El fenómeno de la vida. \APACaddressPublisherDolmen Ediciones. \PrintBackRefs\CurrentBib
  • Vittadello \BBA Stumpf (\APACyear2021) \APACinsertmetastarvittadello2021model{APACrefauthors}Vittadello, S\BPBIT.\BCBT \BBA Stumpf, M\BPBIP.  \APACrefYearMonthDay2021. \BBOQ\APACrefatitleModel comparison via simplicial complexes and persistent homology Model comparison via simplicial complexes and persistent homology.\BBCQ \APACjournalVolNumPagesRoyal Society Open Science810211361. {APACrefDOI} 10.1098/rsos.211361 \PrintBackRefs\CurrentBib
  • Vittadello \BBA Stumpf (\APACyear2022\APACexlab\BCnt1) \APACinsertmetastarvittadello2022group{APACrefauthors}Vittadello, S\BPBIT.\BCBT \BBA Stumpf, M\BPBIP.  \APACrefYearMonthDay2022\BCnt1. \BBOQ\APACrefatitleA group theoretic approach to model comparison with simplicial representations A group theoretic approach to model comparison with simplicial representations.\BBCQ \APACjournalVolNumPagesJournal of Mathematical Biology85548. {APACrefDOI} 10.1007/s00285-022-01807-2 \PrintBackRefs\CurrentBib
  • Vittadello \BBA Stumpf (\APACyear2022\APACexlab\BCnt2) \APACinsertmetastarvittadello2022open{APACrefauthors}Vittadello, S\BPBIT.\BCBT \BBA Stumpf, M\BPBIP.  \APACrefYearMonthDay2022\BCnt2. \BBOQ\APACrefatitleOpen problems in Mathematical Biology Open problems in mathematical biology.\BBCQ \APACjournalVolNumPagesMathematical Biosciences354108926. {APACrefDOI} 10.1016/j.mbs.2022.108926 \PrintBackRefs\CurrentBib
  • Waddington (\APACyear1957) \APACinsertmetastarwaddington1957strategy{APACrefauthors}Waddington, C\BPBIH.  \APACrefYear1957. \APACrefbtitleThe Strategy of the Genes The strategy of the genes. \APACaddressPublisherRoutledge. \PrintBackRefs\CurrentBib
  • Walsh (\APACyear2018) \APACinsertmetastarwalsh2018objectcy{APACrefauthors}Walsh, D\BPBIM.  \APACrefYearMonthDay2018. \BBOQ\APACrefatitleObjectcy and agency: Towards a methodological vitalism Objectcy and agency: Towards a methodological vitalism.\BBCQ \BIn D\BPBIJ. Nicholson \BBA J. Dupré (\BEDS), \APACrefbtitleEverything Flows Everything flows (\BPGS 167–185). \APACaddressPublisherOxford University Press. {APACrefDOI} 10.1093/oso/9780198779636.003.0008 \PrintBackRefs\CurrentBib
  • Wimsatt (\APACyear2007) \APACinsertmetastarwimsatt2007re{APACrefauthors}Wimsatt, W\BPBIC.  \APACrefYear2007. \APACrefbtitleRe-engineering philosophy for limited beings: Piecewise approximations to reality Re-engineering philosophy for limited beings: Piecewise approximations to reality. \APACaddressPublisherHarvard University Press. {APACrefDOI} 10.2307/j.ctv1pncnrh \PrintBackRefs\CurrentBib
  • Yang \BOthers. (\APACyear2003) \APACinsertmetastaryang2003protein{APACrefauthors}Yang, H., Luo, G., Karnchanaphanurach, P., Louie, T\BHBIM., Rech, I., Cova, S.\BDBLXie, X\BPBIS.  \APACrefYearMonthDay2003. \BBOQ\APACrefatitleProtein conformational dynamics probed by single-molecule electron transfer Protein conformational dynamics probed by single-molecule electron transfer.\BBCQ \APACjournalVolNumPagesScience3025643262–266. {APACrefDOI} 10.1126/science.1086911 \PrintBackRefs\CurrentBib

Appendix A Appendix

In this section we discuss in greater detail the pepe-graph I𝒫I_{\mathcal{P}} from Figure 8 in Section 4.3, and prove that the maps ϕ5:I𝒫FA\phi_{5}\colon I_{\mathcal{P}}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}F_{A} and ϕ6:I𝒫A3\phi_{6}\colon I_{\mathcal{P}}\mathrel{\vphantom{\rightarrow}\mathchoice{\ooalign{$\displaystyle\circ\mkern 1.0mu$\cr$\displaystyle\longrightarrow$\cr}}{\ooalign{$\textstyle\circ\mkern 1.0mu$\cr$\textstyle\longrightarrow$\cr}}{\ooalign{$\scriptstyle\circ\mkern 1.0mu$\cr$\scriptstyle\longrightarrow$\cr}}{\ooalign{$\scriptscriptstyle\circ\mkern 1.0mu$\cr$\scriptscriptstyle\longrightarrow$\cr}}}A_{3} are homorheisms. We abbreviate the processes described in I𝒫I_{\mathcal{P}}, FAF_{A}, and A3A_{3} through the correspondence in Table A.1.

Table A.1: Abbreviations for the processes in the pepe-graphs I𝒫I_{\mathcal{P}}, FAF_{A}, and A3A_{3} from Figure 8 in Section 4.3.
Abbreviation Process
CC Covalent Chemistry
PF Protein Folding
IT Ion Transport
MM Membrane Maintenance
NT Nutrient Transport
F Fabrication
A Assembly
I Maintenance of the intracellular milieu
M Metabolic Reaction Network
B Boundary Maintenance
T Transport

We first confirm that each of the edges in I𝒫I_{\mathcal{P}} is a direct enablement:

  • NT \to CC. (Enablement): Once nutrients enter the cell via NT, they are metabolised through the process of CC. The latter would not occur without the former as there would be no reactants for all of the metabolic reactions in CC. (Direct): As nutrients enter the cell they immediately begin interacting with enzymes and other molecules.

  • CC \to MM. (Enablement): CC produces the lipids that are necessary to form the cell membrane. (Direct): Membrane maintenance includes the trafficking of lipids to the membrane. As lipids are synthesised within the cytoplasm they immediately begin to self-assemble and commence their journey to the cell membrane.

  • CC \to PF. (Enablement): CC produces the polypeptides that are needed for protein folding to occur. (Direct): As polypeptides are produced by the ribosome they interact with the local intracellular environment and begin to take on secondary structural features through the process of cotranslational folding.

  • PF \to CC. (Enablement): Folded enzymes are needed to catalyze the reactions in CC. (Direct): Protein folding is a continuous process, since proteins shift between conformations as they diffuse around the cell and participate in metabolic reactions.

  • PF \to IT. (Enablement): Folded ion transporters are needed to transport ions across the membrane. (Direct): Membrane transporters constantly have to maintain their conformation as they interact with the intracellular milieu, the extracellular environment, and the hydrophobic interior of the cell membrane.

  • PF \to NT. (Enablement): Folded membrane transporters (like GLUTs) are needed to transport nutrients across the membrane. (Direct): As for PF \to IT.

  • PF \to MM. (Enablement): Folded enzymes (like flippase and scramblase) are needed to keep the cell membrane regulated and intact. (Direct): As for PF \to IT.

  • MM \to PF. (Enablement): The maintenance of the membrane ensures that ions remain within the cell, thus controlling the composition of the intracellular milieu which mediates PF. Without MM, a protein’s folding environment would be wildly unregulated, likely causing it to fold incorrectly. (Direct): Proteins that fold near the boundary of the cell will interact with ions that have just bounced off the cell membrane.

  • MM \to NT. (Enablement): Without a membrane there is no distinction between the interior and exterior of the cell and therefore transport cannot occur. (Direct): Membrane transport proteins need to exist within the lipid bilayer in order to function properly.

  • MM \to IT. (Enablement): As for MM \to NT. (Direct): As for MM \to NT.

  • MM \to CC. (Enablement): The membrane ensures that all reactants and catalysts involved in CC remain close enough together to form a functional metabolic reaction network. (Direct): The bouncing of particular molecules against the boundary of the cell membrane causes reactants and catalysts to interact.

  • IT \to PF. (Enablement): Ions are needed to ensure a proper environment for protein folding. (Direct): Some protein folding occurs near ion channels, allowing for some interaction between these two processes there.

We will now show that ϕ5\phi_{5} and ϕ6\phi_{6} are homorheisms. Recall that a homorheism is a weak graph homomorphism that preserves closure and reflects closure. Note that ϕ5:V(I𝒫)V(FA)\phi_{5}\colon V(I_{\mathcal{P}})\to V(F_{A}) is the vertex map given by CCFCC\mapsto F, NTFNT\mapsto F, PFAPF\mapsto A, MMAMM\mapsto A, and ITMIT\mapsto M; and ϕ6:V(I𝒫)V(A3)\phi_{6}\colon V(I_{\mathcal{P}})\to V(A_{3}) is the vertex map given by CCMCC\mapsto M, NTTNT\mapsto T, PFMPF\mapsto M, MMBMM\mapsto B, and ITTIT\mapsto T. Both ϕ5\phi_{5} and ϕ6\phi_{6} are weak graph homomorphisms, so it suffices to show that they preserve closure and reflect closure. For reference, all of the cycles in I𝒫I_{\mathcal{P}} are shown in Figure A.1.

CCMMC1C_{1} CCPFC2C_{2} PFITC3C_{3}
PFMMC4C_{4} PFCCMMC5C_{5} PFCCMMC6C_{6}
NTCCMMC7C_{7} NTCCPFC8C_{8} PFITMMC9C_{9}
PFCCMMNTC10C_{10} ITPFCCMMC11C_{11} NTPFCCMMC12C_{12}
NTCCMMITPFC13C_{13}
Figure A.1: The cycles of the pepe-graph I𝒫I_{\mathcal{P}}. The colour coding corresponds to the maps ϕ5\phi_{5} and ϕ6\phi_{6} described in Figure 8 in Section 4.3.
Theorem A.1.

The weak graph homomorphisms ϕ5\phi_{5} and ϕ6\phi_{6} preserve closure.

Proof.

By inspecting the colourings of each cycle in I𝒫I_{\mathcal{P}} from Figure A.1, note that the induced image of every cycle under ϕ5\phi_{5} and ϕ6\phi_{6} contains more than one vertex, except for the graph A3[ϕ6(C2)]A_{3}[\phi_{6}(C_{2})] which has the single vertex MM and a single loop, and except for the graph FA[ϕ5(C4)]F_{A}[\phi_{5}(C_{4})] which has the single vertex AA and a single loop. It therefore follows from Theorem 3.5 that ϕ5\phi_{5} and ϕ6\phi_{6} preserve closure. ∎

Theorem A.2.

The weak graph homomorphisms ϕ5\phi_{5} and ϕ6\phi_{6} reflect closure.

Proof.

Since FA[ϕ5(C13)]=FAF_{A}[\phi_{5}(C_{13})]=F_{A} and A3[ϕ6(C13)]=A3A_{3}[\phi_{6}(C_{13})]=A_{3}, it follows from Theorem 3.10 that ϕ5\phi_{5} and ϕ6\phi_{6} reflect closure. ∎