SBMLBioModels: M - M

Several liver disorders result from perturbations in the metabolism of hepatocytes, and their underlying mechanisms can be outlined through the use of genome-scale metabolic models (GEMs). Here we reconstruct a consensus GEM for hepatocytes, which we call iHepatocytes2322, that extends previous models by including an extensive description of lipid metabolism. We build iHepatocytes2322 using Human Metabolic Reaction 2.0 database and proteomics data in Human Protein Atlas, which experimentally validates the incorporated reactions. The reconstruction process enables improved annotation of the proteomics data using the network centric view of iHepatocytes2322. We then use iHepatocytes2322 to analyse transcriptomics data obtained from patients with non-alcoholic fatty liver disease. We show that blood concentrations of chondroitin and heparan sulphates are suitable for diagnosing non-alcoholic steatohepatitis and for the staging of non-alcoholic fatty liver disease. Furthermore, we observe serine deficiency in patients with NASH and identify PSPH, SHMT1 and BCAT1 as potential therapeutic targets for the treatment of non-alcoholic steatohepatitis. link: http://identifiers.org/pubmed/24419221

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice. link: http://identifiers.org/pubmed/26475342

A wave of apoptosis (programmed cell death) occurs normally in pancreatic beta-cells of newborn mice. We previously showed that macrophages from non-obese diabetic (NOD) mice become activated more slowly and engulf apoptotic cells at a lower rate than macrophages from control (Balb/c) mice. It has been hypothesized that this low clearance could result in secondary necrosis, escalating inflammation and self-antigen presentation that later triggers autoimmune, Type 1 diabetes (T1D). We here investigate whether this hypothesis could offer a reasonable and parsimonious explanation for onset of T1D in NOD mice. We quantify variants of the Copenhagen model (Freiesleben De Blasio et al. 1999 Diabetes 48, 1677), based on parameters from NOD and Balb/c experimental data. We show that the original Copenhagen model fails to explain observed phenomena within a reasonable range of parameter values, predicting an unrealistic all-or-none disease occurrence for both strains. However, if we take into account that, in general, activated macrophages produce harmful cytokines only when engulfing necrotic (but not apoptotic) cells, then the revised model becomes qualitatively and quantitatively reasonable. Further, we show that known differences between NOD and Balb/c mouse macrophage kinetics are large enough to account for the fact that an apoptotic wave can trigger escalating inflammatory response in NOD, but not Balb/c mice. In Balb/c mice, macrophages clear the apoptotic wave so efficiently, that chronic inflammation is prevented. link: http://identifiers.org/pubmed/16608707

Intracellular calcium oscillations, which are oscillatory changes of cytosolic calcium concentration in response to agonist stimulation, are experimentally well observed in various living cells. Simple calcium oscillations represent the most common pattern and many mathematical models have been published to describe this type of oscillation. On the other hand, relatively few theoretical studies have been proposed to give an explanation of complex intracellular calcium oscillations, such as bursting and chaos. In this paper, we develop a new possible mechanism for complex calcium oscillations based on the interplay between three calcium stores in the cell: the endoplasmic reticulum (ER), mitochondria and cytosolic proteins. The majority ( approximately 80%) of calcium released from the ER is first very quickly sequestered by mitochondria. Afterwards, a much slower release of calcium from the mitochondria serves as the calcium supply for the intermediate calcium exchanges between the ER and the cytosolic proteins causing bursting calcium oscillations. Depending on the permeability of the ER channels and on the kinetic properties of calcium binding to the cytosolic proteins, different patterns of complex calcium oscillations appear. With our model, we are able to explain simple calcium oscillations, bursting and chaos. Chaos is also observed for calcium oscillations in the bursting mode. link: http://identifiers.org/pubmed/11004387

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity. link: http://identifiers.org/pubmed/14744999

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity. link: http://identifiers.org/pubmed/14744999

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity. link: http://identifiers.org/pubmed/14744999

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity. link: http://identifiers.org/pubmed/14744999

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity. link: http://identifiers.org/pubmed/14744999

Mitogen-activated protein kinase (MAPK) cascades can operate as bistable switches residing in either of two different stable states. MAPK cascades are often embedded in positive feedback loops, which are considered to be a prerequisite for bistable behavior. Here we demonstrate that in the absence of any imposed feedback regulation, bistability and hysteresis can arise solely from a distributive kinetic mechanism of the two-site MAPK phosphorylation and dephosphorylation. Importantly, the reported kinetic properties of the kinase (MEK) and phosphatase (MKP3) of extracellular signal-regulated kinase (ERK) fulfill the essential requirements for generating a bistable switch at a single MAPK cascade level. Likewise, a cycle where multisite phosphorylations are performed by different kinases, but dephosphorylation reactions are catalyzed by the same phosphatase, can also exhibit bistability and hysteresis. Hence, bistability induced by multisite covalent modification may be a widespread mechanism of the control of protein activity. link: http://identifiers.org/pubmed/14744999

The ability of reactive oxygen species (ROS) to cause molecular damage has meant that chronic oxidative stress has been mostly studied from the point of view of being a source of toxicity to the cell. However, the known duality of ROS molecules as both damaging agents and cellular redox signals implies another perspective in the study of sustained oxidative stress. This is a perspective of studying oxidative stress as a constitutive signal within the cell. In this work, we adopt a theoretical perspective as an exploratory and explanatory approach to examine how chronic oxidative stress can interfere with signal processing by redox signalling pathways in the cell. We report that constitutive signals can give rise to a 'molecular habituation' effect that can prime for a gradual loss of biological function. This is because a constitutive signal in the environment has the potential to reduce the responsiveness of a signalling pathway through the prolonged activation of negative regulators. Additionally, we demonstrate how this phenomenon is likely to occur in different signalling pathways exposed to persistent signals and furthermore at different levels of biological organisation. link: http://identifiers.org/pubmed/29146308

The ability of reactive oxygen species (ROS) to cause molecular damage has meant that chronic oxidative stress has been mostly studied from the point of view of being a source of toxicity to the cell. However, the known duality of ROS molecules as both damaging agents and cellular redox signals implies another perspective in the study of sustained oxidative stress. This is a perspective of studying oxidative stress as a constitutive signal within the cell. In this work, we adopt a theoretical perspective as an exploratory and explanatory approach to examine how chronic oxidative stress can interfere with signal processing by redox signalling pathways in the cell. We report that constitutive signals can give rise to a 'molecular habituation' effect that can prime for a gradual loss of biological function. This is because a constitutive signal in the environment has the potential to reduce the responsiveness of a signalling pathway through the prolonged activation of negative regulators. Additionally, we demonstrate how this phenomenon is likely to occur in different signalling pathways exposed to persistent signals and furthermore at different levels of biological organisation. link: http://identifiers.org/pubmed/29146308

The ability of reactive oxygen species (ROS) to cause molecular damage has meant that chronic oxidative stress has been mostly studied from the point of view of being a source of toxicity to the cell. However, the known duality of ROS molecules as both damaging agents and cellular redox signals implies another perspective in the study of sustained oxidative stress. This is a perspective of studying oxidative stress as a constitutive signal within the cell. In this work, we adopt a theoretical perspective as an exploratory and explanatory approach to examine how chronic oxidative stress can interfere with signal processing by redox signalling pathways in the cell. We report that constitutive signals can give rise to a 'molecular habituation' effect that can prime for a gradual loss of biological function. This is because a constitutive signal in the environment has the potential to reduce the responsiveness of a signalling pathway through the prolonged activation of negative regulators. Additionally, we demonstrate how this phenomenon is likely to occur in different signalling pathways exposed to persistent signals and furthermore at different levels of biological organisation. link: http://identifiers.org/pubmed/29146308

The ability of reactive oxygen species (ROS) to cause molecular damage has meant that chronic oxidative stress has been mostly studied from the point of view of being a source of toxicity to the cell. However, the known duality of ROS molecules as both damaging agents and cellular redox signals implies another perspective in the study of sustained oxidative stress. This is a perspective of studying oxidative stress as a constitutive signal within the cell. In this work, we adopt a theoretical perspective as an exploratory and explanatory approach to examine how chronic oxidative stress can interfere with signal processing by redox signalling pathways in the cell. We report that constitutive signals can give rise to a 'molecular habituation' effect that can prime for a gradual loss of biological function. This is because a constitutive signal in the environment has the potential to reduce the responsiveness of a signalling pathway through the prolonged activation of negative regulators. Additionally, we demonstrate how this phenomenon is likely to occur in different signalling pathways exposed to persistent signals and furthermore at different levels of biological organisation. link: http://identifiers.org/pubmed/29146308

The ability of reactive oxygen species (ROS) to cause molecular damage has meant that chronic oxidative stress has been mostly studied from the point of view of being a source of toxicity to the cell. However, the known duality of ROS molecules as both damaging agents and cellular redox signals implies another perspective in the study of sustained oxidative stress. This is a perspective of studying oxidative stress as a constitutive signal within the cell. In this work, we adopt a theoretical perspective as an exploratory and explanatory approach to examine how chronic oxidative stress can interfere with signal processing by redox signalling pathways in the cell. We report that constitutive signals can give rise to a 'molecular habituation' effect that can prime for a gradual loss of biological function. This is because a constitutive signal in the environment has the potential to reduce the responsiveness of a signalling pathway through the prolonged activation of negative regulators. Additionally, we demonstrate how this phenomenon is likely to occur in different signalling pathways exposed to persistent signals and furthermore at different levels of biological organisation. link: http://identifiers.org/pubmed/29146308

The ability of reactive oxygen species (ROS) to cause molecular damage has meant that chronic oxidative stress has been mostly studied from the point of view of being a source of toxicity to the cell. However, the known duality of ROS molecules as both damaging agents and cellular redox signals implies another perspective in the study of sustained oxidative stress. This is a perspective of studying oxidative stress as a constitutive signal within the cell. In this work, we adopt a theoretical perspective as an exploratory and explanatory approach to examine how chronic oxidative stress can interfere with signal processing by redox signalling pathways in the cell. We report that constitutive signals can give rise to a 'molecular habituation' effect that can prime for a gradual loss of biological function. This is because a constitutive signal in the environment has the potential to reduce the responsiveness of a signalling pathway through the prolonged activation of negative regulators. Additionally, we demonstrate how this phenomenon is likely to occur in different signalling pathways exposed to persistent signals and furthermore at different levels of biological organisation. link: http://identifiers.org/pubmed/29146308

CD4+ T cells orchestrate the adaptive immune response in vertebrates. While both experimental and modeling work has been conducted to understand the molecular genetic mechanisms involved in CD4+ T cell responses and fate attainment, the dynamic role of intrinsic (produced by CD4+ T lymphocytes) versus extrinsic (produced by other cells) components remains unclear, and the mechanistic and dynamic understanding of the plastic responses of these cells remains incomplete. In this work, we studied a regulatory network for the core transcription factors involved in CD4+ T cell-fate attainment. We first show that this core is not sufficient to recover common CD4+ T phenotypes. We thus postulate a minimal Boolean regulatory network model derived from a larger and more comprehensive network that is based on experimental data. The minimal network integrates transcriptional regulation, signaling pathways and the micro-environment. This network model recovers reported configurations of most of the characterized cell types (Th0, Th1, Th2, Th17, Tfh, Th9, iTreg, and Foxp3-independent T regulatory cells). This transcriptional-signaling regulatory network is robust and recovers mutant configurations that have been reported experimentally. Additionally, this model recovers many of the plasticity patterns documented for different T CD4+ cell types, as summarized in a cell-fate map. We tested the effects of various micro-environments and transient perturbations on such transitions among CD4+ T cell types. Interestingly, most cell-fate transitions were induced by transient activations, with the opposite behavior associated with transient inhibitions. Finally, we used a novel methodology was used to establish that T-bet, TGF-β and suppressors of cytokine signaling proteins are keys to recovering observed CD4+ T cell plastic responses. In conclusion, the observed CD4+ T cell-types and transition patterns emerge from the feedback between the intrinsic or intracellular regulatory core and the micro-environment. We discuss the broader use of this approach for other plastic systems and possible therapeutic interventions. link: http://identifiers.org/pubmed/26090929

The kinetics of glyoxalase I [(R)-S-lactoylglutathione methylglyoxal-lyase; EC 4.4.1.5] and glyoxalase II (S-2-hydroxyacylglutathione hydrolase; EC 3.1.2.6) from Saccharomyces cerevisiae was studied in situ, in digitonin permeabilized cells, using two different approaches: initial rate analysis and progress curves analysis. Initial rate analysis was performed by hyperbolic regression of initial rates using the program HYPERFIT. Glyoxalase I exhibited saturation kinetics on 0.05-2.5 mM hemithioacetal concentration range, with kinetic parameters Km 0.53 +/- 0.07 mM and V (3.18 +/- 0.16) x 10(-2) mM.min(-1). Glyoxalase II also showed saturation kinetics in the SD-lactoylglutathione concentration range of 0.15-3 mM and Km 0.32 +/- 0.13 mM and V (1.03 +/- 0.10) x 10(-3) mM.min(-1) were obtained. The kinetic parameters of both enzymes were also estimated by nonlinear regression of progress curves using the raw absorbance data and integrated differential rate equations with the program GEPASI. Several optimization methods were used to minimize the sum of squares of residuals. The best parameter fit for the glyoxalase I reaction was obtained with a single curve analysis, using the irreversible Michaelis-Menten model. The kinetic parameters obtained, Km 0.62 +/- 0.18 mM and V (2.86 +/- 0.01) x 10(-2) mM.min(-1), were in agreement with those obtained by initial rate analysis. The results obtained for glyoxalase II, using either the irreversible Michaelis-Menten model or a phenomenological reversible hyperbolic model, showed a high correlation of residuals with time and/or high values of standard deviation associated with Km. The possible causes for the discrepancy between data obtained from initial rate analysis and progress curve analysis, for glyoxalase II, are discussed. link: http://identifiers.org/pubmed/11453985

A kinetic model for N-(1-deoxy-D-fructos-1-yl)-glycine (DFG) thermal decomposition was proposed. Two temperatures (100 and 120 degrees C) and two pHs (5.5 and 6.8) were studied. The measured responses were DFG, 3-deoxyosone, 1-deoxyosone, methylglyoxal, acetic acid, formic acid, glucose, fructose, mannose and melanoidins. For each system the model parameters, the rate constants, were estimated by non-linear regression, via multiresponse modelling. The determinant criterion was used as the statistical fit criterion. Model discrimination was performed by both chemical insight and statistical tests (Posterior Probability and Akaike criterion). Kinetic analysis showed that at lower pH DFG 1,2-enolization is favoured whereas with increasing pH 2,3-enolization becomes a more relevant degradation pathway. The lower amount observed of 1-DG is related with its high reactivity. It was shown that acetic acid, a main degradation product from DFG, was mainly formed through 1-DG degradation. Also from the estimated parameters 3-DG was found to be the main precursor in carbohydrate fragments formation, responsible for colour formation. Some indication was given that as the reaction proceeded other compounds besides DFG become reactants themselves with the formation among others of methylglyoxal. The multiresponse kinetic analysis was shown to be both helpful in deriving relevant kinetic parameters as well as in obtaining insight into the reaction mechanism. link: http://identifiers.org/pubmed/12873422

Saccharomyces cerevisiae cultures growing exponentially and after starvation are distinctly different, as shown by several studies at the physiological, biochemical, and morphological levels. One group of studies attempted to be mechanistic, characterizing a few molecules and interactions, while another focused on global observations but remained descriptive or at best phenomenological. Recent advances in large-scale molecular profiling technologies, theoretical, and computational biology, are making possible integrative studies of biological systems, where global observations are explained through computational models with solid theoretical bases. A case study of the systems biology approach applied to the characterization of baker's yeast cultures in exponential growth and post-diauxic phases is presented.

Twenty cell cultures of S. cerevisiae were grown under similar environmental conditions. Samples from ten of these cultures were collected 11 hours after inoculation, while samples from the other ten were collected 4 days after inoculation. These samples were analyzed for their RNA and metabolite composition using Affymetrix chips and gas chromatography-mass spectrometry (GC-MS). The data were interpreted with statistical analyses and through the use of computer simulations of a kinetic model that was built by merging two independent models of glycolysis and glycerol biosynthesis. The simulation results agree with the exponential growth phase data, while no model is available for the post-diauxic phase. We discuss the need for expanding the number of kinetic models of S. cerevisiae, combining metabolism and genetic regulation, and covering all of its biochemistry. link: http://identifiers.org/doi/10.2174/1389202043348643

A possible therapeutic strategy for amyloid diseases involves the use of small molecule compounds to inhibit protein assembly into insoluble aggregates. According to the recently proposed Crystallization-Like Model, the kinetics of amyloid fibrillization can be retarded by decreasing the frequency of new fibril formation or by decreasing the elongation rate of existing fibrils. To the compounds that affect the nucleation and/or the growth steps we call true inhibitors. An apparent inhibition mechanism may however result from the alteration of thermodynamic properties such as the solubility of the amyloidogenic protein. Apparent inhibitors markedly influence protein aggregation kinetics measured in vitro, yet they are likely to lead to disappointing results when tested in vivo. This is because cells and tissues media are in general much more buffered against small variations in composition than the solutions prepared in lab. Here we show how to discriminate between true and apparent inhibition mechanisms from experimental data on protein aggregation kinetics. The goal is to be able to identify false positives much earlier during the drug development process. link: http://identifiers.org/pubmed/23232498

Astrocytes are the most abundant cells of the central nervous system; they have a predominant role in maintaining brain metabolism. In this sense, abnormal metabolic states have been found in different neuropathological diseases. Determination of metabolic states of astrocytes is difficult to model using current experimental approaches given the high number of reactions and metabolites present. Thus, genome-scale metabolic networks derived from transcriptomic data can be used as a framework to elucidate how astrocytes modulate human brain metabolic states during normal conditions and in neurodegenerative diseases. We performed a Genome-Scale Reconstruction of the Human Astrocyte Metabolic Network with the purpose of elucidating a significant portion of the metabolic map of the astrocyte. This is the first global high-quality, manually curated metabolic reconstruction network of a human astrocyte. It includes 5,007 metabolites and 5,659 reactions distributed among 8 cell compartments, (extracellular, cytoplasm, mitochondria, endoplasmic reticle, Golgi apparatus, lysosome, peroxisome and nucleus). Using the reconstructed network, the metabolic capabilities of human astrocytes were calculated and compared both in normal and ischemic conditions. We identified reactions activated in these two states, which can be useful for understanding the astrocytic pathways that are affected during brain disease. Additionally, we also showed that the obtained flux distributions in the model, are in accordance with literature-based findings. Up to date, this is the most complete representation of the human astrocyte in terms of inclusion of genes, proteins, reactions and metabolic pathways, being a useful guide for in-silico analysis of several metabolic behaviors of the astrocyte during normal and pathologic states. link: http://identifiers.org/pubmed/28243200

Mutants of Physarum polycephalum can be complemented by fusion of plasmodial cells followed by cytoplasmic mixing. Complementation between strains carrying different mutational defects in the sporulation control network may depend on the signaling state of the network components. We have previously suggested that time-resolved somatic complementation (TRSC) analysis with such mutants may be used to probe network architecture and dynamics. By computer simulation it is now shown how and under which conditions the regulatory hierarchy of genes can be determined experimentally. A kinetic model of the sporulation control network is developed, which is then used to demonstrate how the mechanisms of TRSC can be understood and simulated at the kinetic level. On the basis of theoretical considerations, experimental parameters that determine whether functional complementation of two mutations will occur are identified. It is also shown how gene dosage-effect relationships can be employed for network analysis. The theoretical framework provided may be used to systematically analyze network structure and dynamics through time-resolved somatic complementation studies. The conclusions drawn are of general relevance in that they do not depend on the validity of the model from which they were derived. link: http://identifiers.org/pubmed/12750324

Amyloid protein aggregates are implicated in many neurodegenerative diseases, including Alzheimer's disease and the prion diseases. Therapeutics to block amyloid formation are often tested in vitro, but it is not clear how to extrapolate from these experiments to a clinical setting, where the effective drug dose may be much lower. Here we address this question using a theoretical kinetic model to calculate the growth rate of protein aggregates as a function of the dose of each of three categories of drug. We find that therapeutics which block the growing ends of amyloids are the most promising, as alternative strategies may be ineffective or even accelerate amyloid formation at low drug concentrations. Our mathematical model can be used to identify and optimise an end-blocking drug in vitro. Our model also suggests an alternative explanation for data previously thought to prove the existence of an entity known as protein X. link: http://identifiers.org/pubmed/11152275

Individual ion channels or exchangers are described with a common set of equations for both the sinoatrial node pacemaker and ventricular cells. New experimental data are included, such as the new kinetics of the inward rectifier K+ channel, delayed rectifier K+ channel, and sustained inward current. The gating model of Shirokov et al. (J Gen Physiol 102: 1005-1030, 1993) is used for both the fast Na+ and L-type Ca2+ channels. When combined with a contraction model (Negroni and Lascano: J Mol Cell Cardiol 28: 915-929, 1996), the experimental staircase phenomenon of contraction is reconstructed. The modulation of the action potential by varying the external Ca2+ and K+ concentrations is well simulated. The conductance of I(CaL) dominates membrane conductance during the action potential so that an artificial increase of I(to), I(Kr), I(Ks), or I(KATP) magnifies I(CaL) amplitude. Repolarizing current is provided sequentially by I(Ks), I(Kr), and I(K1). Depression of ATP production results in the shortening of action potential through the activation of I(KATP). The ratio of Ca2+ released from SR over Ca2+ entering via I(CaL) (Ca2+ gain = approximately 15) in excitation-contraction coupling well agrees with the experimental data. The model serves as a predictive tool in generating testable hypotheses. link: http://identifiers.org/pubmed/12877767

Biochemical systems embed complex networks and hence development and analysis of their detailed models pose a challenge for computation. Coarse-grained biochemical models, called reduced-order models (ROMs), consisting of essential biochemical mechanisms are more useful for computational analysis and for studying important features of a biochemical network. The authors present a novel method to model-reduction by identifying potentially important parameters using multidimensional sensitivity analysis. A ROM is generated for the GTPase-cycle module of m1 muscarinic acetylcholine receptor, Gq, and regulator of G-protein signalling 4 (a GTPase-activating protein or GAP) starting from a detailed model of 48 reactions. The resulting ROM has only 17 reactions. The ROM suggested that complexes of G-protein coupled receptor (GPCR) and GAP–which were proposed in the detailed model as a hypothesis–are required to fit the experimental data. Models previously published in the literature are also simulated and compared with the ROM. Through this comparison, a minimal ROM, that also requires complexes of GPCR and GAP, with just 15 parameters is generated. The proposed reduced-order modelling methodology is scalable to larger networks and provides a general framework for the reduction of models of biochemical systems. link: http://identifiers.org/pubmed/16986265

Signal transducer and actuator of transcription (STATs) are a family of transcription factors activated by various cytokines, growth factors and hormones. They are important mediators of immune responses and growth and differentiation of various cell types. The STAT signalling system represents a defined functional module with a pattern of signalling that is conserved from flies to mammals. In order to probe and gain insights into the signalling properties of the STAT module by computational means, we developed a simple non-linear ordinary differential equations model within the 'Virtual Cell' framework. Our results demonstrate that the STAT module can operate as a 'biphasic amplitude filter' with an ability to amplify input signals within a specific intermediate range. We show that dimerisation of phosphorylated STAT is crucial for signal amplification and the amplitude filtering function. We also demonstrate that maximal amplification at intermediate levels of STAT activation is a moderately robust property of STAT module. We propose that these observations can be extrapolated to the analogous SMAD signalling module. link: http://identifiers.org/pubmed/17091582

The phosphatidylinositol (PI) cycle is central to eukaryotic cell signaling. Its complexity, due to the number of reactions and lipid and inositol phosphate intermediates involved makes it difficult to analyze experimentally. Computational modelling approaches are seen as a way forward to elucidate complex biological regulatory mechanisms when this cannot be achieved solely through experimental approaches. Whilst mathematical modelling is well established in informing biological systems, many models are often informed by data sourced from different cell types (mosaic data), to inform model parameters. For instance, kinetic rate constants are often determined from purified enzyme data in vitro or use experimental concentrations obtained from multiple unrelated cell types. Thus they do not represent any specific cell type nor fully capture cell specific behaviours. In this work, we develop a model of the PI cycle informed by in-vivo omics data taken from a single cell type, namely platelets. Our model recapitulates the known experimental dynamics before and after stimulation with different agonists and demonstrates the importance of lipid- and protein-binding proteins in regulating second messenger outputs. Furthermore, we were able to make a number of predictions regarding the regulation of PI cycle enzymes and the importance of the number of receptors required for successful GPCR signaling. We then consider how pathway behavior differs, when fully informed by data for HeLa cells and show that model predictions remain relatively consistent. However, when informed by mosaic experimental data model predictions greatly vary. Our work illustrates the risks of using mosaic datasets from unrelated cell types which leads to over 75% of outputs not fitting with expected behaviors. link:

The microbial community present in the human mouth is engaged in a complex network of diverse metabolic activities. In addition to serving as energy and building-block sources, metabolites are key players in interspecies and host-pathogen interactions. Metabolites are also implicated in triggering the local inflammatory response, which can affect systemic conditions such as atherosclerosis, obesity, and diabetes. While the genome of several oral pathogens has been sequenced, quantitative understanding of the metabolic functions of any oral pathogen at the system level has not been explored yet. Here we pursue the computational construction and analysis of the genome-scale metabolic network of Porphyromonas gingivalis, a gram-negative anaerobe that is endemic in the human population and largely responsible for adult periodontitis. Integrating information from the genome, online databases, and literature screening, we built a stoichiometric model that encompasses 679 metabolic reactions. By using flux balance approaches and automated network visualization, we analyze the growth capacity under amino-acid-rich medium and provide evidence that amino acid preference and cytotoxic by-product secretion rates are suitably reproduced by the model. To provide further insight into the basic metabolic functions of P. gingivalis and suggest potential drug targets, we study systematically how the network responds to any reaction knockout. We focus specifically on the lipopolysaccharide biosynthesis pathway and identify eight putative targets, one of which has been recently verified experimentally. The current model, which is amenable to further experimental testing and refinements, could prove useful in evaluating the oral microbiome dynamics and in the development of novel biomedical applications. link: http://identifiers.org/pubmed/18931137

Given the complexity of developmental networks, it is often difficult to predict the effect of genetic perturbations, even within coding genes. Regulatory factors generally have pleiotropic effects, exhibit partially redundant roles, and regulate highly interconnected pathways with ample cross-talk. Here, we delineate a logical model encompassing 48 components and 82 regulatory interactions involved in mesoderm specification during Drosophila development, thereby providing a formal integration of all available genetic information from the literature. The four main tissues derived from mesoderm correspond to alternative stable states. We demonstrate that the model can predict known mutant phenotypes and use it to systematically predict the effects of over 300 new, often non-intuitive, loss- and gain-of-function mutations, and combinations thereof. We further validated several novel predictions experimentally, thereby demonstrating the robustness of model. Logical modelling can thus contribute to formally explain and predict regulatory outcomes underlying cell fate decisions. link: http://identifiers.org/pubmed/27599298

1. The electrical activity of Cardiac Purkinje fibres was reconstructed using a mathematical model of the membrane current. The individual components of ionic curent were described by equations which wee based as closely as possible on previous experiments using the voltage clamp technique. 2. Membrane action potentials and pace-maker activity were calculated and compared with time course of underlying changes in two functionally distinct outeard currents, iX1 and iK2. 3. The repolarization of the theoretical action potential is triggered by the onset of iX1, which becomes activated over the plateau range of potentials. iK2 also activates during the plateau but does not play a controlling role in the repolarization. Hwever, iK2 does govern the slow pace-maker depolarization through its subsequent deactivation at negative potentials. 4. The individual phases of the calculated action potential and their 'experimental' modifications were compared with published records. The upstroke is generated by a Hodgkin-Huxley type sodium conductance (gNa), and rises with a maximum rate of 478 V/sec, somewhat less than experimentally observed values ( up to 800 V/sec). The discrepancy is discussed in relation to experimental attempts at measuring gNa. 5. The ole of the transient outward chloride current (called igr) was studied in calculations of the rapid phase of repolarization and 'notch' configuration...

link: http://identifiers.org/pubmed/1185607

BACKGROUND: Global demographic changes have stimulated marked interest in the process of aging. There has been, and will continue to be, an unrelenting rise in the number of the oldest old ( >85 years of age). Together with an ageing population there comes an increase in the prevalence of age related disease. Of the diseases of ageing, cardiovascular disease (CVD) has by far the highest prevalence. It is regarded that a finely tuned lipid profile may help to prevent CVD as there is a long established relationship between alterations to lipid metabolism and CVD risk. In fact elevated plasma cholesterol, particularly Low Density Lipoprotein Cholesterol (LDL-C) has consistently stood out as a risk factor for having a cardiovascular event. Moreover it is widely acknowledged that LDL-C may rise with age in both sexes in a wide variety of groups. The aim of this work was to use a whole-body mathematical model to investigate why LDL-C rises with age, and to test the hypothesis that mechanistic changes to cholesterol absorption and LDL-C removal from the plasma are responsible for the rise. The whole-body mechanistic nature of the model differs from previous models of cholesterol metabolism which have either focused on intracellular cholesterol homeostasis or have concentrated on an isolated area of lipoprotein dynamics. The model integrates both current and previously published data relating to molecular biology, physiology, ageing and nutrition in an integrated fashion. RESULTS: The model was used to test the hypothesis that alterations to the rate of cholesterol absorption and changes to the rate of removal of LDL-C from the plasma are integral to understanding why LDL-C rises with age. The model demonstrates that increasing the rate of intestinal cholesterol absorption from 50% to 80% by age 65 years can result in an increase of LDL-C by as much as 34 mg/dL in a hypothetical male subject. The model also shows that decreasing the rate of hepatic clearance of LDL-C gradually to 50% by age 65 years can result in an increase of LDL-C by as much as 116 mg/dL. CONCLUSIONS: Our model clearly demonstrates that of the two putative mechanisms that have been implicated in the dysregulation of cholesterol metabolism with age, alterations to the removal rate of plasma LDL-C has the most significant impact on cholesterol metabolism and small changes to the number of hepatic LDL receptors can result in a significant rise in LDL-C. This first whole-body systems based model of cholesterol balance could potentially be used as a tool to further improve our understanding of whole-body cholesterol metabolism and its dysregulation with age. Furthermore, given further fine tuning the model may help to investigate potential dietary and lifestyle regimes that have the potential to mitigate the effects aging has on cholesterol metabolism. link: http://identifiers.org/pubmed/23046614

Cells must respond specifically to different environmental stimuli in order to survive. The signal transduction pathways involved in sensing these stimuli often share the same or homologous proteins. Despite potential cross-wiring, cells show specificity of response. We show, through modeling, that the physiological response of such pathways exposed to simultaneous and temporally ordered inputs can demonstrate system-level mechanisms by which pathways achieve specificity. We apply these results to the hyperosmolar and pheromone mitogen-activated protein (MAP) kinase pathways in the yeast Saccharomyces cerevisiae. These two pathways specifically sense osmolar and pheromone signals, despite sharing a MAPKKK, Ste11, and having homologous MAPKs (Fus3 and Hog1). We show that in a single cell, the pathways are bistable over a range of inputs, and the cell responds to only one stimulus even when exposed to both. Our results imply that these pathways achieve specificity by filtering out spurious cross-talk through mutual inhibition. The variability between cells allows for heterogeneity of the decisions. link: http://identifiers.org/pubmed/17259986

Heart disease remains the leading cause of death globally. Although reperfusion following myocardial ischemia can prevent death by restoring nutrient flow, ischemia/reperfusion injury can cause significant heart damage. The mechanisms that drive ischemia/reperfusion injury are not well understood; currently, few methods can predict the state of the cardiac muscle cell and its metabolic conditions during ischemia. Here, we explored the energetic sustainability of cardiomyocytes, using a model for cellular metabolism to predict the levels of ATP following hypoxia. We modeled glycolytic metabolism with a system of coupled ordinary differential equations describing the individual metabolic reactions within the cardiomyocyte over time. Reduced oxygen levels and ATP consumption rates were simulated to characterize metabolite responses to ischemia. By tracking biochemical species within the cell, our model enables prediction of the cell's condition up to the moment of reperfusion. The simulations revealed a distinct transition between energetically sustainable and unsustainable ATP concentrations for various energetic demands. Our model illustrates how even low oxygen concentrations allow the cell to perform essential functions. We found that the oxygen level required for a sustainable level of ATP increases roughly linearly with the ATP consumption rate. An extracellular O2 concentration of ∼0.007 mm could supply basic energy needs in non-beating cardiomyocytes, suggesting that increased collateral circulation may provide an important source of oxygen to sustain the cardiomyocyte during extended ischemia. Our model provides a time-dependent framework for studying various intervention strategies to change the outcome of reperfusion. link: http://identifiers.org/pubmed/28487363

A new mechanism is proposed for the apparent breakthrough of HIV that occurs approximately 6 months after the commencement of therapy with zidovudine (AZT). Using a simple mathematical model of the interacting population dynamics of HIV and its major host cell in the circulation (the CD4+ lymphocyte), predicted patterns of HIV plasma viraemia in the weeks following treatment with zidovudine are generated. These are in close agreement with observed patterns despite the fact that the model contains no mechanisms for the development of drug-resistant strains of virus. It is suggested that the patterns of viral abundance observed during the first 6 months after treatment may be the result of non-linearities in the interactions between HIV and CD4+ cells, and that it is only after the first post-treatment burst of viral production that drug resistance plays an important role. link: http://identifiers.org/pubmed/1677807

The electrical response of pancreatic beta-cells to step increases in glucose concentration is biphasic, consisting of a prolonged depolarization with action potentials (Phase 1) followed by membrane potential oscillations known as bursts. We have proposed that the Phase 1 response results from the combined depolarizing influences of potassium channel closure and an inward, nonselective cation current (ICRAN) that activates as intracellular calcium stores empty during exposure to basal glucose (Bertram et al., 1995). The stores refill during Phase 1, deactivating ICRAN and allowing steady-state bursting to commence. We support this hypothesis with additional simulations and experimental results indicating that Phase 1 duration is sensitive to the filling state of intracellular calcium stores. First, the duration of the Phase 1 transient increases with duration of prior exposure to basal (2.8 mM) glucose, reflecting the increased time required to fill calcium stores that have been emptying for longer periods. Second, Phase 1 duration is reduced when islets are exposed to elevated K+ to refill calcium stores in the presence of basal glucose. Third, when extracellular calcium is removed during the basal glucose exposure to reduce calcium influx into the stores, Phase 1 duration increases. Finally, no Phase 1 is observed following hyperpolarization of the beta-cell membrane with diazoxide in the continued presence of 11 mm glucose, a condition in which intracellular calcium stores remain full. Application of carbachol to empty calcium stores during basal glucose exposure did not increase Phase 1 duration as the model predicts. Despite this discrepancy, the good agreement between most of the experimental results and the model predictions provides evidence that a calcium release-activated current mediates the Phase 1 electrical response of the pancreatic beta-cell. link: http://identifiers.org/pubmed/9002424

The generation of off-flavors in soybean homogenates such as n-hexanal via the lipoxygenase (LOX) pathway can be a problem in the processed food industry. Previous studies have examined the effect of using soybean varieties missing one or more of the 3 LOX isozymes on n-hexanal generation. A dynamic mathematical model of the soybean LOX pathway using ordinary differential equations was constructed using parameters estimated from existing data with the aim of predicting how n-hexanal generation could be reduced. Time-course simulations of LOX-null beans were run and compared with experimental results. Model L(2), L(3), and L(12) beans were within the range relative to the wild type found experimentally, with L(13) and L(23) beans close to the experimental range. Model L(1) beans produced much more n-hexanal relative to the wild type than those in experiments. Sensitivity analysis indicates that reducing the estimated K(m) parameter for LOX isozyme 3 (L-3) would improve the fit between model predictions and experimental results found in the literature. The model also predicts that increasing L-3 or reducing L-2 levels within beans may reduce n-hexanal generation.This work describes the use of mathematics to attempt to quantify the enzyme-catalyzed conversions of compounds in soybean homogenates into undesirable flavors, primarily from the compound n-hexanal. The effect of different soybean genotypes and enzyme kinetic constants was also studied, leading to recommendations on which combinations might minimize off-flavor levels and what further work might be carried out to substantiate these conclusions. link: http://identifiers.org/pubmed/21535565

A B S T R A C T The trajectories of the dynamic system which regulates the competition between the populations of malignant cells and immune cells may tend to an asymptotically stable equilibrium in which the sizes of these populations do not vary, which is called tumor dormancy. Especially for lower steady-state sizes of the population of malignant cells, this equilibrium represents a desirable clinical condition since the tumor growth is blocked. In this context, it is of mandatory importance to analyze the robustness of this clinical favorable state of health in the face of perturbations. To this end, the paper presents an optimization technique to determine whether an assigned rectangular region, which surrounds an asymptotically stable equilibrium point of a quadratic systems, is included into the domain of attraction of the equilibrium itself. The biological relevance of the application of this technique to the analysis of tumor growth dynamics is shown on the basis of a recent quadraticmodel of the tumor–immune system competition dynamics. Indeed the application of the proposedmethodology allows to ensure that a given safety region, determined on the basis of clinical considerations, belongs to the domain of attraction of the tumor blocked equilibrium; therefore for the set of perturbed initial conditions which belong to such region, the convergence to the healthy steady state is guaranteed. The proposed methodology can also provide an optimal strategy for cancer treatment. link: http://identifiers.org/doi/10.1016/j.bspc.2008.02.001

We present the quantification and kinetic characterisation of the enzymes of the pentose phosphate pathway in Saccharomyces cerevisiae. The data are combined into a mathematical model that describes the dynamics of this system and allows for the predicting changes in metabolite concentrations and fluxes in response to perturbations. We use the model to study the response of yeast to a glucose pulse. We then combine the model with an existing glycolysis one to study the effect of oxidative stress on carbohydrate metabolism. The combination of these two models was made possible by the standardized enzyme kinetic experiments carried out in both studies. This work demonstrates the feasibility of constructing larger network models by merging smaller pathway models. link: http://identifiers.org/doi/10.7287/peerj.preprints.146v2

We present the quantification and kinetic characterisation of the enzymes of the pentose phosphate pathway in Saccharomyces cerevisiae. The data are combined into a mathematical model that describes the dynamics of this system and allows for the predicting changes in metabolite concentrations and fluxes in response to perturbations. We use the model to study the response of yeast to a glucose pulse. We then combine the model with an existing glycolysis one to study the effect of oxidative stress on carbohydrate metabolism. The combination of these two models was made possible by the standardized enzyme kinetic experiments carried out in both studies. This work demonstrates the feasibility of constructing larger network models by merging smaller pathway models. link: http://identifiers.org/doi/10.7287/peerj.preprints.146v2

link: http://identifiers.org/pubmed/1867714

Seasonal and pandemic influenza A virus (IAV) continues to be a public health threat. However, we lack a detailed and quantitative understanding of the immune response kinetics to IAV infection and which biological parameters most strongly influence infection outcomes. To address these issues, we use modeling approaches combined with experimental data to quantitatively investigate the innate and adaptive immune responses to primary IAV infection. Mathematical models were developed to describe the dynamic interactions between target (epithelial) cells, influenza virus, cytotoxic T lymphocytes (CTLs), and virus-specific IgG and IgM. IAV and immune kinetic parameters were estimated by fitting models to a large data set obtained from primary H3N2 IAV infection of 340 mice. Prior to a detectable virus-specific immune response (before day 5), the estimated half-life of infected epithelial cells is approximately 1.2 days, and the half-life of free infectious IAV is approximately 4 h. During the adaptive immune response (after day 5), the average half-life of infected epithelial cells is approximately 0.5 days, and the average half-life of free infectious virus is approximately 1.8 min. During the adaptive phase, model fitting confirms that CD8(+) CTLs are crucial for limiting infected cells, while virus-specific IgM regulates free IAV levels. This may imply that CD4 T cells and class-switched IgG antibodies are more relevant for generating IAV-specific memory and preventing future infection via a more rapid secondary immune response. Also, simulation studies were performed to understand the relative contributions of biological parameters to IAV clearance. This study provides a basis to better understand and predict influenza virus immunity. link: http://identifiers.org/pubmed/20410284

The B cell response to influenza infection of the respiratory tract contributes to viral clearance and establishes profound resistance to reinfection by related viruses. Numerous studies have measured virus-specific antibody-secreting cell (ASC) frequencies in different anatomical compartments after influenza infection and provided a general picture of the kinetics of ASC formation and dispersion. However, the dynamics of ASC populations are difficult to determine experimentally and have received little attention. Here, we applied mathematical modeling to investigate the dynamics of ASC growth, death, and migration over the 2-week period following primary influenza infection in mice. Experimental data for model fitting came from high frequency measurements of virus-specific IgM, IgG, and IgA ASCs in the mediastinal lymph node (MLN), spleen, and lung. Model construction was based on a set of assumptions about ASC gain and loss from the sampled sites, and also on the directionality of ASC trafficking pathways. Most notably, modeling results suggest that differences in ASC fate and trafficking patterns reflect the site of formation and the expressed antibody class. Essentially all early IgA ASCs in the MLN migrated to spleen or lung, whereas cell death was likely the major reason for IgM and IgG ASC loss from the MLN. In contrast, the spleen contributed most of the IgM and IgG ASCs that migrated to the lung, but essentially none of the IgA ASCs. This finding points to a critical role for regional lymph nodes such as the MLN in the rapid generation of IgA ASCs that seed the lung. Results for the MLN also suggest that ASC death is a significant early feature of the B cell response. Overall, our analysis is consistent with accepted concepts in many regards, but it also indicates novel features of the B cell response to influenza that warrant further investigation. link: http://identifiers.org/pubmed/25171166

The hormone gibberellin (GA) is a key regulator of plant growth. Many of the components of the gibberellin signal transduction [e.g., GIBBERELLIN INSENSITIVE DWARF 1 (GID1) and DELLA], biosynthesis [e.g., GA 20-oxidase (GA20ox) and GA3ox], and deactivation pathways have been identified. Gibberellin binds its receptor, GID1, to form a complex that mediates the degradation of DELLA proteins. In this way, gibberellin relieves DELLA-dependent growth repression. However, gibberellin regulates expression of GID1, GA20ox, and GA3ox, and there is also evidence that it regulates DELLA expression. In this paper, we use integrated mathematical modeling and experiments to understand how these feedback loops interact to control gibberellin signaling. Model simulations are in good agreement with in vitro data on the signal transduction and biosynthesis pathways and in vivo data on the expression levels of gibberellin-responsive genes. We find that GA-GID1 interactions are characterized by two timescales (because of a lid on GID1 that can open and close slowly relative to GA-GID1 binding and dissociation). Furthermore, the model accurately predicts the response to exogenous gibberellin after a number of chemical and genetic perturbations. Finally, we investigate the role of the various feedback loops in gibberellin signaling. We find that regulation of GA20ox transcription plays a significant role in both modulating the level of endogenous gibberellin and generating overshoots after the removal of exogenous gibberellin. Moreover, although the contribution of other individual feedback loops seems relatively small, GID1 and DELLA transcriptional regulation acts synergistically with GA20ox feedback. link: http://identifiers.org/pubmed/22523240

The metabolism of microorganisms is regulated through two main mechanisms: changes of enzyme capacities as a consequence of gene expression modulation (“hierarchical control”) and changes of enzyme activities through metabolite-enzyme interactions. An increasing body of evidence indicates that hierarchical control is insufficient to explain metabolic behaviors, but the system-wide impact of metabolic regulation remains largely uncharacterized. To clarify its role, we developed and validated a detailed kinetic model of Escherichia coli central metabolism that links growth to environment. Metabolic control analyses confirm that the control is widely distributed across the network and highlight strong interconnections between all the pathways. Exploration of the model solution space reveals that several robust properties emerge from metabolic regulation, from the molecular level (e.g. homeostasis of total metabolite pool) to the overall cellular physiology (e.g. coordination of carbon uptake, catabolism, energy and redox production, and growth), while allowing a large degree of flexibility at most individual metabolic steps. These properties have important physiological implications for E. coli and significantly expand the self-regulating capacities of its metabolism. link: http://identifiers.org/doi/10.1371/journal.pcbi.1005396

Overflow metabolism refers to the production of seemingly wasteful by-products by cells during growth on glucose even when oxygen is abundant. Two theories have been proposed to explain acetate overflow in Escherichia coli - global control of the central metabolism and local control of the acetate pathway - but neither accounts for all observations. Here, we develop a kinetic model of E. coli metabolism that quantitatively accounts for observed behaviors and successfully predicts the response of E. coli to new perturbations. We reconcile these theories and clarify the origin, control and regulation of the acetate flux. We also find that, in turns, acetate regulates glucose metabolism by coordinating the expression of glycolytic and TCA genes. Acetate should not be considered a wasteful end-product since it is also a co-substrate and a global regulator of glucose metabolism in E. coli. This has broad implications for our understanding of overflow metabolism. link: http://identifiers.org/doi/10.1101/2020.08.18.255356

BACKGROUND: Solventogenic clostridia offer a sustainable alternative to petroleum-based production of butanol–an important chemical feedstock and potential fuel additive or replacement. C. beijerinckii is an attractive microorganism for strain design to improve butanol production because it (i) naturally produces the highest recorded butanol concentrations as a byproduct of fermentation; and (ii) can co-ferment pentose and hexose sugars (the primary products from lignocellulosic hydrolysis). Interrogating C. beijerinckii metabolism from a systems viewpoint using constraint-based modeling allows for simulation of the global effect of genetic modifications. RESULTS: We present the first genome-scale metabolic model (iCM925) for C. beijerinckii, containing 925 genes, 938 reactions, and 881 metabolites. To build the model we employed a semi-automated procedure that integrated genome annotation information from KEGG, BioCyc, and The SEED, and utilized computational algorithms with manual curation to improve model completeness. Interestingly, we found only a 34% overlap in reactions collected from the three databases–highlighting the importance of evaluating the predictive accuracy of the resulting genome-scale model. To validate iCM925, we conducted fermentation experiments using the NCIMB 8052 strain, and evaluated the ability of the model to simulate measured substrate uptake and product production rates. Experimentally observed fermentation profiles were found to lie within the solution space of the model; however, under an optimal growth objective, additional constraints were needed to reproduce the observed profiles–suggesting the existence of selective pressures other than optimal growth. Notably, a significantly enriched fraction of actively utilized reactions in simulations–constrained to reflect experimental rates–originated from the set of reactions that overlapped between all three databases (P = 3.52 × 10-9, Fisher's exact test). Inhibition of the hydrogenase reaction was found to have a strong effect on butanol formation–as experimentally observed. CONCLUSIONS: Microbial production of butanol by C. beijerinckii offers a promising, sustainable, method for generation of this important chemical and potential biofuel. iCM925 is a predictive model that can accurately reproduce physiological behavior and provide insight into the underlying mechanisms of microbial butanol production. As such, the model will be instrumental in efforts to better understand, and metabolically engineer, this microorganism for improved butanol production. link: http://identifiers.org/pubmed/21846360

Iron is essential for all known life due to its redox properties; however, these same properties can also lead to its toxicity in overload through the production of reactive oxygen species. Robust systemic and cellular control are required to maintain safe levels of iron, and the liver seems to be where this regulation is mainly located. Iron misregulation is implicated in many diseases, and as our understanding of iron metabolism improves, the list of iron-related disorders grows. Recent developments have resulted in greater knowledge of the fate of iron in the body and have led to a detailed map of its metabolism; however, a quantitative understanding at the systems level of how its components interact to produce tight regulation remains elusive. A mechanistic computational model of human liver iron metabolism, which includes the core regulatory components, is presented here. It was constructed based on known mechanisms of regulation and on their kinetic properties, obtained from several publications. The model was then quantitatively validated by comparing its results with previously published physiological data, and it is able to reproduce multiple experimental findings. A time course simulation following an oral dose of iron was compared to a clinical time course study and the simulation was found to recreate the dynamics and time scale of the systems response to iron challenge. A disease state simulation of haemochromatosis was created by altering a single reaction parameter that mimics a human haemochromatosis gene (HFE) mutation. The simulation provides a quantitative understanding of the liver iron overload that arises in this disease. This model supports and supplements understanding of the role of the liver as an iron sensor and provides a framework for further modelling, including simulations to identify valuable drug targets and design of experiments to improve further our knowledge of this system. link: http://identifiers.org/pubmed/24244122

INTRODUCTION: The therapeutic potential of a hemostatic agent can be assessed by investigating its effects on the quantitative parameters of thrombin generation. For recombinant activated factor VII (rFVIIa)–a promising hemostasis-inducing biologic–experimental studies addressing its effects on thrombin generation yielded disparate results. To elucidate the inherent ability of rFVIIa to modulate thrombin production, it is necessary to identify rFVIIa-induced effects that are compatible with the available biochemical knowledge about thrombin generation mechanisms. MATERIALS AND METHODS: The existing body of knowledge about coagulation biochemistry can be rigorously represented by a computational model that incorporates the known reactions and parameter values constituting the biochemical network. We used a thoroughly validated numerical model to generate activated factor VII (FVIIa) titration curves in the cases of normal blood composition, hemophilia A and B blood, blood lacking factor VII, blood lacking tissue factor pathway inhibitor, and diluted blood. We utilized the generated curves to perform systematic fold-change analyses for five quantitative parameters characterizing thrombin accumulation. RESULTS: The largest fold changes induced by increasing FVIIa concentration were observed for clotting time, thrombin peak time, and maximum slope of the thrombin curve. By contrast, thrombin peak height was much less affected by FVIIa titrations, and the area under the thrombin curve stayed practically unchanged. Comparisons with experimental data demonstrated that the computationally derived patterns can be observed in vitro. CONCLUSIONS: rFVIIa modulates thrombin generation primarily by accelerating the process, without significantly affecting the total amount of generated thrombin. link: http://identifiers.org/pubmed/21641634

BACKGROUND: Hypothermia, which can result from tissue hypoperfusion, body exposure, and transfusion of cold resuscitation fluids, is a major factor contributing to coagulopathy of trauma and surgery. Despite considerable efforts, the mechanisms of hypothermia-induced blood coagulation impairment have not been fully understood. We introduce a kinetic modeling approach to investigate the effects of hypothermia on thrombin generation. METHODS: We extended a validated computational model to predict and analyze the impact of low temperatures (with or without concomitant blood dilution) on thrombin generation and its quantitative parameters. The computational model reflects the existing knowledge about the mechanistic details of thrombin generation biochemistry. We performed the analysis for an "average" subject, as well as for 472 subjects in the control group of the Leiden Thrombophilia Study. RESULTS: We computed and analyzed thousands of kinetic curves characterizing the generation of thrombin and the formation of the thrombin-antithrombin complex (TAT). In all simulations, hypothermia in the temperature interval 31°C to 36°C progressively slowed down thrombin generation, as reflected by clotting time, thrombin peak time, and prothrombin time, which increased in all subjects (P < 10(-5)). Maximum slope of the thrombin curve was progressively decreased, and the area under the thrombin curve was increased in hypothermia (P < 10(-5)); thrombin peak height remained practically unaffected. TAT formation was noticeably delayed (P < 10(-5)), but the final TAT levels were not significantly affected. Hypothermia-induced fold changes in the affected thrombin generation parameters were larger for lower temperatures, but were practically independent of the parameter itself and of the subjects' clotting factor composition, despite substantial variability in the subject group. Hypothermia and blood dilution acted additively on the thrombin generation parameters. CONCLUSIONS: We developed a general computational strategy that can be used to simulate the effects of changing temperature on the kinetics of biochemical systems and applied this strategy to analyze the effects of hypothermia on thrombin generation. We found that thrombin generation can be noticeably impaired in subjects with different blood plasma composition even in moderate hypothermia. Our work provides mechanistic support to the notion that thrombin generation impairment may be a key factor in coagulopathy induced by hypothermia and complicated by blood plasma dilution. link: http://identifiers.org/pubmed/23868891

Current mechanistic knowledge of protein interactions driving blood coagulation has come largely from experiments with simple synthetic systems, which only partially represent the molecular composition of human blood plasma. Here, we investigate the ability of the suggested molecular mechanisms to account for fibrin generation and degradation kinetics in diverse, physiologically relevant in vitro systems. We represented the protein interaction network responsible for thrombin generation, fibrin formation, and fibrinolysis as a computational kinetic model and benchmarked it against published and newly generated data reflecting diverse experimental conditions. We then applied the model to investigate the ability of fibrinogen and a recently proposed prothrombin complex concentrate composition, PCC-AT (a combination of the clotting factors II, IX, X, and antithrombin), to restore normal thrombin and fibrin generation in diluted plasma. The kinetic model captured essential features of empirically detected effects of prothrombin, fibrinogen, and thrombin-activatable fibrinolysis inhibitor titrations on fibrin formation and degradation kinetics. Moreover, the model qualitatively predicted the impact of tissue factor and tPA/tenecteplase level variations on the fibrin output. In the majority of considered cases, PCC-AT combined with fibrinogen accurately approximated both normal thrombin and fibrin generation in diluted plasma, which could not be accomplished by fibrinogen or PCC-AT acting alone. We conclude that a common network of protein interactions can account for key kinetic features characterizing fibrin accumulation and degradation in human blood plasma under diverse experimental conditions. Combined PCC-AT/fibrinogen supplementation is a promising strategy to reverse the deleterious effects of dilution-induced coagulopathy associated with traumatic bleeding. link: http://identifiers.org/pubmed/24958246

Acidosis, a frequent complication of trauma and complex surgery, results from tissue hypoperfusion and IV resuscitation with acidic fluids. While acidosis is known to inhibit the function of distinct enzymatic reactions, its cumulative effect on the blood coagulation system is not fully understood. Here, we use computational modeling to test the hypothesis that acidosis delays and reduces the amount of thrombin generation in human blood plasma. Moreover, we investigate the sensitivity of different thrombin generation parameters to acidosis, both at the individual and population level.We used a kinetic model to simulate and analyze the generation of thrombin and thrombin-antithrombin complexes (TAT), which were the end points of this study. Large groups of temporal thrombin and TAT trajectories were simulated and used to calculate quantitative parameters, such as clotting time (CT), thrombin peak time, maximum slope of the thrombin curve, thrombin peak height, area under the thrombin trajectory (AUC), and prothrombin time. The resulting samples of parameter values at different pH levels were compared to assess the acidosis-induced effects. To investigate intersubject variability, we parameterized the computational model using the data on clotting factor composition for 472 subjects from the Leiden Thrombophilia Study. To compare acidosis-induced relative parameter changes in individual ("virtual") subjects, we estimated the probabilities of relative change patterns by counting the pattern occurrences in our virtual subjects. Distribution overlaps for thrombin generation parameters at distinct pH levels were quantified using the Bhattacharyya coefficient.Acidosis in the range of pH 6.9 to 7.3 progressively increased CT, thrombin peak time, AUC, and prothrombin time, while decreasing maximum slope of the thrombin curve and thrombin peak height (P < 10). Acidosis delayed the onset and decreased the amount of TAT generation (P < 10). As a measure of intrasubject variability, maximum slope of the thrombin curve and CT displayed the largest and second-largest acidosis-induced relative changes, and AUC displayed the smallest relative changes among all thrombin generation parameters in our virtual subject group (1-sided 95% lower confidence limit on the fraction of subjects displaying the patterns, 0.99). As a measure of intersubject variability, the overlaps between the maximum slope of the thrombin curve distributions at acidotic pH levels with the maximum slope of the thrombin curve distribution at physiological pH level systematically exceeded analogous distribution overlaps for CT, thrombin peak time, and prothrombin time.Acidosis affected all quantitative parameters of thrombin and TAT generation. While maximum slope of the thrombin curve showed the highest sensitivity to acidosis at the individual-subject level, it may be outperformed by CT, thrombin peak time, and prothrombin time as an indicator of acidosis at the subject-group level. link: http://identifiers.org/pubmed/25839182

The use of prothrombin complex concentrates in trauma- and surgery-induced coagulopathy is complicated by the possibility of thromboembolic events. To explore the effects of these agents on thrombin generation (TG), we investigated combinations of coagulation factors equivalent to 3- and 4-factor prothrombin complex concentrates with and without added antithrombin (AT), as well as recombinant factor VIIa (rFVIIa), in a dilutional model. These data were then used to develop a computational model to test whether such a model could predict the TG profiles of these agents used to treat dilutional coagulopathy.We measured TG in plasma collected from 10 healthy volunteers using Calibrated Automated Thrombogram. TG measurements were performed in undiluted plasma, 3-fold saline-diluted plasma, and diluted plasma supplemented with the following factors: rFVIIa (group rFVIIa); factors (F)II, FIX, FX, and AT (group "combination of coagulation factors" [CCF]-AT); or FII, FVII, FIX, and FX (group CCF-FVII). We extended an existing computational model of TG to include additional reactions that impact the Calibrated Automated Thrombogram readout. We developed and applied a computational strategy to train the model using only a subset of the obtained TG data and used the remaining data for model validation.rFVIIa decreased lag time and the time to thrombin peak generation beyond their predilution levels (P < 0.001) but did not restore normal thrombin peak height (P < 0.001). CCF-FVII supplementation decreased lag time (P = 0.034) and thrombin peak time (P < 0.001) and increased both peak height (P < 0.001) and endogenous thrombin potential (P = 0.055) beyond their predilution levels. CCF-AT supplementation in diluted plasma resulted in an improvement in TG without causing the exaggerated effects of rFVIIa and CCF-FVII supplementation. The differences between the effects of CCF-AT and supplementation with rFVIIa and CCF-FVII were significant for lag time (P < 0.001 and P = 0.005, respectively), time to thrombin peak (P < 0.001 and P = 0.004, respectively), velocity index (P < 0.001 and P = 0.019, respectively), thrombin peak height (P < 0.001 for both comparisons), and endogenous thrombin potential (P = 0.034 and P = 0.019, respectively). The computational model generated subject-specific predictions and identified typical patterns of TG improvement.In this study of the effects of hemodilution, CCF-AT supplementation improved the dilution-impaired plasma TG potential in a more balanced way than either rFVIIa alone or CCF-FVII supplementation. Predictive computational modeling can guide plasma dilution/supplementation experiments. link: http://identifiers.org/pubmed/27541717

We present and analyze a model for the interaction of human immunodeficiency virus type 1 (HIV-1) with target cells that includes a time delay between initial infection and the formation of productively infected cells. Assuming that the variation among cells with respect to this 'intracellular' delay can be approximated by a gamma distribution, a high flexible distribution that can mimic a variety of biologically plausible delays, we provide analytical solutions for the expected decline in plasma virus concentration after the initiation of antiretroviral therapy with one or more protease inhibitors. We then use the model to investigate whether the parameters that characterize viral dynamics can be identified from biological data. Using non-linear least-squares regression to fit the model to simulated data in which the delays conform to a gamma distribution, we show that good estimates for free viral clearance rates, infected cell death rates, and parameters characterizing the gamma distribution can be obtained. For simulated data sets in which the delays were generated using other biologically plausible distributions, reasonably good estimates for viral clearance rates, infected cell death rates, and mean delay times can be obtained using the gamma-delay model. For simulated data sets that include added simulated noise, viral clearance rate estimates are not as reliable. If the mean intracellular delay is known, however, we show that reasonable estimates for the viral clearance rate can be obtained by taking the harmonic mean of viral clearance rate estimates from a group of patients. These results demonstrate that it is possible to incorporate distributed intracellular delays into existing models for HIV dynamics and to use these refined models to estimate the half-life of free virus from data on the decline in HIV-1 RNA following treatment. link: http://identifiers.org/pubmed/9780612

BACKGROUND: Alzheimer's disease (AD) is the most common cause of dementia among the elderly. To clarify pathogenesis of AD, thousands of reports have been accumulating. However, knowledge of signaling pathways in the field of AD has not been compiled as a database before. DESCRIPTION: Here, we have constructed a publicly available pathway map called "AlzPathway" that comprehensively catalogs signaling pathways in the field of AD. We have collected and manually curated over 100 review articles related to AD, and have built an AD pathway map using CellDesigner. AlzPathway is currently composed of 1347 molecules and 1070 reactions in neuron, brain blood barrier, presynaptic, postsynaptic, astrocyte, and microglial cells and their cellular localizations. AlzPathway is available as both the SBML (Systems Biology Markup Language) map for CellDesigner and the high resolution image map. AlzPathway is also available as a web service (online map) based on Payao system, a community-based, collaborative web service platform for pathway model curation, enabling continuous updates by AD researchers. CONCLUSIONS: AlzPathway is the first comprehensive map of intra, inter and extra cellular AD signaling pathways which can enable mechanistic deciphering of AD pathogenesis. The AlzPathway map is accessible at http://alzpathway.org/. link: http://identifiers.org/pubmed/22647208

A type of cancer which originates from the breast tissue is referred to as breast cancer. Globally, it is the most common cause of death in women. Treatments such as radiotherapy, chemotherapy, hormone therapy, immunotherapy, and gene therapy are the main strategies in the fight against breast cancer. The present study aims at investigating the effects of the combined radiotherapy and chemotherapy as a way to treat breast cancer, and different treatment approaches are incorporated into the model. Also, the model is fitted to data on patients with breast cancer in Tanzania. We determine new treatment strategies, and finally, we show that when sufficient amount of chemotherapy and radiotherapy with a low decay rate is used, the drug will be significantly more effective in combating the disease while health cells remain above the threshold. link: http://identifiers.org/pubmed/31611927

The precise coordination of growth and proliferation has a universal prevalence in cell homeostasis. As a prominent property, cell size is modulated by the coordination between these processes in bacterial, yeast, and mammalian cells, but the underlying molecular mechanisms are largely unknown. Here, we show that multifunctional chaperone systems play a concerted and limiting role in cell-cycle entry, specifically driving nuclear accumulation of the G1 Cdk-cyclin complex. Based on these findings, we establish and test a molecular competition model that recapitulates cell-cycle-entry dependence on growth rate. As key predictions at a single-cell level, we show that availability of the Ydj1 chaperone and nuclear accumulation of the G1 cyclin Cln3 are inversely dependent on growth rate and readily respond to changes in protein synthesis and stress conditions that alter protein folding requirements. Thus, chaperone workload would subordinate Start to the biosynthetic machinery and dynamically adjust proliferation to the growth potential of the cell. link: http://identifiers.org/pubmed/30988162

The cholesterol biosynthesis pathway has recently been shown to play an important role in the innate immune response to viral infection with host protection occurring through a coordinate down regulation of the enzymes catalysing each metabolic step. In contrast, statin based drugs, which form the principle pharmaceutical agents for decreasing the activity of this pathway, target a single enzyme. Here, we build an ordinary differential equation model of the cholesterol biosynthesis pathway in order to investigate how the two regulatory strategies impact upon the behaviour of the pathway. We employ a modest set of assumptions: that the pathway operates away from saturation, that each metabolite is involved in multiple cellular interactions and that mRNA levels reflect enzyme concentrations. Using data taken from primary bone marrow derived macrophage cells infected with murine cytomegalovirus or treated with IFNγ, we show that, under these assumptions, coordinate down-regulation of enzyme activity imparts a graduated reduction in flux along the pathway. In contrast, modelling a statin-like treatment that achieves the same degree of down-regulation in cholesterol production, we show that this delivers a step change in flux along the pathway. The graduated reduction mediated by physiological coordinate regulation of multiple enzymes supports a mechanism that allows a greater level of specificity, altering cholesterol levels with less impact upon interactions branching from the pathway, than pharmacological step reductions. We argue that coordinate regulation is likely to show a long-term evolutionary advantage over single enzyme regulation. Finally, the results from our models have implications for future pharmaceutical therapies intended to target cholesterol production with greater specificity and fewer off target effects, suggesting that this can be achieved by mimicking the coordinated down-regulation observed in immunological responses. link: http://identifiers.org/pubmed/22664637

White adipocytes are specialized for energy storage, whereas brown adipocytes are specialized for energy expenditure. Explicating this difference can help identify therapeutic targets for obesity. A common tool to assess metabolic differences between such cells is the Seahorse Extracellular Flux (XF) Analyzer, which measures oxygen consumption and media acidification in the presence of different substrates and perturbagens. Here, we integrate the Analyzer's metabolic profile from human white and brown adipocytes with a genome-scale metabolic model to predict flux differences across the metabolic map. Predictions matched experimental data for the metabolite 4-aminobutyrate, the protein ABAT, and the fluxes for glucose, glutamine, and palmitate. We also uncovered a difference in how adipocytes dispose of nitrogenous waste, with brown adipocytes secreting less ammonia and more urea than white adipocytes. Thus, the method and software we developed allow for broader metabolic phenotyping and provide a distinct approach to uncovering metabolic differences. link: http://identifiers.org/pubmed/29241534

A better understanding of human metabolism and its relationship with diseases is an important task in human systems biology studies. In this paper, we present a high-quality human metabolic network manually reconstructed by integrating genome annotation information from different databases and metabolic reaction information from literature. The network contains nearly 3000 metabolic reactions, which were reorganized into about 70 human-specific metabolic pathways according to their functional relationships. By analysis of the functional connectivity of the metabolites in the network, the bow-tie structure, which was found previously by structure analysis, is reconfirmed. Furthermore, the distribution of the disease related genes in the network suggests that the IN (substrates) subset of the bow-tie structure has more flexibility than other parts. link: http://identifiers.org/pubmed/17882155

A better understanding of human metabolism and its relationship with diseases is an important task in human systems biology studies. In this paper, we present a high-quality human metabolic network manually reconstructed by integrating genome annotation information from different databases and metabolic reaction information from literature. The network contains nearly 3000 metabolic reactions, which were reorganized into about 70 human-specific metabolic pathways according to their functional relationships. By analysis of the functional connectivity of the metabolites in the network, the bow-tie structure, which was found previously by structure analysis, is reconfirmed. Furthermore, the distribution of the disease related genes in the network suggests that the IN (substrates) subset of the bow-tie structure has more flexibility than other parts. link: http://identifiers.org/pubmed/17882155

After whole-genome duplication (WGD), deletions return most loci to single copy. However, duplicate loci may survive through selection for increased dosage. Here, we show how the WGD increased copy number of some glycolytic genes could have conferred an almost immediate selective advantage to an ancestor of Saccharomyces cerevisiae, providing a rationale for the success of the WGD. We propose that the loss of other redundant genes throughout the genome resulted in incremental dosage increases for the surviving duplicated glycolytic genes. This increase gave post-WGD yeasts a growth advantage through rapid glucose fermentation; one of this lineage's many adaptations to glucose-rich environments. Our hypothesis is supported by data from enzyme kinetics and comparative genomics. Because changes in gene dosage follow directly from post-WGD deletions, dosage selection can confer an almost instantaneous benefit after WGD, unlike neofunctionalization or subfunctionalization, which require specific mutations. We also show theoretically that increased fermentative capacity is of greatest advantage when glucose resources are both large and dense, an observation potentially related to the appearance of angiosperms around the time of WGD. link: http://identifiers.org/pubmed/17667951

After whole-genome duplication (WGD), deletions return most loci to single copy. However, duplicate loci may survive through selection for increased dosage. Here, we show how the WGD increased copy number of some glycolytic genes could have conferred an almost immediate selective advantage to an ancestor of Saccharomyces cerevisiae, providing a rationale for the success of the WGD. We propose that the loss of other redundant genes throughout the genome resulted in incremental dosage increases for the surviving duplicated glycolytic genes. This increase gave post-WGD yeasts a growth advantage through rapid glucose fermentation; one of this lineage's many adaptations to glucose-rich environments. Our hypothesis is supported by data from enzyme kinetics and comparative genomics. Because changes in gene dosage follow directly from post-WGD deletions, dosage selection can confer an almost instantaneous benefit after WGD, unlike neofunctionalization or subfunctionalization, which require specific mutations. We also show theoretically that increased fermentative capacity is of greatest advantage when glucose resources are both large and dense, an observation potentially related to the appearance of angiosperms around the time of WGD. link: http://identifiers.org/pubmed/17667951

MOTIVATION: It is a question of whether the supramolecular organization of the protein complex has an impact on its function, or not. In the case of the photosystem II (PSII), water splitting might be influenced by cooperation of the PSIIs. Since PSII is the source of the atmospheric oxygen and because better understanding of the water splitting may contribute to the effective use of water as an alternative energy source, possible cooperation should be analyzed and discussed. RESULTS: We suggest that the dimeric organization of the PSII induces cooperation in the water splitting. We show that the model of monomeric PSII is unable to produce the oxygen after the second short flash (associated with the double turnover of the PSII), in contrast to the experimental data and model of dimeric PSII with considered cooperation. On the basis of this fact and partially from the support from other studies, we concluded that the double turnover of the PSII induced by short flashes might be caused by the cooperation in the water splitting. We further discuss a possibility that the known pathway of the electron transport through the PSII might be incomplete and besides D1-Y161, other cofactor which is able to oxidize the special chlorophyll pair (P680) must be considered in the monomeric PSII to explain the oxygen production after the second short flash. AVAILABILITY: Commented SBML codes (.XML files) of the monomeric and dimeric PSII models will be available (at the time of publication) in the BioModels database (www.ebi.ac.uk/biomodels). link: http://identifiers.org/pubmed/18845578

MOTIVATION: It is a question of whether the supramolecular organization of the protein complex has an impact on its function, or not. In the case of the photosystem II (PSII), water splitting might be influenced by cooperation of the PSIIs. Since PSII is the source of the atmospheric oxygen and because better understanding of the water splitting may contribute to the effective use of water as an alternative energy source, possible cooperation should be analyzed and discussed. RESULTS: We suggest that the dimeric organization of the PSII induces cooperation in the water splitting. We show that the model of monomeric PSII is unable to produce the oxygen after the second short flash (associated with the double turnover of the PSII), in contrast to the experimental data and model of dimeric PSII with considered cooperation. On the basis of this fact and partially from the support from other studies, we concluded that the double turnover of the PSII induced by short flashes might be caused by the cooperation in the water splitting. We further discuss a possibility that the known pathway of the electron transport through the PSII might be incomplete and besides D1-Y161, other cofactor which is able to oxidize the special chlorophyll pair (P680) must be considered in the monomeric PSII to explain the oxygen production after the second short flash. AVAILABILITY: Commented SBML codes (.XML files) of the monomeric and dimeric PSII models will be available (at the time of publication) in the BioModels database (www.ebi.ac.uk/biomodels). link: http://identifiers.org/pubmed/18845578

Glucose addition and subsequent run-out experiments were compared to simulations with a detailed glycolytic model of Lactococcus lactis. The model was constructed largely on bases of enzyme kinetic data taken from literature and not adjusted for the specific simulations shown here. Upon glucose depletion a rapid increase in PEP, inorganic phosphate and a gradual decrease in fructose 1,6-bisphosphate (FBP) were measured and predicted by simulation. The dynamic changes in these and other intermediate concentrations as measured in the experiments were well predicted by the kinetic model. link: http://identifiers.org/pubmed/12241048

Quorum sensing (QS) is an important determinant of bacterial phenotype. Many cell functions are regulated by intricate and multimodal QS signal transduction processes. The LuxS/AI-2 QS system is highly conserved among Eubacteria and AI-2 is reported as a 'universal' signal molecule. To understand the hierarchical organization of AI-2 circuitry, a comprehensive approach incorporating stochastic simulations was developed. We investigated the synthesis, uptake, and regulation of AI-2, developed testable hypotheses, and made several discoveries: (1) the mRNA transcript and protein levels of AI-2 synthases, Pfs and LuxS, do not contribute to the dramatically increased level of AI-2 found when cells are grown in the presence of glucose; (2) a concomitant increase in metabolic flux through this synthesis pathway in the presence of glucose only partially accounts for this difference. We predict that 'high-flux' alternative pathways or additional biological steps are involved in AI-2 synthesis; and (3) experimental results validate this hypothesis. This work demonstrates the utility of linking cell physiology with systems-based stochastic models that can be assembled de novo with partial knowledge of biochemical pathways. link: http://identifiers.org/pubmed/17170762

Modulated immune signal (CD14-TLR and TNF) in leishmaniasis can be linked to EGFR pathway involved in wound healing, through crosstalk points. This signaling network can be further linked to a synthetic gene circuit acting as a positive feedback loop to elicit a synchronized intercellular communication among the immune cells which may contribute to a better understanding of signaling dynamics in leishmaniasis.Network reconstruction with positive feedback loop, simulation (ODE 15s solver) and sensitivity analysis of CD14-TLR, TNF and EGFR was done in SimBiology (MATLAB 7.11.1). Cytoscape and adjacency matrix were used to calculate network topology. PCA was extracted by using sensitivity coefficient in MATLAB. Model reduction was done using time, flux and sensitivity score.Network has five crosstalk points: NIK, IκB-NFκB and MKK (4/7, 3/6, 1/2) which show high flux and sensitivity. PI3K in EGFR pathway shows high flux and sensitivity. PCA score was high for cytoplasmic ERK1/2, PI3K, Atk, STAT1/3 and nuclear JNK. Of the 125 parameters, 20% are crucial as deduced by model reduction.EGFR can be linked to CD14-TLR and TNF through the MAPK crosstalk points. These pathways may be controlled through Ras and Raf that lie upstream of signaling components ERK ½ (c) and JNK (n) that have a high PCA score via a synthetic gene circuit for activating cell-cell communication to elicit an inflammatory response. Also a disease resolving effect may be achieved through PI3K in the EGFR pathway.The reconstructed signaling network can be linked to a gene circuit with a positive feedback loop, for cell-cell communication resulting in synchronized response in the immune cell population, for disease resolving effect in leishmaniasis. link: http://identifiers.org/pubmed/23994140