Self-Organization and Information Processing: From Basic Enzymatic Activities to Complex Adaptive Cellular Behavior

Abstract

One of the main aims of current biology is to understand the origin of the molecular organization that underlies the complex dynamic architecture of cellular life. Here, we present an overview of the main sources of biomolecular order and complexity spanning from the most elementary levels of molecular activity to the emergence of cellular systemic behaviors. First, we have addressed the dissipative self-organization, the principal source of molecular order in the cell. Intensive studies over the last four decades have demonstrated that self-organization is central to understand enzyme activity under cellular conditions, functional coordination between enzymatic reactions, the emergence of dissipative metabolic networks (DMN), and molecular rhythms. The second fundamental source of order is molecular information processing. Studies on effective connectivity based on transfer entropy (TE) have made possible the quantification in bits of biomolecular information flows in DMN. This information processing enables efficient self-regulatory control of metabolism. As a consequence of both main sources of order, systemic functional structures emerge in the cell; in fact, quantitative analyses with DMN have revealed that the basic units of life display a global enzymatic structure that seems to be an essential characteristic of the systemic functional metabolism. This global metabolic structure has been verified experimentally in both prokaryotic and eukaryotic cells. Here, we also discuss how the study of systemic DMN, using Artificial Intelligence and advanced tools of Statistic Mechanics, has shown the emergence of Hopfield-like dynamics characterized by exhibiting associative memory. We have recently confirmed this thesis by testing associative conditioning behavior in individual amoeba cells. In these Pavlovian-like experiments, several hundreds of cells could learn new systemic migratory behaviors and remember them over long periods relative to their cell cycle, forgetting them later. Such associative process seems to correspond to an epigenetic memory. The cellular capacity of learning new adaptive systemic behaviors represents a fundamental evolutionary mechanism for cell adaptation.

Introduction

The cell is a highly self-organized dynamic molecular system that constitutes the basic fundamental unit of all known life forms.

A breakthrough has been made in the knowledge of the cellular molecular components and understanding the basic mechanisms of many of their interactions. However, we still do not know many of the essential aspects of the highly self-ordered dynamics that characterize the biochemical life of cells.

All molecular components of cells exhibit complex chemical transformations shaping an extreme super complex system, where simultaneously the concentrations of millions of ions and molecules change permanently following high-ordered patterns, structured in space and in time.

The consequences of these incessant self-organized and self-regulated dynamics of reactive molecular transformations are adequate cellular growth, the development of all physiological processes to maintain its functional structures, continuous adaptation to the environment, and, lastly, mitosis. Alterations in these self-organized dynamics lead to different pathologies and, on certain occasions, the dynamic order collapses and cell dies.

The limitations in the study of cellular molecular dynamics by traditional methods have been compensated by the efforts of quantitative sciences applied to Biology, especially Systems Biology. Physical-mathematical approaches are making possible to understand the principles underlying the complex cellular molecular dynamics adequately, so concepts such as self-organized dissipative structures, self-regulation, the emergence of functional properties, systemic attractors, and information processing, are part of the new landscape of contemporary biology.

In light of these quantitative advances, each cell is not a mere molecular aggregate but a super complex dynamic system where highly self-ordered dissipative patterns, the main source of molecular order in the cell, do emerge together with self-regulatory behaviors originated by biomolecular information processing, another fundamental source of organization in the cell.

Here, we present an overview of the main sources of biomolecular order and complexity that underline the molecular dynamics of cells. These sources span from the most elementary levels of molecular activity to the emergence of systemic behaviors.

• -

  In the “The Dynamics Originated by the Molecular Turnover and the Role of the Enzymes in These Super Complex Dynamics” section, we have approached two fundamental issues: the dynamics originated by the molecular turnover of cells and the role of the enzymes in these super complex dynamics.

All molecular components of the cell are synthesized and degraded continually following sophisticated interdependent processes that defy the human intellect. This continuous molecular turnover is produced by a super complex dynamic system formed by millions of biochemical reactions, which occur continuously at every moment of cellular life. The continuous recycling and incessant chemical transformations, that encompass practically all molecules, shape a critical scenario, only in which cellular life is possible (De la Fuente, 2015).

The super complex dynamics originated by this huge molecular turnover constitute the fundamental systemic characteristic of all basic life units.

The specialists of the life sciences and those of the quantitative sciences, attracted by the biological fact, should adequately understand the fundamental molecular dynamics of the cell, and the surprising magnitude that they reach at a global level.

Numerous works have been carried out on molecular recycling, mainly of proteins (Mathieson et al., 2018; Ross et al., 2020). However, the study of recycling other molecular components has been much less addressed. Besides, most of these works have not been thoroughly considered comprehensively. Here, for the first time, we present an integral review of the molecular turnover that includes the proteome, lipidome, glycome, metabolome, transcriptome, and main cellular structures.

All reactive molecular transformations occurring in the cell are essentially chemical reactions, that is, modifications on how atoms are bonded. Practically, all these dynamics of chemical changes in the cell are mediated by enzymes, the most extraordinary macromolecular nanomachines, responsible for breaking and joining covalent bonds, the fundamental chemical bonds in biological molecules. Through their ability to decrease the bond activation energy, enzymes permanently modify the way atoms are bonded, making possible the accelerated creation and destruction of molecules. Accordingly, enzymes are fundamental and essential molecules for metabolic life. They are the main focus of attention of this review.

Without understanding the meaning and magnitude of the dynamic molecular turnover that occurs in all living cells, as well as the fundamental role of enzymes in these permanent metabolic dynamics (which means activities developed by enzymes), it is not possible to properly understand what a cell is, nor the essential elements that underlie in its self-ordered and self-regulated functionality.

• -

  Section “Dissipative Self-Organization of Enzymatic Activities, the Main Source of Molecular Organization in the Cell” addresses the dissipative self-organization of enzymatic activities, the main source of molecular organization in the cell.

Enzymes are responsible for molecular recycling processes shaping multienzyme complexes (reversible structures formed by several enzymes of a metabolic pathway), which carry out their activity with autonomy between them playing distinctive and essential roles in cellular physiology. When a multienzyme complex operates far enough from equilibrium dissipative self-organization can emerge (Goldbeter, 2007; De la Fuente, 2010, 2014).

Here, we show the main enzymatic dissipative structures: biomolecular oscillations (metabolic rhythms). Even though self-organized structures are very diverse, and we also cover other different dissipative patterns including circadian rhythms, birhythmicity, multi-stability, and spatial traveling waves.

Numerous physical-mathematical analyses of enzymatic pathways have contributed to a better understanding of the emergence of self-organized processes. They are essential to research the main source of molecular organization in the cell. Most of these functional biochemical studies have been carried out using systems of differential equations e.g., the Krebs cycle (Mogilevskaya et al., 2006), the amino acid biosynthetic pathways (Yang et al., 2005), the oxidative phosphorylation subsystem (Korzeniewski and Zoladz, 2001), the glycolytic subsystem (Bier et al., 1996), the transduction in G-protein enzyme cascade (Kass and Bray, 1996), the gene expression (Gonze et al., 2004), the ascorbate-glutathione cycle in chloroplasts (Valero et al., 2016), and the cell cycle (Tyson, 1991). On the other hand, more complex analyses have been performed using non-ordinary differential equations, for example in the yeast glycolytic pathway (De la Fuente et al., 1995). In these studies, different attractor dynamics were analyzed linked to Hopf bifurcations (De la Fuente et al., 1996a), tangent bifurcations (De la Fuente et al., 1996b), the classical period-doubling cascade preceding chaos (De la Fuente et al., 1999a), persistent behaviors (De la Fuente et al., 1998a, b, 1999b) and the multiplicity of coexisting attractors in the phase space (De la Fuente et al., 1998c; De la Fuente, 1999).

The dissipative structures that emerge in self-organized multi-enzymatic complexes constitute one of the most genuine properties of cells, and the rigorous knowledge of their nature and significance is an essential element in the comprehension of the biological fact at its most basic and elementary levels.

• -

  Dissipative metabolic networks (DMN) and the emergence of the systemic metabolic structure (SMS) are covered in the “Dissipative Metabolic Networks and the Emergence of the Systemic Metabolic Structure” section.

Enzymes form basic functional metabolic subsystems (multienzyme complexes) for turnover activities where dissipative patterns emerge. At a higher level of complexity, they shape different DMN, which originate the emergence of the SMS at a global level.

To research the functionality of the cellular metabolism (the set of activities of all enzymes) at a systemic level, taking into account the self-organized behavior, DMN was created in 1999 (De la Fuente et al., 1999b). Essentially, a DMN is an open system formed by a given set of self-organized multi enzymatic complexes interconnected by biochemical substrate fluxes and three classes of biomolecular regulatory signals: activatory (positive allosteric modulation), inhibitory (negative allosteric modulation), and all-or-nothing type (which correspond to the regulatory enzymes of covalent modulation). The SMS was observed for the first time after performing an exhaustive numerical analysis with several millions of different DMN (De la Fuente et al., 1999b, 2008). Such systemic organization is mainly defined by a small set of different self-organized multienzyme complexes which are always in active states (metabolic core), while the rest of dissipative multienzyme subsystems exhibit on-off dynamic states. Later, this global enzymatic system was verified using flux balance analysis in several prokaryotes and eukaryote cells such as Saccharomyces cerevisiae, Helicobacter pylori, and Escherichia coli (Almaas et al., 2004, 2005; Almaas, 2007).

Adenylate energy system constitutes another global process responsible for the functional enzymatic integration in the SMS at a global level. A primary systemic ratio between ATP, ADP, and AMP concentrations was proposed by Atkinson (Atkinson and Walton, 1967) to calculate the energetic cellular level, which was denominated the Adenylate Energy Charge (AEC).

Cellular energy quantification shows that almost all organisms keep their AEC within narrow values under growth conditions, more specifically between 0.7 and 0.95, despite the extreme complex fluctuations in the adenine nucleotide concentrations.

Further quantitative studies have shown that each cell is a complex non-linearly open system in which there is not a specific energy value that is conserved, but rather dynamic forms of energy change (De la Fuente et al., 2014). Specifically, intensive experimental measurements, which under growth cellular conditions, have shown that AEC changes between 0.7 and 0.95 are invariantly maintained in practically all classes of cells which seems to represent a functional integration feature common to all cellular organisms.

• -

  The “Enzymatic Information Processing, the Second Leading Source of Molecular Order in the Cell” section of this review deals with enzymatic information processing: the second fundamental source of molecular order in the cell.

Enzyme activity during molecular turnover not only originates the emergence of high self-organized dissipative patterns, structured in space and time, but is also capable of producing complex self-regulatory behavior by information processing which highly increases the functional complexity of cellular system. The sophisticated cycles of construction and destruction of the molecular components could not be accomplished adequately without this source of new information, order, and organization.

The first quantitative work on the processing of molecular information in bits, in a multi-enzymatic complex, was performed in 2012 in the yeast glycolysis (De la Fuente and Cortés, 2012), one of the most analyzed dissipative multienzyme system, which was studied using TE.

On the other hand, a dynamic structure of information processing appears in the DMN when TE is considered under systemic activity (De la Fuente et al., 2011). According to these results, the DMN behaves as a highly complex decentralized information processing super-system that generates molecular information flows between the self-organized enzymatic sets, defining flows of biochemical instructions, at each moment, that make every enzymatic activity to evolve with a precise and particular turnover pattern.

Along with these works, different biological examples are displayed in the text showing that unicellular organisms possess a complex self-regulatory behavior originated by biomolecular information processing, the second fundamental source of organization in the cell.

• -

  Finally, in the “Hopfield-Type Systemic Attractors in Dissipative Metabolic Networks. Cellular Systemic Behaviors” section, the emergence of Hopfield-type systemic attractors in DMN is analyzed. It is necessary to underline that these attractors allow the cellular metabolic structure (CMS) to behave as an individual and complete integral system. As a consequence, cellular systemic behaviors occur.

The emergence of Hopfield-type dynamics was quantitatively observed using advanced tools of Statistic Mechanics and techniques of Artificial Intelligence in complex enzymatic dynamics (De la Fuente et al., 2013). Hopfield-like dynamics are characterized by manifesting associative memory (Hopfield, 1982), and the systemic enzymatic study was the first quantitative evidence that an associative memory can emerge in an individual cell. Such memory would be a manifestation of emergent properties underlying the complex dynamics of the systemic cellular metabolic networks when dissipative enzymatic self-organization and molecular information processing act together.

In continuation of this systemic enzymatic study (De la Fuente et al., 2013), and following Pavlov’s methodological approach with dogs (Pavlov, 1927), it was observed that two different unicellular organisms (Amoeba proteus sp. and Metamoeba leningradensis sp.) showed associative learning behaviors (De la Fuente et al., 2019a). To analyze such conditioned behavior in amoebae an electric field was used as a conditioned stimulus and a specific chemotactic peptide as an unconditioned stimulus. The migratory trajectories of more than 700 amoebae under different experimental conditions were studied, and the results showed that the unicellular organisms were able to learn new systemic behaviors, which can be considered as a rudimentary form of associative memory, crucial to govern properly cellular migration. This quantitative study demonstrated, for the first time, that associative memory was also possible in unicellular organisms. Such a type of memory could be the manifestation of the emergent properties underlying the complex dynamics of the systemic metabolic networks and their Hopfield-like attractors.

One of the most important goals of contemporary biology is to understand the fundamental principles and quantitative laws governing the functional metabolic architecture of the cell. In this review, we have mainly focused on the role of dynamic enzymatic processes, the mechanisms implicated in self-organization (dissipative structures), the self-regulation of metabolic dynamics (molecular information processing), and the emergence of the cellular systemic properties (Hopfield-like dynamics). All these issues are essential to correctly understand what a cell is.

To help life scientists who are unfamiliar with quantitative sciences understand these fundamental multidisciplinary issues, we have not used mathematical formulations here. Besides, we have reinforced the fundamental ideas with a considerable number of biological examples.

Systems Biology is being increasingly used as an appropriate methodology to address dynamic metabolic processes and the higher-level properties emerging in the cell from interactions between its elementary molecular parts, forming complex self-organized and self-regulated structures, from basic enzymatic activities to complex cellular adaptive behavior.

Systems biology is fundamental to understand the functional architecture of the most complex, sophisticated, and overwhelming molecular system known in nature, the living cell.

The Dynamics Originated by the Molecular Turnover and the Role of the Enzymes in These Super Complex Dynamics

Cells are dynamic metabolic reactors in which millions of biochemical reactions, tightly interrelated and integrated into sophisticated networks, shape the most complex molecular system in nature.

All molecular components of the cell are in a dynamic state of reactive transformations. The proteome, lipidome, glycome, metabolome, and transcriptome are synthesized and degraded continually in all cell types (see Supplementary Material 01 for further details). Even outside the growth period, the intracellular macromolecular pool is dynamic, and considerable energy is expended in the continuous processes of synthesis and degradation. For instance, only the protein turnover accounts for 38–47% of the total energy produced in every cell (Lahtvee et al., 2014).

Not only molecules but also cellular structures are subjected to cycles of construction and destruction. For instance, the endoplasmic reticulum, the largest endomembrane system, exhibits a permanent turnover; the cytoskeleton is another molecular dynamic system that undergoes continuous reorganization through the synthesis and degradation of its structural elements during the cell cycle; the mitochondria are dynamic organelles that are incessantly undergoing fusion, fission, and molecular destruction; the peroxisome turnover also takes place by autophagy-related mechanism; also, cellular membranes are highly reactive dynamic structures in non-stop recycling, so, in 30 min, an active cell such as a macrophage recycles by endocytosis and exocytosis an amount of plasma membrane that equals its complete plasma membrane (see Supplementary Material 01).

From a molecular perspective, each cell is in a permanent state of self-construction and self-destruction. This is the fundamental systemic characteristic of all basic units of life. The continuous molecular recycling that defines the systemic cellular functionality relies on a super complex dynamic system formed by millions and millions of biochemical reactions which occur continuously, at every moment of the cellular life.

Enzymes are the essential functional molecules of this super complex dynamic reactor. They are responsible for most of the molecular transformations, and such enzymatic activity, when considered globally, is called cellular metabolism (Stryer et al., 2002). Almost all enzymes are proteins; however, some RNAs, called ribozymes, also exhibit catalytic activity (Stryer et al., 2002). As a piece of metabolic activity, some non-enzymatic reactions also occur, which are important for the evolution of metabolic pathways (Keller et al., 2015).

Almost all molecular synthesis and destruction processes need to be catalyzed so that they occur at adequate fast rates to sustain cellular life. These catalytic activities involve the formation and breakdown of covalent bonds, i.e., the creation or destruction of molecules. Enzymes are accelerators of biochemical reactions. They speed up the catalytic processes in different forms, all of which decrease the activation energy of bonds (Stryer et al., 2002). As a result, the rates of enzymatic processes are accelerated from thousands to million-fold, so biochemical reactions that would take years in the absence of catalysis can occur in fractions of seconds inside the cell. For instance, a very fast catalytic process is propitiated by orotidine 5′-phosphate decarboxylase which allows the decarboxylation of orotic acid, an intermediate in the biosynthesis of pyrimidine nucleotides, in milliseconds under cellular conditions, while this reaction takes place with a half-time of 78 million years under raw natural conditions (Radzicka and Wolfenden, 1995).

Metabolism can be classified into two fundamental classes: catabolic processes (the breakdown of molecules that usually results in the release of energy and smaller essential molecular parts) and anabolic processes (the synthesis of molecules such as proteins, lipids, glycans, nucleic acids, etc., from more minor molecular elements). These anabolic and catabolic reactions ensure molecular recycling in the cell.

Enzymes are the biomolecules that perform the most important biological cellular functions. Unlike the rest of the biomolecules, enzymes are not passive. They are the only ones that can carry out a catalytic activity, transforming some molecules into others and developing complex interrelated reactive dynamics (the enzymatic networks, essential to the functionality of cell physiology). These biocatalysts are responsible for the formation of all the ATP in the cell, DNA synthesis, apoptosis, mitosis, cytoskeleton organization, DNA repair, alternative splicing, chromosome regulation, gene expression, post-translational modifications, Golgi apparatus activity, etc. In short, all the main activities in the cell are mediated by enzymes. They constitute the fundamental elements of all essential physiological processes, which allow a cell to grow, multiply, and adapt to the external medium. However, being biologically essential, many issues of enzyme activities remain poorly understood which warrant further investigation.

There are a lot of important biomolecules in the cell, such as ATP, nucleotides, lipids, structural proteins, etc., but enzymes constitute the unique active biomolecules with the capacity to reordering the atoms, i.e., the ability to transform some molecules into others creating or destroying bonds between atoms.

Genes encode information about the primary structure of each enzyme, and this amino acid sequence determines its catalytic specificity. However, the enzymatic functionality, the set of different patterns of activity carried out by every enzyme, cannot be predicted from the primary structure alone. Enzymes are not rigid molecular structures; they have complex dynamic conformations and are continuously undergoing a wide range of conformational fluctuations which modify their catalytic activities. Even in the native state, an enzyme exhibits a range of complex interconverting conformations driven by thermodynamic fluctuations (Ramanathan et al., 2014; Agarwal, 2019). Besides, the dynamic changes in the hydration shell rapidly control these conformational fluctuations of enzymes (they can be described by an energy landscape) that hinder the connection between enzyme sequence and its catalytic function (Fenimore et al., 2004; Dunaway-Mariano, 2008).

Besides, the activities of most enzymes are not constant in biological environments; instead, catalytic activity can be complexly regulated thus emerging sophisticated activity patterns (Figure 1). These enzymatic behaviors depend on different factors, mainly substrate concentration, and any inhibitors or activators. The substrate fluxes and regulatory molecular flows form complex dynamical networks and, as a result, all metabolite concentrations permanently change with time due to collective interactions between them. Therefore, metabolite concentrations are a function of self-organized dynamics that emerge in the enzymatic networks (see below). Such collective connectivity integrated into biochemical networks makes highly variable the rate at which catalytic reactions proceed (see Figure 2), being able to exhibit infinite patterns of activity (De la Fuente, 2014, 2015). It should be noted that a great number of fundamental elements of the enzymatic activity under cellular conditions, continue to be incomplete after decades of research (Ringe and Petsko, 2009).

FIGURE 1

FIGURE 2

Practically all physiological processes rely on the enzymatic activity which shapes the cellular metabolism. This is an essential issue to understand the basic principles of biochemical functionality at all levels of its cellular architecture (Noble et al., 2014). Knowing the cellular metabolism is necessary to understand the laws and principles that determine cellular functions as a whole.

Enzymatic Self-Organization at the Molecular Level

How are enzymatic reactions organized in cells?

Over decades it was widely assumed that enzymes were randomly distributed in the intracellular medium. However, numerous researches on protein-to-protein interactions have demonstrated that homologous or heterologous protein-protein interactions shape complex macromolecular associations (Pang et al., 2008). More precisely, yeast proteome analysis has established that up to 83% of all proteins constitute supra-molecular associations (Gavin et al., 2002), and most of the cellular processes are not carried out by isolated proteins, but rather by multiprotein complexes (Gavin and Superti-Furga, 2003; Kastritis and Gavin, 2018). Such associative protein organization occurs in both eukaryote and prokaryote cells (Ho et al., 2002; Kühner et al., 2009).

Multi Enzymatic Associations

Different approaches have been applied to detect and study protein sets, such as advanced experimental molecular biology procedures, bioinformatics techniques, statistical physics tools, and computational methods (SabziNezhad and Jalili, 2020; Sartori and Leibler, 2020).

Of all protein complexes, enzyme-enzyme complexes are the most essential for biological functionality. The basic concept of multienzyme associations was initially suggested by A. M. Kuzin in the seventies of the last century (Kuzin, 1970). Besides, Paul A. Srere was the first to observe a non-trivial supramolecular enzymatic organization in the citric acid cycle (Srere, 1972), who later called “metabolon” to all these types of multienzyme sets that can catalyze two or more sequences of a metabolic route (Srere, 1985, 1987).

The multienzyme complex is a functional reversible superstructure formed by several enzymes of a metabolic pathway that uses non-covalent interactions such as hydrogen bonding, cation-pi interaction, and hydrophobic force to maintain ensemble several catalytic elements (Figure 1). Such enzymatic complexes allow highly selective reactions, increasing their catalytic efficiency (Conrado et al., 2008; Araiza-Olivera et al., 2013).

Many enzymes belonging to metabolic pathways that shape these complex supramolecular structures have been studied, for instance, tryptophan synthase complex (Hilario et al., 2016), fatty acid synthase (Maier et al., 2006), glycolysis (Kohnhorst et al., 2017), Calvin cycle (Gontero et al., 2018), cellulosome (Dou et al., 2015), pyrimidine biosynthesis (Serre et al., 1998), purine biosynthesis (Hoskins et al., 2004), synthesis of 5′-phosphosulfate (Sun and Leyh, 2006), aldolase dehydrogenase complex (Carere et al., 2011), polyketide synthases (Tsai and Ames, 2009), proteasome (Jung and Grune, 2012), pyruvate dehydrogenase (Izard et al., 1999), carbamoyl phosphate synthetase (Thoden et al., 1997), alpha-ketoglutarate dehydrogenase (Tretter and Adam-Vizi, 2005), urea cycle (Watford, 1989), and aldehyde–alcohol dehydrogenase (Kim et al., 2020).

A multienzyme complex provides efficient advantages over individual enzymes. In many observations, multienzyme sets can transfer the metabolites from one enzyme to the active site of the next consecutive enzyme without diffusion into the bulk medium, avoiding their eventual loss (substrate channeling). The time of diffusion between successive catalytic steps is minimized. This non-covalent direct transfer process of metabolic intermediates originates faster catalysis, making the catalysis more efficient in the multienzyme sets (Ovádi and Srere, 2000; Ovádi and Saks, 2004; Castellana et al., 2014). Substrate channeling may occur in several ways i.e., across enzyme channels or over the electrostatic surface of the protein superstructure (Sweetlove and Fernie, 2018; Svedružić et al., 2020; Zhang and Fernie, 2020).

Also, intracellular membranes and structural proteins may reversibly interact with enzymatic sets originating microcompartments, which represent another advantage of multienzyme complexes (Saks et al., 2007, 2009; Monge et al., 2009). Thus, it has been observed in many cases that enzymes activity is produced in these small specialized microenvironments, limited either by surface or restricted by volume, which is more stable against the biochemical changes in the intracellular medium, allowing more efficient catalytic reactions (Araiza-Olivera et al., 2013; Küken et al., 2018).

Other pieces of evidence have shown that some metabolic compartments are liquid-like, which are formed by phase separation from the cytoplasm. Such liquid-like compartments describe non-membrane-bound restricted spaces as phase-separated, also creating small specialized microenvironments (Hyman et al., 2014).

In the global molecular crowding that constitutes the cell, enzymes rarely function in isolation, but rather as components of catalytic sets and interdependent coupled bioprocesses (Greco and Cristea, 2016).

Dissipative Self-Organization of Enzymatic Activities, the Main Source of Molecular Organization in the Cell

Cells are open dynamic metabolic reactors that exchange matter and energy with the external environment, obtaining fundamental elements for survival from the nutrients through multi-enzymatic reactions. All living organisms exhibit highly ordered enzymatic dynamic processes and use the energy of nutrients to maintain non-equilibrium states, where sophisticated structures, complex biochemical patterns, and collective biomolecular self-organization emerge. When any cell reaches equilibrium with the external medium, its characteristic high functional and molecular order disappears and it dies.

Thermodynamics is a necessary field of scientific knowledge to properly understand some essential principles of biochemical processes (von Stockar, 2010; Kondepudi and Prigogine, 2014). So, entropy is an important concept that allows quantification of uncertainty and disorganization in a system. In isolated systems, entropy never decreases. Conversely, it increases toward a maximum at equilibrium. In open systems, which exchange energy and matter with their environment, entropy can be sustained or decrease, avoiding a state of thermodynamic equilibrium. The negative variation (producing a consequent increase in the entropy of the surroundings) corresponds to a relative order inside the system. This is possible by the dissipation of energy i.e., the energy becomes not only unavailable but irrecoverably lost in the surroundings through irreversible processes (Niebel et al., 2019). The tendency of an open system that consumes energy with low entropy to dissipate energy with higher entropy to its surroundings (thus, the system succeeds in liberating itself of part of its own produced entropy) can allow self-organization of the system (Ebeling and Ulbricht, 1986). This principle is fundamental to understand the main source of molecular order in the cell.

The theoretical basis of self-organization was formulated in 1977 by the Nobel Prize Laureate in Chemistry Ilya Prigogine in his work on dissipative structures (Nicolis and Prigogine, 1977). Specifically, dissipative self-organization at the molecular level constitutes spontaneous highly ordered microscopic structures far from thermodynamic equilibrium which are characterized by coherent spatial and/or temporal patterns. Due to complex non-linear interactions among molecular components and driven by energy dissipation, these self-organized structures emerge increasing the structural and functional complexity of the system. Such emergent self-organized dynamic structures in space and time are called self-organized dissipative structures (Nicolis and Prigogine, 1977; Kondepudi and Prigogine, 2014).

Dissipative self-organized processes use the energy inflow to generate a negative entropy variation in the open system which corresponds to an emergent positive increment in the information self-contained in the system. Such information increases the complexity, being able to produce highly-ordered macrostructures and complex functional dynamic behaviors. Non-linear interactions and irreversible processes may amplify fluctuations leading to a dynamic state, far from the equilibrium, in which the system becomes spatially and temporally self-organized (Klimontovich, 1999; Halley and Winkler, 2008; Misteli, 2009)

One of the most significant scientific discoveries of the 20th century concerns this new source of molecular order in nature, the dissipative self-organization. Such a source is critical to understand the dynamic activities of multienzyme complexes and the systemic functional organization in the cell. The work of Ilya Prigogine represents a profound and original treatment of molecular organization in nature, with deep implications for understanding biological life at its most fundamental levels.

Dissipative Self-Organization of Multienzyme Complexes

How does the functionality of multienzyme complexes in the prevailing conditions inside the cell work?

It is well established that enzymes, in addition to form multicatalytic sets, can shape dissipative structures in which two fundamental types of dynamic self-organization can emerge: temporal rhythms and spatial waves.

When a multienzyme set operates far enough from equilibrium and dissipative self-organization emerges, all the enzymes of the complex do function as a whole, showing reactive coordination between them (long-range coherence) in such a way that all the substrate and product concentrations spontaneously start to oscillate along time (temporal rhythms), see Figure 1. As a consequence, thousands and thousands of molecules and ions that shape the multicatalytic system (substrates, products, protons, regulatory molecules, and other metabolites) experiment massive oscillatory reorganizations in their molecular concentrations. Such dynamics are mainly characterized by collective synchronized behaviors, functional correlations between molecular components separated by macroscopic distances, coherent patterns, and highly coordinated integrative processes (De la Fuente, 2010). Dissipative self-organized catalytic behaviors find their roots in the non-linear enzymatic processes which mainly involve allosteric regulation, autocatalysis, cooperativity, and feed-back interactions (Goldbeter, 2002, 2007, 2018; De la Fuente, 2014).

All metabolite concentrations of multienzyme sets that exhibit dissipative self-organization processes present transitions between complex oscillations and quasi-steady states under cellular conditions (De la Fuente, 2015). The quantification of certain intracellular molecules through nano biosensors in living cells has shown complex quasi-steady states and oscillatory patterns whose dynamics are never constant (Ozalp et al., 2010). Different experimental pieces of evidence suggest that oscillatory behaviors are much more frequent than quasi-steady states in cellular conditions (Lloyd and Murray, 2005, 2006).

Biomolecular oscillations are a genuine behavior of the complex self-organized enzymatic processes in basic units of life. Such oscillatory behavior spontaneously occurs far from thermodynamic equilibrium, when the system is formed by irreversible enzymatic processes with non-linear kinetics, Figure 1. Intracellular molecular oscillations were reported in: ATP, ADP, and AMP nucleotides (De la Fuente et al., 2014), microtubule polymerization (Lange et al., 1988) and other cytoskeletal structures (Kruse and Jülicher, 2005), cyclins (Hungerbuehler et al., 2007), cyclic AMP concentration (Holz et al., 2008), cytokinins (Hartig and Beck, 2005), free fatty acids synthesis (Getty-Kaushik et al., 2005), actin polymerization (Rengan and Omann, 1999; Negrete et al., 2016), biosynthesis of phospholipids (Márquez et al., 2004), urea cycle (Fuentes et al., 1994), proteolysis (Kindzelskii et al., 1998), glycolysis (Dano et al., 1999), metabolism of carbohydrates (Jules et al., 2005), intracellular glutathione concentration (Lloyd and Murray, 2005), Krebs cycle (Wittmann et al., 2005), mitochondrial metabolic processes (Aon et al., 2008), photosynthetic reactions (Smrcinová et al., 1998), CO₂ (Tyson, 2002), protein kinase activities (Chiam and Rajagopal, 2007), transcription factors (Garmendia-Torres et al., 2007), respiratory metabolism (Lloyd et al., 2002), intracellular calcium concentration (Ishii et al., 2006), peroxidase-oxidase reactions (Møller et al., 1998), membrane receptor activities (Placantonakis and Welsh, 2001), membrane lipid oscillation (Nakamura, 2018), ERK/MAPK metabolism (Shankaran et al., 2009), intracellular pH (Sánchez-Armáss et al., 2006), intracellular free amino acid pools (Hans et al., 2003), membrane potential (De Forest and Wheeler, 1999), metabolism of mRNA (Klevecz and Murray, 2001), min-proteins (Wettmann and Karsten, 2018), beta-oxidation of fatty acids (Getty et al., 2000), amino acid transports (Barril and Potter, 1968), insulin secretion (Tyson, 2002), etc.

Moreover, a metabolic orchestration has been observed in cells at a systemic level, by which the entire metabolome and most of the transcriptome oscillate (Klevecz et al., 2004; Lloyd and Murray, 2006; Murray et al., 2007). More specifically, several studies have analyzed oscillatory processes in the transcriptome (Tonozuka et al., 2001; Tian et al., 2005; Chabot et al., 2007). These studies have shown that at least 60% of all transcriptions oscillate with an approximate period of 300 min (Tu et al., 2008). Microarray analysis from S. cerevisiae in continuous synchronous culture shows a genome-wide oscillation in transcription processes coupled to respiration with maxima of transcript levels at intervals of about 40 min (Klevecz et al., 2004).

The oscillatory periods of metabolic rhythms range from milliseconds to minutes and hours (Chance et al., 1973; Berridge and Galione, 1988; Aon et al., 2006; Brodsky, 2006; Roussel et al., 2006). Complex periodic oscillations, including bursting rhythms and deterministic chaotic phenomena, have often been detected (Olsen and Degn, 1985; Dekhuijzen and Bagust, 1996; Almeira and Gusman, 2017; Heltberg et al., 2019).

Dissipative structures are very diverse (Goldbeter, 2018). For instance, besides rhythmic metabolic behaviors, the coexistence between two stable oscillations in a biochemical set (birhythmicity) and the coexistence of several stable steady states (multistability) may also emerge in the same molecular system. Multistability has been observed in different biosystems such as the calcium/calmodulin-dependent protein kinase II (Zhabotinsky, 2000), NF-κB signaling (Pêkalski et al., 2013), apoptosis (Bentele et al., 2004), cytoskeleton control (Byrne et al., 2016), lac operon (Santillán and Mackey, 2008), and signal transduction (Xiong and Ferrell, 2003). The coexistence between two oscillatory regimes and a steady-state with periodic behavior (hard excitation) has also been observed in numerical studies of multienzyme complexes (De la Fuente et al., 1998c; De la Fuente, 1999) and chemical systems (Guttman et al., 1980).

Another class of self-organized biological dissipative processes is circadian rhythms, which exhibit an oscillatory period close to 24 h (dark–light cycle during the Earth’s rotation period). These endogenous autonomous oscillators can adapt their internal metabolism to changes in the external environment (light, temperature, food availability, etc.) during 24-h day/night cycles (Yanling et al., 2019). Circadian rhythms exist in all types of cells from prokaryotes to eukaryotes, and these dissipative patterns regulate a great variety of important physiological processes (Wijnen and Young, 2006). For instance, oscillating genes are also usually circadian (Tuttle et al., 2005). Even, it has been observed in some cells that 80–90% of the transcriptome show a rhythmic gene expression with cycles of 24–26 h (Connor and Gracey, 2011).

Circadian transcription-translation cycles are intertwined with different metabolic activities and respiratory oscillations, which ensures the systematic coordination and integration of cellular physiological processes (Causton, 2020).

Lastly, a fundamental type of dissipative structure in cells is the spatial traveling waves which consist of three-dimensional self-organized coherent oscillation of metabolite concentrations that propagates progressively across the intracellular medium, reminiscent of a wave moving across water. These wave pulses of biochemical activity, moving through subcellular domains over large intracellular distances, are too fast (5–30 μm s^–1 for calcium waves; Jaffe, 2002) and represent an essential mechanism for long-range functional interconnection in the cell. Spatial traveling waves have a crucial role in coordination and synchronization among different metabolic processes, and these dissipative patterns exhibit other characteristics in their shape of oscillation, cellular location, and molecular composition (Guthrie et al., 1999; Scemes and Giaume, 2006; Carsten and Karsten, 2017). Some examples of spatial biochemical oscillations have been observed in calcium ions (Jaffe, 2008), actin dynamics during cell locomotion (Vicker, 2002; Allard and Mogilner, 2013), apoptotic signals (Cheng and Ferrell, 2018), mitochondrial redox (Romashko et al., 1998), NAD(P)H (Kindzelskii and Petty, 2002; Slaby and Lebiedz, 2009), sodium ions waves in astrocytes (Bernardinelli et al., 2004), phosphoprotein waves (Markevich et al., 2006), Cdk1 waves implicated in the cell cycle (Deneke et al., 2016), mitotic waves (Nolet et al., 2020), adenosine triphosphate (Ueda et al., 1990; Newman, 2001), NAD(P)H and protons (Petty et al., 2000), phosphatidylinositol (3,4,5)-trisphosphate (Asano et al., 2008), ROS molecules (Zhou et al., 2010), and mitochondria activity (Kurz et al., 2010).

Each cell is a super dynamic metabolic system in which self-construction and self-destruction of molecules occur following complex dissipative catalytic rhythms, coordinated by complex spatial traveling waves.

Notice that after more than five decades of research many aspects of self-organized processes remain poorly understood.

Molecular Self-Assembly

Apart from dissipative structures, another mechanism of a molecular organization is self-assembly, which can be defined as the spontaneous formation of supramolecular structures usually driven by specific non-covalent interactions, yielding functional and stable three-dimensional macromolecular complexes. The forces that directly generate molecular self-assembly tend to be weak; specifically, they are non-covalent interactions such as van der Waals electric forces, hydrophobic effects, electrostatic attraction, and hydrogen bonding. Under cellular conditions, there are two classes of molecular self-assembly: static and dynamic. Static self-assembly, an important spontaneous mechanism generating molecular order, involves processes that are near thermodynamic equilibrium and do not dissipate energy. Some examples are the formation of the lipid bilayer, viral capsid, protein aggregates to hold the quaternary structure, and some supramolecular polymerization (Sartori and Leibler, 2020). Dynamic self-assembly occurs in non-equilibrium states when the system is dissipating energy. Therefore it is governed by irreversible enzyme kinetics, for example, the self-assembly of actin filaments due to hydrolysis of ATP (Whitesides and Grzybowski, 2002; Rieß et al., 2020). Both classes of self-assembly processes may also occur simultaneously with dissipative self-organization (Halley and Winkler, 2008).

Dissipative self-assembly and dissipative self-organization are the fundamental pillars of the molecular order of all living cells (De la Fuente, 2015).

Dissipative Metabolic Networks and the Emergence of the SMS

Two main principles govern the enzymatic organization in the cell: metabolic segregation and systemic metabolic integration.

Metabolic Segregation

All cellular metabolisms are segregated in multiple autonomous biochemical specializations. This separation of enzymatic roles demands that correlated catalytic processes are grouped. Consequently, the systemic metabolism is functionally segregated at a great variety of specialization levels, originating different types of specific catalytic activities (e.g., the lysosomal metabolism, the oxidative phosphorylation, the signal transduction, etc.).

All modular metabolic networks that shape the cell can be considered as discrete dynamic sets, which perform relatively autonomous activities with specific and coherent enzymatic functionalities. Several scientific works have demonstrated that these elemental networks shape specialized modular units showing the genuine segregation of the metabolic activity (Ravasz et al., 2002; Nurse, 2008; Hao et al., 2012; Kaltenbach and Stelling, 2012; Geryk and Slanina, 2013). Segregation networks have been detected in all principal cellular metabolic activities, e.g., “enzymatic organization of metabolic pathways (Muto et al., 2013), chaperone activities (Korcsmáros et al., 2007), chemotaxis (Postma et al., 2004; Shimizu et al., 2010), apoptosis processes (Harrington et al., 2008), cell cycle (Boruc et al., 2010; Hsu et al., 2011), Golgi apparatus (Nakamura et al., 2012), and kinetochore organization (Petrovic et al., 2014). Specific modular metabolic networks also take part in the transcriptional system such as miRNA regulatory networks (Gennarino et al., 2012), transcription factors networks (Neph et al., 2012), RNA polymerase complexes (Schubert, 2014), and mRNA dynamics (Vlasova-St Louis and Bohjanen, 2011)” (De la Fuente, 2015).

Dissipative multienzyme complexes can also be considered fundamental modular networks (De la Fuente et al., 2008; De la Fuente, 2010, 2014). In this respect, the Metabolic Subsystem concept was suggested in 1999 to define dissipatively structured enzyme complexes in which autonomous catalytic processes with complex quasi-steady-state patterns and molecular rhythms spontaneously emerge (De la Fuente et al., 1999b). Multienzyme complexes carry out specific catalytic activities that also interact reversibly with structural proteins and membranes, originating microcompartments thermodynamically open, in which biomolecular oscillations, substrate channeling, and integrative mechanisms may emerge. These functional integrative processes increase the efficiency of the self-organized multienzyme complexes. Dissipative multienzyme complexes represent highly efficient modular networks capable of performing discrete catalytic activities that allow cellular functional organization. Such self-organized multienzyme sets constitute the essential and basic catalytic elements of the cell (De la Fuente, 2015).

Metabolic Integration

Aside from the metabolic segregation, catalytic activities present another essential property of organization: systemic metabolic integration.

Biocatalytic sets organized into distinct basic modular networks are functionally integrated with the cell through global coordination resulting from the concerted metabolic actions of the different specialized biochemical activities as a whole (Almaas et al., 2004, 2005; Yoon et al., 2007; Buescher et al., 2012; San Román et al., 2014).

While functional segregation expresses the partial interdependence of specialized metabolic nets (Figure 1), metabolic integration is an additional principle shared by the autonomous catalytic process that reflects high deviation activities from the functional autonomy of these modular networks, in such a way that the collective metabolic elements of the cell produce coherent systemic dynamics (Figures 2, 3) that are functionally integrated, shaping a super-complex systemic dynamic structure: the SMS (De la Fuente, 2015).

FIGURE 3

The SMS was observed for the first time in 1999 in an exhaustive numerical analysis with several millions of different DMN (De la Fuente et al., 1999b). Such systemic organization is defined by a small set of different self-organized multienzyme complexes which always present active states, while the rest of dissipative multienzyme subsystems exhibit on-off active states (Figures 2, 3). Under this integrative metabolic activity, each active multicatalytic subsystem generates output responses with complex dynamic transitions between steady states and oscillatory patterns. The set of multienzyme complexes always within an active state are called the metabolic core (De la Fuente et al., 2008).

Later, this global enzymatic system was verified employing flux balance analysis in several prokaryotes and eukaryote cells such as S. cerevisiae, H. pylori, and E. coli. Besides, essential enzymatic reactions for biomass formation were observed in the metabolic core (Almaas et al., 2004, 2005).

Adenylate energy system constitutes another global process responsible for the functional integration of cellular metabolism. Energy is an essential element to keep the metabolic architecture of cells, and adenosine nucleotides couple all bio-energetic enzymatic processes one each other. The dissipative multienzyme complexes regulate their functional activity through changes in ATP, ADP, and AMP levels, along with other biochemical factors.

The adenosine nucleotide levels are determined by the adenylate energy system at a systemic level. Such dynamic energetic structure is shaped by enzymatic reactions which interconvert AMP, ADP, and ATP, as well as energy-consumption processes coupled to the synthesis of ATP (De la Fuente et al., 2014).

Five decades ago, a primary systemic ratio between ATP, ADP, and AMP concentrations was proposed by Atkinson (Atkinson and Walton, 1967) to calculate the energetic cellular level, which was denominated the AEC. Cellular energy quantification shows that almost all organisms appear to keep their AEC within narrow values under growth conditions, more specifically between 0.7 and 0.95, despite the extreme complex fluctuations in the adenine nucleotide concentrations in the cell. This quantitative physiological invariant represents a key property of the integrative mechanisms of energetic and functional dynamics operating at the systemic cellular level that is supported by numerous experimental researches (De la Fuente et al., 2014).

In short, two essential principles govern the different ways of enzymatic organization in the cell: metabolic segregation and systemic metabolic integration. Both are not excluding, but complementary, and the interplay between them represents a fundamental property for the systemic cellular metabolism and its functional global architecture.

On the other hand, in modular and systemic networks, three main issues of metabolic connectivity can also be considered: structural connectivity (structural molecular links), functional connectivity (non-causal statistical dependencies), and effective connectivity (causal dependencies), see an example in Figure 3 (De la Fuente et al., 2011, 2013; De la Fuente and Cortés, 2012).

Enzymatic Information Processing: The Second Primary Source of Molecular Order in the Cell

A significant number of experimental and numerical works have shown that basic units of life exhibit complex self-regulatory behaviors originated by biomolecular information processing, the second essential source of order and organization in cells.

Self-Regulation in Dissipative Multienzyme Complexes

According to the Information Theory, TE permits quantifying, in bits, the information of biomolecular patterns (see Supplementary Material 02 for more details) (Schreiber, 2000). The information processing in yeast glycolysis, a dissipative multienzyme complex, was studied in 2012 using TE (De la Fuente and Cortés, 2012). This biochemical pathway is one of the most analyzed and it has been considered a classical prototype of a metabolic oscillator (see Supplementary Material 02). In Figure 1A, the main irreversible enzymatic processes and the basic structural connectivity of the yeast glycolytic pathway are shown. Such metabolic subsystem under a specific glucose source allows observing a quasi-periodicity route to chaos (Ruelle-Takens-Newhouse route). Figure 1B displays a transition sequence of this route with a stable periodic pattern, quasi-periodic rhythms, complex quasi-periodic oscillations, and deterministic chaos (De la Fuente et al., 1996a).

The molecular information from these glycolytic oscillatory processes was captured by TE. Specifically, the information values in bits were normalized with the maximum TE value obtained, so the value of 1.00 represents the TE value between the functional connection of phosphofructokinase and pyruvate kinase enzymes which is the leading source of information in the glycolytic pathway (De la Fuente and Cortés, 2012). The functional influence obtained ranged from 0.58 ≤ TE ≤ 1.00 (with mean ± SD = 0.79 ± 0.12), which suggests highly effective connectivity in the dissipative multienzyme subsystem. The analysis also indicated that the flows of functional connectivity (contained in the different metabolite patterns) change in all catalytic transitions studied (see Figure 1B), under other inputs of glucose concentrations. The leading source of molecular information coincided with the phosphofructokinase enzyme at the edge of deterministic chaos, when complex quasiperiodic patterns in the glycolytic pathway appeared. This is consistent with other investigations: when a dynamic system works on the border between deterministic chaos and periodic order, its complexity reaches its maximum (Kauffman and Johnsen, 1991; Bertschinger and Natschlager, 2004).

The difference between the total TE input minus the TE output in each irreversible enzyme was calculated to study the total information flow of the metabolic subsystem. From these differences, negative values indicated targets of causality flow, while positive numbers were interpreted as sources. The highest value of total transfer information corresponded to the phosphofructokinase enzyme (0.41), and it occurred with the emergence of complex quasiperiodic patterns. The phosphofructokinase (E2) was the principal source of molecular information in the system, and the pyruvate kinase worked as a sink, which corresponds to a functional target; the irreversible enzyme hexokinase appeared with a flow close to zero and less constrained (Figure 1).

In summary, as a whole, the glycolytic pathway operates as a molecular information processing system, which at every moment redefines metabolic information flows that make each irreversible enzyme evolve with well-defined catalytic activities, under different inputs of glucose concentrations. The glycolytic metabolic subsystem behaves as a catalytic unit highly self-organized (dissipative processes) and self-regulated against the external changes of glucose input concentrations (molecular information processing). Each irreversible enzyme carries out three informative operations simultaneously: molecular information reception, molecular information integration, and source of new informative molecular flows. Thus, the dissipative multienzyme subsystems exhibit, at any time, precise sets of molecular instruction fluxes, and, as a result, each multienzyme complex shows specific and well-defined metabolic activities that permanently self-regulate the system. These dissipative multienzyme complexes can be considered the most fundamental units to accomplish the processing of molecular information in the cell (De la Fuente and Cortés, 2012).

Self-Regulation by Molecular Information Processing in DMN

To understand how different multienzyme complexes functionally operate together, TE was applied for the study of the catalytic activities of DMNs (De la Fuente et al., 2011). In this work, a complex DMN of 18 metabolic subsystems was analyzed which is depicted in Figure 2. Each catalytic unit represents a self-organized multienzyme set (MSb). Figures 3A,B shows the structural organization of substrate fluxes, two classes of allosteric processes (activatory and inhibitory signals), and covalently modifying activities. DMNs are open systems, and S1 and S2 represent substrate external sources to MSb3 and MSb10, and MSb6 activity ends outside the metabolic network.

First, the dynamics of that network were studied with two inputs of substrate sources, S1 and S2. The spontaneous emergence of a SMS was observed under these external conditions, i.e., MSb12 was permanently active working as the metabolic core, the other subsystems presented on-off shifting states, and MSb15 was inactive (Figures 2A,B). All active self-organized enzymatic subsystems presented complex catalytic behaviors with different steady and oscillatory patterns. For instance, in Figures 2C,D, some dissipative catalytic activities of the MSb12 are showed. Specifically, this dissipative multienzyme set exhibited activity cycles with more than 350 transitions between different spontaneous periodic oscillations and steady-states.

Removing S2 and keeping only S1, as the external source of the substrate, a different functional reorganization emerged in the network (Figures 3D–F, 4, 5). However, the SMS was maintained (the same metabolic core MSb12, while the rest of the subsystem, including the MSb15, was eventually active with on-off dynamics). Even more, the network adapted its metabolic activity to the new environmental situation (a single stimulus) through its structural and flux plasticity (Almaas et al., 2005: De la Fuente, 2015).

FIGURE 4

FIGURE 5

The TE of the biochemical network was evaluated by analyzing the amplitude of metabolic patterns (Figures 2E,F). There were just nine dissipative multienzyme sets with significant TE values when two concurrent stimuli were studied (Figure 3E). In this picture, the thickness of the lines was proportionate to the TE values, whilst the directions of arrows indicated the directionality of the effective connectivity. The highest values of molecular information corresponding to the connection between MSb13 and MSb16 were also observed under two stimuli (Figure 3E). The analysis of TE, when only a source of external substrate S1 was maintained is depicted in Figure 3F. These results indicated that the molecular information structure was far more complex: 1° effective connections were found in 10 multienzyme subsystems; 2° the density of informative connections was incremented in the metabolic network (from 10.802 to 13.580); 3° the strength of links was also incremented (from 0.079 to 0.133 in average). 4° the maximum TE value was 0.377 (one stimulus) against 0.179 bits (two stimuli). Quantitative and qualitative alterations of TE values were evident in the metabolic core (see some examples in Figure 4). In particular, it can be observed that the core gets a molecular information signal from MSb8 under two stimuli (with a TE of 0.116, Figure 4A, on the left). However, the TE flow direction is reversed, and the flux of molecular information is sent to MSb8 by the core, under one stimulus (with a TE of 0.153, see Figure 4A, on the right).

According to the TE analysis, it can be observed that, in addition to the molecular network’s structure (Figures 3A,B), a systemic functional structure formed by molecular information fluxes emerges in the network (Figures 3E,F, 5A,B). This SMS defines a complex effective information system that integrates every catalytic activity in a global metabolic entity. Such systemic informative structure is formed by modules, and the dynamic transitions between them shape metabolic switches that enable essential functional changes in the network (Figure 5). Therefore, the functional switches and the modules of effective connectivity are also critical elements of the SMS.

In summary, a complex structure of information processing emerges in the dissipative metabolic system at the global level, similar to parallel computing. “The dissipative network behaves as a highly complex decentralized information processing super-system that generates molecular information flows between the self-organized enzymatic sets, forcing them to be functionally interlocked, i.e., each catalytic set is conditioned to cooperate with others at a specific and precise activity regime in concordance with the activity. According to the overall functional process, the network defines sets of biochemical instructions at each moment that make every metabolic subsystem evolve with precise and particular dynamic catalytic patterns” (De la Fuente, 2015). In short, the network is capable of self-organizing far from equilibrium through dissipative processes (exhibiting complex quasi-stationary and oscillatory catalytic patterns) and self-regulating versus external perturbations through molecular information processing, adapting dissipative metabolite patterns to alterations coming from outside the system. Besides, the information generated in these processes highly increases the complexity of cellular system (De la Fuente et al., 2010, 2011).

One of the best-known examples of information processing at the cellular level is the E. coli chemotaxis network. This is one of the most studied metabolic complexes able to store biochemical patterns by processing molecular information (see Supplementary Material 02 for more details). On the other hand, information processing at a molecular level has been observed in numerous biochemical systems (see Supplementary Material 02).

Despite all investigations, many aspects of molecular information processing remain to be elucidated.

Hopfield-Type Systemic Attractors in DMN. Cellular Systemic Behaviors

“Memory is a fundamental requirement for efficient information processing” (De la Fuente, 2015). In this regard, Hopfield-like dynamics with memory were verified in DMN in 2013 (De la Fuente et al., 2013). In such a study, global information and energy of the emerging biochemical dynamics in DMN (Figure 2) were studied by using methods of Artificial Intelligence and Statistic Mechanics. In addition to the macroscopic properties analyzed, Hopfield-like behaviors were obtained using a Boltzmann machine (Figures 6, 7). As it is well known, these Hopfield-like dynamics are characterized by manifesting associative memory (Hopfield, 1982).

FIGURE 6

FIGURE 7

Specifically, the preliminary analyses of the DMN showed spontaneous self-organization of all enzymatic patterns which lead to the emergence of a SMS characterized by a metabolic core formed by always active enzymatic processes, while the others are eventually active. Also, the dissipative metabolic network was able to self-adjust the internal enzymatic activities according to the external changes by self-regulatory processes (molecular information processing), exhibiting both metabolic structural plasticity (persistent alterations in the state of multi enzymatic subsystems, switches between active, inactive, and varying on-off states) and metabolic flux plasticity (changes in metabolic synapse values leading to a differential catalytic activity); this type of global self-regulation has been observed in several unicellular organisms (Almaas et al., 2004; Almaas, 2007). The quantitative study of the DMN showed the following main results:

First, the Lyapunov function of the biochemical system was calculated. Starting in any initial state under a specific external stimulus, the enzymatic dynamics of the system evolve intending to decrease the energy function. In noise absence, this function will monotonously decrease until it reaches a terminal metabolic state (local minimum corresponding to a stable state of the DMN). There are multiple functional stable states for a DMN.

Second, due to self-organized and self-regulated processes in the DMN, systemic attractors (locally stable states) emerge and globally govern the network. All enzymatic dynamics are stabilized in these systemic attractors. Multiple attractors (and therefore multiple local minima) define the complex landscape on which the systemic enzymatic dynamics rely. As a consequence of the attractor dynamics, the metabolic system works as a fully integrated and individual entity.

Third, “the systemic attractors that emerge in the system correspond to Hopfield-like dynamics which can store metabolic information patterns (Hopfield, 1982). The biochemical information contained in the attractors behaves as a functional metabolic memory, i.e., enzymatic activity patterns stored as stable states, regulates both the permanent catalytic changes in the systemic network and the metabolic responses originated to integrate the perturbations coming from the external environment properly. When an external stimulus pattern is presented to the system, the network states are driven by the intrinsic enzymatic dynamics toward a determined systemic attractor which corresponds to a set of memorized biochemical patterns” (De la Fuente, 2015). As it is generally known, such “pattern-completion” dynamics have long been studied in Hopfield networks, and they are usually referred to as “associative memory” (Hopfield, 1982; Hertz et al., 1991; Amit, 1992; Peretto, 1992).

In sort, Hopfield-like metabolic attractors store systemic enzymatic dynamics which may be properly recovered from external stimuli. Metabolic dynamics are systemically governed by these metabolic attractors which modulate the enzymatic activities, change the connections among dissipative multienzyme associations, and stably store these dynamic connections.

The study with Hopfield networks (De la Fuente et al., 2013) quantitatively evidenced, for the first time, the possibility that a cell can have associative memory. Such Hopfield-like memory would be a cellular type of epigenetic memory (De la Fuente, 2015).

Verification of Associative Memory in Unicellular Organisms

The principles of associative memory were discovered at the beginning of the 20th century by Nobel Laureate Ivan Pavlov in his classic studies with dogs, which opened up a new paradigm in behavioral sciences.

In all organisms with a nervous system, associative learning and memory are essential cognitive characteristics that allow obtaining critical information for adaptation and survival, mastering new behaviors through the association of diverse stimuli. These characteristics are ubiquitous in numerous species from mollusks to humans, but until now they have never been observed in individual cells.

To experimentally verify the emergence of associative memory in cells, the systemic motility patterns of two eukaryotic microorganisms such as A. proteus and M. leningradensis were analyzed under an associative conditioning setting (De la Fuente et al., 2019a). For this purpose, two stimuli were tested following Pavlov’s concept of experiments, a suitable direct-current (DC) electric field, and an N-formylated tripeptide linked to amoeba nutrition.

The most important aspects of this work on cellular associative conditioning are summarized below:

• -

  First, it was found that in the absence of stimuli cellular locomotion directionality displayed a random distribution when amoebae and metamoebae probed almost all the directions of the experimental chamber (Figure 8A).

FIGURE 8

• -

  In the next experimental step, galvanotaxis, amoebae, and metamoebae exhibited an unambiguous systemic reaction consisting of movement to the negative pole when subjected to a strong DC electric field of about 300–600 mV/mm (Figure 8B).

• -

  In the third step, the response of both unicellular organisms was studied during their exposure to a gradient of nFMLP, distributed in the left part of the experimental chamber (chemotaxis, Figure 8C). In these circumstances, most cells showed a stochastic migration with robust directionality toward the peptide.

• -

  In the fourth step, an induction process was performed. To this end, both stimuli, galvanotactic and chemotactic, were applied to amoebae simultaneously for 30 min. The nFMLP peptide solution was added to the anode well, on the left part of the chamber (Figure 8D). In these experiments, only about half of the amoebae and metamoebae migrated to the cathode and the rest moved to the anode site, where the peptide was deposited.

• -

  Finally, to check if the cells conserve their movement pattern to the anode-nFMLP, the cells that showed such migration in the fourth step were exposed once more for 30 min to a single electric field stimulus, omitting peptide at the anode. The individual cell trajectories analyzed after the experiments revealed the migration of most of these cells to the anode in the absence of peptide (Figure 9), so it was shown that after the induction process, amoebae and metamoebae developed a new pattern of locomotion: moving to the anode under galvanotaxis conditions. It is worth remembering that almost all cells in the DC field move to the cathode (see Figure 8B).

FIGURE 9

Under the application of two simultaneous stimuli (induction process), the DC electric field and the peptide related to the nutrition of amoeba, placed at the positive pole, some amoebae seem to relate these stimuli (anode and peptide), and subsequent control tests showed that most of the conditioned A. proteus and M. leningradensis changed their systemic behavior, moving toward the anode where there was no nFMLP peptide.

They acquired a new persistent pattern of cellular motion, described by the move toward the anode (Figure 9), instead of their established propensity to run to the cathode. The conditioned amoebae were able to associate simultaneous but unrelated past events, producing persistent movement (contrary to their natural tendency), which on average could continue for 44 min. Besides, a vast number of strict controls proved the robustness and predictability of this behavior.

These experiments, in which a single cell links two different simultaneous stimuli thus developing a new behavior, are the first published evidence that an associative Pavlovian-like memory is present in single cells (De la Fuente et al., 2019a).

Cellular conditioning may be necessary for critical functions, for example, the regulation of migration. Locomotion movements are fundamental not only to survive and avoid predators, but also in individual cells of multicellular organisms to carry out essential life processes like embryogenesis, organogenesis, immune response, tissue repair, and others. Several severe health conditions, such as cancer, vascular and heart diseases, and intellectual disability, may develop as the consequences of inaccurate control of this systemic motility property. Such migratory cellular abilities are also present in the metastatic process of cancer (De la Fuente and López, 2020).

Integrative Summary

In this review, we have mainly approached four aspects that characterize cellular functioning: the role of dynamic molecular processes, the mechanisms implicated in self-organization (dissipative structures) and the self-regulation of metabolic dynamics (molecular information processing), and the emergence of the cellular systemic properties.

More specifically, the work focuses on the main principles that govern enzymatic activity, under complex dynamic conditions prevailing inside the cell, which are crucial to elucidate the structural and functional architecture of basic units of life.

The scope of Systems Biology applied here is based on integrating concepts and knowledge in Biochemistry, Enzymology, Molecular and Cell Biology, and other biological sciences. We also discussed the results of quantitative research, analyzed by application of physical-mathematical tools such as statistical mechanics, differential calculus, discrete mathematics, artificial intelligence, computing, and others. Nevertheless, favoring the didactic purpose of the review we did not use the mathematical formulations.

The following paragraphs add some conclusive examples and summarize the main aspects of this review in an integrative way.

Cellular Molecular Dynamics

Strictly speaking from a biochemical viewpoint, the cell is a complex molecular reactor, extremely self-organized in a sophisticated manner, and in a permanent biomolecular recycling status (section “The Dynamics Originated by the Molecular Turnover and the Role of the Enzymes in These Super Complex Dynamics”).

Millions of biochemical reactions happen simultaneously in every basic unit of life at any time. This process of reactive transformations also includes subcellular structures. Nothing is inert, all the molecules and structures that make up each cell undergo complex chemical transformations (see Supplementary Material 01). Notice that, despite the efforts, many aspects of cellular recycling remain poorly understood by molecular experimentalists.

The cell is an extremely complex and dynamic system; the endless molecular turnover, which characterizes it, modifies the chemical composition of the whole reactor every second and, as a consequence, the concentration patterns of ions and molecules change permanently over time. These continuous cycles of synthesis and destruction, and incessant chemical molecular transformations shape the essential scenario in which cellular life is possible. When this dynamic process of molecular turnover collapses, the cell dies. This is the main systemic characteristic of all living unicellular organisms.

The consequences of self-regulated turnover dynamics are adequate cellular growth, the development of all physiological processes to maintain its functional structures, continuous adaptation to the environment, and mitosis.

The highly complex organization of cellular turnover dynamics defies the human intellect. Any attempt to synthetically reproduce this global molecular turnover, either in vitro or in silico, has failed so far. In every cell, molecular synthesis and destruction are compensated and harmonized between them, following complex and unrepeatable reactive patterns, whose laws and defining principles are still unknown.

Enzymes, Essential Cellular Molecules

Enzymes are fundamental molecules for metabolic life. Molecular reactive transformations are essentially chemical reactions, modifications on how atoms are bonded. Practically, all reactive transformations in the cell are mediated by enzymes, extraordinary macromolecular nanomachines responsible for breaking and joining covalent bonds, the fundamental chemical links in biological molecules (section “The Dynamics Originated by the Molecular Turnover and the Role of the Enzymes in These Super Complex Dynamics”). Cellular life would be impossible without enzymes; they are key protagonists in all reactive and physiological processes in the cell.

Unlike other cellular molecules, enzymes are active macromolecules responsible for the chemical work. Through their ability to decrease the activation energy of bonds, enzymes permanently modify how atoms are bonded, making possible the accelerated creation and destruction of molecules.

Catalytic enzymatic activities are termed metabolic processes. There are two basic types of enzymatic reactions, catabolism (implicated in molecule destruction) and anabolism (implicated in molecule synthesis). Together, these enzymatic activities are called cellular metabolism. In short, the cell is a complex enzyme-mediated metabolic system where reactive molecular transformations of synthesis and destruction happen unceasingly mediated by catalytic activities.

Dissipative Self-Organization, the Primary Source of Molecular Organization in the Cell

Reactive cellular processes have little to do with the Chemistry of Equilibrium. The fundamental characteristic of cellular biochemistry is that the catalytic reactions, considered globally, occurs far from the Thermodynamic Equilibrium. Under these particular conditions, irreversible enzymatic processes raise the spontaneous emergence of self-organized dissipative structures (section “Dissipative Self-Organization of Enzymatic Activities, the Main Source of Molecular Organization in the Cell”).

Different catalytic mechanisms allow processes far from Thermodynamic Equilibrium, and one of the most important is the non-linear enzyme kinetics of many irreversible enzymes with allosteric regulation (Goldbeter, 2002, 2007).

Two principal types of self-organized dissipative structures emerge in the cell, temporal molecular rhythms, and spatial waves. In temporal rhythms, molecular concentrations fluctuate over time in an oscillatory manner (Figure 1). Every cell appears as a dynamic system in which cellular self-destruction and self-construction happen following complex quasi-stationary patterns and unrepeatable molecular rhythms. On the other hand, the dissipative spatial molecular waves modulate the cellular metabolic activities, synchronizing and regulating functionally different enzymatic processes. All these self-organized behaviors are unlikely for the Chemistry of the Thermodynamic Equilibrium (De la Fuente, 2010).

So, the cell is a metabolic system in which self-construction and self-destruction occur following complex catalytic rhythms, mainly coordinated by spatial molecular waves.

The self-organized dissipative structures were discovered by Nobel Laureate Ilya Prigogine (Nicolis and Prigogine, 1977), and these dynamic processes constitute the primary source of molecular order in the cell. The concept of dissipative order is a fundamental issue to describe the molecular-functional architecture of living cells. Dissipative self-assembling and dissipative self-organization are the two pillars of the molecular order and the emergent functional complexity of cells (De la Fuente, 2015). However, being biologically essential, many issues of cellular dissipative structures remain unclear, which warrant further investigation.

Cellular Enzyme Organization. Cellular Metabolic Structure

At the cellular level, enzymatic activity dramatically depends on the collective functional structure of catalytic processes which is defined by the two fundamental principles that connect the main forms of enzyme organization in the cell: metabolic segregation and systemic metabolic integration, operating simultaneously (section “Dissipative Metabolic Networks and the Emergence of the Systemic Metabolic Structure”).

Enzymes do not work independently of each other in the cellular molecular crowding, but instead shape different types of multienzyme associations, grouped by their metabolic tasks and catalytic roles. Thus all chemical reactions that occur in living cells are segregated functionally (De la Fuente, 2015).

Moreover, these catalytic complexes interact reversibly with structural proteins and membranes originating thermodynamically open microcompartments, in which molecular rhythms can regulate the efficiency of the enzymatic reactions involved (De la Fuente, 2010).

The term Metabolic Subsystem was suggested in 1999 to refer to these self-organized multienzyme sets in which complex quasi-steady-state behavior and molecular oscillations may emerge spontaneously inside the cell (De la Fuente et al., 1999b). These dissipative catalytic complexes (Figure 1) represent highly efficient nanostructures capable to perform autonomous biochemical works which constitute basic catalytic units of the cellular organisms (De la Fuente, 2015).

Multienzyme complexes are organized into distinct modular networks with specific and coherent autonomous activities. At a superior level, the metabolic networks (Figure 2) appear to be integrated (systemic metabolic integration) forming a very complex dynamic super-system, the CMS. Such global structure was observed for the first time in 1999 by analyzing DMN through numerical studies (De la Fuente et al., 1999b, 2008) and was later corroborated by flux balance analysis in prokaryotic and eukaryotic cells (Almaas et al., 2004, 2005; Almaas, 2007; De la Fuente, 2015).

Information Processing. The Second Source of Molecular Order in the Cell

A large number of experimental studies have shown that unicellular organisms possess complex self-regulatory behavior originated by biomolecular information processing, another fundamental source of organization in the cell (section “Enzymatic Information Processing, the Second Leading Source of Molecular Order in the Cell”).

The dissipative multienzyme complexes permanently send information through metabolite patterns between their different enzyme components. As a result of the overall process, these enzymatic assemblies operate as basic molecular information processing units (Figure 1). Thus, each dissipative multienzyme complex defines, at any time, precise sets of molecular instruction fluxes, and as a result, each multienzyme complex exhibits well-defined and specific metabolic activities (De la Fuente and Cortés, 2012).

On the other hand, a complex parallel super-system of information processing appears in the metabolic networks when the systemic activity is considered (Figure 2). CMS behaves as a self-regulated dynamic entity that permanently sends specific regulatory signals to multi-enzyme complexes make them evolve in well-defined and precise metabolic activities by information processing, see Figures 3–5 (De la Fuente et al., 2010, 2011, 2013; De la Fuente, 2015).

Molecular information processing allows the formation in the cellular metabolism of a network of functional links that readjust the catalytic patterns, adapting them to the physiological needs of cells during environmental changes (self-regulation process). Besides, the information generated in these processes highly increases the complexity of cellular system. Physarum polycephalum is a paradigmatic example of the emergence of complex cellular behaviors mediated by information processing in unicellular organisms. This microorganism can discover the minimum-length option between two distant points in a labyrinth (Nakagaki et al., 2000; Nakagaki, 2001). Note that finding the shortest path problem in a maze needs a rigorous mathematical analysis (Miyaji and Ohnishi, 2008). P. polycephalum is also capable of improving the selection of the better route configuration obtained by the shortest Steiner’s minimum tree connections, thus developing adapted strategies to maximize its access to nutrients (Nakagaki et al., 2004a; 2004b). Even more, P. polycephalum achieves to solve difficult problems, for example, finding a high-quality solution to the problem of the traveling salesman, a question which is known to be NP-hard (Aono et al., 2011a, b; Zhu et al., 2011). This multinucleated amoeba is capable of designing an optimal network very similar to the purpose-intended system in the Tokyo railway organization (Tero et al., 2010). During the process of adaptation to different inputs, P. polycephalum succeeds memorizing changes in its environment, recalling them later to adapt its behavior to the new conditions appropriately, for example, anticipating a cold-dry pattern in the environment 1 h before the change happens (Saigusa et al., 2008). It has also been observed that this protist accomplishes complex dilemmas of multi-objective foraging (Bonner, 2010; Dussutour et al., 2010; Latty and Beekman, 2011). Recently it has been shown that P. polycephalum is capable of developing a kind of rudimentary learning (Boisseau et al., 2016). Molecular information processing has been observed in numerous other biochemical systems (see Supplementary Material 02).

Hopfield-Like Dynamics and Associative Memory

Efficient information processing is impossible without memory (section “Hopfield-Type Systemic Attractors in Dissipative Metabolic Networks. Cellular Systemic Behaviors”). The presence of Hopfield-like dynamics characterized by associative memory in dissipative enzymatic networks was numerically verified in 2013 (De la Fuente et al., 2013). Such a memory would be a manifestation of emergent properties underlying the complex dynamics of the systemic cellular metabolic networks when dissipative enzymatic self-organization and molecular information processing act together. This study was the first quantitative evidence that an individual cell can use associative memory, which is supported by epigenetic mechanisms (De la Fuente et al., 2013; De la Fuente, 2015). As a consequence of the emergent Hopfield-like attractors, the systemic metabolism works as a fully integrated and individual entity (Figures 6, 7).

Associative memory enables all organisms with a developed nervous system, from mollusks to humans, to learn and adapt successfully to specific environmental stimuli (Pavlov, 1927; Mackintosh, 1983), but until now this type of learning has never been described in unicellular organisms. Very recently, patterns of locomotion consistent with an associative conditioned behavior have been observed in two microorganisms, A. proteus and M. leningradensis (De la Fuente et al., 2019a). Both amoebae species were capable of linking two simultaneous but unrelated past events, demonstrating a new migration behavior that prevailed for a relatively long time, 44 minutes on average (Figures 8, 9). This finding confirmed in 2019 the prediction from the previous studies of 2013 on the emergence of Hopfield-like attractors with associative memory at the systemic cellular level (De la Fuente et al., 2013).

Along with the investigation of associative conditioning in single cells, a quantitative study on the dynamic characteristics of the systemic motility of A. proteus was carried out recently through advanced analysis of locomotion movements in intact and enucleated (cytoplasts) amoebae. This research also allowed the quantification of the importance of the nucleus during the cellular migration (De la Fuente et al., 2019b). Such study calculated specifically the movement fluctuations along the locomotion trajectories of the amoebae utilizing the “rmsf” method, a known technique in Statistical Mechanics based on a work by Einstein (1905), later elaborated and widely used to analyze time-series. The research revealed that individual cells without nuclei (cytoplasts) and intact cells displayed a similar migratory structure mainly described by non-trivial long-range correlations. “This dynamic memory (non-trivial correlations) represents a key characteristic of A. proteus movements during cell migration. It is worth noting that this temporal persistence (correlations with 41.5 min on average) in cells and cytoplasts matches with the Pavlovian-like dynamic memory (44 min on average)” (De la Fuente et al., 2019a).

The nucleus has been classically considered an essential element governing cell movement, nevertheless direct quantitative evidence was lacking. The migratory analysis of enucleated and non-enucleated amoebae described above (Bringas et al., 2017; De la Fuente et al., 2019b; De la Fuente and López, 2020) also demonstrated that the absence of nucleus does not impact the systemic movements of amoebae significantly. An independent group has reached a similar conclusion using other experimental methods in human cells (Graham et al., 2018).

To summarize, between the extracellular space and the DNA a dissipative self-organized CMS exists (De la Fuente, 2015). This super-complex dynamic structure behaves as a decentralized information processing system, “generating sets of biochemical instructions that drive each enzymatic activity to a particular and precise dynamic of change, allowing self-regulation and adaptation to the external medium. CMS permanently sends a flow of molecular signals to the DNA-associated metabolism, thus contributing to forming the complex transcriptional system. These molecular flows allow the accurate regulation of gene expression, so that only the specific polypeptides necessary for the adaptive maintenance of the CMS are synthesized. Altogether, both informative systems (CMS and DNA) coordinate the physiological development of the cell” (De la Fuente, 2015).

An extensive study performed with DMN (around 15,210,000), to research the mechanism behind the emergence of the CMS, allowed us to observe that this global superstructure is a property common to all SMSs with an high number of dissipative self-organized multienzyme complexes (De la Fuente et al., 2009). The fundamental factor that ensures the spontaneous emergence of CMS is a high multiplicity of dissipative metabolic sets inside the cell.

The high redundancy of the enzymes (with an increased number of copies) and the consequent multiplicity of dissipative self-organized multienzyme complexes that prevail inside the cell would determine the spontaneous formation of the SMS.

A continuous process of molecular synthesis and destruction represents thousands of metabolic reactions that happen permanently and simultaneously in the cell, ensuring the emergence, robustness, and stability of the CMS.

The only possible scenario for cellular life is a dynamic system in permanent self-construction and self-destruction process that guarantee the functionality of a high number of dissipative self-organized multienzyme complexes. There is no other alternative. This sophisticated and massive dynamics of endless molecular transformations encompasses the whole metabolic system and constitutes the most critical dynamic property, essential and definitory, of cellular life.

The word “metabolism” comes from the Greek metabole meaning “change” or “transformation.” This word reflects exactly the essential characteristic of the basic units of life. Rather than a molecular genotheque (box of genes) in evolution, the cell could be considered a singular dynamic metabolic unit capable to self-organize and self-regulate in an endless process of molecular turnover.

“Metabolism can be considered the largest known source of molecular complexity in nature” (De la Fuente, 2015). Roughly 3,700 million years ago (Dodd et al., 2017), an exceptional and singular metabolic dynamic emerged from primeval matter, characterized by a high structural and functional order, improbable for the Chemistry of Equilibrium. This extraordinary molecular organization perpetuates itself by direct transmission after every mitosis, and there is no scientific proof that another parallel metabolic organization has emerged “de novo.”

Across the millennia, many adaptive molecular mechanisms have been acquired in the biological evolution by the CMS. “Through the structural and functional complexities, the metabolic evolution has developed a large diversity of biochemical organizational forms, ranging from sophisticated bacterial metabolism, via the complex enzymatic networks, characteristic of protists, to the extraordinary structures and processes derived from the multicellular eukaryotic organization, such as embryogenesis or the neural networks present in higher mammals. The millions of living metabolic species represent the inexhaustible source of complexity developed by the dynamic forces of the cellular metabolism” (De la Fuente, 2015).

Every unicellular organism results from the direct and uninterrupted transmission of this extraordinary and singular metabolic organization. Over millions of years, sophisticated metabolic dynamics have evolved to originate complex biological supracellular structures, reproducing at different scales the singular intrinsic organization enclosed in itself … up to become self-aware.

References

• 1

  AgarwalK. P. (2019). A biophysical perspective on enzyme catalysis.Biochemistry58(6): 438–449. 10.1021/acs.biochem.8b01004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AllardJ.MogilnerA. (2013). Traveling waves in actin dynamics and cell motility.Curr. Opin. Cell. Biol.25107–115. 10.1016/j.ceb.2012.08.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  AlmaasE.KovacsB.VicsekT.OltvaiZ. N.BarabásiA. L. (2004). Global organization of metabolic fluxes in the bacterium Escherichia coli.Nature427839–843. 10.1038/nature02289

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  AlmaasE.OltvaiZ. N.BarabásiA. L. (2005). The activity reaction core and plasticity of metabolic networks.PLoS Comput. Biol.1:e68. 10.1371/journal.pcbi.0010068

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  AlmaasE. (2007). Biological impacts and context of network theory.J. Exp. Biol.210(Pt. 9)1548–1558. 10.1242/jeb.003731

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  AlmeiraN.GusmanS. R. (2017). Role of transcriptional bursts in cellular oscillations.J. Theor. Biol.749–56. 10.1016/j.jtbi.2017.05.029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  AmitD. J. (1992). Modeling Brain Function: The World of Attractor Neural Networks.Cambridge: Cambridge University Press.

  • Google Scholar

• 8

  AonM. A.CortassaS.O’RourkeB. (2006). The fundamental organization of cardiac mitochondria as a network of coupled oscillators.Biophys. J.914317–4327. 10.1529/biophysj.106.087817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  AonM. A.RousselM. R.CortassaS.O’RourkeB.MurrayD. B.BeckmannM.et al (2008). The scale-free dynamics of eukaryotic cells.PLoS One3:e3624. 10.1371/journal.pone.0003624

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  AonoM.ZhuL.HaraM. (2011a). Amoeba-based neurocomputing for 8-city traveling salesman problem.Int. J. Unconv. Comput.7463–480. 10.1016/j.biosystems.2013.01.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  AonoM.ZhuL.KimS. J.HaraM. (2011b). Performance enhancement of amoeba-based neurocomputer for 8-city traveling salesman problem.Proc. NOLTAA2L-B3, 104–107.

  • Google Scholar

• 12

  Araiza-OliveraD.Chiquete-FelixN.Rosas-LemusM.SampedroJ. G.PeñaA.MujicaA.et al (2013). A glycolytic metabolon in Saccharomyces cerevisiae is stabilized by F-actin.FEBS J.2803887–3905. 10.1111/febs.12387

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  AsanoY.NagasakiA.UyedaT. Q. P. (2008). Correlated waves of actin filaments and PIP3 in Dictyostelium cells.Cell Mot. Cytosk.65923–934. 10.1002/cm.20314

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  AtkinsonD. E.WaltonG. M. (1967). Adenosine triphosphate conservation in metabolic regulation. Rat liver citrate cleavage enzyme.J. Biol. Chem.2423239–3241.

  • Google Scholar

• 15

  BarrilE. F.PotterA. R. (1968). Systematic oscillations of amino acid transport in liver from rats adapted to controlled feeding schedules.J. Nutr.95228–237. 10.1093/jn/132.11.3369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  BenteleM.LavrikI.UlrichM.StösserS.HeermannD. W.KalthoffH.et al (2004). Mathematical modeling reveals threshold mechanism in CD95-induced apoptosis.J. Cell Biol.166839–851. 10.1083/jcb.200404158

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  BernardinelliY.MagistrettiP. J.ChattonJ. Y. (2004). Astrocytes generate Na+-mediated metabolic waves.Proc. Natl. Acad. Sci. U.S.A.10114937–14942. 10.1073/pnas.0405315101

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  BerridgeM. J.GalioneA. (1988). Cytosolic calcium oscillators.FASEB J.23074–3082. 10.1096/fasebj.2.15.2847949

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  BertschingerN.NatschlagerT. (2004). Real-time computation at the edge of chaos in recurrent neural networks.Neural Comput.161413–1436. 10.1162/089976604323057443

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  BierM.TeusinkB.KholodenkoB. N.WesterhoffH. V. (1996). Control analysis of glycolytic oscillations.Biophys. Chem.6215–24.

  • Google Scholar

• 21

  BoisseauR. P.VogelD.DussutourA. (2016). Habituation in non-neural organisms: evidence from slime moulds.Proc. R. Soc. B.28320160446. 10.1098/rspb.2016.0446

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  BonnerJ. T. (2010). Brainless behavior: a myxomycete chooses a balanced diet.Proc. Natl. Acad. Sci. U.S.A.1075267–5268. 10.1073/pnas.1000861107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  BorucJ.Van den DaeleH.HollunderJ.RombautsS.MylleE.HilsonP.et al (2010). Functional modules in the Arabidopsis core cell cycle binary protein-protein interaction network.Plant Cell221264–1280. 10.1105/tpc.109.073635

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  BringasC.MalainaI.Pérez-SamartínA.BoyanoM. D.FedetzM.Pérez-YarzaG.et al (2017). Long-term memory in the migration movements of enucleated Amoeba proteus.BioRxiv[Preprint]10.1101/125054

  • CrossRef
  • Google Scholar

• 25

  BrodskyV. Y. (2006). Direct cell-cell communication: a new approach derived from recent data on the nature and self-organisation of ultradian (circahoralian) intracellular rhythms.Biol. Rev. Camb. Philos. Soc.81143–162. 10.1017/S1464793105006937

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  BuescherJ. M.LiebermeisterW.JulesM.UhrM.MuntelJ.BotellaE.et al (2012). Global network reorganization during dynamic adaptations of Bacillus subtilis metabolism.Science.3351099–1103. 10.1126/science.1206871

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  ByrneK. M.MonsefiN.DawsonJ. C.DegasperiA.Bukowski-WillsJ.-C.VolinskyN.et al (2016). Bistability in the Rac1, PAK, and RhoA signaling network drives actin cytoskeleton dynamics and cell motility switches.Cell Syst.238–48. 10.1016/j.cels.2016.01.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  CarereJ.BakerP.SeahS. Y. K. (2011). Investigating the molecular determinants for substrate channeling in BphI–BphJ, an aldolase–dehydrogenase complex from the polychlorinated biphenyls degradation pathway.Biochemistry508407–8416. 10.1021/bi200960j

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  CarstenB.KarstenK. (2017). Intracellular oscillations and waves.Annu. Rev. Condens. Matt. Phys.8239–264. 10.1146/annurev-conmatphys-031016-025210

  • CrossRef
  • Google Scholar

• 30

  CastellanaM.WilsonM. Z.XuY.JoshiP.CristeaI. M.RabinowitzJ. D.et al (2014). Enzyme clustering accelerates processing of intermediates through metabolic channeling.Nat. Biotechnol.321011–1018. 10.1038/nbt.3018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  CaustonH. C. (2020). Metabolic rhythms: a framework for coordinating cellular function.Eur. J. Neurosci.511–12. 10.1111/ejn.14296

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  ChabotJ. R.PedrazaJ. M.LuitelP.van OudenaardenA. (2007). Stochastic gene expression out-of-steady-state in the cyanobacterial circadian clock.Nature4501249–1252. 10.1038/nature06395

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  ChanceB.PyeE. K.GhoshA. D.HessB. (eds) (1973). Biological and Biochemical Oscillations.New York, NY: Academic Press.

  • Google Scholar

• 34

  ChengX.FerrellJ. E.Jr. (2018). Apoptosis propagates through the cytoplasm as trigger waves.Science361607–612. 10.1126/science.aah4065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  ChiamH. K.RajagopalG. (2007). Oscillations in intracellular signaling cascades.Phys. Rev. E.75(Pt. 1)061901. 10.1103/PhysRevE.75.061901

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  ConnorK. M.GraceyA. Y. (2011). Circadian cycles are the dominant transcriptional rhythm in the intertidal mussel Mytilus californianus.Proc. Natl. Acad. Sci. U.S.A.10816110–16115. 10.1073/pnas.1111076108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  ConradoR. J.VarnerJ. D.DeLisaM. P. (2008). Engineering the spatial organization of metabolic enzymes: mimicking nature’s synergy.Curr. Opin. Biotechnol.19492–499. 10.1016/j.copbio.2008.07.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  DanoS.SorensenP.HynneF. (1999). Sustained oscillations in living cells.Nature.402320–322. 10.1038/46329

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  De ForestM.WheelerC. J. (1999). Coherent oscillations in membrane potential synchronize impulse bursts in central olfactory neurons of the crayfish.J. Neurophysiol.811231–1241. 10.1152/jn.1999.81.3.1231

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  De la FuenteI. M.MartínezL.VeguillasJ. (1995). Dynamic behavior in glycolytic oscillations with phase shifts.Biosystems351–13.

  • Google Scholar

• 41

  De la FuenteI. M.MartínezL.VeguillasJ.AguirregabiriaJ. M. (1996a) Intermittency Route to Chaos. Quasiperiodicity route to chaos in a biochemical system.Biophys. J712375–2379. 10.1016/S0006-3495(96)79431-6

  • CrossRef
  • Google Scholar

• 42

  De la FuenteI. M.MartínezL.VeguillasJ. (1996b). Intermittency route to chaos in a biochemical system.Biosystems3987–92.

  • Google Scholar

• 43

  De la FuenteI. M.MartínezL.AguirregabiriaJ. M.VeguillasJ. (1998a). R/S analysis in strange attractors.Fractals695–100.

  • Google Scholar

• 44

  De la FuenteI. M.MartínezL.BenitezN.VeguillasJ.AguirregabiriaJ. M. (1998b). Persistent behavior in a phase-shift sequence of periodical biochemical oscillations.Bull. Math. Biol.60689–702.

  • Google Scholar

• 45

  De la FuenteI. M.MartínezL.VeguillasJ.AguirregabiriaJ. M. (1998c). Coexistence of multiple periodic and chaotic regimes in biochemical oscillations.Acta Biotheor.4637–51.

  • Google Scholar

• 46

  De la FuenteI. M.BenítezN.SantamaríaA.VeguillasJ.AguirregabiriaJ. M. (1999b). Persistence in metabolic nets.Bull. Mathemat. Biol61573–595. 10.1006/bulm.1999.0103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  De la FuenteI. M. (1999). Diversity of temporal self-organized behaviors in a biochemical system.BioSystems5083–97. 10.1016/s0303-2647(98)00094-x

  • CrossRef
  • Google Scholar

• 48

  De la FuenteI. M.MartínezL.AguirregabiriaJ. M.VeguillasJ.IriarteM. (1999a). Long-range correlations in the phase-shifts of numerical simulations of biochemical oscillations and in experimental cardiac rhythms.J. Biol. Syst.7113–130.

  • Google Scholar

• 49

  De la FuenteI. M.MartínezL.Pérez-SamartínA. L.OrmaetxeaL.AmezagaC.Vera-LópezA. (2008). Global self-organization of the cellular metabolic structure.PLoS One3:e3100. 10.1371/journal.pone.0003100

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  De la FuenteI. M.VadilloF.Pérez-PinillaM.Vera-LópezA.VeguillasJ. (2009). The number of catalytic elements is crucial for the emergence of metabolic cores.PLoS One4:e7510.

  • Google Scholar

• 51

  De la FuenteI. M. (2010). Quantitative analysis of cellular metabolic dissipative, self-organized structures.Int. J. Mol. Sci113540–3599. 10.3390/ijms11093540

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  De la FuenteI. M.VadilloF.Pérez-SamartínA.Pérez-PinillaM.BidaurrazagaJ.Vera-LópezA. (2010). Global self-regulations of the cellular metabolic structure.PLoS One5:e9484. 10.1371/journal.pone.0009484

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  De la FuenteI. M.CortésJ.Pérez-PinillaM.Ruiz-RodríguezV.VeguillasJ. (2011). The metabolic core and catalytic switches are fundamental elements in the self-regulation of the systemic metabolic structure of cells.PLoS One6:e2722. 10.1371/journal.pone.0027224

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  De la FuenteI. M.CortésJ. M. (2012). Quantitative analysis of the effective functional structure in yeast glycolysis.PLoS One7:e30162. 10.1371/journal.pone.0030162

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  De la FuenteI. M.CortesJ. M.PeltaD. A.VeguillasJ. (2013). Attractor metabolic networks.PLoS One8:e58284. 10.1371/journal.pone.0058284

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  De la FuenteI. M. (2014). “Metabolic dissipative structures,” in Systems Biology of Metabolic and Signaling Networks: Energy, Mass and Information Transfer, edsAonM. A.SaksV.SchlattnerU. (New York, NY: Springer Books) 179–212.

  • Google Scholar

• 57

  De la FuenteI. M.CortésJ. M.ValeroE.DesrochesM.RodriguesS.MalainaI.et al (2014). On the dynamics of the adenylate energy system: homeorhesis vs homeostasis.PLoS One9:e108676. 10.1371/journal.pone.0108676

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  De la FuenteI. M. (2015). Elements of the cellular metabolic structure.Front. Mol. Biosc2:16. 10.3389/fmolb.2015.00016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  De la FuenteI. M.BringasC.MalainaI.FedetzM.Carrasco-PujanteJ.MoralesM.et al (2019a). Evidence of conditioned behavior in amoebae.Nat. Commun.10:3690. 10.1038/s41467-019-11677-w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  De la FuenteI. M.BringasC.MalainaI.RegnerB.Pérez-SamartinA.BoyanoM. D.et al (2019b). The nucleus does not significantly affect the migratory trajectories of amoeba in two-dimensional environments.Sci. Rep.9:16369. 10.1038/s41598-019-52716-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  De la FuenteI. M.LópezJ. I. (2020). Cell motility and cancer.Cancers122177. 10.3390/cancers12082177

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  DekhuijzenA.BagustJ. (1996). Analysis of neural bursting: nonrhythmic and rhythmic activity in isolated spinal cord.J. Neurosci. Methods67141–147. 10.1016/S0165-0270(96)00033-7

  • CrossRef
  • Google Scholar

• 63

  DenekeV. E.MelbingerA.VergassolaM.Di TaliaS. (2016). Waves of Cdk1 activity in S phase synchronize the cell cycle in Drosophila embryos.Dev. Cell38399–412. 10.1016/j.devcel.2016.07.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  DoddM. S.PapineauD.GrenneT.SlackJ. F.RittnerM.PirajnoF.et al (2017). Evidence for early life in Earth’s oldest hydrothermal vent precipitates.Nature54360–64. 10.1038/nature21377

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  DouT.-Y.LuanH.-W.GeG.-B.DongM.-M.ZouH.-F.HeY.-Q.et al (2015). Functional and structural properties of a novel cellulosome-like multienzyme complex: efficient glycoside hydrolysis of water-insoluble 7-xylosyl-10-deacetylpaclitaxel.Sci. Rep.5:13768. 10.1038/srep13768

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  Dunaway-MarianoD. (2008). Enzyme function discovery.Structure161599–1600. 10.1016/j.str.2008.10.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  DussutourA.LattyT.BeekmanM.SimpsonS. J. (2010). Amoeboid organism solves complex nutritional challenges.Proc. Natl. Acad. Sci. U.S.A.1074607–4611. 10.1073/pnas.0912198107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  EbelingW.UlbrichtH. (eds) (1986). Self-Organization by Non-Linear Irreversible Processes.Berlin: Springer-Verlag.

  • Google Scholar

• 69

  EinsteinA. (1905). Zum gegenwärtigen Stand des Strahlungsproblesm.Phys. Zeits.10185–193.

  • Google Scholar

• 70

  FenimoreP. W.FrauenfelderH.McMahonB. H.YoungR. D. (2004). Bulk-solvent and hydration-shell fluctuations, similar to alpha- and beta-fluctuations in glasses, control protein motions and functions.P. Natl. Acad. Sci. U.S.A.10114408–14413. 10.1073/pnas.0405573101

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  FuentesJ. M.PascualM. R.SalidoG.SolerJ. A. (1994). Oscillations in rat liver cytosolic enzyme activities of the urea cycle.Arch. Int. Physiol. Biochim Biophys.102237–241. 10.3109/13813459409003936

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  Garmendia-TorresC.GoldbeterA.JacquetM. (2007). Nucleocytoplasmic oscillations of the yeast transcription factor Msn2: evidence for periodic PKA activation.Curr. Biol.171044–1049. 10.1016/j.cub.2007.05.032

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  GavinA. C.BoscheM.KrauseR.GrandiP.MarziochM.BauerA.et al (2002). Functional organization of the yeast proteome by systematic analysis of protein complexes.Nature415141–147. 10.1038/415141a

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  GavinA. C.Superti-FurgaG. (2003). Protein complexes and proteome organization from yeast to man.Curr. Opin. Chem. Biol.721–27. 10.1016/S1367-5931(02)00007-8

  • CrossRef
  • Google Scholar

• 75

  GennarinoV. A.D’AngeloG.DharmalingamG.FernandezS.RussolilloG.SangesR.et al (2012). Identification of microRNA-regulated gene networks by expression analysis of target genes.Genome Res.221163–1172. 10.1101/gr.130435.111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  GerykJ.SlaninaF. (2013). Modules in the metabolic network of E. coli with regulatory interactions.Int. J. Data Min. Bioinform.8188–202. 10.1504/IJDMB.2013.055500

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  GettyL.PanteleonA. E.MittelmanS. D.DeaM. K.BergmanR. N. (2000). Rapid oscillations in omental lipolysis are independent of changing insulin levels in vivo.J. Clin. Invest.106421–430. 10.1172/JCI7815

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  Getty-KaushikL.RichardA. T.CorkeyB. E. (2005). Free fatty acid regulation of glucose-dependent intrinsic oscillatory lipolysis in perifused isolated rat adipocytes.Diabetes54629–637. 10.2337/diabetes.54.3.629

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  GoldbeterA. (2002). Computational approaches to cellular rhythms.Nature420238–245. 10.1038/nature01259

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  GoldbeterA. (2007). Biological rhythms as temporal dissipative structures.Adv. Chem. Phys.135253–295.

  • Google Scholar

• 81

  GoldbeterA. (2018). Dissipative structures in biological systems: bistability, oscillations, spatial patterns and waves.Philos. Trans. R. Soc. Lond. Ser. A Mat. Phys. Eng. Sci.37620170376. 10.1098/rsta.2017.0376

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  GonteroB.LebretonS.GracietE. (2018). “Multienzyme complexes involved in the benson–calvin cycle and in fatty acid metabolism,” in Annual Plant Reviews Book Series, edsAndrewC. A.WilliamA. L.MichaelT. M. (Weinheim: Wiley-VCH), 125–157.

  • Google Scholar

• 83

  GonzeD.HalloyJ.GoldbeterA. (2004). Stochastic models for circadian oscillations: Emergence of a biological rhythm.In.t J. Quantum Chem.98228–238.

  • Google Scholar

• 84

  GrahamD. M.AndersenT.SharekL.UzerG.RothenbergK.HoffmanB. D.et al (2018). Enucleated cells reveal differential roles of the nucleus in cell migration, polarity, and mechanotransduction.J. Cell Biol.217895–914. 10.1083/jcb.201706097

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  GrecoT. M.CristeaI. M. (2016). The biochemical evolution of protein complexes.Trends Biochem. Sci.414–6. 10.1016/j.tibs.2015.11.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  GuthrieP. B.KnappenbergerJ.SegalM.BennettM. V. L.CharlesA. C.KaterS. B.et al (1999). ATP released from astrocytes mediates glial calcium waves.J. Neurosci.19520–528. 10.1523/JNEUROSCI.19-02-00520.1999

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  GuttmanR.LewisS.RinzelJ. (1980). Control of repetitive firing in squid axon membrane as a model for a neuroneoscillator.J. Physiol.305377–395. 10.1113/jphysiol.1980.sp013370

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  HalleyJ. D.WinklerD. A. (2008). Consistent concepts of self-organization and self-assembly.Complexity1410–17. 10.1002/cplx.20235

  • CrossRef
  • Google Scholar

• 89

  HansM. A.HeinzleE.WittmannC. (2003). Free intracellular amino acid pools during autonomous oscillations in Saccharomyces cerevisiae.Biotechnol. Bioeng.82143–151. 10.1002/bit.10553

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  HaoD.RenC.LiC. (2012). Revisiting the variation of clustering coefficient of biological networks suggests new modular structure.BMC Syst. Biol.6:34. 10.1186/1752-0509-6-34

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91

  HarringtonH. A.HoK. L.GhoshS.TungK. C. (2008). Construction and analysis of a modular model of caspase activation in apoptosis.Theor. Biol. Med. Model.526. 10.1186/1742-4682-5-26

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  HartigK.BeckE. (2005). Endogenous cytokinin oscillations control cell cycle progression of tobacco BY-2 cells.Plant. Biol.733–40. 10.1055/s-2004-830474

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 93

  HeltbergM. L.KrishnaS.JensenM. H. (2019). On chaotic dynamics in transcription factors and the associated effects in differential gene regulation.Nat. Commun.1071. 10.1038/s41467-018-07932-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 94

  HertzJ. A.KroghA. S.PalmerR. G. (1991). Introduction to the Theory of Neural Computation. Santa Fe Institute Series.Boca Raton: CRC Press352. 10.1201/9780429499661

  • CrossRef
  • Google Scholar

• 95

  HilarioE.CaulkinsB. G.HuangY.-M. M.YouW.ChangC.-E. A.MuellerL. J.et al (2016). Visualizing the tunnel in tryptophan synthase with crystallography: insights into a selective filter for accommodating indole and rejecting water.Biochim. Biophys. Acta1864268–279. 10.1016/j.bbapap.2015.12.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96

  HoY.GruhlerA.HeilbutA.BaderG. D.MooreL.AdamsS. L.et al (2002). Systematic identification of protein complexes in Saccharomyces cerevisiae by mass spectrometry.Nature415180–183. 10.1038/415180a

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 97

  HolzG. G.HeartE.LeechC. A. (2008). Synchronizing Ca2+ and cAMP oscillations in pancreatic beta cells: a role for glucose metabolism and GLP-1 receptors? Focus on “Regulation of cAMP dynamics by Ca2+ and G protein-coupled receptors in the pancreatic β-cell: a computational approach”.Am. J. Physiol. Cell Physiol.294C4–C6. 10.1152/ajpcell.00522.2007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 98

  HopfieldJ. J. (1982). Neural networks and physical systems with emergent collective computational abilities.Proc. Nat. Acad. Sci. U.S.A.792554–2558. 10.1073/pnas.79.8.2554

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  HoskinsA. A.AnandR.EalickE. E.StubbeJ. (2004). the formylglycinamide ribonucleotide amidotransferase complex from Bacillus subtilis: metabolite-mediated complex formation.Biochemistry4310314–10327. 10.1021/bi049127h

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 100

  HsuJ. T.PengC. H.HsiehW. P.LanC. Y.TangC. Y. (2011). A novel method to identify cooperative functional modules: study of module coordination in the Saccharomyces cerevisiae cell cycle.BMC Bioinform.12:281. 10.1186/1471-2105-12-281

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 101

  HungerbuehlerA. K.PhilippsenP.GladfelterA. S. (2007). Limited functional redundancy and oscillation of cyclins in multinucleated Ashbya gossypii fungal cells.Eukaryot. Cell6473–486. 10.1128/EC.00273-06

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 102

  HymanA. A.WeberC. A.JülicherF. (2014). Liquid-liquid phase separation in biology.Ann. Rev. Cell Dev. Biol.3039–58. 10.1146/annurev-cellbio-100913-013325

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 103

  IshiiK.HiroseK.IinoM. (2006). Ca2+ shuttling between endoplasmic reticulum and mitochondria underlying Ca2+ oscillations.EMBO Rep.7390–396. 10.1038/sj.embor.7400620

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 104

  IzardT.AevarssonA.AllenM. D.WestphalA. H.PerhamR. N.de KokA.et al (1999). Principles of quasi-equivalence and Euclidean geometry govern the assembly of cubic and dodecahedral cores of pyruvate dehydrogenase complexes.Proc. Natl. Acad. Sci. U.S.A.961240–1245. 10.1073/pnas.96.4.1240

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 105

  JaffeL. F. (2002). On the conservation of fast calcium wave speeds.Cell Calcium32217–229. 10.1016/s0143416002001574

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 106

  JaffeL. F. (2008). Calcium waves.Philos. Trans. R. Soc. Lond. B Biol. Sci.3631311–1317. 10.1098/rstb.2007.2249

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 107

  JulesM.FrancoisJ.ParrouJ. L. (2005). Autonomous oscillations in Saccharomyces cerevisiae during batch cultures on trehalose.FEBS J.2721490–1500. 10.1111/j.1742-4658.2005.04588.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 108

  JungT.GruneT. (2012). Structure of the proteasome.Prog. Mol. Biol. Transl. Sci.1091–39. 10.1016/b978-0-12-397863-9.00001-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 109

  KaltenbachH. M.StellingJ. (2012). Modular analysis of biological networks.Adv. Exp. Med. Biol.7363–17. 10.1007/978-1-4419-7210-1_1

  • CrossRef
  • Google Scholar

• 110

  KassL.BrayW. O. (1996). Kinetic model for phototransduction and G-protein enzyme cascade: understanding quantal bumps during inhibition of CaM-KII or PP2B.J. Photochem. Photobiol. B Biol.35105–113.

  • Google Scholar

• 111

  KastritisP. L.GavinA. C. (2018). Enzymatic complexes across scales.Essays Biochem.62501–514. 10.1042/EBC20180008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 112

  KauffmanS.JohnsenS. (1991). Coevolution to the edge of chaos: Coupled fitness landscapes, poised states, and coevolutionary avalanches.J. Theo. Biol.149467–505. 10.1016/S0022-5193(05)80094-3

  • CrossRef
  • Google Scholar

• 113

  KellerM.PiedrafitaG.RalserM. (2015). The widespread role of non-enzymatic reactions in cellular metabolism.Curr. Opin. Biotechnol34153–161. 10.1016/j.copbio.2014.12.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 114

  KimG.YangJ.JangJ.ChoiJ.RoeA. J.OlwynB.et al (2020). Aldehyde-alcohol dehydrogenase undergoes structural transition to form extended spirosomes for substrate channeling.Commun. Biol.3:298. 10.1038/s42003-020-1030-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 115

  KindzelskiiA. L.ZhouM. J.HauglandP. R.BoxerA. L.PettyR. H. (1998). Oscillatory pericellular proteolysis and oxidant deposition during neutrophil locomotion.Biophys. J.7490–97. 10.1016/S0006-3495(98)77770-7

  • CrossRef
  • Google Scholar

• 116

  KindzelskiiA. L.PettyH. R. (2002). Apparent role of traveling metabolic waves in oxidant release by living neutrophils.Proc. Natl. Acad. Sci. U.S.A.999207–9212. 10.1073/pnas.132630999

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 117

  KleveczR. R.MurrayD. B. (2001). Genome wide oscillations in expression. Wavelet analysis of time series data from yeast expression arrays uncovers the dynamic architecture of phenotype.Mol. Biol. Rep.2873–82. 10.1023/A:1017909012215

  • CrossRef
  • Google Scholar

• 118

  KleveczR. R.BolenJ.ForrestG.MurrayD. B. (2004). A genomewide oscillation in transcription gates DNA replication and cell cycle.Proc. Natl. Acad. Sci. U.S.A.1011200–1205. 10.1073/pnas.0306490101

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 119

  KlimontovichY. L. (1999). Entropy and information of open systems.Phys. Uspekhi.42375–384. 10.1515/JNETDY.2008.006

  • CrossRef
  • Google Scholar

• 120

  KohnhorstC. L.KyoungM.JeonM.SchmittD. L.KennedyE. L.RamirezJ.et al (2017). Identification of a multienzyme complex for glucose metabolism in living cells.J. Biol. Chem.2929191–9203. 10.1074/jbc.M117.783050

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 121

  KondepudiD.PrigogineI. (2014). Modern Thermodynamics: From Heat Engines to Dissipative Structures, 2nd Edn. Weinheim: Wiley-VCH.

  • Google Scholar

• 122

  KorzeniewskiB.ZoladzJ. A. (2001). A model of oxidative phosphorylation in mammalian skeletal muscle.Biophys. Chem.9217–34.

  • Google Scholar

• 123

  KorcsmárosT.KovácsI. A.SzalayM. S.CsermelyP. (2007). Molecular chaperones: the modular evolution of cellular networks.J. Biosci.32441–446. 10.1007/s12038-007-0043-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 124

  KruseK.JülicherF. (2005). Oscillations in cell biology.Curr. Opin. Cell Biol.1720–26. 10.1016/j.ceb.2004.12.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 125

  KühnerS.van NoortV.BettsM. J.Leo-MaciasA.BatisseC.RodeM.et al (2009). Proteome organization in a genome-reduced bacterium.Science3261235–1240. 10.1126/science.1176343

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 126

  KükenA.SommerF.Yaneva-RoderL.MackinderL. C.HöhneM.GeimerS.et al (2018). Effects of microcompartmentation on flux distribution and metabolic pools in Chlamydomonas reinhardtii chloroplasts.eLife7:e37960. 10.7554/eLife.37960

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 127

  KurzF. T.AonM. A.O’RourkeB.ArmoundasA. A. (2010). Spatio-temporal oscillations of individual mitochondria in cardiac myocytes reveal modulation of synchronized mitochondrial clusters.Proc. Natl. Acad. Sci. U.S.A.10714315–14320. 10.1073/pnas.1007562107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 128

  KuzinA. M. (1970). Structural – Metabolic Hypothesis in Radiobiology.Moscow: Nauka.

  • Google Scholar

• 129

  LahtveeP.SeimanA.ArikeL.AdambergK.ViluR. (2014). Protein turnover forms one of the highest maintenance costs in Lactococcus lactis.Microbiology160(Pt .7) 1501–1512. 10.1099/mic.0.078089-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 130

  LangeG.MandelkowE. M.JaglaA.MandelkowE. (1988). Tubulin oligomers and microtubule oscillations. antagonistic role of microtubule stabilizers and destabilizers.Eur. J. Biochem.17861–69. 10.1111/j.1432-1033.1988.tb14429.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 131

  LattyT.BeekmanM. (2011). Speed-accuracy trade-offs during foraging decisions in the acellular slime mould Physarum polycephalum.Proc. Biol. Sci.278539–545. 10.1098/rspb.2010.1624

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 132

  LloydD.EshanthaL.SalgadoJ.TurnerM. P.MurrayD. B. (2002). Respiratory oscillations in yeast: clock-driven mitochondrial cycles of energization.FEBS Lett.51941–44. 10.1016/S0014-5793(02)02704-7

  • CrossRef
  • Google Scholar

• 133

  LloydD.MurrayD. B. (2005). Ultradian metronome: timekeeper for orchestration of cellular coherence.Trends Biochem. Sci.30373–377. 10.1016/j.tibs.2005.05.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 134

  LloydD.MurrayD. B. (2006). The temporal architecture of eukaryotic growth.FEBS Lett.5802830–2835. 10.1016/j.febslet.2006.02.066

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 135

  MackintoshN. J. (1983). Conditioning and Associative Learning.Oxford: Clarendon Press.

  • Google Scholar

• 136

  MaierT.JenniS.BanN. (2006). Architecture of mammalian fatty acid synthase at 4.5 A resolution.Science3111258–1262. 10.1126/science.1123248

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 137

  MarkevichN. I.TsyganovM. A.HoekJ. B.KholodenkoB. N. (2006). Long−range signaling by phosphoprotein waves arising from bistability in protein kinase cascades.Mol. Syst. Biol.261. 10.1038/msb4100108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 138

  MárquezS.CrespoP.CarliniV.Garbarino-PicoE.BalerR.CaputtoB. L.et al (2004). The metabolism of phospholipids oscillates rhythmically in cultures of fibroblasts and is regulated by the clock protein PERIOD 1.FASEB J.18519–521. 10.1096/fj.03-0417fje

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 139

  MisteliT. (2009). Self-organization in the genome.Proc. Natl. Acad. Sci. U.S.A.1066885–6886. 10.1073/pnas.0902010106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 140

  MiyajiT.OhnishiI. (2008). Physarum can solve the shortest path problem on riemannian surface mathematically rigourously.Int. J. Pure App. Math.47353–369.

  • Google Scholar

• 141

  MøllerA. C.HauserM. J. B.OlsenL. F. (1998). Oscillations in peroxidase-catalyzed reactions and their potential function in vivo.Biophys. Chem.7263–72. 10.1016/S0301-4622(98)00123-9

  • CrossRef
  • Google Scholar

• 142

  MogilevskayaE.DeminO.GoryaninI. (2006). Kinetic model of mitochondrial krebs cycle: unraveling the mechanism of salicylate hepatotoxic effects.J. Biol. Phys.32245–271.

  • Google Scholar

• 143

  MongeC.GrichineA.RostovtsevaT.SackettP.SaksV. A. (2009). Compartmentation of ATP in cardiomyocytes and mitochondria. kinetic studies and direct measurements.Biophys. J.96:241a. 10.1016/j.bpj.2008.12.1188

  • CrossRef
  • Google Scholar

• 144

  MurrayD.BeckmannM.KitanoH. (2007). Regulation of yeast oscillatory dynamics.Proc. Natl. Acad. Sci. U.S.A.1042241–2246. 10.1073/pnas.0606677104

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 145

  MutoA.KoteraM.TokimatsuT.NakagawaZ.GotoS.KanehisaM. (2013). Modular architecture of metabolic pathways revealed by conserved sequences of reactions.J. Chem. Inf. Model.53613–622. 10.1021/ci3005379

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 146

  NakagakiT.YamadaH.TothA. (2000). Maze-solving by an amoeboid organism.Nature407:470. 10.1038/35035159

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 147

  NakagakiT. (2001). Smart behavior of true slime mold in a labyrinth.Res. Microbiol.152767–770. 10.1016/S0923-2508(01)01259-1

  • CrossRef
  • Google Scholar

• 148

  NakagakiT.KobayashiR.UedaT.NishiuraY. (2004a). Obtaining multiple separate food sources: behavioural intelligence in the Physarum plasmodium.Proc. Biol. Sci.2712305–2310. 10.1098/rspb.2004.2856

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 149

  NakagakiT.YamadaH.HaraM. (2004b). Smart network solutions in an amoeboid organism.Biophys. Chem.1071–5. 10.1016/S0301-4622(03)00189-3

  • CrossRef
  • Google Scholar

• 150

  NakamuraN.WeiJ. H.SeemannJ. (2012). Modular organization of the mammalian Golgi apparatus.Curr. Opin. Cell. Biol.24467–474. 10.1016/j.ceb.2012.05.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 151

  NakamuraY. (2018). Membrane Lipid oscillation: an emerging system of molecular dynamics in the plant membrane.Plant Cell Physiol.59441–447. 10.1093/pcp/pcy023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 152

  NegreteJ.PumirA.HsuH.-F.WestendorfC.TarantolaM.BetaC.et al (2016). Noisy oscillations in the actin cytoskeleton of chemotactic amoeba.Phys. Rev. Lett.117148102. 10.1103/PhysRevLett.117.148102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 153

  NephS.StergachisA. B.ReynoldsA.SandstromR.BorensteinE.StamatoyannopoulosJ. A. (2012). Circuitry and dynamics of human transcription factor regulatory networks.Cell1501274–1286. 10.1016/j.cell.2012.04.040

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 154

  NewmanE. A. (2001). Propagation of intercellular calcium waves in retinal astrocytes and müller cells.J. Neurosci.212215–2223. 10.1523/JNEUROSCI.21-07-02215.2001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 155

  NiebelB.LeupoldS.HeinemannM. (2019). An upper limit on gibbs energy dissipation governs cellular metabolism.Nat. Metab.1125–132. 10.1038/s42255-018-0006-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 156

  NicolisG.PrigogineI. (1977). Self-Organization in Nonequilibrium Systems. From Dissipative Structures to Order Through Fluctuations.New York, NY: Wiley.

  • Google Scholar

• 157

  NobleD.JablonkaE.JoynerM. J.MüllerG. B.OmholtS. W. (2014). Evolution evolves: physiology returns to centre stage.J. Physiol.592(Pt. 11)2237–2244. 10.1113/jphysiol.2014.273151

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 158

  NoletF. E.VanderveldeA.VanderbekeA.PiñerosL.ChangJ. B.GelensL. (2020). Nuclei determine the spatial origin of mitotic waves.eLife9:e52868. 10.7554/eLife.52868

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 159

  NurseP. (2008). Life, logic and information.Nature454424–426. 10.1038/454424a

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 160

  MathiesonT.FrankenH.KosinskiJ.KurzawaN.ZinnN.SweetmanG.et al (2018). Systematic analysis of protein turnover in primary cells.Nat. Commun.9:689. 10.1038/s41467-018-03106-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 161

  OlsenL. F.DegnH. (1985). Chaos in biological systems.Q. Rev. Biophys.18165–225. 10.1017/s0033583500005175

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 162

  OvádiJ.SrereP. A. (2000). Macromolecular compartmentation and channeling.Int. Rev. Cytol.192255–280. 10.1016/S0074-7696(08)60529-X

  • CrossRef
  • Google Scholar

• 163

  OvádiJ.SaksV. (2004). On the origin of the ideas of intracellular compartmentation and organized metabolic systems.Mol. Cell. Biochem.2565–12. 10.1023/B:MCBI.0000009855.14648.2c

  • CrossRef
  • Google Scholar

• 164

  OzalpV. C.PedersenT. R.NielsenL. J.OlsenL. F. (2010). Time-resolved measurements of intracellular ATP in the yeast Saccharomyces cerevisiae using a new type of nanobiosensor.J. Biol. Chem.28537579–37588. 10.1074/jbc.M110.155119

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 165

  PangC. N.KrycerJ. R.LekA.WilkinsM. R. (2008). Are protein complexes made of cores, modules and attachments?Proteomics8425–434. 10.1002/pmic.200700801

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 166

  PavlovI. P. (1927). Conditioned Reflexes: An Investigation of the Physiological Activity of the Cerebral Cortex.Oxford: Oxford Univ. Press.

  • Google Scholar

• 167

  PêkalskiJ.ZukP. J.KochańczykM.JunkinM.KelloggR.TayS.et al (2013). Spontaneous NF-κB activation by autocrine TNFα signaling: a computational analysis.PLoS One8:e78887. 10.1371/journal.pone.0078887

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 168

  PerettoP. (1992). An Introduction to the Modeling of Neural Networks.Cambridge: Cambridge University Press.

  • Google Scholar

• 169

  PetrovicA.MosalagantiS.KellerJ.MattiuzzoM.OverlackK.KrennV. (2014). Modular assembly of RWD domains on the Mis12 complex underlies outer kinetochore organization.Mol. Cell53591–605. 10.1016/j.molcel.2014.01.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 170

  PettyH. R.WorthR. G.KindzelskiiA. L. (2000). Imaging sustained dissipative patterns in the metabolism of individual living cells.Phys. Rev. Lett.842754–2757. 10.1103/PhysRevLett.84.2754

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 171

  PlacantonakisD. G.WelshJ. P. (2001). Two distinct oscillatory states determined by the NMDA receptor in rat inferior olive.J. Physiol.534(Pt. 1)123–140. 10.1111/j.1469-7793.2001.t01-1-00123.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 172

  PostmaM.BosgraafL.LooversH. M.Van HaastertP. J. (2004). Chemotaxis: signalling modules join hands at front and tail.EMBO Rep.535–40. 10.1038/sj.embor.7400051

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 173

  RadzickaA.WolfendenR. (1995). A proficient enzyme.Science26790–93. 10.1126/science.7809611

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 174

  RamanathanA.SavolA.BurgerV.ChennubhotlaC. S.AgarwalP. K. (2014). Protein conformational populations and functionally relevant substates.Acc. Chem. Res.47149–156. 10.1021/ar400084s

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 175

  RavaszE.SomeraA. L.MongruD. A.OltvaiZ. N.BarabásiA. L. (2002). Hierarchical organization of modularity in metabolic networks.Science2971551–1555. 10.1126/science.1073374

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 176

  RenganR.OmannG. M. (1999). Regulation of oscillations in filamentous actin content in polymorphonuclear leukocytes stimulated with leukotriene B4 and platelet-activating factor.Biochem. Biophys. Res. Commun.262479–486. 10.1006/bbrc.1999.1222

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 177

  RießB.GrötschR. K.BoekhovenJ. (2020). The design of dissipative molecular assemblies driven by chemical reaction cycles.Chemistry6552–578. 10.1016/j.chempr.2019.11.008

  • CrossRef
  • Google Scholar

• 178

  RingeD.PetskoG. A. (2009). How enzymes work.Science3201428–1429. 10.1126/science.1159747

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 179

  RomashkoD. N.MarbanE.O’RourkeB. (1998). Subcellular metabolic transients and mitochondrial redox waves in heart cells.Proc. Natl. Acad. Sci. U.S.A.951618–1623. 10.1073/pnas.95.4.1618

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 180

  RossA. B.LangerJ. D.JovanovicM. (2020). Proteome turnover in the spotlight: approaches, applications, and perspectives.Mol. Cell Proteomics.20:00016. 10.1074/mcp.R120.002190

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 181

  RousselM. R.IvlevA. A.IgamberdievA. U. (2006). Oscillations of the internal CO2 concentration in tobacco leaves transferred to low CO2.J. Plant Physiol.341188–1196. 10.1016/j.jplph.2006.08.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 182

  SabziNezhadA.JaliliS. (2020). DPCT: a dynamic method for detecting protein complexes From TAP-aware weighted PPI network.Front. Genet11:567. 10.3389/fgene.2020.00567

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 183

  SaigusaT.TeroA.NakagakiT.KuramotoY. (2008). Amoebae anticipate periodic events.Phys. Rev. Lett.100:018101. 10.1103/PhysRevLett.100.01810

  • CrossRef
  • Google Scholar

• 184

  SaksV. A.MongeC.AnmannT.DzejaP. (2007). “integrated and organized cellular energetic systems: theories of cell energetics, compartmentation and metabolic channeling,” in Molecular System Bioenergetics. Energy for Life, ed.SaksV. A. (Weinheim: Wiley-VCH), 59–110.

  • Google Scholar

• 185

  SaksV. A.MongeC.GuzunR. (2009). Philosophical basis and some historical aspects of systems biology: From Hegel to Noble – Applications for bioenergetic research.Int. J. Mol. Sci.101161–1192. 10.3390/ijms10031161

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 186

  San RománM.CancelaH.AcerenzaL. (2014). Source and regulation of flux variability in Escherichia coli.BMC Syst. Biol.8:67. 10.1186/1752-0509-8-67

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 187

  Sánchez-ArmássS.SennouneS. R.MaitiD.OrtegaF.Martínez-ZaguilaR. (2006). Spectral imaging microscopy demonstrates cytoplasmic pH oscillations in glial cells.Am. J. Physiol. Cell Physiol.290C524–C538. 10.1152/ajpcell.00290.2005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 188

  SantillánM.MackeyM. C. (2008). Quantitative approaches to the study of bistability in the lac operon of Escherichia coli.J. R. Soc. Interf.5(Suppl 1)S29–S39. 10.1098/rsif.2008.0086.focus

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 189

  SartoriP.LeiblerS. (2020). Lessons from equilibrium statistical physics regarding the assembly of protein complexes.P. Natl. Acad. Sci.117114–120. 10.1073/pnas.1911028117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 190

  ScemesE.GiaumeC. (2006). Astrocyte calcium waves: what they are and what they do.Glia54716–725. 10.1002/glia.20374

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 191

  SchreiberT. (2000). Measuring information transfer.Phys. Rev. Lett.85461–464. 10.1103/PhysRevLett.85.461

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 192

  SchubertV. (2014). RNA polymerase II forms transcription networks in rye and Arabidopsis nuclei and its amount increases with endopolyploidy.Cytogenet. Genome Res.14369–77. 10.1159/000365233

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 193

  SerreV.GuyH.LiuX.PenverneB.HervéG.EvansD. (1998). Allosteric regulation and substrate channeling in multifunctional pyrimidine biosynthetic complexes: analysis of isolated domains and yeast-mammalian chimeric proteins.J. Mol. Biol.281363–377. 10.1006/jmbi.1998.1856

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 194

  ShankaranH.IppolitoD. L.ChrislerW. B.ResatH.BollingerN.OpreskoL. K.et al (2009). Rapid and sustained nuclear–cytoplasmic ERK oscillations induced by epidermal growth factor.Mol. Syst. Biol.5:332. 10.1038/msb.2009.90

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 195

  ShimizuT. S.TuY.BergH. C. (2010). A modular gradient-sensing network for chemotaxis in Escherichia coli revealed by responses to time-varying stimuli.Mol. Syst. Biol.6382. 10.1038/msb.2010.37

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 196

  SlabyO.LebiedzD. (2009). Oscillatory NAD(P)H waves and calcium oscillations in neutrophils? A modeling study of feasibility.Biophys. J.96417–428. 10.1016/j.bpj.2008.09.044

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 197

  SmrcinováM.SørensenP. G.KrempaskyJ.BalloP. (1998). Chaotic oscillations in a chloroplast system under constant illumination.Int. J. Bifurcat. Chaos.82467–2470. 10.1142/S0218127498001984

  • CrossRef
  • Google Scholar

• 198

  SrereP. A. (1972). “Is there an organization of krebs cycle enzymes in the mitochondrial matrix?,” in Energy Metabolism and the Regulation of Metabolic Processes in Mitochondria, edsHansonR. W.MehlmanW. A. (New York, NY: Academic Press), 79–91.

  • Google Scholar

• 199

  SrereP. A. (1985). The metabolon.Trends Biochem. Sci.10109–110. 10.1016/0968-0004(85)90266-X

  • CrossRef
  • Google Scholar

• 200

  SrereP. A. (1987). Complexes of sequential metabolic enzymes.Annu. Rev. Biochem.5689–124. 10.1146/annurev.bi.56.070187.000513

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 201

  StryerL.BergJ. M.TymoczkoJ. L. (2002). Biochemistry, 5th Edn. San Francisco, CA: W.H. Freeman.

  • Google Scholar

• 202

  SunM.LeyhT. S. (2006). Channeling in sulfate activating complexes.Biochemistry4511304–11311. 10.1021/bi060421e

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 203

  SvedružićZ. M.OdorèićI.ChangC. H.SvedružićD. (2020). Substrate channeling via a transient protein-protein complex: the case of D-Glyceraldehyde-3-phosphate dehydrogenase and L-lactate dehydrogenase.Sci. Rep.10:10404. 10.1038/s41598-020-67079-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 204

  SweetloveL. J.FernieA. R. (2018). The role of dynamic enzyme assemblies and substrate channelling in metabolic regulation.Nat. commun.9:2136. 10.1038/s41467-018-04543-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 205

  TeroA.TakagiS.SaigusaT.ItoK.BebberD. P.FrickerM. D.et al (2010). Rules for biologically inspired adaptive network design.Science327439–442. 10.1126/science.1177894

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 206

  ThodenJ. B.HoldenH. M.WesenbergG.RaushelF. M.RaymentI. (1997). Structure of carbamoyl phosphate synthetase: a journey of 96 Å from substrate to product.Biochemistry366305–6316. 10.1021/bi970503q

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 207

  TianB.NowakD. E.BrasierA. R. (2005). A TNF-induced gene expression program under oscillatory NF-κB control.BMC Genomics6:137. 10.1186/1471-2164-6-137

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 208

  TonozukaH.WangJ.MitsuiK.SaitoT.HamadaY.TsurugiK. (2001). Analysis of the upstream regulatory region of the GTS1 gene required for its oscillatory expression.J. Biochem.130589–595. 10.1093/oxfordjournals.jbchem.a003023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 209

  TretterL.Adam-ViziV. (2005). Alpha-ketoglutarate dehydrogenase: a target and generator of oxidative stress.Philos. Trans. R. Soc. Lond. B. Biol. Sci.3602335–2345. 10.1098/rstb.2005.1764

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 210

  TsaiS.AmesB. D. (2009). Structural enzymology of polyketide synthases.Methods Enzymol.45917–47. 10.1016/S0076-6879(09)04602-3

  • CrossRef
  • Google Scholar

• 211

  TuY.ShimizuT. S.BergH. C. (2008). Modeling the chemotactic response of Escherichia coli to time-varying stimuli.Proc. Natl. Acad. Sci. U.S.A10514855–14860. 10.1073/pnas.0807569105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 212

  TuttleL. M.SalisH.TomshineJ.KaznessisY. N. (2005). Model-driven designs of an oscillating gene network.Biophys. J.893873–3883. 10.1529/biophysj.105.064204

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 213

  TysonJ. J. (1991). Modeling the cell division cycle: cdc2 and cyclin interactions.PNAS887328–7332.

  • Google Scholar

• 214

  TysonJ. J. (2002). “Biochemical oscillations,” in Computational Cell Biology. Interdisciplinary Applied Mathematics, Vol. 20edsFallC. P.MarlandE. S.WagnerJ. M.TysonJ. J. (New York, NY: Springer), 230–260.

  • Google Scholar

• 215

  UedaT.NakagakiT.YamadaT. (1990). Dynamic organization of ATP and birefringent fibrils during free locomotion and galvanotaxis in the plasmodium of Physarum polycephalum.J. Cell. Biol.1101097–1102. 10.1083/jcb.110.4.1097

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 216

  ValeroE.MaciáH.de la FuenteI. M.HernandezJ. A.González-SánchezM.IGarcía-CarmonaF. (2016). Modeling the under light/dark conditions.BMC Syst. Biol.10:11. 10.1186/s12918-015-0239-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 217

  von StockarU. (2010). Biothermodynamics of live cells: a tool for biotechnology and biochemical engineering.J. Non Equil. Therm.35415–475. 10.1515/JNETDY.2010.024

  • CrossRef
  • Google Scholar

• 218

  VickerM. G. (2002). Eukaryotic cell locomotion depends on the propagation of self-organized reaction-diffusion waves and oscillations of actin filament assembly.Exp. Cel. Res.27554–66. 10.1006/excr.2001.5466

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 219

  Vlasova-St LouisI.BohjanenP. R. (2011). Coordinate regulation of mRNA decay networks by GU-rich elements and CELF1.Curr. Opin. Genet. Dev.21444–451. 10.1016/j.gde.2011.03.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 220

  WatfordM. (1989). Channeling in the urea cycle: a metabolon spanning two compartments.Trends Biochem. Sci.14313–314. 10.1016/0968-0004(89)90156-4

  • CrossRef
  • Google Scholar

• 221

  WettmannL.KarstenK. (2018). The min-protein oscillations in Escherichia coli: An example of self-organized cellular protein waves.Philos. Trans. R. Soc. Lond. B. Biol. Sci37320170111. 10.1098/rstb.2017.0111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 222

  WhitesidesG. M.GrzybowskiB. (2002). Self-assembly at all scales.Science2952418–2421. 10.1126/science.1070821

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 223

  WijnenH.YoungW. M. (2006). Interplay of circadian clocks and metabolic rhythms.Annu Rev Genet.40409–448. 10.1146/annurev.genet.40.110405.090603

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 224

  WittmannC.HansM.Van WindenA. W.RasC.HeijnenJ. J. (2005). Dynamics of intracellular metabolites of glycolysis and TCA cycle during cell-cycle-related oscillation in Saccharomyces cerevisiae.Biotechnol. Bioeng.89839–847. 10.1002/bit.20408

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 225

  XiongW.FerrellJ. E. (2003). A positive-feedback-based bistable “memory module” that governs a cell fate decision.Nature426460–465. 10.1038/nature02089

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 226

  YangC. R.ShapiroB. E.HungS. P.MjolsnessE. D.HatfieldG. W. (2005). A mathematical model for the branched chain amino acid biosynthetic pathways of Escherichia coli K12.J Biol Chem.28011224–11232.

  • Google Scholar

• 227

  YanlingX.QingmingT.GuangjinChMengruX.ShaolingY.JiajiaZ.et al (2019). New insights into the circadian rhythm and its related diseases.Front. Physiol10:682. 10.3389/fphys.2019.00682

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 228

  YoonJ.SiY.NolanR.LeeK. (2007). Modular decomposition of metabolic reaction networks based on flux analysis and pathway projection.Bioinformatics232433–2440. 10.1093/bioinformatics/btm374

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 229

  ZhabotinskyA. M. (2000). Bistability in the Ca(2+) calmodulin-dependent protein kinase-phosphatase system.Biophys. J.792211–2221. 10.1016/S0006-3495(00)76469-1

  • CrossRef
  • Google Scholar

• 230

  ZhangY.FernieA. R. (2020). Metabolons, enzyme–enzyme assemblies that mediate substrate channeling, and their roles in plant metabolism.Plant Commun.2:100081. 10.1016/j.xplc.2020.100081

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 231

  ZhouL.AonM. A.AlmasT.CortassaS.WinslowR. L.O’RourkeB. (2010). A reaction-diffusion model of ros-induced ros release in a mitochondrial network.PLoS Comput. Biol.6:e1000657. 10.1371/journal.pcbi.1000657

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 232

  ZhuL.AonoM.KimS.-J.HaraM. (2011). Problem-sizes scalability of amoeba-based neurocomputer for traveling salesman problem.Proc. NOLTA45108–111. 10.34385/proc.45.A2L-B4

  • CrossRef
  • Google Scholar