Understanding human metabolic physiology: a genome-to-systems approach

Manifold fitting reveals metabolomic heterogeneity and disease associations in UK Biobank populations

NMR-based metabolic biomarkers provide comprehensive insights into human metabolism; however, extracting biologically meaningful patterns from such high-dimensional data remains a significant challenge. In this study, we propose a manifold-fitting-based framework to analyze metabolic heterogeneity within the UK Biobank population, utilizing measurements of 251 NMR biomarkers from 212,853 participants. Initially, our method clusters these biomarkers into seven distinct metabolic categories that reflect the modular organization of human metabolism. Subsequent manifold fitting to each category unveils underlying low-dimensional structures, elucidating fundamental variations from basic energy metabolism to hormone-mediated regulation. Importantly, three of these manifolds clearly stratify the population, identifying subgroups with distinct metabolic profiles and associated disease risks. These subgroups exhibit consistent links with specific diseases, including severe metabolic dysregulation and its complications, as well as cardiovascular and autoimmune conditions, highlighting the intricate relationship between metabolic states and disease susceptibility. Supported by strong correlations with demographic factors, clinical measurements, and lifestyle variables, these findings validate the biological relevance of the identified manifolds. By utilizing a geometrically informed approach to dissect metabolic heterogeneity, our framework enhances the accuracy of population stratification and deepens our understanding of metabolic health, potentially guiding personalized interventions and preventive healthcare strategies.

MetaboRank: network-based recommendation system to interpret and enrich metabolomics results

Motivation

Metabolomics has shown great potential to improve the understanding of complex diseases, potentially leading to therapeutic target identification. However, no single analytical method allows monitoring all metabolites in a sample, resulting in incomplete metabolic fingerprints. This incompleteness constitutes a stumbling block to interpretation, raising the need for methods that can enrich those fingerprints. We propose MetaboRank, a new solution inspired by social network recommendation systems for the identification of metabolites potentially related to a metabolic fingerprint.

Results

MetaboRank method had been used to enrich metabolomics data obtained on cerebrospinal fluid samples from patients suffering from hepatic encephalopathy (HE). MetaboRank successfully recommended metabolites not present in the original fingerprint. The quality of recommendations was evaluated by using literature automatic search, in order to check that recommended metabolites could be related to the disease. Complementary mass spectrometry experiments and raw data analysis were performed to confirm these suggestions. In particular, MetaboRank recommended the overlooked α-ketoglutaramate as a metabolite which should be added to the metabolic fingerprint of HE, thus suggesting that metabolic fingerprints enhancement can provide new insight on complex diseases.

Availability and implementation

Method is implemented in the MetExplore server and is available at www.metexplore.fr. A tutorial is available at https://metexplore.toulouse.inra.fr/com/tutorials/MetaboRank/2017-MetaboRank.pdf.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-scale metabolic modelling common cofactors metabolism in microorganisms

The common cofactors ATP/ADP, NAD(P)(H), and acetyl-CoA/CoA are indispensable participants in biochemical reactions in industrial microbes. To systematically explore the effects of these cofactors on cell growth and metabolic phenotypes, the first genome-scale cofactor metabolic model, icmNX6434, including 6434 genes, 1782 metabolites, and 6877 reactions, was constructed from 14 genome-scale metabolic models of 14 industrial strains. The origin, consumption, and interactions of these common cofactors in microbial cells were elucidated by the icmNX6434 model, and they played important roles in cell growth. The essential cofactor modules contained 2480 genes and 2948 reactions; therefore, improving cofactor biosynthesis, directing these cofactors into essential metabolic pathways, as well as avoiding cofactor utilization during byproduct biosynthesis and futile cycles, are three ways to increase cell growth. The effects of these common cofactors on the distribution and rate of the carbon flux in four universal modes, as well as an optimized metabolic flux, could be obtained by manipulating cofactor availability and balance. Significant changes in the ATP, NAD(H), NADP(H), or acetyl-CoA concentrations triggered relevant metabolic responses to acidic, oxidative, heat, and osmotic stress. Globally, the model icmNX6434 provides a comprehensive platform to elucidate the physiological effects of these cofactors on cell growth, metabolic flux, and industrial robustness. Moreover, the results of this study are a further example of using a consensus genome-scale metabolic model to increase our understanding of key biological processes.

Perspectives and Challenges of Microbial Application for Crop Improvement

Global population increases and climate change pose a challenge to worldwide crop production. There is a need to intensify agricultural production in a sustainable manner and to find solutions to combat abiotic stress, pathogens, and pests. Plants are associated with complex microbiomes, which have an ability to promote plant growth and stress tolerance, support plant nutrition, and antagonize plant pathogens. The integration of beneficial plant-microbe and microbiome interactions may represent a promising sustainable solution to improve agricultural production. The widespread commercial use of the plant beneficial microorganisms will require a number of issues addressed. Systems approach using microscale information technology for microbiome metabolic reconstruction has potential to advance the microbial reproducible application under natural conditions.

Integrative analysis of human omics data using biomolecular networks

High-throughput '-omics' technologies have given rise to an increasing abundance of genome-scale data detailing human biology at the molecular level. Although these datasets have already made substantial contributions to a more comprehensive understanding of human physiology and diseases, their interpretation becomes increasingly cryptic and nontrivial as they continue to expand in size and complexity. Systems biology networks offer a scaffold upon which omics data can be integrated, facilitating the extraction of new and physiologically relevant information from the data. Two of the most prevalent networks that have been used for such integrative analyses of omics data are genome-scale metabolic models (GEMs) and protein-protein interaction (PPI) networks, both of which have demonstrated success among many different omics and sample types. This integrative approach seeks to unite 'top-down' omics data with 'bottom-up' biological networks in a synergistic fashion that draws on the strengths of both strategies. As the volume and resolution of high-throughput omics data continue to grow, integrative network-based analyses are expected to play an increasingly important role in their interpretation.

New strategy for drug discovery by large-scale association analysis of molecular networks of different species

The development of modern omics technology has not significantly improved the efficiency of drug development. Rather precise and targeted drug discovery remains unsolved. Here a large-scale cross-species molecular network association (CSMNA) approach for targeted drug screening from natural sources is presented. The algorithm integrates molecular network omics data from humans and 267 plants and microbes, establishing the biological relationships between them and extracting evolutionarily convergent chemicals. This technique allows the researcher to assess targeted drugs for specific human diseases based on specific plant or microbe pathways. In a perspective validation, connections between the plant Halliwell-Asada (HA) cycle and the human Nrf2-ARE pathway were verified and the manner by which the HA cycle molecules act on the human Nrf2-ARE pathway as antioxidants was determined. This shows the potential applicability of this approach in drug discovery. The current method integrates disparate evolutionary species into chemico-biologically coherent circuits, suggesting a new cross-species omics analysis strategy for rational drug development.

To supplement or not to supplement: a metabolic network framework for human nutritional supplements

Flux balance analysis and constraint based modeling have been successfully used in the past to elucidate the metabolism of single cellular organisms. However, limited work has been done with multicellular organisms and even less with humans. The focus of this paper is to present a novel use of this technique by investigating human nutrition, a challenging field of study. Specifically, we present a steady state constraint based model of skeletal muscle tissue to investigate amino acid supplementation's effect on protein synthesis. We implement several in silico supplementation strategies to study whether amino acid supplementation might be beneficial for increasing muscle contractile protein synthesis. Concurrent with published data on amino acid supplementation's effect on protein synthesis in a post resistance exercise state, our results suggest that increasing bioavailability of methionine, arginine, and the branched-chain amino acids can increase the flux of contractile protein synthesis. The study also suggests that a common commercial supplement, glutamine, is not an effective supplement in the context of increasing protein synthesis and thus, muscle mass. Similar to any study in a model organism, the computational modeling of this research has some limitations. Thus, this paper introduces the prospect of using systems biology as a framework to formally investigate how supplementation and nutrition can affect human metabolism and physiology.

Consensus and conflict cards for metabolic pathway databases

Background

The metabolic network of H. sapiens and many other organisms is described in multiple pathway databases. The level of agreement between these descriptions, however, has proven to be low. We can use these different descriptions to our advantage by identifying conflicting information and combining their knowledge into a single, more accurate, and more complete description. This task is, however, far from trivial.

Results

We introduce the concept of Consensus and Conflict Cards (C₂Cards) to provide concise overviews of what the databases do or do not agree on. Each card is centered at a single gene, EC number or reaction. These three complementary perspectives make it possible to distinguish disagreements on the underlying biology of a metabolic process from differences that can be explained by different decisions on how and in what detail to represent knowledge. As a proof-of-concept, we implemented C₂Cards^Human, as a web application http://www.molgenis.org/c2cards, covering five human pathway databases.

Conclusions

C₂Cards can contribute to ongoing reconciliation efforts by simplifying the identification of consensus and conflicts between pathway databases and lowering the threshold for experts to contribute. Several case studies illustrate the potential of the C₂Cards in identifying disagreements on the underlying biology of a metabolic process. The overviews may also point out controversial biological knowledge that should be subject of further research. Finally, the examples provided emphasize the importance of manual curation and the need for a broad community involvement.

Exploring metabolic pathways that contribute to the stem cell phenotype

BACKGROUND

Stem cells must negotiate their surrounding nutritional and signaling environment and respond accordingly to perform various functions. Metabolic pathways enable these responses, providing energy and biosynthetic precursors for cell proliferation, motility, and other functions. As a result, metabolic enzymes and the molecules which control them are emerging as attractive targets for the manipulation of stem cells. To exploit these targets a detailed characterization of metabolic flux regulation is required.

SCOPE OF REVIEW

Here we outline recent advances in our understanding of metabolism in pluripotent stem cells and adult progenitors. We describe the regulation of glycolysis, mitochondrial metabolism, and the redox state of stem cells, highlighting key enzymes and transcription factors involved in the control of these pathways.

MAJOR CONCLUSIONS

A general description of stem cell metabolism has emerged, involving increased glycolysis, limited oxidative metabolism, and resistance to oxidative damage. Moving forward, the application of systems-based approaches to stem cells will help shed light on metabolic pathway utilization in proliferating and quiescent stem cells.

GENERAL SIGNIFICANCE

Metabolic flux contributes to the unique properties of stem cells and progenitors. This review provides a detailed overview of how stem cells metabolize their surrounding nutrients to proliferate and maintain lineage homeostasis. This article is part of a Special Issue entitled Biochemistry of Stem Cells.

Genome-scale metabolic model of the fission yeast Schizosaccharomyces pombe and the reconciliation of in silico/in vivo mutant growth

Background

Over the last decade, the genome-scale metabolic models have been playing increasingly important roles in elucidating metabolic characteristics of biological systems for a wide range of applications including, but not limited to, system-wide identification of drug targets and production of high value biochemical compounds. However, these genome-scale metabolic models must be able to first predict known in vivo phenotypes before it is applied towards these applications with high confidence. One benchmark for measuring the in silico capability in predicting in vivo phenotypes is the use of single-gene mutant libraries to measure the accuracy of knockout simulations in predicting mutant growth phenotypes.

Results

Here we employed a systematic and iterative process, designated as Reconciling In silico/in vivo mutaNt Growth (RING), to settle discrepancies between in silico prediction and in vivo observations to a newly reconstructed genome-scale metabolic model of the fission yeast, Schizosaccharomyces pombe, SpoMBEL1693. The predictive capabilities of the genome-scale metabolic model in predicting single-gene mutant growth phenotypes were measured against the single-gene mutant library of S. pombe. The use of RING resulted in improving the overall predictive capability of SpoMBEL1693 by 21.5%, from 61.2% to 82.7% (92.5% of the negative predictions matched the observed growth phenotype and 79.7% the positive predictions matched the observed growth phenotype).

Conclusion

This study presents validation and refinement of a newly reconstructed metabolic model of the yeast S. pombe, through improving the metabolic model’s predictive capabilities by reconciling the in silico predicted growth phenotypes of single-gene knockout mutants, with experimental in vivo growth data.

Ruminant metabolic systems biology: reconstruction and integration of transcriptome dynamics underlying functional responses of tissues to nutrition and physiological state

High-throughput ‘omics’ data analysis via bioinformatics is one key component of the systems biology approach. The systems approach is particularly well-suited for the study of the interactions between nutrition and physiological state with tissue metabolism and functions during key life stages of organisms such as the transition from pregnancy to lactation in mammals, ie, the peripartal period. In modern dairy cows with an unprecedented genetic potential for milk synthesis, the nature of the physiologic and metabolic adaptations during the peripartal period is multifaceted and involves key tissues such as liver, adipose, and mammary. In order to understand such adaptation, we have reviewed several works performed in our and other labs. In addition, we have used a novel bioinformatics approach, Dynamic Impact Approach (DIA), in combination with partly previously published data to help interpret longitudinal biological adaptations of bovine liver, adipose, and mammary tissue to lactation using transcriptomics datasets. Use of DIA with transcriptomic data from those tissues during normal physiological adaptations and in animals fed different levels of energy prepartum allowed visualization and integration of most-impacted metabolic pathways around the time of parturition. The DIA is a suitable tool for applying the integrative systems biology approach. The ultimate goal is to visualize the complexity of the systems at study and uncover key molecular players involved in the tissue’s adaptations to physiological state or nutrition.

Constraints-based stoichiometric analysis of hypoxic stress on steroidogenesis in fathead minnows, Pimephales promelas

In this study, an in silico genome-scale metabolic model of steroidogenesis was used to investigate the effects of hypoxic stress on steroid hormone productions in fish. Adult female fathead minnows (Pimephales promelas) were exposed to hypoxia for 7 days with fish sub-sampled on days 1, 3 and 7 of exposure. At each time point, selected steroid enzyme gene expressions and steroid hormone productions were quantified in ovaries. Fold changes in steroid enzyme gene expressions were used to qualitatively scale transcript enzyme reaction constraints (akin to the range of an enzyme’s catalytic activity) in the in silico model. Subsequently, in silico predicted steroid hormone productions were qualitatively compared with experimental results. Key findings were as follows. (1) In silico gene deletion analysis identified highly conserved ‘essential’ genes required for steroid hormone productions. These agreed well (75%) with literature-published genes downregulated in vertebrates (fish and mammal) exposed to hypoxia. (2) Quantification of steroid hormones produced ex vivo from ovaries showed a significant reduction for 17β-estradiol and 17α,20β-dihydroxypregnenone production after 24 h (day 1) of exposure. This lowered 17β-estradiol production was concomitant with downregulation of cyp19a1a gene expression in ovaries. In silico predictions showed agreement with experimentation by predicting effects on estrogen (17β-estradiol and estrone) production. (3) Stochastic sampling of in silico reactions indicated that cholesterol uptake and catalysis to pregnenolone along with estrogen methyltransferase and glucuronidation reactions were also impacted by hypoxia. Taken together, this in silico analysis introduces a powerful model for pathway analysis that can lend insights on the effects of various stressor scenarios on metabolic functions.

Bioinformatics and systems biology

Delivering personalized therapeutic options to cancer patients based on the genetic and molecular aberrations of the tumor offers great promise to improve the outcomes of cancer therapy. Significant progress in biotechnology has allowed the measurement of tens of thousands of "omic" data points across multiple levels (DNA, RNA protein, metabolomics) from a single tumor biopsy sample in a reasonable time frame for making clinical decisions. With this data in hand, the challenge from the bioinformatics and systems biology point of view is how does one convert data into information and knowledge that can improve the delivery of personalized therapy to the patient.

Systems medicine and metabolic modelling

Several complex diseases are caused by the malfunction of human metabolism, and deciphering the underlying molecular mechanisms can elucidate their aetiology. Systems biology is an integrative approach combining experimental and computational biology to identify and describe the molecular mechanisms of complex biological systems. Systems medicine has the potential to elucidate the onset and progression of complex metabolic diseases through the use of computational approaches. Advances in biotechnology have resulted in the provision of high-throughput data, which provide information about different metabolic processes. The systems medicine approach can utilize such data to reconstruct genome-scale metabolic models that can be used to study the function of specific enzymes and pathways in the context of the complete metabolic network. In this review, we outline the importance of genome-scale models in systems medicine and discuss how they may contribute towards the development of personalized medicine.

Critical assessment of human metabolic pathway databases: a stepping stone for future integration

Background

Multiple pathway databases are available that describe the human metabolic network and have proven their usefulness in many applications, ranging from the analysis and interpretation of high-throughput data to their use as a reference repository. However, so far the various human metabolic networks described by these databases have not been systematically compared and contrasted, nor has the extent to which they differ been quantified. For a researcher using these databases for particular analyses of human metabolism, it is crucial to know the extent of the differences in content and their underlying causes. Moreover, the outcomes of such a comparison are important for ongoing integration efforts.

Results

We compared the genes, EC numbers and reactions of five frequently used human metabolic pathway databases. The overlap is surprisingly low, especially on reaction level, where the databases agree on 3% of the 6968 reactions they have combined. Even for the well-established tricarboxylic acid cycle the databases agree on only 5 out of the 30 reactions in total. We identified the main causes for the lack of overlap. Importantly, the databases are partly complementary. Other explanations include the number of steps a conversion is described in and the number of possible alternative substrates listed. Missing metabolite identifiers and ambiguous names for metabolites also affect the comparison.

Conclusions

Our results show that each of the five networks compared provides us with a valuable piece of the puzzle of the complete reconstruction of the human metabolic network. To enable integration of the networks, next to a need for standardizing the metabolite names and identifiers, the conceptual differences between the databases should be resolved. Considerable manual intervention is required to reach the ultimate goal of a unified and biologically accurate model for studying the systems biology of human metabolism. Our comparison provides a stepping stone for such an endeavor.

Interactome networks and human disease

Complex biological systems and cellular networks may underlie most genotype to phenotype relationships. Here, we review basic concepts in network biology, discussing different types of interactome networks and the insights that can come from analyzing them. We elaborate on why interactome networks are important to consider in biology, how they can be mapped and integrated with each other, what global properties are starting to emerge from interactome network models, and how these properties may relate to human disease.

Synthetic biology: putting synthesis into biology

The ability to manipulate living organisms is at the heart of a range of emerging technologies that serve to address important and current problems in environment, energy, and health. However, with all its complexity and interconnectivity, biology has for many years been recalcitrant to engineering manipulations. The recent advances in synthesis, analysis, and modeling methods have finally provided the tools necessary to manipulate living systems in meaningful ways and have led to the coining of a field named synthetic biology. The scope of synthetic biology is as complicated as life itself--encompassing many branches of science and across many scales of application. New DNA synthesis and assembly techniques have made routine customization of very large DNA molecules. This in turn has allowed the incorporation of multiple genes and pathways. By coupling these with techniques that allow for the modeling and design of protein functions, scientists have now gained the tools to create completely novel biological machineries. Even the ultimate biological machinery--a self-replicating organism--is being pursued at this moment. The aim of this article is to dissect and organize these various components of synthetic biology into a coherent picture.

Reverse engineering adverse outcome pathways

The toxicological effects of many stressors are mediated through unknown, or incompletely characterized, mechanisms of action. The application of reverse engineering complex interaction networks from high dimensional omics data (gene, protein, metabolic, signaling) can be used to overcome these limitations. This approach was used to characterize adverse outcome pathways (AOPs) for chemicals that disrupt the hypothalamus-pituitary-gonadal endocrine axis in fathead minnows (FHM, Pimephales promelas). Gene expression changes in FHM ovaries in response to seven different chemicals, over different times, doses, and in vivo versus in vitro conditions, were captured in a large data set of 868 arrays. Potential AOPs of the antiandrogen flutamide were examined using two mutual information-based methods to infer gene regulatory networks and potential AOPs. Representative networks from these studies were used to predict network paths from stressor to adverse outcome as candidate AOPs. The relationship of individual chemicals to an adverse outcome can be determined by following perturbations through the network in response to chemical treatment, thus leading to the nodes associated with the adverse outcome. Identification of candidate pathways allows for formation of testable hypotheses about key biological processes, biomarkers, or alternative endpoints that can be used to monitor an AOP. Finally, the unique challenges facing the application of this approach in ecotoxicology were identified and a road map for the utilization of these tools presented.

In situ to in silico and back: elucidating the physiology and ecology of Geobacter spp. using genome-scale modelling

There is a wide diversity of unexplored metabolism encoded in the genomes of microorganisms that have an important environmental role. Genome-scale metabolic modelling enables the individual reactions that are encoded in annotated genomes to be organized into a coherent whole, which can then be used to predict metabolic fluxes that will optimize cell function under a range of conditions. In this Review, we summarize a series of studies in which genome-scale metabolic modelling of Geobacter spp. has resulted in an in-depth understanding of their central metabolism and ecology. A similar iterative modelling and experimental approach could accelerate elucidation of the physiology and ecology of other microorganisms inhabiting a diversity of environments, and could guide optimization of the practical applications of these species.

Flux coupling analysis of metabolic networks is sensitive to missing reactions

Genome-scale metabolic reconstructions are comprehensive, yet incomplete, models of real-world metabolic networks. While flux coupling analysis (FCA) has proved an appropriate method for analyzing metabolic relationships and for detecting functionally related reactions in such models, little is known about the impact of missing reactions on the accuracy of FCA. Based on an alternative characterization of flux coupling relations using elementary flux modes, this paper studies the changes that flux coupling relations may undergo due to missing reactions. In particular, we show that two uncoupled reactions in a metabolic network may be detected as directionally, partially or fully coupled in an incomplete version of the same network. Even a single missing reaction can cause significant changes in flux coupling relations. In case of two consecutive Escherichia coli genome-scale networks many fully coupled reaction pairs in the incomplete network become directionally coupled or even uncoupled in the more complete reconstruction. In this context, we found gene expression correlation values being significantly higher for the pairs that remained fully coupled than for the uncoupled or directionally coupled pairs. Our study clearly suggests that FCA results are indeed sensitive to missing reactions. Since the currently available genome-scale metabolic models are incomplete, we advise to use FCA results with care.

Integration of metabolic databases for the reconstruction of genome-scale metabolic networks

BACKGROUND

Genome-scale metabolic reconstructions have been recognised as a valuable tool for a variety of applications ranging from metabolic engineering to evolutionary studies. However, the reconstruction of such networks remains an arduous process requiring a high level of human intervention. This process is further complicated by occurrences of missing or conflicting information and the absence of common annotation standards between different data sources.

RESULTS

In this article, we report a semi-automated methodology aimed at streamlining the process of metabolic network reconstruction by enabling the integration of different genome-wide databases of metabolic reactions. We present results obtained by applying this methodology to the metabolic network of the plant Arabidopsis thaliana. A systematic comparison of compounds and reactions between two genome-wide databases allowed us to obtain a high-quality core consensus reconstruction, which was validated for stoichiometric consistency. A lower level of consensus led to a larger reconstruction, which has a lower quality standard but provides a baseline for further manual curation.

CONCLUSION

This semi-automated methodology may be applied to other organisms and help to streamline the process of genome-scale network reconstruction in order to accelerate the transfer of such models to applications.

Data integration and analysis of biological networks

During the past decade, bottom-up and top-down approaches of network reconstruction have greatly facilitated integration and analysis of biological networks, including transcriptional, protein interaction, and metabolic networks. As increasing amounts of multidimensional high-throughput data become available, biological networks have also been upgraded, allowing more accurate understanding of whole cellular characteristics. The network size is constantly expanding as larger volume of information and omics data are further integrated into the biological networks previously built upon a single type of data. Such effort more recently led to the modeling of human metabolic network and prediction of its tissue-specific metabolism, reconstruction of consensus yeast metabolic network, and simulation of mutual interactions among multiple microorganisms. It is expected that this trend will continue, the outcomes of which will allow development of more sophisticated networks integrating diverse omics data, and enhance our understanding of biological systems.

MR connectomics: Principles and challenges

MR connectomics is an emerging framework in neuro-science that combines diffusion MRI and whole brain tractography methodologies with the analytical tools of network science. In the present work we review the current methods enabling structural connectivity mapping with MRI and show how such data can be used to infer new information of both brain structure and function. We also list the technical challenges that should be addressed in the future to achieve high-resolution maps of structural connectivity. From the resulting tremendous amount of data that is going to be accumulated soon, we discuss what new challenges must be tackled in terms of methods for advanced network analysis and visualization, as well data organization and distribution. This new framework is well suited to investigate key questions on brain complexity and we try to foresee what fields will most benefit from these approaches.

Dissection of a complex transcriptional response using genome-wide transcriptional modelling

Modern genomics technologies generate huge data sets creating a demand for systems level, experimentally verified, analysis techniques. We examined the transcriptional response to DNA damage in a human T cell line (MOLT4) using microarrays. By measuring both mRNA accumulation and degradation over a short time course, we were able to construct a mechanistic model of the transcriptional response. The model predicted three dominant transcriptional activity profiles—an early response controlled by NFκB and c-Jun, a delayed response controlled by p53, and a late response related to cell cycle re-entry. The method also identified, with defined confidence limits, the transcriptional targets associated with each activity. Experimental inhibition of NFκB, c-Jun and p53 confirmed that target predictions were accurate. Model predictions directly explained 70% of the 200 most significantly upregulated genes in the DNA-damage response. Genome-wide transcriptional modelling (GWTM) requires no prior knowledge of either transcription factors or their targets. GWTM is an economical and effective method for identifying the main transcriptional activators in a complex response and confidently predicting their targets.

Evidence for the adaptation of protein pH-dependence to subcellular pH

BACKGROUND

The availability of genome sequences, and inferred protein coding genes, has led to several proteome-wide studies of isoelectric points. Generally, isoelectric points are distributed following variations on a biomodal theme that originates from the predominant acid and base amino acid sidechain pKas. The relative populations of the peaks in such distributions may correlate with environment, either for a whole organism or for subcellular compartments. There is also a tendency for isoelectric points averaged over a subcellular location to not coincide with the local pH, which could be related to solubility. We now calculate the correlation of other pH-dependent properties, calculated from 3D structure, with subcellular pH.

RESULTS

For proteins with known structure and subcellular annotation, the predicted pH at which a protein is most stable, averaged over a location, gives a significantly better correlation with subcellular pH than does isoelectric point. This observation relates to the cumulative properties of proteins, since maximal stability for individual proteins follows the bimodal isoelectric point distribution. Histidine residue location underlies the correlation, a conclusion that is tested against a background of proteins randomised with respect to this feature, and for which the observed correlation drops substantially.

CONCLUSION

There exists a constraint on protein pH-dependence, in relation to the local pH, that is manifested in the pKa distribution of histidine sub-proteomes. This is discussed in terms of protein stability, pH homeostasis, and fluctuations in proton concentration.

Bioinformatics in bioinorganic chemistry

Bioinformatics is a central discipline in modern life sciences aimed at describing the complex properties of living organisms starting from large-scale data sets of cellular constituents such as genes and proteins. In order for this wealth of information to provide useful biological knowledge, databases and software tools for data collection, analysis and interpretation need to be developed. In this paper, we review recent advances in the design and implementation of bioinformatics resources devoted to the study of metals in biological systems, a research field traditionally at the heart of bioinorganic chemistry. We show how metalloproteomes can be extracted from genome sequences, how structural properties can be related to function, how databases can be implemented, and how hints on interactions can be obtained from bioinformatics.

Systems biology of the cell cycle of Saccharomyces cerevisiae: From network mining to system-level properties

Following a brief description of the operational procedures of systems biology (SB), the cell cycle of budding yeast is discussed as a successful example of a top-down SB analysis. After the reconstruction of the steps that have led to the identification of a sizer plus timer network in the G1 to S transition, it is shown that basic functions of the cell cycle (the setting of the critical cell size and the accuracy of DNA replication) are system-level properties, detected only by integrating molecular analysis with modelling and simulation of their underlying networks. A detailed network structure of a second relevant regulatory step of the cell cycle, the exit from mitosis, derived from extensive data mining, is constructed and discussed. To reach a quantitative understanding of how nutrients control, through signalling, metabolism and transcription, cell growth and cycle is a very relevant aim of SB. Since we know that about 900 gene products are required for cell cycle execution and control in budding yeast, it is quite clear that a purely systematic approach would require too much time. Therefore lines for a modular SB approach, which prioritises molecular and computational investigations for faster cell cycle understanding, are proposed. The relevance of the insight coming from the cell cycle SB studies in developing a new framework for tackling very complex biological processes, such as cancer and aging, is discussed.