Respiratory oscillations in yeast: clock-driven mitochondrial cycles of energization

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

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

Flavin-based metabolic cycles are integral features of growth and division in single yeast cells

The yeast metabolic cycle (YMC) is a fascinating example of biological organization, in which cells constrain the function of specific genetic, protein and metabolic networks to precise temporal windows as they grow and divide. However, understanding the intracellular origins of the YMC remains a challenging goal, as measuring the oxygen oscillations traditionally associated with it requires the use of synchronized cultures growing in nutrient-limited chemostat environments. To address these limitations, we used custom-built microfluidic devices and time-lapse fluorescence microscopy to search for metabolic cycling in the form of endogenous flavin fluorescence in unsynchronized single yeast cells. We uncovered robust and pervasive metabolic cycles that were synchronized with the cell division cycle (CDC) and oscillated across four different nutrient conditions. We then studied the response of these metabolic cycles to chemical and genetic perturbations, showing that their phase synchronization with the CDC can be altered through treatment with rapamycin, and that metabolic cycles continue even in respiratory deficient strains. These results provide a foundation for future studies of the physiological importance of metabolic cycles in processes such as CDC control, metabolic regulation and cell aging.

Temporal metabolic partitioning of the yeast and protist cellular networks: the cell is a global scale-invariant (fractal or self-similar) multioscillator

Britton Chance, electronics expert when a teenager, became an enthusiastic student of biological oscillations, passing on this enthusiasm to many students and colleagues, including one of us (DL). This historical essay traces BC's influence through the accumulated work of DL to DL's many collaborators. The overall temporal organization of mass-energy, information, and signaling networks in yeast in self-synchronized continuous cultures represents, until now, the most characterized example of in vivo elucidation of time structure. Continuous online monitoring of dissolved gases by direct measurement (membrane-inlet mass spectrometry, together with NAD(P)H and flavin fluorescence) gives strain-specific dynamic information from timescales of minutes to hours as does two-photon imaging. The predominantly oscillatory behavior of network components becomes evident, with spontaneously synchronized cellular respiration cycles between discrete periods of increased oxygen consumption (oxidative phase) and decreased oxygen consumption (reductive phase). This temperature-compensated ultradian clock provides coordination, linking temporally partitioned functions by direct feedback loops between the energetic and redox state of the cell and its growing ultrastructure. Multioscillatory outputs in dissolved gases with 13 h, 40 min, and 4 min periods gave statistical self-similarity in power spectral and relative dispersional analyses: i.e., complex nonlinear (chaotic) behavior and a functional scale-free (fractal) network operating simultaneously over several timescales.

Regulation of the Stress-Activated Degradation of Mitochondrial Respiratory Complexes in Yeast

Repair and removal of damaged mitochondria is a key process for eukaryotic cell homeostasis. Here we investigate in the yeast model how different protein complexes of the mitochondrial electron transport chain are subject to specific degradation upon high respiration load and organelle damage. We find that the turnover of subunits of the electron transport complex I equivalent and complex III is preferentially stimulated upon high respiration rates. Particular mitochondrial proteases, but not mitophagy, are involved in this activated degradation. Further mitochondrial damage by valinomycin treatment of yeast cells triggers the mitophagic removal of the same respiratory complexes. This selective protein degradation depends on the mitochondrial fusion and fission apparatus and the autophagy adaptor protein Atg11, but not on the mitochondrial mitophagy receptor Atg32. Loss of autophagosomal protein function leads to valinomycin sensitivity and an overproduction of reactive oxygen species upon mitochondrial damage. A specific event in this selective turnover of electron transport chain complexes seems to be the association of Atg11 with the mitochondrial network, which can be achieved by overexpression of the Atg11 protein even in the absence of Atg32. Furthermore, the interaction of various Atg11 molecules via the C-terminal coil domain is specifically and rapidly stimulated upon mitochondrial damage and could therefore be an early trigger of selective mitophagy in response to the organelles dysfunction. Our work indicates that autophagic quality control upon mitochondrial damage operates in a selective manner.

The role of alternative Polyadenylation in regulation of rhythmic gene expression

Background

Alternative transcription is common in eukaryotic cells and plays important role in regulation of cellular processes. Alternative polyadenylation results from ambiguous PolyA signals in 3′ untranslated region (UTR) of a gene. Such alternative transcripts share the same coding part, but differ by a stretch of UTR that may contain important functional sites.

Methods

The methodoogy of this study is based on mathematical modeling, analytical solution, and subsequent validation by datamining in multiple independent experimental data from previously published studies.

Results

In this study we propose a mathematical model that describes the population dynamics of alternatively polyadenylated transcripts in conjunction with rhythmic expression such as transcription oscillation driven by circadian or metabolic oscillators. Analysis of the model shows that alternative transcripts with different turnover rates acquire a phase shift if the transcript decay rate is different. Difference in decay rate is one of the consequences of alternative polyadenylation. Phase shift can reach values equal to half the period of oscillation, which makes alternative transcripts oscillate in abundance in counter-phase to each other. Since counter-phased transcripts share the coding part, the rate of translation becomes constant. We have analyzed a few data sets collected in circadian timeline for the occurrence of transcript behavior that fits the mathematical model.

Conclusion

Alternative transcripts with different turnover rate create the effect of rectifier. This “molecular diode” moderates or completely eliminates oscillation of individual transcripts and stabilizes overall protein production rate. In our observation this phenomenon is very common in different tissues in plants, mice, and humans. The occurrence of counter-phased alternative transcripts is also tissue-specific and affects functions of multiple biological pathways. Accounting for this mechanism is important for understanding the natural and engineering the synthetic cellular circuits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3958-1) contains supplementary material, which is available to authorized users.

Elements of the cellular metabolic structure

A large number of studies have demonstrated the existence of metabolic covalent modifications in different molecular structures, which are able to store biochemical information that is not encoded by DNA. Some of these covalent mark patterns can be transmitted across generations (epigenetic changes). Recently, the emergence of Hopfield-like attractor dynamics has been observed in self-organized enzymatic networks, which have the capacity to store functional catalytic patterns that can be correctly recovered by specific input stimuli. Hopfield-like metabolic dynamics are stable and can be maintained as a long-term biochemical memory. In addition, specific molecular information can be transferred from the functional dynamics of the metabolic networks to the enzymatic activity involved in covalent post-translational modulation, so that determined functional memory can be embedded in multiple stable molecular marks. The metabolic dynamics governed by Hopfield-type attractors (functional processes), as well as the enzymatic covalent modifications of specific molecules (structural dynamic processes) seem to represent the two stages of the dynamical memory of cellular metabolism (metabolic memory). Epigenetic processes appear to be the structural manifestation of this cellular metabolic memory. Here, a new framework for molecular information storage in the cell is presented, which is characterized by two functionally and molecularly interrelated systems: a dynamic, flexible and adaptive system (metabolic memory) and an essentially conservative system (genetic memory). The molecular information of both systems seems to coordinate the physiological development of the whole cell.

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.

Oscillation of oxygenation in skeletal muscle at rest and in light exercise

The aim of the present study was to compare the frequency of oxygenation determined in the vastus lateralis by near-infrared spectroscopy (NIRS) in light exercise with that at rest. A subject rested in a recumbent position for 5 min and changed body position to a sitting position on a cycle ergometer for 9 min. Then exercise with low intensity (work rate of 60% of maximal oxygen uptake) was carried out for 30 min. Total hemoglobin and myoglobin (THb/Mb) suddenly decreased after the start of exercise and gradually increased for 6 min. Oxygenated hemoglobin and myoglobin (Hb/MbO2) suddenly decreased and returned to a steady-state after the start of exercise. The difference between Hb/MbO2 and THb/Mb showed a sudden decrease and then a steady-state. This difference was analyzed by fast Fourier transform. The peak frequencies of the power spectrum density (PSD) were 0.0169 ± 0.0076 Hz at rest and 0.0117 ± 0.0042 Hz in exercise. The peak frequency of PSD was significantly decreased in exercise. In exercise, the range of frequencies was expanded. It is concluded that there are oscillations at rest as well as in exercise and that the frequency of peak PSD becomes lower in exercise than at rest.

The yin and yang of yeast transcription: elements of a global feedback system between metabolism and chromatin

When grown in continuous culture, budding yeast cells tend to synchronize their respiratory activity to form a stable oscillation that percolates throughout cellular physiology and involves the majority of the protein-coding transcriptome. Oscillations in batch culture and at single cell level support the idea that these dynamics constitute a general growth principle. The precise molecular mechanisms and biological functions of the oscillation remain elusive. Fourier analysis of transcriptome time series datasets from two different oscillation periods (0.7 h and 5 h) reveals seven distinct co-expression clusters common to both systems (34% of all yeast ORF), which consolidate into two superclusters when correlated with a compilation of 1,327 unrelated transcriptome datasets. These superclusters encode for cell growth and anabolism during the phase of high, and mitochondrial growth, catabolism and stress response during the phase of low oxygen uptake. The promoters of each cluster are characterized by different nucleotide contents, promoter nucleosome configurations, and dependence on ATP-dependent nucleosome remodeling complexes. We show that the ATP:ADP ratio oscillates, compatible with alternating metabolic activity of the two superclusters and differential feedback on their transcription via activating (RSC) and repressive (Isw2) types of promoter structure remodeling. We propose a novel feedback mechanism, where the energetic state of the cell, reflected in the ATP:ADP ratio, gates the transcription of large, but functionally coherent groups of genes via differential effects of ATP-dependent nucleosome remodeling machineries. Besides providing a mechanistic hypothesis for the delayed negative feedback that results in the oscillatory phenotype, this mechanism may underpin the continuous adaptation of growth to environmental conditions.

What yeast and cardiomyocytes share: ultradian oscillatory redox mechanisms of cellular coherence and survival

The coherent and robust, yet sensitively adaptable, nature of organisms is an astonishing phenomenon that involves massive parallel processing and concerted network performance at the molecular level. Unravelling the dynamic complexities of the living state underlines the essential operation of ultradian oscillations, rhythms and clocks for the establishment and maintenance of functional order simultaneously on fast and slower timescales. Non-invasive monitoring of respiration, mitochondrial inner membrane potentials, and redox states (especially those of NAD(P)H, flavin, and the monochlorobimane complex of glutathione), even after more than 50 years research, continue to provide both new insights and biomedical applications. Experiments with yeast and in cardiac cells reveal astonishing parallels and similarities in their dynamic biochemical organization.

The metabolic core and catalytic switches are fundamental elements in the self-regulation of the systemic metabolic structure of cells

BACKGROUND

Experimental observations and numerical studies with dissipative metabolic networks have shown that cellular enzymatic activity self-organizes spontaneously leading to the emergence of a metabolic core formed by a set of enzymatic reactions which are always active under all environmental conditions, while the rest of catalytic processes are only intermittently active. The reactions of the metabolic core are essential for biomass formation and to assure optimal metabolic performance. The on-off catalytic reactions and the metabolic core are essential elements of a Systemic Metabolic Structure which seems to be a key feature common to all cellular organisms.

METHODOLOGY/PRINCIPAL FINDINGS

In order to investigate the functional importance of the metabolic core we have studied different catalytic patterns of a dissipative metabolic network under different external conditions. The emerging biochemical data have been analysed using information-based dynamic tools, such as Pearson's correlation and Transfer Entropy (which measures effective functionality). Our results show that a functional structure of effective connectivity emerges which is dynamical and characterized by significant variations of bio-molecular information flows.

CONCLUSIONS/SIGNIFICANCE

We have quantified essential aspects of the metabolic core functionality. The always active enzymatic reactions form a hub--with a high degree of effective connectivity--exhibiting a wide range of functional information values being able to act either as a source or as a sink of bio-molecular causal interactions. Likewise, we have found that the metabolic core is an essential part of an emergent functional structure characterized by catalytic modules and metabolic switches which allow critical transitions in enzymatic activity. Both, the metabolic core and the catalytic switches in which also intermittently-active enzymes are involved seem to be fundamental elements in the self-regulation of the Systemic Metabolic Structure.

PSK2 coordinates glucose metabolism and utilization to maintain ultradian clock-coupled respiratory oscillation in Saccharomyces cerevisiae yeast

Ultradian clock-coupled respiratory oscillation (UCRO) in an aerobic continuous culture of Saccharomyces cerevisiae S288C is principally regulated by control of certain redox reactions of energy metabolism. It is also modulated by the metabolism of storage carbohydrates during adaptation to environmental change. However, the mechanism of cell sensing and response to environmental nutrients in UCRO is unknown. The purpose of the present study was to determine the role of PSK2 kinase in UCRO in yeast. S. cerevisiae in culture showed oscillation in PSK2 mRNA levels with a definite phase relationship to the respiratory oscillation. Furthermore, inactivation of Psk2 by gene disruption severely affected UCRO and its decline to undetectable levels within 2days. In addition, the extracellular and intracellular glucose concentrations of PSK2 deletion mutants in culture were higher and lower, respectively, than those of the wild type. PSK2 mutant cells showed no alteration in redox state. Furthermore, the levels of storage carbohydrates such as glycogen and trehalose fluctuated in PSK2 mutants with attenuated amplitudes comparable to those in the wild type. The results indicated that PSK2 kinase is important for the uptake of glucose and regulation of storage-carbohydrate synthesis and hence the maintenance of an unperturbed continuously oscillating state.

Quantitative analysis of cellular metabolic dissipative, self-organized structures

One of the most important goals of the postgenomic era is understanding the metabolic dynamic processes and the functional structures generated by them. Extensive studies during the last three decades have shown that the dissipative self-organization of the functional enzymatic associations, the catalytic reactions produced during the metabolite channeling, the microcompartmentalization of these metabolic processes and the emergence of dissipative networks are the fundamental elements of the dynamical organization of cell metabolism. Here we present an overview of how mathematical models can be used to address the properties of dissipative metabolic structures at different organizational levels, both for individual enzymatic associations and for enzymatic networks. Recent analyses performed with dissipative metabolic networks have shown that unicellular organisms display a singular global enzymatic structure common to all living cellular organisms, which seems to be an intrinsic property of the functional metabolism as a whole. Mathematical models firmly based on experiments and their corresponding computational approaches are needed to fully grasp the molecular mechanisms of metabolic dynamical processes. They are necessary to enable the quantitative and qualitative analysis of the cellular catalytic reactions and also to help comprehend the conditions under which the structural dynamical phenomena and biological rhythms arise. Understanding the molecular mechanisms responsible for the metabolic dissipative structures is crucial for unraveling the dynamics of cellular life.

Global self-regulation of the cellular metabolic structure

BACKGROUND

Different studies have shown that cellular enzymatic activities are able to self-organize spontaneously, forming a metabolic core of reactive processes that remain active under different growth conditions while the rest of the molecular catalytic reactions exhibit structural plasticity. This global cellular metabolic structure appears to be an intrinsic characteristic common to all cellular organisms. Recent work performed with dissipative metabolic networks has shown that the fundamental element for the spontaneous emergence of this global self-organized enzymatic structure could be the number of catalytic elements in the metabolic networks.

METHODOLOGY/PRINCIPAL FINDINGS

In order to investigate the factors that may affect the catalytic dynamics under a global metabolic structure characterized by the presence of metabolic cores we have studied different transitions in catalytic patterns belonging to a dissipative metabolic network. The data were analyzed using non-linear dynamics tools: power spectra, reconstructed attractors, long-term correlations, maximum Lyapunov exponent and Approximate Entropy; and we have found the emergence of self-regulation phenomena during the transitions in the metabolic activities.

CONCLUSIONS/SIGNIFICANCE

The analysis has also shown that the chaotic numerical series analyzed correspond to the fractional Brownian motion and they exhibit long-term correlations and low Approximate Entropy indicating a high level of predictability and information during the self-regulation of the metabolic transitions. The results illustrate some aspects of the mechanisms behind the emergence of the metabolic self-regulation processes, which may constitute an important property of the global structure of the cellular metabolism.

The number of catalytic elements is crucial for the emergence of metabolic cores

Background

Different studies show evidence that several unicellular organisms display a cellular metabolic structure characterized by a set of enzymes which are always in an active state (metabolic core), while the rest of the molecular catalytic reactions exhibit on-off changing states. This self-organized enzymatic configuration seems to be an intrinsic characteristic of metabolism, common to all living cellular organisms. In a recent analysis performed with dissipative metabolic networks (DMNs) we have shown that this global functional structure emerges in metabolic networks with a relatively high number of catalytic elements, under particular conditions of enzymatic covalent regulatory activity.

Methodology/Principal Findings

Here, to investigate the mechanism behind the emergence of this supramolecular organization of enzymes, we have performed extensive DMNs simulations (around 15,210,000 networks) taking into account the proportion of the allosterically regulated enzymes and covalent enzymes present in the networks, the variation in the number of substrate fluxes and regulatory signals per catalytic element, as well as the random selection of the catalytic elements that receive substrate fluxes from the exterior. The numerical approximations obtained show that the percentages of DMNs with metabolic cores grow with the number of catalytic elements, converging to 100% for all cases.

Conclusions/Significance

The results show evidence that the fundamental factor for the spontaneous emergence of this global self-organized enzymatic structure is the number of catalytic elements in the metabolic networks. Our analysis corroborates and expands on our previous studies illustrating a crucial property of the global structure of the cellular metabolism. These results also offer important insights into the mechanisms which ensure the robustness and stability of living cells.

Respiratory oscillations in yeasts

Respiratory oscillations in yeasts have been studied in three time domains with periods of (a) about a minute, (b) about 40 min, and (c) about a day. Reactive responses (damped oscillations), rhythms and temperature-compensated clocks have been described for (b) and (c), but a timekeeping clock has not yet been shown for (a). Synchronous populations reveal the time-structure that can only otherwise be studied in single organisms; this is because time-averaging through an asynchronous population conceals its fine structure. Early studies with synchronous cultures made by size selection methods indicated ultradian-clock driven oscillations in respiration, pools of adenylates, total protein, RNA synthesis and many enzyme activities (tau = 40 min in Schizosaccharomyces pombe, 30 min in Candida utilis), and more recently in self-synchronised continuous cultures of Saccharomyces cerevisiae (tau = 48 min). Most detailed understanding comes from the latter system, where continuous, noninvasive real-time monitoring (of 02 uptake, CO2 production, and NAD(P)H redox state) is combined with frequent discrete time samples (for other redox components, including H2S, GSH and cytochromes, metabolites, and mRNA levels). A redox switch lies at the heart of this ultradian clock and a plethora of outputs is optimized to a time-base that is genetically-determined and differs in different organisms. It is suggested that the entire temporal landscape of all eukaryotic organisms and the cells of higher plants and animals is constructed on this basis. A time frame for the coordination and coherence of all intracellular processes and the construction and assembly of cellular structures is provided by the ultradian clock The circadian clock matches these functions to the daily cycle of the external environment.

Collective behavior in gene regulation: the cell is an oscillator, the cell cycle a developmental process

The finding of a genome-wide oscillation in transcription that gates cells into S phase and coordinates mitochondrial and metabolic functions has altered our understanding of how the cell cycle is timed and how stable cellular phenotypes are maintained. Here we present the evidence and arguments in support of the idea that everything oscillates, and the rationale for viewing the cell as an attractor from which deterministic noise can be tuned by appropriate coupling among the many feedback loops, or regulons, that make up the transcriptional-respiratory attractor cycle. The existence of this attractor also explains many of the dynamic macroscopic properties of the cell cycle and appears to be the timekeeping oscillator in both cell cycles and circadian rhythms. The path taken by this primordial oscillator in the course of differentiation or drug response may involve period-doubling behavior. Evidence for a relatively high-frequency timekeeping oscillator in yeast and mammalian cells comes from expression array analysis, and GC/MS in the case of yeast, and primarily from macroscopic measures of phase response to perturbation in the case of mammalian cells. Low-amplitude, genome-wide oscillations, a ubiquitous but often unrecognized attribute of phenotype, may be a source of seemingly intractable biological noise in microarray and proteomic studies. These oscillations in transcript and protein levels and the repeated cycles of synthesis and degradation they require, represent a high energy cost to the cell which must, from an evolutionary point of view, be recovered as essential information. We suggest that the information contained in this genome-wide oscillation is the dynamic code that organizes a stable phenotype from an otherwise passive genome.

A tuneable attractor underlies yeast respiratory dynamics

Our understanding of the molecular structure and function in the budding yeast, Saccharomyces cerevisiae, surpasses that of all other eukaryotic cells. However, the fundamental properties of the complex processes and their control systems have been difficult to reconstruct from detailed dissection of their molecular components. Spontaneous oscillatory dynamics observed in self-synchronized continuous cultures is pervasive, involves much of the cellular network, and provides unique insights into integrative cell physiology. Here, in non-invasive experiments in vivo, we exploit these oscillatory dynamics to analyse the global timing of the cellular network to show the presence of a low-order chaotic component. Although robust to a wide range of environmental perturbations, the system responds and reacts to the imposition of harsh environmental conditions, in this case low pH, by dynamic re-organization of respiration, and this feeds upwards to affect cell division. These complex dynamics can be represented by a tuneable attractor that orchestrates cellular complexity and coherence to the environment.

Hydrogen sulfide: clandestine microbial messenger?

Although the toxicity of hydrogen sulfide (H(2)S) has been substantiated for almost 230 years, its pivotal roles in both aerobic and anaerobic organisms have only recently become evident. In low oxygen environments with millimolar concentrations of H(2)S, it functions as an electron donor and as an energy source in some systems. At micromolar levels, intracellular H(2)S in aerobic organisms has a vital role in redox balancing. At even lower concentrations, H(2)S provides essential signals in yeast, in the brain and in smooth and cardiac muscles. Here, other possible coordinating roles within and between microorganisms are suggested, including the possibility that H(2)S functions as a signalling mediator in prokaryotes. It is expected that future research will uncover a host of novel functions, not only in eukaryotes but also in prokaryotic species.

Modeling the mechanism of metabolic oscillations in ischemic cardiac myocytes

Oscillations in energy metabolism have been observed in a variety of cells under metabolically deprived conditions such as ischemia. In cardiac ventricular myocytes these metabolic oscillations may cause oscillations in the action potential duration, creating the potential for cardiac arrhythmias during ischemia (O'Rourke, 2000). A mathematical model of the mechanism behind metabolic oscillations is developed here. The model consists of descriptions of the mitochondrial components that regulate mitochondrial membrane potential (Psi), mitochondrial inorganic phosphate concentration, mitochondrial magnesium concentration, and cellular NADH and NAD(+) concentrations. Using parameters from the experimental literature, the model produces physiological values for these both under normoxic (steady state) and ischemic (oscillatory) conditions. The model includes the mitochondrial inner membrane anion channel (IMAC), the centum picosiemen channel (mCS), the phosphate carrier (PIC), and the respiration driven proton pumps. The model suggests that these are the essential components for producing oscillations with mCS essential for the rapid depolarization, PIC for the recovery from depolarization, and IMAC for the slow depolarization between depolarization peaks. A decrease of the inner membrane potential due to ischemia or experimental conditions seems to be a triggering factor for the oscillations. The model simulates the experimental observations that high levels of mitochondrial ADP and ATP abolish the oscillations, as does inhibition of electron transport. The model makes predictions on the influence of pH and magnesium levels on metabolic oscillations.

The temporal architecture of eukaryotic growth

Coherence of the time structure of growing organisms depends on a metronome-like orchestration. In a continuously perfused culture of Saccharomyces cerevisiae the redox state of the cell shows a temperature-compensated oscillation manifest in respiratory cycles, which are measured by continuous and non-invasive electrodes of probes such as dissolved oxygen and probes such as fluorometric NAD(P)H. Although the entire transcriptome exhibits low-amplitude oscillatory behaviour, transcripts involved in the vast majority of metabolism, stress response, cellular structure, protein turnover, mRNA turnover, and DNA synthesis are amongst the top oscillators and their orchestration occurs by an intricate network of transcriptional regulators. Therefore cellular auto-dynamism is a function of a large ensemble of excitable intracellular components of that self-organized temporally and spatially that encompasses mitochondrial, nuclear, transcriptional and metabolic dynamics, coupled by cellular redox state.

Effects of uncoupling of mitochondrial energy conservation on the ultradian clock-driven oscillations in Saccharomyces cerevisiae continuous culture

Protonophores have several different perturbative effects on dissolved O2 concentrations in continuous cultures of Saccharomyces cerevisiae. As well as uncoupling energy conservation from mitochondrial electron transport in vivo, they reset ultradian clock-driven respiratory oscillations and produce cell cycle effects. Thus, additions at low concentration (1.25 microM) of either m-chlorocarbonyl-cyanide phenylhydrazone (CCCP) or 5-chloro-3-t-butyl-2-chloro-4(1)-nitrosalicylanilide (S13) led to phase resetting of the 48 min ultradian clock-driven respiratory oscillations. At 2.5 microM CCCP or 4 microM S13, transient inhibition of oscillatory respiration (for 5 h) preceded synchronisation of the cell division cycle seen as a slow (9 h period) wave that enveloped the 48 min oscillation. At still higher concentrations of CCCP (5 microM), the cell division cycle was prolonged by about 7 h, and during this phase, the respiratory oscillation became undetectable. The significance of these observations with respect to the time-keeping functions of the ultradian clock is discussed.

Ultradian metronome: timekeeper for orchestration of cellular coherence

Dynamic intracellular spatial and temporal organization emerges from spontaneous synchronization of a massive array of weakly coupled oscillators; the majority of subcellular processes are implicated in this integrated expression of cellular physiology. Evidence for this view comes mainly from studies of Saccharomyces cerevisiae growing in self-synchronized continuous cultures, in which a temperature-compensated ultradian clock (period of approximately 40 min) couples fermentation with redox state in addition to the transcriptome and cell-division-cycle progression. Functions for ultradian clocks have also been determined in other yeasts (e.g. Schizosaccharomyces pombe and Candida utilis), seven protists (e.g. Acanthamoeba castellanii and Paramecium tetraurelia), as well as cultured mammalian cells. We suggest that ultradian timekeeping is a basic universal necessity for coordinated intracellular coherence.

Deletion ofPDE2, the gene encoding the high-affinity cAMP phosphodiesterase, results in changes of the cell wall and membrane inCandida albicans

A role for the cAMP-dependent pathway in regulation of the cell wall in the model yeast Saccharomyces cerevisiae has recently been demonstrated. In this study we report the results of a phenotypic analysis of a Candida albicans mutant, characterized by a constitutive activation of the cAMP pathway due to deletion of PDE2, the gene encoding the high cAMP-affinity phosphodiesterase. Unlike wild-type strains, this mutant has an increased sensitivity to cell wall and membrane perturbing agents such as SDS and CFW, and antifungals such as amphotericin B and flucytosine. Moreover, the mutant is characterized by an altered sensitivity and a significantly reduced tolerance to fluconazole. The mutant's membrane has around 30% higher ergosterol content and the cell wall glucan was 22% lower than in the wild-type. These cell wall and membrane changes are manifested by a considerable reduction in the thickness of the cell wall, which in the mutant is on average 60-65 nm, compared to 80-85 nm in the wild-type strains as revealed by electron microscopy. These results suggest that constitutive activation of the cAMP pathway affects cell wall and membrane structure, and biosynthesis, not only in the model yeast S. cerevisiae but also in the human fungal pathogen C. albicans.

What makes biochemical networks tick?

In view of the increasing number of reported concentration oscillations in living cells, methods are needed that can identify the causes of these oscillations. These causes always derive from the influences that concentrations have on reaction rates. The influences reach over many molecular reaction steps and are defined by the detailed molecular topology of the network. So-called 'autoinfluence paths', which quantify the influence of one molecular species upon itself through a particular path through the network, can have positive or negative values. The former bring a tendency towards instability. In this molecular context a new graphical approach is presented that enables the classification of network topologies into oscillophoretic and nonoscillophoretic, i.e. into ones that can and ones that cannot induce concentration oscillations. The network topologies are formulated in terms of a set of uni-molecular and bi-molecular reactions, organized into branched cycles of directed reactions, and presented as graphs. Subgraphs of the network topologies are then classified as negative ones (which can) and positive ones (which cannot) give rise to oscillations. A subgraph is oscillophoretic (negative) when it contains more positive than negative autoinfluence paths. Whether the former generates oscillations depends on the values of the other subgraphs, which again depend on the kinetic parameters. An example shows how this can be established. By following the rules of our new approach, various oscillatory kinetic models can be constructed and analyzed, starting from the classified simplest topologies and then working towards desirable complications. Realistic biochemical examples are analyzed with the new method, illustrating two new main classes of oscillophore topologies.

Oscillatory metabolism of Saccharomyces cerevisiae: an overview of mechanisms and models

The budding yeast Saccharomyces cerevisiae displays steady oscillations in continuous cultures under certain conditions. Oscillatory responses are important both metabolically and in process applications. Although much information has become available, a definitive theory to explain and model these oscillations is yet to be formulated. Models of oscillatory cultivation have focussed primarily either on intracellular reactions or on transport processes coupled to substantially lumped intracellular kinetics. This review discusses the development of the models and the directions they provide for a comprehensive model of oscillatory metabolism.

Generation and maintenance of synchrony in Saccharomyces cerevisiae continuous culture

Cultures of Saccharomyces cerevisiae grown continuously produce an autonomous oscillation in many metabolic outputs. The most conveniently measured variable, i.e., dissolved oxygen concentration, oscillates with a period of 40-55 min. Previously we have identified two compounds capable of resetting phase, acetaldehyde and hydrogen sulfide. The phase-response curves constructed for acetaldehyde show a strong (Type 0) response at 3.0 mM and a weak (Type 1) response at 1.0 mM. Ammonium sulfide phase-response curves (pulse injected at 1.0 microM and 3.0 microM) revealed that sulfide is only an effective perturbation agent when endogenous sulfide concentrations are at a maximum. Also only Type 1 phase responses were observed. When the phase-response curve for sulfite (at 3.0 M) was constructed, phase responses were at a maximum at 60 degrees, indicating the possible involvement of sulfite in cell synchronization. It is concluded that endogenously produced acetaldehyde and sulfite tune the oscillation of mitochondrial energization state whereas sulfide mediates population synchrony.

Respiratory oscillations in yeast: mitochondrial reactive oxygen species, apoptosis and time; a hypothesis

Oscillatory metabolic activities occur more widely than is generally realised; detectability requires observation over extended times of single yeast cells or synchrony of individuals to provide a coherent population. Where oscillations in intracellular metabolite concentrations are observed, the phenomenon has been ascribed to sloppy control, energetic optimisation, signalling, temporal compartmentation of incompatible reactions, or timekeeping functions. Here we emphasise the consequences of respiratory oscillations as a source of mitochondrially generated reactive O(2) metabolites. Temporal co-ordination of intracellular activities necessitates a time base. This is provided by an ultradian clock, and one result of its long-term operation is cyclic energisation of mitochondria, and thereby the generation of deleterious free radical species. Our hypothesis is that unrepaired cellular constituents and components (especially mitochondria) eventually lead to cellular senescence and apoptosis when a finite number of respiratory cycles has occurred.

The Gts1 protein stabilizes the autonomous oscillator in yeast

Continuous cultures of Saccharomyces cerevisiae show a robust autonomous temperature compensated oscillation in many metabolic functions. Respiratory activity, a convenient output to measure, oscillates with a period of 40 min. Deletion of GTS1, whose protein product has homology to the circadian per protein, has been implicated in temporal events within yeast, causes a reduction in periodicity to 18 min (wild-type period 40-60 min). The dilution rate was steadily increased from 0.04/h to 0.085/h and the oscillation stabilized after four to six dilutions. However, Gts1p's involvement in the maintenance and generation of metabolic synchrony, and in the central oscillating loop, appear to be minimal, as the mutant oscillation was robust and autonomous. Deletion of GTS1 did cause decreased temperature compensation of the period of the oscillation from Q(10) = 1.07 for the wild-type to Q(10) = 1.6 for the mutant. Also the degree of nutrient compensation observed for the wild-type was not observed in the GTS1-null mutant strain. It is postulated that Gts1p is involved in the mechanism that communicates external conditions, such as temperature, to the central oscillating loop.

Decomposition Analyses Applied to a Complex Ultradian Biorhythm: The Oscillating NADH Oxidase Activity of Plasma Membranes Having a Potential Time-Keeping (Clock) Function

Seasonal decomposition analyses were applied to the statistical evaluation of an oscillating activity for a plasma membrane NADH oxidase activity with a temperature compensated period of 24 min. The decomposition fits were used to validate the cyclic oscillatory pattern. Three measured values, average percentage error (MAPE), a measure of the periodic oscillation, mean average deviation (MAD), a measure of the absolute average deviations from the fitted values, and mean standard deviation (MSD), the measure of standard deviation from the fitted values plus R-squared and the Henriksson-Merton p value were used to evaluate accuracy.Decomposition was carried out by fitting a trend line to the data, then detrending the data if necessary, by subtracting the trend component. The data, with or without detrending, were then smoothed by subtracting a centered moving average of length equal to the period length determined by Fourier analysis. Finally, the time series were decomposed into cyclic and error components. The findings not only validate the periodic nature of the major oscillations but suggest, as well, that the minor intervening fluctuations also recur within each period with a reproducible pattern of recurrence.