Contribution of serine, folate and glycine metabolism to the ATP, NADPH and purine requirements of cancer cells.

Recent observations on cancer cell metabolism indicate increased serine synthesis from glucose as a marker of poor prognosis. We have predicted that a fraction of the synthesized serine is routed to a pathway for ATP production. The pathway is composed by reactions from serine synthesis, one-carbon (folate) metabolism and the glycine cleavage system (SOG pathway). Here we show that the SOG pathway is upregulated at the level of gene expression in a subset of human tumors and that its level of expression correlates with gene signatures of cell proliferation and Myc target activation. We have also estimated the SOG pathway metabolic flux in the NCI60 tumor-derived cell lines, using previously reported exchange fluxes and a personalized model of cell metabolism. We find that the estimated rates of reactions in the SOG pathway are highly correlated with the proliferation rates of these cell lines. We also observe that the SOG pathway contributes significantly to the energy requirements of biosynthesis, to the NADPH requirement for fatty acid synthesis and to the synthesis of purines. Finally, when the PC-3 prostate cancer cell line is treated with the antifolate methotrexate, we observe a decrease in the ATP levels, AMP kinase activation and a decrease in ribonucleotides and fatty acids synthesized from [1,2-(13)C2]-D-glucose as the single tracer. Taken together our results indicate that the SOG pathway activity increases with the rate of cell proliferation and it contributes to the biosynthetic requirements of purines, ATP and NADPH of cancer cells.

Identification of novel bacterial histidine biosynthesis inhibitors using docking, ensemble rescoring, and whole-cell assays.

The rapid spread on multidrug-resistant strains of Staphylococcus aureus requires not just novel treatment options, but the development of faster methods for the identification of new hits for drug development. The exponentially increasing speed of computational methods makes a more extensive use in the early stages of drug discovery attractive if sufficient accuracy can be achieved. Computational target identification using systems-level methods suggested the histidine biosynthesis pathway as an attractive target against S. aureus. Potential inhibitors for the pathway were identified through docking, followed by ensemble rescoring, that is sufficiently accurate to justify immediate testing of the identified compounds by whole-cell assays, avoiding the need for time-consuming and often difficult intermediary enzyme assays. This novel strategy is demonstrated for three key enzymes of the S. aureus histidine biosynthesis pathway, which is predicted to be essential for bacterial biomass productions. Virtual screening of a library of approximately 10(6) compounds identified 49 potential inhibitors of three enzymes of this pathway. Eighteen representative compounds were directly tested on three S. aureus- and two Escherichia coli strains in standard disk inhibition assays. Thirteen compounds are inhibitors of some or all of the S. aureus strains, while 14 compounds weakly inhibit growth in one or both E. coli strains. The high hit rate obtained from a fast virtual screen demonstrates the applicability of this novel strategy to the histidine biosynthesis pathway.

Blueprint for antimicrobial hit discovery targeting metabolic networks.

Advances in genome analysis, network biology, and computational chemistry have the potential to revolutionize drug discovery by combining system-level identification of drug targets with the atomistic modeling of small molecules capable of modulating their activity. To demonstrate the effectiveness of such a discovery pipeline, we deduced common antibiotic targets in Escherichia coli and Staphylococcus aureus by identifying shared tissue-specific or uniformly essential metabolic reactions in their metabolic networks. We then predicted through virtual screening dozens of potential inhibitors for several enzymes of these reactions and showed experimentally that a subset of these inhibited both enzyme activities in vitro and bacterial cell viability. This blueprint is applicable for any sequenced organism with high-quality metabolic reconstruction and suggests a general strategy for strain-specific antiinfective therapy.

The implications of human metabolic network topology for disease comorbidity.

Most diseases are the consequence of the breakdown of cellular processes, but the relationships among genetic/epigenetic defects, the molecular interaction networks underlying them, and the disease phenotypes remain poorly understood. To gain insights into such relationships, here we constructed a bipartite human disease association network in which nodes are diseases and two diseases are linked if mutated enzymes associated with them catalyze adjacent metabolic reactions. We find that connected disease pairs display higher correlated reaction flux rate, corresponding enzyme-encoding gene coexpression, and higher comorbidity than those that have no metabolic link between them. Furthermore, the more connected a disease is to other diseases, the higher is its prevalence and associated mortality rate. The network topology-based approach also helps to uncover potential mechanisms that contribute to their shared pathophysiology. Thus, the structure and modeled function of the human metabolic network can provide insights into disease comorbidity, with potentially important consequences for disease diagnosis and prevention.

Intracellular crowding defines the mode and sequence of substrate uptake by Escherichia coli and constrains its metabolic activity.

The influence of the high intracellular concentration of macromolecules on cell physiology is increasingly appreciated, but its impact on system-level cellular functions remains poorly quantified. To assess its potential effect, here we develop a flux balance model of Escherichia coli cell metabolism that takes into account a systems-level constraint for the concentration of enzymes catalyzing the various metabolic reactions in the crowded cytoplasm. We demonstrate that the model's predictions for the relative maximum growth rate of wild-type and mutant E. coli cells in single substrate-limited media, and the sequence and mode of substrate uptake and utilization from a complex medium are in good agreement with subsequent experimental observations. These results suggest that molecular crowding represents a bound on the achievable functional states of a metabolic network, and they indicate that models incorporating this constraint can systematically identify alterations in cellular metabolism activated in response to environmental change.

Topological units of environmental signal processing in the transcriptional regulatory network of Escherichia coli.

Recent evidence indicates that potential interactions within metabolic, protein-protein interaction, and transcriptional regulatory networks are used differentially according to the environmental conditions in which a cell exists. However, the topological units underlying such differential utilization are not understood. Here we use the transcriptional regulatory network of Escherichia coli to identify such units, called origons, representing regulatory subnetworks that originate at a distinct class of sensor transcription factors. Using microarray data, we find that specific environmental signals affect mRNA expression levels significantly only within the origons responsible for their detection and processing. We also show that small regulatory interaction patterns, called subgraphs and motifs, occupy distinct positions in and between origons, offering insights into their dynamical role in information processing. The identified features are likely to represent a general framework for environmental signal processing in prokaryotes.

The topological relationship between the large-scale attributes and local interaction patterns of complex networks.

Recent evidence indicates that the abundance of recurring elementary interaction patterns in complex networks, often called subgraphs or motifs, carry significant information about their function and overall organization. Yet, the underlying reasons for the variable quantity of different subgraph types, their propensity to form clusters, and their relationship with the networks' global organization remain poorly understood. Here we show that a network's large-scale topological organization and its local subgraph structure mutually define and predict each other, as confirmed by direct measurements in five well studied cellular networks. We also demonstrate the inherent existence of two distinct classes of subgraphs, and show that, in contrast to the low-density type II subgraphs, the highly abundant type I subgraphs cannot exist in isolation but must naturally aggregate into subgraph clusters. The identified topological framework may have important implications for our understanding of the origin and function of subgraphs in all complex networks.

Global organization of metabolic fluxes in the bacterium Escherichia coli.

Cellular metabolism, the integrated interconversion of thousands of metabolic substrates through enzyme-catalysed biochemical reactions, is the most investigated complex intracellular web of molecular interactions. Although the topological organization of individual reactions into metabolic networks is well understood, the principles that govern their global functional use under different growth conditions raise many unanswered questions. By implementing a flux balance analysis of the metabolism of Escherichia coli strain MG1655, here we show that network use is highly uneven. Whereas most metabolic reactions have low fluxes, the overall activity of the metabolism is dominated by several reactions with very high fluxes. E. coli responds to changes in growth conditions by reorganizing the rates of selected fluxes predominantly within this high-flux backbone. This behaviour probably represents a universal feature of metabolic activity in all cells, with potential implications for metabolic engineering.

Experimental determination and system level analysis of essential genes in Escherichia coli MG1655.

Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

Evolutionary conservation of motif constituents in the yeast protein interaction network.

Understanding why some cellular components are conserved across species but others evolve rapidly is a key question of modern biology. Here we show that in Saccharomyces cerevisiae, proteins organized in cohesive patterns of interactions are conserved to a substantially higher degree than those that do not participate in such motifs. We find that the conservation of proteins in distinct topological motifs correlates with the interconnectedness and function of that motif and also depends on the structure of the overall interactome topology. These findings indicate that motifs may represent evolutionary conserved topological units of cellular networks molded in accordance with the specific biological function in which they participate.

Osteopontin contributes to hepatocyte growth factor-induced tumor growth and metastasis formation.

The cytokine hepatocyte growth factor (HGF)/scatter factor-1 and its cognate receptor, Met, are involved in the etiology and progression of many types of cancer. Despite recent advances in understanding the signal transduction pathways activated by HGF, the mechanism by which HGF exerts its tumorigenic effect is not well understood. To identify proteins that may be involved in mediating HGF-induced cell motility, invasiveness, and tumorigenesis, we used two separate differential display screening methods to identify changes in gene expression that are initiated by HGF in an epithelial cell culture system. Among several known and unknown genes whose expression was modified, osteopontin (OPN), a protein previously associated with tumorigenesis, was found to be upregulated within 6 h following HGF stimulation. OPN expression was dependent on activation of the PI-3 kinase pathway. Autocrine secretion of HGF resulted in sustained expression of OPN. Downregulation of opn expression by stable antisense transfection attenuated OPN expression and repressed HGF-induced invasiveness in vitro and decreased HGF-mediated tumor growth and metastasis formation in vivo. Constitutive expression of OPN in itself exerted partial invasiveness in vitro, but its expression itself was not sufficient to initiate tumor growth or metastasis formation in vivo. Thus, together with other molecules, OPN activity contributes to HGF-induced tumor growth and invasiveness.

Hierarchical organization of modularity in metabolic networks.

Spatially or chemically isolated functional modules composed of several cellular components and carrying discrete functions are considered fundamental building blocks of cellular organization, but their presence in highly integrated biochemical networks lacks quantitative support. Here, we show that the metabolic networks of 43 distinct organisms are organized into many small, highly connected topologic modules that combine in a hierarchical manner into larger, less cohesive units, with their number and degree of clustering following a power law. Within Escherichia coli, the uncovered hierarchical modularity closely overlaps with known metabolic functions. The identified network architecture may be generic to system-level cellular organization.

Comparable system-level organization of Archaea and Eukaryotes.

A central and long-standing issue in evolutionary theory is the origin of the biological variation upon which natural selection acts. Some hypotheses suggest that evolutionary change represents an adaptation to the surrounding environment within the constraints of an organism's innate characteristics. Elucidation of the origin and evolutionary relationship of species has been complemented by nucleotide sequence and gene content analyses, with profound implications for recognizing life's major domains. Understanding of evolutionary relationships may be further expanded by comparing systemic higher-level organization among species. Here we employ multivariate analyses to evaluate the biochemical reaction pathways characterizing 43 species. Comparison of the information transfer pathways of Archaea and Eukaryotes indicates a close relationship between these domains. In addition, whereas eukaryotic metabolic enzymes are primarily of bacterial origin, the pathway-level organization of archaeal and eukaryotic metabolic networks is more closely related. Our analyses therefore suggest that during the symbiotic evolution of eukaryotes, incorporation of bacterial metabolic enzymes into the proto-archaeal proteome was constrained by the host's pre-existing metabolic architecture.

The large-scale organization of metabolic networks.

In a cell or microorganism, the processes that generate mass, energy, information transfer and cell-fate specification are seamlessly integrated through a complex network of cellular constituents and reactions. However, despite the key role of these networks in sustaining cellular functions, their large-scale structure is essentially unknown. Here we present a systematic comparative mathematical analysis of the metabolic networks of 43 organisms representing all three domains of life. We show that, despite significant variation in their individual constituents and pathways, these metabolic networks have the same topological scaling properties and show striking similarities to the inherent organization of complex non-biological systems. This may indicate that metabolic organization is not only identical for all living organisms, but also complies with the design principles of robust and error-tolerant scale-free networks, and may represent a common blueprint for the large-scale organization of interactions among all cellular constituents.

Interleukin-3 induces the phosphorylation of a distinct fraction of bcl-2.

Bcl-2-related proteins (i.e. Bcl-2 and Bax) regulate the effector stage of apoptosis and can modulate the entry of quiescent cells into the cell cycle. Phosphorylation of Bcl-2 is presumed to modify its apoptosis-inhibitory function. By utilizing an interleukin-3 (IL-3)-dependent hematopoietic cell line, we examined the structural requirements of Bcl-2 phosphorylation and the correlation of this post-translational modification with its function. In the presence of IL-3, constitutively expressed Bcl-2 was phosphorylated on serine residue(s), and phosphorylated Bcl-2 lost its capacity to heterodimerize with Bax. Whereas the majority of Bcl-2 resided in mitochondria, phosphorylation only affected a minor pool of total Bcl-2 that selectively partitioned into a soluble fraction. Cytosolic targeting of Bcl-2 greatly increased its ratio of phosphorylation. Bcl-2 phosphorylation was reduced during IL-3 deprivation, and its phosphorylation was also delayed after transient cytokine deprivation. This pattern of phosphorylation temporally correlated with the accelerated exit and delayed reentry of Bcl-2-expressing cells into the cell cycle upon transient IL-3 deprivation and subsequent cytokine restimulation. Thus, IL-3-induced phosphorylation of a distinct pool of Bcl-2 may contribute to the inactivation of its antiproliferative function.

Bcl-2 counters apoptosis by Bax heterodimerization-dependent and -independent mechanisms in the T-cell lineage.

The effect of the cell death inhibitor Bcl-2 in relation to its capacity to dimerize with apoptosis promoter Bax or its homologs at their physiological expression levels was explored in the T-cell lineage. Transgenic mice expressing a BH1 mutant Bcl-2 (Bcl-2 mI-3), which fails to heterodimerize with proapoptotic members of the Bcl-2 family, such as Bax, were generated. Bcl-2 mI-3 protected immature CD4+8- thymocytes from spontaneous, glucocorticoid and anti-CD3-induced apoptosis and altered T cell maturation, resulting in increased percentages of CD3(hi) and CD4-8+ thymocytes. In contrast, apoptosis of peripheral T-cells was unaffected by transgene expression. This correlated with their high Bax expression level and insensitivity to the caspase inhibitor, zVAD-fmk, a functional hallmark of Bax-like activity. Thus, within the T-cell lineage Bcl-2 can inhibit apoptosis independent of its association with Bax or its homologs; yet, above a threshold level of their physiologic proapoptotic activity, the capacity of Bcl-2 to heterodimerize with Bax or its homologs appears essential for it to counter cell death.

The ability of Bcl-x(L) and Bcl-2 to prevent apoptosis can be differentially regulated.

Although expression of Bcl-2 has been shown to prevent apoptosis under many circumstances, there are several systems in which Bcl-2 fails to promote cell survival. We have previously reported that Bcl-2 and Bcl-x(L) display differential ability to protect WEHI-231 cells from multiple inducers of apoptosis. A possible explanation for this paradox was provided by the discovery of Bax. Bax is a Bcl-2-related protein which can inhibit the ability of Bcl-2 to enhance the survival of growth factor-dependent cell lines in the absence of growth factor. Consistent with the possibility that Bcl-2 function in WEHI-231 cells is inhibited by Bax, WEHI-231 cells were found to express a high level of Bax. To directly test the effects of Bax expression on Bcl-x(L) function, FL5.12 cells were transfected with both genes. Although Bax overexpression can inhibit Bcl-2 from prolonging cell survival upon growth factor withdrawal, Bax overexpression did not inhibit Bcl-x(L) from preventing apoptosis in this cell line. Although Bcl-2 and Bcl-x(L) were both found to be able to form heterodimers with Bax, the majority of Bax in both cases was not complexed to a partner. Our data suggest that Bcl-x(L) does not function by simply preventing the formation of Bax homodimers which promote cell death. Instead Bax appears to display selectivity in its ability to inhibit Bcl-2 but not Bcl-x(L) from prolonging survival. Furthermore, our data suggest that the abilities of Bcl-2 and Bcl-x(L) to promote cell survival are not identical and can be independently regulated within a cell. Regulation of a cell's apoptotic threshold is likely to result from a complex set of interactions among Bcl-2 family members and other, as yet uncharacterised, regulators of apoptosis.

Multiple Bcl-2 family members demonstrate selective dimerizations with Bax.

A family of Bcl-2-related proteins regulates cell death and shares highly conserved BH1 and BH2 domains. BH1 and BH2 domains of Bcl-2 were required for it to heterodimerize with Bax and to repress apoptosis. A yeast two-hybrid assay accurately reproduced this interaction and defined a selectivity and hierarchy of further dimerizations. Bax also heterodimerizes with Bcl-xL, Mcl-1, and A1. A Gly-159-->Ala substitution in BH1 of Bcl-xL disrupted its heterodimerization with Bax and abrogated its inhibition of apoptosis in mammalian cells. This suggests that the susceptibility to apoptosis is determined by multiple competing dimerizations in which Bax may be a common partner.

Reactive oxygen species and the regulation of cell death by the Bcl-2 gene family.

The maintenance of homeostasis in normal tissues reflects a balance between cell proliferation and cell death. Bcl-2 inaugurated a new category of oncogenes, regulators of cell death. The Bcl-2 gene was identified at the chromosomal breakpoint of t(14;18) bearing B cell lymphomas. Bcl-2 proved unique by blocking programmed cell death rather than promoting proliferation. In adults, Bcl-2 is topographically restricted to progenitor cells and longlived cells but is much more widespread in the developing embryo. Transgenic mice that overexpress Bcl-2 demonstrate extended cell survival, and progress to high grade lymphomas. Bcl-2 has been localized to mitochondria, endoplasmic reticulum and nuclear membranes, also the sites of reactive oxygen species generation. Bcl-2 does not appear to influence the generation of oxygen free radicals but does prevent oxidative damage to cellular constituents including lipid membranes. Bcl-2 deficient mice complete embryonic development but undergo fulminant lymphoid apoptosis of thymus and spleen. Moreover, they demonstrate two unexpected pathologies resulting from cell death, polycystic kidney disease and hair hypopigmentation. The latter is a potential oxidant injury from the melanin biosynthetic pathway. A family of Bcl-2 related genes is emerging that includes Bax, a conserved homolog that heterodimerizes in vivo with Bcl-2 and promotes cell death. The ratio of family members, such as Bcl-2/Bax, determines the survival or death of cells following an apoptotic stimulus.