Free intracellular amino acid pools during autonomous oscillations in Saccharomyces cerevisiae

Cellular Labeling of Phosphatidylserine Using Clickable Serine Probes

Phosphatidylserine (PS) is a key lipid that plays important roles in disease-related biological processes, and therefore, the means to track PS in live cells are invaluable. Herein, we describe the metabolic labeling of PS in Saccharomyces cerevisiae cells using analogues of serine, a PS precursor, derivatized with azide moieties at either the amino (N-l-SerN3) or carbonyl (C-l-SerN3) groups. The conservative click tag modification enabled these compounds to infiltrate normal lipid biosynthetic pathways, thereby producing tagged PS molecules as supported by mass spectrometry studies, thin-layer chromatography (TLC) analysis, and further derivatization with fluorescent reporters via click chemistry to enable imaging in yeast cells. This approach shows strong prospects for elucidating the complex biosynthetic and trafficking pathways involving PS.

High-level Production of Isoleucine and Fusel alcohol by expression of the Feedback Inhibition-insensitive Threonine deaminase in Saccharomyces cerevisiae

A variety of the yeast Saccharomyces cerevisiae with intracellular accumulation of isoleucine (Ile) would be a promising strain for developing a distinct kind of sake, a traditional Japanese alcoholic beverage, because Ile-derived volatile compounds have a great impact on the flavor and taste of fermented foods. In this study, we isolated an Ile-accumulating mutant (strain K9-I48) derived from a diploid sake yeast of S. cerevisiae by conventional mutagenesis. Strain K9-I48 carries a novel mutation in the ILV1 gene encoding the His480Tyr variant of threonine deaminase (TD). Interestingly, the TD activity of the His480Tyr variant was markedly insensitive to feedback inhibition by Ile, but was not upregulated by valine, leading to intracellular accumulation of Ile and extracellular overproduction of 2-methyl-1-butanol, a fusel alcohol derived from Ile, in yeast cells. The present study demonstrated for the first time that the conserved histidine residue located in a linker region between two regulatory domains is involved in allosteric regulation of TD. Moreover, sake brewed with strain K9-I48 contained 2 to 3 times more 2-methyl-1-butanol and 2-methylbutyl acetate than sake brewed with the parent strain. These findings are valuable for the engineering of TD to increase the productivity of Ile and its derived fusel alcohols.

IMPORTANCE Fruit-like flavors of isoleucine-derived volatile compounds, 2-methyl-1-butanol (2MB) and its acetate ester, contribute to a variety of the flavors and tastes of alcoholic beverages. Besides its value as aroma components in foods and cosmetics, 2MB has attracted significant attention as second-generation biofuels. Threonine deaminase (TD) catalyzes the first step in isoleucine biosynthesis and its activity is subject to feedback inhibition by isoleucine. Here, we isolated an isoleucine-accumulating sake yeast mutant and identified a mutant gene encoding a novel variant of TD. The variant TD exhibited much less sensitivity to isoleucine, leading to higher production of 2MB as well as isoleucine than the wild-type TD. Furthermore, sake brewed with a mutant yeast expressing the variant TD contained more 2MB and its acetate ester than that brewed with the parent strain. These findings will contribute to the development of superior industrial yeast strains for high-level production of isoleucine and its related fusel alcohols.

Repair or Degrade: the Thermodynamic Dilemma of Cellular Protein Quality-Control

Life is a non-equilibrium phenomenon. Owing to their high free energy content, the macromolecules of life tend to spontaneously react with ambient oxygen and water and turn into more stable inorganic molecules. A similar thermodynamic picture applies to the complex shapes of proteins: While a polypeptide is emerging unfolded from the ribosome, it may spontaneously acquire secondary structures and collapse into its functional native conformation. The spontaneity of this process is evidence that the free energy of the unstructured state is higher than that of the structured native state. Yet, under stress or because of mutations, complex polypeptides may fail to reach their native conformation and form instead thermodynamically stable aggregates devoid of biological activity. Cells have evolved molecular chaperones to actively counteract the misfolding of stress-labile proteins dictated by equilibrium thermodynamics. HSP60, HSP70 and HSP100 can inject energy from ATP hydrolysis into the forceful unfolding of stable misfolded structures in proteins and convert them into unstable intermediates that can collapse into the native state, even under conditions inauspicious for that state. Aggregates and misfolded proteins may also be forcefully unfolded and degraded by chaperone-gated endo-cellular proteases, and in eukaryotes also by chaperone-mediated autophagy, paving the way for their replacement by new, unaltered functional proteins. The greater energy cost of degrading and replacing a polypeptide, with respect to the cost of its chaperone-mediated repair represents a thermodynamic dilemma: some easily repairable proteins are better to be processed by chaperones, while it can be wasteful to uselessly try recover overly compromised molecules, which should instead be degraded and replaced. Evolution has solved this conundrum by creating a host of unfolding chaperones and degradation machines and by tuning their cellular amounts and activity rates.

The Metabolic Flux Probe (MFP)-Secreted Protein as a Non-Disruptive Information Carrier for 13C-Based Metabolic Flux Analysis

Novel cultivation technologies demand the adaptation of existing analytical concepts. Metabolic flux analysis (MFA) requires stable-isotope labeling of biomass-bound protein as the primary information source. Obtaining the required protein in cultivation set-ups where biomass is inaccessible due to low cell densities and cell immobilization is difficult to date. We developed a non-disruptive analytical concept for ¹³C-based metabolic flux analysis based on secreted protein as an information carrier for isotope mapping in the protein-bound amino acids. This “metabolic flux probe” (MFP) concept was investigated in different cultivation set-ups with a recombinant, protein-secreting yeast strain. The obtained results grant insight into intracellular protein turnover dynamics. Experiments under metabolic but isotopically nonstationary conditions in continuous glucose-limited chemostats at high dilution rates demonstrated faster incorporation of isotope information from labeled glucose into the recombinant reporter protein than in biomass-bound protein. Our results suggest that the reporter protein was polymerized from intracellular amino acid pools with higher turnover rates than biomass-bound protein. The latter aspect might be vital for ¹³C-flux analyses under isotopically nonstationary conditions for analyzing fast metabolic dynamics.

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.

How evolution shapes enzyme selectivity - lessons from aminoacyl-tRNA synthetases and other amino acid utilizing enzymes

Aminoacyl-tRNA synthetases (AARSs) charge tRNA with their cognate amino acids. Many other enzymes use amino acids as substrates, yet discrimination against noncognate amino acids that threaten the accuracy of protein translation is a hallmark of AARSs. Comparing AARSs to these other enzymes allowed us to recognize patterns in molecular recognition and strategies used by evolution for exercising selectivity. Overall, AARSs are 2-3 orders of magnitude more selective than most other amino acid utilizing enzymes. AARSs also reveal the physicochemical limits of molecular discrimination. For example, amino acids smaller by a single methyl moiety present a discrimination ceiling of ~200, while larger ones can be discriminated by up to 10⁵ -fold. In contrast, substrates larger by a hydroxyl group challenge AARS selectivity, due to promiscuous H-bonding with polar active site groups. This 'hydroxyl paradox' is resolved by editing. Indeed, when the physicochemical discrimination limits are reached, post-transfer editing - hydrolysis of tRNAs charged with noncognate amino acids, evolved. The editing site often selectively recognizes the edited noncognate substrate using the very same feature that the synthetic site could not efficiently discriminate against. Finally, the comparison to other enzymes also reveals that the selectivity of AARSs is an explicitly evolved trait, showing some clear examples of how selection acted not only to optimize catalytic efficiency with the target substrate, but also to abolish activity with noncognate threat substrates ('negative selection').

Genome-Scale Analysis of Perturbations in Translation Elongation Based on a Computational Model

Perturbations play an important role both in engineered systems and cellular processes. Thus, understanding their effect on protein synthesis should contribute to all biomedical disciplines. Here we describe the first genome-scale analysis of perturbations in translation-related factors in S. cerevisiae. To this end, we used simulations based on a computational model that takes into consideration the fundamental stochastic and bio-physical nature of translation. We found that the initiation rate has a key role in determining the sensitivity to perturbations. For low initiation rates, the first codons of the coding region dominate the sensitivity, which is highly correlated with the ratio between initiation rate and mean elongation rate (r = −0.95), with the open reading frame (ORF) length (r = 0.6) and with protein abundance (r = 0.45). For high initiation rates (that may rise, for example, due to cellular growth), the sensitivity of a gene is dominated by all internal codons and is correlated with the decoding rate. We found that various central intracellular functions are associated with the sensitivity: for example, both genes that are sensitive and genes that are robust to perturbations are over-represented in the group of genes related to translation regulation; this may suggest that robustness to perturbations is a trait that undergoes evolutionary selection in relation to the function of the encoded protein. We believe that the reported results, due to their quantitative value and genome-wide perspective, should contribute to disciplines such as synthetic biology, functional genomics, comparative genomics and molecular evolution.

Role of phosphatidylserine synthase in shaping the phospholipidome of Candida albicans

Phosphatidylserine (PS) synthase (Cho1p) and the PS decarboxylase enzymes (Psd1p and Psd2p), which synthesize PS and phosphatidylethanolamine (PE), respectively, are crucial for Candida albicans virulence. Mutations that disrupt these enzymes compromise virulence. These enzymes are part of the cytidine diphosphate-diacylglycerol pathway (i.e. de novo pathway) for phospholipid synthesis. Understanding how losses of PS and/or PE synthesis pathways affect the phospholipidome of Candida is important for fully understanding how these enzymes impact virulence. The cho1Δ/Δ and psd1Δ/Δ psd2Δ/Δ mutations cause similar changes in levels of phosphatidic acid, phosphatidylglycerol, phosphatidylinositol and PS. However, only slight changes were seen in PE and phosphatidylcholine (PC). This finding suggests that the alternative mechanism for making PE and PC, the Kennedy pathway, can compensate for loss of the de novo synthesis pathway. Candida albicans Cho1p, the lipid biosynthetic enzyme with the most potential as a drug target, has been biochemically characterized, and analysis of its substrate specificity and kinetics reveal that these are similar to those previously published for Saccharomyces cerevisiae Cho1p.

Diversification of the kinetic properties of yeast NADP-glutamate-dehydrogenase isozymes proceeds independently of their evolutionary origin

In the yeast Saccharomyces cerevisiae, the ScGDH1 and ScGDH3 encoded glutamate dehydrogenases (NADP-GDHs) catalyze the synthesis of glutamate from ammonium and α-ketoglutarate (α-KG). Previous kinetic characterization showed that these enzymes displayed different allosteric properties and respectively high or low rate of α-KG utilization. Accordingly, the coordinated action of ScGdh1 and ScGdh3, regulated balanced α-KG utilization for glutamate biosynthesis under either fermentative or respiratory conditions, safeguarding energy provision. Here, we have addressed the question of whether there is a correlation between the regulation and kinetic properties of the NADP-GDH isozymes present in S. cerevisiae (ScGdh1 and ScGdh3), Kluyveromyces lactis (KlGdh1), and Lachancea kluyveri (LkGdh1) and their evolutionary history. Our results show that the kinetic properties of K. lactis and L. kluyveri single NADP-GDHs are respectively similar to either ScGDH3 or ScGDH1, which arose from the whole genome duplication event of the S. cerevisiae lineage, although, KlGDH1 and LkGDH1 originated from a GDH clade, through an ancient interspecies hybridization event that preceded the divergence between the Saccharomyces clade and the one containing the genera Kluyveromyces, Lachancea, and Eremothecium. Thus, the kinetic properties which determine the NADP-GDHs capacity to utilize α-KG and synthesize glutamate do not correlate with their evolutionary origin.

(1)H NMR metabolomic study of auxotrophic starvation in yeast using Multivariate Curve Resolution-Alternating Least Squares for Pathway Analysis

Disruption of specific metabolic pathways constitutes the mode of action of many known toxicants and it is responsible for the adverse phenotypes associated to human genetic defects. Conversely, many industrial applications rely on metabolic alterations of diverse microorganisms, whereas many therapeutic drugs aim to selectively disrupt pathogens’ metabolism. In this work we analyzed metabolic changes induced by auxotrophic starvation conditions in yeast in a non-targeted approach, using one-dimensional proton Nuclear Magnetic Resonance spectroscopy (¹H NMR) and chemometric analyses. Analysis of the raw spectral datasets showed specific changes linked to the different stages during unrestricted yeast growth, as well as specific changes linked to each of the four tested starvation conditions (L-methionine, L-histidine, L-leucine and uracil). Analysis of changes in concentrations of more than 40 metabolites by Multivariate Curve Resolution – Alternating Least Squares (MCR-ALS) showed the normal progression of key metabolites during lag, exponential and stationary unrestricted growth phases, while reflecting the metabolic blockage induced by the starvation conditions. In this case, different metabolic intermediates accumulated over time, allowing identification of the different metabolic pathways specifically affected by each gene disruption. This synergy between NMR metabolomics and molecular biology may have clear implications for both genetic diagnostics and drug development.

The dynamics of cellular energetics during continuous yeast culture

A plethora of data is accumulating from high throughput methods on metabolites, coenzymes, proteins, and nucleic acids and their interactions as well as the signalling and regulatory functions and pathways of the cellular network. The frozen moment viewed in a single discrete time sample requires frequent repetition and updating before any appreciation of the dynamics of component interaction becomes possible. Even then in a sample derived from a cell population, time-averaging of processes and events that occur in out-of-phase individuals blur the detailed complexity of single cell organization. Continuously-grown cultures of yeast can become spontaneously self-synchronized, thereby enabling resolution of far more detailed temporal structure. Continuous on-line monitoring by rapidly responding sensors (O2 electrode and membrane-inlet mass spectrometry for O2, CO2 and H2S; direct fluorimetry for NAD(P)H and flavins) gives dynamic information from time-scales of minutes to hours. Supplemented with capillary electophoresis and gas chromatography mass spectrometry and transcriptomics the predominantly oscillatory behaviour of network components becomes evident, with a 40 min cycle between a phase of increased respiration (oxidative phase) and decreased respiration (reductive phase). Highly pervasive, this ultradian clock provides a coordinating function that links mitochondrial energetics and redox balance to transcriptional regulation, mitochondrial structure and organelle remodelling, DNA duplication and cell division events. Ultimately, this leads to a global partitioning of anabolism and catabolism and the enzymes involved, mediated by a relatively simple ATP feedback loop on chromatin architecture.

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.

¹³C-based metabolic flux analysis of recombinant Pichia pastoris

Overexpression of a foreign protein may negatively affect several cell growth parameters, as well as cause cellular stress. Central (or core) metabolism plays a crucial role since it supplies energy, reduction equivalents, and precursor molecules for the recombinant product, cell's maintenance, and growth needs. However, the number of quantitative physiology studies of the impact of recombinant protein production on the central metabolic pathways of yeast cell factories has been traditionally rather limited, thereby hampering the application of rational strain engineering strategies targeting central metabolism.The development and application of quantitative physiology and modelling tools and methodologies is allowing for a systems-level understanding of the effect of bioprocess parameters such as growth rate, temperature, oxygen availability, and substrate(s) choice on metabolism, and its subsequent interactions with recombinant protein synthesis, folding, and secretion.Here, we review the recent developments and applications of (13)C-based metabolic flux analysis ((13)C-MFA) of Pichia pastoris and the gained understanding of the metabolic behavior of this yeast in recombinant protein production bioprocesses. We also discuss the potential of multilevel studies integrating (13)C-MFA with other omics analyses, as well as future perspectives on the metabolic modelling approaches to study and design metabolic engineering strategies for improved protein production.

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.

Sphingoid bases and the serine catabolic enzyme CHA1 define a novel feedforward/feedback mechanism in the response to serine availability

Targets of bioactive sphingolipids in Saccharomyces cerevisiae were previously identified using microarray experiments focused on sphingolipid-dependent responses to heat stress. One of these heat-induced genes is the serine deamidase/dehydratase Cha1 known to be regulated by increased serine availability. This study investigated the hypothesis that sphingolipids may mediate the induction of Cha1 in response to serine availability. The results showed that inhibition of de novo synthesis of sphingolipids, pharmacologically or genetically, prevented the induction of Cha1 in response to increased serine availability. Additional studies implicated the sphingoid bases phytosphingosine and dihydrosphingosine as the likely mediators of Cha1 up-regulation. The yeast protein kinases Pkh1 and Pkh2, known sphingoid base effectors, were found to mediate CHA1 up-regulation via the transcription factor Cha4. Because the results disclosed a role for sphingolipids in negative feedback regulation of serine metabolism, we investigated the effects of disrupting this mechanism on sphingolipid levels and on cell growth. Intriguingly, exposure of the cha1Δ strain to high serine resulted in hyperaccumulation of endogenous serine and in turn a significant accumulation of sphingoid bases and ceramides. Under these conditions, the cha1Δ strain displayed a significant growth defect that was sphingolipid-dependent. Together, this work reveals a feedforward/feedback loop whereby the sphingoid bases serve as sensors of serine availability and mediate up-regulation of Cha1 in response to serine availability, which in turn regulates sphingolipid levels by limiting serine accumulation.

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.

Metabolic flux analysis gives an insight on verapamil induced changes in central metabolism of HL-1 cells

Verapamil has been shown to inhibit glucose transport in several cell types. However, the consequences of this inhibition on central metabolism are not well known. In this study we focused on verapamil induced changes in metabolic fluxes in a murine atrial cell line (HL-1 cells). These cells were adapted to serum free conditions and incubated with 4 μM verapamil and [U-¹³C₅] glutamine. Specific extracellular metabolite uptake/production rates together with mass isotopomer fractions in alanine and glutamate were implemented into a metabolic network model to calculate metabolic flux distributions in the central metabolism. Verapamil decreased specific glucose consumption rate and glycolytic activity by 60%. Although the HL-1 cells show Warburg effect with high lactate production, verapamil treated cells completely stopped lactate production after 24 h while maintaining growth comparable to the untreated cells. Calculated fluxes in TCA cycle reactions as well as NADH/FADH₂ production rates were similar in both treated and untreated cells. This was confirmed by measurement of cell respiration. Reduction of lactate production seems to be the consequence of decreased glucose uptake due to verapamil. In case of tumors, this may have two fold effects; firstly depriving cancer cells of substrate for anaerobic glycolysis on which their growth is dependent; secondly changing pH of the tumor environment, as lactate secretion keeps the pH acidic and facilitates tumor growth. The results shown in this study may partly explain recent observations in which verapamil has been proposed to be a potential anticancer agent. Moreover, in biotechnological production using cell lines, verapamil may be used to reduce glucose uptake and lactate secretion thereby increasing protein production without introduction of genetic modifications and application of more complicated fed-batch processes.

Kinetic and chemical mechanisms of homocitrate synthase from Thermus thermophilus

The homocitrate synthase from Thermus thermophilus (TtHCS) is a metal-activated enzyme with either Mg(2+) or Mn(2+) capable of serving as the divalent cation. The enzyme exhibits a sequential kinetic mechanism. The mechanism is steady state ordered with α-ketoglutarate (α-Kg) binding prior to acetyl-CoA (AcCoA) with Mn(2+), whereas it is steady state random with Mg(2+), suggesting a difference in the competence of the E·Mn·α-Kg·AcCoA and E·Mg·α-Kg·AcCoA complexes. The mechanism is supported by product and dead-end inhibition studies. The primary isotope effect obtained with deuterioacetylCoA (AcCoA-d(3)) in the presence of Mg(2+) is unity (value 1.0) at low concentrations of AcCoA, whereas it is 2 at high concentrations of AcCoA. Data suggest the presence of a slow conformational change induced by binding of AcCoA that accompanies deprotonation of the methyl group of AcCoA. The solvent kinetic deuterium isotope effect is also unity at low AcCoA, but is 1.7 at high AcCoA, consistent with the proposed slow conformational change. The maximum rate is pH independent with either Mg(2+) or Mn(2+) as the divalent metal ion, whereas V/K(α-Kg) (with Mn(2+)) decreases at low and high pH giving pK values of about 6.5 and 8.0. Lysine is a competitive inhibitor that binds to the active site of TtHCS, and shares some of the same binding determinants as α-Kg. Lysine binding exhibits negative cooperativity, indicating cross-talk between the two monomers of the TtHCS dimer. Data are discussed in terms of the overall mechanism of TtHCS.

The role of tRNA and ribosome competition in coupling the expression of different mRNAs in Saccharomyces cerevisiae

Protein synthesis translates information from messenger RNAs into functional proteomes. Because of the finite nature of the resources required by the translational machinery, both the overall protein synthesis activity of a cell and activity on individual mRNAs are controlled by the allocation of limiting resources. Upon introduction of heterologous sequences into an organism—for example for the purposes of bioprocessing or synthetic biology—limiting resources may also become overstretched, thus negatively affecting both endogenous and heterologous gene expression. In this study, we present a mean-field model of translation in Saccharomyces cerevisiae for the investigation of two particular translational resources, namely ribosomes and aminoacylated tRNAs. We firstly use comparisons of experiments with heterologous sequences and simulations of the same conditions to calibrate our model, and then analyse the behaviour of the translational system in yeast upon introduction of different types of heterologous sequences. Our main findings are that: competition for ribosomes, rather than tRNAs, limits global translation in this organism; that tRNA aminoacylation levels exert, at most, weak control over translational activity; and that decoding speeds and codon adaptation exert strong control over local (mRNA specific) translation rates.

A quantitative model for mRNA translation in Saccharomyces cerevisiae

Messenger RNA (mRNA) translation is an essential step in eukaryotic gene expression that contributes to the regulation of this process. We describe a deterministic model based on ordinary differential equations that describe mRNA translation in Saccharomyces cerevisiae. This model, which was parameterized using published data, was developed to examine the kinetic behaviour of translation initiation factors in response to amino acid availability. The model predicts that the abundance of the eIF1-eIF3-eIF5 complex increases under amino acid starvation conditions, suggesting a possible auxiliary role for these factors in modulating translation initiation in addition to the known mechanisms involving eIF2. Our analyses of the robustness of the mRNA translation model suggest that individual cells within a randomly generated population are sensitive to external perturbations (such as changes in amino acid availability) through Gcn2 signalling. However, the model predicts that individual cells exhibit robustness against internal perturbations (such as changes in the abundance of translation initiation factors and kinetic parameters). Gcn2 appears to enhance this robustness within the system. These findings suggest a trade-off between the robustness and performance of this biological network. The model also predicts that individual cells exhibit considerable heterogeneity with respect to their absolute translation rates, due to random internal perturbations. Therefore, averaging the kinetic behaviour of cell populations probably obscures the dynamic robustness of individual cells. This highlights the importance of single-cell measurements for evaluating network properties.

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.

Structural basis for L-lysine feedback inhibition of homocitrate synthase

The alpha-aminoadipate pathway of lysine biosynthesis is modulated at the transcriptional and biochemical levels by feedback inhibition. The first enzyme in the alpha-aminoadipate pathway, homocitrate synthase (HCS), is the target of the feedback regulation and is strongly inhibited by l-lysine. Here we report the structure of Schizosaccharomyces pombe HCS (SpHCS) in complex with l-lysine. The structure illustrates that the amino acid directly competes with the substrate 2-oxoglutarate for binding within the active site of HCS. Differential recognition of the substrate and inhibitor is achieved via a switch position within the (alpha/beta)(8) TIM barrel of the enzyme that can distinguish between the C5-carboxylate group of 2-oxoglutarate and the epsilon-ammonium group of l-lysine. In vitro and in vivo assays demonstrate that mutations of the switch residues, which interact with the l-lysine epsilon-ammonium group, abrogate feedback inhibition, as do substitutions of residues within the C-terminal domain that were identified in a previous study of l-lysine-insensitive HCS mutants in Saccharomyces cerevisiae. Together, these results yield new insights into the mechanism of feedback regulation of an enzyme central to lysine biosynthesis.

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.

Amino acid metabolic origin as an evolutionary influence on protein sequence in yeast

The metabolic cycle of Saccharomyces cerevisiae consists of alternating oxidative (respiration) and reductive (glycolysis) energy-yielding reactions. The intracellular concentrations of amino acid precursors generated by these reactions oscillate accordingly, attaining maximal concentration during the middle of their respective yeast metabolic cycle phases. Typically, the amino acids themselves are most abundant at the end of their precursor's phase. We show that this metabolic cycling has likely biased the amino acid composition of proteins across the S. cerevisiae genome. In particular, we observed that the metabolic source of amino acids is the single most important source of variation in the amino acid compositions of functionally related proteins and that this signal appears only in (facultative) organisms using both oxidative and reductive metabolism. Periodically expressed proteins are enriched for amino acids generated in the preceding phase of the metabolic cycle. Proteins expressed during the oxidative phase contain more glycolysis-derived amino acids, whereas proteins expressed during the reductive phase contain more respiration-derived amino acids. Rare amino acids (e.g., tryptophan) are greatly overrepresented or underrepresented, relative to the proteomic average, in periodically expressed proteins, whereas common amino acids vary by a few percent. Genome-wide, we infer that 20,000 to 60,000 residues have been modified by this previously unappreciated pressure. This trend is strongest in ancient proteins, suggesting that oscillating endogenous amino acid availability exerted genome-wide selective pressure on protein sequences across evolutionary time.

Global self-organization of the cellular metabolic structure

BACKGROUND

Over many years, it has been assumed that enzymes work either in an isolated way, or organized in small catalytic groups. Several studies performed using "metabolic networks models" are helping to understand the degree of functional complexity that characterizes enzymatic dynamic systems. In a previous work, we used "dissipative metabolic networks" (DMNs) to show that enzymes can present a self-organized global functional structure, in which several sets of enzymes are always in an active state, whereas the rest of molecular catalytic sets exhibit dynamics of on-off changing states. We suggested that this kind of global metabolic dynamics might be a genuine and universal functional configuration of the cellular metabolic structure, common to all living cells. Later, a different group has shown experimentally that this kind of functional structure does, indeed, exist in several microorganisms.

METHODOLOGY/PRINCIPAL FINDINGS

Here we have analyzed around 2.500.000 different DMNs in order to investigate the underlying mechanism of this dynamic global configuration. The numerical analyses that we have performed show that this global configuration is an emergent property inherent to the cellular metabolic dynamics. Concretely, we have found that the existence of a high number of enzymatic subsystems belonging to the DMNs is the fundamental element for the spontaneous emergence of a functional reactive structure characterized by a metabolic core formed by several sets of enzymes always in an active state.

CONCLUSIONS/SIGNIFICANCE

This self-organized dynamic structure seems to be an intrinsic characteristic of metabolism, common to all living cellular organisms. To better understand cellular functionality, it will be crucial to structurally characterize these enzymatic self-organized global structures.

Suppression of metabolic defects of yeast isocitrate dehydrogenase and aconitase mutants by loss of citrate synthase

Yeast mutants lacking mitochondrial NAD(+)-specific isocitrate dehydrogenase (idhDelta) or aconitase (aco1Delta) were found to share several growth phenotypes as well as patterns of specific protein expression that differed from the parental strain. These shared properties of idhDelta and aco1Delta strains were eliminated or moderated by co-disruption of the CIT1 gene encoding mitochondrial citrate synthase. Gas chromatography/mass spectrometry analyses indicated a particularly dramatic increase in cellular citrate levels in idhDelta and aco1Delta strains, whereas citrate levels were substantially lower in idhDeltacit1Delta and aco1Deltacit1Delta strains. Exogenous addition of citrate to parental strain cultures partially recapitulated effects of high endogenous levels of citrate in idhDelta and aco1Delta strains. Finally, effects of elevated cellular citrate in idhDelta and aco1Delta mutant strains were partially alleviated by addition of iron or by an increase in pH of the growth medium, suggesting that detrimental effects of citrate are due to elevated levels of the ionized form of this metabolite.

Metabolic profiling of Escherichia coli cultivations: evaluation of extraction and metabolite analysis procedures

The quantitative estimation of intracellular metabolite concentrations (metabolic profiling) is a prerequisite for a better understanding of biological processes and thus inevitable for the rational improvement of microbial production strains and process design. Since pool sizes of substrates regulate flux through different enzymes, the accurate determination of intracellular metabolite concentrations is necessary to understand in vivo reaction kinetics. Quantification of intracellular concentrations of glycolytic intermediates in Escherichia coli K12 was achieved by using a novel in situ rapid sampling and quenching procedure. A new extraction procedure using buffered hot water was established. By use of simultaneous multi-substrate feeding with various ratios of glucose, fructose and acetate during continuous cultivations several metabolic states were induced. Metabolic flux analysis and the newly developed metabolic profiling procedure were used to determine in vivo enzyme kinetics as exemplified for fructose 1,6-bisphosphate aldolase and citrate synthase.

Selective Substrate Supply in the Regulation of Yeast de Novo Sphingolipid Synthesis

The heat stress response of Saccharomyces cerevisiae is characterized by transient cell cycle arrest, altered gene expression, degradation of nutrient permeases, trehalose accumulation, and translation initiation of heat shock proteins. Importantly heat stress also induces de novo sphingolipid synthesis upon which many of these subprograms of the heat stress response depend. Despite extensive data addressing the roles for sphingolipids in heat stress, the mechanism(s) by which heat induces sphingolipid synthesis remains unknown. This study was undertaken to determine the events and/or factors required for heat stress-induced sphingolipid synthesis. Data presented indicate that heat does not directly alter the in vitro activity of serine palmitoyltransferase (SPT), the enzyme responsible for initiating de novo sphingolipid synthesis. Moreover deletion of the small peptide Tsc3p, which is thought to maximize SPT activity, specifically reduced production of C(20) sphingolipid species by over 70% but did not significantly decrease overall sphingoid base production. In contrast, the fatty-acid synthase inhibitor cerulenin nearly completely blocked sphingoid base production after heat, indicating a requirement for endogenous fatty acids for heat-mediated sphingoid base synthesis. Consistent with this, genetic studies show that fatty acid import does not contribute to heat-induced de novo synthesis under normal conditions. Interestingly the absence of medium serine also ameliorated heat-induced sphingoid base production, indicating a requirement for exogenous serine for the response, and consistent with this finding, disruption of synthesis of endogenous serine did not affect heat-induced sphingolipid synthesis. Serine uptake assays indicated that heat increased serine uptake from medium by 100% during the first 10 min of heat stress. Moreover treatments that increase serine uptake in the absence of heat including acute medium acidification and glucose treatment also enhanced de novo sphingoid base synthesis equivalent to that induced by heat stress. These data agree with findings from mammalian systems that availability of substrates is a key determinant of flux through sphingolipid synthesis. Moreover data presented here indicate that SPT activity can be driven by several factors that increase serine uptake in the absence of heat. These findings may provide insights into the many systems in which de novo synthesis is increased in the absence of elevated in vitro SPT activity.

Glutamine-dependent NAD+ synthetase. How a two-domain, three-substrate enzyme avoids waste

Glutamine-dependent NAD(+) synthetase, Qns1, utilizes a glutamine aminotransferase domain to supply ammonia for amidation of nicotinic acid adenine dinucleotide (NaAD(+)) to NAD(+). Earlier characterization of Qns1 suggested that glutamine consumption exceeds NAD(+) production by 40%. To explore whether Qns1 is systematically wasteful or whether additional features account for this behavior, we performed a careful kinetic and molecular genetic analysis. In fact, Qns1 possesses remarkable properties to reduce waste. The glutaminase active site is stimulated by NaAD(+) more than 50-fold such that glutamine is not appreciably consumed in the absence of NaAD(+). Glutamine consumption exceeds NAD(+) production over the whole range of glutamine and NaAD(+) substrate concentrations with greatest efficiency occurring at saturation of both substrates. Kinetic data coupled with site-directed mutagenesis of amino acids in the predicted ammonia channel indicate that NaAD(+) stimulates the glutaminase active site in the k(cat) term by a synergistic mechanism that does not require ammonia utilization by the NaAD(+) substrate. Six distinct classes of Qns1 mutants that fall within the glutaminase domain and the synthetase domain selectively inhibit components of the coordinated reaction.

Characterization of the metabolic shift between oxidative and fermentative growth in Saccharomyces cerevisiae by comparative 13C flux analysis

Background

One of the most fascinating properties of the biotechnologically important organism Saccharomyces cerevisiae is its ability to perform simultaneous respiration and fermentation at high growth rate even under fully aerobic conditions. In the present work, this Crabtree effect called phenomenon was investigated in detail by comparative ¹³C metabolic flux analysis of S. cerevisiae growing under purely oxidative, respiro-fermentative and predominantly fermentative conditions.

Results

The metabolic shift from oxidative to fermentative growth was accompanied by complex changes of carbon flux throughout the whole central metabolism. This involved a flux redirection from the pentose phosphate pathway (PPP) towards glycolysis, an increased flux through pyruvate carboxylase, the fermentative pathways and malic enzyme, a flux decrease through the TCA cycle, and a partial relocation of alanine biosynthesis from the mitochondrion to the cytosol. S. cerevisiae exhibited a by-pass of pyruvate dehydrogenase in all physiological regimes. During oxidative growth this by-pass was mainly provided via pyruvate decarboxylase, acetaldehyde dehydrogenase, acetyl-CoA synthase and transport of acetyl-CoA into the mitochondrion. During fermentative growth this route, however, was saturated due to limited enzyme capacity. Under these conditions the cells exhibited high carbon flux through a chain of reactions involving pyruvate carboxylase, the oxaloacetate transporter and malic enzyme. During purely oxidative growth the PPP alone was sufficient to completely supply NADPH for anabolism. During fermentation, it provided only 60 % of the required NADPH.

Conclusion

We conclude that, in order to overcome the limited capacity of pyruvate dehydrogenase, S. cerevisiae possesses different metabolic by-passes to channel carbon into the mitochondrion. This involves the conversion of cytosolic pyruvate either into acetyl CoA or oxaloacetate followed by intercompartmental transport of these metabolites. During oxidative growth mainly the NAD specific isoforms of acetaldehyde dehydrogenase and isocitrate dehydrogenase catalyze the corresponding reactions in S. cerevisiae, whereas NADPH supply under fermentative conditions involves significant contribution of sources other than the PPP such as e. g. NADPH specific acetaldehyde dehydrogenase or isocitrate dehydrogenase.

Metabolic pathways reconstruction by frequency and amplitude response to forced glycolytic oscillations in yeast

The hypothesis that frequency and amplitude response can be used in a complicated metabolic pathway kinetics model for optimal parameter estimation, as speculated by its successful prior usage for a mechanical oscillator and a heterogeneous chemical system, is tested here. Given the complexity of the glycolysis model of yeast chosen, this question is limited to three kinetics parameters of the 87 in the in vitro model developed in the literature. The direct application of the approach, used with the uninformed selection of operating conditions for the oscillation of external glucose concentration, led to miring the data assimilation process in local minima. Application of linear systems theory, however, identified two natural resonant frequencies that, when excited by external forced oscillations of the same frequency, result in the expression of many harmonics in the Fourier spectra, that is, information-rich experiments. A single such information-rich experiment at one of the resonant frequencies was sufficient to break away from the local minima to find the optimum kinetics parameter estimates. The resonant frequencies themselves represent oscillation modes in glycolysis akin to those previously observed. Furthermore, operation of the bioreactor with large amplitude oscillations of glucose feed (25%) leads to enhanced ethanol average yield by 1.6% at the resonant frequency.

Regulatory mechanism of histidine-tagged homocitrate synthase from Saccharomyces cerevisiae. I. Kinetic studies

Homocitrate synthase (HCS) catalyzes one of the regulated steps of the alpha-aminoadipate pathway for lysine biosynthesis in fungi. The kinetic mechanism of regulation of HCS from Saccharomyces cerevisiae by Na+ and the feedback inhibitor lysine was studied by measuring the initial rate in the absence and presence of the effectors. The data suggest that Na+ is an activator at low concentrations and an inhibitor at high concentrations and that these effects occur as a result of the monovalent ion binding to two different sites in the free enzyme. Inhibition and activation by Na+ can occur simultaneously, with the net rate of the enzyme determined by Na+/K(iNa+) and Na+/K(act), where K(iNa+) and K(act) are the inhibition and activation constants, respectively. The inhibition by Na+ was eliminated at high concentrations of acetyl-CoA, the second substrate bound, but the activation remained. Fluorescence binding studies indicated that lysine bound with high affinity to its binding site as an inhibitor. The inhibition by lysine was competitive versus alpha-ketoglutarate and linear in the physiological range of lysine concentrations up to 5 mm. The effects of Na+ and lysine were independent of one another. A model is developed for regulation of HCS that takes into account all of the effects discussed above.

Autonomous oscillations in Saccharomyces cerevisiae during batch cultures on trehalose

We report that autonomous oscillations, which usually happen in aerobic glucose-limited continuous cultures of yeast at low dilution rate, were also observed in trehalose discontinuous cultures of Saccharomyces cerevisiae. This unexpected oscillatory behaviour was therefore examined using fast Fourier transformation of online gas measurements. This robust mathematical analysis underlined the existence of two types of oscillation. The first was found to be linked to the cell cycle because (a) the periodicity corresponded to a fraction of the generation time and (b) the oscillations were accompanied by a transient increase in the budding index, mobilization of storage carbohydrates, and fermentative activity. Moreover, these oscillations occurred in a range of specific growth rates between 0.04 and 0.15 h(-1). All these criteria were consistent with the cell-cycle-related metabolic oscillations observed in the same range of growth rates in glucose-limited continuous cultures. The second type were short-period respiratory oscillations, independent of the specific growth rate. Both types of oscillation were found to take place consecutively and/or simultaneously during batch culture on trehalose. In addition, mobilization of intracellular trehalose emerged as a key parameter for the sustainability of these autonomous oscillations as they were no longer observed in a mutant defective in neutral trehalase activity. We propose that batch culture on trehalose may be an excellent device for further investigation of the molecular mechanisms that underlie autonomous oscillations in yeast.

Dynamics of intracellular metabolites of glycolysis and TCA cycle during cell-cycle-related oscillation inSaccharomyces cerevisiae

In the present work LC-MS/MS was applied to measure the concentrations of intermediates of glycolysis and TCA cycle during autonomous, cell-cycle synchronized oscillations in aerobic, glucose-limited chemostat cultures of Saccharomyces cerevisiae. This study complements previously reported oscillations in carbon dioxide production rate, intracellular concentrations of trehalose and various free amino acids, and extracellular acetate and pyruvate in the same culture. Of the glycolytic intermediates, fructose 1,6-bisphosphate, 2- and 3-phosphoglycerate, and phosphoenolpyruvate show the most pronounced oscillatory behavior, the latter three compounds oscillating out of phase with the former. This agrees with previously observed metabolic control by phosphofructokinase and pyruvate kinase. Although individually not clearly oscillating, several intermediates of the TCA cycle, i.e., alpha-ketoglutarate, succinate, fumarate, and malate, exhibited increasing concentration during the cell cycle phase with high carbon flux through glycolysis and TCA cycle. The average mass action ratios of beta-phosphoglucomutase and fumarase agreed well with previously determined in vitro equilibrium constants. Minor differences resulted for phosphoglucose isomerase and enolase. Together with the observed close correlation of the pool sizes of the involved metabolites, this might indicate that, in vivo, these reactions are operating close to equilibrium, whereby care must be taken due to possible differences between in vivo and in vitro conditions. Combining the data with previously determined intracellular amino acid levels from the same culture, a few clear correlations between catabolism and anabolism could be identified: phosphoglycerate/serine and alpha-ketoglutarate/lysine exhibited correlated oscillatory behavior, albeit with different phase shifts. Oscillations in intracellular amino acids might therefore be, at least partly, following oscillations of their anabolic precursors.

Impact of transamination reactions and protein turnover on labeling dynamics in (13)C-labeling experiments

A dynamic model describing carbon atom transitions in the central metabolism of Saccharomyces cerevisiae is used to investigate the influence of transamination reactions and protein turnover on the transient behavior of (13)C-labeling chemostat experiments. The simulations performed suggest that carbon exchange due to transamination and protein turnover can significantly increase the required time needed for metabolites in the TCA cycle to reach isotopic steady state, which is in agreement with published experimental observations. On the other hand, transamination and protein turnover will speed-up the net rate of incorporation of labeled carbon into some free and protein-bound amino acids. The simulation results indicate that the pattern of labeled carbon incorporation into amino acids obtained from biomass hydrolysate shows significant deviation from the commonly assumed first-order kinetics behavior until after three residence times. These observations suggest that greater caution should be used while also pointing to new opportunities in the design and interpretation of (13)C-labeling experiments.

The quest for the mechanisms of life

The genomic revolution, manifested by the sequencing of the complete genome of many organisms, along with technological advances, such as DNA microarrays and developments in high-throughput analysis of proteins, metabolites, and isotopic tracer distribution patterns, challenged the conventional ways in which questions are approached in the biological sciences: (a) rather than examining a small number of genes and/or reactions at any one time;, we can now analyze gene expression and protein activity in the context of systems of interacting genes and gene products; (b) comprehensive analysis of biological systems requires the integration of all cellular fingerprints: genome sequence, maps of gene expression, protein expression, metabolic output, and in vivo enzymatic activity; and (c) collecting, managing, and analyzing comparable data from various cellular profiles requires expertise from several fields that transcend traditional discipline boundaries. While researchers in systems biology have still to address difficult challenges in both experimental and computational arenas, they possess, for the first time, the opportunity to unravel the mechanisms of life. The enormous impact of these discoveries in diverse areas, such as metabolic engineering, strain selection, drug screening and development, bioprocess development, disease prognosis and diagnosis, gene and other medical therapies, is an obvious motivation for pursuing integrated analyses of cellular systems.