Filling Kinetic Gaps: Dynamic Modeling of Metabolism Where Detailed Kinetic Information Is Lacking

Integrative analysis between dynamical modeling of metabolic networks and data obtained from high throughput technology represents a worthy effort toward a holistic understanding of the link among phenotype and dynamical response. Even though the theoretical foundation for modeling metabolic network has been extensively treated elsewhere, the lack of kinetic information has limited the analysis in most of the cases. To overcome this constraint, we present and illustrate a new statistical approach that has two purposes: integrate high throughput data and survey the general dynamical mechanisms emerging for a slightly perturbed metabolic network.This paper presents a statistic framework capable to study how and how fast the metabolites participating in a perturbed metabolic network reach a steady-state. Instead of requiring accurate kinetic information, this approach uses high throughput metabolome technology to define a feasible kinetic library, which constitutes the base for identifying, statistical and dynamical properties during the relaxation. For the sake of illustration we have applied this approach to the human Red blood cell metabolism (hRBC) and its capacity to predict temporal phenomena was evaluated. Remarkable, the main dynamical properties obtained from a detailed kinetic model in hRBC were recovered by our statistical approach. Furthermore, robust properties in time scale and metabolite organization were identify and one concluded that they are a consequence of the combined performance of redundancies and variability in metabolite participation.In this work we present an approach that integrates high throughput metabolome data to define the dynamic behavior of a slightly perturbed metabolic network where kinetic information is lacking. Having information of metabolite concentrations at steady-state, this method has significant relevance due its potential scope to analyze others genome scale metabolic reconstructions. Thus, I expect this approach will significantly contribute to explore the relationship between dynamic and physiology in other metabolic reconstructions, particularly those whose kinetic information is practically nulls. For instances, I envisage that this approach can be useful in genomic medicine or pharmacogenomics, where the estimation of time scales and the identification of metabolite organization may be crucial to characterize and identify (dis)functional stages

Modeling Core Metabolism in Cancer Cells: Surveying the Topology Underlying the Warburg Effect

BACKGROUND: Alterations on glucose consumption and biosynthetic activity of amino acids, lipids and nucleotides are metabolic changes for sustaining cell proliferation in cancer cells. Irrevocable evidence of this fact is the Warburg effect which establishes that cancer cells prefers glycolysis over oxidative phosphorylation to generate ATP. Regulatory action over metabolic enzymes has opened a new window for designing more effective anti-cancer treatments. This enterprise is not trivial and the development of computational models that contribute to identifying potential enzymes for breaking the robustness of cancer cells is a priority. METHODOLOGY/PRINCIPAL FINDINGS: This work presents a constraint-base modeling of the most experimentally studied metabolic pathways supporting cancer cells: glycolysis, TCA cycle, pentose phosphate, glutaminolysis and oxidative phosphorylation. To evaluate its predictive capacities, a growth kinetics study for Hela cell lines was accomplished and qualitatively compared with in silico predictions. Furthermore, based on pure computational criteria, we concluded that a set of enzymes (such as lactate dehydrogenase and pyruvate dehydrogenase) perform a pivotal role in cancer cell growth, findings supported by an experimental counterpart. CONCLUSIONS/SIGNIFICANCE: Alterations on metabolic activity are crucial to initiate and sustain cancer phenotype. In this work, we analyzed the phenotype capacities emerged from a constructed metabolic network conformed by the most experimentally studied pathways sustaining cancer cell growth. Remarkably, in silico model was able to resemble the physiological conditions in cancer cells and successfully identified some enzymes currently studied by its therapeutic effect. Overall, we supplied evidence that constraint-based modeling constitutes a promising computational platform to: 1) integrate high throughput technology and establish a crosstalk between experimental validation and in silico prediction in cancer cell phenotype; 2) explore the fundamental metabolic mechanism that confers robustness in cancer; and 3) suggest new metabolic targets for anticancer treatments. All these issues being central to explore cancer cell metabolism from a systems biology perspective

A Systems Biology Approach to Drug Targets in Pseudomonas aeruginosa Biofilm

Antibiotic resistance is an increasing problem in the health care system and we are in a constant race with evolving bacteria. Biofilm-associated growth is thought to play a key role in bacterial adaptability and antibiotic resistance. We employed a systems biology approach to identify candidate drug targets for biofilm-associated bacteria by imitating specific microenvironments found in microbial communities associated with biofilm formation. A previously reconstructed metabolic model of Pseudomonas aeruginosa (PA) was used to study the effect of gene deletion on bacterial growth in planktonic and biofilm-like environmental conditions. A set of 26 genes essential in both conditions was identified. Moreover, these genes have no homology with any human gene. While none of these genes were essential in only one of the conditions, we found condition-dependent genes, which could be used to slow growth specifically in biofilm-associated PA. Furthermore, we performed a double gene deletion study and obtained 17 combinations consisting of 21 different genes, which were conditionally essential. While most of the difference in double essential gene sets could be explained by different medium composition found in biofilm-like and planktonic conditions, we observed a clear effect of changes in oxygen availability on the growth performance. Eight gene pairs were found to be synthetic lethal in oxygen-limited conditions. These gene sets may serve as novel metabolic drug targets to combat particularly biofilm-associated PA. Taken together, this study demonstrates that metabolic modeling of human pathogens can be used to identify oxygen-sensitive drug targets and thus, that this systems biology approach represents a powerful tool to identify novel candidate antibiotic targets

Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets

<p>Abstract</p> <p>Background:</p> <p><it>Mycobacterium tuberculosis </it>continues to be a major pathogen in the third world, killing almost 2 million people a year by the most recent estimates. Even in industrialized countries, the emergence of multi-drug resistant (MDR) strains of tuberculosis hails the need to develop additional medications for treatment. Many of the drugs used for treatment of tuberculosis target metabolic enzymes. Genome-scale models can be used for analysis, discovery, and as hypothesis generating tools, which will hopefully assist the rational drug development process. These models need to be able to assimilate data from large datasets and analyze them.</p> <p>Results:</p> <p>We completed a bottom up reconstruction of the metabolic network of <it>Mycobacterium tuberculosis </it>H37Rv. This functional <it>in silico </it>bacterium, <it>iNJ</it>661, contains 661 genes and 939 reactions and can produce many of the complex compounds characteristic to tuberculosis, such as mycolic acids and mycocerosates. We grew this bacterium <it>in silico </it>on various media, analyzed the model in the context of multiple high-throughput data sets, and finally we analyzed the network in an 'unbiased' manner by calculating the Hard Coupled Reaction (HCR) sets, groups of reactions that are forced to operate in unison due to mass conservation and connectivity constraints.</p> <p>Conclusion:</p> <p>Although we observed growth rates comparable to experimental observations (doubling times ranging from about 12 to 24 hours) in different media, comparisons of gene essentiality with experimental data were less encouraging (generally about 55%). The reasons for the often conflicting results were multi-fold, including gene expression variability under different conditions and lack of complete biological knowledge. Some of the inconsistencies between <it>in vitro </it>and <it>in silico </it>or <it>in vivo </it>and <it>in silico </it>results highlight specific loci that are worth further experimental investigations. Finally, by considering the HCR sets in the context of known drug targets for tuberculosis treatment we proposed new alternative, but equivalent drug targets.</p

On the Accessibility of Adaptive Phenotypes of a Bacterial Metabolic Network

The mechanisms by which adaptive phenotypes spread within an evolving population after their emergence are understood fairly well. Much less is known about the factors that influence the evolutionary accessibility of such phenotypes, a pre-requisite for their emergence in a population. Here, we investigate the influence of environmental quality on the accessibility of adaptive phenotypes of Escherichia coli's central metabolic network. We used an established flux-balance model of metabolism as the basis for a genotype-phenotype map (GPM). We quantified the effects of seven qualitatively different environments (corresponding to both carbohydrate and gluconeogenic metabolic substrates) on the structure of this GPM. We found that the GPM has a more rugged structure in qualitatively poorer environments, suggesting that adaptive phenotypes could be intrinsically less accessible in such environments. Nevertheless, on average ∼74% of the genotype can be altered by neutral drift, in the environment where the GPM is most rugged; this could allow evolving populations to circumvent such ruggedness. Furthermore, we found that the normalized mutual information (NMI) of genotype differences relative to phenotype differences, which measures the GPM's capacity to transmit information about phenotype differences, is positively correlated with (simulation-based) estimates of the accessibility of adaptive phenotypes in different environments. These results are consistent with the predictions of a simple analytic theory that makes explicit the relationship between the NMI and the speed of adaptation. The results suggest an intuitive information-theoretic principle for evolutionary adaptation; adaptation could be faster in environments where the GPM has a greater capacity to transmit information about phenotype differences. More generally, our results provide insight into fundamental environment-specific differences in the accessibility of adaptive phenotypes, and they suggest opportunities for research at the interface between information theory and evolutionary biology

Determining the association of medical co-morbidity with subjective and objective cognitive performance in an inner city memory disorders clinic: a retrospective chart review

<p>Abstract</p> <p>Background</p> <p>Medical co-morbidity may be associated with impaired cognitive function based on prior studies. However, no studies to date have determined to what extent this association is linked to medical illness or other factors that may be linked to medical illness (such as education, income levels, depression or subjective memory loss). The present study examined how medical co-morbidity, socioeconomic status (defined as residential SES), education and depression are associated with subjective and objective memory function in a sample of patients recruited from a university affiliated Memory Disorders Clinic located in a large Canadian inner city teaching hospital.</p> <p>Methods</p> <p>Data was collected from 85 consecutive referrals to an Inner City Memory Disorders Clinic including socio-demographic characteristics, cognitive status and medical co-morbidity. Descriptive and correlational analyses were conducted.</p> <p>Results</p> <p>Impaired objective cognitive function correlated significantly with increased medical co-morbidity and partially with education but not with residential SES or depression. Elevated memory complaints correlated significantly with depression, inversely with residential SES and not at all with medical co-morbidity or education.</p> <p>Conclusions</p> <p>Increased medical co-morbidity is significantly associated with impaired cognitive performance but not with subjective memory complaints in an Inner City Memory Clinic sample.</p

Inferring Carbon Sources from Gene Expression Profiles Using Metabolic Flux Models

Background: Bacteria have evolved the ability to efficiently and resourcefully adapt to changing environments. A key means by which they optimize their use of available nutrients is through adjustments in gene expression with consequent changes in enzyme activity. We report a new method for drawing environmental inferences from gene expression data. Our method prioritizes a list of candidate carbon sources for their compatibility with a gene expression profile using the framework of flux balance analysis to model the organism’s metabolic network. Principal Findings: For each of six gene expression profiles for Escherichia coli grown under differing nutrient conditions, we applied our method to prioritize a set of eighteen different candidate carbon sources. Our method ranked the correct carbon source as one of the top three candidates for five of the six expression sets when used with a genome-scale model. The correct candidate ranked fifth in the remaining case. Additional analyses show that these rankings are robust with respect to biological and measurement variation, and depend on specific gene expression, rather than general expression level. The gene expression profiles are highly adaptive: simulated production of biomass averaged 94.84% of maximum when the in silico carbon source matched the in vitro source of the expression profile, and 65.97% when it did not. Conclusions: Inferences about a microorganism’s nutrient environment can be made by integrating gene expression data into a metabolic framework. This work demonstrates that reaction flux limits for a model can be computed which are realistic in the sense that they affect in silico growth in a manner analogous to that in which a microorganism’s alteration of gene expression is adaptive to its nutrient environment.National Institute of Allergy and Infectious Diseases (U.S.) (grant HHSN 2722008000059C)National Institute of Allergy and Infectious Diseases (U.S.) (grant HHSN 26620040000IC)Bill & Melinda Gates Foundation (grant 18651010-37352-A

Genome-Scale Reconstruction of Escherichia coli's Transcriptional and Translational Machinery: A Knowledge Base, Its Mathematical Formulation, and Its Functional Characterization

Metabolic network reconstructions represent valuable scaffolds for ‘-omics’ data integration and are used to computationally interrogate network properties. However, they do not explicitly account for the synthesis of macromolecules (i.e., proteins and RNA). Here, we present the first genome-scale, fine-grained reconstruction of Escherichia coli's transcriptional and translational machinery, which produces 423 functional gene products in a sequence-specific manner and accounts for all necessary chemical transformations. Legacy data from over 500 publications and three databases were reviewed, and many pathways were considered, including stable RNA maturation and modification, protein complex formation, and iron–sulfur cluster biogenesis. This reconstruction represents the most comprehensive knowledge base for these important cellular functions in E. coli and is unique in its scope. Furthermore, it was converted into a mathematical model and used to: (1) quantitatively integrate gene expression data as reaction constraints and (2) compute functional network states, which were compared to reported experimental data. For example, the model predicted accurately the ribosome production, without any parameterization. Also, in silico rRNA operon deletion suggested that a high RNA polymerase density on the remaining rRNA operons is needed to reproduce the reported experimental ribosome numbers. Moreover, functional protein modules were determined, and many were found to contain gene products from multiple subsystems, highlighting the functional interaction of these proteins. This genome-scale reconstruction of E. coli's transcriptional and translational machinery presents a milestone in systems biology because it will enable quantitative integration of ‘-omics’ datasets and thus the study of the mechanistic principles underlying the genotype–phenotype relationship

Context-Specific Metabolic Networks Are Consistent with Experiments

Reconstructions of cellular metabolism are publicly available for a variety of different microorganisms and some mammalian genomes. To date, these reconstructions are “genome-scale” and strive to include all reactions implied by the genome annotation, as well as those with direct experimental evidence. Clearly, many of the reactions in a genome-scale reconstruction will not be active under particular conditions or in a particular cell type. Methods to tailor these comprehensive genome-scale reconstructions into context-specific networks will aid predictive in silico modeling for a particular situation. We present a method called Gene Inactivity Moderated by Metabolism and Expression (GIMME) to achieve this goal. The GIMME algorithm uses quantitative gene expression data and one or more presupposed metabolic objectives to produce the context-specific reconstruction that is most consistent with the available data. Furthermore, the algorithm provides a quantitative inconsistency score indicating how consistent a set of gene expression data is with a particular metabolic objective. We show that this algorithm produces results consistent with biological experiments and intuition for adaptive evolution of bacteria, rational design of metabolic engineering strains, and human skeletal muscle cells. This work represents progress towards producing constraint-based models of metabolism that are specific to the conditions where the expression profiling data is available