Aging as a Loss of Goal-Directedness: An Evolutionary Simulation and Analysis Unifying Regeneration with Anatomical Rejuvenation

1. Introduction

1.1. Theories of Aging

Aging is a system-level, near-ubiquitous biological process characterized by morphological and functional alterations in cellular and extracellular components, resulting in a systematic decline in biological functions [1]. This degenerative process is modulated by genetics, lifestyle, and environmental factors [2]. As cellular functionality diminishes with age, there is a corresponding decrease in physiological capacities and an enhanced vulnerability to a spectrum of diseases, including cardiovascular disorders, cancer, neurodegenerative diseases, type II diabetes, and various infections. These conditions collectively compromise the quality of human life [2,3]. fig:introduction:motivation (A) schematically illustrates the typical life cycle of a complex organism comprising various stages from embryogenesis through maturity, decline, aging, and eventual death.

Several theories have been proposed to explain the aging process, that we can classify into damage-based and programmatic theories [1,4,5], with the former being more extensively studied [1,3,6,7]. Damage-based theories assert that aging results from the accumulation of molecular damage, impacting crucial cellular components like the genome, telomeres, mitochondria, and proteins, driven by imperfect repair mechanisms [6,7,8]. Conversely, programmatic theories suggest that aging is genetically predetermined, or quasi-programmed and encoded within the genome, not merely due to stochastic accumulation of damage [9,10,11]. Here, we formulate and computationally simulate a fundamentally different hypothesis based on the nature of morphogenesis as a homeodynamic system working toward anatomical setpoints: that aging is the result of a proto-cognitive disfunction among the cellular collective, and does not require damage, or evolved antagonistic mechanisms in order to occur. Instead, we propose that, as observed in animal behavior [12], cellular living systems without goals to pursue tend to degrade. Our hypothesis is that aging is the result of intrinsic dynamics of complex cybernetic (goal-directed) systems, not caused by external noise or damage. Nevertheless, it is consistent with other theories of aging [13,14] which propose additional dynamics that can lead to loss of setpoints of anatomical homeostasis, and can be accelerated by various kinds of damage.

Antagonistic pleiotropy posits that certain genes have beneficial effects in early life but detrimental effects in later life. These genes confer advantages such as increased fertility and resilience in youth, which are critical for survival and reproduction. However, as an individual ages, the same genes contribute to aging and decline because their negative effects accumulate. The theory suggests that natural selection favors these genes because of their early-life benefits, despite their harmful subsequent impacts. This results in a trade-off between early benefits and later costs [15]. Another related idea is the hyperfunction theory, which suggests that aging is not primarily a result of accumulating molecular damage or genetic programming, but rather due to a continuous and excessive operation of cellular growth processes that were beneficial during the earlier stages of life [13,14]. Central to this theory is the role of the Target of Rapamycin (TOR) pathway, which is a major regulator of growth and metabolism in cells. During the growth phase of an organism’s life, the TOR pathway promotes cellular and organismal growth by driving processes like protein synthesis, nutrient uptake, and cell division. These functions are vital for development and reproductive success. However, aging arises when these growth-promoting activities are not appropriately down-regulated after the organism has passed its reproductive prime. Instead of adjusting to a maintenance mode that would be more suitable for the non-reproductive phase of life, the TOR pathway continues to stimulate cellular processes actively. This ongoing activity can lead to “hyperfunction”, where cellular components and systems are overstimulated [13,14]. It is not clear, in that scenario, why the efficient processes of morphological homeostasis [16,17] do not orchestrate these active processes toward upkeep of functional anatomy (c.f., fig:introduction:motivation (B,C)).

Our approach is consistent with the antagonistic pleiotropy and hyperfunction models in that we also propose that fundamental morphogenetic processes are involved in aging. However, we do not focus on specific genes with early vs. late roles, but on the role that large-scale morphogenetic plasticity plays in establishing and maintaining bodies and the possible disorders of such cybernetic (goal-driven) processes. We and others have emphasized the indirect link [18,19,20] between genotype and phenotype: not just complexity and pleiotropy, but the fact that the morphogenetic layer has active homeodynamic capabilities that enable collections of cells to be aligned with respect to setpoints in the space of anatomical possibilities [16,21,22] (c.f., Figure 1 (B,C)). Living materials must navigate the space of possible organ and tissue arrangements as needed, starting from diverse initial states – single cells (egg-based development), a different body (in the case of metamorphosis), or a damaged body (in the case of regeneration). In all of these cases, drastic growth and remodeling must occur while failing or no longer needed sub-components must be replaced on-the-fly. This process has been described as a collective intelligence [23,24], because of its problem-solving capacities; like other advanced intelligences, it will be subject to specific failure modes whose origin is not organic disease or damage, but the informational dead end of completing one’s mission and not having a new setpoint toward which to strive.

We formulate, and simulate the consequences of, the hypothesis that the aging phenotype results from morphogenetic processes ceasing their targeted navigation of anatomical morphospace, after development is completed. But adulthood, and the rigors of daily life, are precisely when these processes are needed, in the context of stress, degeneration, and cell replacement, to maintain the Ship of Theseus that is the body. What will bind the cellular components of a body toward common paths in anatomical space after their goals have been reached and the body built? Tissue renewal will continue locally, for a time, but the thing that bound individual cells into a coherent “embryo” – a commitment to a specific navigation path in anatomical space – disappears at the global scale. Damage-based theories of aging are thus hypothesized to be secondary: they induce aging to the extent that they rob cells of clear and consistent anatomical and histological setpoints for which to aim. But the primary new proposal here is that complex goal-driven systems can experience a kind of intrinsic disorganization after they have met their goal. In the absence of strong effort toward targets in morphospace, the cellular collective begins to disband, giving way to systemic disease states such as degeneration and cancer – the inevitable outcome of a body that increasingly becomes a group of individual cells, not a unified higher-order morphogenetic agent [25].

1.2. Knowledge Gaps

Recent research suggests a shift in our understanding of aging, proposing that it results from changes in biological information processing rather than solely from cellular and molecular damage. This concept, linked to a malfunction in the “software of life”, focuses on how aging processes are governed by the epigenetic changes, either bioelectrical or chromatin modifications [5,26,27]. The differing dynamics, capabilities, and failure modes of the molecular layer and the physiological computational layers have long been studied in neuroscience, but are only now beginning to be related to areas such as developmental biology and cancer [28,29,30]. Novel computational frameworks are needed to begin to integrate these ideas in aging and longevity research. Although substantial advancements have been made in extending lifespan in simple model organisms through genetic, dietary, and pharmacological means [2,31,32], the fundamental mechanisms driving aging in humans remain elusive and widely debated. Current life extension strategies show limited effectiveness [33]. Much of the research has adopted a bottom-up approach, focusing on the biological hardware, the cellular components such as genes, molecular pathways, and, more recently, the epigenome. However, managing aging at such a micro-level may not be practical for complex structures like human anatomy [17,34,35]. This suggests that a broader, possibly more integrated approach may be necessary in order to fully understand and effectively intervene in the aging process.

Computational modeling has proven useful in the past to help characterize dynamics guiding complex morphogenesis [36,37,38,39,40,41,42,43,44]. However, to the best of our knowledge, there is currently no computational model in the literature explicating the aging process as a loss of goal-directedness, or directionality within the morphogenetic control system. The field is mostly dominated by experimental approaches and focused on the mathematical and computational modelling of a specific hallmark of aging [6], network approaches to aging or the impact of the environment on aging [45,46,47,48]. One prior study utilizes, similar to our model, a neural cellular automaton (see next section) to perform morphogenesis: with gradient-based deep learning, the self-regulatory interaction dynamics of distributed cellular agents have been optimized to grow a target morphology, fully self-orchestrated. However, without explicitly stabilizing morphogenesis by additionally training the cellular agents on regenerative tasks, anatomies kept growing uncontrollably after development. Thus, the study focused on an adversarial approach to stabilize self-orchestrated morphogenesis rather than studying aging and its causes in such a multi-cellular tissue per se [49].

Existing mathematical models have not yet captured the biological signaling dynamics in aging tissue or revealed the underlying information dynamics that drive aging and the progressive degradation of anatomy over time. Unlike previous approaches, we uniquely address this gap by developing an evolutionary framework that integrates both molecular and systemic aspects of aging within a multi-competency architecture. Our model demonstrates how cells, trained solely through evolution to execute developmental tasks, inherently lack a regenerative goal post-development. By simulating this dynamic, we reveal that aging arises from the absence of targeted navigation in anatomical morphospace after development, offering a novel perspective on the root cause of aging. Moreover, our approach uncovers the potential for dormant regenerative programs to be reactivated, enabling targeted restoration of anatomical structures (c.f., fig:introduction:motivation).

1.3. Biological Architecture of the Model

Morphogenesis is the result of multicellular collectives acting together to achieve a specific endpoint within the space of anatomical possibilities [19] (see fig:introduction:motivation). This process is driven by intercellular communication and information processing, where cells use different communication channels (biochemical, biomechanical or bioelectrical) to exchange information with their neighbors in a tissue [21,50] and implement a process of anatomical homeostasis where deviations from the species-specific target morphology are progressively minimized by active changes in growth and form. This enables cells to adaptively respond in a context-sensitive manner and regulate their behavior through local signaling protocols while collectively following a system-level agenda to grow and maintain an organism’ species-specifis anatomical form in a self-orchestrated way. Thus, the development of embryonic tissues, the growth and remodeling of adult organs, and even the suppression of cancer and aging can be understood as morphological computation [50]. In this sense, morphogenesis can be seen as an intermediate layer of physiological computation between genetic variation and a mature anatomical phenotype (see fig:introduction:motivation (C), significantly effecting the process of evolution [19,41].

Another important aspect of biology, so far neglected in the models of aging is multiscale competency [19,23,51] (see fig:introduction:motivation (B, C)). Biological systems implement an architecture in which context-sensitive problem-solving operates across various levels, from molecular networks to individual cells, tissues, organs, organisms, and even collective groups such as swarms. This structure allows biological entities to solve problems within specific domains and problem spaces relevant to each scale (metabolic, physiological, anatomical, and behavioral state spaces), contributing to the overall adaptability and resilience of life [21,50,52]. At the molecular and cellular level, biological units navigate metabolic spaces to remain within adaptive regions despite changing circumstances. At higher organizational levels, these units form tissues and organs, which tackle more complex and larger-scale physiological and morphological problems. For instance, cells can make autonomous decisions to maintain local balance while contributing to larger-scale biological processes like tissue repair or development and how higher levels of organization constrain and facilitate the behavior of their parts [19,53,54,55,56,57,58]. This enables organisms to adapt to changing conditions and maintain functionality even when faced with disruptions, as observed in regulative development and regeneration where cellular activity is managed toward specific anatomical outcomes despite injury, changes of cell number and size, etc. [22]. This is critical because resistance to both cancer and aging must involve upkeep of complex structures and a resistance to their progressive decline within the Ship of Theseus of the body [59].

A key advantage of this multiscale architecture is that each level is capable of addressing its own problem-spaces without needing constant guidance from higher organizational structures (local homeodynamic competency, see fig:introduction:motivation (C)). This decentralized problem-solving approach ensures that different biological processes, from cell division to behavioral adaptation, occur simultaneously and efficiently. Furthermore, this structure allows biological systems to be both robust and flexible, as failures at one level can often be compensated for by others. Additionally, the interplay between levels within this architecture enables the emergence of complex properties that cannot be explained by examining components in isolation. For example, the collective behavior of cells during processes like regeneration or morphogenesis illustrates how local interactions scale up to contribute to the organism’s overall form and function [23]. Similarly, defects of collective behavior can lead to developmental defects or aging. These interactions are not merely a sum of their parts; rather, they represent emergent capabilities that arise from the coordinated effort of biological subunits. The relative autonomy of layers in biology must surely impact on evolutionary and physiological constraints and drivers of aging, making it essential to study their properties with respect to the lifespan of individuals and evolutionary lineages.

This hierarchical but integrated organization is essential for the evolution of complex life forms. It promotes adaptability, allowing organisms to respond to environmental changes, evolve new features, and develop intricate systems that support survival and reproduction. Moreover, this multiscale competency architecture provides a framework for innovation in fields such as synthetic biology and regenerative medicine [24,34,60], where understanding and harnessing these principles could lead to breakthroughs in tissue engineering, organismal repair and aging [24,26]. This aspect of living forms – the critical layer of computational, problem-solving competency that lies between the genotype and the phenotype is an essential component of any theory of aging that integrates evolutionary and physiological dynamics to understand this complex phenomenon. Thus, we focused our computational model on the role of multiscale morphogenetic competency in the aging process.

1.4. Our Computational Model

Recently [41], we quantitatively investigated in silico the influence of varying levels of cellular competency of a multi-scale competency architecture [19,50] on the evolutionary process. We utilized tools from Artificial Life [61], a field of research focused on developing computational and cybernetic models that exhibit life-like behavior based on ideas derived from biology; prominent examples are Cellular Automata [62,63,64,65] such as Conway’s Game of Life [66]. Specifically, we deployed Neural Cellular Automata [67] (NCAs) as a model for evolving morphogenesis (see Figure 2 (A-C)): NCAs comprise a regular spatial grid of locally interacting cellular agents which maintain a numerical state vector that is regulated by an internal, cell-specific Artificial Neural Network (ANN). In that way, the cells of an NCA can be trained to perform a collective self-orchestrated system-level task, such as morphogenesis [39,40,41] or even collective navigation behavior [51,68], solely via successive perceptions of their local environment, i.e.,, reading neighboring cell state information, and computing corresponding actions to regulate their own states in accordance with a system-level agenda.

Thus, NCAs exhibit striking similarities with the genome- and collective intelligence-based multi-scale competency architectures of biological life [19,20,51,69]: the instructions of how to grow an organism from the bottom up is encoded in its genome, corresponding to the NCA’s (functional) ANN parameters. Moreover, through subsequent localized perception-action cycles of their ANN-based unicellular agents, NCAs model the intercellular communication- and intracellular information processing pathways of the adaptive and flexible internal decision-making mechanisms of biological cells: while biological cells regulate their gene expressions via gene regulatory networks, the cells of an NCA regulate their cell types via their internal ANN, and while biological cells utilize bio-electrical and molecular cues for intercellular signaling, cells in an NCA exchange numerical state information alike. In both cases, such communication pathways are localized between neighboring cells while the collective of cells may still follow a global agenda to self-orchestrate the process of morphogenesis and thereby give rise to the anatomical form of a system-level organism. Here, we use ANNs to model the memory, decision-making, and integrated information processes operating inside of cells [70,71,72,73,74].

Following our previous work [41], we deploy a swarm of virtual unicellular agents on the spatial grid of an NCA. Each agent hosts a decision-making machinery modeled by an ANN (notably, identical ANNs are deployed for each cell throughout a particular NCA instance), enabling the cellular agents to perceive neighboring cell states and independently propose cell-state regulatory actions, as illustrated in fig:introduction:simulation (A, B) and detailed below. Furthermore, we utilize evolutionary algorithms (EAs) [75] to optimize the parameters of NCAs, comprising structural and functional parts [41], such that the cellular agents evolve to self-assemble a particular spatial target pattern of predefined cell states in a fixed number of developmental steps, solely based on local communication rules, see fig:introduction:simulation (C).

Similar to ref. [41], we evolve NCAs to perform morphogenesis and morphostasis tasks, here specifically of a $(16\times 16)$ smiley-face-pattern reminiscent of the bioelectric craniofacial prepattern that establishes the rough outlines of the amphibian face during early development7 [76]. We deploy the evolved NCA solution for much longer times compared to the developmental stage at which the evolutionary selection occurs (c.f., fig:introduction:simulation (C, D)) and study the NCA’s long-term behavior with information-dynamic measures such as Active Information Storage (AIS)[77], Transfer Entropy (TE) [78], and Spatial Entropy (SE) . This parallels a biological an organism’s ability to maintain its physiological integrity via multiscale homeodynamic competencies for much longer than necessary for genetic reproduction. In turn, this may lead to an accumulation of a variety of entropically and self-regulatory favored morphological changes that are not directly fatal but potentially increasingly dysfunctional, changing the setpoint of the unicellular ’agents’ self-regulatory scheme and leading to misdifferentation (c.f., fig:introduction:simulation (D)).

Moreover, based on such virtual organisms, we here investigate in silico the information-dynamic implications of temporally affected competencies or information loss at the cellular level of a multi-cellular collective whose parts pursue a system-level agenda. To systematically interfere with the cellular competencies in our NCA model over time, we deploy four different strategies: We (a) increasing the misdifferentiation rate of cell-state expressions, (b) reduce the cellular agent’s ability to self-regulate their respective states by corrupting the reliability at which cellular actions are actually considered in the NCA’s cell state updates, (c) reduce the connectivity between neighboring cells reminiscent of closing an increasing number of gap junctions (GJs) to neighbors, and (d) accumulate genetic damage to the NCA’s functional parameters by randomly modifying the corresponding ANN parameters with Gaussian noise successively throughout the lifetime of our virtual organisms.

Thus, we relate degrading competencies and temporal information loss at the cellular level to aging at the system level. Finally, based on the hypothesis of the information-dynamical implications of aging caused by decaying cellular competencies, we demonstrate - in an silico minimal model - a reliable way of rejuvenation based on targeted interventions at the cellular- or organ level - of facial features.

We suggest that even in the absence of accumulated cellular defect of information processing and communication or genetic damage, a biological system left without any new goal in the morphospace will start degrading anatomically (and that this situation is allowed by evolution because it typically prioritizes development over long-term morphostasis, leaving organisms without any goal in the morphospace after development). Using an in silico model of homeostatic morphogenesis with a multiscale competency architecture and information dynamics analysis, we find: : (1) Absence of Long-Term Morphostasis: Aging emerges naturally after development due to the lack of an evolved regenerative goal, rather than being caused by specific detrimental properties of developmental programs (e.g., antagonistic pleiotropy or hyperfunction); (2) Acceleration Factors vs. Root Cause: Cellular misdifferenciation, reduced competency, communication failures, and genetic damage all accelerate aging but are not its primary cause; (3) Information Dynamics in Aging: Aging correlates with increased active information storage and transfer entropy, while spatial entropy measures distinguish two dynamics, (i() the loss of structure and (ii) morphological noise accumulation; (4) Dormant Regenerative Potential: Despite organ loss, spatial information persists in our cybernetic tissue. This indicates a memory of lost structures, which can be reactivated for organ restoration through targeted regenerative information; and (5) Optimized Regeneration Strategies: Restoration is most efficient when regenerative information includes differential patterns of affected cells and their neighboring tissue, highlighting strategies and protocols for rejuvenation. These findings provide a novel perspective on aging dynamics with significant implications for longevity research and regenerative medicine.

2. Material And Methods

2.1. Neural Cellular Automaton: A Multi-Agent Model for Morphogenesis, and Aging?

Cellular Automata (CAs) were initially introduced by von Neumann to study self-replicating machines [62]. Since then, they have become widely used as simple models for Artificial Life [61]. The core concept behind CAs revolves around maintaining a discrete spatial grid of cells, each individual cell i being assigned a binary, numerical, or even vector-valued state ${\mathbf{c}}_i\left(t_k\right)$ at every time step $t_k$.

The evolution of these cell states over time follows local update rules. Specifically, we have that ${\mathbf{c}}_i\left(t_{k+1}\right)=u\left({\mathcal{N}}_i\left(t_k\right)\right)$, where $u(·)$ is a function of the current state of cell i, i.e., ${\mathbf{c}}_i\left(t_k\right)\equiv {\mathbf{c}}_{i_0}\left(t_k\right)$, and of the states ${\mathbf{c}}_{i_{\nu }}\left(t_k\right)$ of its $i_{\nu }=i_1,\cdots ,i_N$ neighboring cells, which we collect in the matrix ${\mathcal{N}}_i\left(t_k\right)=\left({\mathbf{c}}_{i_0}\left(t_k\right),\cdots ,{\mathbf{c}}_{i_N}\left(t_k\right)\right)$.

Despite their typically simple and predefined update rules, CAs often exhibit complex dynamics (c.f., Conway’s Game of Life [66]) and have even been employed for universal computation tasks, as in Wolfram’s rule 110[63,64].

Neural Cellular Automata (NCAs) [67] can be seen as an extension to “traditional” CAs. In this context, the local update rule is replaced with a more flexible ANN, $u(·)\to f_{\theta }(·)$, where $\theta$ denotes the set of trainable parameters of the ANN, $f_{\theta }(·)$ (see, e.g., Appendix A in Ref. [41]). By employing Machine Learning (ML) techniques, NCAs have been utilized to perform tasks such as self-orchestrated pattern formation [39], the simultaneous co-evolution of a rigid robot’s morphology and controller [51], and have been proposed as a promising candidate for robust, decentralized controllers of autonomous drug delivery systems [68].

An NCA is essentially a grid of interconnected cells, each equipped with an identical ANN that is capable of perceiving the numerical states of its immediate cellular neighbors, ${\mathcal{N}}_i\left(t_k\right)$, and proposing actions, ${\mathbf{a}}_i\left(t_k\right)=f_{\theta }\left({\mathcal{N}}_i\left(t_k\right)\right)$, to regulate its own cell state, following the equation

$${\mathbf{c}}_i\left(t_{k+1}\right)={\mathbf{c}}_i\left(t_k\right)+{\mathbf{a}}_i\left(t_k\right)+{\xi }_{\mathbf{c}},$$

(1)

where we also account for possible imperfections in the cell-specific state updates via a Gaussian noise term of amplitude ${\xi }_{\mathbf{c}}$.

Thus, the cellular agents of an NCA perceive the numerical states of their immediate, respective neighborhoods, ${\mathcal{N}}_i\left(t_k\right)$, at any given time step, $t_k$, to update their own states (c.f., fig:introduction:simulation (A, B)). Importantly, we do not include any explicit signature of aging in our model, i.e.,, no cell-internal marker or explicit input to the cellular agents other than their self-regulated cellular states that would contain temporal information.

However, being immersed in a dynamical multi-cellular environment where every part has its own agenda, the cellular agents’ state updates can also be utilized for active communicating with their neighbors, following a policy $\pi \left({\mathcal{N}}_i\left(t_k\right)\right)\approx f_{\theta }\left({\mathcal{N}}_i\left(t_k\right)\right)$. From the perspective of Reinforcement Learning [80], an NCA can thus be considered a trainable multi-agent system that needs to utilize local communication rules to achieve a target system-level outcome.

So far, our approach is agnostic to the particular ANN architecture of the update function, $f_{\theta }(·)$. Closely following Ref. [41], we here briefly describe the genuine ANN architecture deployed in our cellular agents, which is schematically depicted in fig:introduction:simulation (B): Inspired by Ref. [81], we partition an NCA’s cell’s ANN into (i) a sensory part, ${\epsilon }_{i_{\nu }}=f_{\theta }^{\left(s\right)}\left({\mathbf{c}}_{i_{\nu }}\left(t_k\right)\right)\in {\mathbb{R}}^s$, preprocessing each neighborhood cell state separately into a respective sensor embedding, which we collect in the matrix $\mathcal{E}\left({\mathcal{N}}_i\left(t_k\right)\right)=\left({\epsilon }_{i_0}\left(t_k\right),\cdots ,{\epsilon }_{i_N}\left(t_k\right)\right)$. These sensor embeddings are then (ii) averaged across all neighbor embeddings into a single context vector ${\mathbf{s}}_i\left(t_k\right)={\displaystyle \frac{1}{N+1}}{\sum }_{\nu =0}^N{\epsilon }_{i_{\nu }}\left(t_k\right)\in {\mathbb{R}}^s$. A subsequent (iii) controller ANN, $f_{\theta }^{\left(c\right)}(·)$, potentially with recurrent feedback connections, eventually outputs the cell’s action, ${\mathbf{a}}_i\left(t_k\right)=f_{\theta }\left({\mathcal{N}}_i\left(t_k\right)\right)=f_{\theta }^{\left(c\right)}\left({\mathbf{s}}_i\left(t_k\right)\right)$ based on the cell-specific context vector ${\mathbf{s}}_i\left(t_k\right)$. Notably, for the controller module we here utilize so-termed RGRN architectures [41] that are inspired by a combination of Recurrent ANNs (RNNs) [82] and Gene Regulatory Networks [83] (see Appendix A in Ref. [41] for details).

Analogous to Ref. [41], we model morphogenesis by employing NCAs on a two-dimensional $N_x\times N_y$ square grid with the objective that all cells of the grid assume a cell-specific predefined target cell type, ${\widehat{g}}_i$, after a fixed number of $t_{\mathrm{D}}$ developmental time steps, starting from an initial cell state configuration ${\mathbf{c}}_i\left(0\right)$. We assign the first $N_{\mathrm{G}}$ elements, ${\mathbf{g}}_i\left(t_k\right)=\left(c_{i1}\left(t_k\right),\cdots ,c_{i{N}_{\mathrm{G}}}\left(t_k\right)\right)$, of an NCA’s $N_C$-dimensional cell state ${\mathbf{c}}_i\left(t_k\right)=\left(c_{i1}\left(t_k\right),\cdots ,c_{i{N}_{\mathrm{C}}}\left(t_k\right)\right)$ as indicators for expressing one of $1,\cdots ,N_G$ discrete cell types, while the remaining $N_H=(N_C-N_G)$ elements of the cell state represent its hidden states, ${\mathbf{h}}_i\left(t_k\right)=\left(c_{i(N_{\mathrm{C}}-N_{\mathrm{H}})}\left(t_k\right),\cdots ,c_{i{N}_{\mathrm{C}}}\left(t_k\right)\right)$, that can be utilized by the NCA for intercellular communication. A cell’s “type”, $g_i\left(t_k\right)$, is now defined as the index, $1,\cdots ,N_G$, of the maximum element of its indicator vector ${\mathbf{g}}_i\left(t_k\right)$:

$$g_i\left(t_k\right)=\underset{\mathbf{g}\in {\mathbb{R}}^{N_G}}{argmax}\left({\mathbf{g}}_i\left(t_k\right)\right).$$

(2)

For an NCA to assemble a predefined target pattern of $N_{\mathrm{j}}=N_x\times N_y$ target cell types $\widehat{\mathbf{g}}=\{{\widehat{g}}_1,\cdots ,{\widehat{g}}_{N_{\mathrm{j}}}\}$, we thus need to find a suitable set of NCA parameters that minimizes the deviation of all actual cell types $g_i\left(t_D\right)$ and the desired ones ${\widehat{g}}_i$ after $t_D$ developmental time steps. Below, we introduce an evolutionary algorithms to evolve suitable sets of NCA parameters that maximize a fitness score based on comparing the “final” cell types of the NCA, $g_i\left(t_D\right)$, after the developmental stage to the predefined target cell types ${\widehat{g}}_i$. In this context, the NCA parameters thus correspond to the virtual organism’s genotype, while the grid of final cell types represent the system-level phenotype that is seen by the selection mechanism of the evolutionary process (c.f., fig:introduction:simulation (C)).

2.2. Neuroevolution of NCAs: An Evolutionary Algorithm Approach to Morphogenesis

Evolutionary Algorithms (EAs) are heuristic optimization techniques designed to maintain and optimize a set, or population, $\mathcal{X}=\{{\mathbf{x}}_1,\cdots ,{\mathbf{x}}_{N_{\mathrm{P}}}\}$, comprising parameters or individuals, ${\mathbf{x}}_j\in {\mathbb{R}}^{\mathrm{X}}$, over successive generations in order to maximize an objective function or fitness score, $r\left({\mathbf{x}}_j\right):{\mathbb{R}}^{\mathrm{X}}\to \mathbb{R}$. Drawing inspiration from the principles of natural selection and biological life’s DNA-based reproduction mechanisms, EAs (i) predominantly select high-fitness individuals for reproduction, and utilize (ii) crossover and (iii) mutation operations to generate new offspring by (ii) merging the genomic material of two high-quality individuals from the current population, ${\mathbf{x}}_o={\mathbf{x}}_j⨁{\mathbf{x}}_k$, and (iii) occasionally mutating offspring genomes, ${\mathbf{x}}_o\to {\mathbf{x}}_o+{\xi }_{\mathbf{x}}$ by adding typically Gaussian noise to the parameters; the ⨁ symbol signifies a genuine merging operation of two genomes, which may vary according to the specific EA implementation. Thus, populations of individuals are directed towards high-fitness regions within the parameter space ${\mathbb{R}}^{\mathrm{X}}$, typically across numerous generations of successive selection and reproduction cycles (i)-(iii) [41].

Closely following our previous work [41], and in contrast to many “traditional” EA use-cases, we here explicitly distinguish between the genotypic parameters, ${\mathbf{x}}_j$, and the corresponding phenotypic realizations, ${\mathbf{p}}_j$. More specifically, we utilize NCAs (c.f., sec:methods:nca) to model the biologically inspired developmental layer [19] in between genotypes and phenotypes, ${\mathbf{x}}_j\to {\mathbf{p}}_j$, akin to the process of morphogenesis: Here, a genotype, ${\mathbf{x}}_j$, corresponds to a particular realization j of an NCA’s parameters, comprising the set of initial cell states ${\mathbf{x}}_j^{\left(\mathrm{S}\right)}={\left\{{\mathbf{c}}_{i=(1,\cdots ,N_j)}\left(0\right)\right\}}_j$ and the corresponding ANN parameters ${\mathbf{x}}_j^{\left(\mathrm{F}\right)}={\theta }_j$, such that

$${\mathbf{x}}_j={\mathbf{x}}_j^{\left(\mathrm{S}\right)}\cup {\mathbf{x}}_j^{\left(\mathrm{F}\right)}=\left({\left\{{\mathbf{c}}_i\left(0\right)\right\}}_j,{\theta }_j\right),$$

(3)

where we explicitly distinguish between structural (S) and functional (F) genes. We then employ eq:nca-update for $t_{\mathrm{D}}$ developmental steps to successively transform the NCA’s cell states from an initial state into a “final” set of cell types ${\left\{g_i\left(t_{\mathrm{D}}\right)\right\}}_j$ on the NCA’s grid. This mature phenotype,

$${\mathbf{p}}_j={\left\{g_i\left(t_{\mathrm{D}}\right)\right\}}_j,$$

(4)

represents a two-dimensional tissue of cells, and is the input to the fitness score based on which the EA selects: As illustrated in fig:introduction:simulation (C), we here optimize for a fitness score in the phenotype- rather than genotype space, $r\left({\mathbf{x}}_j\right)\to r\left({\mathbf{p}}_j\right)$, while still performing genetic recombination and mutation in the genotype space.

Let’s recall: our goal is to achieve morphogenesis of a two-dimensional spatial tissue of cell types, ${\mathbf{p}}_j$, that optimally resembles a predefined target pattern, $\{{\widehat{g}}_1,\cdots ,{\widehat{g}}_{N_{\mathrm{j}}}\}$, of a total of $N_j$ cells on an $N_x\times N_y$ square grid of an NCA. Moreover, we aim to derive a virtual organism with this capability by modeling a biologically inspired evolutionary process with EAs to evolve a suitable set of NCA parameters (c.f., fig:introduction:simulation (C)). Thus, we introduce our phenotype-based fitness score, $r\left({\mathbf{p}}_j\right)$, as [41]

$$r\left({\mathbf{p}}_j\right)=(2n_j^{\left(\mathrm{G}\right)}-N_j)+r_T{n}_j^{\left(\mathrm{T}\right)}-r_S{n}_j^{\left(\mathrm{S}\right)},$$

(5)

where $n_j^{\left(\mathrm{G}\right)}$ counts the number of correctly assumed cell types after $t_{\mathrm{D}}$ developmental steps (where $g_i\left(t_{\mathrm{D}}\right)={\widehat{g}}_i$), $n_j^{\left(\mathrm{T}\right)}$ counts the number of developmental steps at which the target pattern is realized entirely (whenever $g_i(t_k\le t_{\mathrm{D}})={\widehat{g}}_i$ for all i), and $n_j^{\left(\mathrm{S}\right)}$ counts pairs of successive time steps, $(t_s,t_{s+1})\le t_{\mathrm{D}}$, without cell type updates in the entire tissue (where $g_i\left(t_{s+1}\right)=g_i\left(t_s\right)$ for all i).

Thus, we design an evolutionary process that primarily selects for NCAs that assume the correct cell type pattern after (precisely) $t_D$ developmental steps and reinforces NCAs (by a factor of $r_T$ in the fitness score) that prematurely achieve and maintain the target pattern. To avoid developmental stagnation and increase the EA’s performance, we specifically discount the fitness score of NCAs with static behavior by $r_S$. A fitness score of $N_j=N_x\times N_y$ indicates that the problem is solved.

Although the above framework [41] can be utilized in conjunction with any black-box evolutionary or genetic algorithm, we found the well-established Covariance Matrix Adaptation Evolutionary Strategy (CMA-ES) [75] to be a viable choice for our purposes.

2.3. Information-Theoretic Analysis: Active Information Storage, Transfer Entropy and Spatial Entropy on the NCA

Information theory [84] offers valuable insights into the dynamics of complex systems. In order to understand the NCA’s information dynamics during aging, we used three key information-theoretic metrics to examine the information dynamics: active information storage [77], transfer entropy [78], and the spatial entropy.

Active Information Storage (AIS) quantifies the amount of information from an agent’s past that is pertinent to predicting its future state. Specifically, AIS refers to the portion of stored information currently utilized for determining the agent’s next state [77]. Mathematically, the AIS of an agent Q is expressed as the local mutual information between its semi-infinite past $q_n^{\left(k\right)}$ as $k\to \infty$ and its next state $q_{n+1}$ at time step $n+1$:

$$a_Q(n+1)=\underset{k\to \infty }{lim}{log}_2{\displaystyle \frac{p(q_n^{\left(k\right)},q_{n+1})}{p\left(q_n^{\left(k\right)}\right)p\left(q_{n+1}\right)}}$$

(6)

Here, $a_Q(n,k)$ is an approximation with history length k. The time-averaged value, weighted by the distribution of $(q_n^{\left(k\right)},q_{n+1})$, is denoted as $A_Q\left(k\right)=〈a_Q(n,k)〉$. In our study, we calculated the local AIS across the states of the cells.

Transfer Entropy (TE) measures the information transferred from a source agent to a destination agent that is not contained in the past of the destination agent. We employed the local TE concept introduced by Lizier [77]. The local TE from a source agent Z to a destination agent Q is defined as the local mutual information between the previous state of the source $z_n$ and the next state of the destination agent $q_{n+1}$, conditioned on the semi-infinite past of the destination $q_n^{\left(k\right)}$ (as $k\to \infty$):

$$t_{Z\to Q}(n+1)=\underset{k\to \infty }{lim}{log}_2{\displaystyle \frac{p\left(q_{n+1}\right|q_n^{\left(k\right)},z_n)}{p\left(q_{n+1}\right|q_n^{\left(k\right)})}}$$

(7)

The transfer entropy $T_Q(n,k)$ is the local TE averaged over time, denoted as $T_Q\left(k\right)=〈t_Q(n,k)〉$, where $t_Q(n,k)$ is an approximation with history length k. Unlike mutual information, which measures static correlation, transfer entropy captures the dynamic, directional flow of information within the agent network.

Spatial Entropy (SE) measures the randomness or disorder within a spatial distribution of states. It provides insight into the complexity of spatial patterns in a system, such as an NCA with multiple states. In our study, we compute SE by evaluating the entropy of the distribution of cell states across the grid at each time step. Formally, the spatial entropy H at time step n is defined as:

$$H\left(n\right)=-\sum _{s\in S}p\left(s\right){log}_{2}p\left(s\right)$$

(8)

where S is the set of possible states, and $p\left(s\right)$ is the probability of state s occurring in the grid at time step n. This measure provides an indication of the level of unpredictability or complexity in the spatial arrangement of cell states. In this article, we calculate the spatial entropy for the states of the NCA across different time steps.

3. System

Following our previous study [41], we utilize EAs to evolve NCAs to perform morphogenesis tasks. Each cell in our NCA implementation is equipped with a gene regulatory network-inspired recurrent ANN, conceptually substituting the intricate internal decision-making competencies of biological cells (see fig:introduction:simulation (B) for a visualization and Appendix A in [41] for numerical details). In other words, each cell in the NCA is equipped with an information processing pipeline that transforms sensory input of its local environment in the NCA’s grid into control instructions to regularize its own cell type. If done in a coordinated way, these local individual cell-state updates guide the formation of a mature system-level phenotype - in our case of a target tissue - over successive developmental steps.

More specifically, we utilize $N_x\times N_y=16\times 16$ cells on an NCA with fixed boundary conditions (see sec:methods:nca and Appendix C in Ref. [41]) and train it analogously to [41] via evolutionary algorithms (see sec:methods:ea) to perform morphogenesis - i.e.,, self-orchestrated tissue growth - of a $(16\times 16)$ smiley-face pattern in a predefined number of $t_{\mathrm{D}}=35$ developmental steps (see fig:introduction:simulation (C) for an illustration of the evolutionary process).

The target tissue consists of $N_{\mathrm{G}}=3$ distinct cell types: one cell type for the background cells, one for the face cells, and a single cell type for the facial organs, i.e.,., for the eyes and the mouth (c.f., blue-, magenta-, and white-colored cells in fig:introduction:simulation (A)). Moreover, we allow the NCA to utilize $N_{\mathrm{H}}=1$ hidden channel of its $N_{\mathrm{C}}=4$ cell state channels for unconstrained intercellular communication. In cell-state updates, given by eq:nca-update, we utilize a noise level of ${\xi }_{\mathbf{c}}=0.125$, limit the numerical value of the proposed action ${\mathbf{a}}_i\left(t_k\right)$ to $l_{\mathbf{a}}=[-1,1]$, and clip the cell-states post-update ${\mathbf{c}}_i\left(t_{k+1}\right)$ to $l_{\mathrm{c}}=[-3,3]$ (see Appendix C and Ref. [41] for more details on our NCA implementation).

Our NCAs transform an embryonic genotype, ${\mathbf{x}}_j$, into a mature phenotype, ${\mathbf{p}}_j$, and we here quantify an NCA’s performance by evaluating the corresponding fitness score, $r\left({\mathbf{p}}_j\right)$, defined by eq:fitness with respect to the $(16\times 16)$ smiley-face pattern over the NCA’s entire life-span of $t_{\mathrm{D}}$ developmental steps. As illustrated in fig:system:morphogenesis, a successfully evolved NCA maximizes the fitness score to $r_{\max }=256$, exhibiting $16\times 16$ correctly expressed cell types after $t_k⪆20$ time steps, then maintaining the target pattern until $t_k=t_{\mathrm{D}}$ (see sec:methods:ea for details on the evolutionary process).

Notably, the cellular agents do not perceive any information about the current nor the final fitness score of their corresponding NCA. The fitness score is only utilized by the EA to optimize the NCA’s parameters over evolutionary time-scales. Moreover, there is no contribution neither to the fitness score nor to the NCA’s ANN sensory inputs that would provide the system with explicit temporal information indicating the “age” of any part of the virtual organism, or even a notion of “time” in general. Thus, the self-orchestrated pattern formation of the here - and previously [41] - investigated NCAs is fully driven via emergent, i.e.,, evolved, intercellular communication rules and intracellular decision-making.

4. Computational Results

4.1. Aging as a Loss of Goal-Directedness: Organism Learned Development During Evolution, Not to Maintain Anatomical Homeostasis After Development

We aimed to assess whether the self-regulatory behavior of biological systems that drive morphogenesis induces aging-like phenomena at long time-frames. To this end, we use NCAs evolved in silico to model the process of morphogenesis and then deploy their self-regulatory dynamics for much longer times compared to the reproduction stage during evolution (i.e., when selected for reproduction). We monitor and analyze the NCAs morphological integrity throughout their lifetimes and compare such long-term morphological trajectories with biological aging. The evolutionary process described in sec:methods:ea,sec:system explicitly selects for optimal fitness scores after $t_{\mathrm{D}}$ developmental steps, which corresponds to the fertility age of real-world organisms. However, this implies that the behavior of an NCA with a potentially much longer lifetime of $t_{\mathrm{A}}>t_{\mathrm{D}}$, and particularly $t_{\mathrm{A}}\gg t_{\mathrm{D}}$, is not subject to selection and may fundamentally be considered ill-defined or ambiguous.

To demonstrate this, we deploy the NCA solution that has been evolved to self-assemble a $(16\times 16)$ smiley-face pattern during a morphogenesis phase of $t_{\mathrm{D}}=35$ developmental steps (c.f., fig:system:morphogenesis) for a much longer lifetime of up to $t_{\mathrm{A}}=1500$ time steps, notably without any further optimization considering the extended lifetime. From the results presented in fig:result:pleiotropy, we learn that after an initial rapid rise in the fitness score during the developmental phase, the NCA’s long-term fitness score gradually declines, primarily due to the loss of different internal organs such as a single or both eyes and, occasionally, the mouth. Moreover, this is not only a feature of the particular NCA demonstrated here, but is a systematic phenomenon of the presented framework.

Thus, we argue that the biologically-ubiquitous self-regulatory behavior of a substrate’s agential parts may induce system-level dynamics that resemble the process of biological aging if applied on timescales that exceed development. The organisms learned during evolution with development as the primary goal; once they reached that goal, we observed that they don’t regenerate by themselves and slowly show signs of morphological deterioration, e.g., aging. In this sense, aging can be seen primarily as a loss of goal-directedness, the goal of development being different from the goal of maintenance of the anatomy over timescales exceeding development. This view on aging is related to programmatic theories of aging like antagonistic pleiotropy or hyperfunction theory [10,14,15] that state that what is good for development could ultimately cause aging. In our framework, the lack of goal-directedness after development may cause some informational antagonistic pleiotropy (or hyperfunction) as cells lose their developmental policy after reaching the target morphology but it is broader in scope as giving a new goal to the cellular collective like regeneration and therefore by reactivating developmental pathways may cure aging (see regenerative experiments below).

4.2. Impact of Defects of Cellular Information Processing at Different Levels on the Rate of Aging in a Multi-Scale Competency Architecture

To determine how the perception-action cycles of the unicellular agents of a collective multi-cellular organism influence, at least in part, the rate of morphological decline, we here investigate the effects of an NCA’s system parameters on its long-term morphological integrity. Based on the NCA system defined in section III, we identify four key mechanisms involved in the corresponding uni-cellular perception-action cycles affecting morphological integrity, namely cell-specific differentiation reliablility (a), reliable state updates (b), reliable communication between neighbors (c), and persistent genetic encoding (d). Our framework allows us to selectively control these corresponding processes (a-d) and explicitly corrupt them independently over the lifetime of the virtual organism. Thereby we demonstrate the severe effects of these processes on the rate of morphological aging.

In the following sections, we therefore study the factors accelerating aging, not the cause of aging per se that we see primarily as a loss of goal-directedness. We introduce these known aging mechanisms (i-iv), respectively, based on and deployed to the biologically motivated unicellular competencies in our NCA system. We then demonstrate and discuss the corresponding implications on enhanced rates of morphological aging in long-term simulations of such virtual organisms. We illustrate the respective results in a compact form in fig:results:schedules(A-D), and continue below with an information-theoretic analysis of these processes.

4.2.1. Cellular Differentiation

One critical parameter in our system is the noise level ${\xi }_{\mathbf{c}}$ utilized to corrupt the NCA’s update step given by eq:nca-update. While a moderate noise level can help an evolutionary process to identify robust NCA parameters that solve a specific morphogenesis task, chosen too low, the evolutionary process will in our case most likely identify a direct encoding of the target pattern, largely neglecting cellular competencies in the pattern formation task [41]. In turn, a too-large noise level will not only prevent the evolutionary process from finding a good solution to the morphogenesis task but, in general, drastically limit the cell specification capacities of the NCA’s cell states. Therefore, noisy cell state updates in our system could be seen as the accumulation of misdifferentiated cells, ultimately corrupting the state perceptions of neighboring cells on the NCA creating a default of cellular communication (a hallmark of aging) similarly to [6,85]. This mechanism accelerates aging but is not a cause per se.

We argue that an increasing noise level ${\xi }_{\mathbf{c}}\to {\xi }_{\mathbf{c}}\left(t_k\right)$ over the lifetime of such an NCA-akin organism successively affects its morphological integrity. To demonstrate this, we thus utilize the pre-evolved NCA from above (c.f., fig:system:morphogenesis,fig:result:pleiotropy) and deploy it for an extended lifetime of $t_{\mathrm{A}}=1000$ timesteps while successively increasing the noise level linearly from ${\xi }_{\mathbf{c}}(t_k=0)$ to ${\xi }_{\mathbf{c}}\left(t_A\right)=1$. The results depicted in fig:results:schedules (A) show that while the correspondingly affected virtual organisms initially follow their unaffected cousins (c.f., fig:system:morphogenesis,fig:result:pleiotropy and black-dashed line in fig:results:schedules (A)), after a particular duration (specific to the particular NCA solution and noise schedule ${\xi }_{\mathbf{c}}\left(t_k\right)$) the fitness score rapidly drops at $t_k\approx 380-400$, emphasizing an enhanced rate of degradation of the target tissue compared to the baseline with constant noise ${\xi }_{\mathbf{c}}=\begin{array}{c}\hfill named-contentcolortypergbcontent-type0,0,00.125\end{array}$.

Noise influences the dynamical cell states in our NCA experiments. These states are the inputs and outputs of the virtual tissues’ ANN augmented cells, and represent dynamic variables of the self-regulatory NCA dynamics, with attractor states or instabilities. We thus interpret the effect discussed here as cellular misdifferentation-enhanced morphological aging. Understanding such “run-away” arguments in the “biological software” of cell state expressions, or instabilities in the context of competent tissues seems central to understanding aging.

4.2.2. Cellular Competency

Next, we directly manipulate the cellular decision-making competency - the reliability with which uni-cellular agents can regulate their own cell states within a tissue (e.g., regulative development and regenerative tissue allostasis). Here, we investigate whether degrading unicellular competencies throughout the lifetime of an organism enhances morphological aging: while young, healthy cells are highly reliable in their (re)actions to external perturbations, especially during development, the reliability of such cell decisions might be strongly affected over the lifetime of an organism, inducing morphological decline. To this end, and following 6, we define the decision-making as a probability $P_{\mathrm{D}}$ at which proposed actions of unicellular agents are actually forwarded to the environment at a given time step. More specifically, at every instance in time $t_k$ we draw a uniform random number $r_i\approx \mathrm{rand}\in [0,1]$ for every cell $i=1,\cdots ,N$ in the tissue. Only if $r\le P_{\mathrm{D}}$ we consider the cell’s action $a_i\left(t_k\right)$ in the NCA’s update step, described by eq:nca-update, and omit it otherwise by setting $a_i\left(t_k\right)=0$. In that way, we can scale the average reliability of decision-making of all cells in the tissue, i.e.,, while $P_{\mathrm{D}}=1$ corresponds to fully reliable decision-making of every cell at every time step, $P_{\mathrm{D}}=0$ parameterizes fully passive cells which can’t make decisions at all.

Here, we utilize $P_{\mathrm{D}}\to P_{\mathrm{D}}\left(t_k\right)$ as decreasing cellular competency by transforming $P_{\mathrm{D}}(t_k=0)=1$ to $P_{\mathrm{D}}\left(t_{\mathrm{A}}\right)=0$ in $t_{\mathrm{A}}=1000$ time steps following the functional form $P_{\mathrm{D}}\left(t_k\right)={sin}^4\left({\displaystyle \frac{\pi }{2}}(1-{\displaystyle \frac{t_k}{t_{\mathrm{A}}}})\right)$ (see fig:results:schedules (B)). This increasingly corrupts cellular competencies by successively decreasing the reliability at which every unicellular agent in the grid can regularize, and thus “correct” their states against noise or faulty cell state configurations.

The results of fig:results:schedules (B) reveal that the fitness score of NCAs whose cellular competencies are gradually affected remains similar to the baseline fitness scores of the unperturbed case (c.f., fig:system:morphogenesis,fig:result:pleiotropy). But eventually, at $t_k\approx 800$, the small but finite noise of ${\xi }_{\mathbf{c}}=\begin{array}{c}\hfill named-contentcolortypergbcontent-type0,0,00.125\end{array}$ used in our simulations leads to accumulated damage in the NCA’s states that is catastrophic for its morphological integrity, thus destroying the target morphology via diffusive processes.

Thus, decreasing cellular decision-making competencies in a noisy collective organism necessarily lead to morphological aging effects, caused by unreliability in the “biological cell-state reconfiguration software”. However, the rate at which these decreasing cellular competencies lead to significant aging effects is an intricate interplay between the actual noise in the system and the robustness of the unicellular agents to persist against accumulating noise in their respective state configurations.

4.2.3. Cell-Cell Communication

Another necessary feature of systems with decentralized collective decision-making is reliable communication between agents. Similar to biological tissue, the cells of an NCA are integrated into a lattice of functionally homogeneous cells that only diverge in their respective cell state history. However, their neighbor coordination typically remains fixed [39,41,67]. Here, we specifically manipulate the connectivity between the unicellular agents in the NCA’s grid by blocking increasingly many communication channels between neighboring cells as time proceeds. This approach is biologically motivated by cells closing their gap junctions (GJs) to their neighbors, which is known in vivo to shift cells from cooperatively working on large morphogenetic setpoints (organ maintenance) toward individual goals appropriate to unicellular organisms such as proliferation and migration [50,86,87,88]. It is also related to altered cell-cell communication found during aging [6]. Indeed, aging involves deficiencies in neural, neuroendocrine, and hormonal signaling pathways [89] as well as the alterations in the bidirectional communication between human genome and microbiome. Chronic inflammation is also part of this altered cellular communication [6].

We achieve this technically by omitting the corresponding input state ${\mathbf{c}}_{ij}=\mathbf{0}$ of a blocked connection, or GJ from cell j to cell i (see fig:introduction:simulation (B)). More specifically, we utilize a probability $P_{\mathrm{GJ}}\left(t_k\right)$ which defines the fraction of GJs that should remain open for each cell in the tissue at a given time $t_k$. In turn, if for a particular cell i, the fraction of connections to its neighbors exceeds $P_{\mathrm{GJ}}\left(t_k\right)$ at time $t_k$, we henceforth randomly block another one of its $i_{\nu }=i_1,\cdots ,i_N$ still connected GJs (once closed, GJs remain closed in our case); see sec:methods:nca and [41] for more details on our definition of cellular neighborhoods.

By linearly transforming $P_{\mathrm{GJ}}(t_k=0)=1$ to $P_{\mathrm{GJ}}\left(t_{\mathrm{A}}\right)=0$, we successively isolate the unicellular agents from their host tissue by blocking an increasing number of their GJs until, eventually at time $t_{\mathrm{A}}=\begin{array}{c}\hfill named-contentcolortypergbcontent-type0,0,01000\end{array}$, all cells are effectively isolated from their surroundings. In fig:results:schedules (C) we demonstrate the effects of this approach: Again, the self-orchestrated morphogenetic program can withstand the effects of increasingly corrupted competencies for a surprisingly long time. However, at $t_k\approx 500$, corresponding to $P_{\mathrm{GJ}}=0.5$ and thus half the GJs of every cell being closed, the effects of the removed communication channels can be observed by a resulting rapid drop in the phenotypical fitness scores caused by a rapid degradation of the target pattern.

Thus, decreasing the ability of cells to communicate with their neighbors leads to enhanced effects of morphological aging if a threshold of GJs is closed.

4.2.4. Accumulation of Genetic Damage

Another potential source for self-induced aging is accumulating damage to the cells’ genetic material via successive mutations of the functional genome ${\mathbf{x}}^{\left(\mathrm{F}\right)}\to {\mathbf{x}}^{\left(\mathrm{F}\right)}\left(t_k\right)$ over the lifetime of the virtual organism: To this end, we utilize the pre-evolved NCA solution discussed above (c.f., fig:system:morphogenesis,fig:result:pleiotropy), but at every time step $t_k>0$, we successively manipulate the cells’ ANN parameters throughout their lifetime. This directly introduces accumulating genetic mutations to the self-regulatory mechanisms of the NCA, causes potential dysfunctional decision-making behavior at the uni-cellular level, and correspondingly leads to enhanced morphological decline. More specifically, we add zero-centered normal distributed noise $\mathcal{N}(0,{\sigma }_x)$ of standard deviation ${\sigma }_x=\begin{array}{c}\hfill named-contentcolortypergbcontent-type0,0,00.01\end{array}$ to the functional genes of the NCA, ${\mathbf{x}}^{\left(\mathrm{F}\right)}\left(t_{k+1}\right)={\mathbf{x}}^{\left(\mathrm{F}\right)}\left(t_k\right)+{\delta }_{\mathbf{x}}$, where ${\delta }_{\mathbf{x}}\approx \mathcal{N}(0,{\sigma }_x)$. The mutated functional genes ${\mathbf{x}}^{\left(\mathrm{F}\right)}\left(t_k\right)$ are then mapped to the ANN parameters $\theta \to \theta \left(t_k\right)$ of the unicellular agents’ controllers at every time step $t_k$, thus directly affecting the intracellular decision-making machinery of the NCA.

From the results depicted in fig:results:schedules (D) we can learn several things: First, the evolved NCA solutions seem incredibly brittle against accumulating genetic damage in the ANN parameters, which is reflected by the fact that most of the deployed simulations (with statistically independent mutations) retain high-fitness target patterns for the entire simulated lifetime of the NCA of $t_{\mathrm{A}}=1000$ timesteps. Second, and in contrast to sec:results:induced-degeneracy:noise,sec:results:induced-degeneracy:competency,sec:results:induced-degeneracy:cell-communication this leads to large qualitative deviations in individual behavior of an ensemble of statistical independently mutated NCA solutions: some solutions die off early and obtain vanishing fitness scores, others degrade rather continuously with a spread-out fitness score of $\approx 100$, and many solutions still maintain their integrity although accumulating significant genetic damage of ${\sigma }_{\max }^2=t_{\mathrm{A}}{\sigma }_{\mathrm{x}}^2=0.1$ in total.

Thus third, while many mutations appear neutral in our framework, critical mutations to specific genes might induce detrimental modifications to the behavior of the unicellular agents (if they occur at the “correct” moment), typically causing a rapid degradation of the target pattern of the affected NCA as reflected by significant drops in the fitness scores of the corresponding virtual tissue. Notably, due to technical reasons, we mutate the functional genome of every agent homogeneously, thus retaining a functionally homogeneous distributed decision-making in the NCA, which nevertheless deviates functionally from the baseline at $t_k=0$. Potentially intriguing effects based on heterogeneous mutations, and thus heterogeneous agents in the NCA are out of the scope of this work, and we leave a more thorough investigation to a further contribution.

4.3. The Acceleration of Aging Is Linked to Increase in AIS and TE, While Spatial Entropy Revealed Two Different Kind of Aging: Loss of Structure and Proliferation, and Accumulation of Morphological Noise

In order to understand the information dynamics of aging and the cellular information processing in a competent tissue, we applied information-theoretic measures such as (i) Active Information Storage (AIS), and Transfer Entropy (TE) and the Spatial Entropy (SE), e.g., the entropy applied to the states of the organism. (see Figure 6).

We observed in the default behavior case, without any add of the noise or induced loss of competency that both the mean AIS and TE are decreasing during development and aging. For some organisms, we can see that TE can increase during aging, indicating episodes of regeneration in the sense that the morphology is impaired and the organisms try to recover. As we can observe it, the spatial entropy is slowly decreasing during aging, and the anatomy is slowly impaired, with the loss of organs like eyes (see Figure 6).

Interestingly, we observe that for the cases where we increase noise, see sec:results:induced-degeneracy:noise, there is a threshold value for increasing the noise level after which a rapid decline in fitness occurs - which we interpret as rapid aging-/dying. Indeed, at the same moment, the SE is decreasing, showing a complete disappearance of the anatomy towards a quasi-random distribution of the states. In this case, the AIS and TE increase around 600 steps, a bit more of half-life. A high TE can be interpreted as a high entropy of actions, the cells trying to act after the decline and disappearance of the anatomy.

However, the loss of competency of the cells does not impact the AIS and TE immediately as in the noise case. It is delayed and is starting to have a lot of random states in the anatomy with an increase of the AIS and TE at the end of the lifespan of the organism . But when noise increases and loss of competency are associated the increase in AIS and TE during aging is higher than in the case when only noise is added. Maybe, the loss of competency alters the capacity of the cells to understand the decline of the anatomy and their potentials to counteract it. Noise in the other hand clearly creates a lot of cellular actions during lifespan.

When the connectivity of the cells is altered, from 500 steps, we can observe an increase in AIS and more variability on TE (decrease and increase during aging). By closing the gap junctions, the cells are not receiving information from their neighbors gradually, increasing the uncertainty about their actions. We can see with the SE that the cells can counteract the loss of connectivity until 500 steps until an important decrease in spatial entropy and a complete inversion of the states of the anatomy.

The case that is creating the highest variability in anatomical outcome is the one where we altered the ANN parameters. AIS decreases and then increases over time slightly, TE is close to 0 after development and SE shows different types of aging, a loss of structure (organs) and proliferation of a specific states in the tissue (cancer), and a random activation of cell states.

We can conclude that if the states of the organisms do not change over time, if the anatomy reached a steady state, the AIS and TE are low. AIS is low because there is no new information to be gained from observing past states—future states are always the same and known, independent of the past. Essentially, the predictability is absolute, but since the system does not utilize past information to determine this (the state is intrinsically static), the AIS has low values. This indicates that in a non-changing system there is no information flow necessary between past and future states. The future state is always known (it’s the same as the current state), and thus, the storage of information from past states is irrelevant for future predictions. However, when AIS increases, it is associated to an increase in TE too. It happens in the cases of catastrophic changes in the morphology, as related to changes in spatial entropy. This important change implies an important information flow between cells, possibly representing the cellular attempts to recover from the high morphological changes. The absence of information flow in the default behavior can indicate that the systems considered that it reached its goal, e.g., development and morphostasis is not one of its goal, leading to a steady decrease of the anatomy, or in other words, aging. SE measures the entropy of the grid states. When SE is very high it corresponds to a random distribution of the states over the tissue. This random distribution can be found mainly in the cases of the loss of competency, noise and when both perturbations of associated. SE reveals also a different kind of aging corresponding to a loss of structure, with the loss of specific organs or the proliferation of a specific states in the tissue, like in the loss of connectivity and ANN parameters cases, meaning that the different perturbations lead to different kind of morphological degradation.

4.4. Regeneration as the Cure for Aging?

4.4.1. Loss of Organs Does Not Imply the Loss of Information About the Organ

In biological systems, certain cells or tissues retain “memories” of past states, influencing future behaviors. For instance, epigenetic memory in cells allows them to “remember” previous environmental exposures or developmental cues. Even if the external trigger is removed, cells can maintain certain gene expression patterns or epigenetic marks [90]. Similarly, tissues may retain information about past structural organization. For example, some amphibians can regenerate limbs or eyes after they are entirely removed [91,92,93]. Planarian flatworms exhibit a kind of holographic memory where, upon cutting into many pieces, each fragment of the animal regenerates everything that is missing to create a perfect and complete copy [94,95]. This persistence of information about the system’s target morphology, which serves as the setpoint for the process of anatomical homeostasis, can be implemented in a variety of biochemical, bioelectrical, and biomechanical information, thus it is important to study this important question.

In order to analyze the information dynamics during loss of organs, we computed the spatial AIS and TE (see Figure 7). We observed that the loss of the eyes during aging does not imply the loss of the AIS at the location of the eye as we can see at step 800 of Figure 7. The AIS is high and located at the position of the eyes in 100% of the cases. This suggests a memory retention in the system. For example, high AIS values at the “eye” locations in our system parallels the prepattern of the developing amphibian face consisting of regions of specialized cellular resting potential that determine the organs [96]. Even after visible features are lost, underlying patterns (like bioelectric signals or AIS) can persist, reflecting the system’s inherent memory of structure. AIS measures how much information about the past states is retained in a system over time. This retention can create a lingering effect, where regions that previously had distinct structures (like the eyes) continue to show high AIS values. Even if the explicit structure is no longer present (e.g., the “eyes” disappear in the smiley face), the underlying system dynamics may still carry information related to those past configurations. This residual AIS indicates that the system “remembers” these regions as important parts of its past structure, maintaining high information storage due to the prior role they played.

We conclude that the loss of organs does not imply the loss of information about the organ (see example Figure 7). Following this result we decided to explore the implications for rejuvenation by inducing the regeneration of the organs in the system in the next section.

4.4.2. Implication for Rejuvenation: A Simulated Experiment of Organ Restoration

The main reason for the long-term fitness degradation in our experiments - as discussed in sec:results:antagonist-pleiotropy, i.e.,, even without artificially affected cellular competencies - is the loss of specific internal organs, such as the eyes or the mouth. Thus, we here aim to correspondingly counteract the loss of detrimentally affected organs by replacing either the entire organ or specifically reconfiguring the affected cells of the NCA with “regenerative” state information (also erasing their internal working memory by resetting their RNN state information, see appendices A and B for details).

To achieve this, we define an organ-specific score function, $r_{\mathrm{o}}\left(t_k\right)=N_{\widehat{\mathbf{g}}}\left(t_k\right)/{\widehat{N}}_{\mathrm{o}}$, that quantifies for every time-step $t_k$ the fraction of correctly expressed target cells $N_{\widehat{\mathbf{g}}}\left(t_k\right)$ out of the total number of target cells ${\widehat{N}}_o$ of a specific organ, $o\in \{\mathit{LeftEye},\mathit{RightEye},\mathit{Mouth}\}$ (c.f., blue, green, and red cells in the snapshot of the $(16\times 16)$ smiley-face pattern NCA in fig:results:restoration:organ, respectively). Whenever an organ’s score falls below a predefined threshold value of $r_{\mathrm{o}}\left(t_k\right)<r_{\mathrm{th}}$, we perform an intervention to the affected organ, aiming at restoring its morphology locally. To interfere as little as possible with the NCA’s inherent developmental program, we propose utilizing the initial cell states as the regenerative state information of affected cells, ${\mathbf{c}}_i^{*}\left(t_k\right)\to {\mathbf{c}}_i\left(0\right)$ (c.f., sec:methods:nca). This corresponds to reconfiguring the affected cells to their initial (primordial) state at the onset of the developmental program of the NCA .

4.4.3. Less Is More: Organ Restoration Induction by Injecting the Regenerative Information Only to Incorrectly-Positioned Cells

We thus deploy the NCA solution discussed in sec:results:antagonist-pleiotropy (c.f., fig:result:pleiotropy) for an extended life-span of $t_A=1000$ time steps and apply the here discussed rejuvenation interventions whenever an organ’s score $r_{\mathrm{o}}\left(t_k\right)$ falls below a threshold value of $r_{\mathrm{th}}=0.6$. In fig:results:restoration:organ (A) and (B) we respectively present the results of this procedure by either (A) resetting all cells of an affected organ with their initial configuration or (B) only resetting the particular cells of an affected organ that express different cell states than their target types. We present for both (A) and (B) the long-term fitness scores of 10 statistically independent simulations, contrast their mean behavior with the baseline case without interventions (see fig:result:pleiotropy), and visualize the organ-specific scores and the corresponding intervention events for selected individuals.

We see that method (A) is not very sustainable: Although we generally measure significant short-term improvements in a respective organs’ scores after such intervention events, the intervened organs dissolve again quickly so that after $\approx 25$ time steps (which we set as the lower limit of two consecutive interventions, allowing for unperturbed development during this period) the organ needs to be reset again, and again. This periodic need for interventions - also reflected in correspondingly oscillating fitness scores of the entire organism - seems fundamental in our system and is not significantly affected by the delay between two interventions. The resulting fitness scores of the intervened individuals are only slightly improved compared to the baseline case.

Strikingly, if only the wrongly expressed cells of an affected organ are reconfigured to their initial state, as in case (B), very few of such targeted interventions (on the order of $\approx 5-15$, c.f., bar plots in fig:results:restoration:organ (B)) suffice to sustainably restore degrading organs and allow the corresponding virtual organism to maintain its morphological integrity over an extended lifespan. The corresponding long-term fitness scores of the rejuvenated individuals show significant improvement compared to the baseline case and can even be stabilized into a steady state of a constant value of, here, $r_{\mathrm{A}}\approx 200$; notably, the difference to the maximum fitness score of $r_{\max }=256$ right after the developmental stage at $t_{\mathrm{D}}=35$ arises due to the accumulative stagnation loss term $r_{\mathrm{S}}$ in eq:fitness.

4.4.4. Boundaries Matter: Organ Restoration Is More Efficient with the Injection of a Differential Pattern Including the Organ Cell States and Neighboring Cell States

Although the above approach is already promising, we further investigate how including the cells at the boundary to adjacent tissue surrounding the respective organs affects regeneration in our procedure. Thus, we additionally consider the respective boundary cells for the left- and right eyes in the evaluation of the respective organs’ scores but intentionally omit the boundaries of the mouth (c.f., emphasized color-coded cells in fig:results:restoration:socket, similar to fig:results:restoration:organ); notably, the mouth is much more stable than the eyes, and its cells already occupy a substantial fraction of the face.

Deploying the same experiments as above in sub:results:restoration:organ, we see in fig:results:restoration:socket that (A) interventions reconfiguring all the cells of affected organs (now also including their respective surrounding tissue) show similar performance to the corresponding case shown in fig:results:restoration:socket (A). However, the more targeted cell-specific interventions depicted in fig:results:restoration:socket (B) of now manipulating not only the wrongly expressed cells of an affected organ but also the wrongly expressed cells of their respective boundaries greatly improve the performance of restoration interventions compared to the cases fig:results:restoration:organ,fig:results:restoration:socket (A) and fig:results:restoration:socket (B). Here, an average number of as little as $5-6$ interventions (and in some cases even one single or two) suffice to stabilize the morphology of the here investigated $(16\times 16)$ smiley-face pattern organism in long lifespan experiments.

We conclude that boundaries (and more generally relations) between organs matter in our regeneration experiments: including surrounding tissue in detecting affected organs and in subsequent regeneration experiments that specifically target the wrongly expressed cells of the affected organs can be an efficient approach for regeneration and rejuvenation applications in our system.

A different but likewise promising approach is replacing the affected organs or cells directly with their target states (rather than their respective initial states as here). However, while “transplanting” entire organs can be effective, there is little quantitative difference in rejuvenation performance between replacing cells with their target states and reconfiguring cells into their initial states. From a biological perspective, the latter approach of resetting cells to their primordial states strikes us as conceptually and technologically more feasible, as the number of specialized target states might greatly exceed the number of necessary primordial initial states: resetting cells to their initial cell state, or at least to a particular primordial cell state suitable for starting a specific morphogenetic task, could serve as a general fallback mechanism to induce regeneration of the respective (here virtual) organism by restarting the developmental program locally.

At this point, we would like to stress that the developmental program of our NCA has neither been optimized for such long lifespans nor for the here proposed intervention experiments. The latter can thus be seen as minimally invasive ad-hoc cellular reconfiguration interventions, i.e., a dynamic and targeted “hacking” of the organism’s developmental plan by an external agent [19] (i.e.,, us, in this case).

Overall, we conclude that such a procedure of targeted cellular interventions is a reliable approach for rejuvenation in our system.

5. Discussion and Conclusion

In this paper, we aim to address a significant gap in the aging literature which so-far lacks - to the best of our knowledge – computational models of aging considering cellular information processing and multi-scale competency. We seek to address this knowledge gap by developing an evolutionary framework integrating the multi-scale competency architecture of multicellular tissue.

Specifically, we here utilize numerical tools from Artificial Life [61], so-called Neural Celular Automata (NCAs) [39,41,67] which we trained with an evolutionary algorithm to perform self-orchestrated targeted pattern formation, i.e.,, morphogenesis of a $(16\times 16)$ smiley-face pattern. The localized unicellular decision-making centers in our NCA implementation are modeled by Artificial Neural Networks (ANNs) which perceive the numerical states of their direct cellular neighbors on the grid of the NCA to regulate their own cell state individually. In that way, NCAs substitute homeostatic and morphogenetic pathways of “real” multicellular biological tissue, the NCA-specific ANNs and numerical states respectively corresponding to cell-internal decision-making machinery (such as gene regulatory networks) and cell-type expressions.

Analogously to their biological counterparts, the NCA’s cells form their decisions following a global agenda, allowing the collective of cells to assume a predefined target pattern over a predefined developmental time, thus simulating the process of embryogenesis. The evolutionary algorithm selects mature phenotypes directly after the predefined developmental stage (i.e.,, fully grown tissue), specifically quantifying how well the cells of the corresponding NCAs can grow the facial features and organs of the $(16\times 16)$ smiley-face pattern in a fully self-orchestrated and decentralized way. Notably, the NCA’s cells do not in any way explicitly perceive time or their age (in terms of simulation time steps), but their dynamics is fully governed by regulating and exchanging noisy state information on the NCA’s grid. Yet, when deploying a successfully trained NCA for much longer times than necessary for development (at which stage evolutionary selection would have occurred), the tissue of the NCA expresses signs of aging, such as successively blurring boundaries between- or entirely losing a single or several organs. We emphasize that this behavior is emergent and happens without any further intervention.

We propose, that even in the absence of accumulated cellular or genetic damage, a biological system left without any goal in the morphospace will start degrading anatomically. Using an in silico model of homeostatic morphogenesis with a multi-scale competency architecture and information dynamics analysis, we find the following, which provides a novel perspective on aging dynamics with significant implications for longevity research and regenerative medicine: (1) Absence of Long-Term Morphostasis goal: Aging emerges naturally after development due to the lack of an evolved regenerative goal, rather than being caused by specific detrimental properties of developmental programs (e.g., antagonistic pleiotropy or hyperfunction); (2) Acceleration Factors vs. Root Cause: Cellular misdifferentiation, reduced competency, communication failures, and genetic damage all accelerate aging but are not its primary cause; (3) Information Dynamics in Aging: Aging correlates with increased active information storage and transfer entropy, while spatial entropy measures distinguish two dynamics, (i) the loss of structure and (ii) morphological noise accumulation; (4) Dormant Regenerative Potential: Despite organ loss, spatial information persists in our cybernetic tissue, indicating a memory of lost structures, which can be reactivated for organ restoration through targeted regenerative information; and (5) Optimized Regeneration Strategies: Restoration is most efficient when regenerative information includes differential patterns of affected cells and their neighboring tissue, highlighting strategies and protocols for rejuvenation.

While the cells in our system learn to perform morphogenesis, i.e, how to grow an adult morphology successfully, we still observe aging effects in long-term simulations even without artificially perturbing the NCA’s parameters. These results suggest that development and morphostasis, while close in terms of morphological trajectories and tasks still significantly diverge. To counteract aging, it is therefore necessary to (re)activate a developmental/regenerative goal. We explicitly demonstrated this by providing our cybernetic tissue with the necessary information for regeneration, thereby inducing rejuvenation through self-organization at the anatomical level. Our results suggest a roadmap for rejuvenation therapies, particularly that curing aging may imply to provide a new goal to the tissue, i.e., regeneration. In other words, this suggests that solving regeneration solves aging. As currently debated, regenerative species, such as salamanders or planaria, might not age at all. While, already known for planaria, this is the subject of current research for salamander. The latter not only lack senescent phenotypes but it is indeed impossible to define an epigenetic clock for specimen older than four years [97]. One extension to the present work thus will be comparing a system that has learned to regenerate and actively maintain morphostasis after development to an aging system that has only been trained on morphogenesis tasks. Understanding their differences in terms of information dynamics will allow us to form a thorough understanding of the homeostatic loops that allow differently expressed cells in a self-assembling tissue to maintain their integrity against adversarial attacks or thermodynamic noise.

This work also emphasized how to optimize regenerative strategies. Our results show the importance of the information at the level of the borders of integrated tissue or organs for regeneration. This importance can be found in the field of bioelectricity. Indeed, research indicates that mutations in critical neurogenesis genes, such as Notch, can lead to developmental defects [98]. These organ defects can be cured through the overexpression or activation of native HCN2 (hyperpolarization-activated cyclic nucleotide-gated) channels [99,100,101]. As voltage-sensitive ion channels, HCN2 channels respond to their cellular environment by enhancing hyperpolarization in cells that are slightly polarized, while remaining inactive in depolarized cells. This mechanism functions analogously to “contrast enhancers,” which sharpen weak differences in membrane potential (${\mathrm{V}}_{\mathrm{mem}}$) across compartment boundaries, much like a sharpening filter brings clarity to a blurred image. In other words, HCN2 enhances and sharpens the borders in the bioelectrical pattern.

Moreover, our results suggest that a multicellular collective is capable of maintaining the memory about fully degraded features of an organism, as explicitly demonstrated here via organ regeneration experiment and information dynamics analysis. This suggests that this information is dormant and can be “retrieved” by injecting the appropriate signals that trigger regeneration in our system. This finding opens up exciting possibilities for future research, particularly in the realm of regenerative medicine and bioengineering. By understanding how cells store and retrieve hierarchical anatomical information, we could potentially develop new methods to promote regeneration in damaged tissues or organs. Additionally, exploring the mechanisms by which cells communicate and coordinate to maintain this memory could lead to breakthroughs in preventing age-related degeneration. For example, targeting specific bioelectrical patterns or enhancing the function of key ion channels like HCN2 could help maintain tissue integrity and function over time. Furthermore, the here discussed principles targeting the “software” of life could be applied to synthetic biology and rejuvenation, where designing artificial systems with similar regenerative capabilities [102,103] as well as new aging treatments might become feasible [26,34].

While our current framework is capable of self-assembling target patterns of heterogeneously distributed cell types, the ANN-based decision-making centers of our NCA’s cells are implemented homogeneously: all cells of an NCA typically share the same ANN structure and parameters, and thus the same basic functionality. Yet, the cells might be configured differently and thus effectively behave heterogeneously due to cell-specific state expressions and input relations, in short, due to their morphogenetic history. In that way, these artificial cells’ functionality mirrors gene regulatory networks. Thus, further in silico studies could investigate minimally complex high-level interventions necessary to induce robust regenerative effects in affected tissue simply by reconfiguring the cellular collective locally. Furthermore, allowing cell-specific functional mutations during the lifetime of such virtual organisms might be an intriguing way to study morphogenetic and aging-related implications of heterogeneity.

Our results validate aspects of both damage-based and programmatic theories of aging [1,4,5] but in very different ways. On the one hand, accelerated morphological decline (or an increased rate of aging) in our system mainly originates from: increasing cellular misdiffentiation, decreasing cellular competency, defective cell-cell communication, and accumulation of genetic damage. While these four characteristics align with damage-based theories of aging [6], they are only sufficient not necessary, i.e., not really causal, as our system undergoes aging also in their absence. Therefore, our model tends towards programmatic theories of aging, antagonistic pleiotropy or hyperfunction [14,15], or the theories of the software design flaw [5,26], in which the causes of aging have to be found in development. Similarly, but not analogously, we here relate aging to the absence of a developmental/regenerative or morphostatic goal, which nonetheless may imply an abnormal continuation of specific developmental pathways that are not beneficial for anatomical homeostasis. However, our results suggest that developmental programs (that are learned during evolution) can also qualify as regenerative pathways when “reactivated” in aging individuals, thus representing a potential protocol for a cure for aging ref. In a sense, regeneration is learned for free and is just dormant after development. This is an important result and indicates that we should focus our research efforts for curing aging on the reactivation of the regenerative capabilities. And our results suggest an experimental research program for regenerative medicine where we create an embryo-like environment allowing the tissue to regenerate. In turn, when interpreting the NCAs’ numerical states as bioelectric degrees of freedom, instead of transcriptomic cell-state expressions, a promising approach toward innovative no-aging or rejuvenation treatments might be to focus on high-level control layers such as bioelectricity that encode collective anatomical goals [26,102].

Our study showed fundamental computational principles in multicellular aging systems, bridging developmental biology with aging research. We explicitly demonstrate that developmental pathways, initially used for morphogenesis, may be dormant and can be reactivated for regeneration and tissue maintenance later in life. Moreover, we observed that cells collectively retain and manipulate morphological information, indicating that targeted interventions could potentially rejuvenate or maintain tissue integrity as organisms age.

Overall, we focused on the loss of goal-directedness in morphospace as the cause of the degradation of the system in this article. However, it may be possible that this cause extends towards nother levels in biological systems but also to artificial ones. For example, it may be possible that to counteract psychological decline during aging by giving new life goals to organism, which is in line with resistance to cognitive decline that has been found in elderly having high social engagement and new activities [104,105,106]. This can take place by providing new patterns (as occurs in metamorphosis [107,108]) within a body, or via communication of signals between bodies, as recently observed in cross-embryo morphogenetic error correction effects [109]. In terms of learning theory, we could ask: what’s the objective of a goal-directed system once it reached its goal? Our results may point towards the idea that a goal-directed system that achieved its goal becomes an aging system in a broad sense: a system deviating from its goal slowly. For a morphogenetic system, this might be reflected by a slowly degrading anatomy, a slight shift in morphospace [21] ultimately leading to death: an anatomical configuration incompatible with life. For AI systems, a decline in performance, similarly to what we found with ChatGPT or other LLMs [110], could be interpreted as aging or developmental defects after the target morphology has been reached [39]. Overall, our results suggest that reaching a goal by generative means is completely different from maintaining it. Thus, to engineer robust and autonomous AI and intelligent robotic systems or synthetic morphology we should integrate the learning of these two different tasks.

Exploring these principles further could lead to significant advancements in delaying or reversing age-related decline. Integrating computational models with experimental biology might pave the way to uncover new strategies to enhance regeneration through manipulating cellular behavior, potentially extending the healthy lifespan of tissues.

Future research will focus on applying these computational insights to real-world biological systems, aiming to develop therapies that modulate aging processes and improve recovery from injuries, thus optimizing the health span of multicellular organisms. These findings not only deepen our understanding of the aging process but also highlight the potential of using computational approaches to predict and control the dynamics of cellular systems in aging. The insights gained from this study open new avenues for research in aging and regenerative medicine, paving the way for novel interventions that could one day transform our approach to aging and chronic disease treatment.