The ubiquitin-proteasome system in retinal health and disease

The ubiquitin–proteasome system (UPS) is the main intracellular pathway for modulated protein turnover, playing an important role in the maintenance of cellular homeostasis. It also exerts a protein quality control through degradation of oxidized, mutant, denatured, or misfolded proteins and is involved in many biological processes where protein level regulation is necessary. This system allows the cell to modulate its protein expression pattern in response to changing physiological conditions and provides a critical protective role in health and disease. Impairments of UPS function in the central nervous system (CNS) underlie an increasing number of genetic and idiopathic diseases, many of which affect the retina. Current knowledge on the UPS composition and function in this tissue, however, is scarce and dispersed. This review focuses on UPS elements reported in the retina, including ubiquitinating and deubiquitinating enzymes (DUBs), and alternative proteasome assemblies. Known and inferred roles of protein ubiquitination, and of the related, SUMO conjugation (SUMOylation) process, in normal retinal development and adult homeostasis are addressed, including modulation of the visual cycle and response to retinal stress and injury. Additionally, the relationship between UPS dysfunction and human neurodegenerative disorders affecting the retina, including Alzheimer's, Parkinson's, and Huntington's diseases, are dealt with, together with numerous instances of retina-specific illnesses with UPS involvement, such as retinitis pigmentosa, macular degenerations, glaucoma, diabetic retinopathy (DR), and aging-related impairments. This information, though still basic and limited, constitutes a suitable framework to be expanded in incoming years and should prove orientative toward future therapy design targeting sight-affecting diseases with a UPS underlying basis.

Exploring oxidative modifications of tyrosine:an update on mechanisms of formation, advances in analysis and biological consequences

Protein oxidation is increasingly recognised as an important modulator of biochemical pathways controlling both physiological and pathological processes. While much attention has focused on cysteine modifications in reversible redox signalling, there is increasing evidence that other protein residues are oxidised in vivo with impact on cellular homeostasis and redox signalling pathways. A notable example is tyrosine, which can undergo a number of oxidative post-translational modifications to form 3-hydroxy-tyrosine, tyrosine crosslinks, 3-nitrotyrosine and halogenated tyrosine, with different effects on cellular functions. Tyrosine oxidation has been studied extensively in vitro, and this has generated detailed information about the molecular mechanisms that may occur in vivo. An important aspect of studying tyrosine oxidation both in vitro and in biological systems is the ability to monitor the formation of oxidised derivatives, which depends on a variety of analytical techniques. While antibody-dependent techniques such as ELISAs are commonly used, these have limitations, and more specific assays based on spectroscopic or spectrometric techniques are required to provide information on the exact residues modified and the nature of the modification. These approaches have helped understanding of the consequences of tyrosine oxidation in biological systems, especially its effects on cell signalling and cell dysfunction, linking to roles in disease. There is mounting evidence that tyrosine oxidation processes are important in vivo and can contribute to cellular pathology.

Experiment: Acidified seawater impacts sea urchin larvae pH regulatory systems relevant for calcification

Calcifying echinoid larvae respond to changes in seawater carbonate chemistry with reduced growth and developmental delay. To date, no information exists on how ocean acidification acts on pH homeostasis in echinoderm larvae. Understanding acid-base regulatory capacities is important because intracellular formation and maintenance of the calcium carbonate skeleton is dependent on pH homeostasis. Using H(+)-selective microelectrodes and the pH-sensitive fluorescent dye BCECF, we conducted in vivo measurements of extracellular and intracellular pH (pH(e) and pH(i)) in echinoderm larvae. We exposed pluteus larvae to a range of seawater CO(2) conditions and demonstrated that the extracellular compartment surrounding the calcifying primary mesenchyme cells (PMCs) conforms to the surrounding seawater with respect to pH during exposure to elevated seawater pCO(2). Using FITC dextran conjugates, we demonstrate that sea urchin larvae have a leaky integument. PMCs and spicules are therefore directly exposed to strong changes in pH(e) whenever seawater pH changes. However, measurements of pH(i) demonstrated that PMCs are able to fully compensate an induced intracellular acidosis. This was highly dependent on Na(+) and HCO(3)(-), suggesting a bicarbonate buffer mechanism involving secondary active Na(+)-dependent membrane transport proteins. We suggest that, under ocean acidification, maintained pH(i) enables calcification to proceed despite decreased pH(e). However, this probably causes enhanced costs. Increased costs for calcification or cellular homeostasis can be one of the main factors leading to modifications in energy partitioning, which then impacts growth and, ultimately, results in increased mortality of echinoid larvae during the pelagic life stage.

Rôle de la protéine DUSP5 dans l’autophagie des cardiomyocytes

Résumé: L’autophagie est un processus essentiel au maintien de l’homéostasie cellulaire. Elle permet de dégrader et recycler aussi bien des organelles entières que des composants cytoplasmiques non fonctionnels. De plus, l’augmentation d’autophagie en condition de stress constitue une réponse adaptative favorisant la survie cellulaire. Chez les cardiomyocytes, l’autophagie en condition basale est indispensable au renouvellement, entre autres, des mitochondries et des protéines formant les sarcomères. De plus, les stress tels l’ischémie cardiaque ou la carence en nutriments induisent une augmentation de l’autophagie protectrice. Dans certaines conditions extrêmes, il a été suggéré qu’un surcroît d’autophagie puisse toutefois exacerber la pathologie cardiaque en provoquant la mort des cardiomyocytes. Considérant l’importance de ce processus dans la physiopathologie cardiaque, l’identification des mécanismes signalétiques régulant l’autophagie chez les cardiomyocytes a été le sujet de recherches intenses. À cet effet, l’activation des Mitogen-Activated Protein Kinase (MAPK) a été démontrée pour réguler, avec d’autres voies signalétiques, l’autophagie et l’apoptose des cardiomyocytes. Il est donc probable que les Dual-Specificity Phosphatase (DUSP), enzymes clés contrôlant l’activité des MAPK, participent aussi à la régulation de l’autophagie. Afin de vérifier cette hypothèse, nous avons induit l’autophagie chez des cardiomyocytes isolés de rats nouveau-nés en culture. L’analyse de marqueurs d’autophagie par immunobuvardage démontre que l’activation des MAPK ERK1/2 et p38 corrèle avec l’activité autophagique chez les cardiomyocytes. Dans ces conditions, la diminution d’expression de la majorité des ARNm encodant les différentes DUSP retrouvées chez les cardiomyocytes contraste de façon marquée avec l’augmentation d’expression de l’ARNm Dusp5. De plus, nous avons démontré par une étude de gain de fonction que l’activation soutenue de p38 par surexpression d’un mutant MKK6 constitutivement actif stimule l’autophagie chez les cardiomyocytes. De façon surprenante, la perte de fonction de p38 obtenue par surexpression d’un mutant p38 dominant négatif n’altère en rien la réponse autophagique initiatrice dans notre modèle in vitro. Nos résultats suggèrent que les DUSP puissent réguler, via leurs actions sur les MAPK, d’importantes étapes du processus autophagique chez les cardiomyocytes.

DATA DRIVEN APPROACHES TO IDENTIFY DETERMINANTS OF HEART DISEASES AND CANCER RESISTANCE

Cancer and cardio-vascular diseases are the leading causes of death world-wide. Caused by systemic genetic and molecular disruptions in cells, these disorders are the manifestation of profound disturbance of normal cellular homeostasis. People suffering or at high risk for these disorders need early diagnosis and personalized therapeutic intervention. Successful implementation of such clinical measures can significantly improve global health. However, development of effective therapies is hindered by the challenges in identifying genetic and molecular determinants of the onset of diseases; and in cases where therapies already exist, the main challenge is to identify molecular determinants that drive resistance to the therapies. Due to the progress in sequencing technologies, the access to a large genome-wide biological data is now extended far beyond few experimental labs to the global research community. The unprecedented availability of the data has revolutionized the capabilities of computational researchers, enabling them to collaboratively address the long standing problems from many different perspectives. Likewise, this thesis tackles the two main public health related challenges using data driven approaches. Numerous association studies have been proposed to identify genomic variants that determine disease. However, their clinical utility remains limited due to their inability to distinguish causal variants from associated variants. In the presented thesis, we first propose a simple scheme that improves association studies in supervised fashion and has shown its applicability in identifying genomic regulatory variants associated with hypertension. Next, we propose a coupled Bayesian regression approach -- eQTeL, which leverages epigenetic data to estimate regulatory and gene interaction potential, and identifies combinations of regulatory genomic variants that explain the gene expression variance. On human heart data, eQTeL not only explains a significantly greater proportion of expression variance in samples, but also predicts gene expression more accurately than other methods. We demonstrate that eQTeL accurately detects causal regulatory SNPs by simulation, particularly those with small effect sizes. Using various functional data, we show that SNPs detected by eQTeL are enriched for allele-specific protein binding and histone modifications, which potentially disrupt binding of core cardiac transcription factors and are spatially proximal to their target. eQTeL SNPs capture a substantial proportion of genetic determinants of expression variance and we estimate that 58% of these SNPs are putatively causal. The challenge of identifying molecular determinants of cancer resistance so far could only be dealt with labor intensive and costly experimental studies, and in case of experimental drugs such studies are infeasible. Here we take a fundamentally different data driven approach to understand the evolving landscape of emerging resistance. We introduce a novel class of genetic interactions termed synthetic rescues (SR) in cancer, which denotes a functional interaction between two genes where a change in the activity of one vulnerable gene (which may be a target of a cancer drug) is lethal, but subsequently altered activity of its partner rescuer gene restores cell viability. Next we describe a comprehensive computational framework --termed INCISOR-- for identifying SR underlying cancer resistance. Applying INCISOR to mine The Cancer Genome Atlas (TCGA), a large collection of cancer patient data, we identified the first pan-cancer SR networks, composed of interactions common to many cancer types. We experimentally test and validate a subset of these interactions involving the master regulator gene mTOR. We find that rescuer genes become increasingly activated as breast cancer progresses, testifying to pervasive ongoing rescue processes. We show that SRs can be utilized to successfully predict patients' survival and response to the majority of current cancer drugs, and importantly, for predicting the emergence of drug resistance from the initial tumor biopsy. Our analysis suggests a potential new strategy for enhancing the effectiveness of existing cancer therapies by targeting their rescuer genes to counteract resistance. The thesis provides statistical frameworks that can harness ever increasing high throughput genomic data to address challenges in determining the molecular underpinnings of hypertension, cardiovascular disease and cancer resistance. We discover novel molecular mechanistic insights that will advance the progress in early disease prevention and personalized therapeutics. Our analyses sheds light on the fundamental biological understanding of gene regulation and interaction, and opens up exciting avenues of translational applications in risk prediction and therapeutics.

Chronic sustained hypoxia-induced redox remodeling causes contractile dysfunction in mouse sternohyoid muscle

Chronic sustained hypoxia (CH) induces structural and functional adaptations in respiratory muscles of animal models, however the underlying molecular mechanisms are unclear. This study explores the putative role of CH-induced redox remodeling in a translational mouse model, with a focus on the sternohyoid—a representative upper airway dilator muscle involved in the control of pharyngeal airway caliber. We hypothesized that exposure to CH induces redox disturbance in mouse sternohyoid muscle in a time-dependent manner affecting metabolic capacity and contractile performance. C57Bl6/J mice were exposed to normoxia or normobaric CH (FiO2 = 0.1) for 1, 3, or 6 weeks. A second cohort of animals was exposed to CH for 6 weeks with and without antioxidant supplementation (tempol or N-acetyl cysteine in the drinking water). Following CH exposure, we performed 2D redox proteomics with mass spectrometry, metabolic enzyme activity assays, and cell-signaling assays. Additionally, we assessed isotonic contractile and endurance properties ex vivo. Temporal changes in protein oxidation and glycolytic enzyme activities were observed. Redox modulation of sternohyoid muscle proteins key to contraction, metabolism and cellular homeostasis was identified. There was no change in redox-sensitive proteasome activity or HIF-1α content, but CH decreased phospho-JNK content independent of antioxidant supplementation. CH was detrimental to sternohyoid force- and power-generating capacity and this was prevented by chronic antioxidant supplementation. We conclude that CH causes upper airway dilator muscle dysfunction due to redox modulation of proteins key to function and homeostasis. Such changes could serve to further disrupt respiratory homeostasis in diseases characterized by CH such as chronic obstructive pulmonary disease. Antioxidants may have potential use as an adjunctive therapy in hypoxic respiratory disease.

Cellular Mechanisms of Sleep Homeostasis : Studies on Brain Energy Metabolism and Aging

Sleep is governed by a homeostatic process in which the duration and quality of previous wake regulate the subsequent sleep. Active wakefulness is characterized with high frequency cortical oscillations and depends on stimulating influence of the arousal systems, such as the cholinergic basal forebrain (BF), while cessation of the activity in the arousal systems is required for slow wave sleep (SWS) to occur. The site-specific accumulation of adenosine (a by-product of ATP breakdown) in the BF during prolonged waking /sleep deprivation (SD) is known to induce sleep, thus coupling energy demand to sleep promotion. The adenosine release in the BF is accompanied with increases in extracellular lactate and nitric oxide (NO) levels. This thesis was aimed at further understanding the cellular processes by which the BF is involved in sleep-wake regulation and how these processes are affected by aging. The BF function was studied simultaneously at three levels of organization: 1) locally at a cellular level by measuring energy metabolites 2) globally at a cortical level (the out-put area of the BF) by measuring EEG oscillations and 3) at a behavioral level by studying changes in vigilance states. Study I showed that wake-promoting BF activation, particularly with glutamate receptor agonist N-methyl-D-aspatate (NMDA), increased extracellular adenosine and lactate levels and led to a homeostatic increase in the subsequent sleep. Blocking NMDA activation during SD reduced the high frequency (HF) EEG theta (7-9 Hz) power and attenuated the subsequent sleep. In aging, activation of the BF during SD or experimentally with NMDA (studies III, IV), did not induce lactate or adenosine release and the increases in the HF EEG theta power during SD and SWS during the subsequent sleep were attenuated as compared to the young. These findings implicate that increased or continuous BF activity is important for active wake maintenance during SD as well as for the generation of homeostatic sleep pressure, and that in aging these mechanisms are impaired. Study II found that induction of the inducible NO synthase (iNOS) during SD is accompanied with activation of the AMP-activated protein kinase (AMPK) in the BF. Because decreased cellular energy charge is the most common cause for AMPK activation, this finding implicates that the BF is selectively sensitive to the metabolic demands of SD as increases were not found in the cortex. In aging (study III), iNOS expression and extracellular levels of NO and adenosine were not significantly increased during SD in the BF. Furthermore, infusion of NO donor into the BF did not lead to sleep promotion as it did in the young. These findings indicated that the NO (and adenosine) mediated sleep induction is impaired in aging and that it could at least partly be due to the reduced sensitivity of the BF to sleep-inducing factors. Taken together, these findings show that reduced sleep promotion by the BF contributes to the attenuated homeostatic sleep response in aging.

Antizyme Inhibitor 2 : A Novel Regulator of Cellular Polyamine Homeostasis

Polyamines are organic polycations that participate in various physiological functions, including cell proliferation, differentiation and apoptosis. Cellular polyamines originate from endogenous biosynthesis and exogenous sources. Their subcellular pool is under strict control, achieved by regulating their uptake and metabolism. Polyamine-induced proteins called antizymes (AZ) act as key regulators of intracellular polyamine concentration. They regulate both the transport of polyamines and the activity and degradation of ornithine decarboxylase (ODC), the rate-limiting enzyme in polyamine biosynthesis. AZs themselves are negatively regulated by antizyme inhibitor (AZIN). AZIN functions as a positive regulator of cellular polyamine homeostasis, which by binding to AZs reactivates ODC and induces the uptake of polyamines. In various pathological conditions, including cancer, polyamine levels are misregulated. Polyamine homeostasis has therefore become an attractive target for therapeutic interventions and it is thus crucial to characterize the molecular basis underlying the homeostatic regulation. A novel human AZIN-resembling protein was previously identified in our group. The purpose of this study was to elucidate the function and distribution of this protein, termed as an antizyme inhibitor 2 (AZIN2). According to my results, AZIN2 functions as a novel regulator of polyamine homeostasis. It shows no enzymatic activity, but instead it binds AZs and negates their activity, which subsequently leads to reactivation of ODC and inhibition of its degradation. Expression of AZIN2 is restricted to terminally differentiated cells, such as mast cells (MC) and neurosecretory cells. In these actively secreting cell types, AZIN2 localizes to subcellular vesicles or granules where its function is important for the vesicle-mediated secretion. In MCs, AZIN2 localizes to the serotonin-containing subset of MC granules, and its expression is coupled to MC activation. The functional role of polyamines as potential mediators of MC activity was also investigated, and it was observed that the secretion of serotonin is selectively dependent on activation of ODC. In neurosecretory cells, AZIN2-positive vesicles localize mainly to the trans-Golgi network (TGN). Depletion of AZIN2 or cellular polyamines causes selective fragmentation of the TGN and retards secretion of proteins. Since addition of exogenous polyamines reverses these effects, the data indicate that AZIN2 and its downstream effectors, polyamines, are functionally implicated in the regulation of secretory vesicle transport. My studies therefore reveal a novel function for polyamines as modulators of both constitutive and regulated secretion. Based on the results, I propose that the role of AZIN2 is to act as a local in situ activator of polyamine biosynthesis.

Melanocyte homeostasis in vivo tolerates Rb1 loss in a developmentally independent fashion

There has been uncertainty regarding the precise role that the pocket protein Rb1 plays in murine melanocyte homeostasis. It has been reported that the TAT-Cre mediated loss of exon 19 from a floxed Rb1 allele causes melanocyte apoptosis in vivo and in vitro. This is at variance with other findings showing, either directly or indirectly, that Rb1 loss in melanocytes has no noticeable effect in vivo, but in vitro leads to a semi-transformed phenotype. In this study, we show that Rb1-null melanocytes lacking exon 19 do not undergo apoptosis and survive both in vitro and in vivo, irrespective of the developmental stage at which Cre-mediated ablation of the exon occurs. Further, Rb1 loss has no serious long-term ramifications on melanocyte homeostasis in vivo, with Rb1-null melanocytes being detected in the skin after numerous hair cycles, inferring that the melanocyte stem cell population carrying the Cre-mediated deletion is maintained. Consequently, whilst Rb1 loss in the melanocyte is able to alter cellular behaviour in vitro, it appears inconsequential with respect to melanocyte homeostasis in the mouse skin.

Failure of amino acid homeostasis causes cell death following proteasome inhibition

The ubiquitin-proteasome system targets many cellular proteins for degradation and thereby controls most cellular processes. Although it is well established that proteasome inhibition is lethal, the underlying mechanism is unknown. Here, we show that proteasome inhibition results in a lethal amino acid shortage. In yeast, mammalian cells, and flies, the deleterious consequences of proteasome inhibition are rescued by amino acid supplementation. In all three systems, this rescuing effect occurs without noticeable changes in the levels of proteasome substrates. In mammalian cells, the amino acid scarcity resulting from proteasome inhibition is the signal that causes induction of both the integrated stress response and autophagy, in an unsuccessful attempt to replenish the pool of intracellular amino acids. These results reveal that cells can tolerate protein waste, but not the amino acid scarcity resulting from proteasome inhibition.

Unravelling Molecular and Cellular Disease Mechanisms in Infantile Neuronal Ceroid Lipofuscinosis (INCL)

Studying neurodegeneration provides an opportunity to gain insights into normal cell physiology, and not just pathophysiology. In this thesis work the focus is on Infantile Neuronal Ceroid Lipofuscinosis (INCL). It is a recessively inherited lysosomal storage disorder. The disease belongs to the neuronal ceroid lipofuscinoses (NCLs), a group of common progressive neurodegenerative diseases of the childhood. Characteristic accumulation of autofluorescent storage material is seen in most tissues but only neurons of the central nervous system are damaged and eventually lost during the course of the disease leaving most other cell types unaffected. The disease is caused by mutations in the CLN1 gene, but the physiological function of the corresponding protein the palmitoyl protein thioesterase (PPT1) has remained elusive. The aim of this thesis work was to shed light on the molecular and cell biological mechanisms behind INCL. This study pinpointed the localization of PPT1 in axonal presynapses of neurons. It also established the role of PPT1 in early neuronal maturation as well as importance in mature neuronal synapses. This study revealed an endocytic defect in INCL patient cells manifesting itself as delayed trafficking of receptor and non-receptor mediated endocytic markers. Furthermore, this study was the first to connect the INCL storage proteins the sphingolipid activator proteins (SAPs) A and D to pathological events on the cellular level. Abnormal endocytic processing and intracellular re-localization was demonstrated in patient cells and disease model knock-out mouse neurons. To identify early affected cellular and metabolic pathways in INCL, knock-out mouse neurons were studied by global transcript profiling and functional analysis. The gene expression analysis revealed changes in neuronal maturation and cell communication strongly associated with the regulated secretory system. Furthermore, cholesterol metabolic pathways were found to be affected. Functional studies with the knock-out mouse model revealed abnormalities in neuronal maturation as well as key neuronal functions including abnormalities in intracellular calcium homeostasis and cholesterol metabolism. Together the findings, introduced in this thesis work, support the essential role of PPT1 in developing neurons as well as synaptic sites of mature neurons. Results of this thesis also elucidate early events in INCL pathogenesis revealing defective pathways ultimately leading to the neurodegenerative process. These results contribute to the understanding of the vital physiological function of PPT1 and broader knowledge of common cellular mechanisms behind neurodegeneration. These results add to the knowledge of these severe diseases offering basis for new approaches in treatment strategies.

USP1 deubiquitinase: cellular functions, regulatory mechanisms and emerging potential as target in cancer therapy

12 p.

The E2F2 Transcription Factor Sustains Hepatic Glycerophospholipid Homeostasis in Mice

Increasing evidence links metabolic signals to cell proliferation, but the molecular wiring that connects the two core machineries remains largely unknown. E2Fs are master regulators of cellular proliferation. We have recently shown that E2F2 activity facilitates the completion of liver regeneration after partial hepatectomy (PH) by regulating the expression of genes required for S-phase entry. Our study also revealed that E2F2 determines the duration of hepatectomy-induced hepatic steatosis. A transcriptomic analysis of normal adult liver identified "lipid metabolism regulation" as a major E2F2 functional target, suggesting that E2F2 has a role in lipid homeostasis. Here we use wild-type (E2F2(+/+)) and E2F2 deficient (E2F2(-/-)) mice to investigate the in vivo role of E2F2 in the composition of liver lipids and fatty acids in two metabolically different contexts: quiescence and 48-h post-PH, when cellular proliferation and anabolic demands are maximal. We show that liver regeneration is accompanied by large triglyceride and protein increases without changes in total phospholipids both in E2F2(+/+) and E2F2(-/-) mice. Remarkably, we found that the phenotype of quiescent liver tissue from E2F2(-/-) mice resembles the phenotype of proliferating E2F2(+/+) liver tissue, characterized by a decreased phosphatidylcholine to phosphatidylethanolamine ratio and a reprogramming of genes involved in generation of choline and ethanolamine derivatives. The diversity of fatty acids in total lipid, triglycerides and phospholipids was essentially preserved on E2F2 loss both in proliferating and non-proliferating liver tissue, although notable exceptions in inflammation-related fatty acids of defined phospholipid classes were detected. Overall, our results indicate that E2F2 activity sustains the hepatic homeostasis of major membrane glycerolipid components while it is dispensable for storage glycerolipid balance.

On the Dynamics of the Adenylate Energy System: Homeorhesis vs Homeostasis

Biochemical energy is the fundamental element that maintains both the adequate turnover of the biomolecular structures and the functional metabolic viability of unicellular organisms. The levels of ATP, ADP and AMP reflect roughly the energetic status of the cell, and a precise ratio relating them was proposed by Atkinson as the adenylate energy charge (AEC). Under growth-phase conditions, cells maintain the AEC within narrow physiological values, despite extremely large fluctuations in the adenine nucleotides concentration. Intensive experimental studies have shown that these AEC values are preserved in a wide variety of organisms, both eukaryotes and prokaryotes. Here, to understand some of the functional elements involved in the cellular energy status, we present a computational model conformed by some key essential parts of the adenylate energy system. Specifically, we have considered (I) the main synthesis process of ATP from ADP, (II) the main catalyzed phosphotransfer reaction for interconversion of ATP, ADP and AMP, (III) the enzymatic hydrolysis of ATP yielding ADP, and (IV) the enzymatic hydrolysis of ATP providing AMP. This leads to a dynamic metabolic model (with the form of a delayed differential system) in which the enzymatic rate equations and all the physiological kinetic parameters have been explicitly considered and experimentally tested in vitro. Our central hypothesis is that cells are characterized by changing energy dynamics (homeorhesis). The results show that the AEC presents stable transitions between steady states and periodic oscillations and, in agreement with experimental data these oscillations range within the narrow AEC window. Furthermore, the model shows sustained oscillations in the Gibbs free energy and in the total nucleotide pool. The present study provides a step forward towards the understanding of the fundamental principles and quantitative laws governing the adenylate energy system, which is a fundamental element for unveiling the dynamics of cellular life.